Pore Pressure Changes Research Articles

Coupled hydro-mechanical model for underground hydrogen storage undergoing cyclic injection and production

ABSTRACT Excess renewable energy can be used to generate and store hydrogen underground, offering a secure and econimical solution to the intermittency challenges of renewable energy sources. The success of underground storage relies on understanding and controlling the geo-mechanical integrity of storage sites during cyclic loadings. While geological deformation during natural gas and carbon dioxide storage is well studied, few have investigated pore pressure and in-situ stress variations during cyclic loading, which is unique to underground hydrogen storage sites. This paper presents a predictive analytical model for the pore pressure and stress changes during cyclic loading. The model was derived based on continuous point source injection, a basic hydro-mechanical model while applying the superposition to consider the cyclic loading. The developed model was evaluated for hydrogen and natural gas subjected to cyclic injection and production. Comparing hydrogen and natural gas storage cases can offer valuable insights into the behavior of hydrogen since the deformation of natural gas storage is more studied. Finally, analytical solutions for both storage cases were validated by numerical models with discrepancies between the models amounting to less than 1%. Comparative analysis of the results from 3 to 99 cycles of injection and production for both hydrogen and methane demonstrates that rise in cycle numbers leads to a corresponding rise in stress accumulation in the system. Neglecting the stress accumulation during these cycles can lead to significant geomechanical issues, potentially compromising the storage system’s safety and integrity.

Dam impoundment near active faults in areas with high seismic potential: Case studies from Bisri and Mseilha dams, Lebanon

Reservoir induced seismicity is caused by stress changes due to the impoundment of water behind dams. In seismically active areas, the presence of critically located active faults makes the impoundment of water behind dams a seismic safety risk. Dam projects in Lebanon have become a soaring example of complacency and negligence that has overlooked the concerns for seismic safety raised over the projects and their high potential of inducing seismicity. In this paper, we use 2D and 3D fully coupled poroelastic modeling to assess the risk of dam impoundment on seismogenic faults located near dam sites in Lebanon. The coulomb failure stresses are calculated along the faults, and their variations are observed in relation to changes in pore pressures and normal stresses. In addition, the expected maximum earthquake magnitudes are computed along those faults. Our results show a high risk for reservoir induced seismicity on faults that are either underneath the reservoir or hydraulically connected to a fault beneath the reservoir. Consequently, the studied dams would present a serious hazard of induced seismicity in time where the region is already at high risk of destructive earthquakes after the catastrophic seismic events that struck Turkey and Syria on 6 February 2023 on the Eastern Anatolian Fault, which is connected to the Dead Sea Transform Fault that passes through Lebanon.

GEOMECHANICAL MODELING OF RESERVOIR DYNAMICS ASSOCIATED WITH SHALLOW INJECTION AND PRODUCTION IN THE DELAWARE BASIN

The Permian Basin is an area of active hydrocarbon production and saltwater disposal, as well as associated induced seismicity and other geomechanical responses that threaten the surface environment. Among the important but not yet fully answered questions are: (i) what are the temporal and spatial changes in stress and pore pressure in response to fluid injection and production, and (ii) how do these changes relate to slip on pre-existing faults and surface deformation? We simulate stress and pore pressure responses to saltwater injection in the Delaware Mountain Group and production from the Wolfcamp and Bone Spring Formations with the aid of 2D geomechanical simulations that capture key mechanical stratigraphic units and representative faults. Primary loading consists of localized pore pressure changes that represent fluid injection into the Delaware Mountain Group and production from the lowermost portion of the Bone Spring Formation, the entire Wolfcamp A Formation, and the uppermost portion of the Wolfcamp B Formation. Injection leads to pore pressure increase, vertical extension, reduction in mean stress, and increase in differential stress in the Delaware Mountain Group. Production leads to decrease in pore pressure, vertical contraction, increase in mean stress, and increase in differential stress in the Bone Spring and Wolfcamp A and B Formations. The net effect of injection and production in our generalized simulations is normal faulting slip on faults that are relatively shallow in the subsurface, similar to faults that are known to have produced seismogenic rupture. The combination of injection and production reproduces a spatially-variant trend in uplift and subsidence, consistent with regional patterns measured in the study area via InSAR analysis. The modeled scenarios with shallow injection and/or production caused only small (lt;0.2 MPa) stress perturbations to propagate downward to basement, which would be unlikely to cause instability of deep-seated seismogenic faults.

