skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Investigating the role of thermal stresses on induced seismicity
We investigate the influence of thermal stresses on induced seismicity in the context of geothermal energy production at Brawley and Coso in California. We resort to thermo-hydro-mechanical modeling and used measurements of surface deformation to validate our simulations. We find that thermal stresses induced by fluid injections for geothermal energy production play an important role in either triggering or impeding seismicity. We speculate that thermal stresses may also play an important role in inducing fracturing and seismicity in other context of fluid injection for CO2 storage, disposal of wastewater or unconventional oil and gas production for example  more » « less
Award ID(s):
1822214
PAR ID:
10226188
Author(s) / Creator(s):
;
Date Published:
Journal Name:
SEG Technical Program Expanded Abstracts 2020
Page Range / eLocation ID:
1335 to 1338
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We investigate the relative importance of injection and production on the spatial‐temporal distribution of induced seismicity at the Raft River geothermal field. We use time‐series of InSAR measurements to document surface deformation and calibrate a hydro‐mechanical model to estimate effective stress changes imparted by injection and production. Seismicity, located predominantly in the basement, is induced primarily by poroelastic stresses from cold water reinjection into a shallower reservoir. The poroelastic effect of production from a deeper reservoir is minimal and inconsistent with observed seismicity, as is pore‐pressure‐diffusion in the basement and along reactivated faults. We estimate an initial strength excess of ∼20 kPa in the basement and sedimentary cover, but the seismicity rate in the sedimentary cover is four times lower, reflecting lower density of seed‐points for earthquake nucleation. Our modeling workflow could be used to assess the impact of fluid extraction or injection on seismicity and help design or guide operations. 
    more » « less
  2. Abstract Reservoir operations for gas extraction, fluid disposal, carbon dioxide storage, or geothermal energy production are capable of inducing seismicity. Modeling tools exist for seismicity forecasting using operational data, but the computational costs and uncertainty quantification (UQ) pose challenges. We address this issue in the context of seismicity induced by gas production from the Groningen gas field using an integrated modeling framework, which combines reservoir modeling, geomechanical modeling, and stress-based earthquake forecasting. The framework is computationally efficient thanks to a 2D finite-element reservoir model, which assumes vertical flow equilibrium, and the use of semianalytical solutions to calculate poroelastic stress changes and predict seismicity rate. The earthquake nucleation model is based on rate-and-state friction and allows for an initial strength excess so that the faults are not assumed initially critically stressed. We estimate uncertainties in the predicted number of earthquakes and magnitudes. To reduce the computational costs, we assume that the stress model is true, but our UQ algorithm is general enough that the uncertainties in reservoir and stress models could be incorporated. We explore how the selection of either a Poisson or a Gaussian likelihood influences the forecast. We also use a synthetic catalog to estimate the improved forecasting performance that would have resulted from a better seismicity detection threshold. Finally, we use tapered and nontapered Gutenberg–Richter distributions to evaluate the most probable maximum magnitude over time and account for uncertainties in its estimation. Although we did not formally account for uncertainties in the stress model, we tested several alternative stress models, and found negligible impact on the predicted temporal evolution of seismicity and forecast uncertainties. Our study shows that the proposed approach yields realistic estimates of the uncertainties of temporal seismicity and is applicable for operational forecasting or induced seismicity monitoring. It can also be used in probabilistic traffic light systems. 
    more » « less
  3. The development and operation of Enhanced Geothermal Systems (EGS) can induce earthquakes during fluid injections, posing significant hazards for both local populations and the continued operation of EGS plants. An improved understanding of the interplay between the physical mechanisms affecting geothermal induced seismicity is important for current and future EGS projects. One way to analyze seismicity is with the b-value, which represents the slope of earthquake magnitude frequency distributions. Previous EGS studies suggest a pore-pressure driven induced seismicity model where spatially – the b-value is high near the injection point and temporally – the b-value decreases during later stages of injection. We evaluate the pore-pressure driven model by comparing the spatiotemporal b-value evolutions of nine different EGS induced seismicity scenarios including Basel and FORGE, along with injections from Soultzsous-Forêts and Cooper Basin. In our spatiotemporal analyses, we find that a majority of examined sequences have a period where the distal b-value is significantly higher than the b-value near the injection point. These sequences indicate that the pore-pressure driven model is inadequate for describing spatiotemporal b-value evolution, and that additional physical mechanisms, like aseismic slip, could have significant effects on EGS induced seismicity. Further characterization of the b-value in EGS induced seismicity sequences provides an opportunity to better constrain the degrees to which different physical mechanisms influence seismicity. 
    more » « less
  4. Abstract Induced seismicity and surface deformation are common observable manifestations of the geomechanical effect of reservoir operations whether related to geothermal energy production, gas extraction or the storage of carbon dioxide, gas, air or hydrogen. Modelling tools to quantitatively predict surface deformation and seismicity based on operation data could thus help manage such reservoirs. To that effect, we present an integrated and modular modelling framework which combines reservoir modelling, geomechanical modelling and earthquake forecasting. To allow effective computational cost, we assume vertical flow equilibrium, semi-analytical Green's functions to calculate surface deformation and poroelastic stresses and a simple earthquake nucleation model based on Coulomb stress changes. We use the test case of the Groningen gas field in the Netherlands to validate the modelling framework and assess its usefulness for reservoir management. For this application, given the relative simplicity of this sandstone reservoir, we assume homogeneous porosity and permeability and single-phase flow. The model fits the measured pressure well, yielding a root mean square error (RMSE) of 0.95 MPa, and the seismicity observations as well. The pressure residuals show, however, a systematic increase with time that probably reflects groundwater ingression into the depleted reservoir. The interaction with groundwater could be accounted for by implementing a multiphase-flow vertical flow equilibrium (VFE) model. This is probably the major factor that limits the general applicability of the modelling framework. Nevertheless, he modelled subsidence and seismicity fit very well the historical observations in the case of the Groningen gas field. 
    more » « less
  5. Foreshocks are the most obvious signature of the earthquake nucleation stage and could, in principle, forewarn of an impending earthquake. However, foreshocks are only sometimes observed, and we have a limited understanding of the physics that controls their occurrence. In this work, we use high-resolution earthquake catalogs and estimates of source properties to understand the spatiotemporal evolution of a sequence of 11 foreshocks that occurred ~ 6.5 hours before the 2020 Mw 4.8 Mentone earthquake in west Texas.  Elevated pore-pressure and poroelastic stressing from subsurface fluid injection from oil-gas operations is often invoked to explain seismicity in west Texas and the surrounding region. However, here we show that static stresses induced from the initial ML 4.0 foreshock significantly perturbed the local shear stress along the fault and could have triggered the Mentone mainshock. The majority (9/11) of the earthquakes leading up to the Mentone mainshock nucleated in areas where the static shear stresses were increased from the initial ML 4.0 foreshock. The spatiotemporal properties of the 11 earthquakes that preceded the mainshock cannot easily be explained in the context of a preslip or cascade nucleation model. We show that at least 6/11 events are better classified as aftershocks of the initial ML 4.0.  Together, our results suggest that a combination of physical mechanisms contributed to the occurrence of the 11 earthquakes that preceded the mainshock, including static-stressing from earthquake-earthquake interactions, aseismic creep, and stress perturbations induced from fluid injection.  Our work highlights the role of earthquake-earthquake triggering in induced earthquake sequences, and suggests that such triggering could help sustain seismic activity following initial stressing perturbations from fluid injection. 
    more » « less