Search Documents

By T. Feldmann, G. P. Demopoulos

The adsorption of gold complexes on inorganic metal oxides has been mainly discussed in geochemical literature, but it is also of interest to the hydrometallurgical industry. This is so because it could be responsible for valuable metal losses or because it could potentially serve as alternative adsorbent material to activated carbon in certain cases such as from hot and acidic process slurries. This paper constitutes an early investigation into this subject by focusing on the effectiveness of FeOOH, Fe3O4 and Al2O3 as sorbent materials. All tests were done in a batch stirred-tank reactor under atmospheric conditions with synthetic gold(III)-chloride solutions at 23°C and 80°C, in the pH range of 1.6 to 9.0. Detailed kinetic tests were performed at different pH values (3 and 6) on FeOOH and Fe3O4 with solutions having gold concentrations from 5 mg/L to 130 mg/L. The adsorption of gold was found to be best both in terms of kinetics and loading characteristics at 80°C and pH=6. This was explained by a change in Au(III)-Cl speciation from a gold(III)-chloride complex to a gold(III)-chloride-hydroxy complex. By far, magnetite (reagent grade chemical) yielded the best gold adsorption characteristics making it an interesting sorbent material for further development considering that it can be advantageously recovered from slurry streams via magnetic separation. On the other hand, the tendency of gold(III)-chloro complexes to adsorb on metal (hydro) oxides needs to be taken into account to minimize undesirable metal losses in chloride-containing leaching media.

Jan 1, 2015

By Huan-Qiang Zhang, Gang-Hai Huang, Yu-Yong Jiao, Fei Zheng, Zhang. Xiu-Li

An efficient and accurate three-dimensional (3-D) contact detection algorithm is of vital importance in large-scale computing of discontinuous deformation analysis. The contact detection algorithm takes up a large proportion of the whole calculation time because of the increasing calculation number in three-dimensional case. For block systems which consist of convex polyhedrons, the efficiency and accuracy of the contact detection algorithm is a great challenge. The presented contact detection algorithm in discontinuous deformation analysis (DDA) method is carried out in three phases: the global neighbour searching, the local contact type searching and the contact transferring. A contact detection algorithm for 3-D DDA is presented with some improvements on the three phases. A grid-based global neighbour searching algorithm is used to eliminate unnecessary neighbour blocks to the exact number. The local contact type searching includes two steps. In the first step, the common plane method is used to determine the contact type between convex polyhedrons. In the second step, the contact pairs detected by the common plane method are divided into four basic types. In the third phase, the four contact types are transferred into two basic contact types, vertex-to-face contact type and edge-to-edge contact type, which are managed in a uniform procedure. Three verification examples are calculated, the results show that the proposed algorithm based on the three contact detection phases is robust and efficient.

Jan 1, 2015

By S. Heusermann, S. Fahland, J. Behlau, A. Schäfers, J. Hammer

In the course of the deep subsurface use, the Federal Institute for Geosciences and Natural Resources has been carrying out extensive geological and geomechanical research and practical project work on domal salt structures, to prove their suitability e.g. for the disposal of radioactive wastes. Recent work on disposal sites has been taking into account the 3-D geological and geomechanical modelling of the underground structures. The geological models are generated using the special-construction openGEOTM code to improve the visualisation and interpretation of the geological data basis, e.g. borehole, mine, and geophysical data. For the geomechanical analysis of underground structures the JIFE finite-element tool has being developed to analyze large 3-D structures and to consider coupled THMC processes. The results of finite-element calculations, i.e. stresses, strains, displacements, and dilatant zones, are basic values for the evaluation of the stability of the structures and the long-term integrity of the geological barrier. As an example, the results of geological and geomechanical investigations of the Morsleben repository are presented including three-dimensional modelling of the entire geological structure of the central part of the mine as well as three-dimensional geomechanical model calculations to provide the necessary results for a long-term safety prediction.

Jan 1, 2015

By Randy Settgast, Daniel Moos, Hyunil Jo, Ghazal Izadi, Christine Baba

