An objective method for the assessment of fluid injection‐induced seismicity and application to tectonically active regions in central California

Changes in seismicity rates, whether of tectonic or of induced origin, can readily be identified in regions where background rates are low but are difficult to detect in seismically active regions. We present a novel method to identify likely induced seismicity in tectonically active regions based on short-range spatiotemporal correlations between changes in fluid injection and seismicity rates. The method searches through the entire parameter space of injection rate thresholds and determines the statistical significance of correlated changes in injection and seismicity rates. Applying our method to Kern County, central California, we find that most earthquakes within the region are tectonic; however, fluid injection contributes to seismicity in four different cases. Three of these are connected to earthquake sequences with events above M4. Each of these sequences followed an abrupt increase in monthly injection rates of at least 15,000 m^3. The probability that the seismicity sequences and the abrupt changes in injection rates in Kern County coincide by chance is only 4%. The identified earthquake sequences display low Gutenberg-Richter b values of ∼0.6–0.7 and at times systematic migration patterns characteristic for a diffusive process. Our results show that injection-induced pressure perturbations can influence seismic activity at distances of 10 km or more. Triggering of earthquakes at these large distances may be facilitated by complex local geology and faults in tectonically active regions. Our study provides the first comprehensive, statistically robust assessment of likely injection-induced seismicity within a large, tectonically active region.

[1]  Paul Segall,et al.  Induced stresses due to fluid extraction from axisymmetric reservoirs , 1992 .

[2]  Rachel E. Abercrombie,et al.  Triggering of the 1999 MW 7.1 Hector Mine earthquake by aftershocks of the 1992 MW 7.3 Landers earthquake , 2002 .

[3]  T. Kohl,et al.  Change of stress regime during geothermal reservoir stimulation , 2014 .

[4]  Domenico Giardini,et al.  Earthquakes Induced By the Stimulation of an Enhanced Geothermal System Below Basel (Switzerland) , 2009 .

[5]  David M. Evans.  The Denver Area Earthquakes and The Rocky Mountain Arsenal Disposal Well , 1966 .

[6]  W. Ellsworth,et al.  Incorporating induced seismicity in the 2014 United States National Seismic Hazard Model: results of the 2014 workshop and sensitivity studies , 2015 .

[7]  Cliff Frohlich,et al.  Two-year survey comparing earthquake activity and injection-well locations in the Barnett Shale, Texas , 2012, Proceedings of the National Academy of Sciences.

[8]  A. McGarr,et al.  Maximum magnitude earthquakes induced by fluid injection , 2014 .

[9]  R. Downie,et al.  Using Microseismic Source Parameters To Evaluate the Influence of Faults on Fracture Treatments: A Geophysical Approach to Interpretation , 2010 .

[10]  D. Schorlemmer,et al.  Off-fault damage and acoustic emission distributions during the evolution of structurally complex faults over series of stick-slip events , 2014 .

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

[12]  Y. Ogata.  Seismicity Analysis through Point-process Modeling: A Review , 1999 .

[13]  Hans-Peter Harjes,et al.  Injection-induced earthquakes and crustal stress at 9 km depth at the KTB deep drilling site, Germany , 1997 .

[14]  G. Dresen,et al.  Stress Tensor Changes Related to Fluid Injection at The Geysers Geothermal Field, California , 2013 .

[15]  P. Shearer,et al.  Waveform Relocated Earthquake Catalog for Southern California (1981 to June 2011) , 2012 .

[16]  Robert J. Skoumal,et al.  Optimizing multi-station earthquake template matching through re-examination of the Youngstown, Ohio, sequence , 2014 .

[17]  Geoffrey C. P. King,et al.  The accommodation of large strains in the upper lithosphere of the earth and other solids by self-similar fault systems: the geometrical origin of b-Value , 1983 .

[18]  Zoe K. Shipton,et al.  A review of recent developments concerning the structure, mechanics and fluid flow properties of fault zones , 2010 .

[19]  K. Felzer,et al.  Decay of aftershock density with distance indicates triggering by dynamic stress , 2006, Nature.

[20]  Cliff Frohlich,et al.  Did (Or Will) Fluid Injection Cause Earthquakes? - Criteria for a Rational Assessment , 1993 .

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

[22]  T. Jordan,et al.  Seismicity in Deep Gold Mines of South Africa: Implications for Tectonic Earthquakes , 2002 .

[23]  C. Scholz.  Earthquakes and friction laws , 1998, Nature.

[24]  E. Hauksson,et al.  Internal deformation of the southern Sierra Nevada microplate associated with foundering lower lithosphere, California , 2014 .

