Collective properties of injection‐induced earthquake sequences: 1. Model description and directivity bias

Induced seismicity is of increasing concern for oil and gas, geothermal, and carbon sequestration operations, with several M > 5 events triggered in recent years. Modeling plays an important role in understanding the causes of this seismicity and in constraining seismic hazard. Here we study the collective properties of induced earthquake sequences and the physics underpinning them. In this first paper of a two-part series, we focus on the directivity ratio, which quantifies whether fault rupture is dominated by one (unilateral) or two (bilateral) propagating fronts. In a second paper, we focus on the spatiotemporal and magnitude-frequency distributions of induced seismicity. We develop a model that couples a fracture mechanics description of 1-D fault rupture with fractal stress heterogeneity and the evolving pore pressure distribution around an injection well that triggers earthquakes. The extent of fault rupture is calculated from the equations of motion for two tips of an expanding crack centered at the earthquake hypocenter. Under tectonic loading conditions, our model exhibits a preference for unilateral rupture and a normal distribution of hypocenter locations, two features that are consistent with seismological observations. On the other hand, catalogs of induced events when injection occurs directly onto a fault exhibit a bias toward ruptures that propagate toward the injection well. This bias is due to relatively favorable conditions for rupture that exist within the high-pressure plume. The strength of the directivity bias depends on a number of factors including the style of pressure buildup, the proximity of the fault to failure and event magnitude. For injection off a fault that triggers earthquakes, the modeled directivity bias is small and may be too weak for practical detection. For two hypothetical injection scenarios, we estimate the number of earthquake observations required to detect directivity bias.

[1]  J. Rice,et al.  Constitutive relations for fault slip and earthquake instabilities , 1983 .

[2]  J. Rice,et al.  A spectral method for three-dimensional elastodynamic fracture problems , 1995 .

[3]  A. Herrero,et al.  A kinematic self-similar rupture process for earthquakes , 1994, Bulletin of the Seismological Society of America.

[4]  H. Kanamori,et al.  Energy Partitioning During an Earthquake , 2006 .

[5]  Raul Madariaga,et al.  On the relation between seismic moment and stress drop in the presence of stress and strength heterogeneity , 1979 .

[6]  Jean-Paul Ampuero,et al.  Earthquake source characteristics from dynamic rupture with constrained stochastic fault stress , 2007 .

[7]  J. D. Eshelby.  The elastic field of a crack extending non-uniformly under general anti-plane loading , 1969 .

[8]  Bruce E. Shaw,et al.  Relation between stress heterogeneity and aftershock rate in the rate-and-state model , 2006 .

[9]  S. Ge,et al.  High-rate injection is associated with the increase in U.S. mid-continent seismicity , 2015, Science.

[10]  Roland N. Home,et al.  Numerical and analytical modeling of the mechanisms of induced seismicity during fluid injection , 2010 .

[11]  J. Dieterich,et al.  Modeling injection-induced seismicity with the physics-based earthquake simulator RSQSim , 2015 .

[12]  S. A. Miller,et al.  Modeling enhanced geothermal systems and the essential nature of large‐scale changes in permeability at the onset of slip , 2015 .

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

[14]  Alexandre Brouste,et al.  Characterization of Fault Roughness at Various Scales: Implications of Three-Dimensional High Resolution Topography Measurements , 2008, 0810.1109.

[15]  S. Shapiro,et al.  Microseismic rupture propagation imaging , 2015 .

[16]  J. McCloskey,et al.  Modelling fluid induced seismicity on a nearby active fault , 2013 .

[17]  Mark D. Zoback,et al.  Oklahoma’s recent earthquakes and saltwater disposal , 2015, Science Advances.

[18]  E. Dunham.  Dissipative interface waves and the transient response of a three dimensional sliding interface with Coulomb friction , 2005 .

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

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

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

[22]  G. Beroza,et al.  A spatial random field model to characterize complexity in earthquake slip , 2002 .

[23]  J. Rice,et al.  Elastodynamic analysis for slow tectonic loading with spontaneous rupture episodes on faults with rate‐ and state‐dependent friction , 2000 .

[24]  R. Console,et al.  Seismological study of Larderello — Travale geothermal area , 1985 .

[25]  Yehuda Ben-Zion,et al.  Asymmetric distribution of aftershocks on large faults in California , 2011 .

