Refine
Year of publication
Language
- English (28) (remove)
Is part of the Bibliography
- yes (28)
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)
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.
The long-term temporal and spatial changes in statistical, source, and stress characteristics of one cluster of induced seismicity recorded at The Geysers geothermal field (U.S.) are analyzed in relation to the field operations, fluid migration, and constraints on the maximum likely magnitude. Two injection wells, Prati-9 and Prati-29, located in the northwestern part of the field and their associated seismicity composed of 1776 events recorded throughout a 7year period were analyzed. The seismicity catalog was relocated, and the source characteristics including focal mechanisms and static source parameters were refined using first-motion polarity, spectral fitting, and mesh spectral ratio analysis techniques. The source characteristics together with statistical parameters (b value) and cluster dynamics were used to investigate and understand the details of fluid migration scheme in the vicinity of injection wells. The observed temporal, spatial, and source characteristics were clearly attributed to fluid injection and fluid migration toward greater depths, involving increasing pore pressure in the reservoir. The seasonal changes of injection rates were found to directly impact the shape and spatial extent of the seismic cloud. A tendency of larger seismic events to occur closer to injection wells and a correlation between the spatial extent of the seismic cloud and source sizes of the largest events was observed suggesting geometrical constraints on the maximum likely magnitude and its correlation to the average injection rate and volume of fluids present in the reservoir.
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.
We analyze the spatiotemporal distribution of fault geometries from seismicity induced by fluid injection at The Geysers geothermal field. The consistency of these faults with the local stress field is investigated using (1) the fault instability coefficient I comparing the orientation of a fault with the optimal orientation for failure in the assumed stress field and (2) the misfit angle beta between slip vectors observed from focal mechanisms and predicted from stress tensor. A statistical approach is applied to calculate the most likely fault instabilities considering the uncertainties from focal mechanisms and stress inversion. We find that faults activated by fluid injection may display a broad range in orientations. About 72% of the analyzed seismicity occurs on faults with favorable orientation for failure with respect to the stress field. However, a number of events are observed either to occur on severely misoriented faults or to slip in a different orientation than predicted from stress field. These events mostly occur during periods of high injection rates and are located in proximity to the injection wells. From the stress inversion, the friction coefficient providing the largest overall instability is mu = 0.5. About 91% of the events are activated with an estimated excess pore pressure <10 MPa, in agreement with previous models considering the combined effect of thermal and poroelastic stress changes from fluid injection. Furthermore, high seismic activity and largest magnitudes occur on favorably oriented faults with large instability coefficients and low slip misfit angles.
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.
In this paper, an underground experiment at the Aspo Hard Rock Laboratory (HRL) is described. Main goal is optimizing geothermal heat exchange in crystalline rock mass at depth by multistage hydraulic fracturing with minimal impact on the environment, that is, seismic events. For this, three arrays with acoustic emission, microseismicity and electromagnetic sensors are installed mapping hydraulic fracture initiation and growth. Fractures are driven by three different water injection schemes (continuous, progressive and pulse pressurization). After a brief review of hydraulic fracture operations in crystalline rock mass at mine scale, the site geology and the stress conditions at Aspo HRL are described. Then, the continuous, single-flow rate and alternative, multiple-flow rate fracture breakdown tests in a horizontal borehole at depth level 410 m are described together with the monitoring networks and sensitivity. Monitoring results include the primary catalogue of acoustic emission hypocentres obtained from four hydraulic fractures with the in situ trigger and localizing network. The continuous versus alternative water injection schemes are discussed in terms of the fracture breakdown pressure, the fracture pattern from impression packer result and the monitoring at the arrays. An example of multistage hydraulic fracturing with several phases of opening and closing of fracture walls is evaluated using data from acoustic emissions, seismic broad-band recordings and electromagnetic signal response. Based on our limited amount of in situ tests (six) and evaluation of three tests in Avro granodiorite, in the multiple-flow rate test with progressively increasing target pressure, the acoustic emission activity starts at a later stage in the fracturing process compared to the conventional fracturing case with continuous water injection. In tendency, also the total number and magnitude of acoustic events are found to be smaller in the progressive treatment with frequent phases of depressurization.
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.
Rock deformation at depths in the Earth’s crust is often localized in high temperature shear zones occurring at different scales in a variety of lithologies. The presence of material heterogeneities is known to trigger shear zone development, but the mechanisms controlling initiation and evolution of localization are not fully understood. To investigate the effect of loading conditions on shear zone nucleation along heterogeneities, we performed torsion experiments under constant twist rate (CTR) and constant torque (CT) conditions in a Paterson-type deformation apparatus. The sample assemblage consisted of cylindrical Carrara marble specimens containing a thin plate of Solnhofen limestone perpendicular to the cylinder’s longitudinal axis. Under experimental conditions (900 °C, 400 MPa confining pressure), samples were plastically deformed and limestone is about 9 times weaker than marble, acting as a weak inclusion in a strong matrix. CTR experiments were performed at maximum bulk shear strain rates of ~ 2*10-4s-1, yielding peak shear stresses of ~ 20 MPa. CT tests were conducted at shear stresses of ~ 20 MPa resulting in bulk shear strain rates of 1-4*10-4s-1. Experiments were terminated at maximum bulk shear strains of ~ 0.3 and 1.0.Strain was localized within the Carrara marble in front of the inclusion in an area of strongly deformed grains and intense grain size reduction. Locally, evidences for coexisting brittle deformation are also observed regardless of the imposed loading conditions. The local shear strain at the inclusion tipis up to 30 times higher than the strain in the adjacent host rock, rapidly dropping to 5times higher at larger distance from the inclusion. At both investigated bulk strains, the evolution of microstructural and textural parameters is independent of loading conditions. Ourresults suggest that loading conditions do not significantly affect material heterogeneity-induced strain localization during its nucleation and transient stages.
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.
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.