scholarly journals On the effective stress law for rock-on-rock frictional sliding, and fault slip triggered by means of fluid injection

Author(s):  
Ernest Rutter ◽  
Abigail Hackston

Fluid injection into rocks is increasingly used for energy extraction and for fluid wastes disposal, and can trigger/induce small- to medium-scale seismicity. Fluctuations in pore fluid pressure may also be associated with natural seismicity. The energy release in anthropogenically induced seismicity is sensitive to amount and pressure of fluid injected, through the way that seismic moment release is related to slipped area, and is strongly affected by the hydraulic conductance of the faulted rock mass. Bearing in mind the scaling issues that apply, fluid injection-driven fault motion can be studied on laboratory-sized samples. Here, we investigate both stable and unstable induced fault slip on pre-cut planar surfaces in Darley Dale and Pennant sandstones, with or without granular gouge. They display contrasting permeabilities, differing by a factor of 10 5 , but mineralogies are broadly comparable. In permeable Darley Dale sandstone, fluid can access the fault plane through the rock matrix and the effective stress law is followed closely. Pore pressure change shifts the whole Mohr circle laterally. In tight Pennant sandstone, fluid only injects into the fault plane itself; stress state in the rock matrix is unaffected. Sudden access by overpressured fluid to the fault plane via hydrofracture causes seismogenic fault slips. This article is part of the themed issue ‘Faulting, friction and weakening: from slow to fast motion’.

2020 ◽  
Vol 11 (1) ◽  
Author(s):  
François X. Passelègue ◽  
Michelle Almakari ◽  
Pierre Dublanchet ◽  
Fabian Barras ◽  
Jérôme Fortin ◽  
...  

Abstract Modern geophysics highlights that the slip behaviour response of faults is variable in space and time and can result in slow or fast ruptures. However, the origin of this variation of the rupture velocity in nature as well as the physics behind it is still debated. Here, we first highlight how the different types of fault slip observed in nature appear to stem from the same physical mechanism. Second, we reproduce at the scale of the laboratory the complete spectrum of rupture velocities observed in nature. Our results show that the rupture velocity can range from a few millimetres to kilometres per second, depending on the available energy at the onset of slip, in agreement with theoretical predictions. This combined set of observations bring a new explanation of the dominance of slow rupture fronts in the shallow part of the crust or in areas suspected to present large fluid pressure.


2020 ◽  
Vol 110 (5) ◽  
pp. 2328-2349
Author(s):  
Kadek Hendrawan Palgunadi ◽  
Alice-Agnes Gabriel ◽  
Thomas Ulrich ◽  
José Ángel López-Comino ◽  
Paul Martin Mai

ABSTRACT The 15 November 2017 Mw 5.5 Pohang, South Korea, earthquake has been linked to hydraulic stimulation and fluid injections, making it the largest induced seismic event associated with an enhanced geothermal system. To understand its source dynamics and fault interactions, we conduct the first 3D high-resolution spontaneous dynamic rupture simulations of an induced earthquake. We account for topography, off-fault plastic deformation under depth-dependent bulk cohesion, rapid velocity weakening friction, and 1D subsurface structure. A guided fault reconstruction approach that clusters spatiotemporal aftershock locations (including their uncertainties) is used to identify a main and a secondary fault plane that intersect under a shallow angle of 15°. Based on simple Mohr–Coulomb failure analysis and 180 dynamic rupture experiments in which we vary local stress loading conditions, fluid pressure, and relative fault strength, we identify a preferred two-fault-plane scenario that well reproduces observations. We find that the regional far-field tectonic stress regime promotes pure strike-slip faulting, whereas local stress conditions constrained by borehole logging generate the observed thrust-faulting component. Our preferred model is characterized by overpressurized pore fluids, nonoptimally oriented but dynamically weak faults and a close-to-critical local stress state. In our model, earthquake rupture “jumps” to the secondary fault by dynamic triggering, generating a measurable non-double-couple component. Our simulations suggest that complex dynamic fault interaction may occur during fluid-injection-induced earthquakes and that local stress perturbations dominate over regional stress conditions. Therefore, our findings have important implications for seismic hazard in active georeservoir.


