Search Documents

By Ed McConaghy, P. LaRochelle

The Ambatovy Nickel Project comprises the development of mine, hydrometallurgical processing facilities and necessary infrastructure in Madagascar, positioning Ambatovy as one of the world’s largest producers of nickel and cobalt from lateritic ore. SNC-Lavalin successfully provided comprehensive Engineering, Procurement, Construction Management (EPCM) services for this major capital project which, apart from its sheer size, scope, and engineering design, represented significant location, infrastructure, labour, environmental and climatic challenges. Today, the project has reached 90% of its nameplate capacity and is in its final performance optimization stage. This paper will discuss the significant project management effort that was required to successfully address the complexity of this major project. It will cover the main project management control systems and the engineering tools and processes used during the detail engineering, procurement, and construction phases to achieve optimum project execution.

Jan 1, 2015

By X. Xu, J. K. -W. Wong, J. -Y. Delenne, K. Soga

Hydraulic fracturing (HF) is a mulit-scale and multi-physics process which makes numerical modelling challenging. Conventional numerical models focus on multi-physics nature of HF but the multiscale nature is compromised by idealization of reservoir to be continuous, homogeneous and isotropic. However, natural fractures take a significant role in the success of HF for exploration of unconventional resources like Shale gas. Advancement in microseismic monitoring provides useful data for simulation in a much greater details. There is an urgent need to develop a suitable numerical tool for multi-scale threedimensional simulation to better understand the roles of discontinuities and heterogeneity in reservoir and its interaction with HF, especially to address the environmental concerns such as induced seismicity. In this paper, a simple discontinuum numerical method – Lattice Element Method (LEM) is proposed to model HF in a large scale three-dimensional model. A reservoir is modelled as a lattice network composed of 1D Hookean’s spring. Fracturing is modelled simply by removing lattices that meet a specified threshold as determined by the critical energy release rate of the rock. By introducing disorder in the model, mesh dependency in modelling fracture growth is minimized. Fluid flow is simplified as flow in pipe network and the permeability of pipe is related to fracture aperture by cubic flow law. Therefore, fracture flow and rock deformation are fully coupled. In this paper, four simulations of the same configuration except with different degrees of reservoir heterogeneity are presented. Fractures formed are diffusive and two groups of fracture are identified – connected fractures and isolated fracture without clear spatial correlation. The former forms a ‘fracture cloud’ and contributes to fluid flow and modeled by pipe network. The flow of fluid is highly tortuous and pipe network is sparely connected. The degree of reservoir heterogeneity controls the growth of ‘fracture cloud’ and pipe network. Microseismic monitoring often shows diffusive fracturing during HF on site and hence the outcome of the proposed model can be compared to such monitoring results.

Jan 1, 2015

By M. Gimpel, R. Galler, T. Antretter, H. Kargl, R. Tichy, W. Ecker, M. J. Mikl-Resch, G. Pittino

"A computational method based on a continuum-mechanical approach to simulate rock cutting is developed and evaluated. The aim is to obtain the local loading acting at the cutting pick during operation. The explicit FEM is chosen to model the pick-rock interaction. The presented method takes advantage of (i) an advanced material model (ii) a robust contact algorithm and (iii) fast calculations due to parallelization. A damage-plasticity law is chosen and modified with user-subroutines. A user integrated element deletion routine evades critical mesh distortions. Infinite elements at the model boundaries avoid unrealistic deformation wave reflection.INTRODUCTIONExcavation in mining and tunneling is a costly and complex process. Special equipment such as the roadheader (Figure 1) is employed to cut the rock in the excavation process but high material strengths limit their usability (Tockner 1999). Nowadays a challenge is to extend their field of application to harder and more abrasive rock materials. The forces on the pick (Figure 2) are of great importance and determine wear and especially the efficiency of the machine. Figure 2 shows the principle of cutting and the predominant forces. The variation of pick geometry, cutting depth or attack angle will result in an alteration of these forces and determine the efficiency of the process (Tockner 1999). The mechanisms of rock cutting follow the principle of indentation (Kou et al. 1999). Directly under the indenter a hydrostatic stress state evolves causing the material to fail only due to crushing and fragmentation mechanisms. In the close surroundings large strains emerge and tensile cracks are formed but in a further distance the material behaves elastically."