"Analyses of production logs of multiply fractured horizontal wells in unconventional plays show that only a small percent (~30%) of perforation clusters contribute to production. We present in this paper one possible approach to mitigating this problem, based on use of a fully 3-D numerical model of fracture propagation through a heterogeneous medium, by studying the influence both of spacing and of perforation friction on the initiation and propagation of simultaneous hydraulic fractures from a horizontal wellbore. The results indicate that while limited-entry perforating can help to equalize the lengths of fractures simultaneously growing from multiple clusters, it cannot equalize their widths. Both fracture width and fracture length can be equilibrated by combining non-even cluster spacing with limited-entry perforating. This has obvious benefits as it allows equalized placement of proppant across the entire series of perforation clusters, resulting in similar fracture widths and conductivities which enhances the likelihood that all clusters will contribute equally to production.INTRODUCTIONThere is substantial evidence from production logs and other information that during the hydraulic stimulation of horizontal wells, not all perforation clusters within a given stage will produce. One possible cause is the interaction between the hydraulic fractures emanating from each perforation cluster. Typically fractures created at the outer perforation clusters create “stress shadows”, i.e. increases in compressive stress in the minimum principal stress direction, which limit the growth of fractures growing from interior clusters; this is confirmed by modelling [1, 2, 3]. An example of this can be seen in Figure 1 which shows the final configuration of simultaneously propagating fractures from multiple, evenly spaced clusters when there is no perforation friction. The inner fractures are substantially shorter, reproducing the results by Bunger and Peirce [3]. The inner fractures also have substantially smaller widths.One approach to equalize the growth of simultaneously propagating fractures from a horizontal wellbore is by using limited-entry perforating as a method to equilibrate the flow rates into each of the different perforation clusters. Perforation friction is caused by the finite diameter of the entry points into the formation, and leads to a pressure drop across each perforation which is proportional to the flow rate squared. Fractures that take more fluid have a higher fluid velocity through the perforations than fractures receiving less fluid. This causes a larger pressure drop between the wellbore and more rapidly growing fractures, compared to the same pressure drop between the wellbore and the fractures that grow more slowly. Because pressure within the wellbore is approximately the same at each entry point, the larger pressure drop leads to less pressure being delivered to rapidly growing fractures, thereby decreasing their ability to propagate and reducing the amount of fluid they are able to receive. Because the mass rate of injection into the wellbore is constant, this forces more fluid into the fractures that were propagating more slowly, increasing their growth rate. It also increases the pressure required to drive fluid injection. The balance of these effects makes the system largely self-regulating and delivers more balanced fracture growth. While this effect is well known, it is often neglected in current numerical models."

Jan 1, 2015

By V. S. Gischig, G. Preisig

Hydraulic stimulation of deep reservoirs is an essential procedure for developing engineered/enhanced geothermal systems (EGS), but also in other Geo-energy contexts. The goal is enhancing permeability to such a degree that fluid flow rates within the reservoir are sufficiently high for an economically productive reservoir. Two fracture stimulation mechanisms are distinguished that may occur concurrently during the hydraulic treatment: 1) During hydro-fracturing (HF) new tensile fractures are propagated from the borehole by means of a fluid pressure overcoming the minimal principle stress d3 plus the tensile strength T0 of intact rock. In case of a pre-existing fracture oriented normal to d 3 , the fracture can be opened by exceeding d 3 only. Fluid injection is performed via small packed intervals for creating a stack of hydraulic fractures. 2) During hydro-shearing (HS), over-pressure induces slip along pre-existing fractures that are favourably oriented in the stress field for reactivation in shear. Such stimulation is usually performed in a large packed interval or in an open borehole. Which mechanism occurs predominantly during stimulation depends on the rock mass structure and in-situ stress field, but also on the orientation of discontinuities intersecting the open-hole section and thus being pressurized. Both mechanisms also differ in how they affect fracture permeability. While permeability gained by HS is mostly irreversible due to rearrangement of asperity contacts accompanied with shear dilation, permeability enhanced during HF reduces nearly reversible after pressurization unless proppant is used to ensure permanent apertures. In this contribution, we critically discuss the pros and cons of both mechanisms regarding efficiency of enhancing permeability and the associated hazard of induced earthquakes based on literature review and numerical modelling. We review literature that reports on how much permeability required for productive EGS reservoirs and how much permanent permeability can be gained by stimulation. We also compile information regarding the degree of seismicity induced during HFor HS-dominated stimulation procedures, together with conceptual studies that reveal characteristics of seismicity associated with the two mechanisms. Many observations and models indicate that HF may have a higher tendency of being aseismic, while felt seismic events are usually associated with HS. Further, we present a hydro-mechanically coupled fracture flow model that investigates, if and to what degree the stimulation strategy can be designed such that one of the two mechanisms is evoked and dominates over the other. It shows that HS may dominate in presence of larger persistent fractures nearly optimallyoriented in the stress field, even if HF is attempted in short packed intervals. The model further demonstrates that stress transfer during HS promotes the development of permeability across an anisotropic layer of the stimulated rock instead of a large volume. Hence, the study sheds light onto the feasibility of creating productive EGS reservoirs in crystalline rock at several kilometres depths.

