Injection-induced fault slip and associated seismicity in the lab: Insights from source mechanisms, local stress states and fault geometry

Probing source mechanisms of natural and induced earthquakes is a powerful tool to unveil associated rupture kinematics. The source processes of failure and slip instability driven by stress loading are affected by fault geometry, but the source ruptures of injection-induced seismicity in relation to fault structures and local stress states remain poorly understood. We have conducted a series of fault reactivation and slip experiments on sandstone samples containing faults with different surface roughness (smooth saw-cut fault and fractured rough fault). We impose progressive fluid injection to induce fault slip, and simultaneously monitor the associated acoustic emission (AE) activity. Using high-resolution AE recordings, we perform full moment tensor inversion of all located AE sources, and investigate the changes of AE source characteristics associated with induced fault slip and their relation to fault roughness. For the complex and rough fault, we observe significant non-double-couple components of AE sources and a high degree of focal mechanism heterogeneity. The temporal changes of AE mechanisms associated with injection-induced fault slip on the smooth fault reveal increasing proportions of double-couple components and decreasing variability of AE focal mechanisms when approaching the onset of slip events. The observed inconsistency between the nodal planes of AE sources and the macroscopic fault plane orientation is attributed to the development of secondary fracture networks surrounding the principal slip surface. We analyze changes in the magnitude-frequency characteristics and source mechanisms of AEs with fault-normal distance, showing that for the smooth (mature) fault, Gutenberg–Richter b-value of on-fault seismicity is lower and focal mechanisms are less heterogeneous, compared to off-fault seismicity. Our results emphasize the important role of roughness-related changes in local fault geometry and associated stress heterogeneity for source mechanisms and rupture kinematics of injection-induced seismicity

Experimental Investigation on Static and Dynamic Bulk Moduli of Dry and Fluid-Saturated Porous Sandstones

Knowledge of pressure-dependent static and dynamic moduli of porous reservoir rocks is of key importance for evaluating geological setting of a reservoir in geo-energy applications. We examined experimentally the evolution of static and dynamic bulk moduli for porous Bentheim sandstone with increasing confining pressure up to about 190 MPa under dry and water-saturated conditions. The static bulk moduli (K-s) were estimated from stress-volumetric strain curves while dynamic bulk moduli (K-d) were derived from the changes in ultrasonic P- and S- wave velocities (similar to 1 MHz) along different traces, which were monitored simultaneously during the entire deformation. In conjunction with published data of other porous sandstones (Berea, Navajo and Weber sandstones), our results reveal that the ratio between dynamic and static bulk moduli (K-d/K-s) reduces rapidly from about 1.5 - 2.0 at ambient pressure to about 1.1 at high pressure under dry conditions and from about 2.0 - 4.0 to about 1.5 under water-saturated conditions, respectively. We interpret such a pressure-dependent reduction by closure of narrow (compliant) cracks, highlighting thatK(d)/K(s)is positively correlated with the amount of narrow cracks. Above the crack closure pressure, where equant (stiff) pores dominate the void space,K-d/K(s)is almost constant. The enhanced difference between dynamic and static bulk moduli under water saturation compared to dry conditions is possibly caused by high pore pressure that is locally maintained if measured using high-frequency ultrasonic wave velocities. In our experiments, the pressure dependence of dynamic bulk modulus of water-saturated Bentheim sandstone at effective pressures above 5 MPa can be roughly predicted by both the effective medium theory (Mori-Tanaka scheme) and the squirt-flow model. Static bulk moduli are found to be more sensitive to narrow cracks than dynamic bulk moduli for porous sandstones under dry and water-saturated conditions

Analysis of microseismicity framing ML > 2.5 earthquakes at The Geysers geothermal field, California

Preparatory mechanisms accompanying or leading to nucleation of larger earthquakes have been observed at both laboratory and field scales, but conditions favoring the occurrence of observable preparatory processes are still largely unknown. In particular, it remains a matter of debate why some earthquakes occur spontaneously without noticeable precursors as opposed to events that are preceded by an extended failure process. In this study, we have generated new high‐resolution seismicity catalogs framing the occurrence of 20 ML > 2.5 earthquakes at The Geysers geothermal field in California. To this end, a seismicity catalog of the 11 days framing each large event was created. We selected 20 sequences sampling different hypocentral depths and hydraulic conditions within the field. Seismic activity and magnitude frequency distributions displayed by the different earthquake sequences are correlated with their location within the reservoir. Sequences located in the northwestern part of the reservoir show overall increased seismic activity and low b values, while the southeastern part is dominated by decreased seismic activity and higher b values. Periods of high injection coincide with high b values and vice versa. These observations potentially reflect varying differential and mean stresses and damage of the reservoir rocks across the field. About 50% of analyzed sequences exhibit no change in seismicity rate in response to the large main event. However, we find complex waveforms at the onset of the main earthquake, suggesting that small ruptures spontaneously grow into or trigger larger events