Jan 1, 2015

By Y. Guglielmi, C. Gout, P. Henry, P. Dick

The key questions about fault reactivation in shales concern the potential enhanced fluid displacement through a previously low-permeability aseismic shale formation. Here we selected a critically-stressed fault cross-cutting the underground research facility of Tournemire (France) which is located in Toarcian shales. Several inclined boreholes were drilled across the 8.5 m thick fault zone characterized by a fractured damage zone (FDZ) and a central 2.4 m thick complex core (FC) made of multiple shear bands with gouge intercalated with highly fractured lenses. The fault was pressurized using a straddle packer system defining a 2.4m long injection chamber normal to the fault surface. Three fault zone intervals were tested, two in the fault damage zone and one in the fault core. A three-dimensional displacement sensor set in the chamber allowed continuous synchronous coupled monitoring of fault movements, injection pressure and flow rate. An anticlockwise 0.22-to-0.05 mm fault slip was captured associated with a 10-to-100 hydraulic aperture increase of the fault. Preliminary analyses of the experiments show that leakage in the fault zone is controlled by strength and stress magnitude contrasts between the fault core and the fault damage zone. Fault damage zone discontinuities display a much higher peak static friction angle than the fault core shear zones, respectively of 34 and 22° that could be related to complex interactions between healing and sealing processes. Fault core leakage is much more dependent on strain than on effective stress variations. Discontinuous, pulse-like leakage, develops following complex failure processes that may simultaneously take place in the gouge zones and along the large pre-existing discontinuities. Such results may be of importance in evaluating fault seal integrity in oil and gas exploration, production or CO2 storage or underground excavations monitoring.

Jan 1, 2015

By J. Arzúa, L. R. Alejano, P. Rodríguez

"After six years of successful excavation with a support pattern, some roof failures have happened in a room-and-pillar underground mine excavated in a stratified rock mass. In this paper the authors give explanation to the failure mechanisms and redesign the support in order to solve the recent problems of roof collapse.The mine magnesite bed is approximately 14 m thick, and runs in the E-W direction with an average dip of 18º. Due to mine constraints and the inexpensiveness of the ore, the excavation is performed in the mineral itself, avoiding the need for mining dumps and waste excavation. Rooms and drifts are oriented in the mineral bed to form an angle of 19º to the direction of the strata, reducing the roadway gradients to less than 10%. Rooms are 11 m wide and pillars are 7 m wide. They are designed to leave a 1 meter thick magnesite bed in the roof of the rooms to act as a supporting beam. More often than not there is a layer of very poor quality marly shale up to 1 m thick resting on this magnesite beam. Overlying the poor quality marly shale there is a better quality, self-supporting sandy shale.First, a characterization of the involved rock masses (magnesite, marly shale and sandy shale) is performed, collecting the geotechnical data from previous studies and performing new in situ geotechnical surveys to check the soundness of the previous data. A distinction is made between different quality zones inside the mine by means of the RMR classification, associated to the occurrence of faults.Back-analysis of some of the failed rooms and data from exploration boreholes were used to determine the failure mechanisms (associated to Voussoir type failure or occurrence of sub-horizontal joints crossing the magnesite beam). Based on an understanding of these mechanisms, the method used to design the support and to analyse the stability of the different possible thicknesses of the magnesite beams of the rooms and drifts of the mine were improved.Numerical models of the most relevant conditions found in the mine were performed in order to better understand the behaviour of the roofs and to simulate the supporting solution obtained in the previous analyses.The most relevant conclusions with regard to the main elements causing the instability of the roof include the appearance of sub-horizontal discontinuities crossing the magnesite bed, which causes total instability of the roof and usually leads to collapse; excessive thinness of the magnesite beam in the problematic rooms; and the relaxation of the deepest (around 270 m in depth) zone of the mine because of existing faults."

