Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pore pressure”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulation study of chemo-mechanical impacts of CO 2 injection in morrow b sandstone reservoir

Here, this study investigates the impacts of geomechanical and geochemical changes on carbon storage in a partially depleted oil reservoir, using results from four different coupled simulation models. Models were used to examine the relative importance of storage mechanisms, and how changing reservoir parameters might affect these mechanisms through time. The study uses data from a Morrowan sandstone reservoir in the Farnsworth Unit (FWU), Ochiltree County, Texas which is currently undergoing CO 2 enhanced oil recovery (EOR). Partially depleted oil reservoirs such as the FWU offer attractive carbon utilization and/or storage targets because of existing infrastructure and economic benefits from incremental oil recovery as well as tax credits. However, prediction of storage capacity or long-term fluid migration in these fields can be difficult because of the wide variation in formation fluids and operational histories that may have undergone. CO 2 injection can cause complex geomechanical and geochemical responses in a reservoir as a result of interplay between dynamic changes in pore pressure, reservoir temperature, fluid composition, and interactions between formation fluids, CO 2 , and reservoir rock. Thus, multiple coupled numerical simulation models must be developed and used to more precisely understand what CO 2 storage mechanisms are most significant, as well as the long-term fate of the stored CO 2 . Our study used results from hydrodynamic, coupled hydro-geomechanical, coupled hydro-geochemical, and coupled hydro-geomechanical-geochemical models to examine how changes in geomechanical and geochemical properties can impact the injectivity or storage capacity of CO 2 . Models simulated historical field operations and then forward-modeled a water-alternate-gas (WAG) operation for 20 years, followed by a 1000-year post-injection monitoring. The work demonstrates that in this specific reservoir, geomechanical impacts appear to be more significant than any geochemical processes, which showed minimal impact on reservoir properties. Furthermore, these simulation scenarios corroborated other studies showing more than 90% of CO 2 injected are being stored and will remain stored in the reservoir. The study was able to draw on a rich set of data from this currently operational project for the work. Because the Morrow is a widespread oil-bearing reservoir in the midcontinent of the United States and is fairly typical of many sandstone reservoirs elsewhere, the study has significance beyond the boundaries of our field area.

58 GEOSCIENCES↗
DOE OSTI2021

Multi-scale Coupled Processes Modeling of Fractures as Porous, Interfacial and Granular Systems from Rock Images with the Numerical Manifold Method

Abstract The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

58 GEOSCIENCES↗
DOE OSTI2025

Thermophysical properties of the Ghareb formation relevant for nuclear waste disposal

The Ghareb formation, a shallowly buried porous chalk in Israel, is currently a candidate for nuclear waste disposal. The potential repository is somewhat unique for its host rock and emplacement in shallow (500 m) engineered large diameter boreholes. Herein, the thermal properties of the Ghareb are determined to support design and performance assessment; the relevant properties measured are thermal conductivity, specific heat capacity, thermal diffusivity, and thermal expansion coefficient along with their relationship with varying temperature. For the temperature range of 40 to 275 °C, the thermal conductivity ranges from 0.30 to 1.10 W/m·K, the thermal diffusivity ranges from 0.20 to 0.72 mm 2 /s, and volumetric heat capacity ranges from 0.86 to 2.00 MJ/m 3 ·K. Thermal strain measurements were used to estimate the linear thermal expansion coefficient to be 6·10 –4 –9·10 –2 °C -1 from 40 to 300 °C. These measured properties were used in a thermomechanical model to estimate near-field stresses an hour and 10 years after waste emplacement; the borehole was found to be stable. Thermal loading after 10 years was predicted to elevate local pore pressures by 1–1.5 MPa. The laboratory measurements coupled with analyses are the first attempts at performance assessment characterization for this first of its kind potential repository setting with this chosen host rock.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Investigation of methane mass transfer and sorption in Marcellus shale under variable net-stress

