ERE5.4 | Coupled Thermo-Hydro-Mechanical-Chemical (THMC) processes in geo-reservoirs with emphasis on clay behavior
EDI
Coupled Thermo-Hydro-Mechanical-Chemical (THMC) processes in geo-reservoirs with emphasis on clay behavior
Co-organized by EMRP1
Convener: Silvia De SimoneECSECS | Co-conveners: Reza TaherdangkooECSECS, Emmanuel Gaucher, Iman Rahimzadeh KiviECSECS, Christoph Butscher, Alexandros DaniilidisECSECS, Keita Yoshioka
Orals
| Mon, 15 Apr, 16:15–18:00 (CEST)
 
Room 0.51
Posters on site
| Attendance Mon, 15 Apr, 10:45–12:30 (CEST) | Display Mon, 15 Apr, 08:30–12:30
 
Hall X4
Orals |
Mon, 16:15
Mon, 10:45
Geological media are a strategic resource for the forthcoming energy transition and constitute an important ally in the fight to mitigate the adverse effects of climate change. Several energy and environmental processes in the subsurface involve fluid circulation into geo-reservoirs, which activates a series of multi-physical interactions in the porous and fractured rock. Changes in pore pressure, stress and temperature, rock deformation and chemical reactions occur simultaneously and impact each other. Forecasts are bounded to the adequate understanding of field data associated with thermo-hydro-mechanical-chemical (THMC) processes and predictive capabilities heavily rely on the quality of the integration between the input data (laboratory and field evidence) and the mathematical models describing the evolution of the multi-physical systems. The relationship between the observations and the modelled changes is, however, often ubiquitous, which challenges the interpretation of the observations in regard of the physical processes in play within the reservoir. Certain geological settings involve an additional complexity, especially when clay materials are present, which have a versatile role acting as both assets and challenges in energy extraction.
This session is dedicated to studies investigating all or part of these THMC interactions by means of experimental, analytical, numerical, multi-scale, data-driven and artificial intelligence methods. Studies focusing on laboratory characterization and on monitoring and interpreting in-situ geological and geophysical evidence to validate or calibrate the multi-physical behavior of rocks and clays are also welcome. Applications in carbon capture and storage (CCS), radioactive waste storage, gas storage, energy storage, mining, geothermal systems, reservoir monitoring and management and geotechnical applications would be relevant.

Session assets

Orals: Mon, 15 Apr | Room 0.51

Chairpersons: Silvia De Simone, Emmanuel Gaucher, Keita Yoshioka
16:15–16:20
16:20–16:40
|
EGU24-8330
|
solicited
|
Virtual presentation
Maartje Boon

Large scale subsurface gas storage in porous reservoirs can play an important role in the energy transition. Geological storage of carbon dioxide will mitigate CO2 emissions while underground energy storage, for example in the form of hydrogen gas, can be used to balance out the renewable energy production and demand. To investigate the feasibility of large scale subsurface gas storage in porous reservoirs, simulation models are needed that accurately capture the multi-phase flow behaviour in porous rock. Important input parameters for reservoir simulators are relative permeability and capillary pressure which highly depend on the wettability of the system. In this presentation, we will show results of different experimental techniques to characterize and visualize gas transport in porous rock including a novel experimental device to characterize wettability under the impact of different driving forces.

How to cite: Boon, M.: Experimental characterization of multi-phase flow in porous rock relevant for subsurface gas storage, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8330, https://doi.org/10.5194/egusphere-egu24-8330, 2024.

16:40–16:50
|
EGU24-15531
|
ECS
|
On-site presentation
Francesco Lazari, Gabriel Meyer, and Marie Violay

In geothermal reservoirs, pore fluid chemistry can affect mechanical and hydraulic properties of rocks inducing mineral dissolution, precipitation, weakening and alteration. Moreover, with increasing pressure and temperature deformation mode can change from localized to ductile (diffused), leading to a major decrease in permeability, thus affecting the exploitability of the reservoir. The effect of fluid chemistry on the transition between localized and ductile deformation of rocks is still marginally understood.

To investigate the effect of water presence and fluid chemistry on localized and ductile deformation, two sets of triaxial experiments (at 20 and 100 MPa effective confinement pressure, in the localized and ductile field respectively) were performed on a porous silicate sandstone (Adamswiller sandstone), dry, with deionized water and with a 6 M NaCl solution (Na+ rich solution), with a 0.1 M HCl solution (pH 1 solution) and a 0.1 M NaOH solution (pH 13 solution). To complement the mechanical properties, complex spectral electrical conductivity was measured during deformation to monitor pore fluid ion content; pore fluid was collected at the beginning and at the end of the experiments and analyzed with ICP-MS; post-mortem microstructural analyses were performed.

In the localized domain, both water presence and pore fluid chemistry had a marginal effect on the strength of the rock, leading to a 5/10% strength reduction over dry rock strength indipendently of the pore fluid composition. In the ductile domain, de-ionized water weakens the rock by 25%, a Na+-rich or pH 1 fluid leads to a 35% weakening and a pH 13 fluid weakens the rock by 40% over dry rock strength. Spectral electrical conductivity does not change during localized deformation, while it increases by 2 to 5 times during ductile deformation when the rock is saturated with deionized water; conductivity does not change with the Na+-rich fluid regardless of the deformation mode and conductivity decreases by half an order of magnitude with both pH 1 and pH 13 fluids both with localized and ductile deformation.

 
 

How to cite: Lazari, F., Meyer, G., and Violay, M.: Effects of pore fluid chemistry  on the localized to ductile transition of sandstone., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-15531, https://doi.org/10.5194/egusphere-egu24-15531, 2024.

16:50–17:00
|
EGU24-19876
|
ECS
|
On-site presentation
Hadrien Rattez, Alexandre Sac-Morane, and Manolis Veveakis

Many natural processes and energy-related project in the subsurface involve the interactions between the porous rock and the fluids inside the pores. In particular, the fluid can interact with the solid matter by various dissolution-precipitation reactions that can modify the microstructure geometry and thus change the mechanical behavior of the rock and its failure potential or directly induce a creep of the rock. These chemo-mechanical couplings can have important implications for storage applications and induced seismicity. In this contribution, we will show discrete element simulations performed at the microstructural scale of porous matter. First, we will show that the dissolution of the cement in sedimentary rocks strongly influence the lateral earth pressure coefficient. The value of this coefficient tends to an attractor by increasing the degree of dissolution, which can lead to stress redistribution at the reservoir scale and promote faulting or induced seismicity. Secondly, we will show a numerical framework coupling discrete element with a phase field model allowing to capture grains shape changes due to local precipitation or dissolution. This model is applied to study the phenomenon of intergranular pressure-solution and allows to reproduce the creep behavior of the material in compaction. It enables also to study the competition between grain rearrangement and pressure solution in fault gouges to induce a rate dependency of the fault mechanical behavior.

