Spatial and temporal multiplet analysis for identification of dominant fluid migration path at The Geysers geothermal field, California

Multiplet analysis is based on the identification of seismic events with very similar waveforms which are used then to enhance seismological analysis e.g. by precise relocation of sources. In underground fluid injection conditions, it is a tool frequently used for imaging of subsurface fracture system. We identify over 150 repeatedly activated seismic sources within seismicity cluster induced by fluid injection in NW part of The Geysers geothermal field (California). Majority of multiple events (ME) occur along N–S oriented planar structure which we interpret as a fault plane. Remaining ME are distributed along structures interpreted as fractures, forming together a system of interconnected cracks enabling fluid migration. Temporal analysis reveals that during periods of relatively low fluid injection the proportion of ME to non-multiple events is higher than during periods of high injection. Moreover, ME which occur within the fault differ in activity rate and source properties from ME designating the fractures and non-multiple events. In this study we utilize observed differences between ME occurring within various structures and non-multiple events to describe hydraulic conditions within the reservoir. We show that spatial and temporal analysis of multiplets can be used for identification and characterization of dominant fluid migration paths.

Long Memory in Earthquake Time Series: The Case Study of the Geysers Geothermal Field

The present study aims at proving the existence of long memory (or long-range dependence) in the earthquake process through the analysis of time series of induced seismicity. Specifically, we apply alternative statistical techniques borrowed from econometrics to the seismic catalog of The Geysers geothermal field (California), the world’s largest geothermal field. The choice of the study area is essentially guided by the completeness of the seismic catalog at smaller magnitudes (a drawback of conventional catalogs of natural seismicity). Contrary to previous studies, where the long-memory property was examined by using non-parametric approaches (e.g., rescaled range analysis), we assume a fractional integration model for which the degree of memory is defined by a real parameter d, which is related to the best known Hurst exponent. In particular, long-memory behavior is observed for d > 0. We estimate and test the value of d (i.e., the hypothesis of long memory) by applying parametric, semi-parametric, and non-parametric approaches to time series describing the daily number of earthquakes and the logarithm of the (total) seismic moment released per day. Attention is also paid to examining the sensitivity of the results to the uncertainty in the completeness magnitude of the catalog, and to investigating to what extent temporal fluctuations in seismic activity induced by injection operations affect the value of d. Temporal variations in the values of d are analyzed together with those of the b-value of the Gutenberg and Richter law. Our results indicate strong evidence of long memory, with d mostly constrained between 0 and 0.5. We observe that the value of d tends to decrease with increasing the magnitude completeness threshold, and therefore appears to be influenced by the number of information in the chain of intervening related events. Moreover, we find a moderate but significant negative correlation between d and the b-value. A negative, albeit weaker correlation is found between d and the fluid injection, as well as between d and the annual number of earthquakes.

An application of induced event interferometry approach at The Geysers Geothermal Field, California, USA

<p>While the classical tomography approaches, e.g., P-, S-, and/or surface-wave traveltime tomography, provide a general structure of the Earth’s interior, new developments in signal processing of interferometry approaches are needed to obtain a high-resolution velocity structure. If the number of earthquakes is adequate, the virtual seismometer method may be a solution in regions with sparse instrumental coverage. Theoretically, the empirical Green’s functions between a pair of events can be retrieved using earthquake’s cross-correlations. Here, an event interferometry approach was used on a very small scale around Prati-9 and Prati-29 injection wells in the NW of The Geysers Geothermal Field. The study region experienced intense injection-induced seismicity. We selected all events with location uncertainties less than 50 m in a cuboid of the horizontal side ~1 × ~2 km and the vertical edge at depths between 1.0 and 2.0 km. The cuboid was cut into 100m thick layers, and we applied to events from each layer criteria enabling a quasi 2D approach. After calculating the Rayleigh wave group velocity dispersion curves, further processing was performed at a 0.2s period, selected based on the sensitivity kernel criterion. Finally, the relative velocity model of each layer at the depth z was obtained by subtracting the velocity model of the just overlying layer (at the depth z-100m) from the model of this layer. Our resultant velocity model in the study area indicated four low-velocity anomalies. The first one can be linked by the two layers interface topography variation at the top of the cuboid (depth 1000 m). The secondary faults can cause the second low-velocity anomaly. The other two anomalies look to result from fluid injection into Prati-9 and Prati-29 wells. <br>This work was supported under the S4CE: "Science for Clean Energy" project, which has received funding from the European Union’s Horizon 2020 research and innovation program, under grant agreement No 764810.</p>

