In this study, we analyze a deep hydraulic stimulation of a fracture zone that was conducted as part of the classical Fenton Hill Hot Dry Rock Program in the 1970s. At the time, it was suggested that a pre-existing fracture or multiple fractures within the fracture zone were jacked open by injection-induced increase in pressure. In this study, we analyze the same stimulation experiment, but investigate the possibility of an alternative mechanism of shear reactivation of pre-existing fractures. We conduct modeling that accounts for both jacking (or elastic-fracture opening) and shear-slip dilation and demonstrate that injection-induced shear reactivation (or hydroshearing) could have occurred simultaneously with seismic events of magnitudes lower than what can be felt by humans. In fact, simulations considering shear reactivation seem to better match observed fluid recovery after multiple injection cycles. Shear reactivation and shear dilation results in locked-open fractures, especially near the injection well that provides permeability of higher flow recovery. We then investigate the sensitivity of the proposed model by varying some of the critical parameters such as maximum aperture, dilation angle, as well as fracture density. Interestingly, none of the simulated cases resulted in a large event that could have been felt by humans, but did result in a cumulative seismic magnitude of less than 1 for each given stimulation step. These results suggest that a permanent irreversible permeability increase of several orders of magnitude can be obtained by hydroshearing in a “seismically” safe manner.

Understanding the injection-induced triggering mechanism is a fundamental step towards controlling the seismicity generated by deep underground exploitation. Here we propose a modeling approach based on coupling the TOUGH2 simulator with a geomechanical-stochastic model. The hydro-mechanical-stochastic model provides a good representation of different mechanisms influencing each other during and after the injection phase. Each mechanism affects the induced seismicity in a different way and at different times during the reservoir stimulation, confirming that a complex interaction is in place, and that more sophisticated and physics-based approaches coupled with statistical model are required to explain such a complex interaction. In addition to previous statistical and hybrid models, our approach accounts for a full 3D formulation of both stresses and fluid flow, further including all the TOUGH2 capabilities. Furthermore, it includes interactions between triggered seismic events through calculation of static stress transfer. In this work, we present the main capabilities of TOUGH2-SEED and apply the model to the Basel EGS case, successfully reproducing the injection pressure as well as the evolution of the seismicity.

Carbon dioxide (CO2) is likely to reach the bottom of injection wells at a colder temperature than that of the storage formation, causing cooling of the rock. This cooling, together with overpressure, tends to open up fractures, which may enhance injectivity. We investigate cooling effects on injectivity enhancement by modeling the In Salah CO2 storage site and a theoretical, long-term injection case. We use stress-dependent permeability functions that predict an increase in permeability as the effective stress acting normal to fractures decreases. Normal effective stress can decrease either due to overpressure or cooling. We calibrate our In Salah model, which includes a fracture zone perpendicular to the well, obtaining a good fitting with the injection pressure measured at KB-502 and the rapid CO2 breakthrough that occurred at the observation well KB-5 located 2 km away from the injection well. CO2 preferentially advances through the fracture zone, which becomes two orders of magnitude more permeable than the rest of the reservoir. Nevertheless, the effect of cooling on the long-term injectivity enhancement is limited in pressure dominated storage sites, like at In Salah, because most of the permeability enhancement is due to overpressure. However, thermal effects enhance injectivity in cooling dominated storage sites, which may decrease the injection pressure by 20%, saving a significant amount of compression energy all over the duration of storage projects. Overall, our simulation results show that cooling has the potential to enhance injectivity in fractured reservoirs.

The spectral boundary integral method is popular for simulating fault, fracture, and frictional processes at a planar interface. However, the method is less commonly used to simulate off-fault dynamic fields. Here we develop a spectral boundary integral method for poroelastodynamic solid. The method has two steps: first, a numerical approximation of a convolution kernel and second, an efficient temporal convolution of slip speed and the appropriate kernel. The first step is computationally expensive but easily parallelizable and scalable such that the computational time is mostly restricted by computational resources. The kernel is independent of the slip history such that the same kernel can be used to explore a wide range of slip scenarios. We apply the method by exploring the short-time dynamic and static responses: first, with a simple source at intermediate and far-field distances and second, with a complex near-field source. We check if similar results can be attained with dynamic elasticity and undrained pore-pressure response and conclude that such an approach works well in the near-field but not necessarily at an intermediate and far-field distance. We analyze the dynamic pore-pressure response and find that the P-wave arrival carries a significant pore pressure peak that may be observed in high sampling rate pore-pressure measurements. We conclude that a spectral boundary integral method may offer a viable alternative to other approaches where the bulk is discretized, providing a better understanding of the near-field dynamics of the bulk in response to finite fault ruptures.

Thermo-mechanical effects are important in geologic carbon storage because CO2 will generally reach the storage formation colder than the rock, inducing thermal stresses. Capillary functions, i.e., retention and relative permeability curves, control the CO2 plume shape, which may affect overpressure and thus, caprock stability. To analyze these thermal and capillary effects, we numerically solve non-isothermal injection of CO2 in deformable porous media considering the In Salah, Algeria, CO2 storage site. We find that changes in the capillary functions have a negligible effect on overpressure and thus, caprock stability is not affected by capillary effects. However, we show that for the strike slip stress regime prevalent at In Salah, stability decreases in the lowest parts of the caprock during injection due to cooling-induced thermal stresses. Simulations show that shear slip along pre-existing fractures may take place in the cooled region, whereas tensile failure is less likely to occur. Indeed, only the injection zone and the lowest tens of meters of the 900-m-thick caprock at In Salah might be affected by cooling effects, which would thus not jeopardize the overall sealing capacity of the caprock. Furthermore, faults are likely to remain stable far away from the injection well because outside the cooled region the injection-induced stress changes are not sufficient to exceed the anticipated shear strength of minor faults. Nevertheless, we recommend that thermal effects should be considered in the site characterization and injection design of future CO2 injection sites to assess caprock stability and guarantee a permanent CO2 storage.