Jan 1, 2015

By Andrea Lisjak Bradley

Clay shales are currently being assessed as possible host rock formations for the deep geological disposal of radioactive waste. However, one main concern is that the favourable long-term isolation properties of the intact rock mass could be negatively affected by the formation of an excavation damaged zone (EDZ) around the underground openings. The goal of this thesis was to develop and validate a mechanically-based model of the deformation and failure process of a clay shale, namely Opalinus Clay (OPA), with particular focus on the influence of anisotropy on the short-term response of circular tunnels. To achieve this goal, a hybrid continuum-discontinuum numerical approach was developed in combination with new field measurements from the Mont Terri underground research laboratory. The response of OPA during the excavation of a full-scale emplacement (FE) test tunnel was characterized by geodetic monitoring of wall displacements, radial extensometers and longitudinal inclinometers. The deformation measurements indicated strong directionality induced by the combined effect of in-situ stress field and bedding planes striking parallel to the tunnel axis, with the most severe deformation occurring in the direction approximately perpendicular to the material layering. Computer simulations were conducted using a newly-extended combined finite-discrete element method (FEM/DEM), a numerical technique which allows the explicit simulation of brittle fracturing and associated seismicity. The numerical experimentation firstly focused on the laboratory-scale analysis of failure processes (e.g., acoustic activity) in brittle rocks, and on the role of strength and modulus anisotropy in the failure behaviour of OPA in tension and compression. The fracturing behaviour of unsupported circular excavations in laminated rock masses was then analyzed under different in-situ stress conditions. Lastly, the proposed modelling methodology was applied to the aforementioned FE tunnel to obtain unique insights into the EDZ formation process around emplacement tunnels for nuclear waste. The calibrated numerical model suggested delamination along bedding planes and subsequent extensional fracturing as key mechanisms of the damage process potentially leading to buckling and spalling phenomena. Overall, the research findings are already having a direct impact on the preliminary design process of an underground repository for nuclear waste in Switzerland.

Jan 1, 2015

By S. Perez

Core drilling is a fundamental process in the mineral exploration industry. As deep exploration in hard rocks is accelerated due to the depletion of near ground mineral resources, deep-drilling operations need to reduce their reliance in operators experience and begin to implement new technologies in order to increase efficiency and cost-effectiveness of the overall process. Currently, lack of real-time information on the bit/rock interaction is a major problem. Therefore the mechanism behind the acoustic emission, which is a proven technique for remote monitoring, generation in mechanical rock cutting must be understood in order to provide real time information about the rock. This paper aims to identify the influence of mineral composition of granitic rocks on the acoustic emission (AE) during rock abrasivity test with a single contact tool and to evaluate the feasibility of AE as a means to infer Cerchar Abrasivity Index (CAI). Single contact tool will provide precise and accurate description of the source of AE. To provide constant grooving speed Cerchar apparatus was modified to monitor AE during rock/tool interaction. Tests on the modified Cerchar apparatus are carried out under the static load of 70N and cutting speed of 1 mm/s. A Scanning Electron Microscope (SEM) is utilized to visualize the groove in the rock and to identify the mineral composition through Back-Scattered Electrons (BSE) reflected from the rock sample. Based on the microscope images, it was found that biotite shows a larger groove area and plate-like fractures accompanied with burst-type AE signal, whereas plagioclase and quartz, which are characterized by similar hardness and BSE brightness, possess shattered-like fractures and more continuous-type AE signals. In addition, a-value of the Gutenberg-Richter scale has been found to well correlate with Cerchar Abrasivity Index (CAI).

Jan 1, 2015

By Bob Drew, Goran Vukovic, Bojan Zivanovic, Klaus Gamweger

