Real-time monitoring of seismic moment and radiated energy

Author(s):  
Davide Scafidi ◽  
Daniele Spallarossa ◽  
Matteo Picozzi ◽  
Dino Bindi

<p>Understanding the dynamics of faulting is a crucial target in earthquake source physics (Yoo et al., 2010). To study earthquake dynamics it is indeed necessary to look at the source complexity from different perspectives; in this regard, useful information is provided by the seismic moment (M0), which is a static measure of the earthquake size, and the seismic radiated energy (ER), which is connected to the rupture kinematics and dynamics (e.g. Bormann & Di Giacomo 2011a). Studying spatial and temporal evolution of scaling relations between scaled energy (i.e., e = ER/M0) versus the static measure of source dimension (M0) can provide valuable indications for understanding the earthquake generation processes, single out precursors of stress concentrations, foreshocks and the nucleation of large earthquakes (Picozzi et al., 2019). In the last ten years, seismology has undergone a terrific development. Evolution in data telemetry opened the new research field of real-time seismology (Kanamori 2005), which targets are the rapid determination of earthquake location and size, the timely implementation of emergency plans and, under favourable conditions, earthquake early warning. On the other hand, the availability of denser and high quality seismic networks deployed near faults made possible to observe very large numbers of micro-to-small earthquakes, which is pushing the seismological community to look for novel big data analysis strategies. Large earthquakes in Italy have the peculiar characteristic of being followed within seconds to months by large aftershocks of magnitude similar to the initial quake or even larger, demonstrating the complexity of the Apennines’ faults system (Gentili and Giovanbattista, 2017). Picozzi et al. (2017) estimated the radiated seismic energy and seismic moment from P-wave signals for almost forty earthquakes with the largest magnitude of the 2016-2017 Central Italy seismic sequence. Focusing on S-wave signals recorded by local networks, Bindi et al. (2018) analysed more than 1400 earthquakes in the magnitude ranges 2.5 ≤ Mw ≤ 6.5 of the same region occurred from 2008 to 2017 and estimated both ER and M0, from which were derived the energy magnitude (Me) and Mw for investigating the impact of different magnitude scales on the aleatory variability associated with ground motion prediction equations. In this work, exploiting first steps made in this direction by Picozzi et al. (2017) and Bindi et al. (2018), we derived a novel approach for the real-time, robust estimation of seismic moment and radiated energy of small to large magnitude earthquakes recorded at local scales. In the first part of the work, we describe the procedure for extracting from the S-wave signals robust estimates of the peak displacement (PDS) and the cumulative squared velocity (IV2S). Then, exploiting a calibration data set of about 6000 earthquakes for which well-constrained M0 and theoretical ER values were available, we describe the calibration of empirical attenuation models. The coefficients and parameters obtained by calibration were then used for determining ER and M0 of a testing dataset</p>

Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 755-762 ◽  
Author(s):  
Arild Buland ◽  
Martin Landrø

The impact of prestack time migration on porosity estimation has been tested on a 2-D seismic line from the Valhall/Hod area in the North Sea. Porosity is estimated in the Cretaceous chalk section in a two‐step procedure. First, P-wave and S-wave velocity and density are estimated by amplitude variation with offset (AVO) inversion. These parameters are then linked to porosity through a petrophysical rock data base based on core plug analysis. The porosity is estimated both from unmigrated and prestack migrated seismic data. For the migrated data set, a standard prestack Kirchhoff time migration is used, followed by simple angle and amplitude corrections. Compared to modern high‐cost, true amplitude migration methods, this approach is faster and more practical. The test line is structurally fairly simple, with a maximum dip of 5°; but the results differ significantly, depending on whether migration is applied prior to the inversion. The maximum difference in estimated porosity is of the order of 10% (about 50% relative change). High‐porosity zones estimated from the unmigrated data were not present on the porosity section estimated from the migrated data.


Author(s):  
Daniele Spallarossa ◽  
Matteo Picozzi ◽  
Davide Scafidi ◽  
Paola Morasca ◽  
Chiara Turino ◽  
...  

