scholarly journals Crustal Fault Zones (CFZ) as Geothermal Power Systems: A Preliminary 3D THM Model Constrained by a Multidisciplinary Approach

Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-24
Author(s):  
Hugo Duwiquet ◽  
Laurent Guillou-Frottier ◽  
Laurent Arbaret ◽  
Mathieu Bellanger ◽  
Théophile Guillon ◽  
...  

The Pontgibaud crustal fault zone (CFZ) in the French Massif Central provides an opportunity to evaluate the high-temperature geothermal potential of these naturally permeable zones. Previous 2D modeling of heat and mass transfer in a fault zone highlighted that a subvertical CFZ concentrates the highest temperature anomalies at shallow depths. By comparing the results of these large-scale 2D numerical models with field data, the depth of the 150°C isotherm was estimated to be at a depth of 2.5 km. However, these results did not consider 3D effects and interactions between fluids, deformation, and temperature. Here, field measurements are used to control the 3D geometry of the geological structures. New 2D (thin-section) and 3D (X-ray microtomography) observations point to a well-defined spatial propagation of fractures and voids, exhibiting the same fracture architecture at different scales (2.5 μm to 2 mm). Moreover, new measurements on porosity and permeability confirm that the highly fractured and altered samples are characterized by large permeability values, one of them reaching 10-12 m2. Based on a thermoporoelastic hypothesis, a preliminary 3D THM numerical model is presented. A first parametric study highlights the role of permeability, stress direction, and intensity on fluid flow. In particular, three different convective patterns have been identified (finger-like, blob-like, and double-like convective patterns). The results suggest that vertical deformation zones oriented at 30 and 70° with respect to the maximum horizontal stress direction would correspond to the potential target for high-temperature anomalies. Finally, a large-scale 3D numerical model of the Pontgibaud CFZ, based on THM coupling and the comparison with field data (temperature, heat flux, and electrical resistivity), allows us to explore the spatial geometry of the 150°C isotherm. Although simplified hypotheses have been used, 3D field data have been reproduced.

2019 ◽  
Vol 7 (1) ◽  
Author(s):  
Hugo Duwiquet ◽  
Laurent Arbaret ◽  
Laurent Guillou-Frottier ◽  
Michael J. Heap ◽  
Mathieu Bellanger

Abstract The present study aims to understand the potential of a new and novel type of geothermal play system for high temperature and electricity production: crustal fault zones (CFZ). According to geological and geophysical data, the Pontgibaud fault zone (French Massif Central) is suspected to host an active hydrothermal system at a depth of a few kilometers. The deep geometry of the fault zone and the permeability distribution are the main unknown parameters that are required to assess the geothermal potential of the Pontgibaud site. Structural and thin-section observations, laboratory permeability and connected porosity measurements and X-ray micro-tomography observations suggest that the hydrothermal system behaves like a double matrix-fracture permeability reservoir. Numerical modeling in which we varied the fault dip and the ratio between the fault zone permeability and host rock, R, was performed. Results indicate that three main convective regimes can be identified (weak convection, single cellular-type convection and bicellular convection). For a sufficiently high fault zone permeability (> 1 × 10−15 m2), buoyancy-driven flow creates a positive thermal anomaly of several tens of °C at a depth of 2–5 km. For a vertical fault zone, the thermal anomaly is larger for higher R values. Numerical models, then applied to the geologically constrained Pontgibaud fault zone, show that a temperature of 150 °C at a depth of 2500 m can be obtained for a fault zone permeability of 1.6 × 10−14 m2. Based on a multi-disciplinary approach, this work establishes a potential predictive tool for future high-temperature geothermal operations within basement rocks hosting large-scale fault systems.


Author(s):  
Qianyu Shi ◽  
Zhijian Wang ◽  
Hui Tang ◽  
Qi Li