How to cite: Rattez, H., Sac-Morane, A., and Veveakis, M.: Weathering and creep in rocks modelled at the microscale, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-19876, https://doi.org/10.5194/egusphere-egu24-19876, 2024.

17:00–17:10
|
EGU24-17626
|
ECS
|
On-site presentation
Lior Suchoy, Adriana Paluszny, and Robert W Zimmerman

The utilisation of geothermal heat as a form of clean energy is experiencing global growth. Many sites designated for geothermal energy extraction are in regions with elevated heat gradients, such as Iceland and Japan. However, the repercussions of the migration of cold fluids in rock formations at high temperatures, and the subsequent fracturing of the host rock, which is relevant in the context of deep geothermal energy systems, remains insufficiently understood.  In this study, we explore the three-dimensional evolution of fractures in 50×9.8×1 mm3  homogeneous slabs of various brittle rocks, subjected to a thermal shock of ΔT=580°C, through numerical simulations over 10 seconds. Initially, we validate our numerical approach using a benchmark of Al2O3  slab, and subsequently, we examine fracture development in granite, basalt, and shales. Our numerical methodology employs a three-dimensional finite-element-based simulator to model thermo-mechanical deformation. The in-house code is a fully coupled THM code which considers damage to predict fracture initiation. Fracture growth is predicted per fracture tip using stress intensity factors. The code uses adaptive meshing and NURBS for fracture surfaces to facilitate mesh-independent fracture growth. In our simulations, we apply triangles and tetrahedra elements to discretise surface and volume elements, respectively. Our results show the development of dozens of fractures in bi- or tri-modality, which penetrate up to 85% of the depth of the slab, for the various simulations. These results demonstrate both qualitative and quantitative agreement between the simulated slab and the benchmark by reproducing the same intertwined short-long fracture patterns and modal distribution of fracture lengths. Furthermore, they illustrate how the fracturing rates (ranging between 1-100 mm/sec), fracture length distribution (unimodal, bimodal or trimodal), and penetration depth of fractures in front of the shock front vary among the different brittle rock types.

How to cite: Suchoy, L., Paluszny, A., and Zimmerman, R. W.: Fracture growth using fully coupled thermo-mechanical model in brittle rocks during thermal shock and resulting network patterns, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17626, https://doi.org/10.5194/egusphere-egu24-17626, 2024.

17:10–17:20
|
EGU24-16640
|
ECS
|
On-site presentation
Thanushika Gunatilake, Antonio Pio Rinaldi, Alba Zappone, Yingqi Zhang, Dominik Zbinden, Marco Mazzotti, and Stefan Wiemer

The escalating global temperature necessitates immediate action to mitigate greenhouse gas emissions. Geological Carbon Storage (GCS) technology has emerged as a promising solution, specifically by injecting carbon dioxide (CO2) into geological formations, particularly deep saline aquifers. However, finding ideal geological reservoir conditions, including caprock stability and storage capacity, is a rare occurrence.

This study comprehensively assesses the potential for CO2 storage in the Triemli saline aquifer in Zurich, Switzerland. The goal is not only to demonstrate the feasibility of CO2 storage in Switzerland but also to emphasize the viability of decentralized storage with multiple small injection point, for emitters like medium -sized citiesin regions with geologically challenging subsurface conditions. Through numerical simulations, we explore CO$_2$ injection, migration, and long-term reservoir stability to bridge the gap between theoretical estimates and practical feasibility.  Our findings underscore the potential of deep saline aquifers for CO2 storage in Switzerland, particularly in the Swiss Molasse Basin and the adjacent Folded Jura, identified as crucial regions for effective CO2 storage. Employing advanced methods and strategic injection techniques, such as multiple vertical or horizontal injection points along a single well, could optimize this storage capacity to approximately 3 million tons of CO2 over the same period.

How to cite: Gunatilake, T., Pio Rinaldi, A., Zappone, A., Zhang, Y., Zbinden, D., Mazzotti, M., and Wiemer, S.: Decentralized CO2 storage in unfavorable conditions: An example from Switzerland, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16640, https://doi.org/10.5194/egusphere-egu24-16640, 2024.

17:20–17:30
|
EGU24-2591
|
On-site presentation
Peter Fokker, Loes Buijze, Maarten Pluymaekers, Chris Maaijwee, Harmen Mijnlieff, Job Mos, Bouko Vogelaar, Sjoukje de Vries, and Mark Vrijlandt

The safe and effective deployment of geothermal energy and storage of carbon dioxide requires an assessment of potential induced seismicity and fracturing through the seal above the reservoir. To aid such assessment we have built the SRIMA tool (Seal and Reservoir Integrity through Mechanical Analysis) and we have made it available online. The tool can be used in the Standard extended Seismic Hazard Analysis, which is part of the seismic hazard and risk assessment for geothermal projects in the Netherlands. SRIMA is a fast semi-analytical tool that provides a scenario-based analysis of pressure and temperature changes around an injection well, the resulting stress changes on nearby faults, reactivated fault area, the maximum credible earthquake magnitude, the resulting PGV distribution and an estimate of damage. SRIMA also computes the potential for development of tensile fracture in the seal and base. SRIMA has been designed to give first-order estimates of these results. The speed of the calculations facilitate them to be performed in a stochastic framework, which allows the assessment of failure probabilities.

 

SRIMA is based on semi-analytical expressions for the fast calculation of temperatures, pressures, and induced poro-elastic and thermo-elastic stresses due to the injection of cold fluid. The expressions for flow and induced stresses have been developed for a homogeneous, isotropic layer cake model under radial symmetry. In the injection layer the flow is assumed to be fully developed and temperature transfer is in an advective way. In the bounding seal and base layers, the pressure and temperature dynamics are assumed diffusive. The derived expressions capture the first-order characteristics of the pressure, temperature and stress changes. Validation of the expressions has been achieved through comparison with finite-difference and finite-element codes for temperature, pressure, and stress changes around an injection well. A fault without offset cutting through the seal, reservoir, and base can be specified within the model space. Poro-elastic and thermo-elastic stress changes are transformed to fault stresses and fault criticality. The fault area over which stresses are critical (i.e. fault reactivation occurs) is used to estimate the magnitude of the largest credible earthquake for each model scenario, assuming that the entire reactivated fault area participates in a single event, slip cannot propagate beyond the reactivated area, and all assumed slip over the fault area is seismic slip. An ensemble of magnitudes is converted to exceedance curves of peak ground velocities (PGV), using nationwide developed Ground Motion Prediction Equations. The resulting PGV distribution as a function of epicentral distance serves as input for calculating the probability of exceeding Damage State 1 using empirical fragility functions for unreinforced masonry buildings. This contribution will summarize details and assumptions behind each step.