This work investigates the injection-induced seismic response of a heterogeneous fault plane, featuring low-permeability asperities embedded into a high-permeability damage zone. We simulate the pressure evolution with a hydrogeological simulator, accounting for the heterogeneous fault plane. Seismicity occurs then on the asperities, represented as unstable patches reactivating by means of the Mohr-Coulomb criterion. The hydrological and seismic modules are implicitly coupled to account for effects of asperity reactivation on the permeability. Results show that permeability changes may cause at a later time a change in seismicity propagation. We also investigated such effects by varying the density of asperities.

Induced seismicity due to natural gas production is observed at different sites worldwide. Common understanding states that the pressure drop caused by gas production leads to compaction, which affects the stress field in the reservoir and the surrounding rock formations and hence reactivates preexisting faults and induces earthquakes. In this study, we show that the multiphase fluid flow involved in natural gas extraction activities should be included. We use a fully coupled fluid flow and geomechanics simulator, which accounts for stress-dependent permeability and linear poroelasticity, to better determine the conditions leading to fault reactivation. In our model setup, gas is produced from a porous reservoir, divided into two compartments that are offset by a normal fault. Results show that fluid flow plays a major role in pore pressure and stress evolution within the fault. Fault strength is significantly reduced due to fluid flow into the fault zone from the neighboring reservoir compartment and other formations. We also analyze scenarios for minimizing seismicity after a period of production, such as (i) well shut-in and (ii) gas reinjection. In the case of well shut-in, a highly stressed fault zone can still be reactivated several decades after production has ceased, although on average the shut-in results in a reduction in seismicity. In the case of gas reinjection, fault reactivation can be avoided if gas is injected directly into the compartment under depletion. However, gas reinjection into a neighboring compartment does not stop the fault from being reactivated.

In 2013, fluid injection during the St. Gallen deep geothermal project, Switzerland, induced hundreds of seismic events, including a ML 3.5 earthquake on a fault hundreds of meters away from the well. Recent studies have suggested the direct pressure effect through permeable hydraulic connections and poroelastic effects as possible mechanisms for inducing seismicity on distant faults. In St. Gallen, operational, seismic, and earthquake data are available to investigate the underlying physical mechanisms using a numerical model. The results show that Coulomb stress changes at the fault can be 3 orders of magnitude greater when a hydraulic connection is present. Combining this with several field observations, we conclude that the direct pressure effect was more likely the predominant mechanism behind the seismicity induced in St. Gallen. The detection of hydraulic connections may be important for future projects as pressure can be driven far from the well and reactivate remote faults.

Ground deformation, commonly observed in storage projects, carries useful information about processes occurring in the injection formation. The Krechba gas field at In Salah (Algeria) is one of the best-known sites for studying ground surface deformation during geological carbon storage. At this first industrial-scale on-shore CO2 demonstration project, satellite-based ground-deformation monitoring data of high quality are available and used to study the large-scale hydrological and geomechanical response of the system to injection. In this work, we carry out coupled fluid flow and geomechanical simulations to understand the uplift at three different CO2 injection wells (KB-501, KB-502, KB-503). Previous numerical studies focused on the KB-502 injection well, where a double-lobe uplift pattern has been observed in the ground-deformation data. The observed uplift patterns at KB-501 and KB-503 have single-lobe patterns, but they can also indicate a deep fracture zone mechanical response to the injection. The current study improves the previous modeling approach by introducing an injection reservoir and a fracture zone, both responding to a Mohr-Coulomb failure criterion. In addition, we model a stress-dependent permeability and bulk modulus, according to a dual continuum model. Mechanical and hydraulic properties are determined through inverse modeling by matching the simulated spatial and temporal evolution of uplift to InSAR observations as well as by matching simulated and measured pressures. The numerical simulations are in agreement with both spatial and temporal observations. The estimated values for the parameterized mechanical and hydraulic properties are in good agreement with previous numerical results. In addition, the formal joint inversion of hydrogeological and geomechanical data provides measures of the estimation uncertainty.

Abstract. In July 2013, the city of St. Gallen conducted a deep geothermal project that aimed to exploit energy for district heating and generating power. A few days after an injection test and two acid stimulations that caused only minor seismicity, a gas kick forced the operators to inject drilling mud to combat the kick. Subsequently, multiple earthquakes were induced on a fault several hundred meters away from the well, including a ML 3.5 event that was felt throughout the nearby population centers. Given the occurrence of a gas kick and a felt seismic sequence with low total injected fluid volumes (~ 1200 m3), the St. Gallen deep geothermal project represents a particularly interesting case study of induced seismicity. Here, we first present a conceptual model based on seismic, borehole and seismological data suggesting a hydraulic connection between the well and the fault. The overpressurized gas, which is assumed to be initially sealed by the fault, may have been released due to the stimulations before entering the well via the hydraulic connection. We test this hypothesis with a numerical model calibrated against the borehole pressure of the injection test. We successfully reproduce the gas kick and the temporal and spatial characteristics of the main seismicity sequence that followed the well control operation. The results indicate that the gas may have destabilized the fault during and after the injection operations and could have enhanced the resulting seismicity. This study may have important implications for future deep hydrothermal projects conducted in similar geological conditions.