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 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

Comprehensive study of micro-seismicity by using an automatic monitoring platform

A modern digital seismic network, with many stations optimally distributed on the earthquake causative seismic zone, enables detection of very low magnitude earthquakes and determination of their source parameters. It is essential to associate to such kind of networks procedures to analyze the huge amount of continuously recorded data for monitoring the space-time-magnitude evolution of natural and/or induced seismicity. Hence, the demand for near-real-time, automated data collection and analysis procedures for assisting seismic network operators in carrying out microearthquake monitoring is growing. In response to this need, we designed a computational software platform, TREMOR, for fast and reliable detection and characterization of seismicity recorded by a dense local seismic network. TREMOR integrates different open-source seismological algorithms for earthquake signal detection, location, and source characterizations in a fully automatic workflow. We applied the platform in play-back mode to the continuous waveform data recorded during 1 month at the Japanese Hi-net seismic network in the Nagano region (Japan) and compared the resulting catalog with the Japan Meteorological Agency bulletin in terms of number of detections, location pattern and magnitudes. The results show that the completeness magnitude of the new seismic catalog decreased by 0.35 units of the local magnitude scale and consequently the number of events increased by about 60% with respect to the available catalog. Moreover, the fault plane solutions resulted coherent with the stress regime of the region, and the Vp/Vs ratio well delineated the main structural features of the area. According to our results, TREMOR has shown to be a valid tool for investigating and studying earthquakes, especially to identify and monitor natural or induced micro-seismicity

4D travel-time tomography as a tool for tracking fluid-driven medium changes in offshore oil-gas exploitation areas

The monitoring of rock volume where offshore exploitation activities take place is crucial to assess the corresponding seismic hazard. Fluid injection/extraction operations generate a pore fluid pressure perturbation into the volume hosting the reservoir which, in turn, may trigger new failures and induce changes in the elastic properties of rocks. Our purpose is to evaluate the feasibility of reconstructing pore pressure perturbation diffusion in the host medium by imaging the 4D velocity changes using active seismic. We simulated repeated active offshore surveys and imaged the target volume. We constructed the velocity model perturbed by the fluid injection using physical modeling and evaluated under which conditions the repeated surveys could image the velocity changes. We found that the induced pressure perturbation causes seismic velocity variations ranging between 2-5% and 15-20%, depending on the different injection conditions and medium properties. So, in most cases, time-lapse tomography is very efficient in tracking the perturbation. The noise level characterizing the recording station sites is a crucial parameter. Since we evaluated the feasibility of the proposed 4D imaging strategy under different realistic environmental and operational conditions, our results can be directly applied to set up and configure the acquisition layout of surveys aimed at retrieving fluid-induced medium changes in the hosting medium. Moreover, our results can be considered as a useful starting point to design the guidelines to monitor exploitation areas

A Bayesian Method for Real-Time Earthquake Location Using Multiparameter Data

A primary task of a network-based, earthquake early warning system is the prompt event detection and location, needed to assess the magnitude of the event and its potential damage through the predicted peak ground shaking amplitude using empirical attenuation relationships. Most of real-time, automatic earthquake location methods ground on the progressive measurement of the first P-wave arrival time at stations located at increasing distances from the source but recent approaches showed the feasibility to improve the accuracy and rapidity of the earthquake location by using the additional information carried by the P-wave polarization or amplitude, especially with unfavorable seismic network lay-outs. Here, we propose an evolutionary, Bayesian method for the real-time earthquake location which combines the information derived from the differential P-wave arrival times, amplitude ratios and back-azimuths measured at a minimum of two stations. As more distant stations record the P-wave, the posterior pdf is updated and new earthquake location parameters are determined along with their uncertainty. To validate the location method, we performed a retrospective analysis of mainshocks (M > 4.5) occurred during the 2016–2017 Central Italy earthquake sequence by simulating the typical acquisition layouts of in-land, coastal and linear array of stations. Results show that with the combined use of the three parameters, 2–4 s after the first P-wave detection, the method converges to stable and accurate determinations of epicentral coordinates and depth even with a nonoptimal coverage of stations. The proposed methodology can be generalized and adapted to the off-line analysis of seismic records collected by standard local networks

Decade-long monitoring of seismic velocity changes at the Irpinia fault system (southern Italy) reveals pore pressure pulsations

Several studies have monitored crustal seismic velocity changes and attempted to relate them to the stress state and physical properties in volume embedding fault systems. The aim is to provide constraints on fault system dynamics and earthquake triggering mechanisms. Here, we reconstruct the spatiotemporal (4D) seismic velocity images of volume embedding the Irpinia fault system (IFS, South Italy), which originated the 1980 Ms 6.9 multi-segmented ruptures. By inverting data from more than ten years of continuous seismicity monitoring, we retrieved time-constant velocity anomalies, whose shapes correlate well with crustal lithology, while time-changing (up to 20%) velocity anomalies are mapped in the central region. Here, the Vp-to-Vs changes at depths of 1–5 km and 8–12 km correlate well with groundwater recharge and geodetic displacement during the same time interval. This correlation provides evidence for the existence of pulsating, pore pressure changes induced by groundwater recharge processes in a deep volume (8–12 km of depth), fractured and saturated with a predominant gas phase (likely CO2). We suggest that tomographic measurements of the Vp-to-Vs spatiotemporal changes are a suitable proxy to track the pore pressure evolution at depth in highly sensitive regions of fault systems