We used the 3DEC code (46), a distinct element method (47), to simulate the interaction between fluid flow and fault slip evolution, including hydromechanical coupling and rate- and state-dependent friction. The model incorporates the full coupling between fluid pressure diffusion and effective stress- and strain-dependent permeability. The method has been previously used to understand the hydromechanical behavior of fractured rocks during fluid injection (48) and to study earthquake rupture and off-fault fracture response (49).

Our 3D model considers a fluid injection into a fault (dip angle of 70°) in a homogeneous elastic and impervious medium (Fig. 1C). The fault geometry reproduces the conditions of the in situ experiment (6). The remote normal (σn) and shear stress (τ) resolved on the fault plane are constant. During injection, the fluid pressure is increased into the fault in a point source using the loading path applied during the in situ experiment (blue curve in Fig. 1B). The initial value of stress (σno = 4.25 MPa, τo = 1.65 MPa) and fluid pressure (po = 0) into the fault represents the conditions of the in situ experiment (6). At time t = 0, water is injected into the fault so that the effective normal stress within the fault is defined by σn′ = σnp.

This model is based on the cubic law to describe the coupling between the fluid pressure and fault normal displacement (50)Q(t)=(bho+un)3wΔp12μfL(5)where Q(t) is the flow rate (in m3/s), w denotes the fault width (in m), bho is the initial hydraulic aperture of the fault (in m), un is the fault normal displacement (in m), μf is the fluid dynamic viscosity (in Pa·s), L is the contact length of the fault surface (in m), and Δp is the fluid pressure change (in Pa). The fluid pressure perturbation follows a diffusivity equationΔp.=TS2p(6)where T is the transmissivityT=(bh+un)312μf(7)and S is the storativityS=(bhKf+1K+4/3G)(8)where bh is the hydraulic aperture (bho + un) (in m), Kf is the fluid bulk modulus (in Pa), and K and G are the rock bulk and shear moduli (in Pa), respectively.

The magnitude of the pressure perturbation and the hydraulic properties affect the dynamics of the transient. We use an initial hydraulic aperture (bho = 9 × 10−6 m) consistent with an initial fault permeability of 7 × 10−12 m2 estimated in situ (4), Kf = 2 GPa, K = 20 GPa, G = 7.5 GPa, g = 9.81 m/s2, ρf = 1000 kg/m3, and μf = 1 × 10−3 Pa·s.

For the frictional behavior of the fault, the model is based on the rate-and-state friction law (Eq. 2) derived from laboratory experiments (1, 20, 36). The state evolution law is described by the aging law (20)dθdt=1vθdc(9)

Consistent with our laboratory measurements at low slip velocity (v < 10 μm/s), here, we assume a rate-weakening fault with μo = 0.6, a = 0.001, b = 0.003, and dc = 10 μm. The frictional resistance of the fault is then given byτ=μ(v,θ)(σnp)(10)

Note: The content above has been extracted from a research article, so it may not display correctly.

Please log in to submit your questions online.
Your question will be posted on the Bio-101 website. We will send your questions to the authors of this protocol and Bio-protocol community members who are experienced with this method. you will be informed using the email address associated with your Bio-protocol account.

We use cookies on this site to enhance your user experience. By using our website, you are agreeing to allow the storage of cookies on your computer.