Joint Inversion Based P/S Wave Velocity Data Processing to Test a New Rock Physical Model Describing Acoustic Hysteresis

Author(s):  
J. Somogyine Molnar ◽  
A. Kiss ◽  
M. Dobroka
Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. KS11-KS22 ◽  
Author(s):  
Nicola Piana Agostinetti ◽  
Alberto Malinverno

We use teleseismic P-to-S converted waves from a permanent station to estimate the uncertainties in a 1D elastic model of the shallow crust (0–7 km depth) obtained from the inversion of receiver function (RF) data. Our earth model consists of layers with a constant S-wave velocity [Formula: see text] and P- to S-wave velocity ratio ([Formula: see text]). We apply a Bayesian formulation and transdimensional Monte Carlo sampling to compute the posterior uncertainties of the earth model. The model uncertainties rely on a realistic representation of the data uncertainties, and we estimate directly from the stacking of the teleseismic data, a full-error covariance matrix. To explore the effect of the number of teleseismic events and the RF frequency content, we compare the results of inverting a single RF computed for a cut-off filter frequency of 4 Hz with the joint inversion of four RFs computed from independent ensembles in a larger pool of events for cut-off frequencies of 0.5, 1, 2, and 4 Hz. The inversion results are compared with the lithostratigraphy and sonic-log measurements from a 7 km deep borehole drilled near the seismic station. The inversion of a single RF results in larger uncertainties in the recovered [Formula: see text] profile and in the depth to seismic discontinuities compared with the multifrequency inversion. Moreover, the multifrequency inversion predicts more accurately the depth to a velocity inversion at approximately 6 km below the surface and matches more closely the borehole sonic-log data. Our results indicate that RF data can be used to map shallow (3–5 km depth) crustal interfaces with uncertainties in the order of 300–500 m, whereas uncertainties are consistently smaller (<300 m) for interfaces in the top kilometer.


Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 1095-1107 ◽  
Author(s):  
Ilya Tsvankin ◽  
Leon Thomsen

In anisotropic media, the short‐spread stacking velocity is generally different from the root‐mean‐square vertical velocity. The influence of anisotropy makes it impossible to recover the vertical velocity (or the reflector depth) using hyperbolic moveout analysis on short‐spread, common‐midpoint (CMP) gathers, even if both P‐ and S‐waves are recorded. Hence, we examine the feasibility of inverting long‐spread (nonhyperbolic) reflection moveouts for parameters of transversely isotropic media with a vertical symmetry axis. One possible solution is to recover the quartic term of the Taylor series expansion for [Formula: see text] curves for P‐ and SV‐waves, and to use it to determine the anisotropy. However, this procedure turns out to be unstable because of the ambiguity in the joint inversion of intermediate‐spread (i.e., spreads of about 1.5 times the reflector depth) P and SV moveouts. The nonuniqueness cannot be overcome by using long spreads (twice as large as the reflector depth) if only P‐wave data are included. A general analysis of the P‐wave inverse problem proves the existence of a broad set of models with different vertical velocities, all of which provide a satisfactory fit to the exact traveltimes. This strong ambiguity is explained by a trade‐off between vertical velocity and the parameters of anisotropy on gathers with a limited angle coverage. The accuracy of the inversion procedure may be significantly increased by combining both long‐spread P and SV moveouts. The high sensitivity of the long‐spread SV moveout to the reflector depth permits a less ambiguous inversion. In some cases, the SV moveout alone may be used to recover the vertical S‐wave velocity, and hence the depth. Success of this inversion depends on the spreadlength and degree of SV‐wave velocity anisotropy, as well as on the constraints on the P‐wave vertical velocity.


2002 ◽  
Vol 10 (03) ◽  
pp. 331-347 ◽  
Author(s):  
QIZHEN DU ◽  
HUIZHU YANG ◽  
YUAN DONG

The paper presents estimates of the S-wave velocity and the crack density at which fractured reservoirs begin to play an important role in oil exploration. Transverse isotropy with a horizontal axis of symmetry (HTI) is the simplest azimuthally anisotropic model used to describe fractured reservoirs that contain parallel vertical cracks. A double profile concept is used to develop an equation for the P-S wave normal-moveout (NMO) velocity. The azimuthal NMO velocities of the P- and P-S waves can then be used to estimate the velocities of the S-waves and Thomsen's coefficient, γ. For multilayered media, a recursive equation is developed for the NMO velocity in each layer. The numerical results indicate that the S-wave NMO velocity can be accurately estimated using the P- and P-S wave NMO velocities in HTI media. An important parameter of fracture systems that can be measured from seismic data is the crack density which can be estimated using the NMO velocities of the P- and S-waves from horizontal reflectors. Therefore, fractures can be completely characterized by the joint inversion of P-waves and converted P-S waves in HTI media.


2021 ◽  
Vol 11 (15) ◽  
pp. 6712
Author(s):  
Chao Zhang ◽  
Ting Lei ◽  
Yi Wang

Surface-wave dispersion and the Z/H ratio are important parameters used to resolve the Earth’s structure, especially for S-wave velocity. Several previous studies have explored using joint inversion of these two datasets. However, all of these studies used a 1-D depth-sensitivity kernel, which lacks precision when the structure is laterally heterogeneous. Adjoint tomography (i.e., full-waveform inversion) is a state-of-the-art imaging method with a high resolution. It can obtain better-resolved lithospheric structures beyond the resolving ability of traditional ray-based travel-time tomography. In this study, we present a systematic investigation of the 2D sensitivities of the surface wave phase and Z/H ratio using the adjoint-state method. The forward-modeling experiments indicated that the 2D phase and Z/H ratio had different sensitivities to the S-wave velocity. Thus, a full-waveform joint-inversion scheme of surface waves with phases and a Z/H ratio was proposed to take advantage of their complementary sensitivities to the Earth’s structure. Both applications to synthetic data sets in large- and small-scale inversions demonstrated the advantage of the joint inversion over the individual inversions, allowing for the creation of a more unified S-wave velocity model. The proposed joint-inversion scheme offers a computationally efficient and inexpensive alternative to imaging fine-scale shallow structures beneath a 2D seismic array.


2019 ◽  
Vol 124 (12) ◽  
pp. 12997-13020 ◽  
Author(s):  
Christian Baillard ◽  
William S.D. Wilcock ◽  
Adrien F. Arnulf ◽  
Maya Tolstoy ◽  
Felix Waldhauser

Sign in / Sign up

Export Citation Format

Share Document