Jan 1, 2015

By S. Gheibi, H. S. B. Duzgun, N. Dadashzadeh

It is widely recognized that there are various uncertainties in the analysis of rock slope stability owing to inadequate information for site characterization and inherent variability and measurement errors in geological and corresponding parameters. Therefore, reliability-based approaches that allow the systematic and quantitative treatment of these uncertainties have become a topic of increasing interest in rock slope engineering. There has been extensive reliability analysis of rock slope stability in the literature. However, all of these analyses are based on the First-Order Reliability Method (FORM). The probability of failure of a rock slope in the First-Order Reliability Method is provided by a linear approximation of the failure function. The FORM is most commonly used because it is efficient; its accuracy, however, deteriorates when the nonlinearity of limit-state function increases. A better approximation of the failure function can establish a more accurate value of the probability of failure. The Second-Order Reliability Method (SORM) overcomes this drawback with a cost of lower efficiency in terms of the number of function calls and computation expenses. Despite considerable studies of the application of the SORM in other engineering fields, a few attempts have been made in the slope stability studies. In this study, a recently published Second-Order Reliability Method with First-Order Efficiency (SORM-FOE) is implemented on a selected rock slope in Amasya, Turkey. The SORM-FOE overcomes the low efficiency of the SORM by decreasing the number of function calls. The probability of failure of both FORM and SORM-FOE are obtained and the results are compared. It is shown that, in high value of failure probability, SORM’s do not significantly affect the analysis and design. However, in very low failure probability, the discrepancy between FORM and SORM-FOE becomes pronounced. It is also shown that, where the second order methods are required, the SORM-FOE is more efficient than the conventional curvature fitting methods in terms of computational costs.

Jan 1, 2015

By Tanai Marin

Modelling of pyrometallurgical processes can become extremely challenging when the process involves complex chemical reactions, phase change, multiphase flows, heat and mass transfer and fluid flow. Traditionally modelling involving complex thermodynamic equilibrium and thermochemistry are treated separately from problems in which transport phenomena are dominant, such as CFD and heat transfer studies. The separation is also necessary due to the lack of modelling tools that can effectively couple both solution thermodynamic and transport phenomena or even multiphysics problems. However, there are, quite often, problems in high temperature process metallurgy in which it is not possible to simplify this two modelling approaches. Typical cases of this are for example, the modelling of slag freezing on furnace walls (electric or flash furnaces) to create a protective layer in front of the refractories hot wall, the formation of build up from the matte or metal phases at the hearth of the furnace, calcine smelting from charging piles into slag phase, etc. The problem is even more challenging in nonferrous pyrometallurgy as most, if not all, of the condense phases are highly non ideal solutions involving several components and exhibiting complex crystallization paths and miscibility gaps, such as the Fe-Ni-Cu-S matte, FeO-Fe2O3-SiO2, FeO-MgO-SiO2 slags, etc. In this paper, an approach to model coupled problems involving multiphysics and solution thermodynamic is presented. The approach consists in developing in-house a library containing the Gibbs energy minimization routines and using the local equilibrium concept for process modelling from a multiphysics commercial software (COMSOL). Simple test cases, such as building a non-ideal binary phase diagram on a CFD package is presented, the limitations and of other classic approaches such as equilibrium or Scheil cooling are discussed and compared to this model.

Jan 1, 2015

By D. Yuan, M. Sitter, S. St. Jean, M. Collins

