Search Documents

By Shahé Shnorhokian

A study was conducted on the seismic activities and associated rockburst events that took place at an underground Canadian metal mine over a period of five years. In the first stage, a back analysis of the microseismic database led to the identification of clusters of major events. The locations of these clusters were identified with respect to the local geology and the planned mining sequence. In the second stage, a numerical model was constructed using the finite difference code FLAC3D to simulate the in-situ stress, mine geometry, local geology, and stope sequence. Mining-induced stress fields were compared against two rockburst events that took place during this period, as well as against the various geologic units present. It was shown that by combining microseismic analysis and numerical modeling, it was possible to provide a basis for evaluating the influence of the in-situ stress regimes and stope sequencing on mining-induced seismic response.

May 1, 2011

By S. Stefanizzi

In the present paper, numerical simulations based on a finite/discrete element approach are used to study deformation and failure process during standard rock mechanics. FEM/DEM modeling has been used as available with the ELFEN code, which explicitly model the transition from continuum to discontinuum. The transition is simulated through a crack nucleation and propagation that obeys to the Griffith?s failure criterion. If the fracture criterion within the intact rock (represented by FEM) is met, a new crack (represented by DEM) is initiated. Remeshing allows the fracture process through the FEM mesh to be tracked and visualized, thus contact properties can be assigned to pre-existing fractures and newly generated fractures. In this study, compressive tests of rock samples are simulated. The applicability of ELFEN to properly and efficiently simulate typical laboratory rock mechanics tests, with reference to the rock fracture process, is described in detail. In addition, an indirect validation of the extensional strain criterion proposed by Stacey (1981) is given.

May 1, 2009

By A. P. E. Dirige

Extensive field study has been performed to assess the capabilities of thin sprayon liners (TSLs) and conventional spray-on support systems for preventing rock and support material damage that often results due to mining-induced rockbursts. Support performance was studied using field scale explosive detonation trials to simulate dynamic failure effects that are known to develop during typical rockburst events. Multiple seismic and high speed photographic monitoring techniques were used to provide detailed information concerning rock motion, surface fracturing, ejected fragment motion, and support liner survivability characteristics. Results of this study have demonstrated that TSL?s and variant layer combinations may be as effective as or better than conventional support materials for mitigating rockburst or like damage in highly stressed mine environments. In this study, field results have been verified using numerical modeling procedures to better understand the support behaviour of TSLs when subjected to highly stressed mine environments and mining-induced rockbursts, a factor that is essential when designing any rock support system. A series of numerical modeling assessments were conducted using FLAC3D modeling of a half circular TSL-lined tunnel, influenced by anisotropic stresses and simulated mining-induced rockbursts, to investigate the support potential of TSLs under prototype conditions. The effect of a mine rockburst was simulated using the three-dimensional dynamic analysis option of FLAC3D. Results of TSL-lined tunnels in highly stressed and mine rockburst conditions indicated exceptional effectiveness of TSLs for suppressing rock deterioration resulting from stresses and rockbursts. The rock support potential of TSLs on tunnel surfaces was compared with that of thin shotcrete linings, and results indicated that thin spray-on lining products currently available may be equally as effective in support as shotcrete materials.

May 1, 2009

By R. Pakalnis, D. Tenney, B. Lang

"This paper describes three case histories whereby numerical methods were employed as a tool in the identification of potential instability. Simple two dimensional (PCBEM) and pseudo three-dimensional (BESOL, MINTAB) codes were employed to identify zones of induced stress which were subsequently related to extensometer, stress meter and observational data as the stope was excavated. The following case histories will be discussed. MacLellan Mine - sill pillar stability; Ruttan Mine-Hanging Wall Stability; and Detour Lake Mine-narrow vein mining.The MacLellan case history showed how modelling and instrumentation enabled mining to be conducted safely and under controlled conditions. MINTAB and PCBEM were employed to predict sill pillar failure. The sill pillar failure was associated, with an unloading of stress which was recorded by instrumentation (stress meter) and a decrease in the over-all rock mass strength. The yielded pillar, however, maintained its over-all configuration enabling access through the pillar.The Ruttan case history is an example of hangingwall and footwall design whereby the existence of adverse structure coupled with a zone of relaxation results in dilution. BESOL and PCBEM were employed to model a test stope and to verify the observations and measurements through instrumentation . Empirical methods of dilution estimates are employed to determine the level of walls lough with respect to the rock mass rating, stope dimensions and rate of excavation.The Detour Lake case history in an example of open stoping a narrow steeply dipping deposit. Stope dimensions are 135m in strike length, 47 m in vertical height with an average width (F/W to H/W) of3.6 m. Stress meters coupled with visual observations are related to modelled estimates of stress and dilution."

