TY - JOUR A1 - Wang, Lei A1 - Kwiatek, Grzegorz A1 - Rybacki, Erik A1 - Bonnelye, Audrey A1 - Bohnhoff, Marco A1 - Dresen, Georg T1 - Laboratory study on fluid-induced fault slip behavior: the role of fluid pressurization rate JF - Geophysical research letters : GRL N2 - Understanding the physical mechanisms governing fluid-induced fault slip is important for improved mitigation of seismic risks associated with large-scale fluid injection. We conducted fluid-induced fault slip experiments in the laboratory on critically stressed saw-cut sandstone samples with high permeability using different fluid pressurization rates. Our experimental results demonstrate that fault slip behavior is governed by fluid pressurization rate rather than injection pressure. Slow stick-slip episodes (peak slip velocity < 4 mu m/s) are induced by fast fluid injection rate, whereas fault creep with slip velocity < 0.4 mu m/s mainly occurs in response to slow fluid injection rate. Fluid-induced fault slip may remain mechanically stable for loading stiffness larger than fault stiffness. Independent of fault slip mode, we observed dynamic frictional weakening of the artificial fault at elevated pore pressure. Our observations highlight that varying fluid injection rates may assist in reducing potential seismic hazards of field-scale fluid injection projects.
Plain Language Summary Human-induced earthquakes from field-scale fluid injection projects including enhanced geothermal system and deep wastewater injection have been documented worldwide. Although it is clear that fluid pressure plays a crucial role in triggering fault slip, the physical mechanism behind induced seismicity still remains poorly understood. We performed laboratory tests, and here we present two fluid-induced slip experiments conducted on permeable Bentheim sandstone samples crosscut by a fault that is critically stressed. Fault slip is then triggered by pumping the water from the bottom end of the sample at different fluid injection rates. Our results show that fault slip is controlled by fluid pressure increase rate rather than by the absolute magnitude of fluid pressure. In contrast to episodes of relatively rapid but stable sliding events caused by a fast fluid injection rate, fault creep is observed during slow fluid injection. Strong weakening of the dynamic friction coefficient of the experimental fault is observed at elevated pore pressure, independent of fault slip mode. These results may provide a better understanding of the complex behavior of fluid-induced fault slip on the field scale. KW - fault slip KW - fluid injection KW - induced seismicity KW - fluid pressurization KW - rate KW - stick-slip KW - fault creep Y1 - 2020 U6 - https://doi.org/10.1029/2019GL086627 SN - 0094-8276 SN - 1944-8007 VL - 47 IS - 6 PB - Wiley CY - Hoboken, NJ ER - TY - JOUR A1 - Richter, Gudrun A1 - Hainzl, Sebastian A1 - Dahm, Torsten A1 - Zöller, Gert T1 - Stress-based, statistical modeling of the induced seismicity at the Groningen gas field BT - the Netherlands JF - Environmental earth sciences N2 - Groningen is the largest onshore gas field under production in Europe. The pressure depletion of the gas field started in 1963. In 1991, the first induced micro-earthquakes have been located at reservoir level with increasing rates in the following decades. Most of these events are of magnitude less than 2.0 and cannot be felt. However, maximum observed magnitudes continuously increased over the years until the largest, significant event with ML=3.6 was recorded in 2014, which finally led to the decision to reduce the production. This causal sequence displays the crucial role of understanding and modeling the relation between production and induced seismicity for economic planing and hazard assessment. Here we test whether the induced seismicity related to gas exploration can be modeled by the statistical response of fault networks with rate-and-state-dependent frictional behavior. We use the long and complete local seismic catalog and additionally detailed information on production-induced changes at the reservoir level to test different seismicity models. Both the changes of the fluid pressure and of the reservoir compaction are tested as input to approximate the Coulomb stress changes. We find that the rate-and-state model with a constant tectonic background seismicity rate can reproduce the observed long delay of the seismicity onset. In contrast, so-called Coulomb failure models with instantaneous earthquake nucleation need to assume that all faults are initially far from a critical state of stress to explain the delay. Our rate-and-state model based on the fluid pore pressure fits the spatiotemporal pattern of the seismicity best, where the fit further improves by taking the fault density and orientation into account. Despite its simplicity with only three free parameters, the rate-and-state model can reproduce the main statistical features of the observed activity. KW - induced seismicity KW - modeling KW - statistical seismology KW - forecast Y1 - 2020 U6 - https://doi.org/10.1007/s12665-020-08941-4 SN - 1866-6280 SN - 1866-6299 VL - 79 IS - 11 PB - Springer CY - New York ER - TY - JOUR A1 - Wang, Lei A1 - Kwiatek, Grzegorz A1 - Rybacki, Erik A1 - Bohnhoff, Marco A1 - Dresen, Georg T1 - Injection-induced seismic moment release and laboratory fault slip BT - implications for fluid-induced seismicity JF - Geophysical research letters N2 - Understanding the relation between injection-induced seismic moment release and operational parameters is crucial for early identification of possible seismic hazards associated with fluid-injection projects. We conducted laboratory fluid-injection experiments on permeable sandstone samples containing a critically stressed fault at different fluid pressurization rates. The observed fluid-induced fault deformation is dominantly aseismic. Fluid-induced stick-slip and fault creep reveal that total seismic moment release of acoustic emission (AE) events is related to total injected volume, independent of respective fault slip behavior. Seismic moment release rate of AE scales with measured fault slip velocity. For injection-induced fault slip in a homogeneous pressurized region, released moment shows a linear scaling with injected volume for stable slip (steady slip and fault creep), while we find a cubic relation for dynamic slip. Our results highlight that monitoring evolution of seismic moment release with injected volume in some cases may assist in discriminating between stable slip and unstable runaway ruptures. KW - induced seismicity KW - seismic moment release KW - fluid injection KW - stick slip KW - fault creep KW - acoustic emission Y1 - 2020 U6 - https://doi.org/10.1029/2020GL089576 SN - 0094-8276 SN - 1944-8007 VL - 47 IS - 22 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Martínez-Garzón, Patricia A1 - Kwiatek, Grzegorz A1 - Bohnhoff, Marco A1 - Dresen, Georg T1 - Volumetric components in the earthquake source related to fluid injection and stress state JF - Geophysical research letters N2 - We investigate source processes of fluid-induced seismicity from The Geysers geothermal reservoir in California to determine their relation with hydraulic operations and improve the corresponding seismic hazard estimates. Analysis of 869 well-constrained full moment tensors (M-w 0.8-3.5) reveals significant non-double-couple components (>25%) for about 65% of the events. Volumetric deformation is governed by cumulative injection rates with larger non-double-couple components observed near the wells and during high injection periods. Source mechanisms are magnitude dependent and vary significantly between faulting regimes. Normal faulting events (M-w<2) reveal substantial volumetric components indicating dilatancy in contrast to strike-slip events that have a dominant double-couple source. Volumetric components indicating closure of cracks in the source region are mostly found for reverse faulting events with M-w>2.5. Our results imply that source processes and magnitudes of fluid-induced seismic events are strongly affected by the hydraulic operations, the reservoir stress state, and the faulting regime. KW - non-double-couple components KW - induced seismicity KW - geothermal KW - stress state KW - tensile opening KW - pore pressure Y1 - 2017 U6 - https://doi.org/10.1002/2016GL071963 SN - 0094-8276 SN - 1944-8007 VL - 44 IS - 2 SP - 800 EP - 809 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Kwiatek, Grzegorz A1 - Martinez-Garzon, Patricia A1 - Plenkers, K. A1 - Leonhardt, Maria A1 - Zang, Arno A1 - von Specht, Sebastian A1 - Dresen, Georg A1 - Bohnhoff, Marco T1 - Insights into complex subdecimeter fracturing processes occurring during a water injection experiment at depth in Aspo Hard Rock Laboratory, Sweden JF - Journal of geophysical research : Solid earth N2 - We investigate the source characteristics of picoseismicity (M-w < -2) recorded during a hydraulic fracturing in situ experiment performed in the underground Aspo Hard Rock Laboratory, Sweden. The experiment consisted of six stimulations driven by three different water injection schemes and was performed inside a 28-m-long, horizontal borehole located at 410-m depth. The fracturing processes were monitored with a variety of seismic networks including broadband seismometers, geophones, high-frequency accelerometers, and acoustic emission sensors thereby covering a wide frequency band between 0.01 and 100,000Hz. Here we study the high-frequency signals with dominant frequencies exceeding 1000 Hz. The combined seismic network allowed for detection and detailed analysis of 196 small-scale seismic events with moment magnitudes M-W < -3.5 (source sizes of decimeter scale) that occurred solely during the stimulations and shortly after. The double-difference relocated hypocenter catalog as well as source parameters were used to study the physical characteristics of the induced seismicity and then compared to the stimulation parameters. We observe a spatiotemporal migration of the picoseismic events away and toward the injection intervals in direct correlation with changes in the hydraulic energy (product of fluid injection pressure and injection rate). We find that the total radiated seismic energy is extremely low with respect to the product of injected fluid volume and pressure (hydraulic energy). The radiated seismic energy correlates well with the hydraulic energy rate. The obtained fault plane solutions for particularly well-characterized events signify the reactivation of preexisting rock defects under influence of increased pore fluid pressure on fault plane orientations in good correspondence with the local stress field orientation. KW - induced seismicity KW - fracking KW - picoseismicity KW - seismomechanics KW - source parameters KW - maximum magnitude Y1 - 2018 U6 - https://doi.org/10.1029/2017JB014715 SN - 2169-9313 SN - 2169-9356 VL - 123 IS - 8 SP - 6616 EP - 6635 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Hofmann, Hannes A1 - Zimmermann, Günter A1 - Farkas, Márton Pál A1 - Huenges, Ernst A1 - Zang, Arno A1 - Leonhardt, Maria A1 - Kwiatek, Grzegorz A1 - Martinez-Garzon, Patricia A1 - Bohnhoff, Marco A1 - Min, Ki-Bok A1 - Fokker, Peter A1 - Westaway, Rob A1 - Bethmann, Falko A1 - Meier, Peter A1 - Yoon, Kern Shin A1 - Choi, JaiWon A1 - Lee, Tae Jong A1 - Kim, Kwang Yeom T1 - First field application of cyclic soft stimulation at the Pohang Enhanced Geothermal System site in Korea JF - Geophysical journal international N2 - Large-magnitude fluid-injection induced seismic events are a potential risk for geothermal energy developments worldwide. One potential risk mitigation measure is the application of cyclic injection schemes. After validation at small (laboratory) and meso (mine) scale, the concept has now been applied for the first time at field scale at the Pohang Enhanced Geothermal System (EGS) site in Korea. From 7 August until 14 August 2017 a total of 1756 m(3) of surface water was injected into Pohang well PX-1 at flow rates between 1 and 10 l s(-1), with a maximum wellhead pressure (WHP) of 22.8 MPa, according to a site-specific cyclic soft stimulation schedule and traffic light system. A total of 52 induced microearthquakes were detected in real-time during and shortly after the injection, the largest of M-w 1.9. After that event a total of 1771 m(3) of water was produced back from the well over roughly 1 month, during which time no larger-magnitude seismic event was observed. The hydraulic data set exhibits pressure-dependent injectivity increase with fracture opening between 15 and 17 MPa WHP, but no significant permanent transmissivity increase was observed. The maximum magnitude of the induced seismicity during the stimulation period was below the target threshold of M-w 2.0 and additional knowledge about the stimulated reservoir was gained. Additionally, the technical feasibility of cyclic injection at field scale was evaluated. The major factors that limited the maximum earthquake magnitude are believed to be: limiting the injected net fluid volume, flowback after the occurrence of the largest induced seismic event, using a cyclic injection scheme, the application of a traffic light system, and including a priori information from previous investigations and operations in the treatment design. KW - Cyclic soft stimulation (CSS) KW - induced seismicity KW - risk mitigation KW - enhanced geothermal systems (EGS) KW - granite KW - Pohang (Korea) Y1 - 2019 U6 - https://doi.org/10.1093/gji/ggz058 SN - 0956-540X SN - 1365-246X VL - 217 IS - 2 SP - 926 EP - 949 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Bentz, Stephan A1 - Martinez-Garzon, Patricia A1 - Kwiatek, Grzegorz A1 - Dresen, Georg A1 - Bohnhoff, Marco T1 - Analysis of Microseismicity Framing M-L > 2.5 Earthquakes at The Geysers Geothermal Field, California JF - Journal of geophysical research : Solid earth N2 - 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 M-L > 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. KW - induced seismicity KW - earthquake nucleation KW - The Geysers KW - earthquake sequences Y1 - 2019 U6 - https://doi.org/10.1029/2019JB017716 SN - 2169-9313 SN - 2169-9356 VL - 124 IS - 8 SP - 8823 EP - 8843 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Dahm, Torsten A1 - Cesca, Simone A1 - Hainzl, Sebastian A1 - Braun, Thomas A1 - Krüger, Frank T1 - Discrimination between induced, triggered, and natural earthquakes close to hydrocarbon reservoirs: A probabilistic approach based on the modeling of depletion-induced stress changes and seismological source parameters JF - Journal of geophysical research : Solid earth N2 - Earthquakes occurring close to hydrocarbon fields under production are often under critical view of being induced or triggered. However, clear and testable rules to discriminate the different events have rarely been developed and tested. The unresolved scientific problem may lead to lengthy public disputes with unpredictable impact on the local acceptance of the exploitation and field operations. We propose a quantitative approach to discriminate induced, triggered, and natural earthquakes, which is based on testable input parameters. Maxima of occurrence probabilities are compared for the cases under question, and a single probability of being triggered or induced is reported. The uncertainties of earthquake location and other input parameters are considered in terms of the integration over probability density functions. The probability that events have been human triggered/induced is derived from the modeling of Coulomb stress changes and a rate and state-dependent seismicity model. In our case a 3-D boundary element method has been adapted for the nuclei of strain approach to estimate the stress changes outside the reservoir, which are related to pore pressure changes in the field formation. The predicted rate of natural earthquakes is either derived from the background seismicity or, in case of rare events, from an estimate of the tectonic stress rate. Instrumentally derived seismological information on the event location, source mechanism, and the size of the rupture plane is of advantage for the method. If the rupture plane has been estimated, the discrimination between induced or only triggered events is theoretically possible if probability functions are convolved with a rupture fault filter. We apply the approach to three recent main shock events: (1) the M-w 4.3 Ekofisk 2001, North Sea, earthquake close to the Ekofisk oil field; (2) the M-w 4.4 Rotenburg 2004, Northern Germany, earthquake in the vicinity of the Sohlingen gas field; and (3) the M-w 6.1 Emilia 2012, Northern Italy, earthquake in the vicinity of a hydrocarbon reservoir. The three test cases cover the complete range of possible causes: clearly human induced, not even human triggered, and a third case in between both extremes. KW - induced seismicity KW - probabilistic discrimination KW - hydrocarbon field KW - triggered earthquake KW - seismic hazard KW - earthquake Y1 - 2015 U6 - https://doi.org/10.1002/2014JB011778 SN - 2169-9313 SN - 2169-9356 VL - 120 IS - 4 SP - 2491 EP - 2509 PB - American Geophysical Union CY - Washington ER -