The Ambatovy Joint Venture comprises an operating lateritic nickel mine, ore slurry pipeline, pressure acid leach plant and nickel refinery in Madagascar. The hydrometallurgical process selected by Sherritt for treatment of the Ambatovy ore is a long-term, sustainable approach to producing nickel. Sulphur, used to prepare hydrogen sulphide, which is used to produce mixed sulphide intermediate, is recovered in the refinery as ammonium sulphate fertilizer. Additional sulphur is used to generate sulphuric acid for treatment of the ore in the high pressure acid leach (HPAL) step; fossil fuel consumption is minimized by using the heat generated from the exothermic reaction between sulphur and oxygen in the acid plant to produce power and steam. Steam consumption is minimized by design of efficient flashing and preheat systems for HPAL, while the high temperature of the HPAL step serves to minimize acid consumption. Fossil fuel consumption by the operation as a whole is further minimized by locating the leach plant and refinery near the port of Toamasina and utilizing a 220 km slurry pipeline to transport ore from the mine to the leach plant. In addition, Ambatovy has demonstrated a strong commitment to sustainability programs such as the biodiversity management program at the mine and local procurement in Madagascar. In this paper, the hydrometallurgical process for nickel and cobalt recovery from Ambatovy ore is described, along with highlights of Ambatovy’s approach to sustainability in Madagascar.

Jan 1, 2015

By V. Prasad, A. Afacan, K. Popli, S. Dubljevic

Froth flotation is one of the most common and widely used separation processes in the mineral processing industry; however, its complexity means that complete understanding of the process has not yet been achieved. A dynamic fundamental model is needed to enhance the understanding of the process and to improve process control. It is challenging to capture all the phenomena involved; hence, we propose a multiscale modeling approach with a framework ranging from macroscopic (equipment scale) to microscopic (attachment, detachment) models along with distribution models for ore particle sizes. Real-time froth properties are obtained using a vision-based sensing system. These measurements can be used for developing a soft sensor for grade using statistical methods such as PCR (principal component regression) and PLS (partial least squares) regression. The fundamental model is updated in real-time using extended Kalman filtering. Attachment dynamics are studied for different particle size distributions. The modeling and soft sensing is tested on batch flotation cell for the pure mineral galena with different particle size distributions using Xanthate and MIBC as collector and frother, respectively. The model, online measurements and soft sensor are useful in developing monitoring methods for froth flotation.

Jan 1, 2015

By Marc-André Villiard, C. St-Onge, V. Silla, J-C Milot

Rio Tinto Fer et Titane (RTFT) operates a metallurgical complex in Sorel-Tracey, Canada. In the Ore Preparation Plant (OPP), ilmenite ore is roasted through four parallel rotary kilns generating 95% of the plant wide sulphur dioxide (SO2) emissions. RTFT proactively set an annual emission cap for SO2 to be achieved by 2015. To achieve this targeted SO2 reductions, a technology evaluation and selection exercise was carried out, followed by detailed engineering and implementation. After a detailed evaluation of SO2 scrubbing technologies, the Circulating Dry Scrubber (CDS) system was selected as the most appropriate desulphurization system for this application. The system was implemented and commissioned from 2013 to 2015. General performance testing conducted in May 2015 demonstrated that the system achieved and surpassed the anticipated design values in terms of SO2 removal efficiency and hydrated lime utilisation. This paper summarizes the SO2 technology selection process, design features of the selected technology, project implementation and final performance testing of the SO2 abatement system.

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 Xueyuan Nie, Xuezhi Zhang, Henry Hu, Shen Xinwei, Jimi Tjong, Li Fang

With the demand for weight reduction in vehicles and airplanes, researchers have shown great interested in casting wrought aluminum alloy 7075 because of its potential advantages such as high strength-to-density ratio and stress-reduction capabilities. It is therefore essential to conduct a fundamental investigation on solidification behavior and grain microstructure of alloy 7075. In this study, a commercially-available wrought aluminum alloy 7075 was used to study the solidification behavior of alloy 7075. A thermal analysis was performed based on the cooling curve and its corresponding first and second derivative curves. The nucleation of the primary a-Al phase and eutectic phases were analyzed. The comparison between the microstructures of the wrought and cast 7075 alloys shows that coarse microstructure and porosity are present in the cast alloy. The techniques for grain refinement in casting, such as the addition of the nucleating particles and applied pressure, should be employed to reduce the dendrite arm spacing and the porosity.

