Description: The influence of elevated temperature on injection-induced fault slip is poorly constrained. In this study, at steady-state elevated temperatures, triaxial shear-flow experiments on a sawcut fault in granite were conducted to simulate injection-induced slip of a critically stressed fault. Our results suggest that an elevated temperature favors a more uniform fluid pressure distribution over the fault surface mainly by reducing water viscosity. At temperatures above ambient, a larger perturbation force from the injected fluid is required to reactivate the fault primarily because of the enhanced thermally activated fault healing processes, resulting in a faster fault slip rate upon failure. This study may partially explain the causal link between higher reservoir temperature and higher maximum magnitude of injection-induced earthquakes in geothermal systems, and the observation that larger magnitude seismic events concentrate near the deeper part of the reservoir, where temperature is higher.

Description: Earthquakes associated with fluid injection in various geo-energy settings, such as shale gas and deep geothermal energy, have shelved many projects with great potential. However, the injection-rate dependence of earthquake nucleation length, i.e., the slowly slipping (creeping) fault length in preparation for a subsequent earthquake (Kaneko & Lapusta, 2008), remains elusive. In this study, we take a step towards this issue by performing fluid injection experiments on low-permeability granite samples containing a critically stressed sawcut fault at different local injection rates (0.2 mL/min and 0.8 mL/min) and confining pressures (31 MPa and 61 MPa) (c. f., Ji & Wu, 2017; Wang et al., 2020). An array of local strain gauges and acoustic emission (AE) hypocenter locations were used to monitor the precursory slip of critically stressed faults before injection-induced stick-slip failure (c. f., Passelègue et al., 2020; Wang et al., 2020). The nucleation length was determined for each injection-induced stick-slip event, and its dependence on effective normal stress and injection rate was explored. Herein, we compile the processed data obtained from the experiments in four Excel worksheets. The full description of the methods is provided in Ji et al. (2022).

Description: The utilisation of geothermal energy in the Netherlands is primarily focused on deep sedimentary aquifers, which are often intersected by major faults. Geothermal operations (i.e. fluid production and injection) may alter the effective stress state along these faults and trigger induced seismic events. Pore pressure perturbations have been generally considered the main driver of injection-induced seismicity. However, thermal stresses caused by temperature gradients between the re-injected cold fluid and the reservoir rock may also contribute to the triggering of earthquakes in geothermal reservoirs. While existing geothermal power plants operating in sandstone reservoirs did not produce any major induced seismicity, it is a matter of debate whether a reduction in the temperature of the re-injected fluid could increase the seismic hazard potential. In this study, we applied modified Gutenberg–Richter statistics based on frictional Coulomb stress variations implemented in a coupled thermo-hydro-mechanical model to estimate the seismic hazard caused by the operation of a geothermal doublet. We conducted a systematic parametric study to assess and rank the impact of different intrinsic (geological) and extrinsic (operational) parameters on the induced seismic hazard potential. We identified a competing mechanism between induced variations in pore pressure and thermal stress within the reservoir in controlling induced seismicity. We found that stress changes induced by pore pressure variations are the main cause of seismic hazard, although thermally induced stresses also contribute significantly. The results indicate that by optimising the operational parameters it is possible to increase production efficiency while maintaining a long-term control over the fluid injection-induced seismicity.

Description: This study investigates numerically several hydraulic fracturing experiments that were performed on intact cubic Pocheon granite samples applying different injection protocols. The goal of the laboratory experiments is to test the concept of cyclic soft stimulation which aims to increase permeability sustainably among others. The Irazu 2D numerical code is used to simulate explicitly coupled hydraulic diffusion and fracturing processes under bi-axial stress conditions. Using the hybrid finite-discrete element modelling approach, we test two injection schemes, constant-rate continuous injection and cyclic progressive injection on homogeneous and heterogeneous samples. Our study focuses on the connection between the geometry of hydraulic fractures, fracturing mechanisms and the permeability increase after injection. The models capture several characteristics of the hydraulic fracturing tests using a time-scaling approach. The numerical simulation results show good agreement with the laboratory experiments in terms of pressure evolution characteristics and fracture pattern. Based on the simulation results, the constant-rate continuous and cyclic progressive injection schemes applied to heterogeneous rock sample with pre-existing fractures show the highest hydraulic aperture increase, and thus permeability enhancement.