Variations in Temperature and Pressure in the “Reservoir–Well” System Triggered by Blasting Recovery of Iron Ore at the Kursk Magnetic Anomaly

This paper presents the results of precise measurements of temperature and pore pressure in the “reservoir–well” system during the development of iron ore deposits of the Kursk Magnetic Anomaly (KMA) via blasting. For the observation period from October 2021 to June 2024, variations in compressibility, permeability and temperature in the upper Albian-Cenomanian confined aquifer, which is used for district water supply, were determined. The general trend in a decrease in water temperature was traced (from 12 °C to 11.4 °C). It was accompanied by an increase in the hydrostatic head (from 3.7 m to 7.4 m). Water temperature in the upper aquifer was measured for 9 industrial explosions in the mine and for 30 explosions in the quarry. For one explosion in the mine and five explosions in the quarry the coseismic changes in water temperature with amplitudes of 0.06–0.09 °C were established, while changes in pore pressure in the “reservoir–well” system were 0.4–2.2 kPa. Local changes in the permeability of the reservoir in the vicinity of the well (the skin effect) are considered to be the main factor that controls the coseismic response of temperature during industrial explosions. As the reservoir permeability increases, the water temperature in the “reservoir–well” system can decrease and vice versa. The same pattern was observed according to regime measurements performed in 2022–2023. The recorded coseismic responses of water temperature in the upper aquifer in the high-frequency range are similar to the effects observed during propagation of seismic waves originated from earthquakes in the low-frequency range at different sites all over the world for the seismic energy density of 0.05–0.45 J/m3. The observed variations in aquifer temperature in the “reservoir–well” system under episodic dynamic impacts are of particular interest from the point of view of activating hydrogeochemical processes that accompany the development of iron ore deposits.

A numerical analysis of Thermo–Hydro–Mechanical behavior in the FE experiment at Mont Terri URL: Investigating capillary effects in bentonite on the disposal system

We investigated thermo–hydro–mechanical (T-H–M) coupled behavior observed during the full-scale heater emplacement experiment at the Mont-Terri underground research laboratory conducted in the Opalinus clay as part of the DECOVALEX-2023 Task C project. Utilizing the OGS-FLAC simulator, we created a three-dimensional model to simulate multiphase flow in the experiment, applying extended Philip and de Vries’ model and incorporating the anisotropic T–H–M properties of the Opalinus clay. The simulation, which included a ventilation process, spanned five years of heating experiments and successfully replicated the measured temperature, pore pressure, displacement, and relative humidity results in bentonite and host rock during the experiment. The analysis revealed that capillary pressure significantly influenced the pore pressure change in the host rock near the tunnel, while thermal pressurization became dominant with increasing distance. Consequently, we conducted a sensitivity analysis on a simplified model to evaluate the effect of capillary pressure on the disposal system. Capillarity is a dominant factor for the multiphase flow depending on the distance from the heat. Variations in capillary pressure were observed depending on the gas entry pressure and water retention model, indicating that the capillarity of unsaturated bentonite could inherently affect the T–H–M results within the disposal system.

Study on the Effect of Liquefiable Overburden Foundations of Rockfill Dams Based on a Pore Pressure Model