Abstract Large scale molten salt storage tanks are widely used in the solar thermal power systems. For these tanks, buckling is a primary failure mode because of its features such as large scale, thinned wall and high temperature. Suffering high temperature condition is a major distinction between molten salt storage tanks and other water or oil tanks. High temperature can cause large thermal deformation for large scale structures which may have an effect on the safety assessment, especially on buckling assessment. Meanwhile, the hydrostatic pressure of molten salt can also cause the change of tank’s configuration. In this paper, a typical large molten salt storage tank has been studied. The critical buckling loads of the tank roof are obtained using nonlinear buckling analysis considering thermal loads and hydrostatic pressure. The results are discussed and some conclusions are proposed for engineering design.


2003 ◽  
Vol 17 (04n06) ◽  
pp. 528-533 ◽  
Author(s):  
RICCARDO TEBANO ◽  
RENATA MELE ◽  
VINCENZO BOFFA ◽  
FEDOR GÖMÖRY ◽  
FRANTISEK STRYCEK ◽  
...  

Reduction of AC losses for large-scale applications of superconductors is a critical issue. Therefore, the quantitative evaluation of AC losses is important for the development of superconductors and their applications to electrical power systems. The development of numerical models that simulate the electromagnetic phenomena inside superconductors allows to understand the electromagnetic behavior of superconductors and to evaluate the AC loss properties. Following an approach proposed by Brandt in several papers, a numerical model was developed in order to study the AC properties of superconducting tapes in different geometrical arrangements and with time dependent current and magnetic field. Here we show an example for simple single-layer model cables to show how this rather simple and versatile numerical approach allows optimizing configurations for actual operational conditions.


1986 ◽  
Vol 1 (20) ◽  
pp. 143
Author(s):  
H.E. Klatter ◽  
J.M.C. Dijkzeul ◽  
G. Hartsuiker ◽  
L. Bijlsma

This paper discusses the application of two-dimensional tidal models to the hydraulic research for the storm surge barrier in the Eastern Scheldt in the Netherlands. At the site of the barrier local energy losses dominate the flow. Three methods are discussed for dealing with these energy losses in a numerical model based on the long wave equations. The construction of the storm surge barrier provided extensive field data for various phases of the construction of the barrier and these field data are used as a test case for the computation at methods developed. One method is preferred since it gives good agreement between computations and field data. The two-dimensional flow patterns, the discharge and the head-difference agree well,, The results of scale model tests were also available for comparison. This comparison demonstrated that depth-averaged velocities, computed by a two-dimensional numerical model, are as accurate as values obtained from a large physical scale model. Even compicated flow patterns with local energy losses and sharp velocity gradients compared well.


2012 ◽  
Vol 1 (33) ◽  
pp. 64 ◽  
Author(s):  
Haiqing Liu Kaczkowski ◽  
Timothy W Kana

Nags Head, located at the northeastern part of North Carolina in the U.S., has sustained chronic erosion over the past 50 years. In 2005, Coastal Science & Engineering (CSE) was retained by the town of Nags Head to develop an interim beach restoration plan. Profile volume change was used in the planning and preliminary design of the project, and longshore and cross-shore numerical models were used in the final design to refine the preliminary nourishment plan and increase potential longevity of the project. This paper focuses on the key factors of the longshore numerical model setup for the project. These include model selection, input data and parameters, model calibration, and applications under different design alternatives. The Generalized Model for Simulating Shoreline Changes (GENESIS) was used in this study to evaluate shoreline evolution under normal wave conditions during various stages of the design life following the beach nourishment project. The model was used to identify the potential occurrence of erosional hotspots and to optimize the nourishment design so that the effects of such hotspots could be avoided or minimized where possible. Model results were also used to evaluate the impact of borrow area dredging on longshore transport in the project area and the impact of nourishment on shoaling in the adjacent inlet. The project encompasses 10.11 miles (mi) (16.28 kilometers-km) of ocean shoreline, and the design nourishment volume is based on the total permitted volume of 4 million cubic yards (cy) (3 million cubic meters-m³). [Note: As-built length was 10.0 mi and volume was 4.615 million cubic yards.] The final design has fill densities varying from north to south in relation to historical erosion rates and model projections. The average fill density is 75 cubic yards per foot (cy/ft) (188 m³/m) and ranges from 38 cy/ft to 150 cy/ft (95 m³/m to 375 m³/m). In conclusion, it is shown that the numerical model selected in this study was capable of predicting the overall performance of the large scale beach nourishment project in Nags Head as well as the performance at a particular location within or adjacent to the project, and its design methods can offer guidance to future projects.