Natural gas in shale exists as free and adsorbed gas, subject to prevailing pore pressures and stress conditions. Accordingly, to accurately estimate/predict the shale gas recovery potential, a central requirement is to represent gas transport and sorption behavior under varying stress conditions. The objective of this work is to facilitate the interpretation of laboratory-scale experiments, at relevant conditions, in an attempt to bridge the gap in scales between laboratory- and field-scale observations. We have conducted a series of high-pressure experiments on a full-diameter core sample from the Marcellus shale. These include gas loading (pressure-decay) and depletion (production) experiments with pure methane (CH 4 ) at variable stress conditions to characterize transport and sorption behavior under reservoir-relevant conditions. Here we have formulated and applied a novel integral model for mass transfer and storage in multi-porosity shale systems that allows us to effectively investigate transport and sorption phenomena: We delineate gas transport by interpreting helium (He) pressure-decay experiments and demonstrate how to use the information gained to calculate the relevant transport coefficients of CH 4 and other gases. A separate measurement of the CH 4 sorption isotherm on a smaller sample (a shale cube) was interpreted and combined with the transport description to predict the production behavior of CH 4 from the experiments with the full-diameter core. Our experiments demonstrate that the representation of sorption hysteresis is crucial for predicting and guiding shale gas production: At the end of both gas production experiments, approximately 20% of the initial gas in place remained in the core. Without accounting for sorption hysteresis, our modeling demonstrates that the CH 4 production could be overestimated by 10%. We demonstrate that our integral, triple-porosity model provides an effective approach for the interpretation and prediction of gas transport and sorption behavior during loading and production experiments on shale cores under variable net-stress conditions. In summary, our work combines measurements and modeling of mass transfer and sorption in shales at different scales to validate a characterization approach that facilitates an improved understanding of shale gas production. Furthermore, the triple-porosity model utilized in our work defines a potential pathway for the translation of laboratory-scale experimentation to larger-scale applications.

58 GEOSCIENCES↗
DOE OSTI2023

The physics of desiccation cracks 2: Modeling and prediction of the crack patterns

Here, this paper extends the experimental results of our companion paper by modeling and predicting the onset and pattern formation of desiccation cracks in geomaterials (Ruoyu et al., 2023). Thin-layer silt samples in controlled atmospheric conditions were tested to obtain the surface strain maps with the digital image correlation (DIC) method during dehydration. Support experiments, including consolidation and displacement-controlled triaxial tests, were conducted for the properties of geomaterials. These experimental results were used to validate a viscoplastic theoretical model by comparing cracks and singularities locations that distribute following the Cnoidal wave pattern. A critical value exists in the viscoplastic model that determines the number of singularities, which accurately predicts the number of cracks in the experiment. The critical values contain two crucial parameters: rate sensitivity and λ. Rate sensitivity describes the rate-dependent stress–strain relation, while λ shows the ability of pore pressure redistribution under the external mechanical loading rate. These results provide a new view to analyze the desiccation cracks considering the rate-dependent viscoplasticity.

42 ENGINEERING↗
DOE OSTI2021

Inference of Rock Flow and Mechanical Properties from Injection-Induced Microseismic Events During Geologic CO 2 Storage

Monitoring microseismic activities during CO 2 injection into geologic formations is important for ensuring the safety of the storage operations. The resulting data provide insight into the response of the storage formation to CO 2 injection and can be used to infer the underlying rock flow and mechanical properties. In this paper, assimilation of microseismic data is performed for dynamic characterization of the storage formation by using a stochastic simulation model to forecast the microseismic response of a geologic formation during CO 2 injection. Two modeling approaches are adopted to predict the space-time distribution of the injection-induced microseismicity. The first model is based on pore pressure relaxation assumption, while the second model uses coupled flow and geomechanics simulation to establish the complex physical relation between the storage formation properties and the corresponding microseismic responses during CO 2 injection. The stochastic predictive models in each case are used in ensemble data assimilation frameworks to estimate rock properties from the observed microseismic data. Two data assimilation methods are considered: (i) a new ensemble-based stochastic point process filter (EnPPF) that can directly integrate discrete microseismic events, and (ii) a variant of ensemble smoother, known as the ensemble smoother with multiple data assimilation (ES-MDA), which requires continuous representation of microseismic events for assimilation. The two methods are successfully applied to a geologically realistic model of the Farnsworth Field in Texas, with complex geologic flow units and interacting fault systems.