High injection rates decrease the probability of creating undesired, far-reaching fluid migration pathways at Cooper Basin geothermal field

<p>Enhanced Geothermal Systems apply the pressurized fluid injection to fracture impermeable rocks to form pathways in which water circulates. The cold water under high pressure is pumped into the hot subsoil, where it heats up and returns to the surface. However, the induced fractures may coalesce into unwanted paths that allow the fluids to reach pre-existing faults, triggering major seismic events.</p><p>This work investigates the relationship between injection and a degree of disordering of sources, ZZ, at Cooper Basin geothermal field in Australia, following the methodology developed and applied to study The Geysers geothermal field case (Lasocki & Orlecka-Sikora, 2020). The parameter ZZ quantifies the potential of seismicity to build pathways for fluid migration. It is the average distance between the seismic events in the eight-dimensional parameter space consisting of three hypocentral coordinates, T- and P-axis plunges, T-axis trend, and polar and azimuthal angles in the spherical system of coordinates beginning at the open hole of an injection well. A decrease of ZZ indicates an increasing hazard of forming far-reaching migration pathways. In The Geysers case, ZZ turned out to be highly correlated with the injection rate.</p><p>Here we focus on the case of Habanero 4 well stimulation from 17 - November 30, 2012 (data access, see: IS EPOS, 2020). We processed 489 seismic events with known focal mechanisms. The events moment magnitude varies between 0.8 and 3.1.  </p><p>Our analysis shows that ZZ is significantly correlated with both the injection rate and the wellhead pressure. The higher the injection rate / the wellhead pressure was, the less probable was the creation of undesired fluid migration pathways. The Cooper Basin’s and The Geyser’s reservoir rocks are vastly different, the former – granite, the latter – greywacke sandstone, likewise the stimulation techniques applied in these two reservoirs. However, in both cases, ZZ was positively correlated with injection rate; thus, the potential to build unwanted paths for fluids was negatively correlated. These results suggest that such correlation may be a global feature of rock fracturing caused by pressurized fluid injections.</p><p><em>This work has been supported by S4CE (Science for Clean Energy) project, funded from the European Union’s Horizon 2020 - Framework Programme, under grant agreement No 764810 and by PRIN-MATISSE (20177EPPN2) project funded by Italian Ministry of Education and Research.</em></p><p> </p><p><strong>References:</strong></p><p>IS EPOS (2020), Episode: COOPER BASIN, https://tcs.ah-epos.eu/#episode:COOPER_BASIN, doi:10.25171/InstGeoph_PAS_ISEPOS-2020-001</p><p>Lasocki, S., & Orlecka-Sikora, B. (2020). High injection rates counteract formation of far-reaching fluid migration pathways at The Geysers geothermal field. Geophysical Research Letters, 47, e2019GL086212. https://doi.org/10.1029/2019GL086212</p>

Double-difference seismic attenuation tomography method and its application to The Geysers geothermal field, California

SUMMARY Knowledge of attenuation structure is important for understanding subsurface material properties. We have developed a double-difference seismic attenuation (DDQ) tomography method for high-resolution imaging of 3-D attenuation structure. Our method includes two main elements, the inversion of event-pair differential ${t^*}$ ($d{t^*}$) data and 3-D attenuation tomography with the $d{t^*}$ data. We developed a new spectral ratio method that jointly inverts spectral ratio data from pairs of events observed at a common set of stations to determine the $d{t^*}$ data. The spectral ratio method cancels out instrument and site response terms, resulting in more accurate $d{t^*}$ data compared to absolute ${t^*}$ from traditional methods using individual spectra. Synthetic tests show that the inversion of $d{t^*}$ data using our spectral ratio method is robust to the choice of source model and a moderate degree of noise. We modified an existing velocity tomography code so that it can invert $d{t^*}$ data for 3-D attenuation structure. We applied the new method to The Geyser geothermal field, California, which has vapour-dominated reservoirs and a long history of water injection. A new Qp model at The Geysers is determined using P-wave data of earthquakes in 2011, using our updated earthquake locations and Vp model. By taking advantage of more accurate $d{t^*}$ data and the cancellation of model uncertainties along the common paths outside of the source region, the DDQ tomography method achieves higher resolution, especially in the earthquake source regions, compared to the standard tomography method using ${t^*}$ data. This is validated by both the real and synthetic data tests. Our Qp and Vp models show consistent variations in a normal temperature reservoir that can be explained by variations in fracturing, permeability and fluid saturation and/or steam pressure. A prominent low-Qp and Vp zone associated with very active seismicity is imaged within a high temperature reservoir at depths below 2 km. This anomalous zone is likely partially saturated with injected fluids.

