Fluid-injection into the subsurface in the frame of reservoir-engineering activities for hydrocarbon and geothermal energy production has resulted in the reactivation of previously unknown critically stressed and thus hazardous faults. In Oklahoma and Kansas (US), a dramatic increase of induced seismicity including four M>5 earthquakes has been seen over the last years. This increased the seismic hazard of the nearby settlements and key-infrastructure, triggering significant public awareness and acceptance issues. SAIDAN aims at investigating the earthquake physics and quantifying the proportion of seismic versus aseismic slip release in fluid-injection environments and investigating the resulting maximum magnitude of the induced earthquakes and the subsequent hazard. To do that, geothermal environments such as The Geysers and Salton Sea are investigated, as well as regions where fluid-induced seismicity is prominent, such as in southern Kansas.
Some of the highlights from our recent studies are the following:
Earthquake clustering properties are investigated in relation to fluid balance H(t) (the difference of fluid injection and production rates) using about nine years of data from The Geysers (both the entire field and a local subset), Coso, and Salton Sea geothermal fields in California. Individual earthquake clusters are identified and classified using the nearest-neighbor approach of Zaliapin and Ben-Zion (2013a, 2013b). These are used to calculate nine complementary cluster statistics as time series with a step of about one month. Three alternative techniques (moving window correlation, analysis of variance, and regression) are employed to assess the relations between (possibly nonstationary) time series of cluster statistics and H(t). A total of 108 pairwise relations between cluster statistics and H(t) are analyzed to clarify effects of fluid activities on seismicity in different places. The seismic clustering response to the fluid balance differs among the examined fields. The Geysers and Salton Sea areas display the highest and lowest clustering responses, respectively. The proportion of clusters consisting of a single event with no offspring (singles) is correlated significantly with H(t) at all examined data sets, with a lower proportion of singles during periods of high fluid balance. This may reflect increased susceptibility to earthquake triggering in time intervals with high injection rates. The background seismicity rates significantly increase with H(t) at the Geysers and Coso, while an opposite relation holds at the Salton Sea. This could be related to the high structural and tectonic complexity at the Salton Sea compared to the other two geothermal fields.
We investigated source mechanisms of 80 earthquakes at the Salton Sea geothermal field in southern California. The performance of two moment tensor inversion techniques (P-wave first-motion amplitudes vs full waveforms) was compared. Only a combination of P-wave amplitude inversion and moment tensor refinement techniques enables a robust determination of the source types. We find that most shallow earthquakes display statistically significant non-double-couple components. Largest volumetric components are inferred for events located near active injection and production wells and primarily correspond to a normal-faulting regime. In contrast, strike-slip events are dominated by shear motion, and their locations extend well beyond the reservoir.
Analysis of 869 well-constrained full moment tensors of fluid-induced seismicity (Mw 0.8–3.5) from The Geysers geothermal reservoir in northern California were investigated revealed 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. 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.
Spatio-temporal trends in statistical, source, and stress characteristics of 1776 induced seismic events from The Geysers geothermal field were analyzed for a period of 7 years. The temporal, spatial, and source characteristics were attributed to changes in the fluid injection rates and the fluid migration toward greater depths, involving increasing pore pressure. 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 magnitude and its correlation to the injection rate and volume of fluid present in the reservoir.