Publications
Journal Articles
2027
- Three dimensional simulation of fluid-driven frictional and tensile ruptures on existing discontinuitiesBrice Lecampion , Sylvain Brisson , Antareep Sarma , Ankit Gupta , Alexis Sáez , and Regina FakhretdinovaComputer Methods in Applied Mechanics and Engineering, 2027
We present an implicit, fully-coupled hydro-mechanical solver for the three-dimensional simulation of fluid-driven rupture propagation along pre-existing discontinuities. The solver simultaneously handles frictional slip and tensile failure along arbitrary intersecting fractures and faults in a linearly elastic and impermeable rock matrix. Spatial discretization combines a displacement discontinuity boundary element method with a Galerkin finite element method for pore-fluid pressure diffusion. Frictional and tensile failure are governed by a poro-elastoplastic interface law incorporating slip-weakening friction, dilatancy, and tensile strength degradation. Block preconditioning of the coupled tangent system ensures robustness across a wide range of fracture behaviors, including friction and tensile hydraulic failure. Solver accuracy is verified – for the first time – against a comprehensive suite of semi-analytical rupture propagation solutions of increasing complexity: fluid-driven frictional ruptures, dilatant ruptures with permeability changes, and penny-shaped hydraulic fractures spanning the viscosity-to-toughness transition. Two multi-fracture examples further demonstrate the solver’s capabilities: injection into intersecting fractures, and a hydraulic fracture intersecting a strike-slip fault. These highlight the ability of the algorithm to capture frictional slip, dilatancy, permeability evolution, and tensile opening within a unified framework, making it well suited for fluid-driven rupture simulation in faulted and fractured rocks.
2025
- Maximum size and magnitude of injection-induced slow slip eventsAlexis Sáez , Francois Passelègue , and Brice LecampionScience Advances, 2025
Fluid injections can induce aseismic slip, resulting in stress changes that may propagate faster than pore pressure diffusion, potentially triggering seismicity at substantial distances from injection wells. Constraining the maximum extent of these aseismic ruptures is, thus, important for better delineating the influence zone of injections concerning their seismic hazard. Here, we derive a scaling relation based on rupture physics for the maximum size of aseismic ruptures, accounting for fluid injections with arbitrary flow rate histories. Moreover, on the basis of mounting evidence that the moment release during these operations is often predominantly aseismic, we derive a scaling relation for the maximum magnitude of aseismic slip events. Our theoretical predictions are consistent with observations over a broad spectrum of event sizes, from laboratory to real-world cases, indicating that fault zone storativity, background stress change, and injected fluid volume are key determinants of the maximum size and magnitude of injection-induced slow slip events. Scaling laws reveal how subsurface fluid injections control the maximum size and magnitude of induced aseismic slip events.
2024
- Fluid-driven slow slip and earthquake nucleation on a slip-weakening circular faultAlexis Sáez , and Brice LecampionJournal of the Mechanics and Physics of Solids, 2024
We investigate the propagation of fluid-driven fault slip on a slip-weakening frictional interface separating two identical half-spaces of a three-dimensional elastic solid. Our focus is on axisymmetric circular shear ruptures as they capture the most essential aspects of the dynamics of unbounded ruptures in three dimensions. In our model, fluid-driven aseismic slip occurs in two modes: as an interfacial rupture that is unconditionally stable, or as the quasi-static nucleation phase of an otherwise dynamic rupture. Unconditionally stable ruptures progress through four stages. Initially, ruptures are diffusively self-similar and the interface behaves as if it were governed by a constant friction coefficient equal to the static friction value. Slip then accelerates due to frictional weakening while the cohesive zone develops. Once the latter gets properly localized, a finite amount of fracture energy emerges along the interface. The rupture dynamics is then governed by an energy balance of the Griffith’s type. In this stage, fault slip transitions from a large-toughness to a small-toughness regime due to the diminishing effect of the fracture energy in the near-front energy balance. Ultimately, self-similarity is recovered and the fault behaves again as having a constant friction coefficient, but this time equal to the dynamic friction value. This condition is equivalent to a fault interface operating to leading order with zero fracture energy. When slow slip is the result of a frustrated dynamic instability, slip also initiates self-similarly at a constant peak friction coefficient. The maximum size that aseismic ruptures can reach before becoming unstable can be as small as a critical nucleation radius (shear modulus divided by the slip-weakening rate) and as large as infinity when faults are close to a well-defined limit that separates the two modes of aseismic sliding. In the former case, earthquake nucleation occurs unaffected by the dynamic friction coefficient. In contrast, the latter case exhibits fracture-mechanics behavior, characterized by a finite influx of elastic strain energy being supplied to and dissipated at the rupture front. We provide analytical and numerical solutions for the problem solved over its full dimensionless parameter space, including expressions for relevant length and time scales characterizing the transition between different stages and regimes. Due to its three-dimensional nature, the model enables quantitative comparisons with field observations as well as preliminary engineering design of hydraulic stimulation operations. Existing laboratory and in-situ experiments of fluid injection into simulated and natural faults are briefly discussed in light of our results.
2023
- Post-injection aseismic slip as a mechanism for the delayed triggering of seismicityAlexis Sáez , and Brice LecampionProceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2023
Injection-induced aseismic slip plays an important role in a broad range of human-made and natural systems, from the exploitation of geo-resources to the understanding of earthquakes. Recent studies have shed light on how aseismic slip propagates in response to continuous fluid injections. Yet much less is known about the response of faults after the injection of fluids has stopped. In this work, we investigate via a hydro-mechanical model the propagation and ultimate arrest of aseismic slip during the so-called post-injection stage. We show that after shut-in, fault slip propagates in pulse-like mode. The conditions that control the propagation as a pulse and notably when and where the ruptures arrest are fully established. In particular, critically stressed faults can host rupture pulses that propagate for several orders of magnitude the injection duration and reach up to nearly double the size of the ruptures at the moment of shut-in. We consequently argue that the persistent stressing of increasingly larger rock volumes caused by post-injection aseismic slip is a plausible mechanism for the triggering of post-injection seismicity—a critical issue in the geo-energy industry. Our physical model shows quantitative agreement with field observations of documented cases of post-injection induced seismicity.
2022
- Three-dimensional fluid-driven stable frictional rupturesAlexis Sáez , Brice Lecampion , Pathikrit Bhattacharya , and Robert C. ViescaJournal of the Mechanics and Physics of Solids, 2022
We investigate the quasi-static growth of a fluid-driven frictional shear crack that propagates in mixed mode (II+III) on a planar fault interface that separates two identical half-spaces of a three-dimensional solid. The fault interface is characterized by a shear strength equal to the product of a constant friction coefficient and the local effective normal stress. Fluid is injected into the fault interface and two different injection scenarios are considered: injection at constant volume rate and injection at constant pressure. We derive analytical solutions for circular ruptures which occur in the limit of a Poisson’s ratio ν=0 and solve numerically for the more general case in which the rupture shape is unknown (ν≠0). For an injection at constant volume rate, the fault slip growth is self-similar. The rupture radius (ν=0) expands as R(t)=λL(t), where L(t) is the nominal position of the fluid pressure front and λ is an amplification factor that is a known function of a unique dimensionless parameter T. The latter is defined as the ratio between the distance to failure under ambient conditions and the strength of the injection. Whenever λ>1, the rupture front outpaces the fluid pressure front. For ν≠0, the rupture shape is quasi-elliptical. The aspect ratio is upper and lower bounded by 1/(1−ν) and (3−ν)/(3−2ν), for the limiting cases of critically stressed faults (λ≫1, T≪1) and marginally pressurized faults (λ≪1, T≫1), respectively. Moreover, the evolution of the rupture area is independent of the Poisson’s ratio and grows simply as Ar(t)=4παλ2t, where α is the fault hydraulic diffusivity. For injection at constant pressure, the fault slip growth is not self-similar: the rupture front evolves at large times as ∝(αt)(1−T)/2 with T between 0 and 1. The frictional rupture moves at most diffusively (∝αt) when the fault is critically stressed, but in general propagates slower than the fluid pressure front. Yet in some conditions, the rupture front outpaces the fluid pressure front. The latter will eventually catch the former if injection is sustained for sufficient time. Our findings provide a basic understanding on how stable (aseismic) ruptures propagate in response to fluid injection in 3-D. Notably, since aseismic ruptures driven by injection at constant rate expand proportionally to the squared root of time, seismicity clouds that are commonly interpreted to be controlled by the direct effect of fluid pressure increase might be controlled by the stress transfer of a propagating aseismic rupture instead. We also demonstrate that the aseismic moment M0 scales to the injected fluid volume V as M0∝V3/2.
Articles Under Review
2026
- Earthquake sequence synchronization and alternation driven by quasistatic elastic interactions across faultsA. Rodriguez-Padilla , A. Sáez , and J.-P. Avouac2026Under review in Journal of Geophysical Research: Solid Earth
- Conditions for run-away ruptures during reservoir operations, implication for the magnitude-frequency distribution of induced earthquakes and application to GroningenQ. Shi , A. Sáez , L. Li , K. Im , S. Bourne , J. Elk , and J.-P. Avouac2026Under review in Proceedings of the National Academy of Sciences of the United States of America
- Recent magmatism sustains liquid water in a rapidly sealing Martian crustJ. Hernández-Montenegro , A. Sáez , J. Dong , and P. Asimow2026Under review in Nature
- Earthquakes induced by hydraulic fracturing, a semi-analytical model with application to the Permian BasinS. Marty , J.-P. Avouac , W. Curry , A. Sáez , and J. Silva2026Under review in Journal of Geophysical Research: Solid Earth
Proceedings/Conference Articles
2023
- Shearing and opening of a pre-existing discontinuity in response to fluid injectionB. Lecampion , A. Sáez , and A. GuptaIn U.S. Rock Mechanics/Geomechanics Symposium , Jun 2023
We investigate the fluid-driven growth ofa shear crack along a frictional discontinuity and its transition to hydraulic fracturing (sometimes referred to as hydraulic jacking) under plane-strain conditions. We focus on the case of a constant friction coefficient and account for the permeability changes associated with the fracture opening. By combining the scaling analysis and numerical simulations, we examine the evolution of both the shear and opening fronts as a function of the hydro-mechanical properties of the pre-existing discontinuity, in-situ stress state, and the fluid injection conditions. Further, we derive an approximate analytical solution for the relation between the positions of the slip and opening fronts at large times. We notably show that the ratio between the slip and opening fronts converges to a constant value at late times which only depends on the ratio between the shear stress and shear strength acting initially along the discontinuity. We compare this approximate solution against numerical simulations and demonstrate its usage to serve as a benchmark solution in the development of coupled hydro-mechanical numerical solvers for frictional fluid-driven fractures.Hydraulic stimulation of pre-existing fractures is used in the geothermal development in order to increase reservoir permeability and achieve economical flow rates - with mixed success (Jung, 2013; McClure and Horne, 2014). Although the primary idea is to shear dilate these pre-existing discontinuities via injection, in a number of field tests (Guglielmi et al., 2020), a large increase of permeability is only observed when fracture opening has been reached (sometimes denoted as hydraulic jacking). Shearing of pre-existing discontinuities can also occur during more traditional hydraulic fracturing operations in oil and gas reservoirs, either by direct fluid pressurization or via stress transfer from the main fractures. In this contribution, we investigate the fluid-driven growth of a shear crack along a frictional discontinuity and its transition to hydraulic fracturing. We focus on the case of a constant friction coefficient and account for permeability changes associated with fracture opening. We use a fully coupled implicit numerical scheme for the solution of this non-linear hydro-mechanical problem. Notably, the frictional interface behavior is modeled using an elastoplastic constitutive relation with a non-associated flow rule.
Thesis
2023
- Three-dimensional fluid-driven frictional ruptures: theory and applicationsAlexis SáezÉcole Polytechnique Fédérale de Lausanne (EPFL) , 2023
Fluid-driven frictional ruptures are important in a broad range of subsurface engineering technologies and natural earthquake-related phenomena. Some examples of subsurface operations where borehole fluid injections can induce frictional slip are deep geothermal energy, CO2 and hydrogen geological storage, and wastewater disposal from oil and gas production, among many others. On the other hand, natural phenomena such as seismic swarms, aftershock sequences, and slow earthquakes, are often associated with transients of pore-fluid pressure and fault slip. Motivated by the aforementioned applications and phenomena, this doctoral thesis aims to mechanistically understand how pre-existing geological structural discontinuities such as fractures and faults slide due to the pressurization of pore fluids, in a three-dimensional configuration that can be used for quantitative comparisons with field observations and preliminary engineering designs. We first develop a numerical solver for frictional slip and fault opening along pre-existing networks of discontinuities in 3D elastic media. We then examine a model with Coulomb’s friction to reproduce the initiation, propagation, and arrest of fluid-driven stable frictional ruptures. We show that a dimensionless number containing information about the initial stress state of the fault and the intensity of the injection governs the dynamics and shape of the ruptures in all its stages. Next, we extend the model to account for a friction coefficient that weakens with slip. This results in a broader range of fault slip behaviors, from unconditionally stable slip to dynamic instabilities. We quantify the conditions controlling both the propagation of stable slip and the nucleation of earthquakes. It is shown that the Coulomb’s friction model is both an early-time and late-time asymptotic solution of the more general slip-weakening model. After, we generalize some important fault rupture regimes to account for fairly arbitrary fluid injections. In particular, the connection between the expansion rate of the slipping surface and the history of injection volume rate is established. Using field observations, we then apply and explore the implications of our modeling results to injection-induced seismicity –a critical concern in the geo-energy industry. We show how the history of injection rate may control the migration patterns of micro-seismicity, and how these patterns may contain important information about in-situ conditions such as the fault stress state. Further, we elaborate on how post-injection pulses of stable slip can continue triggering seismicity due to stress-transfer effects, long after fluid injections stop. Finally, we propose a theoretical scaling relation for the maximum magnitude of injection-induced slow slip events, which sheds light on how aseismic motions release potential energy during injection operations. This doctoral research offers for the first time a quantitative and conceptual framework for understanding various subsurface and earthquake-science problems associated with frictional ruptures induced by fluid injections. While some implications of our findings are already explored in detail here, this work opens the possibility of understanding a much broader range of field observations. Moreover, our research paves the way for designing laboratory experiments, which is the next step to ground the theory we develop.