Refine
Has Fulltext
- no (18)
Document Type
- Article (18) (remove)
Language
- English (18)
Is part of the Bibliography
- yes (18) (remove)
Keywords
- induced seismicity (6)
- fault creep (2)
- fluid injection (2)
- fluid-induced seismicity (2)
- geothermal (2)
- maximum magnitude (2)
- pore pressure (2)
- reservoir characterization (2)
- source parameters (2)
- 3D geomechanical numerical model (1)
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.
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.
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.
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.
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.
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 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.
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.