Jan 1, 2015

By W. Lin, Z. Z. Dong, C. L. Zhou, Y. C. Kuang

The inversion of rock mechanical properties based on well-logging data is a common method to drillbit selection in petroleum engineering, the basic principle of which is to use deep-well logging data and experimental regression model to calculate rock mechanical properties. The profile of rock mechanical properties can be established. According to the profile, the type of drillbit can be selected in different formations, which will improve the economic benefit of bit use indicators. Compressive strength of rock is one of the most important parameters. Conventional method to establish regression models of acoustic velocity and uniaxial compressive strength is on the basis of rock core sonic experiments and uniaxial compressive experiments. However, with the increase of deep drilling operation, triaxial compressive strength prediction model will play a more vital role in drillbit optimization selection in deep wells and improving ROP (Rate of Penetration). Triaxial compressive strength is an important parameter to evaluate the formation strength property. In this paper, to the sandstone commonly encountered in oil drilling, combined the rock mechanical experimental results with the changing porosity experimental results under confining pressure, a new interpretation model of triaxial compressive strength based on well logging data is proposed. Firstly, three types of sandstone samples are collected. Then measured components, uniaxial compressive strength, triaxial compressive strength and porosity. Experimental results show that: confining pressure has a significant influence on sandstone strength, with the increase of confining pressure, triaxial compressive strength increase; furthermore, the triaxial compressive strength of different porosity sandstones have different change rules when confining pressures changes; sandstone porosity decreases with the increase of confining pressure and it obeys logarithmic function. According to the results, regression models of sandstone, in-situ uniaxial compressive strength model and uniaxial compressive strength model under normal pressure are generated. Meanwhile, using well logging data in the acoustic velocity, density and gamma value can derive in-situ rock formations uniaxial compressive strength and in-situ porosity. Eventually, a new model established in this paper, containing in porosity, uniaxial compressive strength, lateral pressure and power exponent of regression, can predict triaxial compressive strength which will guide drillbit selection.

Jan 1, 2015

By Tassos Grammatikopoulos, Vladimiros Papangelakis, Douglass Duffy, Cheryl Washer, Srinath Garg

Flotation operations in Sudbury, Ontario, seeking to separate pentlandite and chalcopyrite from mixed sulphide deposits, produce pyrrhotite tailings that have been disposed of in tailings ponds for over 50 years. These reactive tailings are viewed as waste, despite containing ~0.8 wt% Ni, and represent an environmental liability because of their acid mine drainage (AMD) potential. Given the large quantities of pyrrhotite currently stored in Sudbury, and declining ore grades elsewhere, these tailings have shifted from being viewed as a waste to a potential resource. Recent advances in the bioleaching of low-grade sulphide ores have raised the prospect of bioleaching as a potential low-cost processing option. However optimization of this proposed process requires a thorough understanding of the mineralogy of the pyrrhotite tailings that would be leached. The work presented here addresses the characterization of two distinct Sudbury pyrrhotite tailings, namely from the Vale and Glencore operations, using electron microprobe, XRD and QEMSCAN analysis. Particular focus is given to how the mineralogy may affect bioleaching as a processing option, through the deportment of nickel and the various phases present

Jan 1, 2015

By Z. Zhong, Y. J. Hu, D. Elsworth

