ENHANCED GEOTHERMAL SYSTEMS (EGS) - NUMERICAL PREDICTION OF THE MODE AND LOCATION OF FRACTURE INITIATION

Mohamad Zaarour
2020
Student project

Abstract

Geo-energy is a comprehensive term used to describe any form of energy that comes from the Earth. This includes hydrocarbons such as gas, oil, and coal, but also geothermal energy (shallow and deep). The focus of this thesis is on Enhanced Geothermal Systems (EGS), which are renewable sources of energy. In EGS, hydraulic stimulation is the technique employed to enhance the permeability e.g. of the granite basements at great depths, which might exceed 5 kilometers. Rocks have natural fractures which affect the behavior of the rock mass. The mechanical goal of hydraulic stimulation is to open these pre-existing fractures. Due to the generally great depths, only indirect information such as in-situ stress, pumping records, and micro-seismicity, are available. Many issues in the field of reservoir geothermal (but also petroleum) production, are attributed to the lack of proper understanding of how these fractures interact with each other. The main goal of this study is to provide an insight into how the pre-existing fractures in rocks interact with each other. Specifically, this research answers the question of how pre-existing flaws of different geometrical parameters contained in rock specimens subjected to various loading configurations, affect the mode and the location of the mechanism of crack initiation. We develop a finite element (FE) linearly elastic MATLAB Partial Differential Equation (PDE) code. This numerical model predicts the mode and the location of the initiation of the local fractures at the inner tips of two interacting pre-existing flaws in Barre granite in the context of EGS. We investigate different combinations of double-flaw geometries (flaw inclination angle β, bridging angle α, and ligament length L) contained in specimens subjected to various loading scenarios (Combination of Vertical Load VL, Lateral Load LL, and Internal Flaw Pressure FP) likely to induce first local tensile and shear fractures, and we determine their corresponding predicted locations around the inner flaw tips. For a given double-flaw pair geometry, we subject the specimen to a certain loading configuration (Combination of VL, LL, and FP). We then retrieve the state of stress (maximum and minimum principal stresses 𝜎1 and 𝜎3, and maximum shear stresses 𝜏𝑚𝑎𝑥) at the Finite Element (FE) corner nodes on the boundary of the inner tips, from which we plot the corresponding Mohr’s circles. We study the interaction of these Mohr’s circles with the Griffith-Coulomb failure envelope. Since the main objective is to predict the first local failure’s mechanism and location, we iteratively vary the magnitude of the loading component until the first Mohr’s circle touches the failure envelope. Depending on where it does, the mechanism of the local fracture is inferred. If the maximum principal stress (𝜎1) of this Mohr’s circle is equal to the tensile strength of the granite (𝜎𝑇=7.5𝑀𝑃𝑎), a local tensile fracture is numerically predicted. If this Mohr’s circle touches the linear shear Coulomb portion of the envelope, it is a predicted local shear failure. The corresponding location of the FE node is then determined at the inner tip of the flaw, and this is where for the given geometry and loading scenarios, the first local fracture’s initiation is numerically expected. We follow this procedure for several double-flaw geometrical configurations. We compare the numerical predictions with the previously obtained experimental results in the MIT CEE Rock Mechanics Laboratory. This allows us to validate our set-up numerical model. 4 This study provides a solid ground to the understanding of the initiation and the interaction of pre-existing fractures in EGS. This thesis comes as an answer to a problem of major practical importance in the field of rock mechanics, and beyond. Typical geotechnical problems, such as rock slope stability and tunnel support, can also be better addressed when clearly understanding the process of crack initiation in rocks. This study also enables the correlation between numerical models and experimental findings, a challenge that has always existed in the research world. This thesis shows how our developed FE code can be used as a tool to predict the location and the mode of fractures in the laboratory, but also in the real world.

Official source

https://infoscience.epfl.ch/record/280102?ln=en

About this result

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related concepts

Related publications

Mohamad Zaarour
2020
Student project

Abstract

Official source

https://infoscience.epfl.ch/record/280102?ln=en

About this result

Related concepts (20)

Related concepts

Related publications (11)

Related publications

Typologie des eaux souterraines de la molasse entre Chambéry et Linz (France, Suisse, Allemagne, Autriche)

