TY - JOUR A1 - Wang, Lei A1 - Rybacki, Erik A1 - Bonnelye, Audrey A1 - Bohnhoff, Marco A1 - Dresen, Georg T1 - Experimental investigation on static and dynamic bulk moduli of dry and fluid-saturated porous sandstones JF - Rock mechanics and rock engineering N2 - 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. KW - Porous sandstone KW - Static bulk modulus KW - Dynamic bulk modulus KW - Narrow KW - (compliant) cracks KW - Equant (stiff) pores Y1 - 2020 U6 - https://doi.org/10.1007/s00603-020-02248-3 SN - 0723-2632 SN - 1434-453X VL - 54 IS - 1 SP - 129 EP - 148 PB - Springer CY - Wien ER - TY - JOUR A1 - Bentz, Stephan A1 - Kwiatek, Grzegorz A1 - Martinez-Garzon, Patricia A1 - Bohnhoff, Marco A1 - Dresen, Georg T1 - Seismic moment evolution during hydraulic stimulations JF - Geophysical research letters N2 - Analysis of past and present stimulation projects reveals that the temporal evolution and growth of maximum observed moment magnitudes may be linked directly to the injected fluid volume and hydraulic energy. Overall evolution of seismic moment seems independent of the tectonic stress regime and is most likely governed by reservoir specific parameters, such as the preexisting structural inventory. Data suggest that magnitudes can grow either in a stable way, indicating the constant propagation of self-arrested ruptures, or unbound, for which the maximum magnitude is only limited by the size of tectonic faults and fault connectivity. Transition between the two states may occur at any time during injection or not at all. Monitoring and traffic light systems used during stimulations need to account for the possibility of unstable rupture propagation from the very beginning of injection by observing the entire seismicity evolution in near-real time and at high resolution for an immediate reaction in injection strategy. Plain Language Summary Predicting and controlling the size of earthquakes caused by fluid injection is currently the major concern of many projects associated with geothermal energy production. Here, we analyze the magnitude and seismic moment evolution with injection parameters for prominent geothermal and scientific projects to date. Evolution of seismicity seems to be largely independent of the tectonic stress background and seemingly depends on reservoir specific characteristics. We find that the maximum observed magnitudes relate linearly to the injected volume or hydraulic energy. A linear relation suggests stable growth of induced ruptures, as predicted by current models, or rupture growth may no longer depend on the stimulated volume but on tectonics. A system may change between the two states during the course of fluid injection. Close-by and high-resolution monitoring of seismic and hydraulic parameters in near-real time may help identify these fundamental changes in ample time to change injection strategy and manage maximum magnitudes. Y1 - 2020 U6 - https://doi.org/10.1029/2019GL086185 SN - 0094-8276 SN - 1944-8007 VL - 47 IS - 5 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Wang, Lei A1 - Dresen, Georg A1 - Rybacki, Erik A1 - Bonnelye, Audrey A1 - Bohnhoff, Marco T1 - Pressure-dependent bulk compressibility of a porous granular material modeled by improved contact mechanics and micromechanical approaches BT - effects of surface roughness of grains JF - Acta materialia N2 - The change of the mechanical properties of granular materials with pressure is an important topic associated with many industrial applications. In this paper we investigate the influence of hydrostatic pressure (P-e) on the effective bulk compressibility (C-eff) of a granular material by applying two modified theoretical approaches that are based on contact mechanics and micromechanics, respectively. For a granular material composed of rough grains, an extended contact model is developed to elucidate the effect of roughness of grain surfaces on bulk compressibility. At relatively low pressures, the model predicts that the decrease of bulk compressibility with pressure may be described by a power law with an exponent of -1/2 (i.e., C-eff proportional to P-e(1/2) ), but deviates at intermediate pressures. At elevated pressures beyond full contact, bulk compressibility remains almost unchanged, which may be roughly evaluated by continuum contact mechanics. As an alternative explanation of pressure-dependent bulk compressibility, we suggest a micromechanical model that accounts for effects of different types of pore space present in granular materials. Narrow and compliant inter-granular cracks are approximated by three-dimensional oblate spheroidal cracks with rough surfaces, whereas the equant and stiff pores surrounded by three and four neighboring grains are modeled as tubular pores with cross sections of three and four cusp-like corners, respectively. In this model, bulk compressibility is strongly reduced with increasing pressure by progressive closure of rough-walled cracks. At pressures exceeding crack closure pressure, deformation of the remaining equant pores is largely insensitive to pressure, with almost no further change in bulk compressibility. To validate these models, we performed hydrostatic compression tests on Bentheim sandstone (a granular rock consisting of quartz with high porosity) under a wide range of pressure. The relation between observed microstructures and measured pressure-dependent bulk compressibility is well explained by both suggested models. KW - Bulk compressibility KW - Granular materials KW - Roughness KW - Micromechanical model KW - Contact model Y1 - 2020 U6 - https://doi.org/10.1016/j.actamat.2020.01.063 SN - 1359-6454 SN - 1873-2453 VL - 188 SP - 259 EP - 272 PB - Elsevier CY - Amsterdam ER - 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 - 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 - Durand, Virginie A1 - Bentz, Stephan A1 - Kwiatek, Grzegorz A1 - Dresen, Georg A1 - Wollin, Christopher A1 - Heidbach, Oliver A1 - Martinez-Garzon, Patricia A1 - Cotton, Fabrice A1 - Nurlu, Murat A1 - Bohnhoff, Marco T1 - A two-scale preparation phase preceded an M-w 5.8 earthquake in the sea of marmara offshore Istanbul, Turkey JF - Seismological research letters N2 - We analyze the spatiotemporal evolution of seismicity during a sequence of moderate (an M-w 4.7 foreshock and M-w 5.8 mainshock) earthquakes occurring in September 2019 at the transition between a creeping and a locked segment of the North Anatolian fault in the central Sea of Marmara, northwest Turkey. To investigate in detail the seismicity evolution, we apply a matched-filter technique to continuous waveforms, thus reducing the magnitude threshold for detection. Sequences of foreshocks preceding the two largest events are clearly seen, exhibiting two different behaviors: a long-term activation of the seismicity along the entire fault segment and a short-term concentration around the epicenters of the large events. We suggest a two-scale preparation phase, with aseismic slip preparing the mainshock final rupture a few days before, and a cascade mechanism leading to the nucleation of the mainshock. Thus, our study shows a combination of seismic and aseismic slip during the foreshock sequence changing the strength of the fault, bringing it closer to failure. Y1 - 2020 U6 - https://doi.org/10.1785/0220200110 SN - 0895-0695 SN - 1938-2057 VL - 91 IS - 6 SP - 3139 EP - 3147 CY - Boulder 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 - Kwiatek, Grzegorz A1 - Saarno, Tero A1 - Ader, Thomas A1 - Blümle, Felix A1 - Bohnhoff, Marco A1 - Chendorain, Michael A1 - Dresen, Georg A1 - Heikkinen, Pekka A1 - Kukkonen, Ilmo A1 - Leary, Peter A1 - Leonhardt, Maria A1 - Malin, Peter A1 - Martinez-Garzon, Patricia A1 - Passmore, Kevin A1 - Passmore, Paul A1 - Valenzuela, Sergio A1 - Wollin, Christopher T1 - Controlling fluid-induced seismicity during a 6.1-km-deep geothermal stimulation in Finland JF - Science Advances N2 - We show that near-real-time seismic monitoring of fluid injection allowed control of induced earthquakes during the stimulation of a 6.1-km-deep geothermal well near Helsinki, Finland. A total of 18,160 m(3) of fresh water was pumped into crystalline rocks over 49 days in June to July 2018. Seismic monitoring was performed with a 24-station borehole seismometer network. Using near-real-time information on induced-earthquake rates, locations, magnitudes, and evolution of seismic and hydraulic energy, pumping was either stopped or varied-in the latter case, between well-head pressures of 60 and 90 MPa and flow rates of 400 and 800 liters/min. This procedure avoided the nucleation of a project-stopping magnitude M-W 2.0 induced earthquake, a limit set by local authorities. Our results suggest a possible physics-based approach to controlling stimulation-induced seismicity in geothermal projects. Y1 - 2019 U6 - https://doi.org/10.1126/sciadv.aav7224 SN - 2375-2548 VL - 5 IS - 5 PB - American Association for the Advancement of Science CY - Washington ER - TY - JOUR A1 - Martinez-Garvon, Patricia A1 - Bohnhoff, Marco A1 - Mencin, David A1 - Kwiatek, Grzegorz A1 - Dresen, Georg A1 - Hodgkinson, Kathleen A1 - Nurlu, Murat A1 - Kadirioglu, Filiz Tuba A1 - Kartal, Recai Feyiz T1 - Slow strain release along the eastern Marmara region offshore Istanbul in conjunction with enhanced local seismic moment release JF - Earth & planetary science letters N2 - We analyze a large transient strainmeter signal recorded at 62.5 m depth along the southern shore of the eastern Sea of Marmara region in northwestern Turkey. This region represents a passage of stress transfer from the Izmit rupture to the Marmara seismic gap. The strain signal was recorded at the Esenkoy site by one of the ICDP-GONAF (International Continental Drilling Programme - Geophysical Observatory at the North Anatolian Fault) strainmeters on the Armutlu peninsula with a maximum amplitude of 5 microstrain and lasting about 50 days. The onset of the strain signal coincided with the origin time of a M-w 4.4 earthquake offshore Yalova, which occurred as part of a seismic sequence including eight M-w >= 3.5 earthquakes. The Mw 4.4 event occurred at a distance of about 30 km from Esenkoy on June 25th 2016 representing the largest earthquake in this region since 2008. Before the event, the maximum horizontal strain was subparallel to the regional maximum horizontal stress derived from stress inversion of local seismicity. During the strain transient, we observe a clockwise rotation in the local horizontal strain field of about 20 degrees. The strain signal does not correlate with known environmental parameters such as annual changes of sea level, rainfall or temperature. The strain signal could indicate local slow slip on the Cinarcik fault and thus a transfer of stress to the eastern Marmara seismic gap. KW - strain transient KW - slow slip KW - transform faults KW - seismic hazard KW - strainmeter data KW - Sea of Marmara Y1 - 2019 U6 - https://doi.org/10.1016/j.epsl.2019.01.001 SN - 0012-821X SN - 1385-013X VL - 510 SP - 209 EP - 218 PB - Elsevier CY - Amsterdam 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 -