Reaction bath temperature is a major factor in processing liquid metals and in many applications where a balance between process duration and temperature has to be found to run a procedure efficiently. Also extremes of temperature can decrease the output of a process if too cold or even make it impossible, conversely overheating can be detrimental to both the product and furnace integrity. With the help of the new online temperature measurement system very accurate monitoring of the process is possible in contrast to the current routine practice where the temperature is measured periodically with dip thermocouples. The system works with a digital, fiber optic pyrometer, which is embedded in a refractory housing. Installation in the lining is possible below liquid level as nitrogen is blown to keep the system open and to ensure continuous sight from the optical head of the pyrometer to the molten metal through the tuyere. This permanent temperature measurement gives optimal process control with the potential of taking appropriate actions at an early stage enabling proactive process control. Therefore energy savings and even process endpoint detections can be realized.

Jan 1, 2015

By Q. Li, X. -L. Lei

Accompanying global energy shortage and high energy prices, shale gas has been the focus of exploration and it has attracted great attention from both academia and industry. The combination of horizontal drilling and hydraulic fracturing enables the extraction of huge quantities of shale gas from impermeable shale formations. During the process of hydraulic fracturing, great amounts of water, sand and chemicals must be injected into the shale formation at a high pressure, which may trigger seismicity. Such events could occur due to the presence of faults or the reactivation of these potential faults. Shale is the main component of source rock in the Sichuan Basin, China, and has the greatest contribution to the natural gas of the basin. In order to investigate the fracturing mechanism of shale in the Sichuan Basin, China, a series of rock fracture tests under triaxial compression were carried out in the laboratory in order to investigate the hydraulic fracturing behavior in relation to the seismicity. Based on the high multichannel waveform recording system, the detailed spatio-temporal distribution of acoustic emission due to micro-cracking was used to examine the deformation and fracture behavior of shale. In addition, the effects of confining pressure on the deformation and fracture in shale under triaxial compression by acoustic emission monitoring and X-ray computed tomography were discussed in the paper. The experimental results indicate that the shale samples demonstrate brittle fracture behavior. The confining pressure has a great effect on the complexity of the fracture planes and the spatial distribution of acoustic emissions, but few contributions to the location of the major hydraulic fracturing plane. The foliation of shale plays a dominant role in the rupture process and the geometry of the major hydraulic fracturing plane. The high confining pressure increases the strength and restrains the deformation of shale during compression process. These results shed some light on the understanding of onset of hydraulic fracture initiation not only for shale gas exploration and production, but also for CO2 geological storage and underground energy storage.

Jan 1, 2015

By R. Simon, A. Khosravi

"Tensile strength of rock samples is measured by direct and indirect tensile tests in accordance with standard guidelines. Due to the difficulties in performing direct tensile test, the indirect tensile (Brazilian) test is conventionally carried out to obtain the ultimate tensile strength. A limited number of studies have been conducted to investigate the role of sample size (length to diameter ratio = L/D) in the stress distribution of the Brazilian test using numerical and analytical methods. These investigations showed that the tensile stress distribution along the loading diameter is not uniform. They also indicated that the difference between the tensile stress at the middle part and at the sides of the sample with higher L/D ratio is greater than that in the other samples.The objective of this study is to show the effect of tensile stress distribution on the fracture surface created by the Brazilian test from a morphological point of view on the samples with different L/D ratio. Hence, a series of fracture surfaces created using the Brazilian tests on the rock samples with L/D of 0.5, 0.75, 1, 1.25 and 1.5 were investigated.In this study, all fracture surfaces were scanned by a 3D laser profilommeter. A number of roughness profiles parallel and perpendicular to loading axis have been selected from the samples with different L/D ratio and spacing of 1 mm. The results show that the roughness of the profiles perpendicular to the loading axis increases with an increase of L/D ratio, while the roughness of the profiles parallel to the loading axis decreases. Furthermore, some roughness profiles parallel to the loading axis were chosen and divided into three segments. The results revealed that the segment in adjacent to the upper loading platen possesses higher roughness than the other segments."

Jan 1, 2015

By K. K. Panthi, P. K. Shrestha

