An official website of the United States government
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
Earthquake Nucleation Characteristics Revealed by Seismicity Response to Seasonal Stress Variations Induced by Gas Production at Groningen
Abstract Deterministic earthquake prediction remains elusive, but time‐dependent probabilistic seismicity forecasting seems within reach thanks to the development of physics‐based models relating seismicity to stress changes. Difficulties include constraining the earthquake nucleation model and fault initial stress state. Here, we analyze induced earthquakes from the Groningen gas field, where production is strongly seasonal, and seismicity began 3 decades after production started. We use the seismicity response to stress variations to constrain the earthquake nucleation process and calibrate models for time‐dependent forecasting of induced earthquakes. Remarkable agreements of modeled and observed seismicity are obtained when we consider (a) the initial strength excess, (b) the finite duration of earthquake nucleation, and (c) the seasonal variations of gas production. We propose a novel metric to quantify the nucleation model's ability to capture the damped amplitude and the phase of the seismicity response to short‐timescale (seasonal) stress variations which allows further tightening the model's parameters.
more »
« less
- Award ID(s):
- 1822214
- PAR ID:
- 10536156
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 50
- Issue:
- 19
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Gas extraction from the Groningen gas field resulted in significant induced seismicity. We analyze the magnitude‐frequency distribution of these earthquakes in space, time and in view of stress changes calculated based on gas production and reservoir properties. Previous studies suggested variations related to reservoir geometry and stress. While we confirm the spatial variations, we do not detect a clear sensitivity of b‐value to Coulomb stress changes. However, we find that b‐value correlates positively with the rate of Coulomb stress changes. This correlation is statistically significant and robust to uncertainties related to stress change calculation. This study thus points to a possible influence of stress change rate on the probability of the magnitude of induced earthquakes.more » « less
-
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
-
Abstract In areas of induced seismicity, earthquakes can be triggered by stress changes due to fluid injection and static deformation from fault slip. Here we present a method to distinguish between injection‐driven and earthquake‐driven triggering of induced seismicity by combining a calibrated, fully coupled, poroelastic stress model of wastewater injection with interpretation of a machine learning algorithm trained on both earthquake catalog and modeled stress features. We investigate seismicity from Paradox Valley, Colorado as an ideal test case: a single, high‐pressure injector that has induced thousands of earthquakes since 1991. Using feature importance analysis, we find that injection‐driven earthquakes are approximately 225% of the total catalog but act as background events that can trigger subsequent aftershocks. Injection‐driven events also have distinct spatiotemporal clustering properties with a larger b‐value, closer proximity to the well, and earlier occurrence in the injection history. Generalization of our technique can help characterize triggering processes in other regions where induced seismicity occurs.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1029/2023gl105455