Description: The effect of normal stress variations on fault frictional strength has been extensively characterized in laboratory experiments and modelling studies based on a rate-and-state-dependent fault friction formalism. However, the role of pore pressure changes during injection-induced fault reactivation and associated frictional phenomena is still not well understood. We apply rate-and-state friction (RSF) theory in finite element models to investigate the effect of fluid pressurization rate on fault (re)activation and on the resulting frictional slip characteristics at the laboratory scale. We consider a stepwise injection scenario where each fluid injection cycle consists of a fluid pressurization phase followed by a constant fluid pressure phase. We first calibrate our model formulation to recently published laboratory results of injection-driven shear slip experiments. In a second stage, we perform a parametric study by varying fluid pressurization rates to cover a higher dimensional parameter space. We demonstrate that, for high permeability laboratory samples, the energy release rate associated with fault reactivation can be effectively controlled by a stepwise fluid injection scheme, i.e. by the applied fluid pressurization rate and the duration of the constant pressure phase between each successive fluid pressurization phase. We observe a gradual transition from fault creep to slow stick–slip as the fluid pressurization rate increases. Furthermore, computed peak velocities for an extended range of fluid pressurization rate scenarios (0.5 MPa/min to 10 MPa/min) indicate a non-linear (power-law) relationship between the imposed fluid pressurization rate and the peak slip velocities, and consequently with the energy release rate, for scenarios with a fluid pressurization rate higher than a critical value of 4 MPa/min. We also observe that higher pressurization rates cause a delay in the stress release by the fault. We therefore argue that by adopting a stepwise fluid injection scheme with lower fluid pressurization rates may provide the operator with a better control over potential induced seismicity. The implications for field-scale applications that we can derive from our study are limited by the high matrix and fault permeability of the selected sample and the direct hydraulic connection between the injection well and the fault, which may not necessarily represent the conditions typical for fracture dominated deep geothermal reservoirs. Nevertheless, our results can serve as a basis for further laboratory experiments and field-scale modelling studies focused on better understanding the impact of stepwise injection protocols on fluid injection-induced seismicity.

Description: We conducted fluid injection experiments on cylindrical low-permeability granite samples with a critically stressed sawcut fault at local injection rates of 0.2 and 0.8 mL/min and confining pressures of 31 and 61 MPa. A local array of six strain gauges attached close to the faults allows us to estimate the nucleation length of each injection-induced dynamic slip event (i.e., laboratory earthquake). We find nucleation lengths decrease from approximately 90% to 〈15% of the fault length with higher injection rate and increased effective normal stress. Injection-induced laboratory earthquakes with smaller nucleation lengths show generally higher peak slip rates and larger fault slip displacements, signifying an intensified seismic hazard. Our results also indicate that initially stable fault patches may be reactivated to slip seismically by increasing injection rates. This study systematically demonstrates that higher injection rates constitute dynamic loading, which increase the seismic hazard by shrinking the earthquake nucleation length.

Description: Fluids injection for hydraulic stimulation and fracturing, typical in the development of enhanced geothermal systems (EGS) in granites, can reactivate deep faults and induce seismicity. Such faults typically contain chlorite coatings as an alteration product that may impact styles of deformation—aseismic through seismic. We performed low velocity shear experiments on simulated granite fault gouges under conditions typifying a geothermal reservoir at ∼4-km depth with a confining pressure of 110 MPa, a temperature of 150°C, fluid pressures of 21–80 MPa, and chlorite contents of 0–100%, to investigate the influence of variation in effective stress and mineral composition on fault strength and stability. Our results show a transition from velocity-strengthening to velocity-weakening behavior in simulated granite gouge when the effective confining pressure was reduced from 89 to 30 MPa, characterized by a transition from fault compaction to dilation—as revealed by microstructural observations—with implications in enabling unstable failure. Conversely, increasing chlorite content stabilizes slip but reduces frictional strength. The microstructures of these mixed gouges exhibit shear localized on chlorite-enriched planes and promoting fault sliding. These results suggest that earthquake ruptures occurring during fluid injection can be facilitated by effective stress variations and that both controlling fluid overpressures (effective stresses) and being aware of the presence of alteration minerals are both important controls in mitigating such injection-induced seismic risks.