Search Documents

By M. B. Dusseault, A. Pirayehgar

"There is concern that hydraulic fracturing (HF) near an existing critically stressed fault could trigger slip and generate a seismic event of significance. To this end, a mathematical discontinuum model representing a naturally fractured rock (NFR) mass containing a fault is created, and fault slip response to fluid injection is studied. UDECTM is used as a Discrete Element Method (DEM) approach to define the seismic moment, and also to assess the influence of injection and slip on further shear slip propagation. In most of the simulations, slip and dilation take place along a number of oriented fracture surfaces, so it is reasonable to assume that in a real case, microseimic emissions would accompany each slip “event”. The moment magnitude scale (Mw) is used to measure the size of the microseismic event in terms of the energy released. Seismic moment is a measure of the total energy (work) released during a seismic event and is used to measure its size. From another aspect, the injection energy is calculated as the total energy required for achieving the hydraulic fracturing; this is viewed as energy to open fractures against the ambient stress field and much of it is potential stored energy that can be released and create discontinuities (work of fracture), deformations (W=F·d), heat and radiated seismic energy. A small amount of the total released energy is transformed to radiated seismic energy and another small amount results in tensile fractures. Fracture energy is used to calculate the work required to create fracture surfaces. Comparing the energy release by fracturing or seismic events indicates their contribution to the total amount of released energy or total work. The amount of energy associated with the stick-slip event is small compared with the overall work done to generate the distortions (i.e. increases in aperture by forcing them open) in the rock mass, and much of the energy in slip is lost as heat from the plastic deformation, rather than as radiated seismic energy.Microseismic monitoring shows the spatial distribution and magnitude of seismicity associated with slip of bedding planes and natural and incipient fractures. Work calculations have been undertaken to look at the various sources of energy storage (elastic strain energy) or dissipation (slip), and likely the energy losses from viscous dissipation could be addressed as well, helping to clarify HF processes in NFRs."

Jan 1, 2015

By U. Caliskan, S. Miskovic

Gravity separation is a technique that depends on the difference in density of minerals for their separation. Although this method is one of the oldest separation techniques, it is still preferred in the industry because of its low capital and operating cost while it is environmental-friendly. In this paper, we present a new numerical methodology that can be used to investigate dynamics of dense slurry flows in gravity separation systems. A simple channel with a set of riffles is selected as a test platform in order to observe the effect of the riffle shape on segregation profile of particles with different densities. It is believed that the riffle shape would affect the slurry flow and possibly cause one class of particles to segregate within the riffle zone. The turbulence, as a disturbance effect within the riffle zone, causes heavy particles to displace the light, that is, the heavy particles will be trapped in riffles, while the lighter particles are carried by water to the end of the channel, which is an example to uninterrupted suspension. MP-PIC model in OpenFOAM is used to study the effect of riffle shape and slurry percentage of solids in this simple gravity separator on particle segregation profile of four discrete mineral density classes, which are Quartz/Clay (SG 2.6), Pyrite (SG 4.2), Chalcopyrite (SG 5.1), and Gold (SG 19.32). Gold particles are found to be readily captured and retained in the riffles because of their high density, while particles with the intermediate densities are more influenced by the flow. The developed model framework has been proven robust and capable of providing physical results even for the flows with a very high percentage of solids, characteristic for mineral processing operations. INTRODUCTION Gravity-based ore concentration techniques make use of differences in density of minerals within an ore to enable separation of valuables from the gangue. As some of the oldest techniques used for ore beneficiation, they are known to offer a range of benefits such as low specific energy consumption and low capital and operating costs (Burt, 1999; Falconer, 2003; Holland-Batt, 1998), however, they suffer from low selectivity, poor separation efficiency, and low capacity. Gravity separation techniques are common practices in the mining industry and are employed widely for the precious mineral concentration, iron ore recovery and coal preparation (Kroll-Rabotin, Bourgeois, & Climent, 2013; Laplante & Spiller, 2002; Williams, 1997). For example, gravity separation systems are highly applicable for precious metal enrichment operations by being included in primary grinding circuits, shaking table tailing separation, pre-concentration of leach feed and flotation concentrate enrichment for both flash and conventional flotation purposes (Burt, 1999; Falconer, 2003; Holland-Batt, 1998; Laplante & Spiller, 2002; McAlister, 1992). Gravity separation can be employed in the copper processing circuits by being integrated into the pre-flotation circuit to capture the gravity recoverable gold (GRG) and any free copper aiming to enhance the gold and copper recovery of the plant. When native copper and GRGs are present in the deposit, this integration not only facilitates plant overall performance via removal of the valuables early in the circuit, but it also reduces downstream operational expenditures due to less retention time of coarser copper particles in the system.

Jan 1, 2019

By S. Pietruszczak, G. Su, T. S. Nguyen, E. Haghighat