How to cite: Fokker, P., Buijze, L., Pluymaekers, M., Maaijwee, C., Mijnlieff, H., Mos, J., Vogelaar, B., de Vries, S., and Vrijlandt, M.: SRIMA: A fast tool to assess seismicity and seal integrity related to fluid injection., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2591, https://doi.org/10.5194/egusphere-egu24-2591, 2024.

17:30–17:40
|
EGU24-11143
|
ECS
|
On-site presentation
Clara Fraile, Thomas Kohl, and Emmanuel Gaucher

In Central Europe, a substantial emitter of CO2 in the energy sector corresponds to the thermal energy required for heating and cooling. Seasonal underground heat storage presents a viable option for storing excess heat generated during the summer months for usage in winter, reducing the need for conventional sources of energy. Today, high-temperature aquifer thermal energy storage (HT-ATES) systems are attracting large interest as they represent a sustainable means of meeting heat demand.

In HT-ATES systems, hot water is injected into a reservoir during summer, while exchanged cold water is injected over the winter season. These fluctuations in temperature and pressure have an impact on the geomechanical and thermo-hydraulic properties of both the reservoir and the surrounding layers. Monitoring the changes in the reservoir properties is a critical aspect of running a heat storage system safely and efficiently. We try to determine whether active seismic imaging could be a suitable method to characterize the temporal and spatial evolution of the reservoir.

With view on designing future geophysical assessment and monitoring systems, we first perform thermo-hydro-mechanical (THM) modelling to estimate the variations in the poroelastic properties due to the geothermal processes. Our modeling is based on the characteristics of the DeepStor demonstrator, currently under development in the north of Karlsruhe (Germany), at the Karlsruhe Institute of Technology (KIT). The three layers model includes different mechanical properties with one borehole. The simulation of cyclic hot water injection and production over time allows to quantify its effect on the underground material properties. In addition to assessing the expected operational parameters of the DeepStor demonstrator, we test additional injection schemes with varying underground properties to simulate the different ranges of porosity changes and look at their effects on the elastic properties.

Linking the THM model parameters to seismic sensitive variables such as velocities and impedances, through empirical equations, allow us to determine the conditions under which active seismic surveys could effectively detect these changes. This approach provides a valuable tool for evaluating the potential of active seismic imaging for monitoring HT-ATES systems.

How to cite: Fraile, C., Kohl, T., and Gaucher, E.: Estimation of seismic velocity changes in a HT-ATES system using THM modelling, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11143, https://doi.org/10.5194/egusphere-egu24-11143, 2024.

17:40–17:50
|
EGU24-16823
|
ECS
|
On-site presentation
Gilles Corman and Frédéric Collin

In the field of radioactive waste confinement, the question of gas transfers in clay formations is a crucial issue. A certain amount of gas, such as Hydrogen may be generated during the exploitation phase in the nearfield by the deterioration of the metal components of the system. As the host medium is characterised by a very low permeability, various gas transport processes could occur as a function of gas accumulation and pressurization [1], including the development of preferential gas pathways through the sound rock mass, which could lead to undesirable changes in the favourable containment properties of the host rock.

There is a growing body of experimental evidences [2, 3] that natural heterogeneities and pre-sxisting fractures in clay-rich materials represent preferred weaknesses for the process of opening discrete gas-filled pathways. Capturing the related transport mechanisms therefore requires to go from macroscopic to microscopic scale. Hence, a multi-scale modelling approach that models the micro-scale effects explicitely on their specific length scale and couples their homogenized effects to the macro-scale is proposed in the present work. Based on a periodicity assumption of the microstructure, a relevant Representative Element Volume (REV) is defined based on experimental data, which makes it possible to idealise the flow behaviour of the material microstructure with different families of discontinuities, and an assembly of tubes substituting the porous matrix blocks. This complete hydraulic constitutive model is solved at the scale of the microstructural constituents, and is directly affected by the mechanical effects tackled at the macroscopic scale, which makes the whole model hydro-mechanically coupled.

This model has been subsequently applied to simulate a gas injection test parallel and perpendicular to the bedding of initially saturated samples of Boom Clay [3]. This analysis provides a rather good agreement with the experimental results in terms of pressure response, outflow volume and average axial strain. In addition, it allows to simulated the creation of a preferential flow pathway along the sample axis (Figure 1b, top), which serves as basis to numerically reproduce the development of random pathways through the sample in plane strain state (Figure 1b, bottom), and aims to improve the mechanistic understanding of the gas transport processes at play in clayey barriers.

[1] P. Marschall, S. Horseman, and T. Gimmi. Characterisation of Gas Transport Properties of the Opalinus Clay, a Potential Host Rock Formation for Radioactive Waste Disposal. Oil & Gas Science and Technology Rev. IFP, 60(1):121-139, 2005. 

[2] Harrington, J.F., Milodowski, A.E., Graham, C.C., Rushton, J.C., & Cuss, R.J. (2012). Evidence for gas-induced pathways in clay using a nanoparticle injection technique. Mineralogical Magazine, 76(8):3327–3336.

[3] Gonzalez-Blanco, L., Romero, E., Jommi, C., Li, X. & Sillen, X. (2016). Gas migration in a Cenozoic clay: Experimental results and numerical modelling. Geomechanics for Energy and the Environment, 6:81–100.

How to cite: Corman, G. and Collin, F.: A multi-scale model to study gas transport processes in clay materials, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16823, https://doi.org/10.5194/egusphere-egu24-16823, 2024.

17:50–18:00
|
EGU24-18283
|
On-site presentation
Identification of self-sealing processes in claystone-bentonite mixtures by digital image correlation
(withdrawn)
Nina Stoppe-Struck, Gerrit Langhorst, and Ben Laurich

Posters on site: Mon, 15 Apr, 10:45–12:30 | Hall X4

Display time: Mon, 15 Apr 08:30–Mon, 15 Apr 12:30
Chairpersons: Iman Rahimzadeh Kivi, Reza Taherdangkoo
X4.76
|
EGU24-5213
Reza Taherdangkoo, Thomas Nagel, and Christoph Butscher

