Refine
Language
- English (20)
Is part of the Bibliography
- yes (20)
Keywords
- induced seismicity (4)
- strain localization (3)
- Fluid Flow (2)
- Ion-Beam (FIB) (2)
- central south island (2)
- fluid-induced seismicity (2)
- frictional properties (2)
- geothermal (2)
- hanging wall (2)
- internal structure (2)
Institute
- Institut für Geowissenschaften (20) (remove)
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.
Localization processes in the viscous lower crust generate ductile shear zones over a broad range of scales affecting long‐term lithosphere deformation and the mechanical response of faults during the seismic cycle. Here we use centimeter‐scale numerical models in order to gain detailed insight into the processes involved in strain localization and rheological weakening in viscously deforming rocks. Our 2‐D Cartesian models are benchmarked to high‐temperature and high‐pressure torsion experiments on Carrara marble samples containing a single weak Solnhofen limestone inclusion. The models successfully reproduce bulk stress‐strain transients and final strain distributions observed in the experiments by applying a simple, first‐order softening law that mimics rheological weakening. We find that local stress concentrations forming at the inclusion tips initiate strain localization inside the host matrix. At the tip of the propagating shear zone, weakening occurs within a process zone, which expands with time from the inclusion tips toward the matrix. Rheological weakening is a precondition for shear zone localization, and the width of this shear zone is found to be controlled by the degree of softening. Introducing a second softening step at elevated strain, a high strain layer develops inside the localized shear zone, analogous to the formation of ultramylonite bands in mylonites. These results elucidate the transient evolution of stress and strain rate during inception and maturation of ductile shear zones.
The spatiotemporal, kinematic, and source characteristics of induced seismicity occurring at different fluid injection rates are investigated to determine the predominant physical mechanisms responsible for induced seismicity at the northwestern part of The Geysers geothermal field, California. We analyze a relocated hypocenter catalog from a seismicity cluster where significant variations of the stress tensor orientation were previously observed to correlate with injection rates. We find that these stress tensor orientation changes may be related to increased pore pressure and the corresponding changes in poroelastic stresses at reservoir depth. Seismic events during peak injections tend to occur at greater distances from the injection well, preferentially trending parallel to the maximum horizontal stress direction. In contrast, at lower injection rates the seismicity tends to align in a different direction which suggests the presence of a local fault. During peak injection intervals, the relative contribution of strike-slip faulting mechanisms increases. Furthermore, increases in fluid injection rates also coincide with a decrease in b values. Our observations suggest that regardless of the injection stage, most of the induced seismicity results from thermal fracturing of the reservoir rock. However, during peak injection intervals, the increase in pore pressure may likewise be responsible for the induced seismicity. By estimating the thermal and hydraulic diffusivities of the reservoir, we confirm that the characteristic diffusion length for pore pressure is much greater than the corresponding length scale for temperature and also more consistent with the spatial extent of seismicity observed during different injection rates.
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.
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.
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.
On 27 December 2007, a 1.9 seismic event occurred within a dyke in the deep-level Mponeng Gold Mine, South Africa. From the seismological network of the mine and the one from the Japanese-German Underground Acoustic Emission Research in South Africa (JAGUARS) group, the hypocentral depth (3,509 m), focal mechanism and aftershock location were estimated. Since no mining activity took place in the days before the event, dynamic triggering due to blasting can be ruled out as the cause. To investigate the hypothesis that stress transfer, due to excavation of the gold reef, induced the event, we set up a small-scale high-resolution three-dimensional (3D) geomechanical numerical model. The model consisted of the four different rock units present in the mine: quartzite (footwall), hard lava (hanging wall), conglomerate (gold reef) and diorite (dykes). The numerical solution was computed using a finite-element method with a discretised mesh of approximately elements. The initial stress state of the model is in agreement with in situ data from a neighbouring mine, and the step-wise excavation was simulated by mass removal from the gold reef. The resulting 3D stress tensor and its changes due to mining were analysed based on the Coulomb failure stress changes on the fault plane of the event. The results show that the seismic event was induced regardless of how the Coulomb failure stress changes were calculated and of the uncertainties in the fault plane solution. We also used the model to assess the seismic hazard due to the excavation towards the dyke. The resulting curve of stress changes shows a significant increase in the last in front of the dyke, indicating that small changes in the mining progress towards the dyke have a substantial impact on the stress transfer.
The Alpine Fault is capable of generating large (MW > 8) earthquakes and is the main geohazard on South Island, NZ, and late in its 250–291-year seismic cycle. To minimize its hazard potential, it is indispensable to identify and understand the processes influencing the geomechanical behavior and strength-evolution of the fault. High-resolution microstructural, mineralogical and geochemical analyses of the Alpine Fault's core demonstrate wall rock fragmentation, assisted by mineral dissolution, and cementation resulting in the formation of a fine-grained principal slip zone (PSZ). A complex network of anastomosing and mutually cross-cutting calcite veins implies that faulting occurred during episodes of dilation, slip and sealing. Fluid-assisted dilatancy leads to a significant volume increase accommodated by vein formation in the fault core. Undeformed euhedral chlorite crystals and calcite veins that have cut footwall gravels demonstrate that these processes occurred very close to the Earth's surface. Microstructural evidence indicates that cataclastic processes dominate the deformation and we suggest that powder lubrication and grain rolling, particularly influenced by abundant nanoparticles, play a key role in the fault core's velocity-weakening behavior rather than frictional sliding. This is further supported by the absence of smectite, which is reasonable given recently measured geothermal gradients of more than 120 °C km−1 and the impermeable nature of the PSZ, which both limit the growth of this phase and restrict its stability to shallow depths. Our observations demonstrate that high-temperature fluids can influence authigenic mineral formation and thus control the fault's geomechanical behavior and the cyclic evolution of its strength.
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.
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.