[25]  Shamita Das,et al.  Why large earthquakes do not nucleate at shallow depths , 1983, Nature.

[26]  Katie M. Keranen,et al.  Enhanced Remote Earthquake Triggering at Fluid-Injection Sites in the Midwestern United States , 2013, Science.

[27]  Mark E. J. Newman,et al.  Power-Law Distributions in Empirical Data , 2007, SIAM Rev..

[28]  J. Hardebeck.  Homogeneity of Small-Scale Earthquake Faulting, Stress, and Fault Strength , 2006 .

[29]  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 .

[30]  P. Segall,et al.  Poroelastic stressing and induced seismicity near the Lacq gas field, southwestern France , 1994 .

[31]  E. Brodsky,et al.  Anthropogenic Seismicity Rates and Operational Parameters at the Salton Sea Geothermal Field , 2013, Science.

[32]  T. Utsu.  Representation and Analysis of the Earthquake Size Distribution: A Historical Review and Some New Approaches , 1999 .

[33]  K. Felzer,et al.  Testing the stress shadow hypothesis , 2005 .

[34]  M. Wyss,et al.  Mapping the frequency-magnitude distribution in asperities: An improved technique to calculate recurrence times? , 1997 .

[35]  Qiang Wang,et al.  Earthquakes induced by water injection at ∼3 km depth within the Rongchang gas field, Chongqing, China , 2008 .

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

[37]  Danijel Schorlemmer,et al.  Acoustic emissions document stress changes over many seismic cycles in stick‐slip experiments , 2013 .

[38]  C. Frohlich,et al.  Two-year survey of earthquakes and injection/production wells in the Eagle Ford Shale, Texas, prior to the MW4.8 20 October 2011 earthquake , 2013 .

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

[40]  E. Brodsky,et al.  Connecting near-field and far-field earthquake triggering to dynamic strain , 2010 .

[41]  S. Wiemer,et al.  Influence of pore‐pressure on the event‐size distribution of induced earthquakes , 2012 .

[42]  T. Teng,et al.  Microearthquakes and water flooding in Los Angeles , 1973, Bulletin of the Seismological Society of America.

[43]  Paul Segall,et al.  Earthquakes triggered by fluid extraction , 1989 .

[44]  J. G. Williams,et al.  A Reservoir Analysis of the Denver Earthquakes: A Case of Induced Seismicity (Paper 80B 1095) , 1981 .

[45]  James P. Evans,et al.  Hydrogeologic Controls on Induced Seismicity in Crystalline Basement Rocks Due to Fluid Injection into Basal Reservoirs , 2013, Ground water.

[46]  D. Schorlemmer,et al.  Identifying fault heterogeneity through mapping spatial anomalies in acoustic emission statistics , 2011 .

[47]  D. Schorlemmer,et al.  Seismic event distributions and off-fault damage during frictional sliding of saw-cut surfaces with pre-defined roughness , 2014 .

[48]  J. Bredehoeft,et al.  An Experiment in Earthquake Control at Rangely, Colorado , 1976, Science.

[49]  Won-Young Kim,et al.  Induced seismicity associated with fluid injection into a deep well in Youngstown, Ohio , 2013 .

[50]  D. F. Sumy,et al.  Observations of static Coulomb stress triggering of the November 2011 M5.7 Oklahoma earthquake sequence , 2014 .

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

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

[53]  William L. Ellsworth,et al.  The 2001–Present Induced Earthquake Sequence in the Raton Basin of Northern New Mexico and Southern Colorado , 2014 .

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

[55]  S. Shapiro,et al.  Estimating the crust permeability from fluid-injection-induced seismic emission at the KTB site , 1997 .

[56]  Stefan Wiemer,et al.  Statistical analysis of the induced Basel 2006 earthquake sequence: introducing a probability-based monitoring approach for Enhanced Geothermal Systems , 2011 .

[57]  G. Dresen,et al.  Spatiotemporal changes, faulting regimes, and source parameters of induced seismicity: A case study from The Geysers geothermal field , 2014 .

[58]  Andrew J. Michael,et al.  Community online resource for statistical seismicity analysis , 2011 .

[59]  E. Brodsky,et al.  The Absence of Remotely Triggered Seismicity in Japan from 1997 to 2002 , 2003 .

[60]  H. Kanamori,et al.  A slow earthquake in the Santa Maria basin, California , 1992, Bulletin of the Seismological Society of America.

[61]  E. Hauksson,et al.  The tectonic crustal stress field and style of faulting along the Pacific North America Plate boundary in Southern California , 2012 .