China’s southwestern region boasts abundant hydropower resources. However, the area is prone to frequent strong earthquakes. The areas surrounding dam sites typically have deep overburden, and the liquefaction of saturated sand foundations by earthquakes poses significant safety risks to the construction of high dams in the southwest. The effects of liquefaction and reinforcing measures on the foundations of rockfill dams on liquefiable overburden under seismic action are currently the subject of somewhat unsystematic investigations. The paper utilizes the total stress and effective stress methods, based on the equivalent linear model, to perform numerical simulations on the overburden foundations of rockfill dams. The study explores how factors such as dam height, overburden thickness, liquefiable layer depth, liquefiable layer thickness, ground motion intensity, and seismic wave characteristics affect the liquefaction of the overburden foundations. Additionally, it examines how rockfill dams impact the dynamic response, considering the liquefaction effects in the overburden. The results show that although the total stress method, which ignores the cumulative evolution of pore pressure during liquefaction, can reveal the basic response trend of the dam, its results in predicting the acceleration response are significantly biased compared to those of the effective stress method, which comprehensively considers the cumulative changes in liquefaction pore pressure. Specifically, when the effect of soil liquefaction is considered, the predicted acceleration response is reduced compared to that when liquefaction is not considered, with the reduction ranging from 4% to 30%; with increases in the thickness and burial depth of the liquefiable layer, the effective stress method considering liquefaction significantly reduces the predicted peak acceleration; the effect of liquefiable soil on the attenuation of the speed response is more sensitive to the low-frequency portion of the seismic wave. The study’s findings are a significant source of reference for the planning and building of rockfill dams on liquefiable overburden.

Mechanical integrity analysis of caprock during the CO2 injection phase

Abstract CO2 geological storage is one of the important means to mitigate the greenhouse effect. The safe underground storage of CO2 largely depends on the mechanical integrity of the caprock. This paper establishes a fluid-solid coupling model for CO2 geological storage to study the changes in pore pressure, vertical displacement, and effective stress in the caprock during the CO2 injection process. Combined with the Mohr-Coulomb criterion, the study determines whether mechanical failure occurs in the caprock. The results indicate that, at the beginning of CO2 injection, significant changes occur in the pore pressure, vertical displacement, and effective stress at the bottom of the caprock near the injection well, which then tend to stabilize; the maximum pore pressure at the bottom of the caprock reaches 36.08 MPa; the caprock near the injection well is considered the most critical area, where the risk of mechanical failure is highest; at the end of CO2 injection, the stress state does not reach the limit, and the caprock remains stable.

Analysis of stress field in the head area of the Three Gorges Reservoir based on coupled fluid-solid theory

The Three Gorges Reservoir, one of the largest water conservation system in the world, has been of keen interest to scientists globally since its impoundment. After construction of the dam, there has been a significant increase in seismic activity in the head area of the reservoir. It is generally accepted that earthquakes in this region are predominantly caused by the Jiuwanxi and Xiannvshan faults. This study focused on the stress changes occurring in the research area. A three-dimensional finite element model of the reservoir area was constructed using the geological structure and digital ground elevation data of the reservoir area. The fluid-solid coupling theory was applied to calculate the dynamic spatial changes in pore pressure and Coulomb stress in the faults and surrounding rocks during reservoir impoundment. The findings indicated that the added head pore water pressure at the bottom of the reservoir had a maximum impact range of approximately −2800 m on the surrounding rock, whereas the Xiannvshan and Jiuwanxi faults had a maximum diffusion range of approximately −4300 m. Rock permeability also played a significant role in the water storage process. During the 1 56 m water impoundment stage, owing to rapid water storage activity, stress could not be transmitted to both sides in a timely manner, resulting in the formation of an extreme stress change zone at −4000 m inside the fault. This may have been the reason for the frequent earthquakes during this stage. The 17 5 m cycle water storage stage also exhibited a significant degree of seismicity, potentially attributable to the long-term infiltration of reservoir water and accumulation of stress in the previous stage. The stress in the study area at the four stages are in a process of accumulation-release-accumulation-release.

Semi-analytical solutions for land subsidence due to groundwater withdrawal