42 ENGINEERING↗
DOE OSTI2022

Assessment of oil and gas fields in California as potential CO 2 storage sites

California's total annual greenhouse gas (GHG) emissions (425.3 MtCO 2 e) in 2018 were about 6.4% of the US total (6,677 MtCO 2 e) and around 1% of global emissions. About 39% of 2018 GHG emissions in California were from the industrial and electrical sectors. Many of these emissions were from large stationary point sources and were suitable for carbon capture retrofit with subsequent storage of the captured carbon dioxide (CO 2 ) in geological formations. Previous studies of California found suitable geology and CO 2 storage resource. This study refines and furthers prior work using a three-stage screening process of oil fields, gas fields, and underground natural gas storage (UGS) sites by combining criteria from previous studies while excluding sites that pose technical risk or are located in regions with surface restrictions including sensitive habitats and dense populations. In the first stage, 129 CO 2 storage sites in California were identified using qualification criteria based upon formation properties including geological conditions and pore pressure. The second stage identified sensitive sites by applying conservative screens including seismic activity, faulting, population density, restricted lands, and sensitive habitats. During the third stage, 61 CO 2 potential storage sites were identified by subtraction of stage 2 areas from stage 1. The potential storage volume in the third stage ranged from 1.0 to 2.0 GtCO 2 . Finally, we applied a scoring system with seven parameters to rank the 61 potential sites based on subsurface technical criteria. The scored sites are classified as high priority, medium priority, and sites for future study. Prospective CO 2 storage sites with high and moderate priority were selected and linked to CO 2 sources. There are 14 prospective sites (above 20 MtCO 2 storage resource per site) with a total storage resource of 1024 MtCO 2 distributed in Northern and Southern California. Of these sites, there are 9 potential CO 2 -EOR sites and 1 depleted oil field with a total estimated CO 2 storage volume of ~800 MtCO 2 in the Southern San Joaquin and Ventura Basin. These 10 prospective sites with a storage resource greater than 20 MtCO 2 could potentially deliver more than 20 years of storage with an average injection rate of 40 MtCO 2 /year. The remaining 4 highly prospective sites are in Northern California. Additionally, study results suggest that saline formations should be re-evaluated in concert with storage in oil, gas, and natural gas storage reservoirs.

58 GEOSCIENCES↗
DOE OSTI2021

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗
DOE OSTI2024

De-risking fault leakage risk and containment integrity for subsurface storage applications

The subsurface is pivotal in the energy transition, for the sequestration of CO 2 and energy storage. It is crucial to understand to what extent geological faults may form leakage pathways that threaten the containment integrity of these projects. Fault flow behavior has been studied in the context of hydrocarbon development, supported by observations from wells drilled through faults, but such observations are rare in geoenergy projects. Focusing on mechanical behavior as early indicator of potential leakage risks, a probabilistic Coulomb Failure Stress workflow is developed and demonstrated using data from the Decatur CO 2 sequestration project to rank faults based on their containment risk. The analysis emphasizes the importance of fault throw relative to reservoir thickness and pore pressure change in assessing reactivation risks. Integrating this mechanical assessment with geological and dynamic fault analyses contributes to derisking fault containment for geoenergy applications, providing valuable insights for the successful development of subsurface storage projects.

58 GEOSCIENCES↗
DOE OSTI2022

Relationships between fault characteristics and seismic responses in a large lab-scale tri-axial injection test conducted on a faulted Castlegate Sandstone