Abstract We present Rapid Assessment of MOmeNt and Energy Service (RAMONES), a service for disseminating through a web interface, the estimates of seismic moment (M0) and radiated energy (ER) for earthquakes occurring in central Italy with local magnitudes above 1.7. The service is based on a fully-automatic procedure developed for downloading and processing open seismological data from the European Integrated Data Archive, Italian Civil Protection repository, and Incorporated Research Institutions for Seismology (IRIS). In its actual configuration, RAMONES uses the seismic catalog generated through the event webservice of the Italian Institute of Geophysics and Volcanology (compliant with International Federation of Digital Seismograph Networks standards) to guide the data download. The concept of RAMONES is to estimate M0 and ER from features extracted directly from recordings, namely the S-wave peak displacement (PDS) and the integral of the squared velocity (IV2S) evaluated over the S-wave window at local distances. A data set composed of 6515 earthquakes recorded in central Italy between 2008 and 2018 was used to calibrate the attenuation models relating M0 to PDS and ER to IV2S, including station corrections. The calibration values for M0 and ER were extracted from the source spectra obtained by applying a decomposition approach to the Fourier amplitude spectra known as the generalized inversion technique. To test the capabilities of RAMONES, we validate the attenuation models by performing residual analysis over about 60 earthquakes occurring in 2019 that were used for the spectral decomposition analysis but not considered in the calibration phase. Since January 2020, a testing operational phase has been running, and RAMONES has analyzed about 800 earthquakes by September 2020. The distribution of the source parameters and their relevant scaling relationships are automatically computed and disseminated in the form of maps, parametric tables, figures, and reports available through the RAMONES web interface.


Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. KS31-KS41 ◽  
Author(s):  
Wenjie Jiao ◽  
Michael Davidson ◽  
Arcangelo Sena ◽  
Bradley L. Bankhead ◽  
Yu Xia ◽  
...  

We investigated the method of estimating seismic moment and moment magnitude for microseismic events. We determined that the [Formula: see text] defined by Bowers and Hudson is the proper scalar moment to be used in microseismic studies for characterizing the size of an event and calculating its moment magnitude. For non-double-couple sources, the proportional relationship between body-wave amplitude and seismic moment in the Brune model breaks down. So under such situations, the Brune model is not an appropriate way to estimate the seismic moment and magnitude. Moreover, the S-wave alone is not sufficient for determining the total seismic moment. Instead, the P-wave must be analyzed. An example Barnett Shale data set was studied, and the results concluded that the magnitudes estimated with the Brune model could be off by as much as 1.92, with an absolute average of 0.35. The moment magnitudes based on the scalar moment [Formula: see text] also gave a significantly different event size distribution and b-value estimation. Finally, attenuation also played a role in estimating the moment magnitude. With a typical average attenuation factor of [Formula: see text], the average magnitude correction for our field data set was on the order of 0.15. However, it could reach 0.3 for events far away from the monitoring well.


Geophysics ◽  
2000 ◽  
Vol 65 (5) ◽  
pp. 1446-1454 ◽  
Author(s):  
Side Jin ◽  
G. Cambois ◽  
C. Vuillermoz

S-wave velocity and density information is crucial for hydrocarbon detection, because they help in the discrimination of pore filling fluids. Unfortunately, these two parameters cannot be accurately resolved from conventional P-wave marine data. Recent developments in ocean‐bottom seismic (OBS) technology make it possible to acquire high quality S-wave data in marine environments. The use of (S)-waves for amplitude variation with offset (AVO) analysis can give better estimates of S-wave velocity and density contrasts. Like P-wave AVO, S-wave AVO is sensitive to various types of noise. We investigate numerically and analytically the sensitivity of AVO inversion to random noise and errors in angles of incidence. Synthetic examples show that random noise and angle errors can strongly bias the parameter estimation. The use of singular value decomposition offers a simple stabilization scheme to solve for the elastic parameters. The AVO inversion is applied to an OBS data set from the North Sea. Special prestack processing techniques are required for the success of S-wave AVO inversion. The derived S-wave velocity and density contrasts help in detecting the fluid contacts and delineating the extent of the reservoir sand.


Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. B183-B195 ◽  
Author(s):  
K. De Meersman ◽  
J.-M. Kendall ◽  
M. van der Baan