2021 ◽  
Author(s):  
Tom Manzocchi

<p>Faults can control the large-scale properties of rock volumes through their behaviour as flow conduits and/or barriers or by localising geomechanical effects. Hence, often the fidelity of a numerical model of faulted site relies on the accuracy with which the fault zone is represented.  There are two distinct factors that must be considered in a modelling study: first, does the model contain the most relevant characteristics of the fault that influence the behaviour of interest; and second, are these characteristics assigned realistic and representative values that capture both their natural variability and the uncertainty with which they can be determined for the specific case of interest. These two factors are contained in the conceptual fault model and choice of modelling proxy-properties, respectively.</p><p>In recent years, two classes of conceptual fault zone model have dominated the description of fault zones, broadly characterised by either a continuous or a discrete approach. Continuous fault zone properties (e.g. fault core and damage zone thickness, displacement partitioning statistics) often show high variability which many modelling studies attempt to capture by running multiple model containing property values sampled from the distribution. Discrete descriptions focus on the presence of individual fault zone elements (e.g. shale smears, relay zones), and models guided by a discrete conceptual model attempt to place representative frequencies of elements. A single discrete model might contain the same property distributions as an ensemble of continuous models yet, because it contains a representative frequency of different elements, its behaviour might lie beyond the extreme behaviour of the continuous ensemble. Hence, the manner in which a geologist’s conceptual model is represented in a modeller’s numerical model can be hugely important for the outcome of the study, and it is in the interest of both modellers and geologists to ensure that they have a correct understanding of the other’s part of the process.</p>


2004 ◽  
Vol 175 (2) ◽  
pp. 147-156 ◽  
Author(s):  
Christian Le Carlier De Veslud ◽  
Paul Alexandre ◽  
Michel Cuney ◽  
Gilles Ruffet ◽  
Alain Cheilletz ◽  
...  