[26]  Olivier Lengliné,et al.  Rupture directivity of microearthquake sequences near Parkfield, California , 2011 .

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

[28]  D. J. Andrews,et al.  A stochastic fault model: 1. Static case , 1980 .

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

[30]  William W Rubey,et al.  ROLE OF FLUID PRESSURE IN MECHANICS OF OVERTHRUST FAULTING I. MECHANICS OF FLUID-FILLED POROUS SOLIDS AND ITS APPLICATION TO OVERTHRUST FAULTING , 1959 .

[31]  Y. Ben‐Zion.  Dynamic ruptures in recent models of earthquake faults , 2001 .

[32]  Jean-Paul Ampuero,et al.  Macroscopic Asymmetry of Dynamic Rupture on a Bimaterial Interface With Velocity- Weakening Friction , 2006 .

[33]  J. Rutqvist,et al.  Seismic rupture and ground accelerations induced by CO2 injection in the shallow crust , 2012 .

[34]  J. Dieterich,et al.  Implications of fault constitutive properties for earthquake prediction. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[35]  Yihe Huang,et al.  Collective properties of injection‐induced earthquake sequences: 2. Spatiotemporal evolution and magnitude frequency distributions , 2016 .

[36]  A Model for Tracking Fronts of Stress-Induced Permeability Enhancement , 2013, Transport in Porous Media.

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

[38]  Christian Pelties,et al.  On the initiation of sustained slip-weakening ruptures by localized stresses , 2015 .

[39]  Peter M. Shearer,et al.  Rupture directivity of small earthquakes at Parkfield , 2013 .

[40]  P. Mai,et al.  Hypocenter locations in finite-source rupture models , 2005 .

[41]  Albert Genter,et al.  Seismic response of the fractured and faulted granite of Soultz-sous-Forêts (France) to 5 km deep massive water injections , 2009 .

[42]  J. Rice.  Spatio‐temporal complexity of slip on a fault , 1993 .

[43]  G. Backus.  Interpreting the seismic glut moments of total degree two or less , 1977 .

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

[45]  Jonny Rutqvist,et al.  Impact of CO2 geological sequestration on the nucleation of earthquakes , 2011 .

[46]  M. Hubbert,et al.  ROLE OF FLUID PRESSURE IN MECHANICS OF OVERTHRUST FAULTING A REPLY , 1959 .

[47]  E. Dunham,et al.  Earthquake Ruptures with Strongly Rate-Weakening Friction and Off-Fault Plasticity, Part 1: Planar Faults , 2011 .

[48]  Raul Madariaga,et al.  On the Self-Healing Fracture Mode , 2003 .

[49]  Thomas H. Jordan,et al.  Predominance of Unilateral Rupture for a Global Catalog of Large Earthquakes , 2002 .

[50]  Shaun D. Fitzgerald,et al.  A note on induced stress changes in hydrocarbon and geothermal reservoirs , 1998 .

[51]  J. Weeks,et al.  Roughness and wear during brittle faulting , 1988 .

[52]  J. Boatwright,et al.  The Persistence of Directivity in Small Earthquakes , 2007 .

[53]  C. Marone.  LABORATORY-DERIVED FRICTION LAWS AND THEIR APPLICATION TO SEISMIC FAULTING , 1998 .

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

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

[56]  H. Freyer.  Preliminary15N studies on atmospheric nitrogenous trace gases , 1978 .

[57]  T. Heaton,et al.  Constraining fault constitutive behavior with slip and stress heterogeneity , 2008 .

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

[59]  Chien-Ching Ma,et al.  Dynamic fracture analysis of a finite crack subjected to an incident horizontally polarized shear wave , 1997 .

[60]  S. Day,et al.  Comparison of finite difference and boundary integral solutions to three‐dimensional spontaneous rupture , 2005 .

[61]  D. Wald,et al.  Spatial and temporal distribution of slip for the 1992 Landers, California, earthquake , 1994, Bulletin of the Seismological Society of America.

[62]  Daniel Moos,et al.  Numerical modeling of injection, stress and permeability enhancement during shear stimulation at the Desert Peak Enhanced Geothermal System , 2015 .

[63]  J. Molinari,et al.  Linear Elastic Fracture Mechanics Predicts the Propagation Distance of Frictional Slip , 2014, Tribology Letters.