Metre-scale damage zone characterization using S-coda waves from active ultrasonic transmission measurements in the STIMTEC project, URL Reiche Zeche, Germany

Studies of controlled hydraulic stimulation experiments with active and passive seismic monitoring conducted in Underground Research Laboratories (URLs) benefit from specific knowledge of hydraulic parameters, close by microseismic monitoring revealing structural details of the rock mass, and detailed evolution of seismicity in response to injection operations. Microseismic monitoring is commonly used to characterize a stimulated reservoir volume, for example, in terms of damage evolution of the rock mass. Since seismic attenuation is affected by damage of the rock volume, active seismic sources covering sizes from the centimetre to decimetre scale may help us to investigate space–time varying attenuation properties in a reservoir. This may allow us to monitor damage evolution of the stimulated rock volume in more detail, also since active seismic sources produce stronger signals leading to a broader frequency range that can be analysed compared to passive seismic signals. Within the STIMTEC project in the URL Reiche Zeche (URL-RZ) in Freiberg (Germany), more than 300 active Ultrasonic Transmission (UT) measurements were performed before and after hydraulic stimulations in two boreholes in the targeted rock volume, an anisotropic metamorphic gneiss. The signal-frequency content ranges between 1 and 60 kHz. Assuming scattering attenuation to dominate over intrinsic attenuation, we here apply the single isotropic scattering model. S-coda waves of 88 spatially representative UT measurements are used to estimate the coda quality factor (QC). We obtain stable QC estimates for centre frequencies of octave-width frequency bands between 3 and 21 kHz. We group neighbouring UT measurements to stabilize the observations and form eight UT groups in total, covering different depth intervals in three boreholes and four different time periods to investigate scattering attenuation changes in a spatiotemporal manner. Our final mean QC (⁠QC¯¯¯¯¯¯¯⁠) estimates show characteristic frequency-dependence as observed at the field scale in geological reservoirs. We find temporal variations of QC are strongly connected to hydraulic stimulation, and these variations are more significant than those resolved from velocity changes. QC¯¯¯¯¯¯¯ estimates at frequencies above 15 kHz indicate healing of injection-induced small-scale fractures during a two-months post-stimulation phase. Larger fractures, mostly sampled by lower frequencies (<15 kHz), seem to be more persistent with time (over 15 months). We observe spatial differences of QC¯¯¯¯¯¯¯ values near the mine galleries (driftway and vein drift) and relate these observations to different extents and characteristics of the galleries’ excavation damage zones. Our results further support previous assumptions based on borehole televiewer logs and mapped structures of an existing fault with larger damage zone that crosses the stimulated rock volume NW-SE between the galleries. We conclude that the coda analysis of active UT measurements complements established imaging methods used during experiments in URLs. In particular, coda analysis is a powerful tool for the detection of damage zones and for monitoring local fracture networks with immediate application for imaging georeservoirs considered for exploitation or underground storage of gases and liquids

Unsupervised clustering of catalogue-driven features for characterizing temporal evolution of labquake stress

Earthquake forecasting poses significant challenges, especially due to the elusive nature of stress states in fault systems. To tackle this problem, we use features derived from seismic catalogues obtained from acoustic emission (AE) signals recorded during triaxial stick-slip experiments on natural fractures in three Westerly granite samples. We extracted 47 physically explainable features from AE data that described spatio-temporal evolution of stress and damage in the vicinity of the fault surface. These features are then subjected to unsupervised clustering using the K-means method, revealing three distinct stages with a proper agreement with the temporal evolution of stress. The recovered stages correspond to the mechanical behaviour of the rock, characterized as initial stable (elastic) deformation, followed by a transitional stage leading to an unstable deformation prior to failure. Notably, AE rate, clustering-localization features, fractal dimension, b-value, interevent time distribution, and correlation integral are identified as significant features for the unsupervised clustering. The systematically evolving stages can provide valuable insights for characterizing preparatory processes preceding earthquake events associated with geothermal activities and waste-water injections. In order to address the upscaling issue, we propose to use the most important features and, in case of normalization challenge, removing non-universal features, such as AE rate. Our findings hold promise for advancing earthquake prediction methodologies based on laboratory experiments and catalogue-driven features

What controls the presence and characteristics of aftershocks in rock fracture in the lab?

