Search Documents

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

By C. C. Li

Traditionally, strong and stiff devices, such as fully encapsulated rebar bolts and concrete lining and arches, are preferred for underground rock support. The attention is paid to the strength of the devices in such support designs. It has been observed in engineering practice that support systems composed of stiff devices are often subjected to premature failure in high stress rock conditions. In the past decade, yield support devices have been employed to improve the performance of rock support systems. A rock support system usually consists of both internal and external devices that are installed within the rock mass and on the free surface of the rock, respectively. Great efforts have been made to make support devices yieldable in the past years, but the efforts are focused on different aspects in the tunneling and mining branches. In the tunneling branch, the focus has been put on the development of external yield devices, such as Lining Stress Controller (LSC), Honeycomb element and compressible concrete, which are embedded in shotcrete liners, while in the mining industry the efforts has been on the development of yield rockbolts, such as the cone bolt and the D-Bolt. The internal and external support devices are not compatible in their deformation capacities in the current rock support systems even though yield devices have been employed. All the devices in a support system should be deformable and compatible in deformation capacity. With such a support system, premature failure of individual support devices could be avoided and the service life of the support system is prolonged. It is a trend at present to take into account the energy absorption capacity, instead of simply the strength, of support devices in rock support design.

Jan 1, 2015

By C. Grisanti, E. Prina Howald

The major of Swiss territory is occupied by mountains. Natural events as floods, debris flows, landslides, rockfalls, and avalanches regularly threaten the safety of people and properties. On the one hand, extreme weather events cause great damages, on the other hand, the enhancement of territory’s exploitation increases the risk. The Swiss Government makes many efforts to deal with the prevention of natural hazards, especially through measures of land management. As stated by federal law, Swiss cantons are obliged to establish hazard maps and to take the risk into account in the land management. Hazard maps are based on the hazard matrix, whose abscissa and ordinate are probability and intensity of the event on site. This project focuses on rockfalls. It consists of three main parts: the creation of a database to collect data on site, the creation of a catalog with all type of measures, and the development of a methodology for the analysis of the effectiveness of protective measures. The creation of the database allows us to collect information: about the site, about protective structures implemented on it, and about failures. With many multiple choice fields and few free text fields, it can collect data in a homogeneous way, without redundancy. The catalog has to be a tool that gives a basis to make decisions. It collects information about all types of protective measures implemented on Canton Vaud’s territory or available on the Swiss market. Free text fields allow for data storage of specifics, to give a complete view of structure’s characteristics and potentials. Methodology consists of four steps. The core of the procedure is the structure’s evaluation following a score system; it allows for calculation of reduction factors that have to be applied to intensity and probability parameters of the structure. Then it is necessary to convert the reduced structure’s intensity value into a residual site’s intensity value. Although the first step enables us to choose suitable measure for a specific site, it is not necessary to perform it. Second step allows us to evaluate structures under scenario 0, where structures are in best conditions. Third step allows us to evaluate structures under six scenarios of possible failings. Fourth step allows us to convert structure’s reduced intensity into event’s residual intensity and to establish the new point on hazard matrix. While the procedure has already been defined, details of score system and conversion have to be developed during the second phase of the project.

Jan 1, 2015

By M. L. F. Paullo, S. Cristian Mejia, D. Roehl

"Rock fracture simulation most often presents non-linear behaviour due to loss of material integrity and strain concentration in the fracture zone. These strains generate local instabilities that propagate to the global system of equilibrium equations, generating convergence problems of the numerical solution. Often, in such cases, analyses in the post-critical regime are unfeasible. To overcome the numerical difficulties associated with problems of stiffness loss in the post-critical regime control methods, also called continuation methods, are employed. Among them are the arc-length method, the energy control method and the indirect displacement control method. These methods are employed in conjunction with a Newton- Raphson scheme for the solution of non-linear systems of equations. However, depending on the system instability, convergence in not guaranteed even with very small increments. In the present work these control techniques are investigated in combination with cohesive elements to simulate rock fracture. The constitutive model employed with the cohesive elements represents stiffness degradation through a damage law, which leads to serious convergence difficulties as reported in the literature. Here the continuation methods mentioned are applied to a Mode I rock fracture problem and to hydraulic fracture simulations. The effectiveness of the different continuation methods is compared.INTRODUCTIONThe process of rock fracture induced by the fluid pressure is called hydraulic fracturing. The propagation of a hydraulic fracture is a complex process, which is basically defined by the mechanical deformation of the rock, the fluid flow within the fracture and the fracture propagation itself (Mokryakov, 2011). In the works of Bendezu et al. (2013), Chen (2012) and Carrier et al. (2012), vertical hydraulic fracture propagation is modeled by cohesive elements with traction-separation law. The softening behavior of rock in post-critical regime is demonstrated by experimental testing and is characterized mainly by high strains and low stresses (Crowder and Bawden, 2004). This nonlinear behaviour can arise due to material nonlinearity, such as micro-cracking and damage, leading to softening behaviour (Chandrakant, S.D., 2012). The cohesive fracture model considers in its formulation the softening behaviour that most rocks present in post-critical regime with propagation criteria."

Jan 1, 2015

By M. Hozayen, E. S. Radwan, P. A. Nawrocki, A. Jebbouri

"A major geomechanical study was conducted for one of the UAE onshore oil fields with the purpose of studying reservoir rock properties and stress distribution, based on which a 3D coupled reservoir geomechanical model was built. This model was later used for wellbore stability analysis and for analyzing the effect of fractures on future reservoir production performance, including fracture deformation and permeability evolution occurring due to reservoir depletion and resulting changes in local stress field. The subject of this study was a complex fractured Upper Cretaceous carbonate reservoir characterized by high degree of heterogeneity, both vertically and laterally. Wellbore instability related drilling problems were encountered while drilling the wells in this field. The problems were mainly in drilling the horizontal wells, including breakouts and mud losses. Experimental results of laboratory tests are summarized in this paper whereas wellbore stability analysis and formulation of a 3D geomechanical model that was built for that purpose is summarized in a sister paper (Part II).A suite of tests reported includes tensile strength, uniaxial and triaxial compression and permeability under stress at simulated reservoir conditions. All tested rocks have been grouped into twelve test units and reservoir geology, rock types and challenges related to designing a viable sampling and testing program for a big reservoir are discussed first followed by sample selection criteria and details of the testing program. All test units showed brittle behavior at zero confinement and some samples showed similar behavior at low confining stresses. With confining stress increase, the behavior gradually turned to plastic but some samples showed features of ductile behavior even at low confinement.Obtained experimental data have been summarized in the form of strength criteria and correlations. The Mohr-Coulomb strength equations were obtained for each test unit and also for the averaged values. It was noticed that triaxial compressive strength (TCS) increases with decreasing porosity or increasing density but their influence on TCS is getting smaller with confining pressure increase. Good correlation of porosity and density with Young’s modulus, Poisson’s ratio, unconfined compressive strength (UCS), friction angle, cohesion and tensile strength has been noticed and there was an increasing trend between the confining pressure and TCS. Rock type (mud or grain supported) was not influencing the TCS as much as porosity and density and the two lithologies encountered affected TCS approximately in the same way. These results and correlations are summarized and explained in this paper. They were further used in the numerical modelling part (Part II) of this paper, Hozayen et al. (2015). It is also concluded that rock lithology is not a decisive factor and mechanical properties of tested rocks should be rather considered dependent on porosity and permeability. This helped greatly in reducing the number of tests that had to be performed."

Jan 1, 2015