"In this paper, the mechanical properties of Tournemire argillite are modeled using a mathematical framework which employs the microstructure tensor approach (Pietruszczak & Mroz, 2000). The experimental data incorporates the results of a series of triaxial tests (Abdi & Evgin, 2013) conducted at different orientation of bedding planes relative to the direction of loading. First, a plasticity framework is presented which involves the specification of the conditions at failure as well as the description of the irreversible deformation process. Later, a systematic procedure is outlined for identification of material functions/parameters involved. The paper is concluded by presenting the results of numerical simulations of the triaxial tests. It is demonstrated that the model can adequately reproduce the primary effects of anisotropy in the stable deformation regime as well as the onset of fracture.INTRODUCTIONThe understanding of the mechanical behaviour of argillaceous rocks is important in many types of geotechnical projects, including petroleum extraction, carbon dioxide sequestration as well as deep geological disposal of radioactive waste. In all these cases, the primary concern is the onset and propagation of damage due to excavation, transport of pore fluids and/or elevated temperature. The Canadian Nuclear Safety Commission (CNSC), Canada’s nuclear regulator, will be responsible for the licensing of future deep geological repositories. In order to evaluate the safety of deep geological repositories, CNSC staff keep abreast of scientific developments by collaborating with national and international partners on experimental and theoretical research on the long-term performance of those engineered and natural barriers. Through that collaboration, the CNSC has access to data resulting from experiments performed at different Underground Research Laboratories (URL). Such an URL is situated in Tournemire, France, in an argillaceous rock formation, the Tournemire argillite. The Tournemire argillite is characterized by the presence of closely spaced bedding planes and exhibits a strong anisotropic behaviour, so that its strength as well as deformation properties are directionally dependent."

Jan 1, 2015

By H. Zniber El Mouhabbis

Ore dilution or stope overbreak increases the cost of production and ultimately reduces the profitability of a mining operation. The adverse economic impact of ore dilution is due to the added costs associated with the mucking, haulage, crushing, hoisting, milling, and treatment of waste rock or low grade ore in the hangingwall and footwall. One factor that could have a significant influence on stope overbreak is the undercutting of the stope into the walls. Undercutting stope walls is often necessary, in particular in narrow vein mines. Moreover, equipment size has a role in drift sizing where stope undercutting may become a must in order to fit the mobile equipment while respecting legal clearances.

May 1, 2009

By N. Manopulo, J. Maier, L. Tong, C. Becker, P. Hora

The contribution presents a new finite element approach in which the common updated Lagrangian formulation is replaced by a spatial flow formulation applying a special arbitrary Lagrangian-Eulerian for-mulation (ALE). The description separates spatially the process in a die inside zone and in a profile zone. For the description of the inside zone, where the filling process occurs, a special Euler time independent mesh is applied. By this formulation the mesh has to be generated only once, which simplifies the generation of spatially optimized meshes. For mapping of the material flow during the filling process the method uses a flow surface traction method as applied in casting simulations. A further strength of the new method is the computational efficiency which allows the modelling of extrusion process with the highest complexity with a response time of 12 hours, which is important for the industrial optimization process.

Jan 1, 2015

By N. Golchinfar

Tunnels in deep mines are subjected to high stresses. Fracture and failure of the rock mass around an opening can occur when mining-induced stresses are high. In burst-prone grounds, mining-induced seismicity can cause additional dynamic loading which can further increase the stress around the tunnels. When the dynamic stress reaches the rock mass strength, fracturing can occur and rockburst may happen (depending on loading stiffness around the failed rock mass). Rock support installed in the tunnel must be capable of dissipating dynamic energy and holding the failed rocks. Hence, it is important to estimate the depth of failure for rock support design. This paper focuses on how to model depth of failure under dynamic loading. In the modeling, two scenarios of rock failure are considered. In the first scenario, rock failure occurs under static loading and subsequent dynamic loading further increases the depth of failure. In the second scenario, no failure occurs under static loading but in the subsequent dynamic loading, failure around the tunnel is created.

Aug 1, 2013

By M. N. Vu, G. Armand, A. Noiret, D. Seyedi

Extensive experimental observations have been performed around drifts and shafts at the Meuse Haute-Marne (North-Eastern France) Underground Research Laboratory (URL) to assess the extent and pattern of the induced fracture networks in Callovo-Oxfordian claystone. Two types of fractures can be distinguished, namely; shear (mode II) and extension unloading (mode I) fractures. The main goal of the present paper is to provide insights on the extent of damaged zone around the drifts by numerical simulations. A two-step approach is considered. First, 2D elastic calculations are performed for simulating the excavation of a drift following the direction of the major in situ horizontal stress. The evolution of the stress and strain tensors around drifts are calculated by applying an elliptic convergence function on the drift wall. Two failure criteria are then verified for the induced damage. An equivalent strain and an elastic energy based damage criteria are considered. The obtained results show that the employed simplified approach provides rather good results for the extent of the damaged zone for drifts parallel to major horizontal stress. However, post-processing elastic results leads to overestimating the vertical extent of the damage zone and underestimating its lateral extent, where the in situ stress state is quasi-isotropic in a section of drift. Nevertheless, the obtained results show that an energy based damage criterion could be considered as a potential way for reproducing the extent of the damaged zones if a convergence function can be defined for the drift.

Jan 1, 2015

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