This paper introduces novel semi-analytical models tailored for estimating land subsidence resulting from groundwater extraction in confined aquifers. These models offer high scalability, allowing them to be applied to various well configurations and pumping schedules. Their development involves the numerical integration of two key analytical solutions: the “nucleus of strain” (NoS) (Mindlin and Chen, 1950), which represents a localised zone within the aquifer where a unit change in pore pressure leads to deformation and subsequent surface displacement, and the classic Theis equation (Theis, 1935) for the pore pressure changes induced by a constant-rate well pumping from a laterally unbounded aquifer. These integrations yield surface displacement components, both horizontal and vertical, expressed as functions of two dimensionless spatial–temporal variables, which encompass aquifer depth, thickness, well placement, pumping schedules, and critical hydro-geomechanical parameters like hydraulic conductivity, porosity, vertical compressibility, and water compressibility. Proposed are two distinct modelling approaches: one employing a lookup table (LT) derived from numerical integration results, and the other providing direct closed-form surface displacement solutions by fitting LT data with “hinge models”, which use piecewise-linear functions linked by sigmoidal curves for computational efficiency. In both cases, surface displacement components are estimated by plugging in the dimensionless variables. Conditions of variable pumping from multiple wells can be addressed by applying superposition of solutions. In essence, these semi-analytical models offer swift computational capabilities for understanding and forecasting land subsidence dynamics. Their scalability makes them adaptable to a wide array of well configurations and scheduling scenarios, rendering them valuable for numerous applications. They are particularly significant for providing preliminary estimates of the impacts of groundwater development, conducting “what-if” tests, and performing sensitivity analyses to identify key factors affecting land subsidence risk. This underscores the importance of these models in sustainable groundwater resource management and in mitigating land subsidence and its associated consequences.

Pore pressure-minimum horizontal stress coupling estimation in Tunu Field, Lower Kutai Basin, Indonesia

Quantifying the change of minimum horizontal stress (Shmin) in response to pore pressure is vital, particularly at the development stage of a hydrocarbon field. However, estimates of the coupling ratio between these parameters, either at basin or reservoir scales, are not well defined in Indonesia’s sedimentary basin, such as Lower Kutai Basin. This study investigates the coupling ratio between Shmin and pore pressure in one of the gas fields within the basin, Tunu Field, using available evidence from an extensive data set of direct pressure measurement, leak-off test, extended LOT, Mini-Frac, and lost circulation. Previous studies have reported overpressured and depleted reservoirs. Nevertheless, the evolution of Shmin due to the change in pore pressure has yet to be well established. This study reveals a Shmin and vertical stress ratio of 0.87 which may relate to the stiffening of the shales via clay diagenesis/cementation. Coupling ratios between Shmin and pore pressure are 0.3 for the overpressured zone at the basin scale and 0.42 for the depletion at the reservoir scale. The approach used has been shown to produce similar coupling ratio (or stress path) results to other basins globally, thus providing a useful estimate of the Shmin and pore pressure relationship for the Tunu Field, and potentially other fields in the Lower Kutai Basin. For the application to infill drillings within the field, this study indicates Shmin decreases in response to pore pressure depletion due to production, ranging from 1.7 to 2.3 ppg within the hydrostatic zone and from 2.5 to 3.7 within the overpressured zone, equivalent to a Pore Pressure (PP)-Shmin window decrease of approximately 28% to 37% and approximately 40% to 88%, respectively. In addition, strike-slip faulting is suggested as the likely in situ stress regime in the Tunu Field.

Dynamic Evolution Law of Production Stress Field in Fractured Tight Sandstone Horizontal Wells Considering Stress Sensitivity of Multiple Media

Inter-well frac-hit has become an important challenge in the development of unconventional oil and gas resources such as fractured tight sandstone. Due to the presence of hydraulic fracturing fractures, secondary induced fractures, natural fractures, and other seepage media in real formations, the acquisition of stress fields requires the coupling effect of seepage and stress. In this process, there is also stress sensitivity, which leads to unclear dynamic evolution laws of stress fields and increases the difficulty of the staged multi-cluster fracturing of horizontal wells. The use of a multi-stage stress-sensitive horizontal well production stress field prediction model is an effective means of analyzing the influence of natural fracture parameters, main fracture parameters, and multi-stage stress sensitivity coefficients on the stress field. This article considers multi-stage stress sensitivity and, based on fractured sandstone reservoir parameters, establishes a numerical model for the dynamic evolution of the production stress field in horizontal wells with matrix self-supporting fracture-supported fracture–seepage–stress coupling. The influence of various factors on the production stress field is analyzed. The results show that under constant pressure production, for low-permeability reservoirs, multi-stage stress sensitivity has a relatively low impact on reservoir stress, and the amplitude of principal stress change in the entire fracture length direction is only within the range of 0.27%, with no significant change in stress distribution; The parameters of the main fracture have a significant impact on the stress field, with a variation amplitude of within 2.85%. The ability of stress to diffuse from the fracture tip to the surrounding areas is stronger, and the stress concentration area spreads from an elliptical distribution to a semi-circular distribution. The random natural fracture parameters have a significant impact on pore pressure. As the density and angle of the fractures increase, the pore pressure changes within the range of 3.32%, and the diffusion area of pore pressure significantly increases, making it easy to communicate with the reservoir on both sides of the fractures.