2021 ◽  
Author(s):  
Georg Dresen ◽  
Lei Wang ◽  
Grzegorz Kwiatek ◽  
Erik Rybacki ◽  
Audrey Bonnelye ◽  
...  

<p>Fluid injection causes fault slip that is partitioned in aseismic and seismic moment release. EGS stimulation campaigns have shown that in addition to total fluid volume injected also the rates of injection and fluid pressure increase affect seismic moment release. We investigate the effect of injection rate on slip characteristics, strain partitioning and energy budget in laboratory fluid injection experiments on reservoir sandstone samples in a triaxial deformation apparatus equipped with a 16-channel acoustic emission (AE) recording system. We injected fluid in sawcut samples containing a critically stressed fault at different pressurization rates. In general, fluid-induced fault deformation is dominantly aseismic. We find slow stick-slip events are induced at high fluid pressurization rate while steady fault creep occurs in response to low fluid pressurization rate. The released total seismic moment is found to be related to total injected volume, independent of fault slip behavior. Seismic moment release rate of AE is related to measured fault slip velocity. Total potential energy change and fracture energy release rate are defined by fault stiffness and largely independent of injection rate. Breakdown power density scales with slip rate and is significantly higher for fast injection and pressurization rates. The relation between moment release and injected volume is affected by fault slip behavior, characterized by a linear relation for slip at constant rate and fault creep while a cubic relation is evident for unstable and dynamic slip. Our experimental results allow separating a stable pressure-controlled injection phase with low rate of energy dissipation from a run-away phase, where breakdown power is high and cumulative moment release with injected volume is non-linear.</p>


2019 ◽  
Vol 5 (3) ◽  
pp. eaau4065 ◽  
Author(s):  
Frédéric Cappa ◽  
Marco Maria Scuderi ◽  
Cristiano Collettini ◽  
Yves Guglielmi ◽  
Jean-Philippe Avouac

Faults can slip seismically or aseismically depending on their hydromechanical properties, which can be measured in the laboratory. Here, we demonstrate that fault slip induced by fluid injection in a natural fault at the decametric scale is quantitatively consistent with fault slip and frictional properties measured in the laboratory. The increase in fluid pressure first induces accelerating aseismic creep and fault opening. As the fluid pressure increases further, friction becomes mainly rate strengthening, favoring aseismic slip. Our study reveals how coupling between fault slip and fluid flow promotes stable fault creep during fluid injection. Seismicity is most probably triggered indirectly by the fluid injection due to loading of nonpressurized fault patches by aseismic creep.


2021 ◽  
Vol 134 ◽  
pp. 104120
Author(s):  
Yijin Zeng ◽  
Qinghua Lei ◽  
Zineng Wang ◽  
Shidong Ding ◽  
Kui Liu ◽  
...  

2021 ◽  
Vol 9 (1) ◽  
Author(s):  
Daniel Bohnsack ◽  
Martin Potten ◽  
Simon Freitag ◽  
Florian Einsiedl ◽  
Kai Zosseder

