Induced seismicity closed-form traffic light system for actuarial decision-making during deep fluid injections

The rise in the frequency of anthropogenic earthquakes due to deep fluid injections is posing serious economic, societal, and legal challenges to many geo-energy and waste-disposal projects. Existing tools to assess such problems are still inherently heuristic and mostly based on expert elicitation (so-called clinical judgment). We propose, as a complementary approach, an adaptive traffic light system (ATLS) that is function of a statistical model of induced seismicity. It offers an actuarial judgement of the risk, which is based on a mapping between earthquake magnitude and risk. Using data from six underground reservoir stimulation experiments, mostly from Enhanced Geothermal Systems, we illustrate how such a data-driven adaptive forecasting system could guarantee a risk-based safety target. The proposed model, which includes a linear relationship between seismicity rate and flow rate, as well as a normal diffusion process for post-injection, is first confirmed to be representative of the data. Being integrable, the model yields a closed-form ATLS solution that is both transparent and robust. Although simulations verify that the safety target is consistently ensured when the ATLS is applied, the model from which simulations are generated is validated on a limited dataset, hence still requiring further tests in additional fluid injection environments.

[1]  Thomas Göbel A comparison of seismicity rates and fluid-injection operations in Oklahoma and California: Implications for crustal stresses , 2015 .

[2]  Austin A. Holland,et al.  Earthquakes Triggered by Hydraulic Fracturing in South‐Central Oklahoma , 2013 .

[3]  A. Mignan Functional shape of the earthquake frequency‐magnitude distribution and completeness magnitude , 2012 .

[4]  A. Mignan Modeling aftershocks as a stretched exponential relaxation , 2015, 1510.01180.

[5]  Markus Häring,et al.  Characterisation of the Basel 1 enhanced geothermal system , 2008 .

[6]  L. Devroye Non-Uniform Random Variate Generation , 1986 .

[7]  Nicholas Deichmann,et al.  High-precision relocation and focal mechanism of the injection-induced seismicity at the Basel EGS , 2014 .

[8]  J. Bommer,et al.  Control of hazard due to seismicity induced by a hot fractured rock geothermal project , 2006 .

[9]  Stefan Wiemer,et al.  Bayesian Estimation of the Spatially Varying Completeness Magnitude of Earthquake Catalogs , 2010 .

[10]  Kristy F. Tiampo,et al.  Surface uplift and time-dependent seismic hazard due to fluid injection in eastern Texas , 2015, Science.

[11]  D. Giardini,et al.  Optimization of a large-scale microseismic monitoring network in northern Switzerland , 2011 .

[12]  Colin Camerer,et al.  Recent developments in modeling preferences: Uncertainty and ambiguity , 1992 .

[13]  Serge A. Shapiro,et al.  Decay rate of fluid-induced seismicity after termination of reservoir stimulations , 2010 .

[14]  G. Dresen,et al.  Effects of long‐term fluid injection on induced seismicity parameters and maximum magnitude in northwestern part of The Geysers geothermal field , 2015 .

[15]  J. N. Breunese,et al.  Induced seismicity of the Groningen gas field: History and recent developments , 2015 .

[16]  Peter A. W. Lewis,et al.  Statistical Analysis of Non-Stationary Series of Events in a Data Base System , 1976, IBM J. Res. Dev..

[17]  A. McGarr Maximum magnitude earthquakes induced by fluid injection , 2014 .

[18]  W. Grove,et al.  Clinical versus statistical prediction: the contribution of Paul E. Meehl. , 2005, Journal of clinical psychology.

[19]  Arnaud Mignan,et al.  Revisiting the 1894 Omori Aftershock Dataset with the Stretched Exponential Function , 2016 .

[20]  Jonny Rutqvist,et al.  Status of the TOUGH-FLAC simulator and recent applications related to coupled fluid flow and crustal deformations , 2011, Comput. Geosci..

[21]  R. Dawes,et al.  Heuristics and Biases: Clinical versus Actuarial Judgment , 2002 .

[22]  Domenico Giardini,et al.  Geothermal quake risks must be faced , 2009, Nature.