Here to investigate mechanisms causing microseismicity (M w < 2) at a CO 2 injection site, a large-scale triaxial block experiment was carried out on a faulted (saw-cut) cubic-meter of Castlegate Sandstone. The experiment consisted of injection tests at varying differential stresses, while monitoring and recording pore pressure and acoustic emissions (AEs). During the experiment, ~33,000 AEs and ~14 mm of horizontal displacement/slip, like a strike-slip fault movement, occurred. To understand the AE responses and ascertain fault characteristics near the located AEs, we modeled the topography of the fault surface, fault aperture, and fault-gouge thickness using pre- and post-experiment laser scans of the fault surface on each half of the block. Additionally, we characterized surface roughness parallel and perpendicular to slip. Models show crushing and flattening of the fault surfaces can be linked to the spatiotemporal distribution of AEs within 50 mm of the fault surface. Approximately 65% of AEs were in areas with small aperture (≤300 μm); thicker fault gouge was observed in adjacent areas with wider aperture and shows a two-fold reduction in grain size relative to unaltered Castlegate Sandstone. This work provides a conceptual understanding on fault surface evolution, which can be applied towards modeling of seismic slip.

58 GEOSCIENCES↗
DOE OSTI2022

Porosity evolution during heating of copper made from powder by friction extrusion

Here, friction extrusion was used to compact and extrude solid copper rod from feedstock powders. Following extrusion, considerable porosity was observed throughout the extrudate cross section due to the entrained porosity in the feedstock material and the extrusion process. The thermal stability of the extrudate was investigated via a series of heat treatments. Porosity evolution exhibited three distinct stages—an unchanged plateau (0°C–300°C) followed by an abrupt increase (400°C–500°C) and ultimately a reduction, as the temperature increases (>500°C). The peak porosity measured was ~25%. The underlying driving force for pore evolution is described as the competition between the internal pore pressure, material strength, and sintering kinetics, as a function of temperature. The observed porosity evolution and driving force are not expected to be limited to copper. Thus, this manuscript reveals an important consideration regarding elevated-temperature thermal stability as advanced manufacturing methods involving direct powder extrusion are explored.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Alteration in micro-mechanical characteristics of wellbore cement fracture surfaces due to fluid exposure

The stiffness of wellbore cement fracture surfaces was measured after exposing to the advective flow of nitrogen, silicone oil, and medium sweet dead crude oil for different exposure periods. The test specimens were extracted from fractured cement cylinders, where the cement fracture surfaces were exposed to the different fluids up to 15 weeks. A nanoindenter with a Berkovich indenter tip was used to measure load-indentation depth data, which was used to extract the elastic modulus (E) and nano-hardness (H) of the cement fracture surfaces. A reduction in the elastic modulus compared with an unexposed specimen were observed in all the specimens. Both elastic modulus and nano-hardness for the specimens exposed to silicone oil were lower than specimens exposed to nitrogen gas and varied with the period of exposure. The elastic modulus and nano-hardness of the specimens exposed to crude oil were the lowest with a significant decrement with the exposure period. The frequency distribution of the nanoindentation measurements shows that the volume-fraction ratio of the two types of cement hydrated nanocomposites for both the unexposed and test specimens is about 70:30%. Phase transformation beneath the indenter is observed for all of the specimens, with more obvious plastic deformation in specimens exposed to crude oil. Analytical measurements (SEM, EDS, FT-IR, and XRD) on exposed cement fracture surfaces reveal different levels of physical and chemical alteration that are consistent with the reduction in stiffness measured by nanoindentation. Furthermore, the study suggests that cement stiffness will decrease due to crude oil exposure, and the fracture will be sensitive to stress and pore pressure with time.

42 ENGINEERING↗
DOE OSTI2024

Finite domain solution of a KGD hydraulic fracture in the viscosity-dominated regime

This paper describes a numerical algorithm for solving the classic problem of a plane strain (KGD) fracture propagating in an impermeable elastic medium with zero toughness. The method, which takes advantage of the self-similar nature of the solution, combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation. This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects, noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method. To enhance the efficiency and accuracy of the numerical scheme, the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening, while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation. The proposed technique yields results that closely match the analytical solution, even with a coarse mesh. This approach offers potential for addressing more complex hydraulic fracturing problems in the future.

Domain-based method↗
DOE OSTI2025

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN↗
DOE OSTI2021