"Displacement of rock mass around a tunnel opening varies according to the rock mass properties, insitu stress conditions and applied support in the tunnel. Such displacement is also altered by presence of bands of rock mass that differ in properties. This paper focuses on the deformation behavior of rock mass around tunnel openings at four tunnel sections in Kaligandaki headrace tunnel where extensive instrumentation was done. Based on the actually measured tunnel displacement records and actually measured support pressure; back analysis has been done to estimate rock mass properties and insitu stresses. Numerical analysis has been done to analyze and assess the deformation behavior of the rock mass around the tunnel sections.INTRODUCTION Instrumentation and monitoring of tunnel convergences have significant importance in evaluation of tunnel stability. One of the key aspects of monitored data is that early convergences can be used in estimation of long term convergence. In addition, these data can also be used in evaluation of rock mass parameters and in-situ stresses around the tunnel. Convergence around a tunnel periphery may vary according to rock mass quality. Presence of bands of different rock mass may hinder the displacement pattern in the tunnel. Particularly, displacement characteristic of tunnel periphery in schistose, foliated and deformed rock mass is altered by presence of bands of comparatively stronger rock mass, e.g., quartz veins or presence of shear seams or faults. Similar observations were also recorded in monitored tunnel sections in Kaligandaki headrace tunnel, located in the Nepal Himalayas. Extensive monitoring was done by installing multiple point borehole extensometers (MPBX) at three different locations of the tunnel periphery at four different tunnel sections. In addition, convergence measurements by tape extensometers and installation of pressure cells to record support pressures were also carried out at each of the tunnel sections. Observations showed that the tunnel convergences varied according to rock mass quality, in-situ stress and effective support pressure.Rock mass parameters and stresses have significant role in the extent of tunnel deformation, but these parameters are often unknown or difficult to predict. Back analysis based on actual tunnel displacement records and known support pressure is believed to be an important method in such endeavor. This paper performs back analysis to evaluate rock mass parameters and discusses correlation of rock mass behavior with tunnel deformation and applied support. For this purpose, the rock mass parameters are estimated from actual tunnel displacement records, mapped rock mass quality records and laboratory tested data. Then, numerical modelling is conducted to correlate post peak behavior of the rock mass with actual tunnel displacement."

Jan 1, 2015

By J. Bensabat, B. Figueiredo, A. Niemi, J. Rutqvist, C. F. Tsang

"A hydro-mechanical model was developed to study the interaction between mechanical deformation and fluid flow in the fault zones during CO2 injection and storage at the Heletz site which was chosen as a test site for CO2 injection experiments in the framework of the EU-MUSTANG project. The potential reservoir for CO2 storage at the Heletz site consists of three sand layers separated by lowpermeability shale layers. The potential consequences caused by fault reactivation, focusing on the shortterm integrity of CO2 storage, are evaluated. The difference in the results obtained by considering the three-layer as an equivalent layer storage formation is analyzed.INTRODUCTIONLarge-scale storage of CO2 in deep underground reservoirs may cause considerable pore-pressure perturbation as well as flux migration and concern has been raised over whether shear slip displacements may be induced in nearby faults, as a consequence of their reactivation. A state-of-the-art review can be found in Rutqvist (2012), which describes the various possible effects resulting from CO2 injection into a deep sedimentary basin, including those involved in reservoir stress–strain changes, micro-seismicity, caprock sealing performance, and potential fault reactivation. Moreover, several studies have recently shown how CO2 injection may produce seismic events both on major faults (i.e., large faults with large initial shear offset) (Cappa & Rutqvist, 2011a, b, 2012), and on minor faults that might have gone undetected during the site characterization, e.g., faults less than 1 km long and without detected initial offset (Mazzoldi et al., 2012).The present paper is concerned with the Heletz formation which was chosen as a test site for CO2 injection experiments in the framework of the EU-MUSTANG project. The Heletz site located at the southern Mediterranean coastal plain of Israel and the structure is an anticline fold with a crest of about 4 km by 2 km. The structure is gently dipping to the east, truncated by a pinch-out line to the west and subdivided into a number of blocks by transversal normal faults with small displacements. The water table is located at 300 m below the land surface. An areal sketch of the site with the elevation of the caprock and the location of the borehole that has been chosen for the CO2 injection experiment is presented in Figure 1."

Jan 1, 2015

By Arthur Barnes, Mika Muinonen, M. J. Lavigne

