skip to main content


Title: Fluid-induced aseismic fault slip outpaces pore-fluid migration
Earthquake swarms attributed to subsurface fluid injection are usually assumed to occur on faults destabilized by increased pore-fluid pressures. However, fluid injection could also activate aseismic slip, which might outpace pore-fluid migration and transmit earthquake-triggering stress changes beyond the fluid-pressurized region. We tested this theoretical prediction against data derived from fluid-injection experiments that activated and measured slow, aseismic slip on preexisting, shallow faults. We found that the pore pressure and slip history imply a fault whose strength is the product of a slip-weakening friction coefficient and the local effective normal stress. Using a coupled shear-rupture model, we derived constraints on the hydromechanical parameters of the actively deforming fault. The inferred aseismic rupture front propagates faster and to larger distances than the diffusion of pressurized pore fluid.  more » « less
Award ID(s):
1653382
NSF-PAR ID:
10096034
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Science
Volume:
364
Issue:
6439
ISSN:
0036-8075
Page Range / eLocation ID:
464 to 468
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Fluid injection stimulates seismicity far from active tectonic regions. However, the details of how fluids modify on‐fault stresses and initiate seismic events remain poorly understood. We conducted laboratory experiments using a biaxial loading apparatus with a 3 m saw‐cut granite fault and compared events induced at different levels of background shear stress. Water was injected at 10 mL/min and normal stress was constant at 4 MPa. In all experiments, aseismic slip initiated on the fault near the location of fluid injection and dynamic rupture eventually initiated from within the aseismic slipping patch. When the fault was near critically stressed, seismic slip initiated only seconds after MPa‐level injection pressures were reached and the dynamic rupture propagated beyond the fluid pressure perturbed region. At lower stress levels, dynamic rupture initiated hundreds of seconds later and was limited to regions where aseismic slip had significantly redistributed stress from within the pressurized region to neighboring locked patches. We found that the initiation of slow slip was broadly consistent with a Coulomb failure stress, but that initiation of dynamic rupture required additional criteria to be met. Even high background stress levels required aseismic slip to modify on‐fault stress to meet initiation criteria. We also observed slow slip events prior to dynamic rupture. Overall, our experiments suggest that initial fault stress, relative to fault strength, is a critical factor in determining whether a fluid‐induced rupture will “runaway” or whether a fluid‐induced rupture will remain localized to the fluid pressurized region.

     
    more » « less
  2. Abstract

    There is a growing recognition that subsurface fluid injection can produce not only earthquakes, but also aseismic slip on faults. A major challenge in understanding interactions between injection-related aseismic and seismic slip on faults is identifying aseismic slip on the field scale, given that most monitored fields are only equipped with seismic arrays. We present a modeling workflow for evaluating the possibility of aseismic slip, given observational constraints on the spatial-temporal distribution of microseismicity, injection rate, and wellhead pressure. Our numerical model simultaneously simulates discrete off-fault microseismic events and aseismic slip on a main fault during fluid injection. We apply the workflow to the 2012 Enhanced Geothermal System injection episode at Cooper Basin, Australia, which aimed to stimulate a water-saturated granitic reservoir containing a highly permeable ($$k = 10^{-13} - 10^{-12}$$k=10-13-10-12$$\hbox {m}{^2}$$m2) fault zone. We find that aseismic slip likely contributed to half of the total moment release. In addition, fault weakening from pore pressure changes, not elastic stress transfer from aseismic slip, induces the majority of observed microseismic events, given the inferred stress state. We derive a theoretical model to better estimate the time-dependent spatial extent of seismicity triggered by increases in pore pressure. To our knowledge, this is the first time injection-induced aseismic slip in a granitic reservoir has been inferred, suggesting that aseismic slip could be widespread across a range of lithologies.

     
    more » « less
  3. Abstract

    To better quantify how injection, prior seismicity, and fault properties control rupture growth and propagation of induced earthquakes, we perform a finite‐fault slip inversion on aMw4.0 earthquake that occurred in April 2015, the largest earthquake in an induced sequence near Guthrie, Oklahoma. The slip inversion reveals a rupture with slip patches that are anti‐correlated to the locations of prior seismicity. The prior seismicity driven by low pore pressure changes and static stress changes occurred on weaker portions of the fault, while theMw4.0 earthquake likely ruptured relatively stronger portions of the fault. To resolve if pore pressure changes or the initial underlying stress distribution and fault strength controlled the final slip distribution of the GuthrieMw4.0 earthquake, we compare strike‐slip events of similar magnitude from tectonically active regions and previously inactive regions. Earthquakes on reactivated faults exhibit different slip distributions than active regions, they have more prominent and well separated slip patches, a behavior often associated with faults of lower fault maturity. Pore pressure shows little effect on the distributions. These observations suggest that the initial underlying stress distribution and fault strength of reactivated faults in low deformation regions is the primary controlling factor of the slip distribution with pore pressure perturbations and earthquake interactions being secondary. Therefore, GuthrieMw4.0 earthquakes slip distribution was enhanced by pore‐pressure perturbations and earthquake interactions by creating an optimal stress state for its failure, but the slip distribution itself is controlled by its fault's initial stress and strength state.

     
    more » « less
  4. Abstract

    Microseismicity associated with fluid pressurization in the subsurface occurs during fluid injection but can also be triggered after injection shut‐in. Understanding the extent and duration of the post‐injection microseismicity is critical to limit the risk of fluid‐induced seismicity and insure the safe utilization of the subsurface. Using theoretical and numerical techniques, we investigated how aseismic slip on a fault plane evolves and stops after a fluid pressurization event. We found that the locking mechanisms controlling the arrest of aseismic slip highly depend on the initial fault stress criticality and the pressurization duration. The absolute arrest time of fault aseismic slip after injection shut‐in is proportional to the pressurization duration and increases significantly with the initial fault stress criticality. Given that microseismicity can be triggered by aseismic slip, these results provide insights into the mechanics controlling the arrest of microseismicity after fluid pressurization as a milestone toward induced seismicity mitigation strategies.

     
    more » « less
  5. Abstract

    It is widely recognized that fluid injection can trigger aseismic fault slip. However, the processes by which the fluid‐rock interactions facilitate or inhibit slip are poorly understood and some are oversimplified in most models of injection‐induced slip. In this study, we perform a 2D anti‐plane shear investigation of aseismic slip that occurs in response to fluid injection into a permeable fault governed by rate‐and‐state friction. We account for porosity and permeability changes that accompany slip, including dilatancy, and quantify how these processes affect pore pressure diffusion, which couples to aseismic slip. Fault response to injection has two phases. In the first phase, slip is negligible and pore pressure closely follows the standard linear diffusion model. Pressurization eventually triggers aseismic slip close to the injection site. In the second phase, aseismic slip front expands outward and dilatancy causes pore pressure to depart from the linear diffusion model. We quantify how prestress, injection rate, permeability and other fluid transport properties affect the slip front migration rate, finding rates ranging from 10 to 1,000 m/day for typical parameters. The migration rate is strongly influenced by the fault's closeness to failure and injection rate. The total slip on the fault, on the other hand, is primarily determined by the injected volume, with minimal sensitivity to injection rate. Additionally, we show that when dilatancy is neglected, slip front migration rate and total slip can be several times higher. Our modeling demonstrates that porosity and permeability evolution, especially dilatancy, fundamentally alters how faults respond to fluid injection.

     
    more » « less