The accurate determination of hydraulic conductivity in unsaturated bentonite is important for the modeling of subsurface thermo-hydro-mechanical and chemical processes. This study introduces a new hybrid approach, employing a constrained CatBoost algorithm coupled with a genetic algorithm for hyperparameter tuning. We benchmarked the effectiveness of the constrained CatBoost model against various data-driven regression models, including lasso, elastic net, polynomial regression, k-nearest neighbors, decision tree, bagging tree, random forest, and standard CatBoost. Our findings demonstrate that the constrained CatBoost model excels in providing accurate estimations of the hydraulic conductivity of compacted bentonite during the wetting phase. The model adequately captures the U-shaped correlation between hydraulic conductivity and suction and reflects the influence of temperature changes on hydraulic conductivity. Furthermore, the bootstrapping analysis, conducted across 800 iterations, confirms the stability and robustness of the constrained CatBoost model. This work provides a reliable tool for predicting hydraulic conductivity in diverse environmental and engineering contexts.

How to cite: Taherdangkoo, R., Nagel, T., and Butscher, C.: A  constraint-enhanced machine learning model for predicting hydraulic conductivity of unsaturated bentonite, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5213, https://doi.org/10.5194/egusphere-egu24-5213, 2024.

X4.77
|
EGU24-8326
Yeongkyoo Kim and Yeongjun Jang

Cobalt is often released into the natural environment through industrial waste from alloying industries and acid mine drainage. Additionally, it exists as a radionuclide (60Co) contributing to high-level radioactive waste. Smectite is a mineral that can be useful for adsorption and isolation of this element. In this investigation, Cheto-type montmorillonite (Cheto-MM), a source clay mineral of The Clay Mineral Society's (CMS) with well-established characteristics, was used as the primary material. The study aimed to assess how cobalt adsorption is affected by the adsorption site in the presence of interlayer water and after subsequent dehydration through heating. Adsorption kinetics and adsorption isotherm models were employed to explore the cobalt adsorption mechanism on Cheto-MM.

Results demonstrated notable variations in adsorption characteristics post-dehydration and subsequent shrinkage. Approximately 38% of cobalt was found to adsorb at the edge of Cheto-MM, while about 62% was adsorbed at the interlayer site, indicating the significant influence of the interlayer on cobalt adsorption in Cheto-MM. Adsorption kinetic models showed that the cobalt adsorption kinetics on Cheto-MM can be explained by a pseudo-second-order model. Moreover, isotherm experimental result was best represented by the Langmuir isotherm adsorption model. This study provides fundamental insights into cobalt adsorption characteristics on montmorillonite, emphasizing distinct adsorption sites. Such findings are instrumental in predicting smectite's adsorption behavior in high-level radioactive waste disposal sites in the future.

How to cite: Kim, Y. and Jang, Y.: Cobalt adsorption on montmorillonite: investigating the influence of dehydration on adsorption properties, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8326, https://doi.org/10.5194/egusphere-egu24-8326, 2024.

X4.78
|
EGU24-12464
|
ECS
Dmitry Bublik, Sebastian Mulder, and Johannes Miocic

Fluid extraction from geological formations for purposes of subsurface utilization leads to a pore pressure drop in reservoirs and potentially to compaction and seismicity. The geomechanical behaviour, and thus production-related compaction, of siliciclastic reservoirs is governed by the composition of reservoir sandstones, which includes porosity, grain size distribution, and detrital and authigenic mineralogy.

One siliciclastic reservoir which is undergoing compaction related to fluid extraction is the Rotliegend of the Groningen gas field. In this research, we investigate potential approaches to upscale the sandstone composition from the micro-scale (thin sections, plugs) to the well and reservoir scale eventually.

Previous research has highlighted that among the petrographic properties, the presence of authigenic clays tends to affect the geomechanical behaviour the most. Intragranular clays can affect the overall stiffness of the individual grains while the presence of the pore-filling clays  and clay rims may result in grain slip or pressure solution in loaded rock samples. In particular, we are defining the fractions of the clay-coating minerals (illite, kaolinite, etc.) and their effect on the inelastic deformation of reservoir sandstone as well as paying close attention to the presence of chlorite as its distribution corresponds to the areas associated with increased subsidence and seismicity.

We employed Short-Wave Infrared (SWIR) spectroscopy, a non-destructive and time-efficient technique, to obtain the mineralogical composition of core slabs from the Groningen Gas field. The SWIR data, based on a resolution of 200 µm pixels, allows for a detailed analysis of compositional variation within the Upper Rotliegend Group. Verification of the SWIR results against a comprehensive petrographic dataset, including X-ray diffraction, thin section descriptions, and modal point count analysis highlights that SWIR data primarily captures qualitative variations in mineralogy rather than providing precise numerical values. Combination of sedimentary facies, SWIR spectroscopy data and conventional petrographic studies allows to generate descriptive mineralogical trends for each of the studied wells.

Ultimately, the results of our research will serve as a foundation for selecting the samples, designing geomechanical experiments to test the proposed hypotheses, and as the means to select the upscaling and modelling approaches to develop a detailed model of the Dutch subsurface that matches well the existing heterogeneous structures. The derived 3D reservoir composition model of the Groningen gas field will be combined with the results of the deformation experiments to link reservoir composition to geomechanical behaviour. This will enable an updated, more realistic, 3D geomechanical model of the Groningen gas field that can be utilised by other researchers to better predict future compaction and subsidence.

How to cite: Bublik, D., Mulder, S., and Miocic, J.: From pore to well - identifying clay mineral distribution in sandstones using SWIR spectroscopy, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12464, https://doi.org/10.5194/egusphere-egu24-12464, 2024.

X4.79
|
EGU24-14166
Magdalena Krol, Farhana Chowdhury, Sifat Papry, Md Abdullah Asad, Pulin Mondal, Tarek Rashwan, and Ian Molnar

The use of bentonite clay in industrial applications is widespread: it is used as an engineered barrier for long-term management of radioactive wastes, CO2 storage, landfill liners, and contaminant containment. These applications have diverse environmental conditions ranging from various temperatures, pHs, saline contents, and ionic concentrations. Since bentonite is a low permeability clay, anion transport is diffusion dominated but geochemical reactions can also play a significant role and transport will be affected by environmental conditions. In this study, anion transport (bisulfide) under various conditions was examined using experimental and numerical techniques to understand the various geochemical and surface mediated reactions that are occurring in the bentonite. The case study presented is for the use of bentonite in long term storage of nuclear waste but can be extended to other applications.

First, diffusion experiments were performed to examine the transport and reactive nature of bisulfide (HS-) through bentonite compacted at dry density of 1090-1330 kg m-3. Experimental data of bisulfide transport were fitted using the inverse solution technique of Hydrus-1D model and different fitting parameters (e.g., diffusion, sorption, and reaction sink). Simulation results suggest that the HS- sorption/reaction affecting itsdiffusive transport through bentonite can be modeled using a simple nonlinear adsorption process.

