scholarly journals Ground water-induced changes in velocities of P and S waves (Vp and Vs) measured using an accurately controlled seismic source

2021 ◽  
Vol 73 (1) ◽  
Author(s):  
Rina Suzuki ◽  
Koshun Yamaoka ◽  
Shuhei Tsuji ◽  
Toshiki Watanabe

AbstractWe analyze temporal variations in the travel times of both P and S waves (Vp and Vs) for 14 months at Toyohashi (central Japan) with a continuously operating vibration source that could efficiently produce both P and S waves. Seismic waves produced by the source, which is called the accurately controlled routinely operated signal system (ACROSS), are recorded by three nearby seismic stations, and the travel time variation at each station is estimated using the transfer function calculated from the recorded data. Long-term variations in Vp and Vs are observed and can be interpreted by the change in fluid saturation and crack density of the subsurface rocks. The variation in fluid saturation and crack density are consistent with that in the groundwater level, which is measured at the station nearest to the ACROSS. Short-term responses to rainfalls are observed at the station nearest to the ACROSS system; the interpretation of the changes in crack density and saturation is inconsistent with the ground water observation, partly owing to the initial response to rainfall. This can be interpreted as an air–water mixture within pores or cracks on a fine scale.


2020 ◽  
Author(s):  
Rina Suzuki ◽  
Koshun Yamaoka ◽  
Shuhei Tsuji ◽  
Toshiki Watanabe

Abstract We analyzed the temporal variation in the travel times of both the P and S waves (Vp and Vs) for 14 months at Toyohashi (central Japan) with a continuously operating vibration source that could produce both P and S waves efficiently. The seismic waves produced by the source, which is named ACROSS (accurately-controlled routinely-operated signal system), were recorded by three nearby seismic stations, and the travel time variation at each station was calculated using the transfer function calculated from the recorded data. We observed the seasonal variations in the Vp and Vs for all the stations—which can be interpreted using the change in the fluid saturation and crack density of subsurface rocks—are consistent with the variation in the ground water level. The short-term responses to rainfall are observed at the nearest station; the interpretation of the changes in crack density and saturation is inconsistent with the ground water observation partly due the initial response to rainfall. This can be interpreted as an air-water mixture within pores or cracks on a fine scale.



Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. N17-N29 ◽  
Author(s):  
Shahpoor Moradi ◽  
Kristopher A. Innanen

Anelastic properties of reservoir rocks are important and sensitive indicators of fluid saturation and viscosity changes due (for instance) to steam injection. The description of seismic waves propagating through viscoelastic continua is quite complex, involving a range of unique homogeneous and inhomogeneous modes. This is true even in the relatively simple theoretical environment of amplitude variation with offset (AVO) analysis. For instance, a complete treatment of the problem of linearizing the solutions of the low-loss viscoelastic Zoeppritz equations to obtain an extended Aki-Richards equations (one that is in accord with the appropriate complex Snell’s law) is lacking in the literature. Also missing is a clear analytical path allowing such forms to be reconciled with more general volume scattering pictures of viscoelastic seismic wave propagation. Our analysis, which provides these two missing elements, leads to approximate reflection and transmission coefficients for the P- and type-I S-waves. These involve additional, complex terms alongside those of the standard isotropic-elastic Aki-Richards equations. The extra terms were shown to have a significant influence on reflection strengths, particularly when the degree of inhomogeneity was high. The particular AVO forms we evaluated were finally shown to be special cases of potentials for volume scattering from viscoelastic inclusions.



1995 ◽  
Vol 11 (2) ◽  
pp. 233-247 ◽  
Author(s):  
Javier F. Pacheco ◽  
Shri Krishna Singh

The Valley of Mexico is exposed to seismic risk from normal-faulting, large intermediate-depth earthquakes. We explore two approaches to estimate future ground motions from such events at CU, a hill-zone site in the valley. In the first we obtain parameters of an ω2 seismic source model and determine amplification of seismic waves due to local site effects at CU. This permits estimation of Fourier spectrum of expected ground motion at CU from postulated earthquakes. We find that the S-waves suffer an amplification of 2.5 between 0.2 to 3.0 Hz. This amplification is similar to that observed from deep teleseismic events but differs from that obtained from shallow coastal events. In the second approach the available recordings at CU are used as empirical Green's functions (EGF) to synthesize motions from future large earthquakes. This approach is very powerful if the smaller event is truly an empirical Green's function for the postulated earthquake.





2018 ◽  
Author(s):  
Hossein Sahour ◽  
◽  
Mohamed Sultan ◽  
Karem Abdelmohsen ◽  
Sita Karki ◽  
...  


2014 ◽  
Vol 57 (3) ◽  
Author(s):  
Andrea Bizzarri

<p>The focus on the present study is on the point-source approximation of a seismic source. First, we compare the synthetic motions on the free surface resulting from different analytical evolutions of the seismic source (the Gabor signal (G), the Bouchon ramp (B), the Cotton and Campillo ramp (CC), the Yoffe function (Y) and the Liu and Archuleta function (LA)). Our numerical experiments indicate that the CC and the Y functions produce synthetics with larger oscillations and correspondingly they have a higher frequency content. Moreover, the CC and the Y functions tend to produce higher peaks in the ground velocity (roughly of a factor of two). We have also found that the falloff at high frequencies is quite different: it roughly follows ω<span><sup>−2</sup></span> in the case of G and LA functions, it decays more faster than ω<span><sup>−2</sup></span> for the B function, while it is slow than ω<span><sup>−1</sup></span> for both the CC and the Y solutions. Then we perform a comparison of seismic waves resulting from 3-D extended ruptures (both supershear and subshear) obeying to different governing laws against those from a single point-source having the same features. It is shown that the point-source models tend to overestimate the ground motions and that they completely miss the Mach fronts emerging from the supershear transition process. When we compare the extended fault solutions against a multiple point-sources model the agreement becomes more significant, although relevant discrepancies still persist. Our results confirm that, and more importantly quantify how, the point-source approximation is unable to adequately describe the radiation emitted during a real world earthquake, even in the most idealized case of planar fault with homogeneous properties and embedded in a homogeneous, perfectly elastic medium.</p>



