Search Documents

By Paolo Farina

Slope monitoring radar has emerged in the last ten years as a leading edge tool for monitoring movements in natural and engineered slopes, thanks to its ability to rapidly measure slope movements with millimetric accuracy over wide areas in any weather conditions, obviating the need to install artificial reflectors. In recent years, slope monitoring radar has experienced significant improvements in its performances, mainly thanks to the use of different radar technology, such as Synthetic Aperture Radar (SAR), able to overcome some of the limitations of the first generation of radar and to provide higher spatial resolution, longer working distances and faster acquisition time. One of the leading applications of this technology is today represented by slope monitoring within the surface mining industry. Prominent international mining groups currently employ slope monitoring radar, based either on Real Aperture Radar (RAR) and SAR technology, for active and background monitoring of open-pit walls. The main technical features of typical slope monitoring radar based on SAR technology along with their possible applications are presented in this paper along with recent case studies of successful slope monitoring in open pit mines.

Jun 1, 2012

By Hamidreza Soltanzadeh

The injection of CO2 into a reservoir or aquifer can potentially result in deformation, induced fracturing and fault reactivation in the injection zone and/or the rocks that surround it. This paper reviews the data collection and geomechanical analyses performed to assess these potential hazards for the Heartland Area Redwater Project (HARP), in which CO2 sequestration is being considered in a carbonate (Leduc Formation) reef. A three-dimensional mechanical earth model (3D MEM) was developed using available geophysical log data from several wells. The geomechanical response to CO2 injection was studied using a finite element numerical model based on the 3D MEM and pore pressure distributions predicted using a reservoir simulator. The modeling results suggest that the potential for induced fracturing and fault reactivation is low and the maximum predicted ground surface heave is less than 2 cm. Sensitivity analyses conducted with a 2D geomechanical model that was extracted from the 3D model confirm that this conclusion holds true even when a reasonable measure of uncertainty in rock mechanical properties, ambient pore pressures and horizontal stress magnitudes is accounted for.

Jun 1, 2012

By SeyedBijan Mahbaz

Development of Geomechanical Earth Models (GEM) is a prerequisite to reservoir management, involving issues such as stress path analysis, prediction of induced seismicity and sanding initiation, and coupled flow-stress modeling that leads to more realistic predictions of oil rates and rock response. Part of GEM development is the stipulation of the 3-D lithostratigraphic disposition (geometry) which is used to choose the geomechanical units. The digital geological model is then populated with quantitative and semi-quantitative data on rock properties, fluid properties, saturations, initial conditions (p, T, [_]) and other information. Both fluid and solid properties are needed this requires integration of geomechanics skills with geophysical and rock physics skills, and this article deals with a subset of this integration process, with an example from two wells in southwest Iran.Rock physics approaches are used to predict reservoir rock and fluid properties with wellbore and 3-D seismic data, which are also used to interpolate rock mechanics properties to the inter-well regions. The first step in a generic rock physics approach is diagnostic and involves the introduction of a suitable Rock Physics model.Acoustic impedance (AI), shear impedance (SI), and density (_) are the fundamental rock properties usually derived using AVO equations, and these parameters, along with attenuation coefficients, can then be correlated to static rock mechanics properties that define deformability and perhaps even strength. In this study, Poisson Impedance is introduced as a new seismic attribute used for lithology and fluid differentiation in the FahliyanFormation (reservoir rock) in two oil wells in the southwest of Iran. To do this, a novel rock physics model was developed for this naturally fractured carbonate reservoir.Currently, these data are also being used to correlate to a limited rock mechanics data base to develop a more extensive GEM to serve as input to more general geomechanical and reservoir modeling tasks.

Jun 1, 2012

By Jeffrey D. Oke

