Refine
Language
- English (26) (remove)
Is part of the Bibliography
- yes (26)
Keywords
- induced seismicity (5)
- strain localization (3)
- Fluid Flow (2)
- Ion-Beam (FIB) (2)
- central south island (2)
- fault creep (2)
- fluid injection (2)
- fluid-induced seismicity (2)
- frictional properties (2)
- geothermal (2)
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.
We study triggering processes in triaxial compression experiments under a constant displacement rate on sandstone and granite samples using spatially located acoustic emission events and their focal mechanisms. We present strong evidence that event-event triggering plays an important role in the presence of large-scale or macrocopic imperfections, while such triggering is basically absent if no significant imperfections are present. In the former case, we recover all established empirical relations of aftershock seismicity including the Gutenberg-Richter relation, a modified version of the Omori-Utsu relation and the productivity relation-despite the fact that the activity is dominated by compaction-type events and triggering cascades have a swarmlike topology. For the Gutenberg-Richter relations, we find that the b value is smaller for triggered events compared to background events. Moreover, we show that triggered acoustic emission events have a focal mechanism much more similar to their associated trigger than expected by chance.
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.
Natural earthquakes often have very few observable foreshocks which significantly complicates tracking potential preparatory processes. To better characterize expected preparatory processes before failures, we study stick-slip events in a series of triaxial compression tests on faulted Westerly granite samples. We focus on the influence of fault roughness on the duration and magnitude of recordable precursors before large stick-slip failure. Rupture preparation in the experiments is detectable over long time scales and involves acoustic emission (AE) and aseismic deformation events. Preparatory fault slip is found to be accelerating during the entire pre-failure loading period, and is accompanied by increasing AE rates punctuated by distinct activity spikes associated with large slip events. Damage evolution across the fault zones and surrounding wall rocks is manifested by precursory decrease of seismic b-values and spatial correlation dimensions. Peaks in spatial event correlation suggest that large slip initiation occurs by failure of multiple asperities. Shear strain estimated from AE data represents only a small fraction (< 1%) of total shear strain accumulated during the preparation phase, implying that most precursory deformation is aseismic. The relative contribution of aseismic deformation is amplified by larger fault roughness. Similarly, seismic coupling is larger for smooth saw-cut faults compared to rough faults. The laboratory observations point towards a long-lasting and continuous preparation process leading to failure and large seismic events. The strain partitioning between aseismic and observable seismic signatures depends on fault structure and instrument resolution.
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.