Jan 1, 1991

By O. K. Mahabadi

The scope of this study is to simulate the behaviour of a layered rock sample under standard laboratory Brazilian Disc test using an innovative combined finite-discrete element method (FEM/DEM) research code. The influence of layering and the direction of loading on the samples behaviour are studied through various layering and loading configurations. A parametric study including the effect of loading rate on the tensile strength of samples has been carried out. This paper demonstrates the suitability of the numerical approach to explicitly model rock deformation and failure.

May 1, 2009

By Mohamed Sassi, I. Kolo, R. L. Sousa, R. K. Abu Al-Rub, Manhal Sirat

"Fractures induced by folds are pivotal in hydrocarbon exploration, ground water transport and harnessing of geothermal energy. This is because of the need to predict and understand fracture propagation in reservoirs which have folded rock formations. Despite its prominence, there is a lack of reliable numerical models for simulating fractures resulting from rock folding. This is partially due to the difficulty involved in fracture modelling. As a contribution to fold-fracture modelling, this work applies a coupled plasticity-damage model to rock fracturing and anticlinal folds. Parametric studies on a single folded layer give insight into the fold formation and behaviour. The anisotropic continuum damage model which considers both tension and compression is formulated using the power damage evolution law. For ease in formulation, strain equivalence hypothesis is adopted whereby the strain is the same for both damaged and undamaged configurations. Lubliner plasticity yield criterion is adopted for plastic deformation. The model is coded in Abaqus user subroutine UMAT and is applicable to quasi-brittle materials.INTRODUCTIONThe understanding of fractures in the earth’s crust helps in the prediction, evaluation and characterization of fractured reservoirs, whether they are petroleum, geothermal or underground water reservoirs. This is partly because reservoir permeability and stability are dependent on fracture properties (Nelson, 2001). To achieve this, various approaches have been developed to analyze reservoir rock fractures. These could be grouped into (a) discrete approaches – where the rock mass is represented by a finite number of well-defined components and (b) continuum approaches – where the mathematical assumption of an infinitesimal element is made (Jing, 2003). Ivanova (Ivanova, 1998) comprehensively studied fracture systems in nature and highlighted folds as one of the major geologic settings that produce fractures. Naturally occurring folds in reservoirs sometimes compound fracture analysis. It might not be clear whether the fractures formed before, during or after rock folding. Fractures formed during the folding process are expected to show stress orientations related to folding and are thus of primary concern in rock folding simulation. It could be argued that there are no numerical models that as"

Jan 1, 2015

By D. J. Reddish

The paper describes a case study where numerical modelling has been applied to predict potential risks associated with inflows of water in active mine workings from an overlying sandstone aquifer. The objective of the modelling was to predict the development of strata fractures and strains as the longwall panels were extracted and investigate their relationship to an sandstone aquifer. The numerical modelling study also included the investigation of the effect of mining on a fault that existed close to the start position of one of the panels. The potential effects of mining on the fault were undertaken by parametric studies adjusting properties for the fault and the relative position of extractions, which allowed the mine planner to position the panel layout accordingly.

May 1, 2003

By C. V. Deutsch, O. Babak, J. G. Manchuk

High levels of geological heterogeneity occurring over relatively short distances pose a challenge to the unconventional open-pit techniques used in oil sands mining. There is strong motivation to mine selectively and minimize dilution and loss; however, the large equipment and high mining rate used to meet production targets tend to reduce selectivity. A numerical model is developed to account for geological complexity, mining equipment, regulatory criteria, and other operational parameters. This framework is used for more accurate forecasts, optimization of equipment, and other operational decisions. An example demonstrates the approach in light of Alberta Energy Regulator criteria.

