[1] Elevated pore pressure can lead to reactivation of slip on pre-existing fractures and faults when the static Coulomb failure is reached locally. As the pressurized region spreads diffusively, slip can accumulate quasi-statically (paced by the pore fluid diffusion) or dynamically. In this work, we consider a prestressed fault with a locally peaked, diffusively spreading pore pressure field to study (1) conditions leading to the escalation of slip and nucleation of dynamic rupture and (2) rupture run-out distance before it is arrested. Nucleation appears in this model when the fault friction decreases from its peak value with slip, while arrest of dynamic propagation is imminent on aseismic faults (i.e., such that prestress t b is less than the residual fault strength t r at ambient conditions). When fluid overpressure is a small-to-moderate fraction of the ambient value of normal effective stress (and prestress is large enough for fault slip to be activated by overpressure), dynamic rupture always nucleates, and the nucleation length increases with decreasing prestress practically independently of the overpressure value. Transition from the ultimately unstable (t b > t r ) to the ultimately stable (t b < t r ) fault loading is marked by a strong increase of the nucleation length (∝1/(t b À t r ) 2 ) as t b approaches t r from above. For aseismic faults (t b < t r ), no dynamic rupture is nucleated at large fluid overpressures for all but the smallest values of prestress. The largest run-out distances of dynamic slip on aseismic faults correspond to overpressure/prestress just sufficient for slip activation. In such cases, the dynamically accumulated slip can lead to enhanced, dynamic fault weakening, resulting in a sustained dynamic rupture and generating a large earthquake. This is consistent with field observations when the largest injection-induced seismicity occurred after fluid injection ended.

No abstract

The focus of this paper is on constructing the solution for a semi-infinite hydraulic crack for arbitrary toughness, which accounts for the presence of a lag of a priori unknown length between the fluid front and the crack tip. First, we formulate the governing equations for a semi-infinite fluid-driven fracture propagating steadily in an impermeable linear elastic medium. Then, since the pressure in the lag zone is known, we suggest a new inversion of the integral equation from elasticity theory to express the opening in terms of the pressure. We then calculate explicitly the contribution to the opening from the loading in the lag zone, and reformulate the problem over the fluid-filled portion of the crack. The asymptotic forms of the solution near and away from the tip are then discussed. It is shown that the solution is not only consistent with the square root singularity of linear elastic fracture mechanics, but that its asymptotic behavior at infinity is actually given by the singular solution of a semi-infinite hydraulic fracture constructed on the assumption that the fluid flows to the tip of the fracture and that the solid has zero toughness. Further, the asymptotic solution for large dimensionless toughness is derived, including the explicit dependence of the solution on the toughness. The intermediate part of the solution (in the region where the solution evolves from the near tip to the far from the tip asymptote) of the problem in the general case is obtained numerically and relevant results are discussed, including the universal relation between the fluid lag and the toughness. [S0021-8936(00)02401-6]

Models for hydraulic fracturing–induced earthquakes in shales typically ascribe fault activation to elevated pore pressure or increased shear stress; however, these mechanisms are incompatible with experiments and rate-state frictional models, which predict stable sliding (aseismic slip) on faults that penetrate rocks with high clay or total organic carbon. Recent studies further indicate that the earthquakes tend to nucleate over relatively short injection time scales and sufficiently far from the injection zone that triggering by either poroelastic stress changes or pore pressure diffusion is unlikely. Here, we invoke an alternative model based on recent laboratory and in situ experiments, wherein distal, unstable regions of a fault are progressively loaded by aseismic slip on proximal, stable regions stimulated by hydraulic fracturing. This model predicts that dynamic rupture initiates when the creep front impinges on a fault region where rock composition favors dynamic and slip rate weakening behavior.

This paper is concerned with an analysis of the near-tip region of a fluid-driven fracture propagating in a permeable saturated rock. The analysis is carried out by considering the stationary problem of a semi-infinite fracture moving at constant speed V. Two basic dissipative processes are taken into account: fracturing of the rock and viscous flow in the fracture, and two fluid balance mechanisms are considered – leak-off and storage of the fracturing fluid in the fracture. It is shown that the solution is characterized by a multiscale singular behaviour at the tip, and that the nature of the dominant singularity depends both on the relative importance of the dissipative processes and on the scale of reference. This solution provides a framework to understand the interaction of representative physical processes near the fracture tip, as well as to track the changing nature of the dominant tip process(es) with the tip velocity and its impact on the global fracture response. Furthermore, it gives a universal scaling of the near-tip processes on the scale of the entire fracture and sets the foundation for developing efficient numerical algorithms relying on accurate modelling of the tip region.