Second, batch experiments were performed to understand the maximum allowable sorption that could take place under key geochemical conditions, including temperature, pH, and ionic strength. The results of batch sorption experiments performed suggest that HS- sorption increases with increasing temperature but decreases with increasing pH and ionic strength.

Lastly, since transport and reactive processes are interconnected, the results of these experiments were incorporated into a 1D transport COMSOL model to understand which geochemical process governs bisulfide transport through bentonite. Various processes were examined including linear and non-linear sorption, reactive transport, and anion exclusion. The model was validated using the experiments and showed that HS- was retained in the bentonite due to reactive processes and anion exclusion effects.

How to cite: Krol, M., Chowdhury, F., Papry, S., Asad, M. A., Mondal, P., Rashwan, T., and Molnar, I.: Anion Transport Through Bentonite Under Various Geochemical Conditions, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-14166, https://doi.org/10.5194/egusphere-egu24-14166, 2024.

X4.80
|
EGU24-4757
|
ECS
|
Erdi Avci and Atiye Tugrul

In the aggregate quarries in Çiftalan and Ağaçlı (Eyüpsultan, Turkey) regions, grovak-shale type rocks belonging to the Thracian Formation are used as concrete and asphalt aggregates. After crushing, these rocks are classified into 0-5 mm, 5-12 mm and 12-24 mm sizes and marketed as concrete and asphalt aggregates. Aggregates with sizes larger than 5 mm can be used directly in concrete and asphalt applications without any washing process, while aggregates with sizes smaller than 5 mm are called fine aggregates and are marketed after the washing process since the methylene blue value is too high for concrete and aggregate. After the washing process, grains smaller than approximately 100 microns in size (especially clay minerals) are stockpiled in waste sites. Within the scope of this study, it is aimed to produce sinter aggregates with strength that can be used in concrete from these waste materials. It was determined that the waste materials with a grain diameter of less than 100 micron consisted of 28.6% quartz, 21.9% albite, 17% muscovite, 24.9% chlorite and 3.6% calcite minerals. In addition to minerals, 3.4% organic matter was detected. In addition, the main oxide compositions of these wastes were determined, and it was found that 54.51% SiO2, 18.1% Al2O3, 19.66% sum of fluxes (CaO, Na2O, Fe2O3, K2O, MgO) and 5.8% loss on ignition. Within the scope of determining the thermal properties of the wastes, DTA analyzes were carried out and it was determined that two different endothermic reactions took place in the range of 500-600 °C and 700-800 °C, thus dehydroxylation reactions took place in these temperature ranges. In the range of 1000-1100 °C, it was determined that an exothermic reaction, that is, a new phase or several phases were realized. Considering these thermal properties of the wastes, 1100 °C was determined as the sintering temperature. For the sintering process, the samples were first dried at 105 °C and then milled and powdered again. The powdered samples were placed in a metal cell with a diameter of 5 cm and a height of 10 cm and then compressed with a load of 200 kg per cm2. The compressed samples were heated to 1100 °C in a high temperature furnace with a temperature increase of 10 °C per minute and kept at 1100 °C for 30 minutes and then allowed to cool in the furnace. Since the height of these sintered cylindrical specimens was less than 10 cm, the Brazilian test method was used to determine their strength. As a result of the experiment, it was determined that the Brazilian test strengths of 3 different specimens reached 8.8 MPa, 9.6 MPa and 15.7 MPa. Considering the strength values obtained, it was determined that ideal products for concrete and asphalt aggregate can be produced from the dust wastes released in these aggregate quarries. 

How to cite: Avci, E. and Tugrul, A.: Utilization of clay-containing aggregate sludge waste as structural concrete and asphalt aggregate, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4757, https://doi.org/10.5194/egusphere-egu24-4757, 2024.

X4.81
|
EGU24-16781
|
ECS
Hangbiao Song and Frédéric Collin

Deep underground geological disposal is widely accepted as one of the most appropriate ways for the long-term safety and management of radioactive waste. The heat generated by the nuclear waste decay brings elevated temperature increase, which may affect the thermo-hydro-mechanical (THM) behaviour of the host rock. A correct evaluation of the thermal impacts on the host rock behaviour is important for the design of the underground geological disposal.

With the temperature increases, thermal pressurization is observed both in the small (laboratory) and large (in-situ) scale tests. Physically, the overpressure induced by the discrepancy of the thermal expansion coefficient between the solid and fluid phases may potentially induce fracture re-opening and propagation. The host rock located in the middle of two adjacent cells may suffer shear or tensile failure, which is dependent on the intensity of the thermal power and the distance between the neighbouring cells. Some research work also shows that soil characteristics like cohesion, elastic modulus and water viscosity are influenced by the rise in temperature [1]. To investigate the thermally induced change on the mechanical property of host rock, triaxial compression tests were conducted at the University of Lorraine at different temperatures (20, 40, 60, 80, 100 and 150 °C), confining pressures (0, 4 and 12 MPa) and samples orientations (parallel and perpendicular to the bedding plane). The results showed the transitory overpressure induced by the thermal dilation during the initial heating, and the degradation of the mechanical strength of the host rock with the increase in temperature [2].

Based on the experimental observations, the triaxial compression tests are represented in a two-dimensional axisymmetric coupled THM model. The modelling is composed of the two steps: isotropic loading (increase confining stress and temperature), and shear process (increase axial loading). The numerical FEM code is LAGAMINE from the University of Liège. The Callovo-Oxfordian (COx) claystone, relying on its low permeability and good plasticity, has been selected as the host rock for the underground geological disposal in Meuse/Haute-Marne in France. The objective of this study is to introduce thermal-mechanical modelling involved with thermal plasticity. The cohesion of the host rock is defined as a function of the temperature to describe the thermally induced change of mechanical behaviour of the host rock. This model will then be validated against experimental observations in the laboratory and further applied to the large-scale heating test.

ACKNOWLEDGEMENT

This study was performed in the framework of the European Joint Programme on Radioactive Waste Management (EURAD). EURAD has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 847593.

REFERENCES

[1] Laloui, L. and Cekerevac, C., 2003. Thermo-plasticity of clays: an isotropic yield mechanism. Computers and Geotechnics, 30(8), pp.649-660. Doi : https://doi.org/10.1016/j.compgeo.2003.09.001.

[2] Gbewade, C.A.F., Grgic, D., Giraud, A. and Schoumacker, L., 2023. Experimental study of the effect of temperature on the mechanical properties of the Callovo-Oxfordian claystone. Rock Mechanics and Rock Engineering, pp.1-22. Doi : https://doi.org/10.1007/s00603-023-03630-7.