"A standard triaxial apparatus was set up to perform sliding tests. Slide-hold-slide (SHS) tests, with hold times from 10 to 3000s and at invariant effective confining pressures of 1MPa, incremented to 3MPa and then 5MPa, were conducted on both dry and wet split limestone samples to study the evolution of fracture friction under holding and sliding processes. Limestone fracture shows time-dependent behaviors in SHS tests, friction drops (creep) gradually during hold periods, and recovers (healing) after re-slide. Frictional healing and creep increment are proportional with hold times, while are inversely proportional with effective confining pressures. Moreover, the wet sample displays heavier healing than the dry one, whereas creep increment of the wet sample is lighter than the dry one. We conclude that the evolution of limestone fracture friction is significantly controlled by mechanisms of fluid-rock interactions.INTRODUCTIONThe evolution of the mechanical properties of fractures and faults exerts a critical influence on the behavior of fault zones and on the earthquake cycle. In addition to natural processes, this behavior has broad implications in induced seismicity and influences oil and gas production, the sequestration of carbon dioxide and the disposal of nuclear wastes.Fault friction is typically measured by slide-hold-slide experiments, where fault gouge is deformed at low velocities (<0.1 mms-1) and over small shear offsets (<30 mm) (Dieterich, 1979, 1981; Mair et al., 1999; Marone, 1998, 2004) or alternatively with rotary-shear devices, that enable higher slip velocities over very large shear offsets (Hirose & Shimamoto, 2005; Mizoguchi et al., 2006). Typical experiments (Dieterich, 1972; Marone, 1998; Olsen et al., 1998; Bos & Spiers, 2002; Tenthorey et al., 2003) show significant creep during the holding period and frictional healing (strength gain) upon reloading. Fault creep and healing are often time-dependent with longer hold times typically resulting in greater magnitudes of creep and in strength gains (Tenthorey et al., 2003; Yasuhara et al., 2005; Carpenter et al., 2012). Frictional healing may arise from growth of grain contacts (potentially induced by pressure solution) (Rabinowicz, 1951; Dieterich & Kilgore, 1994; Yasuhara et al., 2005), strengthening of the contacts (correlates with cementation) (Hickman & Evans, 1992; Beeler et al., 1994; Karner et al., 1997; Olsen et al., 1998; Muhuri et al., 2003; Tenthorey et al., 2003; Li et al., 2011) or via other changes in state including changes in porosity and rearrangement of force chains. Experimental evidence (Dieterich, 1979, 1981; Johnson, 1981; Marone & Scholz, 1988; Marone, 1998; Saffer & Marone, 2003; Moore & Lockner, 2007; Ikari et al., 2007; Ikari et al., 2011; Carpenter et al., 2012) quantifies both the mode and magnitude of second-order changes in the velocity-dependent (including velocity- strengthening and weakening) evolution of frictional strength."

Jan 1, 2015

By V. Shumila, T. Butler

Microseismic monitoring systems listen to a rockmass and record ground vibrations associated with brittle rock fracture. Modern systems record continuously and can identify as many as one thousand small and possibly larger magnitude events occurring each day from induced stress change and redistribution. This paper presents examples of the real time use of microseismic information in mining applications. It shows how alerts can quickly identify large magnitude events or unexpected clusters of seismicity and how seismic re-entry protocol can help to improve safety decisions. Additionally, newly developed advanced analysis techniques applied to mine seismic data sets are shown including the use of clustering methods to see into a cloud of seismic data to identify significant activity and quantify possible planar features. The use of source parameter trend analysis is shown to be a beneficial method for interpreting changes in relative stress levels. Seismic moment magnitude distributions are analyzed using the TGR (Truncated Gutenberg Richter) and TMD (Tapered Moment Distribution) methods to estimate bvalue and expected maximum magnitude. Seismic hazard analysis is shown as a useful way to simplify a dataset and highlight seismic hot spots. The rockmass response to excavation is a complicated spatial and temporal process and microseismic monitoring is a method that can be used to help understand the location and mechanics of failure occurring in the rock ultimately helping to improve safety and productivity at a mine site.

Jan 1, 2015

By M. Bustin, A. Hosseinian, G. Preisig