Jan 1, 2014

By T. Saksala

In this paper, the embedded discontinuity approach is applied to finite element modeling of rock fracture. A rate-dependent constitutive model based on the embedded displacement discontinuity theory is developed to describe the mode I, mode II and mixed mode fracture of rock in tension and compression. The bulk material is described as linear elastic until reaching the elastic limit. Beyond this limit, a ratedependent exponential softening law governs the evolution of the embedded displacement jump. The present approach incorporates the rock heterogeneity by random description of the mineral texture of rock. Moreover, the inherent initial microcrack populations of natural rocks are accounted for as randomly oriented embedded discontinuities. Numerical examples demonstrate the model behavior in uniaxial compression and tension. The effect of loading rate and confining pressure is tested as well in 2D numerical simulations. These simulations show that the model captures the main features of rock in confined compression and uniaxial tension. The developed method has the computational efficiency of continuum plasticity models. However, it has an important advantage of accounting for the orientation of introduced microcracks.

Jan 1, 2015

By F. Liu, H. A. Meier, D. M. Valiveti, P. A. Gordon

Hydraulic fracture is a frequently used technique in the Oil and Gas industry to create high permeability pathways for extraction of hydrocarbons. Numerical simulation of this method poses challenges such as representation of 3D fractures, propagation of randomly oriented fractures in 3D media, enabling fluid leak-off from fracture to poro-elastic medium. We have developed a numerical model using the Generalized Finite Element Method (GFEM) to simulate fluid driven fracture in a poro-elastic medium that addresses these challenges. In the GFEM framework, special enrichment functions are used to augment the solution space in the displacement and pore-pressure fields in order to accurately capture the response in the vicinity of the fracture. Adaptive refinement of the finite element mesh is also employed to facilitate accurate and efficient solution of the model. Our model fully couples fluid flow in the fracture, represented with the lubrication approximation, to the geo-mechanical behavior of the bulk poroelastic rock, comprised of an elastic rock skeleton saturated with fluid. Fracturing fluid exerts traction on the fracture surface and is also allowed to leak-off into the poroelastic medium using a transfer function approach. Stress intensity factors (SIFs) are evaluated based on the Contour Integral Method, accounting for the fluid pressure on the fracture surface. The SIFs are used to govern a full 3-D evolution of the fracture. Verification of the model is performed for known behaviors for the poroelastic response and fluid flow in the fracture.

Jan 1, 2015

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 T. Miura, Y. Mori, H. Aoki

In a flash smelting furnace, easily combustible concentrate particles excessively oxidized and melt. The excessively oxidized molten particles then collide with less oxidized solid ones to create uniform oxidized molten particles. This collision phenomenon plays a very important role in achieving high smelting performance. The authors developed a mathematical model describing above combustion phenomena. This model incorporated fluid flow, heat and mass transfer, chemical reactions and collisions of concentrate particles. Both gas flow and particle motion were calculated using the Eulerian method. Copper concentrate was assumed to consist of chalcopyrite(CuFeS2). Considered in the concentrate reactions were the decomposition of CuFeS2, and the oxidation of the resulting sulfur(S) and pyrrhotite(FeS) to produce magnetite(Fe30 4) and sulfur dioxide(S02). Moreover a reaction due to concentrate particle collisions, in which Fe304 in the excessively oxidized molten particle reacts with FeS in the less oxidized particle, was also considered. The particle collisions were assumed to occur according to a collision probability rule obtained by computational experiments. Particle growth in diameter was described as a change in the volume fractions of particle phases. When applied to a commercial flash smelting furnace, this model reproduced the particle growth and component transfer as observed in the actual furnace.

Jan 1, 2007

By S. H. Seyedein

A mathematical model was developed to study fluid flow and solidification in a single-roll casting system. Fluid flow coupled with energy equations were solved numerically using a finite volume method. A body-fitted coordinate transformation technique was adopted to handle the irregularly shaped domain. In order to couple the continuity and momentum equations, a SIMPLEC algorithm was used. A fixed-grid enthalpy method was considered to model the solidification in the system. The turbulence effect on flow field was modeled using a low-Re k-e model The model developed in this work was further run under various conditions to optimize the operating parameters of the process.