How to cite: Song, H. and Collin, F.: A thermal-mechanical constitutive modelling for Callovo-Oxfordian Claystone in the context of nuclear waste disposal, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16781, https://doi.org/10.5194/egusphere-egu24-16781, 2024.

X4.82
|
EGU24-8331
Mohamad Omidi and Thomas Baumann

Scaling refers to the accumulation of solid precipitates on the surfaces of pipes, heat exchangers, and other equipment across various industrial processes, notably within geothermal systems. This process can lead to decreased efficiency and increased maintenance needs, highlighting the importance of accurate prediction and control methods.

Hydrogeochemical models tend to overestimate the extent of scalings, especially if the scalings are caused by a disruption of the lime-carbonic acid equilibrium due to degassing in the geothermal fluid. This is because the models are not capable to describe diffusion-limited crystal growth, partial volume effects, and local saturation states adequately. In this study we introduce a novel approach by coupling a multiphase CFD-model (OpenFOAM) for the description of the gas phases in the produced geothermal fluids with PhreeqC to simulate the hydrochemical effects of the stripping of CO2 by the gas phase. This results in a model with high spatial and temporal resolution which allows to quantitatively differentiate between processes in the fluid and the processes taking place at the solid interfaces (pipe walls or matrix).

The code is validated with published experiments in bubble columns. The coupling results in a highly flexible model which can account for different hydrochemical conditions, different matrix, and varying gas composition using a well established thermodynamic database. Transfer to other hydrochemical conditions is therefore facilitated.

How to cite: Omidi, M. and Baumann, T.: Coupling OpenFOAM and PhreeqC to quantify local disruptions of hydrochemical equilibria due to bubble formation and stripping of CO2, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8331, https://doi.org/10.5194/egusphere-egu24-8331, 2024.

X4.83
|
EGU24-20645
|
ECS
Lucy Tweed, Jerome Neufeld, and Mike Bickle

Depleted oil and gas reservoirs are attractive sites for CO2 sequestration. However, the injection of CO2 into depleted reservoirs carries the potential for significant Joule-Thomson cooling, when dense, supercritical COis injected into a low-pressure reservoir. The resulting low temperatures around the well-bore risk causing thermal fracturing and/or freezing of pore waters or precipitation of gas hydrates which would reduce injectivity and jeopardise near-well stability. Injection into reservoirs at subcritical pressure also leads to a phase transition from liquid to vapour CO2. This is accompanied by cooling, due to the latent heat of vaporisation, and dramatic changes in fluid properties including density, compressibility and viscosity.

We present models of non-isothermal flow of CO2 in the near well-bore region, which demonstrate the controls on cooling and constrain the different pressure-temperature regimes that can emerge. We show that during radial injection, with fixed injection rate, transient Joule-Thomson cooling can be described by similarity solutions at early times. The positions of the CO2 and thermal fronts are described by self-similar scaling relations. The scaling analysis here identifies the parametric dependence of Joule-Thomson cooling. We present a sensitivity analysis which demonstrates that the primary controls on the degree of cooling are the injection rate and Joule-Thomson coefficient. The analysis presented provides a computationally efficient approach to assessing the degree of Joule-Thomson cooling expected during injection start-up, providing a complement to full numerical simulations.

How to cite: Tweed, L., Neufeld, J., and Bickle, M.: Joule-Thomson cooling and phase transitions during CO2 injection in depleted reservoirs, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-20645, https://doi.org/10.5194/egusphere-egu24-20645, 2024.

X4.84
|
EGU24-17117
|
ECS
|
Craig Allsop, Matteo Pedrotti, and Alessandro Tarantino

The successful deployment of carbon dioxide (CO2) geological sequestration in porous media is reliant on the sealing efficiency of the overlying, clay-rich caprock to act as a physical barrier. Clay-rich caprock formations are considered favorable materials to act as a seal due to them characteristically consisting of small pores providing high capillary entry pressures, hence preventing the intrusion of a non-wetting fluid.

The juxtaposition and availability of deep seated (buried) caprock-reservoir systems to carbon capture and storage clusters may not be available. Therefore, the assessment of shallow seated, weakly consolidated caprock-reservoir systems (e.g., Sleipner) will be required. Our experimental campaign tests analogous caprock geomaterials which have relatively high compressibility, representative of shallow seated (buried or less indurated) clay-rich caprocks.

Past experimental campaigns demonstrate that CO2 breakthrough is dominated by the creation of very localized channels across the sealing barrier which occur at pressures far lower than the one predicted by the Laplace’s equation [1]. However, limited data characterizing these pathways exists. Furthermore, the physical indicators of susceptibility which underly the micro-mechanisms of failure (e.g., fracturing), are still only postulated for clay-rich geomaterials.

where Pc* is the capillary breakthrough pressure [kPa], ψ, reflects pore shape [-], Ts, interfacial tension between water and gas (e.g., CO2), and θ, represents wettability [°].

An innovative experimental set-up which allowed for the onset of surface crack formation to be captured during gas injection (representing the non-wetting fluid in CO2 geological sequestration) into intact clay-rich geomaterials is presented. This allowed for the investigation of physical indicators of susceptibility to gas breakthrough via localized pathways.

Results on different fracture patterns when non-wetting gas (i.e., air) is injected into consolidated clay show the formation of large fractures that nucleate from within the sample. Upon air pressurization, before fracture formation, the sample undergoes volumetric deformation (i.e., consolidation), as the resulting action of the vertical stress applied at the air-water interface (menisci). Once a fracture forms deformation stops and breakthrough occurred at lower pressures than traditionally recorded. The mechanisms of air intrusion are expected to be of a similar nature as CO2 intrusion. Post-mortem assessment of the internal nature of these localized pathways was then visualized using xCT imaging.

As a continuum mechanics framework will not predict fracture formation under our test conditions, it appears that the experimental evidence support the underlying hypothesis that disjoining pressure governs the mechanisms that ultimately control fracture formation and thus, eventually CO2 breakthrough. The disjoining pressure is governed by the electrostatic double-layer interactions, van der Waal’s dispersion forces, structural forces, and solvation forces.

If the pore size distribution is such that high gas pressures are required to overcome capillarity, the gas pressure will force single clay particles apart, displacing water form adjacent interparticle spaces. This represents a localized failure mechanism at the clay interface, resulting in fracture nucleation. It is expected that clay displaying large swelling pressures will subsequently display high gas entry pressure, termed “pathway dilation”. Therefore, pressurized gas will enter a dilated pathway at lower pressures than anticipated.

