SH-SH wave inversion for S-wave velocity and density

Geophysics ◽  
2022 ◽  
pp. 1-59
Author(s):  
Fucai Dai ◽  
Feng Zhang ◽  
Xiangyang Li

SS-waves (SV-SV waves and SH-SH waves) are capable of inverting S-wave velocity ( VS) and density ( ρ) because they are sensitive to both parameters. SH-SH waves can be separated from multicomponent data sets more effectively than the SV-SV wave because the former is decoupled from the PP-wave in isotropic media. In addition, the SH-SH wave can be better modeled than the SV-SV wave in the case of strong velocity/impedance contrast because the SV-SV wave has multicritical angles, some of which can be quite small when velocity/ impedance contrast is strong. We derived an approximate equation of the SH-SH wave reflection coefficient as a function of VS and ρ in natural logarithm variables. The approximation has high accuracy, and it enables the inversion of VS and ρ in a direct manner. Both coefficients corresponding to VS and ρ are “model-parameter independent” and thus there is no need for prior estimate of any model parameter in inversion. Then, we developed an SH-SH wave inversion method, and demonstrated it by using synthetic data sets and a real SH-SH wave prestack data set from the west of China. We found that VS and ρ can be reliably estimated from the SH-SH wave of small angles.

Geophysics ◽  
2003 ◽  
Vol 68 (1) ◽  
pp. 185-198 ◽  
Author(s):  
Arild Buland ◽  
Henning Omre

A new linearized AVO inversion technique is developed in a Bayesian framework. The objective is to obtain posterior distributions for P‐wave velocity, S‐wave velocity, and density. Distributions for other elastic parameters can also be assessed—for example, acoustic impedance, shear impedance, and P‐wave to S‐wave velocity ratio. The inversion algorithm is based on the convolutional model and a linearized weak contrast approximation of the Zoeppritz equation. The solution is represented by a Gaussian posterior distribution with explicit expressions for the posterior expectation and covariance; hence, exact prediction intervals for the inverted parameters can be computed under the specified model. The explicit analytical form of the posterior distribution provides a computationally fast inversion method. Tests on synthetic data show that all inverted parameters were almost perfectly retrieved when the noise approached zero. With realistic noise levels, acoustic impedance was the best determined parameter, while the inversion provided practically no information about the density. The inversion algorithm has also been tested on a real 3‐D data set from the Sleipner field. The results show good agreement with well logs, but the uncertainty is high.


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 ◽  
1993 ◽  
Vol 58 (5) ◽  
pp. 713-719 ◽  
Author(s):  
Ghassan I. Al‐Eqabi ◽  
Robert B. Herrmann

The objective of this study is to demonstrate that a laterally varying shallow S‐wave structure, derived from the dispersion of the ground roll, can explain observed lateral variations in the direct S‐wave arrival. The data set consists of multichannel seismic refraction data from a USGS-GSC survey in the state of Maine and the province of Quebec. These data exhibit significant lateral changes in the moveout of the ground‐roll as well as the S‐wave first arrivals. A sequence of surface‐wave processing steps are used to obtain a final laterally varying S‐wave velocity model. These steps include visual examination of the data, stacking, waveform inversion of selected traces, phase velocity adjustment by crosscorrelation, and phase velocity inversion. These models are used to predict the S‐wave first arrivals by using two‐dimensional (2D) ray tracing techniques. Observed and calculated S‐wave arrivals match well over 30 km long data paths, where lateral variations in the S‐wave velocity in the upper 1–2 km are as much as ±8 percent. The modeled correlation between the lateral variations in the ground‐roll and S‐wave arrival demonstrates that a laterally varying structure can be constrained by using surface‐wave data. The application of this technique to data from shorter spreads and shallower depths is discussed.


2020 ◽  
Vol 91 (6) ◽  
pp. 3622-3633
Author(s):  
Rufus D. Catchings ◽  
Mark R. Goldman ◽  
Jamison H. Steidl ◽  
Joanne H. Chan ◽  
Amir A. Allam ◽  
...  

Abstract The 2019 Ridgecrest, California, earthquake sequence included Mw 6.4 and 7.1 earthquakes that occurred on successive days beginning on 4 July 2019. These two largest earthquakes of the sequence occurred on orthogonal faults that ruptured the Earth’s surface. To better evaluate the 3D subsurface fault structure, (P- and S-wave) velocity, 3D and temporal variations in seismicity, and other important aspects of the earthquake sequence, we recorded aftershocks and ambient noise using up to 461 three-component nodal seismographs for about two months, beginning about one day after the Mw 7.1 mainshock. The ∼30,000Mw≥1 earthquakes that were recorded on the dense arrays provide an unusually large volume of data with which to evaluate the earthquake sequence. This report describes the recording arrays and is intended to provide metadata for researchers interested in evaluating various aspects of the 2019 Ridgecrest earthquake sequence using the nodal data set.


2019 ◽  
Vol 491 (4) ◽  
pp. 5238-5247 ◽  
Author(s):  
X Saad-Olivera ◽  
C F Martinez ◽  
A Costa de Souza ◽  
F Roig ◽  
D Nesvorný

