Statistical Method for Early Detection of Changes in Seismic Rate Associated with Wastewater Injections

It is well documented that wastewater disposal by underground injection can increase seismic activity. However, it is challenging to predict whether and when this effect occurs, as it can manifest itself even many years after the beginning of disposal operations. Continuous instrumental monitoring of the seismic activity is a first step enabling early detection of symptoms of change, updating of site and regional risk assessment, and prompt response. To avoid misdetection, however, seismic records have to be rationally analyzed, properly modeling the statistical features of the earthquake occurrence process. In this study, we develop a statistical approach to detect increments in seismic rate, accounting for model uncertainty (which is particularly acute when the monitoring period is short) and interdependence among events. The approach is composed of two steps: (1) stochastic earthquake declustering identifies mainshocks, and (2) the hypothesis of a constant rate of mainshocks is statistically tested. The method is applied to the analysis of the Oklahoma injection region, demonstrating that it is able to detect an increment in seismic rate before the change is large enough to produce severe consequences. We also investigate the statistical power of the method using synthetic data simulated for a wide range of scenarios. We believe that the method can be employed as a preliminary decision-support tool to identify areas with trends in seismicity that require further investigation and to guide adaptive management and operation of disposal wells.

[1]  Won-Young Kim,et al.  A Fluid-Injection-Triggered Earthquake Sequence in Ashtabula, Ohio: Implications for Seismogenesis in Stable Continental Regions , 2004 .

[2]  M. Christie,et al.  Masking of earthquake triggering behavior by a high background rate and implications for epidemic‐type aftershock sequence inversions , 2011 .

[3]  François Renard,et al.  Changes in effective stress during the 2003–2004 Ubaye seismic swarm, France , 2011 .

[4]  P. Talwani,et al.  Pore pressure diffusion and the mechanism of reservoir-induced seismicity , 1984 .

[5]  L. Knopoff,et al.  Is the sequence of earthquakes in Southern California, with aftershocks removed, Poissonian? , 1974, Bulletin of the Seismological Society of America.

[6]  R. Horne,et al.  Investigation of injection-induced seismicity using a coupled fluid flow and rate/state friction model , 2011 .

[7]  Stefan Wiemer,et al.  Geomechanical modeling of induced seismicity source parameters and implications for seismic hazard assessment , 2013 .

[8]  P. Reasenberg Second‐order moment of central California seismicity, 1969–1982 , 1985 .

[9]  Richard R. Hillis,et al.  Estimating fault stability and sustainable fluid pressures for underground storage of CO2 in porous rock , 2004 .

[10]  Yosihiko Ogata,et al.  Detecting fluid signals in seismicity data through statistical earthquake modeling , 2005 .

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

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

[13]  R. Saltus,et al.  Precambrian basement structure map of the continental United States - An interpretation of geologic and aeromagnetic data , 2008 .

[14]  D. Vere-Jones,et al.  Stochastic Declustering of Space-Time Earthquake Occurrences , 2002 .

[15]  R. Matsu’ura,et al.  A Point-process Analysis of the Matsushiro Earthquake Swarm Sequence: The Effect of Water on Earthquake Occurrence , 2005 .

[16]  Mark Naylor,et al.  Statistical Modeling of the 1997–1998 Colfiorito Earthquake Sequence: Locating a Stationary Solution within Parameter Uncertainty , 2014 .

[17]  G. Abers,et al.  Potentially induced earthquakes in Oklahoma, USA: Links between wastewater injection and the 2011 Mw 5.7 earthquake sequence , 2013 .

[18]  Dimitrios Karvounis,et al.  Towards a Real-Time Forecast of Induced Seismicity for Enhanced Geothermal Systems , 2014 .

[19]  Y. Ogata,et al.  The Centenary of the Omori Formula for a Decay Law of Aftershock Activity , 1995 .

[20]  William Walden Rubey,et al.  The Denver EarthquakeS. , 1968, Science.

[21]  Roger Fletcher,et al.  A Rapidly Convergent Descent Method for Minimization , 1963, Comput. J..

[22]  E. Lehmann Testing Statistical Hypotheses , 1960 .

[23]  Valentin Gischig,et al.  A stochastic model for induced seismicity based on non-linear pressure diffusion and irreversible permeability enhancement , 2013 .

[24]  S. Wiemer,et al.  Assessing the Quality of Earthquake Catalogues: Estimating the Magnitude of Completeness and Its Uncertainty , 2005 .

[25]  Stephen P. Horton,et al.  Disposal of Hydrofracking Waste Fluid by Injection into Subsurface Aquifers Triggers Earthquake Swarm in Central Arkansas with Potential for Damaging Earthquake , 2012 .

[26]  Warner Marzocchi,et al.  On the Increase of Background Seismicity Rate during the 1997-1998 Umbria-Marche, Central Italy, Sequence: Apparent Variation or Fluid-Driven Triggering? , 2010 .

[27]  Stephen H. Hickman,et al.  Introduction to Special Section: Mechanical Involvement of Fluids in Faulting , 1995 .

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

[29]  J. Campbell,et al.  Geologic Provinces of Oklahoma , 1996 .

[30]  A. Michael,et al.  Modeling Earthquake Rate Changes in Oklahoma and Arkansas: Possible Signatures of Induced Seismicity , 2013 .