Jan 1, 2004

By Y. Jin, H. You, H. Pang, B. Lin

Successful application of steam assisted gravity drainage (SAGD) on the recovery of land facies ultra-heavy oil in Xinjiang oilfield depends largely on shortening the preheating period with the aid of the micro-fracturing technology, which is conducted through controlled water injection to a pair of horizontal wells to create homogenously distributed microcracks in the zone between the two wells. The numerical simulation of this mechanical process is yet to be investigated. Up to date, there has been research work on rock mechanics during the SAGD production period of marine deposited ultra-heavy oilsands, but the counterpart research on land facies ultra-heavy oilsands of Xinjiang oilfield is left undone. In this regard, this study is attempted to perform numerical simulation of the fracturing process on land facies unconsolidated oilsands. With this objective, the study first investigates geological and physical identification of the field collected cores, then carries out mechanical analysis on laboratory triaixal test results, based on which numerical simulation of the fractured zone is performed. The simulation is applied through a finite element platform adopting extended Drucker-Prager constitutive model, with the parameters determined by experimental work. It is verified that hydraulic fracturing of unconsolidated oilsands generates a microcrack feature in terms of a dilative zone instead of linear apertures. The extent to which the fractured zone propagates is visually demonstrated based on the magnitude distribution of equivalent plastic strain or volumetric dilation. Evaluation of the connection of injection and production wells can possibly be determined by the expansion of the saturated region around each well. Finally, the implications on the efficiency of hydraulic fracturing in stimulation of the pay zone are investigated.

Jan 1, 2015

By K. Hashiba, Y. Z. Liang, K. Fukui, M. Koizumi

"Efficiency and penetration rate during percussive rock drilling significantly affect the speed of tunneling and excavation. Many researchers and engineers have continuously engaged in the advancement of rock drilling performance. Development and modification of the rock drills based on test data are timeconsuming and costly. Hence it is essential to make effective use of numerical simulation. However, most of the numerical models proposed in previous studies were too simple in comparison with the actual percussive rock drilling in which many factors are involved. The objective of our research project is to construct a new realistic numerical model of continuous percussive drilling, which is based on the previous Okubo-Nishimatsu’s model. In this study, the parts of the previous model involved in stress wave propagation and attenuation were modified. The two models were prepared for HD712 hydraulic rock drill mounted on crawler drills and HD210 hydraulic rock drill mounted on drilling jumbos. Stress wave propagation was computed with the discretized one-dimensional wave equation. The calculated results were compared to the stress waves measured on a rod during percussive drilling. In both models for HD712 and HD210, the measured stress waves can be reproduced with the one-dimensional model of the same acoustic impedance as the actual piston, shank rod and rod combined with the CI+spring model as a rod joint. The calculated results indicate that the spring constant in the model of a rod joint needs to be set appropriately depending on drilling conditions. INTRODUCTIONEfficiency and penetration rate during percussive rock drilling significantly affect the speed of tunneling and excavation. Many researchers and engineers have continuously engaged in the advancement of rock drilling performance. Components related to the percussive drilling are a drill body, a piston, a shank rod, rods, rod joints, a bit and rock. The piston in the drill body is reciprocated by hydraulic pressure and collides with the shank rod a few thousand times per minute. Stress waves generated by the collision propagate and attenuate in the rods and rod joints, and then reach the bit. The bit penetrates and crushes rock. During percussive drilling, the components mutually affect each other and self-induced oscillation occurs. The number of blows and the rate of penetration are determined by the practical drilling conditions such as hydraulic pressure, thrust force, number of rods and mechanical properties of rock. Development and modification of the rock drills based on test data are time-consuming and costly. Hence it is essential to make effective use of numerical simulation. However, most of the numerical models proposed in previous studies were too simple in comparison with the actual percussive rock drilling in which many factors are involved."

Jan 1, 2015

By T. G. Sitharam, Resmi Sebastian