Using geodetic data in geothermal areas

Geodetic observations, often in conjunction with other data, provide a cost-effective means for identifying and characterizing geothermal resources. A review of the various methods reveals how the technology for measuring deformation has advanced considerably in the past few decades. Currently, interferometric synthetic aperture radar is the method of choice for monitoring deformation at a geothermal field. A discussion of geodetic monitoring at The Geysers geothermal field, California, illustrates some of the progress made and the challenges that remain.

Induced earthquake potential in geothermal reservoirs: Insights from The Geysers, California

Geothermal reservoir production and associated induced seismicity may experience pronounced attention in the near future, given the ambitious plans for reducing greenhouse gas emissions toward a carbon-neutral economy and society. At some geothermal sites, the occurrence of hazard- and risk-prone induced earthquakes caused by or associated with reservoir stimulation has resulted in project shutdown (e.g., Pohang, South Korea, and Basel Deep Heat Mining, Switzerland). At other geothermal sites, the maximum event magnitudes were successfully maintained below a threshold defined by local authorities (e.g., Helsinki St1 Deep Heat project in Helsinki, Finland). In this study, we review some of our results from seismological and geomechanical reservoir characterization at The Geysers geothermal reservoir in California, USA, the largest producing geothermal field worldwide. We relate our findings to other geothermal sites to better understand the variability of reservoir behavior. In particular, we obtain a constant and relatively low seismic injection efficiency at The Geysers, which is interpreted to be related to the large energy dissipation through thermal processes and additional dissipation through aseismic slip, the latter now being considered to play a fundamental role in earthquake nucleation. We discuss some characteristics of the seismicity from The Geysers that suggest stable reservoir seismic injection efficiency and possibly low potential to rupture into large induced earthquakes, reducing the associated seismic hazard.

Magnitude distribution complexity and variation at The Geysers geothermal field

SUMMARY Earthquake magnitude (size) distribution is a major component required for seismic hazard assessment and therefore, the accurate determination of its functional shape and variation is a task of utmost importance. Although often considered as stationary, the magnitude distribution at particular sites may significantly vary over time and space. In this study, the well-known Gutenberg–Richter (GR) law, which is widely assumed to describe earthquake magnitude distribution, is tested for a case study of seismicity induced by fluid injection at The Geysers (CA, USA) geothermal field. Statistical tests are developed and applied in order to characterize the magnitude distribution of a high quality catalogue comprising seismicity directly associated with two injection wells, at the north western part of The Geysers. The events size distribution variation is investigated with respect to spatial, temporal, fluid injection and magnitude cut-off criteria. A thorough spatio-temporal analysis is performed for defining seismicity Clusters demonstrating characteristic magnitude distributions which significantly differ from the ones of the nearby Clusters. The magnitude distributions of the entire seismic population as well as of the individual Clusters are tested for their complexity in terms of exponentiality, multimodal and multibump structure. Then, the Clusters identified are further processed and their characteristics are determined in connection to injection rate fluctuations. The results of the analysis clearly indicate that the entire magnitude distribution is definitely complex and non-exponential, whereas subsequent periods demonstrating significantly diverse magnitude distributions are identified. The regional seismicity population is divided into three major families, for one of which exponentiality of magnitude distribution is clearly rejected, whereas for the other two the GR law b-value is directly proportional to fluid injection. In addition, the b-values of these Families seem to be significantly magnitude dependent, a fact that is of major importance for seismic hazard assessment implementations. To conclude, it is strongly suggested that magnitude exponentiality must be tested before proceeding to any b-value calculations, particularly in anthropogenic seismicity cases where complex and time changeable processes take place.