[23]  Baruch Fischhoff,et al.  The realities of risk-cost-benefit analysis , 2015, Science.

[24]  Serge A. Shapiro,et al.  Inter event times of fluid induced earthquakes suggest their Poisson nature , 2011 .

[25]  A. Mignan Reply to “Comment on ‘Revisiting the 1894 Omori Aftershock Dataset with the Stretched Exponential Function’ by A. Mignan” by S. Hainzl and A. Christophersen , 2016 .

[26]  W. Ellsworth Injection-Induced Earthquakes , 2013, Science.

[27]  Arnaud Mignan Static behaviour of induced seismicity , 2015 .

[28]  Yosihiko Ogata,et al.  Statistical Models for Earthquake Occurrences and Residual Analysis for Point Processes , 1988 .

[29]  Pierre Henry,et al.  Seismicity triggered by fluid injection–induced aseismic slip , 2015, Science.

[30]  Comment on “Revisiting the 1894 Omori Aftershock Dataset with the Stretched Exponential Function” by A. Mignan , 2016 .

[31]  G. Atkinson,et al.  “Did You Feel It?” Intensity Data: A Surprisingly Good Measure of Earthquake Ground Motion , 2007 .

[32]  Serge A. Shapiro,et al.  Seismotectonic state of reservoirs inferred from magnitude distributions of fluid-induced seismicity , 2012, Journal of Seismology.

[33]  Ted Urbancic,et al.  Towards an improved understanding of induced seismicity associated with hydraulic fracturing , 2016 .

[34]  Morgan T. Page,et al.  Induced earthquake magnitudes are as large as (statistically) expected , 2015 .

[35]  J. Vrijling,et al.  An overview of quantitative risk measures for loss of life and economic damage. , 2003, Journal of hazardous materials.

[36]  Joshua A. White,et al.  Assessing induced seismicity risk at CO2 storage projects: Recent progress and remaining challenges , 2016 .

[37]  Justin L. Rubinstein,et al.  Myths and Facts on Wastewater Injection, Hydraulic Fracturing, Enhanced Oil Recovery, and Induced Seismicity , 2015 .

[38]  Hans-Peter Harjes,et al.  Source parameters of injection-induced microearthquakes at 9 km depth at the KTB Deep Drilling site, Germany , 1998, Bulletin of the Seismological Society of America.

[39]  G. Shedler,et al.  Simulation of Nonhomogeneous Poisson Processes by Thinning , 1979 .

[40]  Lisa V. Block,et al.  Deep-Injection and Closely Monitored Induced Seismicity at Paradox Valley, Colorado , 2005 .

[41]  S. Maxwell,et al.  Geomechanical modeling of induced seismicity resulting from hydraulic fracturing , 2015 .

[42]  S. Ge,et al.  Sharp increase in central Oklahoma seismicity since 2008 induced by massive wastewater injection , 2014, Science.

[43]  H. Akaike A new look at the statistical model identification , 1974 .

[44]  Susan Petty,et al.  Improving Geothermal Project Economics with Multi-zone Stimulation: Results from the Newberry Volcano EGS Demonstration , 2013 .

[45]  Fabio Moia,et al.  A survey of the induced seismic responses to fluid injection in geothermal and CO2 reservoirs in Europe , 2012 .

[46]  K. Aki 17. Maximum Likelihood Estimate of b in the Formula logN=a-bM and its Confidence Limits , 1965 .

[47]  Yehuda Ben-Zion,et al.  Theoretical limits on detection and analysis of small earthquakes , 2016 .

[48]  Susan E. Hough,et al.  Shaking from Injection‐Induced Earthquakes in the Central and Eastern United States , 2014 .

[49]  Stefan Wiemer,et al.  Induced seismicity risk analysis of the 2006 Basel, Switzerland, Enhanced Geothermal System project: Influence of uncertainties on risk mitigation , 2015 .

[50]  Wei Huang,et al.  A Three-Dimensional Variational Data Assimilation System for MM5: Implementation and Initial Results , 2004 .

[51]  Álvaro González,et al.  The 2013 September–October seismic sequence offshore Spain: a case of seismicity triggered by gas injection? , 2014 .