"Wave propagation studies across jointed rocks have significance in the fields of rock mechanics and geophysics. The passage of a seismic wave induces strains in rocks and often, these strains are of low magnitudes except at the vicinity of source. These low magnitude strains are generally in the elastic range and hence elastic constants necessary for the design and analysis of strength of rock mass, can be derived from wave velocities. Measurement of travel time velocities and measurement of resonance frequencies are the two methods usually adopted for obtaining wave propagation velocities.Resonant column apparatus measures wave velocities from the resonant frequencies obtained during testing. Depending on the type vibration applied, torsional or flexural, shear and compression wave velocities are secured. By varying the amplitude of vibrations applied, various strain levels can be introduced and corresponding resonant frequencies and wave velocities can be calculated. The resonant column experiments were performed on a model material, hardened Plaster of Paris (POP) samples, in laboratory. The paper explains the preparation of samples and the working mechanism of resonant column apparatus along with the estimation of wave velocity from resonant frequency. The numerical simulation of resonant column apparatus using distinct element method and its validation with the experimental results are presented in the paper.A distinct element code is usually used to simulate the movement of materials with intersecting discontinuities. Finite rotations and displacements of individual blocks representing the material are permitted and new contacts between the blocks are identified as the calculation progresses. Three dimensional Distinct Element Code (3DEC), the software intended for the numerical simulation of rock engineering projects, can well simulate the behavior of jointed rocks under static and dynamic loading conditions. The numerical simulation of resonant column apparatus has been performed with 3DEC considering its efficiency in modeling jointed rocks.The numerical modeling of the resonant column apparatus with application of torsional and flexural vibrations to the sample is elaborated in the paper. The response of samples to a range of frequencies is obtained and the frequency that induces maximum displacement in the sample is measured as the resonant frequency, for any particular strain level. The displacements measured at various points of the sample during application of vibrations are compared and presented. The numerical modeling of resonant column apparatus with an intact sample (without joint) and a sample with a joint incorporated are explained at a particular strain level. The resonant frequencies obtained from laboratory experiments and numerical simulations are compared and illustrated. It has been found that resonant frequencies obtained using numerical and laboratory experiments are in good comparison and can be extended for various strain levels and sample conditions."

Jan 1, 2015

"A numerical procedure based on Smoothed Particle Hydrodynamics (SPH) is developed for rock blasting phenomena involving multiple intersecting joints. In this procedure, an elasto-plastic damage constitutive model has been applied in the framework of SPH to simulate yielding response of the surrounding rock medium. In this paper, a square shape rock medium with centrally located emulsion explosive in a borehole is numerically blasted with and without the presence of two diagonal rock joints. The model is simulated by assuming all four sides of the rock medium as free surfaces. The paper elaborates the propagation of stress wave, expansion dynamics of gaseous products, and development of fracture, fragmentation of the medium and dissipation of energy due to damage.INTRODUCTIONDamage prediction of rock mass and its mean fragmentation size are of important to design blast parameters for a given explosive. Current blasting practice solely depends on empirical relationships between fragmentation sizes to those blast design parameters while mathematical descriptions of the physical laws governing the process of rock blasting are limited. A constitutive relationship based on properties of explosives and rock mass is of special interest because of its capabilities in describing the key features of rock fracture and fragmentation, resulting in understanding of the physical process involved in rock blasting. It is well known that blast-induced damage and fragmentation of rock mass occur due to detonation induced stress wave and product gas driven fracture propagation. Nevertheless, the fundamental mechanism of the failure process and the behavior of rock characteristics under these two different loading are not well understood in many cases, which are closely related to growth of fractures and interaction of cracks under different loading conditions.Continuum damage mechanics is considered to be appropriate to describe the failure process under such high-amplitude stress wave. Based on this approach, several theoretical models have been developed to analyze the dynamic damage of brittle material such as rock and concrete with microstructure (Chen, 1999; Taylor et al., 1986; Wu et al., 2004; Yazdchi et al., 1996). These models are mostly based on the theory of damage mechanics or fracture mechanics with some statistical treatment to incorporate random distribution of microcracks. In the application of fracture mechanics theory, sometimes it is difficult to estimate the micromechanical parameters such as crack density and fracture toughness in different modes. In this study, damage in rock material is evaluated according to the earlier model developed by Grady and Kipp (1980), which is a fracture model coupled with a material description for stress wave propagation to predict quantitatively fracture and fragmentation under explosive loading conditions."