How to cite: Allsop, C., Pedrotti, M., and Tarantino, A.: Insight into the micromechanisms of gas breakthrough in water-saturated clay-rich geomaterials – Implications for CO2 sequestration, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17117, https://doi.org/10.5194/egusphere-egu24-17117, 2024.

X4.85
|
EGU24-18636
|
ECS
Kalliopi Tzoufka, Guido Blöcher, Mauro Cacace, Daniela Pfrang, and Kai Zosseder

Concepts of High-Temperature Aquifer Thermal Energy Storage (HT-ATES) (> 50 °C) are investigated for system application in the Upper Jurassic reservoir (Malm aquifer) of the German Molasse Basin (North Alpine Foreland Basin). The karstified and fractured carbonate rocks exhibit favorable conditions for conventional geothermal exploitation of the hydrothermal resource. Here, to further assess the sustainability of HT-ATES development in the Upper Jurassic reservoir, a physics-based numerical analysis is performed. With an estimated heating capacity of ca. 21 MW over half a year, our approach aims at determining numerically the efficiency of heat storage under the in-situ Upper Jurassic reservoir conditions and locally feasible operation parameters.

The numerical models build upon datasets from three operating geothermal sites at depths of ca. 2000-3000 m TVD located in a subset of the reservoir which is dominated by karst-controlled fluid fluxes. Commonly considered as a single homogeneous unit, the 500 m thick reservoir is subdivided into three discrete layers based on field tests and borehole logs from the three considered sites. This introduced vertical heterogeneity with associated layer-specific enhanced permeabilities allows to examine potentially arising favorable heat transfer, and in combination with the facilitated high operation flow rates to evaluate thermal recoveries in the multilayered reservoir.

Computation results reveal that the reservoir layering induces preferential fluid and heat migration primarily into the high-permeability zone, while thermal front propagation into the lower permeable rock matrix is restricted. The simulations further display the gradual temperature increase in the warm wellbore and its surrounding host rock, and the consequent progressive improvement in the heat recovery efficiencies. Despite the elevated permeability that may trigger advective heat losses, heat recovery factor values range from ca. 0.7 over the first year of operation to over 0.85 after 10 years of operation. An additional scenario is examined with fluid injection solely in the high permeable zone, in order to quantify potential improvement in the recovery efficiency by omitting fluid injection in the lower-permeability layers where heat propagation is diminished. This is due to the geometrical shape of the thermally perturbed rock volume as heat losses occur at the interface between thermal front and adjacent reservoir rock. Consequently, conclusions on the performance of the two different system designs under this layered reservoir setting are derived. All simulations account for density and viscosity variation through the IAPWS (International Association for the Properties of Water and Steam) thermodynamic property formulations. Results show that density-induced buoyant fluxes which would considerably decrease thermal efficiencies are inhibited, and thus the prevailing heat transport mechanism is forced convection.

How to cite: Tzoufka, K., Blöcher, G., Cacace, M., Pfrang, D., and Zosseder, K.: Physics-based numerical evaluation of High-Temperature Aquifer Thermal Energy Storage (HT-ATES) in the Upper Jurassic reservoir of the German Molasse Basin, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18636, https://doi.org/10.5194/egusphere-egu24-18636, 2024.

X4.86
|
EGU24-21553
Guido Blöcher, Liang Pei, Stefan Kranz, Christian Cunow, Lioba Virchow, and Ali Saadat

Surplus heat as stored in an ATES (Aquifer Thermal Energy Storage) system in summer could partly meet the increasing demand of energy in winter. Better assessment on the wellbore integrity permits sustainable operation of such systems. Therefore, an artesian flow test was conducted in a research well located in Berlin, Germany. In this test, artesian flow of 16.8°C from Jurassic sand at depths from 220 m to 230 m was produced at 14°C and at a flow rate of 3.6m3/h from the annular space between the production casing and the anchor casing. The depth-resolved temperature at the production casing as monitored using the distributed temperature sensing (DTS) technique manifested the depths of the artesian aquifer. A hydro-thermal coupled numerical model for the artesian flow was calibrated by matching the simulated flow rate to the wellhead-measured values. The simulated and the DTS-monitored temperatures suggested that the heating-up in the near-wellbore materials by the artesian flow was hindered by the deployment-related inclusion of water behind the anchor casing, and the cooling in these materials in the shut-in test stage was enhanced by such inclusion. 

How to cite: Blöcher, G., Pei, L., Kranz, S., Cunow, C., Virchow, L., and Saadat, A.: Analysis of Wellbore Integrity using DTS Monitoring and Numerical Modelling in the Practice of ATES, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-21553, https://doi.org/10.5194/egusphere-egu24-21553, 2024.

X4.87
|
EGU24-19122
|
ECS
Aurelio Andy, Felix Schölderle, Daniela Pfrang, and Kai Zosseder

In order to gain insights into the hydraulic and thermal long-term behavior of deep geothermal boreholes during operation, a fiber-optic monitoring (FOM) system was implemented at the Stadtwerke Munich's (SWM) geothermal plant 'Schäftlarnstrasse' in the national project Geothermal-Alliance Bavaria (GAB). The permanently installed fiber optic cables enable continuous measuring of the temperature in boreholes with high spatial and temporal resolution via Distributed Temperature Sensing (DTS) and of acoustics/vibration via Distributed Acoustic Sensing (DAS). In 2019, one production and one injection well were equipped with FOM in different setups and have been collecting data since then. In the injector, the cable was cemented behind the anchor pipe in the upper section of the borehole, whereas the producer was equipped with a fiber-optic cable until the total depth of the reservoir. Additionally, a fiber-optic gauge provides pressure and temperature at the top of the reservoir in this well.

Among other things, precise inflow profiling in the reservoir section of the producer could be carried out using the temperature data, which helps deepen the understanding of the hydraulic behavior of the reservoir. Changes over time in the temperature of the produced thermal water can be traced back to changes in the hydraulically active zones or inflow temperatures in the reservoir. Other factors that influence the wellhead temperature, such as borehole heat losses to the surrounding rock, can be quantified using the DTS data.

At the beginning of 2024, a third fiber-optic cable will be installed in an injection well at the Schäftlarnstrasse geothermal site. A similar setup as in the production well will allow for continuous measuring in the entire borehole until total depth (TD) and pressure and temperature from p/T gauges at the top of the reservoir and at TD. Therefore, both production and injection conditions in the reservoir and wellbore will be continuously monitored at all depths using DTS and DAS for the first time.