Abstract The eastern part of the “Seuil du Poitou” area had been selected by the French National radioactive waste management agency (ANDRA) as a potential site for building an underground laboratory in granitic rocks. 17 cored bore holes, completed by petrographical, geochemical [Cuney et al., 2001], geophysical [Virlogeux et al., 1999] and structural [Gros and Genter, 1999] studies, have provided a detailed knowledge of the Charroux-Civray complex, hidden under a Jurassic sedimentary cover. Three main types of magmas were distinguished : medium-K calk-alkaline, high-K calk-alkaline and peraluminous. The first two types are largely predominant and belong to the “Limousin Tonalitic line” (LTL) [Peiffer, 1985 ; 1986]. They were emplaced between 356 ± 5 Ma to 349 ± 5 Ma from U-Pb dating of zircon [Bertrand et al., 2001], at a depth of 14 ± 2 km [Freiberger et al., 2001]. This work aims to reconstruct the thermochronological evolution of the Charroux-Civray complex from 40Ar/39Ar dating of biotite and amphibole, combined with available U/Pb ages [Bertrand et al., 2001] and thermo-barometric data derived from mineral paragenesis and fluid inclusions [Freiberger et al., 2001]. This reconstruction may provide interesting constraints on a stage not well understood in the evolution of the French Massif Central : the emplacement of the LTL granitoids. The datings were performed on alteration-free, single grain of biotite and amphibole from the main petrologic types, according to the procedure described by Ruffet et al. [1991] and [1995]. The closure temperatures of the isotopic systems have been assumed to be 300 ± 30 oC for biotites, 500 ± 50 oC for amphiboles, and 850 ± 50 oC for zircons [Villa, 1998 ; McDougall and Harrison, 1999]. Six samples were dated : two tonalites (samples 112 and 212t), a monzogranite (sample 106), a monzogabbro-diorite (sample 115), a monzodiorite (sample 104), and a granodiorite (sample 105). Some of the analyses have been performed twice to test the reproducibility of the 40Ar/39Ar measurements. The 14 age spectra obtained may be divided into four groups : plateau ages, which provide robust ages for the amphiboles of samples 104, 106, 112, and biotites from samples 106, 115 and 212t ; pseudo-plateaux ages : three biotites (samples 104, 105 and 112) display spectrum shapes that could be interpreted as resulting from 39Ar recoil, related to an incipient chloritisation [Ruffet et al., 1991 ; McDougall and Harrison, 1999]. The most reliable ages are therefore close to the apparent ages given by the high temperature steps ; 40Ar* excess, as suggested by the spectrum shape of the amphibole from sample 212t [McDougall and Harrison, 1999]. The preferred age is defined with 83 % of released gas, and has been confirmed by a duplicate analysis ; a meaningless spectrum has been obtained on the amphiboles from sample 115. A duplicate analysis provided an approximate age of 347 ± 1 Ma, calculated on a relatively flat segment of the age spectrum. These results show that : (1) the closure of the isotopic system of the amphiboles occurred at approximately the same time over the entire complex (about 348 Ma) ; (2) the closure of the isotopic system of the biotites occurred slightly after the closure of the amphiboles, but spread over a larger time interval (350–343 Ma), (3) all the samples display high temperature gradients between 500 and 300 oC (> 40 oC.my-1). These results are in good agreement with mineralogical and fluid inclusion thermo-barometric data [Freiberger et al., 2001]. Two scénarios may be invoked to explain such high temperature gradients : a fast exhumation episode (several mm/y) during the 350–340 Ma period. This model is not acceptable because it is incompatible with pre- and post-intrusion conditions constrained by thermo-barometric data ; a fast thermal equilibration of the complex with surrounding rocks at the end of a succession of nearly-synchronous emplacement of calk-alkaline intrusions. First-order numerical models were used to simulate the thermal equilibration of the intrusive bodies with surrounding rocks, assuming a purely conductive heat regime [Carslaw and Jaeger, 1959]. These models show that according to the size of intrusions, the thermal equilibrium with surrounding rocks is reached in less than 5 to 10 m.y. The calculated temperature gradients derived from these models are compatible with those deduced from 40Ar/39Ar ages. These data confirm the existence of a major calk-alkaline magmatic event on the Seuil du Poitou, at about 355–350 Ma, which would be synchronous with the emplacement of the large peraluminous Guéret-type granodiorites in the northern Limousin. The 40Ar/39Ar biotite ages indicate that a regional temperature of 250–300 oC was reached at ca. 340 Ma at a depth of about 9 km.


2021 ◽  
Vol 1 ◽  
pp. 165-167
Author(s):  
Tuanny Cajuhi ◽  
Jobst Maßmann ◽  
Gesa Ziefle