The geomechanical response of the Gulf of Mexico Green Canyon 955 reservoir to gas hydrate dissociation: A model based on sediment properties with and without gas hydrate

Gas hydrate is commonly found concentrated in sand- and silt-rich reservoirs and these are thought to have great potential as an energy resource. To evaluate this potential, it is critical to understand their in-situ stress state and their behavior during production. We explore the geomechanical behavior of sandy silt sediments from a hydrate reservoir at Green Canyon GC 955, deep-water Gulf of Mexico without hydrate and compare these results to previous work on the same sandy silty sediments with hydrate present. Previous results from hydrate-bearing sediments at GC 955 using intact specimens show that the ratio of lateral to axial effective stress under uniaxial strain (K0) is strongly time-dependent and over long timescales approaches 1 (isostatic stress state). New results from sediments without hydrate using reconstituted specimens reveal that K0 is 0.51 and independent of time and stress. The reservoir without hydrate is approximately three times more compressible than with hydrate. In addition, we find that hydrate-free sediments have a normalized creep rate of Cα/Cc of about 0.028 at all stress levels, whereas hydrate bearing sediments exhibit a decrease in Cα/Cc with stress. We use these observations to present a new geomechanical model that describes how stress and porosity evolve during production of the hydrate reservoir. The key insight in our proposed model is that the sediment skeleton and gas hydrate phases share any applied external load in proportion to their relative stiffnesses. The model results predict that hydrate-bearing sediments are stiffer than gas hydrate-free sediments and that K0 reduces dramatically after gas hydrate dissociation, as is observed. At GC 955, a 5 MPa depressurization is needed to cause gas hydrate dissociation. This change in pore pressure increases the vertical and horizontal effective stresses from 3.8 to 8.8 MPa and 3.8–4.4 MPa, respectively. The dramatic vertical effective stress increase will cause significant compaction (7%) and perturb stresses near the wellbore. These experimental observations and the load-sharing model presented have the potential to underpin a new generation of hydrate reservoir simulation models, which may drive more effective production approaches.

Rolling table centrifuge modelling of partially saturated granular material to inform on instability during solid bulk cargo transport

AbstractThe response of partially saturated granular cargoes during maritime transportation has resulted in the capsize and sinking of 27 bulk carriers at sea and the loss over 90 seafarers’ lives in the last decade. The partially saturated granular material response to energy imparted during cargo loading, ship engine vibrations and vessel rolling motions from sea states causes a change in state of the granular cargo, which can lead to vessel instability and ultimately capsize. However, the mechanisms driving the response of partially saturated granular cargos within a bulk carrier hold are not well understood. This paper presents results from an experimental study of rolling table centrifuge model tests on a partially saturated silica sand to contribute to improved understanding of the response of granular cargoes during maritime transport. Observed settlement, pore pressure, moisture content and density changes during and/or following a sequence of large amplitude rolling motions are presented. The results indicate that for the conditions considered, a progressive upwards migration of pore water during rolling led to creation of a free surface of water above the granular sample that was left in a denser, lower moisture content sample compared to its initial state. Sloshing of free water on top of even a competent cargo during rolling motions of the ship can contribute to loss of ship stability and ultimately capsize.

Application of Compressible Carbon in Cement to Mitigate Cement Pore Pressure Reduction and Improve Cement Bonding