While significant progress in underground construction engineering has occurred in the past few decades, a number of fundamental deficiencies persist regarding the recognition, prediction, and mitigation of underground construction hazards associated with complex geological conditions in rock and soil works, such as tunnelling. Many of these shortcomings could be addressed by improving geomechanical models which, in turn, would provide more realistic frameworks for determining input parameters for design analysis. This would improve predictions of construction efficiency for tunnelling and the formulation of more accurate modelling techniques (2D, 3D and 4D) required to optimize the design of engineered rock structures. This research focuses primarily on the issues related to forepole temporary tunnel support within weak rock masses. This research attempts to determine how best to numerically model forepole support. The behaviour and mechanics of forepole support has not been quantified and/or optimized. As well, there is a lack of in-situ and/or laboratory forepole and weak rockmass interaction parameters available. In terms of numerical modelling, there is a lack of understanding of the use of proper interaction parameters associated with the forepole and the rockmass. The parametric study that has been conducted by the authors demonstrates the influence of each of the forepole-rockmass interaction parameters and their sensitivity. Introduction

Jun 1, 2012

By Kathy S. Kalenchuk

Damages have been observed in shafts, backfill infrastructure and the slick line at the Leeville Mine. Shaft damage (shearing and spalling of the concrete liner) is attributed to mine-induced rock mass deformations concentrated where regional faults intersect the shafts. Aggregate bin, mixer station and conveyor drive damage (shotcrete shearing, heave of concrete floor slabs and bending of a vertical steel beams) is attributed to rock mass squeezing induced by elevated magnitudes of the vertically orientated major principal stresses. Separation of the slick line has occurred do to extension of the stressshadowed rock mass immediately above the mined out areas. Three-dimensional numerical models have been developed and calibrated to simulated rock mass deformations with magnitudes similar to those believed to have produced the observed damages. By stepping through the mine sequence, models reproduce deformations at the same mining stage as when damage has occurred. Calibrated models have been used to forward project further damage over the remaining mine life.

Jun 1, 2012

By Karina Forster

Open stope design depends on the ore body geometry, loading conditions and rock mass properties to estimate stope dilution and ground support requirements. Rock mass conditions for design primarily rely upon estimates of well established rock classification systems. Many sources of data exist for obtaining estimates for these rock mass classification values which include core logging, geotechnical mapping, geological mapping as well as linking one classification system to another through the use of empirically derived equations. The extent of mapping coverage can also vary significantly based on the data collection technique. This paper proposes an approach where individual components of rock mass classification systems are collected by the different available data collection techniques. Data from different rock types are also considered separately based on their alteration characteristics and general observation products. Data from the Cameco Corporation Eagle Point Mine is presented and an approach is outlined for obtaining average representative classification values for stope hanging wall characterization. The effectiveness of this approach is demonstrated through the analysis of a data base of hanging wall dilution information.

Jun 1, 2012

By Somchet Vongpaisal

This paper presents the R&D results on the development of a new knowledge-based system for predicting stope instability and associated unplanned dilution in deep underground mining operations. A new prototype was developed based on databases generated from 3D non-linear elasto-plastic finite element stress analysis, heuristic rules, and case histories. This prototype aims to assist mine operators to make decisions on strategic mine planning and feasibility in a quick and efficient manner thereby minimizing risk on ground failure and optimizing mining costs.

Jun 1, 2012

By Will S. Pettitt

Fracture Network Engineering (FNE) is the engineering design of rock mass disturbance through the use of advanced techniques to model fractured rock masses numerically, and then correlate field observations with simulated fractures generated within the models. Microseismic (MS) monitoring is a standard tool for evaluating the geometry and evolution of the fracture network induced during a hydraulic treatment, principally by source locating MS hypocenters and visualizing these with respect to the treatment volume and infrastructure (Figure 1). The integrated use of Synthetic Rock Mass (SRM) modeling of the hydrofracturing with Enhanced Microseismic Analysis (EMA) provides a feedback loop in which the SRM is constrained by the information provided by the MS data and the in-situ behavior of the fracture network is better understood, leading to informed decisions on future field operations. This paper discusses the technologies used in FNE and some of the developmental challenges we face in order to provide a more efficient and robust application of the approach. One particular challenge is the accurate representation of MS source data within the models (locations, magnitudes and mechanisms), and the ability to correlate these data effectively with field results. This challenge is being overcome by the development of more efficient numerical modeling algorithms, providing simulations in three dimensions at reservoir scale, combined with validation of these algorithms using laboratory and field case studies.

Jun 1, 2012

By Mohamed Hamoud