We relocate 303 microseismic events recorded in 1998 by sensors in a single borehole in the North Sea Valhall oil field. A semiautomated array analysis method repicks the P- and S-wave arrival times and P-wave polarizations, which are needed to locate these events. The relocated sources are confined predominantly to a [Formula: see text]-thick zone just above the reservoir, and location uncertainties are half those of previous efforts. Multiplet analysis identifies 40 multiplet groups, which include 208 of the 303 events. The largest group contains 24 events, and five groups contain 10 or more events. Within each multiplet group, we further improve arrival-time picking through crosscorrelation, which enhances the relative accuracy of the relocated events and reveals that more than 99% of the seismic activity lies spatially in three distinct clusters. The spatial distribution of events and wave-form similarities reveal two faultlike structures that match well with north-northwest–south-southeast-trending fault planes interpreted from 3D surface seismic data. Most waveform differences between multiplet groups located on these faults can be attributed to S-wave phase content and polarity or P-to-S amplitude ratio. The range in P-to-S amplitude ratios observed on the faults is explained best in terms of varying source mechanisms. We also find a correlation between multiplet groups and temporal variations in seismic anisotropy, as revealed by S-wave splitting analysis. We explain these findings in the context of a cyclic recharge and dissipation of cap-rock stresses in response to production-driven compaction of the underlying oil reservoir. The cyclic nature of this mechanism drives the short-term variations in seismic anisotropy and the reactivation of microseismic source mechanisms over time.


Author(s):  
Manudul Pahansen de Alwis ◽  
Karl Garme

The stochastic environmental conditions together with craft design and operational characteristics make it difficult to predict the vibration environments aboard high-performance marine craft, particularly the risk of impact acceleration events and the shock component of the exposure often being associated with structural failure and human injuries. The different timescales and the magnitudes involved complicate the real-time analysis of vibration and shock conditions aboard these craft. The article introduces a new measure, severity index, indicating the risk of severe impact acceleration, and proposes a method for real-time feedback on the severity of impact exposure together with accumulated vibration exposure. The method analyzes the immediate 60 s of vibration exposure history and computes the severity of impact exposure as for the present state based on severity index. The severity index probes the characteristic of the present acceleration stochastic process, that is, the risk of an upcoming heavy impact, and serves as an alert to the crew. The accumulated vibration exposure, important for mapping and logging the crew exposure, is determined by the ISO 2631:1997 vibration dose value. The severity due to the impact and accumulated vibration exposure is communicated to the crew every second as a color-coded indicator: green, yellow and red, representing low, medium and high, based on defined impact and dose limits. The severity index and feedback method are developed and validated by a data set of 27 three-hour simulations of a planning craft in irregular waves and verified for its feasibility in real-world applications by full-scale acceleration data recorded aboard high-speed planing craft in operation.


2015 ◽  
Vol 3 (1) ◽  
pp. SF43-SF54 ◽  
Author(s):  
Shelby L. Peterie ◽  
Richard D. Miller

Tunnel locations are accurately interpreted from diffraction sections of focused mode converted P- to S-wave diffractions from a perpendicular tunnel and P-wave diffractions from a nonperpendicular (oblique) tunnel. Near-surface tunnels are ideal candidates for diffraction imaging due to their small size relative to the seismic wavelength and large acoustic impedance contrast at the tunnel interface. Diffraction imaging algorithms generally assume that the velocities of the primary wave and the diffracted wave are approximately equal, and that the diffraction apex is recorded directly above the scatterpoint. Scattering phenomena from shallow tunnels with kinematic properties that violate these assumptions were observed in one field data set and one synthetic data set. We developed the traveltime equations for mode-converted and oblique diffractions and demonstrated a diffraction imaging algorithm designed for the roll-along style of acquisition. Potential processing and interpretation pitfalls specific to these diffraction types were identified. Based on our observations, recommendations were made to recognize and image mode-converted and oblique diffractions and accurately interpret tunnel depth, horizontal location, and azimuth with respect to the seismic line.


Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1172-1180 ◽  
Author(s):  
W. Scott Leaney ◽  
Colin M. Sayers ◽  
Douglas E. Miller