"The carbothermic reduction of chromite ore to produce high carbon ferrochromium is among the most energy intensive metal extraction processes performed, due not only to the highly endothermic nature of the reduction reactions but also the very high operating temperatures required for the smelting operation to separate the chrome-containing alloy from the discard slag.KWG Resources and XPS have collaborated over the last 3 years in developing a patented (Winter) processing alternative to the traditional smelting operation in which all processing takes place in the solid state at temperatures many hundreds of degrees lower than those required for liquid processing. Amongst the attractions of the process are much lower capital costs, as well as reduced operating costs, particularly of energy.A techno-economic study on the new process identifies a reduction in overall energy consumption of 70% against conventional processing and approximately 40% lower than the current best practice. By utilising more natural gas than coal based energy sources, the KWG process is able to show an overall 50% reduction in greenhouse gas emissions compared with the most energy efficient current practice. The impact of the new process on future processing is therefore regarded as highly significant, with global energy reductions equivalent to the effect of completely eliminating energy demand from a country the size of Italy.The study concluded that the potential for this process to completely revolutionise the global ferrochrome industry should not be underestimated. Its impact reaches far beyond exploitation of a regional natural resource."

Jan 1, 2015

By T. Banerjee, D. Felsher, G. Leavoy

"Lake Shore Gold is engaged in commercial gold production at both the Timmins West and Bell Creek mines, located in Timmins, Ontario. Run of mine ore is processed at the Bell Creek mill which consists of a conventional gold mill flow sheet involving crushing, grinding, gravity concentration and carbon-in-leach (CIL) for gold recovery. Since 2009, the operating capacity of the mill has been upgraded from 2000 tpd tonnes per day (tpd) to over 3,000 tpd. The most recent expansion involved the installation of a low aspect ratio SAG mill to replace the secondary crusher. Without operational data on the low aspect mill, combined with a site team that was new to SAG milling, there were many challenges in achieving the target throughput from the circuit. Through the implementation of AwaRE advanced controllers, running in closed loop with PRC fragmentation cameras as well as embedded PARC modules within the Allen Bradley PLC, the site was able to rapidly develop a best operating practise and has been successful in creating a stable and consistent operational environment. The process, lessons learned and eventual solution will be reviewed in this paper.INTRODUCTIONLake Shore Gold is engaged in commercial gold production at both the Timmins West and Bell Creek mines, located in Timmins, Ontario. Run of mine ore is processed at the Bell Creek mill which consists of a conventional gold mill flow sheet involving crushing, grinding, gravity concentration and carbon-in-leach (CIL) for gold recovery. Since 2009, the operating capacity of the mill has been upgraded from 2,000 tpd tonnes per day (tpd) to 3,150 tpd. The most recent expansion involved the installation of a low aspect ratio SAG mill to replace the secondary crusher. This endeavor was split into two phases: the expansion of the back-end (Phase 1) and the front-end (Phase 2) of the mill. Phase 1 of the expansion was completed at the end of 2012 and involved enhancements and capacity addition to the leaching and adsorption circuits, which increased the mill’s capacity to 2,500 tpd. Phase 2 of the expansion commenced at the end of July 2013 and involved the expansion of the crushing and grinding circuits. The expanded front-end of the mill now has a capacity of 5,500 tpd, and processes ore from both Bell Creek and Timmins West mines. The updated flow-sheet, as shown in Figure 1, illustrates the replacement of the two-stage crushing and ball milling circuit with a single stage crusher and a low aspect ratio SAG mill."

Jan 1, 2015

By D. Moronkeji, R. Maharidge, D. A. Curry, B. J. Davis, U. Prasad, J. Franquet