AbstractIn geothermal reservoir systems, changes in pore pressure due to production (depletion), injection or temperature changes result in a displacement of the effective stresses acting on the rock matrix of the aquifer. To compensate for these intrinsic stress changes, the rock matrix is subjected to poroelastic deformation through changes in rock and pore volume. This in turn may induce changes in the effective pore network and thus in the hydraulic properties of the aquifer. Therefore, for the conception of precise reservoir models and for long-term simulations, stress sensitivity of porosity and permeability is required for parametrization. Stress sensitivity was measured in hydrostatic compression tests on 14 samples of rock cores stemming from two boreholes of the Upper Jurassic Malm aquifer of the Bavarian Molasse Basin. To account for the heterogeneity of this carbonate sequence, typical rock and facies types representing the productive zones within the thermal reservoir were used. Prior to hydrostatic investigations, the hydraulic (effective porosity, permeability) and geomechanical (rock strength, dynamic, and static moduli) parameters as well as the microstructure (pore and pore throat size) of each rock sample were studied for thorough sample characterization. Subsequently, the samples were tested in a triaxial test setup with effective stresses of up to 28 MPa (hydrostatic) to simulate in-situ stress conditions for depths up to 2000 m. It was shown that stress sensitivity of the porosity was comparably low, resulting in a relative reduction of 0.7–2.1% at maximum effective stress. In contrast, relative permeability losses were observed in the range of 17.3–56.7% compared to the initial permeability at low effective stresses. Stress sensitivity coefficients for porosity and permeability were derived for characterization of each sample and the different rock types. For the stress sensitivity of porosity, a negative correlation with rock strength and a positive correlation with initial porosity was observed. The stress sensitivity of permeability is probably controlled by more complex processes than that of porosity, where the latter is mainly controlled by the compressibility of the pore space. It may depend more on the compaction of precedented flow paths and the geometry of pores and pore throats controlling the connectivity within the rock matrix. In general, limestone samples showed a higher stress sensitivity than dolomitic limestone or dolostones, because dolomitization of the rock matrix may lead to an increasing stiffness of the rock. Furthermore, the stress sensitivity is related to the history of burial diagenesis, during which changes in the pore network (dissolution, precipitation, and replacement of minerals and cements) as well as compaction and microcrack formation may occur. This study, in addition to improving the quality of input parameters for hydraulic–mechanical modeling, shows that hydraulic properties in flow zones largely characterized by less stiff, porous limestones can deteriorate significantly with increasing effective stress.


2021 ◽  
Author(s):  
Hariharan Ramachandran ◽  
Andreia Plaza-Faverola ◽  
Hugh Daigle ◽  
Stefan Buenz

<p>Evidences of subsurface fluid flow-driven fractures (from seismic interpretation) are quite common at Vestnesa Ridge (around 79ºN in the Arctic Ocean), W-Svalbard margin. Ultimately, the fractured systems have led to the formation of pockmarks on the seafloor. At present day, the eastern segment of the ridge has active pockmarks with continuous methane seep observations in sonar data. The pockmarks in the western segment are considered inactive or to seep at a rate that is harder to identify. The ridge is at ~1200m water depth with the base of the gas hydrate stability zone (GHSZ) at ~200m below the seafloor. Considerable free gas zone is present below the hydrates. Besides the obvious concern of amount and rates of historic methane seeping into the ocean biosphere and its associated effects, significant gaps exist in the ability to model the processes of flow of methane through this faulted and fractured region. Our aim is to highlight the interactions between physical flow, geomechanics and geological control processes that govern the rates and timing of methane seepage.</p><p>For this purpose, we performed numerical fluid flow simulations. We integrate fundamental mass and component conservation equations with a phase equilibrium approach accounting for hydrate phase boundary effects to simulate the transport of gas from the base of the GHSZ through rock matrix and interconnected fractures until the seafloor. The relation between effective stress and fluid pressure is considered and fractures are activated once the effective stress exceeds the tensile limit. We use field data (seismic, oedometer tests on calypso cores, pore fluid pressure and temperature) to constrain the range of validity of various flow and geomechanical parameters in the simulation (such as vertical stress, porosity, permeability, saturations).</p><p>Preliminary results indicate fluid overpressure greater than 1.5 MPa is required to initiate fractures at the base of the gas hydrate stability zone for the investigated system. Focused fluid flow occurs through the narrow fracture networks and the gas reaches the seafloor within 1 day. The surrounding regions near the fracture network exhibit slower seepage towards the seafloor, but over a wider area. Advective flux through the less fractured surrounding regions, reaches the seafloor within 15 years and a diffusive flux reaches within 1200 years. These times are controlled by the permeability of the sediments and are retarded further due to considerable hydrate/carbonate formation during vertical migration. Next course of action includes constraining the methane availability at the base of the GHSZ and estimating its impact on seepage behavior.</p>


Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-23 ◽  
Author(s):  
Rob Westaway ◽  
Neil M. Burnside

