An official website of the United States government
Injection-induced seismicity and aseismic slip often involve the reactivation of long-dormant faults, which may have extremely low permeability prior to slip. In contrast, most previous models of fluid-driven aseismic slip have assumed linear pressure diffusion in a fault zone of constant permeability and porosity. Slip occurs within a frictional shear crack whose edge can either lag or lead pressure diffusion, depending on the dimensionless stress-injection parameter that quantifies the prestress and injection conditions. Here, we extend this foundational work by accounting for permeability enhancement and dilatancy, assumed to occur instantaneously upon the onset of slip. The fault zone ahead of the crack is assumed to be impermeable, so fluid flow and pressure diffusion are confined to the interior, slipped part of the crack. The confinement of flow increases the pressurization rate and reduction of fault strength, facilitating crack growth even for severely understressed faults. Suctions from dilatancy slow crack growth, preventing propagation beyond the hydraulic diffusion length. Our new two-dimensional and three-dimensional solutions can facilitate the interpretation of induced seismicity data sets. They are especially relevant for faults in initially low permeability formations, such as shale layers serving as caprock seals for geologic carbon storage, or for hydraulic stimulation of geothermal reservoirs. This article is part of the theme issue ‘Induced seismicity in coupled subsurface systems’.
more »
« less
Effect of Porosity and Permeability Evolution on Injection‐Induced Aseismic Slip
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
- Award ID(s):
- 1947448
- 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
-
-
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
-
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 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 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
