Earthquakes Triggered by Hydraulic Fracturing in South‐Central Oklahoma

In January 2011, a sequence of earthquakes occurred in close proximity to a well, which was being hydraulically fractured in south‐central Oklahoma. The hydraulic fracturing of the Picket Unit B Well 4–18 occurred from 16 January 2011 18:43 through 22 January 16:54 UTC. This vertical well penetrated into the mature Eola‐Robberson oil field. Earthquakes were identified by cross correlating template waveforms from manually identified earthquakes and cross correlating these templates through the entire operation period of the Earthscope USArray Transportable Array (TA) station X34A. This produced a series of 116 earthquakes, which occurred from 17 January 2011 19:06 through 23 January 3:13 UTC with no other similar earthquakes identified at other times prior to or post‐hydraulic fracturing. The identified earthquakes range in local magnitude ( M L) from 0.6 to 2.9, with 16 earthquakes M L 2 or greater and a b ‐value of 0.98. There is a strong temporal correlation between hydraulic fracturing and earthquakes. This correlation is strengthened because hydraulic fracturing operations ceased for ∼2 days due to bad weather, and earthquakes can be observed to cease during this period and resume after hydraulic fracturing had resumed. Earthquakes were relocated using cross‐correlated phase arrivals and bootstrap iterations of hypoDD. Locations were well constrained for 86 earthquakes. These earthquake locations clearly delineate a fault which strikes ∼166°, subparallel to the mapped minor fault sets in the area, and dips steeply to the west. The earthquakes appear to have occurred at shallow depths from ∼2 to 3 km and within ∼2.5  km horizontally of the well. The first earthquake occurred ∼24 hrs after hydraulic fracturing began at the well. This delay is consistent with the diffusion of pore pressure in the subsurface over a distance of ∼2  km. Online Material: Results from bootstrap hypoDD relocations using cross‐correlation phase arrivals.

[1]  J. McCloskey,et al.  What controls an earthquake's size? Results from a heterogeneous cellular automaton , 1998 .

[2]  Julian J. Bommer,et al.  Induced seismicity associated with Enhanced Geothermal Systems , 2007 .

[3]  J. G. McCaskill.  Multiple Stratigraphic Indicators of Major Strike-slip Movement along the Eola Fault, Subsurface Arbuckle Mountains, Oklahoma , 1998 .

[4]  Ted Urbancic,et al.  Microseismic Imaging of Hydraulic Fracture Complexity in the Barnett Shale , 2002 .

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

[6]  G. R. Keller,et al.  THE ROLE OF RIFTING IN THE TECTONIC DEVELOPMENT OF THE MIDCONTINENT, U.S.A. , 1983 .

[7]  L. Ceranna,et al.  Probing the Crust to 9-km Depth: Fluid-Injection Experiments and Induced Seismicity at the KTB Superdeep Drilling Hole, Germany , 2002 .

[8]  Felix Waldhauser,et al.  Waveform Cross-Correlation-Based Differential Travel-Time Measurements at the Northern California Seismic Network , 2005 .

[9]  J. Tanner.  Wrench Fault Movements Along Washita Valley Fault, Arbuckle Mountain Area, Oklahoma: GEOLOGICAL NOTES , 1967 .

[10]  Alexander Y. Rozhko,et al.  Role of seepage forces on seismicity triggering , 2010 .

[11]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[12]  R. Wheeler,et al.  Known and suggested quaternary faulting in the midcontinent United States , 2001 .

[13]  C. Langston,et al.  Empirical Distance Attenuation and the Local-Magnitude Scale for the Central United States , 2007 .

[14]  D. L. Anderson,et al.  Theoretical Basis of Some Empirical Relations in Seismology by Hiroo Kanamori And , 1975 .

[15]  Lion Krischer,et al.  ObsPy: A Python Toolbox for Seismology , 2010 .

[16]  Marco Tedesco,et al.  Snowmelt detection over the Greenland ice sheet from SSM/I brightness temperature daily variations , 2007 .

[17]  B. Gutenberg,et al.  Frequency of Earthquakes in California , 1944, Nature.

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

[19]  N. R. Warpinski,et al.  Microseismic Monitoring: Inside and Out , 2009 .

[20]  Felix Waldhauser,et al.  hypoDD -- A Program to Compute Double-Difference Hypocenter Locations , 2001 .

[21]  S. Maxwell,et al.  Fault Activation During Hydraulic Fracturing , 2010 .

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

[23]  C. Dinske,et al.  Probability of a given‐magnitude earthquake induced by a fluid injection , 2007 .

[24]  B. H. Harlton.  Tectonic Framework of Eola and Southeast Hoover Oil Fields and West Timbered Hills Area, Garvin and Murray Counties, Oklahoma , 1964 .

[25]  J. Granath.  Structural evolution of the Ardmore Basin, Oklahoma: Progressive deformation in the foreland of the Ouachita Collision , 1989 .

[26]  Serge A. Shapiro,et al.  Microseismic monitoring of borehole fluid injections: Data modeling and inversion for hydraulic properties of rocks , 2003 .

[27]  Robert L. Wesson,et al.  Earthquake hazard associated with deep well injection: a report to the U.S. Environmental Protection Agency , 1990 .

[28]  P. Talwani,et al.  Seismogenic permeability, ks , 2007 .

[29]  S. Maxwell,et al.  A comparison of passive seismic monitoring of fracture stimulation from water and CO2 injection , 2010 .

[30]  D. Wyborn,et al.  Induced Seismicity during the Stimulation of a Geothermal hfr Reservoir in the Cooper Basin, Australia , 2006 .

[31]  Travis E. Oliphant,et al.  Python for Scientific Computing , 2007, Computing in Science & Engineering.

[32]  Jens Havskov,et al.  A Computer Program for Locating Earthquakes Both Locally and Globally , 1995 .

[33]  F. Ringdal,et al.  The detection of low magnitude seismic events using array-based waveform correlation , 2006 .

[34]  B. Lienert,et al.  HYPOCENTER: An earthquake location method using centered, scaled, and adaptively damped least squares , 1986 .

[35]  Pascal Audigane,et al.  Large-scale in situ permeability tensor of rocks from induced microseismicity , 1999 .

[36]  D. Nychka,et al.  Spatial patterns of probabilistic temperature change projections from a multivariate Bayesian analysis , 2007 .

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

[38]  F. Waldhauser,et al.  A Double-Difference Earthquake Location Algorithm: Method and Application to the Northern Hayward Fault, California , 2000 .

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

[40]  Richard M. Allen,et al.  The deterministic nature of earthquake rupture , 2005, Nature.

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