Relative permeability curves represent an important element of special core analyses, but challenges remain in obtaining valid results and the extrapolation and use of existing permeability results is problematic. This situation is compounded in cases where core properties are altered from changes in effective stress. During the production of oil and gas reservoirs, the effective stress changes as a result of changes in pore pressure or temperature. Therefore, changes in rock properties such as porosity and permeability will occur and will consequently influence reservoir performance calculations. An experimental program was carried out to understand the behavior of dense, reconstituted sand samples at low effective stress states. Drained triaxial compression tests were set up to study the effect of stress alteration on relative permeability data through porous media. Initially, a saturation stage was performed followed by isotropic loading to bring the core sample to reservoir fluid and stress conditions. Then, while axial (deviatoric stress) loading was applied to the sample, a series of absolute and relative permeability tests was performed under steady state flow at different levels of strain during pre and post peak shearing stress. The stress-strain behavior results from the study revealed that the absolute permeability decreased during isotropic loading (compaction) and increased during shearing (dilation). However, the shape of relative permeability curve changed from contraction to dilation and showed higher initial kro and the lower initial kro respectively. The goal of this project was to understand the behavior of both absolute and relative permeability at different stresses. It is expected that geomechanical processes will be adopted for special core analyses to support the current integration of geomechanics into reservoir engineering. This will be required if, during SAGD, conventional and/or unconventional reservoir simulation, reservoir engineers implement a proper set of relative permeability data at each in situ state of stress, instead of constant curves throughout the field?s life.

Jun 1, 2012

By S. Balideh

A good understanding of the behavior of broken rock helps to make more informed decisions regarding slope stability and waste dump issues in mines. The cohesion coefficient of rock has little influence on the strength of broken rock; therefore, a criterion based on friction and cohesion of intact materials does not provide an accurate analysis for the strength of broken rock. This paper illustrates the effective strength of broken rock through a brief literature review about failure of intact rock and rock masses. The complete stress- strain curve and some important points of residual strength of rock are summarized in this paper.

Jun 1, 2012

By Navid Bahrani