The geochemical typology of natural subsurface waters of the Molasse basin is part of the AQUITYP project, initiated by the Geology Laboratory (GEOLEP) of the Swiss Federal Institute of Technology in Lausanne (EPFL). This study provides an overview of the variation in space of the hydrochemical parameters, especially of trace elements, in order to determine the geological origin of their aquifers. Our data base is mainly made up of two observation networks: a subsurface water network and a deep groundwater network. Complete analyses of all samples have been executed using standard methods for anions and silicium, and high resolution inductively coupled mass spectrometry (HR-ICP-MS) for major cations and trace elements down to the detection limit of about 0.2 µg/l. Subsurface water network It is the main observation network of this investigation and comprises 112 water samples from different springs, collected all over the Molasse basin between Chambéry (France) and Linz (Austria), but mainly in the Swiss Molasse basin. For this region, four groups are generally distinguished, according to the sediment deposition in a continental or marine environment, in ascending stratigraphic order: (1) the Lower Marine Molasse (UMM and North Helvetic "flysch" p.p.), (2) the Lower Freshwater Molasse (USM), (3) the Upper Marine Molasse (OMM) and (4) the Upper Freshwater Molasse (OSM). Our water samples originate from typical sedimentary deposits which we divided in about 30 different sub-type aquifers according to their geological origin and their depositional environment. They are defined on the basis of lithological and structural homogeneity. Therefore, only catchments that completely drain a single, lithologically homogenous, Molasse terrain unit are considered here. Spring catchments selected were generally not affected by runoff from towns, or paved roads. Agricultural activities have been avoided as well as possible. Such spring water should be free of major human influences and only in few cases contain a certain amount of fertilizer products (p.ex: nitrates). Leachate tests of crushed rocks and comparison with snow and soil water compositions allowed to decipher the geogene origin of the dissolved ions in natural groundwaters. The water withdrawn can be generally classified as a calcium-magnesium-hydrogenocarbonate (Ca-Mg-HCO3)-type water. The trace analyses allowed to distinguish several subtypes within this homogenous water type group. The general spacial distribution of the trace elements show, that spring waters in the Bavarian Molasse basin are caracterized by elements as lithium and uranium linked to acid rocks, while those of the Swiss Molasse basin are preferentially influenced by elements derived from mafic or ultramafic rocks, such as chromium and cobalt. The following elements show specific links between aquifer rock type and groundwater. Sulfates: Total mineralisation of spring waters issued from the gypse-bearing Molasse (USM) and the "Glimmersande" (OSM) are highly influenced by high sulfate contents. It can be demonstrated that sulfate ions are produced from two different sources within this two sub-type aquifers: from the dissolution of gypsum in water, producing high calcium contents as well, and from the oxydation of sulfurous ore minerals, especially pyrite, which are abundant in the "Glimmersande". The "Glimmersand" aquifers in NE Switzerland are furtheron caracterized by high lithium, molybdenum and uranium concentrations, reflecting the granitic composition of this sandstones. High barium contents (>150 µg/l) have been found in spring waters issued from aquifers with a dominant fissure permeability ("Grès de la Cornalle", "Gibloux-Delta" and "Rigi-Schüttung"). Its origin is supposed to be the mineral barite, which tends to precipitate in fractures. Therefore, high barium concentrations may indicate the presence of fracture flow within the aquifer. Considerably high silicon contents are linked to the quartzitic and fine grained aquifers of the brackish marine molasse in Bavaria. Both, petrographical composition and slow flow rate are favouring the water-rock interaction. Caracteristic lithium footprints are found in several waters, as in the aquifers composed of the "Glimmersande", the gypse-bearing molasse, the Rigi-conglomerates, the brackish and the estuarine marine molasse. Chromium strongly marks the waters originating from the Swiss marine Molasse aquifers. It is probably related to ophiolithic detritus derived from the flysch of the Prealpes. A couple of waters show weak anomalies of cobalt, which are nevertheless intimately linked to the presence of ophiolithic components in the sediments. Cobalt might be used as a hydrostratigraphical tool to identify greenstone inputs. Analysed colloidal aluminum may serve as a hydrodynamical parameter, indicating microturbidity within the aquifer. High amounts of aluminum varying in time, have been found in estuarin deposits of the OMM. Generally speaking, we propose three main origins of trace elements in the spring water of the Molasse basin: presence of certain petrographical components within the sediments such as basic and ultrabasic rocks, heavy minerals and evaporites. a small flow rate through interstitial porosity of fine grained sediments, favouring water-rock interaction. the abundance of colloidal material in the waters which depends on the petrographical composition of the aquifer and on sample treatment, e.g. filtration techniques. Deep groundwater network A second observation network composed of 21 borehole waters has been added for comparison. Compared to the subsurface waters, deep water chemistry might differ. The following reasons seem to be responsible: Long residence times of water in the underground favours slow reaction kinetics. Chemistry of very deep water can be influenced by mixing with older formation waters. Most of the deep boreholes in the Swiss Molasse basin are situated in the OMM and therefore represent a likely homogenous lithology. Borehole waters, because of different aquifer rock types and unknown surface and borehole influences, are least likely to represent typical Molasse water composition.