2001 ◽  
Vol 09 (04) ◽  
pp. 1407-1416 ◽  
Author(s):  
GIULIANA ROSSI ◽  
ALDO VESNAVER

Converted waves can play a basic role in the traveltime inversion of seismic waves. The sought velocity fields of P and S waves are almost decoupled, when considering pure P and S arrivals: their only connection are the possible common reflecting interfaces in the Earth. Converted waves provide new equations in the linear system to be inverted, which directly relates the two velocity fields. Since the new equations do not introduce additional unknowns, they increase the system rank or its redundancy, so making its solutions better constrained and robust.



Geophysics ◽  
2001 ◽  
Vol 66 (5) ◽  
pp. 1519-1527 ◽  
Author(s):  
Robert Sun ◽  
George A. McMechan

Reflected P‐to‐P and P‐to‐S converted seismic waves in a two‐component elastic common‐source gather generated with a P‐wave source in a two‐dimensional model can be imaged by two independent scalar reverse‐time depth migrations. The inputs to migration are pure P‐ and S‐waves that are extracted by divergence and curl calculations during (shallow) extrapolation of the elastic data recorded at the earth’s surface. For both P‐to‐P and P‐to‐S converted reflected waves, the imaging time at each point is the P‐wave traveltime from the source to that point. The extracted P‐wave is reverse‐time extrapolated and imaged with a P‐velocity model, using a finite difference solution of the scalar wave equation. The extracted S‐wave is reverse‐time extrapolated and imaged similarly, but with an S‐velocity model. Converted S‐wave data requires a polarity correction prior to migration to ensure constructive interference between data from adjacent sources. Synthetic examples show that the algorithm gives satisfactory results for laterally inhomogeneous models.



Author(s):  
Maryam Safarshahi ◽  
Igor B. Morozov

ABSTRACT Empirical models of geometrical-, Q-, t-star, and kappa-type attenuation of seismic waves and ground-motion prediction equations (GMPEs) are viewed as cases of a common empirical standard model describing variation of wave amplitudes with time and frequency. Compared with existing parametric and nonparametric approaches, several new features are included in this model: (1) flexible empirical parameterization with possible nonmonotonous time or distance dependencies; (2) joint inversion for time or distance and frequency dependencies, source spectra, site responses, kappas, and Q; (3) additional constraints removing spurious correlations of model parameters and data residuals with source–receiver distances and frequencies; (4) possible kappa terms for sources as well as for receivers; (5) orientation-independent horizontal- and three-component amplitudes; and (6) adaptive filtering to reduce noise effects. The approach is applied to local and regional S-wave amplitudes in southeastern Iran. Comparisons with previous studies show that conventional attenuation models often contain method-specific biases caused by limited parameterizations of frequency-independent amplitude decays and assumptions about the models, such as smoothness of amplitude variations. Without such assumptions, the frequency-independent spreading of S waves is much faster than inferred by conventional modeling. For example, transverse-component amplitudes decrease with travel time t as about t−1.8 at distances closer than 90 km and as t−2.5 beyond 115 km. The rapid amplitude decay at larger distances could be caused by scattering within the near surface. From about 90 to 115 km distances, the amplitude increases by a factor of about 3, which could be due to reflections from the Moho and within the crust. With more accurate geometrical-spreading and kappa models, the Q factor for the study area is frequency independent and exceeds 2000. The frequency-independent and Q-type attenuation for vertical-component and multicomponent amplitudes is somewhat weaker than for the horizontal components. These observations appear to be general and likely apply to other areas.



2020 ◽  
Vol 224 (2) ◽  
pp. 1443-1463
Author(s):  
Victor Mensah ◽  
Arturo Hidalgo

SUMMARY The accurate numerical modelling of 3-D seismic wave propagation is essential in understanding details to seismic wavefields which are, observed on regional and global scales on the Earth’s surface. The diffusive-viscous wave (DVW) equation was proposed to study the connection between fluid saturation and frequency dependence of reflections and to characterize the attenuation property of the seismic wave in a fluid-saturated medium. The attenuation of DVW is primarily described by the active attenuation parameters (AAP) in the equation. It is, therefore, imperative to acquire these parameters and to additionally specify the characteristics of the DVW. In this paper, quality factor, Q is used to obtain the AAP, and they are compared to those of the visco-acoustic wave. We further derive the 3-D numerical schemes based on a second order accurate finite-volume scheme with a second order Runge–Kutta approximation for the time discretization and a fourth order accurate finite-difference scheme with a fourth order Runge–Kutta approximation for the time discretization. We then simulate the propagation of seismic waves in a 3-D fluid-saturated medium based on the derived schemes. The numerical results indicate stronger attenuation when compared to the visco-acoustic case.



Sign in / Sign up

Export Citation Format

Share Document