Hydraulic fracturing and hydraulic shearing techniques are essential for the development of unconventional oil & gas and geothermal resources. Such techniques will also play a key role in the development of CO2 sequestration, and management of high stresses encountered with deep mining. However, the optimal design of hydraulic injections is often limited by continuum treatments of the rock mass without consideration of the influence of the natural fractures present. These may interact with the hydraulic fracture to promote offsets or arrest, or may cause the natural fractures to undergo shear and dilation. Given the limited access to the rock afforded by deep boreholes common to oil & gas and geothermal projects, advanced numerical modeling is critical for understanding key mechanisms and interactions. In this work, we first present a coupled hydromechanical distinct-element model focusing on the propagation of a hydraulic fracture. Normal and shear displacements arising from the generated hydraulic fractures are subsequently integrated into a fracture flow model for detailed analysis of specific flow behaviour and fracture permeability. For this, a scanned 3-D profile of a natural fracture is used to correlate surface roughness to dilation as a function of shearing. Fracture flow is computed at different time steps corresponding to different fracture propagation states and geometries. Modeling results indicate that the opening of a tensile hydraulic fracture is accompanied by millimetre-scale shear displacements associated with slip and wedging between adjacent blocks. This may leads to a mismatch between the upper and lower walls of the fracture if asperities are irregular and strong enough. However, fluid-induced tensile fractures have low Joint Roughness Coefficients (JRC), and associated millimetre-scale shear displacements are too small to ensure a significant offset of asperities that prevent a fully-elastic closure when the hydraulic fracture is depressurized. Injection of proppant is thus critical to keeping hydraulic fracture apertures open. Future work is planned that will see the models calibrated against field data including microseismic event locations, magnitudes and source mechanisms, tiltmeter deformations, stress change, and pore pressures.

Jan 1, 2015

By L. J. Han, Y. J. Zong, S. Y. Wen, J. J. Wei

"Stress strain curves characteristics, deformation and strength properties, and failure process and characteristics of sandstone samples under different stress states were analyzed through triaxial compression tests carried on using TAW-2000 microcomputer control electro-hydraulic servo rock triaxial test system. The experimental results reveal that peak strength, residual strength, elasticity modulus and deformation modulus increase linearly with confining pressure, and failure models transform from fragile failure under low confining pressure to ductility failure under high confining pressure. Macroscopic failure forms of samples under uniaxial compression was split failure paralleling to the axis of samples, while macroscopic failure forms under uniaxial compression was shear failure, the shear failure angle of which decreased linearly with confining pressure. There were significant volume dilatation properties in the loading process of sandstone under different confining pressure, and then based on acoustic emission damage and volumetric dilatation damage the damage evolution properties of samples were analyzed, and damage constitutive model were established, realizing the real-time quantitative evaluation of samples damage state in loading process. INTRODUCTIONDamage is a phenomenon that micro defects in material under monotonic loading or reloading will lead to progressive decrease of the cohesion and damage of volume units. Many scholars at home and abroad have carried out systematic researchers on damage evolution characteristics and constitutive model of rock, and have achieved a series of remarkable results in this field. Lu Y.D. (2004) studied the complete stress-strain curve’s characteristics of marble under triaxial compression, and established the bilinear elastic-linear strain softening-residual ideal plastic damage constitutive model. Ren J.X. (2001) studied the damage evolution laws of coal and rock mass by computerized tomography triaxial loading system, realizing the quantitative evaluation of damage state. Zhang Z.L. (2011) studied the deformation and failure mechanism of strong weathered sandstone by triaxial compression test, and analyzed the damage evolution processes and established the damage evolution equations based on the density method. Li X.J. (2012) studied the microscopic damage characteristics of siltstone under triaxial compression by scanning electron microscope and digital image technology, and analyzed the statistical distribution characteristics of cracks azimuth angle, length, width. Zhou J.W. (2012) studied the strength, deformation and fracture damage characteristics of sandstone by uniaxial cyclic loading and unloading tests, and defined damage variable based on linear damage mechanics theory and acoustic emission (AE). Cao W.G. (2006) established a new damage model and a statistical damage constitutive model of rock under specific confining pressure according to the force of damaged part and undamaged part in rock on the base of statistical damage theory. Li S.C. (2014) studied the damage evolution characteristics in the whole uniaxial compression process of sandstone by electrical resistivity and AE, and proposed a status qualitative criterion for rock damage."

Jan 1, 2015

By R. Cammack