"Proper shale characterization is crucial in shale gas drilling, exploration, completion and production. Unlike porous sandstone and carbonate reservoirs, the challenge with shales comes from their micro- to nano-darcy permeability ranges and the anisotropy of their mechanical and hydraulic properties. Most porosity-dependent logs are not sensitive enough to capture the true behavior of shales. Consequently, mineralogy and other strength and stiffness-related properties must be determined to characterize magnitudes and variances in shale mechanical and flow properties. In addition, these properties must pass quality assurance/quality control (QA/QC) practices to give confidence in their use in subsequent log calibrations, during geomechanical modeling, or in fracture design models. This paper presents a technique to integrate shale properties that are measured by various methods, and a methodology to perform QA/QC on the data before application.Shale properties originate from various sources, including downhole and laboratory measurements. Physical samples for laboratory testing are usually limited to a few depths. Lab tests are often costly, and replicating in-situ conditions can be difficult or impossible. In addition, results frequently show wide scatter, reflecting inherent depth variability and lateral heterogeneity. Wireline logs provide nearcontinuous data, but rarely yield direct measurements of rock strength properties. Although geomechanical and fracture models use log data, calibrating it with reliable and previously verified lab-measured properties is important. This paper presents a process to perform QA/QC on lab-measured data using various crossplots.Ten mechanical, microstructural and physical properties are graphically displayed to visually characterize shale strength behavior in ambient and in-situ conditions, whether measured or modeled statically or dynamically. These properties are often sourced from different departments. The properties are integrated but are individually subject to the previously mentioned uncertainties. Pairs of property values are crossplotted and compared with correlations from the literature or previously established with qualified data from related formations. Conformance with these correlations gives confidence in the data, whereas outliers are flagged as questionable. Repeating the process with questionable data crossplotted against other parameters for the same depth highlights the property with the greatest uncertainty. The determination of its value is investigated to locate the source of the inconsistency and establish a more reliable value. Where this is not possible, the user has a rational basis for rejecting that property value."

Jan 1, 2015

By K. E. Waters, N. Brodusch, R. Gauvin, A. Jordens, Teng. C., H. Demers

Rare earth elements (REEs) become increasingly significant for renewable energy devices and high-technology electronic products because of their large-scale popularity and applications. In this project, a field-emission scanning electron microscope was used to characterize the composition and distribution of REEs-bearing minerals. X-ray elemental maps were acquired with a large area x-ray detector (60 and 80 mm2 collecting area) and annular detector (60 mm2 collecting area), which allow fast x-ray acquisition at high count rate. These elemental maps were transformed into phase maps, which are used to identify the composition of the mineral phases and the REEs distribution over a large area. The phase maps are an easy way to distinguish the phases with and without REEs and compare the efficiency of different REEs extraction processes. Different methods (manual and automatic) for producing phase map were studied by comparing the distribution of REEs obtained from the same area.

Jan 1, 2015

When underground extractions occur beneath or in close proximity to valleys and other forms of irregular surface topographic features, the observed vertical subsidence at the base of the valley is less than that would be expected in flat terrain, while the valley sides are observed to move inwards towards the centre of the valley. The reduction in subsidence at the bottom of the valley is referred to as upsidence and the convergence of two sides of the valley is termed valley closure. Research coupled with mining experience over the past few decades in Australia has identified valley closure subsidence as being very real and quite significant part of non-conventional subsidence effects, which always have adverse effects on natural features, with risks to water flow and aquatic ecosystems. In this paper, we describe a threedimensional numerical modelling approach to investigate the geomechanical mechanisms which cause the valley closure subsidence behaviour, using distinct element computer package 3DEC. The calibrated models can accurately predict subsidence and horizontal displacements, and they are subsequently extended to specifically study the effect of a series of mining geological and topographical factors on valley closure subsidence phenomenon, by developing a range of hypothetical scenarios. Previous analyses carried out using a two-dimensional numerical modelling technique have identified the potential parameters that affect the valley closure subsidence, which mainly include longwall position, horizontal stresses, valley sloping angle, and valley shape. Further three-dimensional analyses regarding the relationships between the major factors and valley closure subsidence are discussed in this paper. The results have significant implications with regard to both the underlying mechanisms that lead to this nonconventional subsidence behaviour, and how these should be incorporated in future valley closure subsidence prediction and mitigation of its impacts on natural features.

Jan 1, 2015

By H. M. Kim, J. Rutqvist, D. W. Ryu