EPFL1995

Nuclear materials have been successfully employed in commercial fission reactors for many decades. Research and material development for advanced fission system and future thermo-nuclear fusion reactors have acquired maturity while great challenges remain to be taken up. The goal of the material scientist is to understand the way neutron irradiation and aggressive reactor environment alter the original material properties and to safely manage the operation conditions for each particular reactor component. In the framework of surveillance programme, the need has arisen to design experimental methods to be used with a limited amount of material for Post Irradiation Examination (PIE) of reactor components to test the integrity of those on a regular time basis. Having to deal with radioactive material, it is particularly important to devise experimental testing procedures as simple as possible but from which reliable information can be drawn. Within materials development programme, irradiations have continuously been performed either in reactors or with particle accelerators. In both cases, the irradiation volumes are usually limited. The concern was to optimize the volume of the irradiated samples in order to increase their number to obtain a good statistical representation of their mechanical behavior after irradiation, and to have a very definite and homogeneous damage dose and irradiation temperature. The concerns explained above originated the development of small-scale specimen testing methods. Advantage is taken of the fact that, within size and geometry limitations, some mechanical properties of a sample are independent of its absolute size. In this way it is possible to perform mechanical tests on samples of reduced volume that represents big advantages to face the aforementioned challenges. Non-standard test on very small specimen have also been proposed to extract mechanical properties of the bulk. If these non-standard testing methods are to be employed for PIE, it becomes necessary to develop modeling tool to understand under what circumstances they can provide useful information and specially how to obtain it. Finite element modeling (FEM) of the non-standard tests has played a key role in this sense, providing a link between standard test analysis and that of the non-standard ones. The goal of this work was to implement finite element models for non-standard tests to ultimately study the mechanical properties of Zircaloy cladding tube and tempered martensitic steels in the unirradiated and irradiated condition. One model was developed for ring tensile tests. It was first validated with a quite ductile iron-chromium ferritic alloy, for which ring and uniaxial tensile test were performed at room temperature. Using the constitutive behavior determined from the uniaxial tensile tests, the load – displacement curves of the ring tensile tests were very well reproduced. Another model was developed for the so-called small ball punch test, which consists in deforming a 3 mm diameter disk clamped between two dies with a ball puncher and to record the load – displacement curve. In order to validate the model, a series of ball punch tests and of tensile tests were carried out on the austenitic steel 316L at room temperature. In this last case, the ball punch test curves were also well reconstructed with the model. Based upon a simple analysis and simulation results, we showed that it is possible to obtain, by ring tensile tests, a quick and accurate assessment of the constitutive behavior of an unknown material, simply by processing the experimental force-elongation data in the same way as for uniaxial tests, taking half the ring perimeter as initial length and twice the ring cross section as the sample cross section. The ring tensile test method was used to determine the hoop mechanical properties of Zircaloy tubes used in fission reactors. Tensile tests were also performed in the axial direction of the Zircaloy tubes. Differences between the yield stress as well as the post-yield behaviors in the axial direction and those in the circumferential direction of the tube could be measured reflecting some component of the anisotropy of the tube. Ball punch specimens were cut also out of the tube in the radial direction and tested. Using the FE model for punch test, in which the anisotropy of the mechanical properties are taken into account through the Hill's theory, the values of the Hill's yield function parameters were estimated. The irradiation hardening of neutron irradiated Zircaloy tube was finally determined. The ball punch test model was extensively used for numerical investigations, to determine how the constitutive behavior of the materials (uniaxial true stress – true strain) employed as input acts in combination with other experimental factors (specimen thickness, friction coefficient) in mediating the shape of the experimental curves. In the light of the results, we derived calibrations between parameters of the non-standard and standard tests. For instance we proposed a new approach to determine the yield stress from the ball punch test. We also showed how the slope of the punch tests scales with the average flow stress of the material. The calibrations proposed in the past by authors and published in the open literature have been revisited and critically reviewed. Based upon the findings of the numerical studies, we propose a given methodology to unambiguously evaluate the material constitutive behavior of a material in a case where the only available data comes from various non-standard tests. Experimental ball punch test results obtained on proton-irradiated specimen (0.5 dpa and Tirr = 523 K), of the tempered martensitic steel, EUROFER97, have been obtained. The irradiation hardening was found equal to about 150 MPa, consistently with previous data on similar steels and irradiation conditions. From the simulations of the ball punch test curves of the irradiated specimens, it was concluded that the irradiation does not affect the strain hardening at relatively low fluence.

