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.

Cobalt is often released into the natural environment through industrial waste from alloying industries and acid mine drainage. Additionally, it exists as a radionuclide (⁶⁰Co) 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.

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.

The use of bentonite clay in industrial applications is widespread: it is used as an engineered barrier for long-term management of radioactive wastes, CO₂ 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.

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.

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.

Depleted oil and gas reservoirs are attractive sites for CO₂ sequestration. However, the injection of CO₂ into depleted reservoirs carries the potential for significant Joule-Thomson cooling, when dense, supercritical CO₂ is 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 CO₂. 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 CO₂ 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 CO₂ 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.

The successful deployment of carbon dioxide (CO₂) 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 CO₂ 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 [-], T_s, interfacial tension between water and gas (e.g., CO₂), 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 CO₂ 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 CO₂ 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 CO₂ 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.

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.

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.6m³/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.

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.

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.

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.

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 (1^st 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 2^nd 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.