"A fundamental component of the Hoek–Brown empirical rock failure criterion is assessing an appropriate mi constant. In developing their empirical criterion Hoek and Brown normalised triaxial test results by dividing the principal stresses at failure by the uniaxial compressive strength to provide mi values for typical rock types. Experience shows that Hoek triaxial tests for individual lithology types commonly have a large inherent scatter. In addition, tests are often conducted on core from random locations, which although being of one parent rock type, display variations in density as well as minor variations in grain size and alteration.This paper presents experience from projects the authors have been involved in which large sample populations of triaxial and UCS tests on similar rock types were available. Geophysical logging data was also available for some projects. The methods used to estimate the compressive strength for each sample and the effects on normalising the data sets are described. The variability of the data and uncertainties and assumptions involved in the process of assigning a mi value from the normalised data are also discussed."

Jan 1, 2015

By Y. J. Yu, C. D. Zhang, Y. S. Zhao, W. G. Liang, Y. D. Geng, J. L. Bi

"In this study, the coefficient of thermal expansion (CTE) and thermal conductivity (TC) of shale in two perpendicular directions, vertical and parallel to the bedding of the rock, are firstly measured by experiment. The experimental results show that the CTE of the shale in the vertical direction is 1.69 times as much as that in the parallel direction, while the TC in the parallel direction is 1.83 times larger than that of the vertical direction. Based on the no initial strain assumptions, numerical simulations for the thermal cracking in shale are carried out with the obtained coefficients. Using finite element method and thermal-mechanical coupling theory, two fracture mechanical parameters, KI and J-integral of the shale are also analyzed for the shale thermal fracturing. The simulation results show that both KI and J-integral reach the maximum when the fracturing angle equals to 0° in the shale, namely parallel to the bedding direciton. It is also found that the two fracture parameters decrease with the fracturing angle increasing. Additionally, the thermal fracturing effect of the shale is evaluated under a given scenario. The study is significant for shale oil or gas recovery by in situ pyrolyzing method.INTRODUCTIONOil shale is regarded as a kind of unconventional resource that can contribute to production of shale oil and hydrocarbon gases (Rahm & Dianne, 2011). As is known that there is a large deposit of oil shale in China, and present technology for exploiting oil shale still concentrates on open-pit mining, obtaining shale oil by retorting shale on the ground (Li, Tang, Xue, Zheng, & Du, 2014; Qian, Yin, 2008). However, this technology has a bad effect on the environment and leads to severe resource waste for utilizing oil shale. Exploiting oil shale by in-situ method for ‘green mining’ is widely necessary in China (Kang, Wang, Cao, Liang, Chang, & Liu, 2013; Waters, Dean, Downie, Kerrihard, Austbo, & McPherson, 2009; Zoback, Kohli, Das, & McClure, 2012). One of effective method to exploit shale oil is heating oil shale in-situ underground. For instance, the in-situ conversion process (ICP technique) from Shell Company etc. used a direct heating method, and convectional heating oil shale underground by injecting superheated steam from Mining Technology Institute, Taiyuan University of Technology and constructing horizontal fracture in vertical wellbore from Chevron Company use indirect heating method (Chaudhary, Ehlig, & Wattenbarger, 2011; Hazra, Lee, Economides, & Moridis, 2013; Zhao et al., 2010). The convectional heating method contains complicated procedures with phase transition in double-media of fracture-matrix of oil shale subjected to multi-physical fields coupling including stress, fluid and temperature (Zhao et al., 2010). The characteristics of thermodynamics, solid mechanics and seepage of oil shale have been investigated separately (Hopkins, Rosen, & Hill,1998; Kang, Zhao, & Yang, 2008; Yang, Zhao, Kang, 2010; Zhao, Yang, Feng, Liang, & Kang, 2008; Zhao, Cao, Zhao, Lin, & Wang, 2008). However, it is quite difficult to solve constrained problem of each physical parameter in multi-physical fields coupling for oil shale pyrolysis."

Jan 1, 2015