Summary A quality cement job is essential to ensure long-term zonal isolation in oil/gas and carbon storage wells. This may be affected by the cement hydration process, during which water is consumed, causing pore pressure reduction and shrinkage, which increases the risk of fluid/gas invasion, formation of microannuli, and even the risk of shear bond strength reduction. In this paper, we present the application of compressible carbon (CC) particles as a cement additive to mitigate cement pore pressure reduction and improve zonal isolation. We designed cement formulations with CC and evaluated them using API standard tests and industry-wide recognized tests to ensure they met functional requirements for well cementing, such as mixability, thickening time, rheology, compressive strength development, fluid loss, and free water. The concept is as follows: When cement is pumped into an annulus, the CC particles will compress under hydrostatic pressure. When the cement is hydrating, the CC will expand to mitigate the pore pressure reduction. To validate the proof of concept, we developed a benchtop high-pressure, high-temperature setup and a pilot-scale test setup. A cement formulation with 9% CC by volume was found to be stable and have controllable performance properties such as thickening time, fluid loss, and free water. We also designed and tested a control slurry without CC for comparison. Both the benchtop and pilot-scale tests demonstrated that adding CC into the cement formulation mitigated the pore pressure reduction during cement hydration. This may reduce the risk of fluid/gas invasion that could result in migration in the set cement. Comparing the cement containing carbon with the control system without carbon, the gas permeability was reduced. The bond strength obtained from shear bond tests improved significantly, which may reduce the risk of debonding and be an indication of the reduction of shrinkage. In addition, a better hydraulic seal was found in one test, but further investigation is needed. Similar mechanical properties such as compressive strength were measured for both cements with and without CC. The data showed that adding CC does not have a negative impact on the hardened cement properties. In fact, adding CC decreased Poisson’s ratio slightly. In this paper, we present the results of proof-of-concept testing for the novel application of CC in cement to improve zonal isolation by mitigating pore pressure reduction. We also present new benchtop and pilot-scale experimental setups to measure pore pressure changes during cement hydration.

Geomechanics analysis of caprock integrity in injection and production operations of aquifer gas storage

ABSTRACT The safety of aquifer gas storage (AGS) depends on the mechanical stability of the caprock to a great extent. The injection and production operation of AGS will cause the disturbance of the in-situ stress due to the change of pore pressure in the reservoir and caprock, which may lead to the failure of caprock. The con AGS is established based on geological structure and hydrogeological information. By using the finite element method, the factors such as boundary condition, properties of reservoir and caprock, and in-situ stress coefficient are considered, respectively, and the variation of formation pressure and disturbed stress in the process of injection and production are analyzed. The results show that the boundary conditions are of significant influence on the stress field. The shear failure risk of reservoir–caprock interface under closed boundary is the greatest. The in-situ stress coefficient has the greatest influence on the failure risk of the caprock, followed by the permeability of the reservoir and caprock.

Vulnerability of firn to hydrofracture: poromechanics modeling

Abstract On the Greenland Ice Sheet, hydrofracture connects the supraglacial and subglacial hydrologic systems, coupling surface runoff dynamics and ice velocity. In recent decades, the growth of low-permeability ice slabs in the wet snow zone has expanded Greenland's runoff zone, but observations suggest that surface-to-bed connections are rare, because meltwater drains through crevasses into the porous firn beneath ice slabs. However, there is little quantitative evidence confirming the absence of surface-to-bed fracture propagation. Here, we use poromechanics to investigate whether water-filled crevasses in ice slabs can propagate vertically through an underlying porous firn layer. Based on numerical simulations, we develop an analytical estimate of the water injection-induced effective stress in the firn given the water level in the crevasse, ice slab thickness, and firn properties. We find that the firn layer substantially reduces the system's vulnerability to hydrofracture because much of the hydrostatic stress is accommodated by a change in pore pressure, rather than being transmitted to the solid skeleton. This result suggests that surface-to-bed hydrofracture will not occur in ice slab regions until all pore space proximal to the initial flaw has been filled with solid ice.

Theoretical Analysis of Drilling Unloading and Pile-Side Soil Pressure Recovery of Nonsqueezing Pipe Piles Installed in K0-Consolidated Soils