We present the first interpretations of the data collected from the newly installed FOM in the injection well as well as insights into the long-term monitoring of the hydraulic and thermal dynamics in the production well to underline the importance of permanent downhole monitoring to ensure efficient and sustainable use of the geothermal reservoir.

How to cite: Andy, A., Schölderle, F., Pfrang, D., and Zosseder, K.: Operational Monitoring of Thermal Dynamics in Deep Geothermal Production and Injection Wells with Fiber Optics from the Surface to the Reservoir, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-19122, https://doi.org/10.5194/egusphere-egu24-19122, 2024.

X4.88
|
EGU24-6340
Sandro Andrés, David Santillán, Miguel Ángel Santoyo, and Luis Cueto-Felgueroso

Induced seismicity has attracted increasing research interest in recent times. This phenomenon is generally associated with fluid injection or extraction wells, in energy industry activities such as hydrocarbon extraction/injection, CO2 sequestration, geothermal energy or underground storage of green hydrogen. However, there are other human activities that can induce or trigger seismic events, such as oscillations in the surface water level due to the construction of hydraulic infrastructures.

Surface water level oscillations, whether natural or anthropogenic, can alter the pore pressure regime in the subsurface. Although these alterations usually have a moderate magnitude, they can destabilize faults that were already close to overcoming their slip resistance and eventually trigger an earthquake. On the other hand, the fault slip itself alters the pressure regime due to the undrained response of the porous medium, caused by the sudden deformation of the fault surrounding area. This undrained pressure, which can take up to weeks to dissipate, can alter the stress state of other nearby fractures and trigger new earthquakes or aftershocks.

In this work we present numerical simulations of the underground area around the Itoiz dam (Spain). We simulate the subsurface as a saturated poroelastic medium, with a fully coupled scheme between fluid flow and solid deformation in the porous medium. Through a seismological analysis we start from series of recent earthquakes to locate geological faults. With our numerical model we study how the undrained effect produced by this series of events can destabilize other nearby faults. We also add the effects of the oscillations in the reservoir level following the historical series of the last years.

Our numerical simulations indicate that in the case of the Itoiz dam, the most recent seismic swarm could have been triggered by the stress transfer from the previous events, which, together with the filling of the reservoir, may have destabilized faults that were critically stressed for failure. Our simulations can contribute to explore how the combined effect of the undrained pressure by the fault slip and the oscillations of the reservoir can trigger faults that, due to the natural state of stress, are already close to slip. Also to clarify if the most relevant mechanism is the oscillation of the reservoir level or the undrained pressure trigger, which according to bibliographic analysis can be of the same magnitude. This could apply to other cases of seismicity induced by hydraulic infrastructures or in general by oscillations in the surface water level.

Acknowledgements

This research Project has been funded by the Comunidad de Madrid through the call Research Grants for Young Investigators from Universidad Politécnica de Madrid under grant APOYO-JOVENES-21-6YB2DD-127-N6ZTY3, RSIEIH project, research program V PRICIT.

How to cite: Andrés, S., Santillán, D., Santoyo, M. Á., and Cueto-Felgueroso, L.: Stress transfer and poroelastic mechanisms to elucidate seismicity triggered by reservoirs, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6340, https://doi.org/10.5194/egusphere-egu24-6340, 2024.

X4.89
|
EGU24-10763
|
ECS
Tatia Sharia, Andreas Rietbrock, Birgit Müller, and Thomas Niederhuber

Understanding perturbations caused by underground fluid injection and extraction is essential for steady long-lasting operations of subsurface energy systems (e.g. geothermal systems). The systems are often too complex to obtain precise analytical solutions and computationally expensive to have fine numerical results. Simplified settings give the benefit of understanding the role of driving geological parameters as well as examining the limits of each approach. In this study, we focus on the influence of permeability on pore pressure and present analytical and numerical solutions of spatial and temporal evolutions of pore pressure in an elastic homogeneous porous media with isotropic and anisotropic permeabilities. We use COMSOL Multiphysics to build a 3D finite element model with injection/production wells and investigate where the numerical solutions of the spatially limited volume coincide and diverge from corresponding analytical solutions of pore pressure in infinite media.

How to cite: Sharia, T., Rietbrock, A., Müller, B., and Niederhuber, T.: Pore pressure evolution in media with isotropic and anisotropic permeability - analytical and numerical solutions, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10763, https://doi.org/10.5194/egusphere-egu24-10763, 2024.

X4.90
|
EGU24-2963
|
ECS
Baibing Yang, Qingfeng Meng, and Fang Hao

The gas-bearing Yinggehai Basin in the northern South China Sea is characterized by high pressure and temperature conditions. The Miocene strata in the eastern slope exhibit intense overpressure with a pressure coefficient exceeding 2.2. Understanding the characteristics of overpressure and clarifying its causes are very important for natural gas exploration and development in reservoirs with intense overpressure. This study identifies and characterizes the pressure structure in the eastern slope of Yinggehai Basin by utilizing formation pressure data from formation tests (MDT) and drill pipe tests (DST) in combination with conventional logging data (acoustic wave, density, and resistivity). The pressure distribution exhibits a distinct three-stage ladder structure. Overpressure initiates in the Lower Pliocene Yinggehai Formation, with a pressure coefficient of 1.2. The formation pressure increases sharply in the underlying Upper Miocene strata (1st member of the Huangliu Formation), with a pressure coefficient ranging from 1.4 to 1.8. The pressure coefficient increases to between 2.1 and 2.3 in the 2nd member of the Huangliu Formation.

We utilized the logging curve combination method, Bowers effective stress method, and Bowers acoustic-density crossplot method to differentiate between two types of overpressure origins: loading and unloading type, and further evaluated their contribution rate to overpressure. Loading overpressure primarily accounts for the overpressure observed in the Yinggehai Formation, contributing to pore pressure with a range from 51% to 93%. In contrast, the overpressure in the Huangliu Formation is predominantly of the unloading type, contributing to pore pressure from 35% to 46%. Additionally, we analyzed the genetic mechanism of overpressure, utilizing lithology, subsidence rate, organic matter maturity, and seismic attribute data. We find that the loading overpressure in the Yinggehai Formation is attributed to unbalanced compaction of mudstone due to rapid sedimentation (sedimentation rate > 500m/ Ma). The unloading intense overpressure in the Miocene strata is most likely caused by the vertical transmission of overpressured fluid along faults and fractures. Our findings provide implications for complex pressure structure, overpressure evaluation, and genetic mechanisms in sedimentary basins.

How to cite: Yang, B., Meng, Q., and Hao, F.: Intense fluid overpressure in the eastern slope of the Yinggehai Basin, South China Sea, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2963, https://doi.org/10.5194/egusphere-egu24-2963, 2024.