Jan 1, 2015

By B. Panton, J. Carvalho, E. T. Brown

"A major consideration in the design of high capacity tiedown anchors for dam, bridge and tower foundations is the tensile resistance of the rock mass to pullout, typically as a result of overturning moments or hydrostatic uplift. Rock mass pullout capacity for installed anchors is developed from the tensile strength and fracture propagation properties of intact rock, the orientation and physical properties of the discontinuities and anchor confinement at depth. The typically assumed, but generally conservative, design approach is to calculate anchor pullout capacity using the dead weight of a uniformly shaped inverted “cone” with an assumed initiation point and breakout angle. In some cases, some tensile or shear strength is considered across the area of the assumed pullout cone, however, these strength estimates are based primarily on based on limited historic testing data, empirical relationships for rock mass strength developed for other applications, and engineering judgement.As an alternative to the current foundation anchor design method, this paper presents a new approach to the rock cone pullout problem using Discrete Fracture Networks (DFN) combined with numerical simulations. The simulations presented in the paper investigate the influence of fractures in a synthetic rock mass on ultimate anchor strength. In addition, the global influence on fracture network statistics on rock mass pullout failure is also discussed and it is contended that some of the aleatory and epistemic uncertainty inherent in the design of rock anchors can be reduced by considering fracture statistics in anchor design."

Jan 1, 2015

By Yang Shuxin, Guo Qiliang, He Ping, Lu Lizhi, Yao Rui

"We combined stepwise regression and constrained regression statistical method to obtain the best boundary conditions for the calculation of in-situ stress field in Wushaoling tunnel using finite element method. Our numerical results reveal that: the main geologic factors contributing to the initial background stress are gravity and tectonics; the principal stress values would drop in fault zones; the angle between the maximal horizontal stress and the axis of Wushaoling tunnel is acute; Wushaoling tunnel tends to become stable when this tunnel has 2:1 ratio of the height to its width; the possibility of rock-burst exists under high stress.INTRODUCTIONStress state of engineering region is an important factor affecting the stability of project. Geological disasters in the project area are closely related to their stress environment. Wushaoling tunnel which is 20,050 meters long, lying along Lanzhou-Wuwei segment of Lan-Xin Railway in China, is the key project of Lan-Wu double-track railway. The region surrounding this tunnel is under a high stress environment. Accurate identification of the stress environment of the engineering area and the interaction between geomaterials and geotechnical engineering construction holds significant scientific and practical value for engineering design, construction, and operation.Many scholars explored and discussed the present-day tectonic stress through investigations of tectonic evolution history (Ding, & Lu, 1986), crustal deformation (Zoback, 1992), focal mechanisms (Yang et al., 2003, Xie, Cui, & Zhao, 2004), and so forth. However, the above methods would only provide the approximate direction of the principal stress in a broad region, and could not give the absolute stress values. In-situ stress measurement is one of the primary means to directly gauge the current crustal stress. However, in-situ stress measurement can only be applied in small number limited by measurement conditions."

Jan 1, 2015

By S. Kawai

Natural convective flow of the electrolyte solution plays an important role in several electrochemical systems. Copper refining electrolysis is the most important industrial process to produce high purity copper on large scale. Turbulent natural convection develops along vertical plane copper electrodes installed in an unstirred electrolyte accompanying with the electrochemical deposition and dissolution of copper during the electrolysis. It introduces the electrolyte stratification phenomenon. Thus, the electrolyte is circulated in a special way depending on each electrorefining tankhouse. The theoretical discussion is indispensable to understand the coupling phenomenon between the ionic mass transfer rate and the morphological variations of electrodeposited copper. In present study, the mass transfer rate of both Cu2+ and H+ ions as well as the momentum transfer rate is numerically analyzed accompanying with the turbulent natural convection along 1 m high Cu metal electrode surface.

Jan 1, 2007