Aftershock cascades are a characteristic feature of natural seismicity, but underlying mechanisms remain debated. Here, we experimentally explore the presence or absence of aftershocks during failure of intact rock and slip on newly created laboratory faults. We show that the overall activity increase and spatial localization of acoustic emission (AE) events during fracture nucleation occurs without temporal (Omori-type) correlations. Our analysis shows that this absence of aftershock sequences occurs even beyond peak stress and also when a macroscopic fracture has formed post peak-stress and propagates. Instead, aftershock triggering does occur during post-fracture stress relaxation along the newly created lab-fault and in the presence of large-scale stress heterogeneities, for example, imposed by a saw-cut notch. The detected aftershocks in these cases can be described by standard seismological relationships such as a modified Omori-Utsu relation and its associated inter-event time distribution and productivity relation. Moreover, AE within all experiments follow the Gutenberg-Richter relation, with smaller b-values for triggered events compared to non-triggered events. Performing full-moment tensor inversions, we find that seismic events with significant isotropic, compaction components play an important role for aftershock triggering. The resulting triggered events tend to have focal mechanisms similar to their trigger. Seismic events with predominant tensile components, on the other hand, show little evidence for aftershock triggering. This opens up a new perspective on aftershocks, going beyond not only the rate-and-state paradigm limited to purely frictional sliding or shear events but also the mainshock attribute paradigm where mainshock attributes control aftershock patterns

Stress drop-magnitude dependence of acoustic emissions during laboratory stick-slip

Earthquake source parameters such as seismic stress drop and corner frequency are observed to vary widely, leading to persistent discussion on potential scaling of stress drop and event size. Physical mechanisms that govern stress drop variations arc difficult to evaluate in nature and are more readily studied in controlled laboratory experiments. We perform two stick-slip experiments on fractured (rough) and cut (smooth) Westerly granite samples to explore fault roughness effects on acoustic emission (AE) source parameters. We separate large stick-slip events that generally saturate the seismic recording system from populations of smaller AE events which are sensitive to fault stresses prior to slip. AE event populations show many similarities to natural seismicity and may be interpreted as laboratory equivalent of natural microseismic events. We then compare the temporal evolution of mechanical data such as measured stress release during slip to temporal changes in stress drops derived from Alis using the spectral ratio technique. We report on two primary observations: (1) In contrast to most case studies for natural earthquakes, we observe a strong increase in seismic stress drop with AE size. (2) The scaling of stress drop with magnitude is governed by fault roughness, whereby the rough fault shows a more rapid increase of the stress drop magnitude relation with progressing large stick-slip events than the smooth fault. The overall range of AE sizes on the rough surface is influenced by both the average grain size and the width of the fault core. The magnitudes of the smallest AE events on smooth faults may also be governed by grain size. However, AEs significantly grow beyond peak roughness and the width of the fault core. Our laboratory tests highlight that source parameters vary substantially in the presence of fault zone heterogeneity (i.e. roughness and narrow grain size distribution), which may affect seismic energy partitioning and static stress drops of small and large AE events

Strain Partitioning and Frictional Behavior of Opalinus Clay During Fault Reactivation

The Opalinus Clay (OPA) formation is considered a suitable host rock candidate for nuclear waste storage. However, the sealing integrity and long-term safety of OPA are potentially compromised by pre-existing natural or artificially induced faults. Therefore, characterizing the mechanical behavior and microscale deformation mechanisms of faults and the surrounding rock is relevant for predicting repository damage evolution. In this study, we performed triaxial tests using saw-cut samples of the shaly and sandy facies of OPA to investigate the influence of pressure and mineral composition on the deformation behavior during fault reactivation. Dried samples were hydrostatically pre-compacted at 50 MPa and then deformed at constant strain rate, drained conditions and confining pressures (pc) of 5–35 MPa. Mechanical data from triaxial tests was complemented by local strain measurements to determine the relative contribution of bulk deformation and fault slip, as well as by acoustic emission (AE) monitoring, and elastic P-wave velocity measurements using ultrasonic transmissions. With increasing pc, we observe a transition from brittle deformation behavior with highly localized fault slip to semi-brittle behavior characterized by non-linear strain hardening with increasing delocalization of deformation. We find that brittle localization behavior is limited by pc at which fault strength exceeds matrix yield strength. AEs were only detected in tests performed on sandy facies samples, and activity decreased with increasing pc. Microstructural analysis of deformed samples revealed a positive correlation between increasing pc and gouge layer thickness. This goes along with a change from brittle fragmentation and frictional sliding to the development of shear zones with a higher contribution of cataclastic and granular flow. Friction coefficient at fault reactivation is only slightly higher for the sandy (µ ~ 0.48) compared to the shaly facies (µ ~ 0.4). Slide-hold-slide tests performed after ~ 6 mm axial shortening suggest stable creeping and long-term weakness of faults at the applied conditions. Our results demonstrate that the mode of fault reactivation highly depends on the present stress field and burial history