Time Domain Source Parameter Estimation of Natural and Man-Induced Microearthquakes at the Geysers Geothermal Field

Water injection in geothermal areas is the preferential strategy to sustain the natural production of geothermal resources. In this context, monitoring microearthquakes is a fundamental tool to track changes in the reservoirs in terms of soil composition, response to injections, and resource exploit...

Full description

Bibliographic Details
Main Authors: Valeria Longobardi, Sahar Nazeri, Simona Colombelli, Raffaele Rea, Grazia De Landro, Aldo Zollo
Format: Article
Language:English
Published: MDPI AG 2023-01-01
Series:Energies
Subjects:
Online Access:https://www.mdpi.com/1996-1073/16/3/1121
_version_ 1797624726397386752
author Valeria Longobardi
Sahar Nazeri
Simona Colombelli
Raffaele Rea
Grazia De Landro
Aldo Zollo
author_facet Valeria Longobardi
Sahar Nazeri
Simona Colombelli
Raffaele Rea
Grazia De Landro
Aldo Zollo
author_sort Valeria Longobardi
collection DOAJ
description Water injection in geothermal areas is the preferential strategy to sustain the natural production of geothermal resources. In this context, monitoring microearthquakes is a fundamental tool to track changes in the reservoirs in terms of soil composition, response to injections, and resource exploitation with space and time. Therefore, refined source characterization is crucial to better estimate the size, source mechanism, and rupture process of microearthquakes, as they are possibly related to industrial activities, and to identify any potential variation in the background seismicity. Standard approaches for source parameter estimation are ordinarily based on the modelling of Fourier displacement spectra and its characteristic parameters: the low-frequency spectral level and corner frequency. Here, we apply an innovative time domain technique that uses the curves of P-wave amplitude vs. time along the seismogram. This methodology allows estimation of seismic moment, source radius, and stress release from the plateau level and the corner time of the average logarithm of P-wave displacement versus time with the assumption of a triangular moment rate function, uniform rupture speed, and a constant/frequency-independent Q-factor. In the current paper, this time domain methodology is implemented on a selected catalog of microearthquakes consisting of 83 events with a moment magnitude ranging between 1.0 and 1.5 that occurred during a 7-year period (2007–2014) of fluid extraction/injection around Prati-9 and Prati-29 wells at The Geysers geothermal field. The results show that the time domain technique provides accurate seismic moment (moment magnitude) and rupture duration/radius estimates of microearthquakes down to the explored limit (M 1) while accounting for the anelastic attenuation effect in the radiated high-frequency wavefield. The retrieved source radius vs. moment scaling is consistent with a self-similar, constant stress drop scaling model, which proves an appropriate attenuation correction and the validity of the assumed, triangular moment rate function for microearthquake ruptures. Two alternative mechanical models are proposed to explain the observed difference (about two orders of magnitude) in the retrieved average stress release estimates between the time and frequency domain methods. We argue that the two quantities may not refer to the same physical quantity representing the stress release of earthquake ruptures. Either the smaller stress release values from the time domain method may indicate a larger fracture area (by a factor of 20) radiating the observed P-waveforms than the one estimated from the corner frequencies, or the frequency domain estimate is a proxy for dynamic stress release while the time domain is more representative of the static release. The latter is associated with a much lower dynamic friction value than static friction value at the fault during the rupture process.
first_indexed 2024-03-11T09:46:37Z
format Article
id doaj.art-0a2dbdf7efee44ad8988736b8ade5990
institution Directory Open Access Journal
issn 1996-1073
language English
last_indexed 2024-03-11T09:46:37Z
publishDate 2023-01-01
publisher MDPI AG
record_format Article
series Energies
spelling doaj.art-0a2dbdf7efee44ad8988736b8ade59902023-11-16T16:32:56ZengMDPI AGEnergies1996-10732023-01-01163112110.3390/en16031121Time Domain Source Parameter Estimation of Natural and Man-Induced Microearthquakes at the Geysers Geothermal FieldValeria Longobardi0Sahar Nazeri1Simona Colombelli2Raffaele Rea3Grazia De Landro4Aldo Zollo5Department of Physics “E. Pancini”, University of Naples “Federico II”, 80126 Naples, ItalyDepartment of Physics “E. Pancini”, University of Naples “Federico II”, 80126 Naples, ItalyDepartment of Physics “E. Pancini”, University of Naples “Federico II”, 80126 Naples, ItalyDepartment of Physics “E. Pancini”, University of Naples “Federico II”, 80126 Naples, ItalyDepartment of Physics “E. Pancini”, University of Naples “Federico II”, 80126 Naples, ItalyDepartment of Physics “E. Pancini”, University of Naples “Federico II”, 80126 Naples, ItalyWater injection in geothermal areas is the preferential strategy to sustain the natural production of geothermal resources. In this context, monitoring microearthquakes is a fundamental tool to track changes in the reservoirs in terms of soil composition, response to injections, and resource exploitation with space and time. Therefore, refined source characterization is crucial to better estimate the size, source mechanism, and rupture process of microearthquakes, as they are possibly related to industrial activities, and to identify any potential variation in the background seismicity. Standard approaches for source parameter estimation are ordinarily based on the modelling of Fourier displacement spectra and its characteristic parameters: the low-frequency spectral level and corner frequency. Here, we apply an innovative time domain technique that uses the curves of P-wave amplitude vs. time along the seismogram. This methodology allows estimation of seismic moment, source radius, and stress release from the plateau level and the corner time of the average logarithm of P-wave displacement versus time with the assumption of a triangular moment rate function, uniform rupture speed, and a constant/frequency-independent Q-factor. In the current paper, this time domain methodology is implemented on a selected catalog of microearthquakes consisting of 83 events with a moment magnitude ranging between 1.0 and 1.5 that occurred during a 7-year period (2007–2014) of fluid extraction/injection around Prati-9 and Prati-29 wells at The Geysers geothermal field. The results show that the time domain technique provides accurate seismic moment (moment magnitude) and rupture duration/radius estimates of microearthquakes down to the explored limit (M 1) while accounting for the anelastic attenuation effect in the radiated high-frequency wavefield. The retrieved source radius vs. moment scaling is consistent with a self-similar, constant stress drop scaling model, which proves an appropriate attenuation correction and the validity of the assumed, triangular moment rate function for microearthquake ruptures. Two alternative mechanical models are proposed to explain the observed difference (about two orders of magnitude) in the retrieved average stress release estimates between the time and frequency domain methods. We argue that the two quantities may not refer to the same physical quantity representing the stress release of earthquake ruptures. Either the smaller stress release values from the time domain method may indicate a larger fracture area (by a factor of 20) radiating the observed P-waveforms than the one estimated from the corner frequencies, or the frequency domain estimate is a proxy for dynamic stress release while the time domain is more representative of the static release. The latter is associated with a much lower dynamic friction value than static friction value at the fault during the rupture process.https://www.mdpi.com/1996-1073/16/3/1121earthquake source parameterstime domain techniqueinduced seismicityThe Geysers geothermal field
spellingShingle Valeria Longobardi
Sahar Nazeri
Simona Colombelli
Raffaele Rea
Grazia De Landro
Aldo Zollo
Time Domain Source Parameter Estimation of Natural and Man-Induced Microearthquakes at the Geysers Geothermal Field
Energies
earthquake source parameters
time domain technique
induced seismicity
The Geysers geothermal field
title Time Domain Source Parameter Estimation of Natural and Man-Induced Microearthquakes at the Geysers Geothermal Field
title_full Time Domain Source Parameter Estimation of Natural and Man-Induced Microearthquakes at the Geysers Geothermal Field
title_fullStr Time Domain Source Parameter Estimation of Natural and Man-Induced Microearthquakes at the Geysers Geothermal Field
title_full_unstemmed Time Domain Source Parameter Estimation of Natural and Man-Induced Microearthquakes at the Geysers Geothermal Field
title_short Time Domain Source Parameter Estimation of Natural and Man-Induced Microearthquakes at the Geysers Geothermal Field
title_sort time domain source parameter estimation of natural and man induced microearthquakes at the geysers geothermal field
topic earthquake source parameters
time domain technique
induced seismicity
The Geysers geothermal field
url https://www.mdpi.com/1996-1073/16/3/1121
work_keys_str_mv AT valerialongobardi timedomainsourceparameterestimationofnaturalandmaninducedmicroearthquakesatthegeysersgeothermalfield
AT saharnazeri timedomainsourceparameterestimationofnaturalandmaninducedmicroearthquakesatthegeysersgeothermalfield
AT simonacolombelli timedomainsourceparameterestimationofnaturalandmaninducedmicroearthquakesatthegeysersgeothermalfield
AT raffaelerea timedomainsourceparameterestimationofnaturalandmaninducedmicroearthquakesatthegeysersgeothermalfield
AT graziadelandro timedomainsourceparameterestimationofnaturalandmaninducedmicroearthquakesatthegeysersgeothermalfield
AT aldozollo timedomainsourceparameterestimationofnaturalandmaninducedmicroearthquakesatthegeysersgeothermalfield