Search Documents

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

By Allan R. Punkkinen

Safe extraction of the valued ore resource in the deep levels of VALE Creighton Mine has progressed in a pillar-less fashion using vertical retreat, slot-slash and uppers retreat methods since the cessation of ?cut and fill? mining in the late 1980?s. Sequencing of the panel extraction and filling has been the key to successful exploitation of high stress orebodies at depth. In late 2010, discussions between the Operating Department and Creighton Ground Control revealed an opportunity to evaluate a trial primarysecondary sequence in the west end of the 400 Orebody at the 2340m horizon. To manage the trial primarysecondary project, a high level risk analysis and management of change was carried out to mitigate potential danger. As part of the requirements, an instrumentation project was conducted in the hanging wall and footwall bottomsills of the secondary panel and within the secondary panel itself. The extraction sequence was modeled with Map3D® software, confirming high stress levels with the advance of the primary panel. Enhanced ground support was installed based on modelling results. A key objective of the instrumentation project was to verify support system response to mining and to obtain knowledge in regards to the behavior of the stope. This paper describes the first time that this type of new innovative instrumented rockbolt technology had been used in a hard rock mining environment. Understanding of the physical deformation mechanisms of ground support in response to high stress conditions is the basis for an effective ground control strategy. The data played an important role in subsequent enhanced ground support system design at the mine and importantly, established the fact that the secondary pillar had not totally failed. Subsequent cavity monitor results helped establish that fact.

Jun 1, 2012

By Cheowchan Leelasukseree

Mae Moh mine, operated by Electricity Generation Authority of Thailand (EGAT), is the largest open pit mine in Southeast Asia. Its production is approximately 15 million lignite tons and 90 million m3 of waste, annually. The bottom of the final pit limit is designed to be approximately 500 m in depth. Presently, the pit depth is almost 250 m below original surface. The complexity of geological structures, such as faults, has caused a number of geotechnical difficulties. In a challenging lowwall area, called Area 4.1, an inclined weak bedding plane, G1, will be daylighted along the lowwall as following by mine plans. The dimensions of Area 4.1 are nearly 300 m wide and 90 m high. The friction and dip angles of the G1 interface are reportedly 12° and 15°, respectively. A three million m3 claystone block, resting above G1, would potentially slide down, and possibly cause a massive failure. EGAT engineers are acutely aware of the potential gigantic sliding failure and are modifying the mine plan in Area 4.1 to improve mine safety and minimize the displacement of the potential sliding block, reduce the possibility of failure and lower the additional mining cost as much as possible. The cut and fill mining technique was selected and studied to achieve those goals for mining in Area 4.1. With this mining technique, G1 will be partially daylighted. Therefore, an arch formation is induced with a sliding block above the mining area, which could perform as a temporary support undercutting a certain span. Initial studies of the behavior of the partially undercut slope are presented in this paper, along with a review of the similar failure history in Mae Moh mine and a 3-dimensional numerical model of the mine. The failure history review provides an estimation of undercut span without failure. The three-dimensional numerical model results, using 3DEC, gave an initial approximation of the displacement of the sliding block in Area 4.1. The next exciting study of the behavior of partially undercut slope is a trial undercut slope field experiment in Area 4.1. Currently, a trial undercut slope field experiment in Area 4.1 is under construction. In the field experiment, a number of monitoring programs will be installed in Area 4.1 to monitor and obtain information to allow a better understanding of the behavior.

Jun 1, 2012

By Melanie Grob

Two commonly estimated fractal dimensions, called b and D values, represent statistical characteristics in the distribution of magnitude sizes (b) and spatial hypocenter locations (D) of microseismic events, respectively. We postulate that the values of b and D are related to specific stress regimes. Through different case studies of fluid injection into an oil/gas reservoir or geothermal field, we infer that the measured temporal variation in fractal dimension b is most likely due to significant changes in the local stress regime. The fracture dimension D indicates the spatial clustering of events. The distribution of these events can vary with time because of stress perturbations that (re)activate fractures whose direction favours their movement under a specific stress regime. Thus a combined statistical study of these two dimensions can help infer the variations of local stresses in a reservoir.

Jun 1, 2012

By Jaime Meléndez Martínez

Recent interest in unconventional reservoirs motivates our work in laboratory measurements of seismic anisotropy. Seismic anisotropy is the variation in speed of a wave as a function of its direction of propagation. Analyzing anisotropy in unconventional reservoirs is important since anisotropy leads, for example, to differential stresses upon loading and could affect hydraulic stimulation (Reinicke et al., 2010, Zimmermann and Reinicke, 2009); in this sense, laboratory measurements are an important tool to study seismic anisotropy at small scales which can aid in the characterization of larger formations. The purpose of this work is to describe the laboratory methodology used to simultaneously measure ultrasonic P and S-waves in three different directions from a single core sample. Elastic anisotropy of a Proterozoic sedimentary rock from south-west Alberta is investigated. In metamorphic rocks seismic anisotropy is caused by preferred mineralogical alignment or foliation. Assuming a transversally isotropic medium, arrays of compressional and shear piezoelectric ceramic transducers were mounted on a sample trimmed from a core to measure travel times of ultrasonic P-waves and S-waves as a function of confining pressure at directions perpendicular, parallel, and oblique to the plane of foliation (bedding). Results show strong anisotropy agreeing with the assumption of a transversely isotropic medium. Dependence of travel times and hysteresis effects can be observed when pressurizing and depressurizing.

Jun 1, 2012