Drilling with prestressed concrete (DPC) pipe pile is a nonsqueezing pile sinking technology, employing drilling, simultaneous pile sinking, a pipe pile protection wall, and pile side grouting. The unloading induced by drilling, the pipe pile supporting effect, and the dissipation of the negative excess pore-water pressure after pile sinking, all of which have significant effects on the recovery of soil pressure on the pile side, are the main concerns of this study, which aim to establish a method to reasonably evaluate the timing selection of pile side grouting. The theoretical solutions for characterizing the unloading and dissipation of the negative excess pore-water pressure are presented based on the cylindrical cavity contraction model and the separated variable method. By inverse-analyzing the measured initial pore pressure change data from borehole unloading, initial soil pressures on the pile side of each soil layer are determined using the presented theoretical solutions. Then, the presented theoretical solutions were verified through a comparative analysis with the corresponding measured results. Moreover, by introducing time-dependent coefficients αt1 and αt2 to characterize the pore pressure dissipation and rheology effects, the effects of the negative excess pore-water pressure dissipation on the pile-side soil pressure recovery are discussed in detail.

Further development of distinct lattice spring model with Finn model for liquefaction analysis

The hazards associated with sand liquefaction induced by dynamic events are significant. The study of the dynamic stability in hydraulic structures presents an interdisciplinary challenge encompassing both geotechnical engineering and engineering seismology. This study was based on an actual hydraulic engineering project. To effectively predict changes in pore pressure caused by earthquakes, we integrated the Finn constitutive equations into a distinct lattice spring model (DLSM). In this study, the code was customized to accommodate multiple materials simultaneously participating in the calculations, thus simplifying the solution to complex engineering problems. Initially, we validated the DLSM’s liquefaction equation by comparing it with the finite difference method. Subsequently, we conducted a comparative analysis of liquefaction in an engineering project of sluice using the enhanced DLSM. Our analysis indicates that untreated sand has a severe risk of trending toward liquefaction, presenting a hazard to hydraulic engineering. The incorporation of a gridded concrete framework significantly mitigated the seismic-induced pore pressure accumulation and irreversible deformations caused by vibrations. A comparative study showed that concrete retaining walls with concrete supports are more effective at reducing liquefaction hazards and minimizing irreversible deformations in engineering structures.

Cyclic Injection Leads to Larger and More Frequent Induced Earthquakes under Volume-Controlled Conditions

Abstract As carbon storage technologies advance globally, methods to understand and mitigate induced earthquakes become increasingly important. Although the physical processes that relate increased subsurface pore pressure changes to induced earthquakes have long been known, reliable methods to forecast and control induced seismic sequences remain elusive. Suggested reservoir engineering scenarios for mitigating induced earthquakes typically involve modulation of the injection rate. Some operators have implemented periodic shutdowns (i.e., effective cycling of injection rates) to allow reservoir pressures to equilibrate (e.g., Paradox Valley) or shut-in wells after the occurrence of an event of concern (e.g., Basel, Switzerland). Other proposed scenarios include altering injection rates, actively managing pressures through coproduction of fluids, and preinjection brine extraction. In this work, we use 3D physics-based earthquake simulations to understand the effects of different injection scenarios on induced earthquake rates, maximum event magnitudes, and postinjection seismicity. For comparability, the modeled injection considers the same cumulative volume over the project’s operational life but varies the schedule and rates of fluid injected. Simulation results show that cyclic injection leads to more frequent and larger events than constant injection. Furthermore, with intermittent injection scenario, a significant number of events are shown to occur during pauses in injection, and the seismicity rate remains elevated for longer into the postinjection phase compared to the constant injection scenario.

A thermo-hydromechanical damage model and its application to a deep geological radioactive repository

This study conducts a numerical analysis of long-term thermo-hydromechanical (THM) processes, with a particular focus on deep geological disposal of radioactive waste. It emphasizes the modeling of damage zones resulting from excavation activities, as well as changes in pore pressure and temperature. The numerical model presented in this paper delineates the fundamental relations of thermo-poro-elasticity. It also introduces a double phase field model specifically formulated for rock materials under compressive stress, taking into account THM coupling and time-dependent behavior. The proposed model has been implemented in a finite element code and applied to numerical analysis of short and long terms responses around a horizontal borehole in the context of the French High-Level Waste (HLW) disposal concept. In particular, displacements, temperature changes, pore fluid pressure variations as well as evolution of induced damage zones are presented and analyzed for different periods including excavation, operation and post-closure. It is found that the occurrence of damage zones is clearly related to the spacing distance between two neighboring disposal boreholes. The time-dependent deformation of host rock plays an significant role in the long-term responses.