Abstract. Salt, crystalline and clay formations are under discussion as potential host rocks for storage of heat-generating radioactive waste. Each of these rocks has a different structure and composition, and consequently a different material behavior. The latter needs to be studied and evaluated with respect to the main aim: to find a place to store the waste in a safe and sustainable manner. Several requirements in the context of the safety of a repository need to be fulfilled, concerning the long-term as well as the operational phase. One key point in this matter is the integrity, which refers to retention of the isolating rock zone's containment capabilities. With the focus on some experimental and numerical investigations on the excavation influenced near-field behavior of Opalinus clay (OPA), this contribution aims to illustrate an example for the role of numerical modeling in safety assessment. Once, e.g. anthropogenic action such as excavation starts, the natural state of equilibrium in the formation is disturbed. Trying to restore it, the rock deforms (convergence) and/or releases energy in other ways such as cracking. This could lead to loss of integrity since crack nucleation and propagation can affect the mechanical stability and create paths to transport contaminants. During operation in the excavated rock, environmental changes, e.g. temperature and humidity, further affect its behavior. The understanding of these dynamic phenomena ideally needs to occur at the in situ scale; however, performing an experiment in the spatial and time scales of interest is not always possible. For this reason, the in situ problem needs to be formulated, abstracted and mathematically modeled. The interpretation of the results must take place with simplifying assumptions and complementary laboratory scale experiments can be used to improve understanding of the system. The real problem is approached stepwise, each step associated to the size of the model and its complexity. The gradually obtained knowledge is necessary to achieve a better understanding of the process and to evaluate the capacities and limitations of the models. This contribution aims at showing the basic practical steps for numerical modeling with particular focus on the preparation and interpretation of the models and results, e.g. model calibration, verification and validation. As an example, the OPA at the Mont Terri site is chosen. The material parameters are obtained either experimentally or from the literature. We choose and perform laboratory scale simulations that are related to nearly the same mechanism as in the in situ scale. To have a first impression on the latter, a simplified, large-scale numerical model is prepared. The mechanism in study is drying and wetting, which is associated with shrinkage and swelling. We analyze the pore pressure and stress development in both scales. Thus, hydraulic mechanically coupled approaches are essential. The concept of effective stress is used, which combines the contributions of the solid and fluid phases (gas and liquid). In the current modeling approach, the gas pressure remains constant (atmospheric pressure) and during drying, the liquid pressure induces capillary pressure development and decrease of saturation. The laboratory scale simulation is important to evaluate the model of choice and to assess potential numerical problems. Furthermore, it can be used to perform a sensitivity study of material and numerical parameters. This step is necessary during the development or extension of numerical models as well as to evaluate their applicability on new research questions. The simplified in situ scale numerical model is then extended. In this phase the numerical model is evaluated once again, especially with respect to its complexity. Furthermore, specific questions related to this scale are posed: overall behavior of the rock, influence of the excavation, seasonal and long-term effects. In this contribution we deal with the long-term cyclic deformation (CD-A) experiment. The CD-A experiment has been taking place in the Mont Terri Rock Laboratory since October 2019. It consists of twin niches, a closed and an open niche, subjected to either high air humidity or seasonal humidity changes leading to saturation/desaturation during summer/winter in the OPA, respectively. Several parameters are periodically or continuously measured, including relative air humidity, convergence and crack development. We attempt to transfer the knowledge and numerical models developed in the small scale to the large scale and to evaluate the possibilities and limitations of the chosen approaches by comparing the numerical and experimental results.


To minimize groundwater pollution and suggest appropriate remedial actions, sound numerical models must be developed to predict the fate, transport and biodegradation of pollutants in partially treated or untreated industrial wastewater. Phenol is an aromatic organic compound produced on a large scale and is also released as major organic pollutant from several industrial wastewater (pharmaceutical, petroleum, coal refineries etc.,). Biodegradation of phenol in soil is generally rapid especially in presence of nutrients and acclimated microbes which are discharged along with the partially treated wastewater. A numerical model has been developed to predict the fate of phenol from industrial wastewater discharged on to a porous unsaturated soil media. The transport processes of advection, dispersion, and biodegradation process using Haldane growth and inhibition have been incorporated in the numerical model. The results suggest that acclimated microbes in the wastewater has a potential to degrade phenol up to 1500 mg/L at a bacterial concentration of 0.1 mg/L and soil depth of 50 cm. The results also show that phenols desorb at a depth of 100cm from 12th day and are simultaneously acted upon by the increased microbial concentration. In essence, high microbial concentration significantly decreases the phenol movement in the unsaturated zone, particularly at a larger depth and at higher time levels which eventually affects the groundwater quality.


2018 ◽  
Author(s):  
Matthias May ◽  
Kira Rehfeld

Greenhouse gas emissions must be cut to limit global warming to 1.5-2C above preindustrial levels. Yet the rate of decarbonisation is currently too low to achieve this. Policy-relevant scenarios therefore rely on the permanent removal of CO<sub>2</sub> from the atmosphere. However, none of the envisaged technologies has demonstrated scalability to the decarbonization targets for the year 2050. In this analysis, we show that artificial photosynthesis for CO<sub>2</sub> reduction may deliver an efficient large-scale carbon sink. This technology is mainly developed towards solar fuels and its potential for negative emissions has been largely overlooked. With high efficiency and low sensitivity to high temperature and illumination conditions, it could, if developed towards a mature technology, present a viable approach to fill the gap in the negative emissions budget.<br>


Sign in / Sign up

Export Citation Format

Share Document