Multioffset vertical seismic profile (VSP) experiments, commonly referred to as walkaways, enable anisotropy to be measured reliably in the field. The results can be fed into modeling programs to study the impact of anisotropy on velocity analysis, migration, and amplitude versus offset (AVO). Properly designed multioffset VSPs can also provide the target AVO response measured under optimum conditions, since the wavelet is recorded just above the reflectors of interest with minimal reflection point dispersal. In this paper, the multioffset VSP technique is extended to include multioffset azimuths, and a multiazimuthal multiple VSP data set acquired over a carbonate reservoir is analyzed for P-wave anisotropy and AVO. Direct arrival times down to the overlying shale and reflection times and amplitudes from the carbonate are analyzed. Data analysis involves a three‐term fit to account for nonhyperbolic moveout, dip, and azimuthal anisotropy. Results indicate that the overlying shale is transversely isotropic with a vertical axis of symmetry (VTI), while the carbonate shows 4–5% azimuthal anisotropy in traveltimes. The fast direction is consistent with the maximum horizontal stress orientation determined from break‐out logs and is also consistent with the strike of major faults. AVO analysis of the reflection from the top of the carbonate layer shows a critical angle reduction in the fast direction and maximum gradient in the slow direction. This agrees with modeling and indicates a greater amplitude sensitivity in the slow direction—the direction perpendicular to fracture strike. In principle, 3-D surveys should have wide azimuthal coverage to characterize fractured reservoirs. If this is not possible, it is important to have azimuthal line coverage in the minimum horizontal stress direction to optimize the use of AVO for fractured reservoir characterization. This direction can be obtained from multiazimuthal walkaways using the azimuthal P-wave analysis techniques presented.


2021 ◽  
Author(s):  
Fatma Sevil Malcıoğlu ◽  
Hakan Süleyman ◽  
Eser Çaktı

Abstract An MW 4.5 earthquake took place on September 24, 2019 in the Marmara Sea. Two days after, on September 26, 2019, Marmara region was rattled by an MW5.7 earthquake. With the intention of compiling an ample strong ground motion data set of recordings, we have utilized the stations of Istanbul Earthquake Rapid Response and Early Warning System operated by the Department of Earthquake Engineering of Boğaziçi University and of the National Strong Motion Network operated by AFAD. All together 438 individual records are used to calculate the source parameters of events; namely, corner frequency, radius, rupture area, average source dislocation, source duration and stress drop. Some of these parameters are compared with empirical relationships and discussed extensively. Duration characteristics are analyzed in two steps; first, by making use of the time difference between P-wave and S-wave onsets and then, by considering S-wave durations and significant durations. It is observed that they yield similar trends with global models. PGA, PGV and SA values are compared with three commonly used ground motion prediction models. At distances closer than about 60 km observed intensity measures mostly conform with the GMPE predictions. Beyond 60 km their attenuation is clearly faster than those of GMPEs. Frequency-dependent Q models are developed for both events. Their consistency with existing regional models are confirmed.


2019 ◽  
Vol 218 (3) ◽  
pp. 1873-1891 ◽  
Author(s):  
Farbod Khosro Anjom ◽  
Daniela Teodor ◽  
Cesare Comina ◽  
Romain Brossier ◽  
Jean Virieux ◽  
...  

SUMMARY The analysis of surface wave dispersion curves (DCs) is widely used for near-surface S-wave velocity (VS) reconstruction. However, a comprehensive characterization of the near-surface requires also the estimation of P-wave velocity (VP). We focus on the estimation of both VS and VP models from surface waves using a direct data transform approach. We estimate a relationship between the wavelength of the fundamental mode of surface waves and the investigation depth and we use it to directly transform the DCs into VS and VP models in laterally varying sites. We apply the workflow to a real data set acquired on a known test site. The accuracy of such reconstruction is validated by a waveform comparison between field data and synthetic data obtained by performing elastic numerical simulations on the estimated VP and VS models. The uncertainties on the estimated velocity models are also computed.


Sign in / Sign up

Export Citation Format

Share Document