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
Effect of Porosity and Permeability Evolution on Injection‐Induced Aseismic Slip
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- Award ID(s):
- 1947448
- NSF-PAR ID:
- 10446201
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 126
- Issue:
- 7
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Fault-zone fluids control effective normal stress and fault strength. While most earthquake models assume a fixed pore fluid pressure distribution, geologists have documented fault valving behavior, that is, cyclic changes in pressure and unsteady fluid migration along faults. Here we quantify fault valving through 2-D antiplane shear simulations of earthquake sequences on a strike-slip fault with rate-and-state friction, upward Darcy flow along a permeable fault zone, and permeability evolution. Fluid overpressure develops during the interseismic period, when healing/sealing reduces fault permeability, and is released after earthquakes enhance permeability. Coupling between fluid flow, permeability and pressure evolution, and slip produces fluid-driven aseismic slip near the base of the seismogenic zone and earthquake swarms within the seismogenic zone, as ascending fluids pressurize and weaken the fault. This model might explain observations of late interseismic fault unlocking, slow slip and creep transients, swarm seismicity, and rapid pressure/stress transmission in induced seismicity sequences.
-
more » « less
Abstract Fluids influence fault zone strength and the occurrence of earthquakes, slow slip events, and aseismic slip. We introduce an earthquake sequence model with fault zone fluid transport, accounting for elastic, viscous, and plastic porosity evolution, with permeability having a power‐law dependence on porosity. Fluids, sourced at a constant rate below the seismogenic zone, ascend along the fault. While the modeling is done for a vertical strike‐slip fault with 2D antiplane shear deformation, the general behavior and processes are anticipated to apply also to subduction zones. The model produces large earthquakes in the seismogenic zone, whose recurrence interval is controlled in part by compaction‐driven pressurization and weakening. The model also produces a complex sequence of slow slip events (SSEs) beneath the seismogenic zone. The SSEs are initiated by compaction‐driven pressurization and weakening and stalled by dilatant suctions. Modeled SSE sequences include long‐term events lasting from a few months to years and very rapid short‐term events lasting for only a few days; slip is ∼1–10 cm. Despite ∼1–10 MPa pore pressure changes, porosity and permeability changes are small and hence fluid flux is relatively constant except in the immediate vicinity of slip fronts. This contrasts with alternative fault valving models that feature much larger changes in permeability from the evolution of pore connectivity. Our model demonstrates the important role that compaction and dilatancy have on fluid pressure and fault slip, with possible relevance to slow slip events in subduction zones and elsewhere.
-
more » « less
Abstract Inversions of InSAR ground deformation in the Delaware Basin have revealed an aseismic slip on semi‐optimally oriented normal faults located close to disposal wells. The slip, occurring over 3–5 years, extends approximately 1 km down‐dip, over 10 km along strike, and reaches 25 cm. We develop and calibrate 2D and pseudo‐3D coupled pore pressure diffusion and rate‐state models with velocity‐strengthening friction tailored to this unique height‐bounded fault geometry. Pressure diffusion is limited to a high‐permeability fault damage zone, and the net influx of fluid is adjusted to match the observed slip. A 1–2 MPa pressure increase initiates slip, with ∼5 MPa additional pressure increase required to produce ∼20 cm slip. Most slip occurs at approximately constant friction. Fault zone permeability must exceed ∼10−13 m2to match the along‐strike extent of slip. Models of the type developed here can be used to operationally manage injection‐induced aseismic slip.
-
more » « less
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}$$ $$\hbox {m}{^2}$$