The grain-based model (GBM) has resolved the difficulty of matching both the directtension (sc) and unconfined-compressive (st strengths of typical brittle rock by providing a synthetic material that mimics deformable, breakable polygonal grains cemented along their adjoining sides with deformable, breakable interfaces modeled as smooth-joint contacts. In this paper, a GBM utilizing a parallel-bonded base material is calibrated to match both st and sc of Lac du Bonnet granite and then used to simulate a Brazilian test. The Brazilian strength (sB) is slightly overestimated and the failure mode includes ductile post-peak response that is not representative of what occurs in laboratory Brazilian tests. By reducing the grain tensile strength, grain failure occurs during the Brazilian test and a more realistic failure mode that includes brittle post-peak response is obtained. However, the unconfinedcompressive strength of the modified material is too low. Sensitivity analysis on the grain and grainboundary strength properties demonstrates that the grain tensile strength controls the sB/sc ratio. Because the grain strength is controlled by the parallel-bond properties and the lowest achievable st/sc ratio in a parallel-bonded material is ~0.25, it is concluded thatsc, st and sB for a typical brittle rock cannot be matched with a GBM that utilizes a parallel-bonded base material.

Jun 1, 2012

By Gautier Njiekak

The current study provides an experimental approach that aims at documenting the effects of different phases of CO2 (gas, liquid, supercritical fluid) on rock seismic responses. This has implications on obtaining quantitative estimates of fluid saturations for geological storage. The ultrasonic pulse transmission technique is the underlying methodology used in this work. Unlike previous rock physics studies involving CO2 as pore fluid, care is taken to fully saturate the studied samples and also to separate pore fluid effects from pore pressure build-up effects during the experiments. This allows to record seismic effects that are solely reflective of changes in the pore fluid properties. Measurements were done on a CO2 saturated Berea sandstone. Laboratory results show P- and S- waveform variations (velocity, amplitude) in response to CO2's varying phase state. Overall, CO2 phase changes (gas to liquid and gas to supercritical fluid) are marked by a drop in velocities, with the drop being smaller and less immediate over the gas-supercritical fluid transition than over the gas-liquid transition.

Jun 1, 2012

By Francois Charette

Corrosion is a significant factor affecting the capacity and working life of ground support. The paper explores a relatively simple method to predict the load capacity of a support system using a Monte Carlo stochastic approach. The analysis uses published field data on existing ground support to evaluate the applicability of the method and its predictive performance. Discussion on corrosion rates found in the field and predicted from laboratory measurement is provided.

Jun 1, 2012

By Azadeh Riahi

The main goal of this paper is to investigate the effects of permeability anisotropy in Large Open Pit (LOP) mines for problems that exhibit homogeneous anisotropic characteristics. The most important effect of the groundwater flow in open pit mines is on the stability of the pit slopes. This study will investigate the effects of having increased permeability in both horizontal and vertical directions on three aspects of LOP behavior, location of phreatic surface, pore pressure distribution and stability of slopes. The hydromechanical coupling is simulated using a one-way coupling approach. The Shear Strength Reduction method (SSR) is applied to assess the stability of slopes. The studies performed in this paper suggest that anisotropy in permeability leads to a significantly different regime of pore pressure dissipation depending on two aspects, ratio of anisotropy and the relative magnitude of permeability to simulated time. As a result, variation in safety factors and the location and depth of the failure surface must be expected.

Jun 1, 2012

By Renato Macciotta

The argillaceous formation known as the Opalinus Clay is host to the Mont Terri Underground Research Laboratory. The URL was developed to investigate the hydrogeological, geochemical and geotechnical properties of the formation. Underground excavations at the site are characterized by non-elastic and time-dependent deformations. Excavations trending perpendicular to the strike of the bedding presented relatively minor deformations, however, tunnels parallel to the strike of the bedding showed much larger non-elastic deformations. This paper presents an analysis of the measured response of a tunnel trending parallel to the strike of the bedding in Opalinus Clay. The measured convergence is compared to tunnel experience in weak rocks to assess the tunnel stability.

Jun 1, 2012

By Kristine Haug

Although most in situ operations in Alberta do not encounter significant problems during resource extraction, there have been a few documented occurrences of oil or steam releases as a result of compromised caprock seals. The objective of the Oil Sands Caprock Integrity Project (OSCRIP) is to characterize the units above and below bitumen zones to determine which geological factors affect the quality of the caprock seal. In addition, geomechanical properties and stress measurements of the bitumen reservoir and overlying units will be examined to provide a better picture of the containment conditions in which in situ operations occur. The information from this study will aid regulatory staff and industry in determining the zone-specific geological and geomechanical factors that could affect the design, operation, and feasibility of in situ recovery projects. The OSCRIP study area (Townships 87 to 99 and Ranges 1 to 13, West 4th Meridian) encompasses the surface mineable area, the northeast extent of the Wabiskaw-McMurray bitumen deposit and the Clearwater Formation, as well as a wide band of the Devonian salt-dissolution zone. Regional geological mapping of the Cretaceous strata focuses on the spatial extent and thickness of the caprock seal. Mapping of the underlying Devonian strata is targeted at better understanding the potential effects of salt dissolution and collapse of overlying strata on the integrity of the caprock seal. Examination of LiDAR (Light Detection And Ranging) bare earth DEM in the study area has revealed localized depressions that may represent sinkholes rooted in subsurface collapse features in Devonian carbonate and evaporite strata. Mapping the location and understanding the timing of salt dissolution-related features and deformation in the Devonian strata will help address questions about the possible effect that geological structures, such as sinkholes, have on the integrity of the caprock, as well as the potential for localized variations in the stress regime of the bitumen zone and overlying units. The stress regime in units above 350 m is such that any induced fractures are believed to be horizontal, and therefore, there would appear to be no mechanisms that would allow these fractures to reach the surface. However, recent incidents have shown that vertical pathways can occur, linking the subsurface to the surface. The final phase of the project will evaluate the changes that may occur in the stress regime under thermal operating conditions in the study area.

Jun 1, 2012

By Bruce Macdonald

Operational reservoir management activities such as cyclical steam, water or CO2 injections frequently cause surface micro-deformation. InSAR monitoring from space based sensors of this surface deformation has become a valuable approach to confirm and monitor extraction and injection activities. Sequestration of CO2 into deep saline aquifers is discussed as an important long term solution to slow future climate change. This method requires reliable continuous monitoring to assure that captured CO2 is sequestered in the substrate for ?all time? (typical minimum requirement by law is 20-30 years). We have shown that combining flow and deformation modeling with spaceborne synthetic aperture RADAR interferometry (InSAR) measurement of ground uplift caused by the underground release and spreading of the CO2 can be forged into a powerful tool to monitor sequestration. We used a novel InSAR approach to resolve the surface motion over the In Salah, Algeria CO2 sequestration site. Our approach exploits all information contained in an interferometric stack by exploiting both persistent distributed as well as persistent point scatterers. Displacement time series are obtained with a network inversion approach applied after topographic error and preliminary atmospheric correction of the network (dual-scale technique). After inverting the original network data to common master we then use an advanced spatio-temporal filtering scheme to remove the remaining atmospheric phase components down to scales of 1 km or less. By comparing final displacement series from multiple interleaving InSAR stacks we show that the accuracy of our refined solution to detect uplift outside of significant temporal data gaps in the InSAR stack is significantly better than 1 mm per year. We further show that accuracy and robustness to data gaps are both significantly improved by co-inverting multiple datasets with the same and different sensors. This co-inversion method also allows the seamless continuation of the monitoring from one sensor to another.

Jun 1, 2012

By Erik C. Westman

Knowledge of the orientation and relative magnitude of in-situ stresses at the mine scale is important for safe and efficient production from an underground mine. With this knowledge, stope location and scheduling can be optimized. Passive seismic tomography was used for imaging horizontal stress distribution at the mine scale. Microseismic measurements from the Creighton mine were used as input for the tomography calculations. The double-difference tomography algorithm was used, which provides the most precise results. The results clearly show a low velocity zone in the previously mined areas bounded by two high-velocity lobes located adjacent to the mined openings at locations consistent with those expected due to the horizontal stresses at the mine. The overarching goal of the research is to provide a tool that can augment both numerical modeling results and underground geotechnical measurements to allow the mine operator to produce in the safest and most efficient manner possible.

Jun 1, 2012

By Mohammad. M. Khajeh

Rock mechanical properties are important parameters in geomechanical simulation. Considering stochastic heterogeneous models rather than conventional homogeneous layer cake models leads to more realistic simulation results and uncertainty assessment. These heterogeneous models are usually fine scale models generated by geostatistical techniques. Although considering high resolution property models is desirable, geomechanical simulation with a fine scale model is not computationally feasible for full?field, coupled geomechanical-flow simulation. Upscaling of property models is necessary to move from fine scale/high resolution property models to coarser scale/low resolution models appropriate for simulation. Local numerical upscaling was proposed by Khajeh et al. (2012) to describe the macroscopic elastic behavior of complex heterogeneous media and is expanded upon in this work. This numerical technique is not restricted to specific geologies as typical analytical techniques are. Moreover, transversely isotropic deformation is considered rather than the usual assumption of isotropic deformation. To compare the accuracy of the upscaling methodology, in this study shear strain is considered as the geomechanical response of a synthetic facies model and the response of the coarse scale model is compared to other conventional analytical averaging techniques. It is found that the numerical upscaling technique better approximates the geomechanical response of the fine scale model.

Jun 1, 2012

By A. K. Alzo?ubi

Block toppling in natural rock slopes has been under investigation since the 1970?s. De Freitas and Watters (1973) introduced the term ?toppling? to describe the movements of rock slopes by rotation in steeply dipping beds. The Lornex pit case history from British Columbia (Canada) was chosen to study the effect of mining on the stability of toppling rock slopes. A discrete element approach that allows fracturing of intact material was applied successfully to model the composite rock slide-rock topple movement. Three main movement modes were observed in the field and the numerical model; the sliding of the top of the slope exposed normal scarps, the toppling at the middle part of the slope exposed obsequent scarps facing opposite to the slope face, and dilation of the rock caused by fracturing at the toe formed loose material. The development of the rupture surface at the toe of the slope was traced and did not evolve into a continuous rupture surface throughout the slope both in the field and the numerical simulation. Good agreement was achieved between the deformations recorded in the field and the numerical simulation.

Jun 1, 2012