"Demand of high pressurized gas storage in underground rock caverns is increasing. Compressed natural gas (CNG) and compressed air energy storage (CAES) are one of the promising examples after since lined rock cavern (LRC) concept has been verified through in-situ pilot tests. Two key essential concerns of high pressured gas storage in underground rock caverns are mechanical stability of the cavern against high inner storage pressure and the risk of gas leakage from the caverns.In this study, we introduce our results of numerical simulations so as to investigate the potential of monitoring for gas leakage at a pilot facility for CAES in underground LRCs. The pilot facility was constructed at a comparatively shallow depth of 100 m, compared to the previous cases with several hundreds of meters, aiming at validating the feasibility of the LRC storage concept in which inner sealing liners, consisting of steel plate and concrete linings, are installed inside the excavated rock caverns (KIGAM, 2011). Compressed air may leak through either initial defects of inner containment liners such as imperfect welding regions and in-situ construction joints of the liners, or through structurally damaged points of the liners during the CAES operation where compression and decompression of the caverns are necessarily repeated and their structural robustness can be deteriorated.For this perspective, we carried out a coupled thermodynamic and geomechanical deformation analysis using TOUGH-FLAC simulator, and investigated the effectiveness of strain monitoring method in detecting mechanical failure of liners which may provide major pathways of air leakage. We therefore used an equivalent continuum model to relate the localized tangential strain to fracture permeability changes, which was developed in our early work (Rutqvist et al, 2012). In detecting and characterizing the leakage location, we tested a cross correlation technique using the synthetic pressure responses monitored in case of the failure of the inner liners. From the investigations, we showed that tangential strain, rather than radial, monitoring especially at the inner face of the sealing liners can be more effective in detecting the mechanical failure of the liners.In localizing the leakage point, we noted that the interface between the inner liners and surrounding rock may not be perfectly bonded so as to increase permeability through which any leaking air pressure may migrate more dominantly. We also demonstrated that the cross correlation of the time histories of pressure measurement especially at these interfaces can be useful in identifying the leakage point. These results may provide useful information in overall design of monitoring and alarming system of our upcoming CAES pilot facilities."

Jan 1, 2015

By Habib Zughbi, Paul Zulli

"A key step towards sustainability in iron and steelmaking is to optimise the use of raw materials and energy consumption, and consequently minimise CO2 emissions. Models based on a Process Integration (PI) approach play a key role in evaluating various scenarios and the impact of introducing changes in raw materials and operational conditions on the overall energy consumption and CO2 emissions. BlueScope uses a PI based comprehensive model for its Port Kembla Steelworks (PKSW), called the Integrated Steelworks Energy and Emissions Model (ISEEM). ISEEM is tightly coupled to the configuration of PKSW. It is based on mass and energy balances as well as thermodynamic relationships. ISEEM is a steady-state model and is used to predict overall and detailed unit-by-unit energy consumptions and CO2 emissions. It is also used for plant-based scenario modelling. This paper investigates the impact of major changes to PKSW configuration, including the shutting down of a blast furnace and cokemaking battery, on the energy consumption and CO2 emissions of the plant. It discusses some of the issues involved in updating the model and some identified opportunities."

Jan 1, 2015

By J. Lee, K. B. Min, J. Rutqvist

Shear slip of fractures in fractured rock masses is a critical issue for many geological engineering applications, such as the design of underground nuclear waste repositories, enhanced geothermal systems (EGS), and geosequestration of CO2. This paper presents new development of the TOUGH-UDEC simulator for analysis of coupled thermal-hydraulic-mechanical (THM) processes in fractured porous media. As an alternative to developing a single coupled THM code in discontinuum, two existing wellestablished codes, namely TOUGH2 and UDEC, are coupled to cover coupled THM analysis of multiphase fluid flow, heat transfer, and deformation in fractured porous media. TOUGH2 is widely used for modeling heat transfer and multiphase, multicomponent fluid flow, and UDEC is a well-known discrete element code for rock mechanics, which have been verified and used by many groups all over the world. Although a sequential coupling of two codes is less efficient than having a single code in principal, a big advantage with coupling of TOUGH2 and UDEC is that both codes are well tested and widely applied in their respective fields. When linking two codes, the coupled equations cannot be solved simultaneously; they should be solved sequentially, with coupling parameters passed to each equation at specific intervals. After solving thermal-hydraulic equations through the TOUGH2 model, pressure (P) and temperature (T) information is delivered to UDEC model. Pressure information is directly stored at the center of each UDEC zone and domain, however, temperature information is interpolated from zone-centers to the corners nodes. After solving the mechanical equation through the UDEC model—considering new pressure and temperature—calculated fracture apertures are converted to equivalent permeability and porosity for TOUGH2 flow analysis. The solution was calculated by iteratively following an explicit sequence for numerical efficiency. Verification examples are presented to demonstrate the capabilities of the coupled TOUGH-UDEC simulator.

Jan 1, 2015