ABSTRACT We characterize the radii and masses of the star and planets in the Kepler-59 system, as well as their orbital parameters. The star parameters are determined through a standard spectroscopic analysis, resulting in a mass of $1.359\pm 0.155\, \mathrm{M}_\odot$ and a radius of $1.367\pm 0.078\, \mathrm{R}_\odot$. The obtained planetary radii are $1.5\pm 0.1\, R_\oplus$ for the inner and $2.2\pm 0.1\, R_\oplus$ for the outer planet. The orbital parameters and the planetary masses are determined by the inversion of Transit Timing Variations (TTV) signals. We consider two different data sets: one provided by Holczer et al. (2016), with TTVs only for Kepler-59c, and the other provided by Rowe et al. (2015), with TTVs for both planets. The inversion method applies an algorithm of Bayesian inference (MultiNest) combined with an efficient N-body integrator (Swift). For each of the data set, we found two possible solutions, both having the same probability according to their corresponding Bayesian evidences. All four solutions appear to be indistinguishable within their 2-σ uncertainties. However, statistical analyses show that the solutions from Rowe et al. (2015) data set provide a better characterization. The first solution infers masses of $5.3_{-2.1}^{+4.0}~M_{\mathrm{\oplus }}$ and $4.6_{-2.0}^{+3.6}~M_{\mathrm{\oplus }}$ for the inner and outer planet, respectively, while the second solution gives masses of $3.0^{+0.8}_{-0.8}~M_{\mathrm{\oplus }}$ and $2.6^{+0.9}_{-0.8}~M_{\mathrm{\oplus }}$. These values point to a system with an inner super-Earth and an outer mini-Neptune. A dynamical study shows that the planets have almost co-planar orbits with small eccentricities (e < 0.1), close to the 3:2 mean motion resonance. A stability analysis indicates that this configuration is stable over million years of evolution.


Geophysics ◽  
1995 ◽  
Vol 60 (3) ◽  
pp. 796-809 ◽  
Author(s):  
Zhong‐Min Song ◽  
Paul R. Williamson ◽  
R. Gerhard Pratt

In full‐wave inversion of seismic data in complex media it is desirable to use finite differences or finite elements for the forward modeling, but such methods are still prohibitively expensive when implemented in 3-D. Full‐wave 2-D inversion schemes are of limited utility even in 2-D media because they do not model 3-D dynamics correctly. Many seismic experiments effectively assume that the geology varies in two dimensions only but generate 3-D (point source) wavefields; that is, they are “two‐and‐one‐half‐dimensional” (2.5-D), and this configuration can be exploited to model 3-D propagation efficiently in such media. We propose a frequency domain full‐wave inversion algorithm which uses a 2.5-D finite difference forward modeling method. The calculated seismogram can be compared directly with real data, which allows the inversion to be iterated. We use a descents‐related method to minimize a least‐squares measure of the wavefield mismatch at the receivers. The acute nonlinearity caused by phase‐wrapping, which corresponds to time‐domain cycle‐skipping, is avoided by the strategy of either starting the inversion using a low frequency component of the data or constructing a starting model using traveltime tomography. The inversion proceeds by stages at successively higher frequencies across the observed bandwidth. The frequency domain is particularly efficient for crosshole configurations and also allows easy incorporation of attenuation, via complex velocities, in both forward modeling and inversion. This also requires the introduction of complex source amplitudes into the inversion as additional unknowns. Synthetic studies show that the iterative scheme enables us to achieve the theoretical maximum resolution for the velocity reconstruction and that strongly attenuative zones can be recovered with reasonable accuracy. Preliminary results from the application of the method to a real data set are also encouraging.


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.


2021 ◽  
Vol 11 (8) ◽  
pp. 3571
Author(s):  
Genggeng Wen ◽  
Kuiyuan Wan ◽  
Shaohong Xia ◽  
Huilong Xu ◽  
Chaoyan Fan ◽  
...  

The detailed studies of converted S-waves recorded on the Ocean Bottom Seismometer (OBS) can provide evidence for constraining lithology and geophysical properties. However, the research of converted S-waves remains a weakness, especially the S-waves’ inversion. In this study, we applied a travel-time inversion method of converted S-waves to obtain the crustal S-wave velocity along the profile NS5. The velocities of the crust are determined by the following four aspects: (1) modelling the P-wave velocity, (2) constrained sediments Vp/Vs ratios and S-wave velocity using PPS phases, (3) the correction of PSS phases’ travel-time, and (4) appropriate parameters and initial model are selected for inversion. Our results show that the vs. and Vp/Vs of the crust are 3.0–4.4 km/s and 1.71–1.80, respectively. The inversion model has a similar trend in velocity and Vp/Vs ratios with the forward model, due to a small difference with ∆Vs of 0.1 km/s and ∆Vp/Vs of 0.03 between two models. In addition, the high-resolution inversion model has revealed many details of the crustal structures, including magma conduits, which further supports our method as feasible.


Sign in / Sign up

Export Citation Format

Share Document