EPFL2005

Theoretical and Experimental Study of a Thermally Driven Heat Pump Based on a Double Organic Rankine Cycle

Jonathan Demierre

Nowadays, one of the main alternatives for a more rational use of energy in heating applications is the heat pumping technologies. The market is dominated by two kinds of heat pump systems: the electrically driven vapor compression heat pumps (EHP), which are the most widely used in residential heating applications and the thermally driven heat pumps (TDHP), that are usually based on a sorption process. In this thesis, theoretical and experimental investigations of a concept of thermally driven heat pump, based on the coupling of a vapor compression heat pump cycle and an organic Rankine cycle, are presented. The system is named here as ORC-ORC. The studied concept uses a single stage centrifugal compressor directly coupled to a single stage radial inflow turbine. The shaft is rotating on gas bearings, which allows the system to be oil-free. Like most of the other TDHP's, this system has the advantage to work with a variety of fuels or heat sources like wood pellets, natural gas, solar heat, geothermal heat or waste heat. The concept studied in this work is a gas red system for space heating and domestic hot water production in small residential buildings. A systematic approach has been used to evaluate, in term of energy efficiency, the potential of ORC-ORC systems and to propose optimal design specifications at different levels of the design process. The method is based on the optimization which allows to identify the best configurations at each design step with respect to the designer choices. This approach is divided in three steps. In the first step, a model of the complete system has been developed based on a process integration approach. This step allows to quickly determine whether the system is potentially attractive or not, for given conditions, before going deeper in the design process. The results show that, in general, it is advantageous for the ORC evaporation to be supercritical. In the second step, based on the process integration results and heuristic rules, a suitable system heat exchanger network has been generated and modeled. The predicted values of the coefficient of performance (COP) at design conditions, for the different situations, are in the range 1.30 and 2.19. In a last step, the off-design characteristics of the compressor-turbine unit are considered, in order to evaluate the COP as well as the heating capacity range that can be covered with a specific compressor-turbine unit design. For this purpose, a model of the compressor-turbine unit has been developed. Experimental investigations have been carried out with R134a as working fluid. A first simple test rig has been developed to test the selected ORC evaporator in supercritical conditions. The measurements have allowed to calibrate and to validate an in-house supercritical evaporator simulation tool. An ORC-ORC prototype has been developed and built. Commercially available equipment has been used, except for the compressor-turbine unit that has been designed especially for this application. The targeted operating point, with an HP evaporation temperature of 0 °C and a condensing temperature of 35 °C, has not been reached, because of experimental difficulties which caused some delay. A maximum temperature difference of 20 °C between the HP evaporation and the condensation has been achieved, which corresponds to a compressor pressure ratio of 1.9. Measurements have been performed with compressor-turbine unit rotational speeds up to 171 krpm. The theoretical investigations shows that the proposed ORC-ORC concept is an interesting alternative of thermally actuated heating device for small residential buildings. Concerning the experimental investigations, although a number of problems have been encountered during the tests, it has been demonstrated that the proposed concept is feasible with today's technology.