Potential Link Between 2020 Mentone, West Texas M5 Earthquake and Nearby Wastewater Injection: Implications for Aquifer Mechanical Properties

Abstract The M5 Mentone earthquake that occurred on March 26, 2020, was the largest event recorded over the last 2 decades in West Texas within the Delaware Basin, a U.S. major petroleum‐producing area. Also, numerous hydrofracturing and wastewater disposal wells are spread across this region. Within a 30 km distance to mainshock, eight class‐II injection wells for industrial wastewater disposal target the deep porous Ellenburger aquifer at an average rate of 1.36 × 10 6 barrel (BBL) per month during 2012–2020. Poroelastic models of fluid diffusion show these nearby injectors collectively imparted up to 80.5 kPa of Coulomb stress at the mainshock location, capable of triggering this M5 event. Assuming the Mentone event occurs when pore‐pressure increase is maximum, the time delay between peak injection and the M5 occurrence corresponds with an optimal permeability of 6.76 × 10 ‐14 m 2 for the Ellenburger aquifer layer, in agreement with independent estimates.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Seismicity and Stress Associated With a Fluid-Driven Fracture: Estimating the Evolving Geometry

A coupled approach, combining the theory of rate- and state-dependent friction and methods from poroelasticity, forms the basis for a quantitative relationship between displacements and fluid leak-off from a growing fracture and changes in the rate of seismic events in the region surrounding the fracture. Poroelastic Green's functions link fracture aperture changes and fluid flow from the fracture to changes in the stress field and pore pressure in the adjacent formation. The theory of rate- and state-dependent friction provides a connection between Coulomb stress changes and variations in the rate of seismic events. Numerical modeling indicates that the Coulomb stress changes can vary significantly between formations with differing properties. The relationship between the seismicity rate changes and the changes in the formation stresses and fluid pressure is nonlinear, but a transformation produces a quantity that is linearly related to the aperture changes and fluid leak-off from the fracture. The methodology provides a means for mapping changes in seismicity into fracture aperture changes and to image an evolving fracture. An application to observed microseismicity associated with a hydrofracture reveals asymmetric fracture propagation within two main zones, with extended propagation in the upper zone. Furthermore, the time-varying volume of the fracture agrees with the injected volume, given by the integration of rate changes at the injection well, providing validation of the estimated aperture changes.

58 GEOSCIENCES↗
DOE OSTI2021

Aquifer Monitoring Using Ambient Seismic Noise Recorded With Distributed Acoustic Sensing (DAS) Deployed on Dark Fiber

Abstract Groundwater is a critical resource for human activities worldwide, and a vital component of many natural ecosystems. However, the state and dynamics of water‐bearing aquifers remain uncertain, mostly due to the paucity of subsurface data at high spatial and temporal resolution. Here, we show that analysis of infrastructure‐generated ambient seismic noise acquired on distributed acoustic sensing (DAS) arrays has potential as a tool to track variations in seismic velocities (d v / v ) caused by groundwater level fluctuations. We analyze 5 months of ambient noise acquired along an unused, 23 km‐long telecommunication fiber‐optic cable in the Sacramento Valley, CA, a so‐called “dark fiber." Three array subsections, ∼6 km apart, are processed and the stretching technique is applied to retrieve daily d v / v beneath each location. Near the Sacramento river, d v / v variations in the order of 2%–3% correlate with precipitation events and fluctuations in river stage of ∼1.5 m. In contrast, regions away (2.5 km) from the river do not experience large d v / v variations. These observations reveal short‐scale spatial variability in aquifer dynamics captured by this approach. Dispersion analysis and surface wave inversion of noise gathers reveal that seismic velocity perturbations occur at depths of 10–30 m. Rock physics modeling confirms that observed d v / v are linked to pore pressure changes at these depths, caused by groundwater table fluctuations. Our results suggest that DAS combined with ambient noise interferometry provides a means of tracking aquifer dynamics at high spatial and temporal resolutions at local to regional scales, relevant for effective groundwater resource management.

58 GEOSCIENCES↗
DOE OSTI2021

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