The November 2017 MW 5.5 Pohang earthquake is one of the largest and most damaging seismic events to have occurred in the Korean peninsula over the last century. Its close proximity to an Enhanced Geothermal System (EGS) site, where hydraulic injection into granite had taken place over the previous two years, has raised the possibility that it was anthropogenic; if so, it was by far the largest earthquake caused by any EGS project worldwide. However, a potential argument that this earthquake was independent of anthropogenic activity considers the delay of two or three months before its occurrence, following the most recent injection into each of the wells. A better understanding of the physical and chemical processes that occur following fluid injection into granite is thus warranted. We show that hydrochemical changes occurring while surface water, injected into granite, reequilibrates chemically with its subsurface environment, can account for time delays for earthquake occurrence of such duration, provided the seismogenic fault was already critically stressed, or very close to the condition for slip. This candidate causal mechanism counters the potential argument that the time delay militates against an anthropogenic cause of the Pohang earthquake and can account for its relatively large magnitude as a consequence of a relatively small-volume injection. The resulting analysis places bounds on combinations of physical and chemical properties of rocks, injected volume, and potential postinjection time delays for significant anthropogenic seismicity during future EGS projects in granite.


1998 ◽  
Vol 41 (2) ◽  
Author(s):  
G. Selvaggi ◽  
F. D'Ajello Caracciolo

We analysed the one-year-long seismic swarm at the Alban Hills volcano which occurred during 1989-1990. We portray spatial distribution of seismic moment release, better delineating the activated volume during the swarm. The seismic structure is imaged as a 7-km long, 3-km wide, and 3-km thick volume, located between 2 and 5 km depth, and NW-SE striking. Fault plane solutions and scalar seismic moments for the largest earthquakes provide the description of the average strain rate tensor. The principal strain rate axes show a dominant extension in NE-SW direction, a SE-NW direction of compression and a negligible thickening rate. P and T axes direction of the smaller earthquakes suggests that the same mode of deformation is distributed all over the activated volume. These results are discussed in terms of seismic deforming processes active at the Alban Hills volcano, in the frame of magmatic inflation recently invoked to explain the rapid vertical uplift affecting part of the volcano. The observed average deformation is consistent with shear failures occurring on faults connecting stress-oriented dykes in response to an increasing fluid pressure.


2021 ◽  
Vol 40 (6) ◽  
pp. 413-417
Author(s):  
Chunfang Meng ◽  
Michael Fehler

As fluids are injected into a reservoir, the pore fluid pressure changes in space and time. These changes induce a mechanical response to the reservoir fractures, which in turn induces changes in stress and deformation to the surrounding rock. The changes in stress and associated deformation comprise the geomechanical response of the reservoir to the injection. This response can result in slip along faults and potentially the loss of fluid containment within a reservoir as a result of cap-rock failure. It is important to recognize that the slip along faults does not occur only due to the changes in pore pressure at the fault location; it can also be a response to poroelastic changes in stress located away from the region where pore pressure itself changes. Our goal here is to briefly describe some of the concepts of geomechanics and the coupled flow-geomechanical response of the reservoir to fluid injection. We will illustrate some of the concepts with modeling examples that help build our intuition for understanding and predicting possible responses of reservoirs to injection. It is essential to understand and apply these concepts to properly use geomechanical modeling to design geophysical acquisition geometries and to properly interpret the geophysical data acquired during fluid injection.


Sign in / Sign up

Export Citation Format

Share Document