Simultaneous inversion of PP and PS seismic data

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. R1-R10 ◽  
Author(s):  
Helene Hafslund Veire ◽  
Martin Landrø
Keyword(s):  
Seismic Data ◽  
Reservoir Characterization ◽  
Model Building ◽  
Joint Inversion ◽  
P Wave ◽  
Inversion Method ◽  
System Of Equations ◽  
S Wave ◽  
Data Set ◽  
Simultaneous Inversion

Elastic parameters derived from seismic data are valuable input for reservoir characterization because they can be related to lithology and fluid content of the reservoir through empirical relationships. The relationship between physical properties of rocks and fluids and P-wave seismic data is nonunique. This leads to large uncertainties in reservoir models derived from P-wave seismic data. Because S- waves do not propagate through fluids, the combined use of P-and S-wave seismic data might increase our ability to derive fluid and lithology effects from seismic data, reducing the uncertainty in reservoir characterization and thereby improving 3D reservoir model-building. We present a joint inversion method for PP and PS seismic data by solving approximated linear expressions of PP and PS reflection coefficients simultaneously using a least-squares estimation algorithm. The resulting system of equations is solved by singular-value decomposition (SVD). By combining the two independent measurements (PP and PS seismic data), we stabilize the system of equations for PP and PS seismic data separately, leading to more robust parameter estimation. The method does not require any knowledge of PP and PS wavelets. We tested the stability of this joint inversion method on a 1D synthetic data set. We also applied the methodology to North Sea multicomponent field data to identify sand layers in a shallow formation. The identified sand layers from our inverted sections are consistent with observations from nearby well logs.

scholarly journals Quadri-joint inversion: Method and application to the Big Sky 9C 3D data set in northern Montana

2018 ◽  
Vol 6 (4) ◽  
pp. SN101-SN118 ◽  
Author(s):  
Vincent Clochard ◽  
Bryan C. DeVault ◽  
David Bowen ◽  
Nicolas Delépine ◽  
Kanokkarn Wangkawong
Keyword(s):  
Reservoir Characterization ◽  
Joint Inversion ◽  
Wave Impedance ◽  
P Wave ◽  
Seismic Survey ◽  
Time Shift ◽  
Inversion Method ◽  
S Wave ◽  
Data Set ◽  
Long Term Storage

The Kevin Dome [Formula: see text] storage project, located in northern Montana, attempted to characterize the Duperow Formation as a potential long-term storage zone for injected [Formula: see text]. A multicomponent (9C) seismic survey was acquired for the Big Sky Carbon Sequestration Partnership over a portion of the Kevin Dome using P- and S-wave sources. Prestack migrated PP, PS, SH, and SV data sets were generated. We then applied several stratigraphic inversion workflows using one or several kinds of seismic wavefield at the same time resulting in joint inversions of each data set. The aim of our study is to demonstrate the benefits of doing quadri-joint inversion of PP-, PS-, SH-, and SV-wavefields for the recovery of the elastic earth parameters, especially the S-wave impedance and density. These are crucial parameters because they can help determine lithology and porefill in the reservoir characterization workflow. Because the inversion workflow always uses the original seismic data recorded in its own time domain, it is necessary to compute registration laws between PP-PS-, PP-SH-, and PP-SV-wavefields using a time shift computation procedure (warping) based on inverted S-wave impedances from inversion of a single wavefield. This generated a significant improvement over methods that rely on attempting to match trace waveforms that may have a different phase, frequency content, and polarity. Finally, we wanted to investigate the reliability of the quadri-joint inversion results in the Bakken/Banff Formations, which have less lateral geologic variation than the underlying Duperow target. This interval shares many of the geophysical characterization challenges common to shale reservoirs in other North American basins. We computed geomechanical parameters, such as Poisson’s ratio and Young’s modulus, which are a proxy for brittleness. Comparison of these results with independent laboratory measurements in the Bakken interval demonstrates the superiority of the quadri-joint inversion method to the traditional inversion using P-wave data only.

Simultaneous inversion for model geometry and elastic parameters

Geophysics ◽  
1999 ◽  
Vol 64 (1) ◽  
pp. 182-190 ◽  
Author(s):  
Yanghua Wang
Keyword(s):  
Joint Inversion ◽  
Real Data ◽  
Structural Effect ◽  
P Wave ◽  
Inversion Method ◽  
S Wave ◽  
Elastic Parameters ◽  
Simultaneous Inversion ◽  
Avo Analysis ◽  
The North

Both traveltimes and amplitudes in reflection seismology are used jointly in an inversion to simultaneously invert for the interface geometry and the elastic parameters at the reflectors. The inverse problem has different physical dimensions in both data and model spaces. Practical approaches are proposed to tackle the dimensional difficulties. In using the joint inversion, which may properly take care of the structural effect, one potentially improves the estimates of the subsurface elastic parameters in the traditional analysis of amplitude variation with offset (AVO). Analysis of the elastic parameters estimated, using the ratio of s-wave to P-wave velocity contrasts and the deviation of this parameter from a normal background trend, promises to have application in AVO analysis. The inversion method is demonstrated by application to real data from the North Sea.

Joint inversion of logging-while-drilling multipole acoustic data to determine P- and S-wave velocities in unconsolidated slow formations

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. D553-D560 ◽  
Author(s):  
Yuan-Da Su ◽  
Can Jiang ◽  
Chun-Xi Zhuang ◽  
Song Xu ◽  
Xiao-Ming Tang
Keyword(s):  
Data Processing ◽  
Joint Inversion ◽  
P Wave ◽  
Inversion Method ◽  
Dispersion Characteristics ◽  
S Wave ◽  
Data Set ◽  
Acoustic Data ◽  
Wave Velocities ◽  
Logging While Drilling

We have developed a joint inversion method for logging-while-drilling (LWD) multipole acoustic data processing to simultaneously determine the formation of P- and S-wave velocities. The presence of the LWD tool strongly influences the dispersion characteristics of quadrupole and monopole leaky-P-waves, especially in unconsolidated slow formations. We have verified that an equivalent-tool theory can be adequately used to model the LWD multipole wave dispersion characteristics and can therefore be used to do forward modeling for the inversion. A major advantage of jointly inverting the multipole data sets, as compared with separately inverting each individual data set, is the reduction of uncertainties in the estimated formation of P- and S-wave velocities. We have applied the method to field data processing. The results found that the method not only corrected the dispersion effect in the quadrupole and leaky-P-wave data but also simultaneously obtained the formation of P- and S-wave velocities.

Simultaneous inversion of prestack seismic data for rock properties using simulated annealing

Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1877-1885 ◽  
Author(s):  
Xin‐Quan Ma
Keyword(s):  
Simulated Annealing ◽  
Seismic Data ◽  
Low Frequency ◽  
Optimization Procedure ◽  
P Wave ◽  
Rock Properties ◽  
Inversion Method ◽  
Inversion Algorithm ◽  
Simultaneous Inversion ◽  
Reflection Seismic

A new prestack inversion algorithm has been developed to simultaneously estimate acoustic and shear impedances from P‐wave reflection seismic data. The algorithm uses a global optimization procedure in the form of simulated annealing. The goal of optimization is to find a global minimum of the objective function, which includes the misfit between synthetic and observed prestack seismic data. During the iterative inversion process, the acoustic and shear impedance models are randomly perturbed, and the synthetic seismic data are calculated and compared with the observed seismic data. To increase stability, constraints have been built into the inversion algorithm, using the low‐frequency impedance and background Vs/Vp models. The inversion method has been successfully applied to synthetic and field data examples to produce acoustic and shear impedances comparable to log data of similar bandwidth. The estimated acoustic and shear impedances can be combined to derive other elastic parameters, which may be used for identifying of lithology and fluid content of reservoirs.

Converted-wave model building and imaging based on common-focus-point methodology

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. U139-U149
Author(s):  
Hongwei Liu ◽  
Mustafa Naser Al-Ali ◽  
Yi Luo
Keyword(s):  
Elastic Wave ◽  
Wave Velocity ◽  
Seismic Data ◽  
Model Building ◽  
P Wave ◽  
Rock Properties ◽  
S Wave ◽  
Converted Wave ◽  
Seismic Images ◽  
Focus Point

Seismic images can be viewed as photographs for underground rocks. These images can be generated from different reflections of elastic waves with different rock properties. Although the dominant seismic data processing is still based on the acoustic wave assumption, elastic wave processing and imaging have become increasingly popular in recent years. A major challenge in elastic wave processing is shear-wave (S-wave) velocity model building. For this reason, we have developed a sequence of procedures for estimating seismic S-wave velocities and the subsequent generation of seismic images using converted waves. We have two main essential new supporting techniques. The first technique is the decoupling of the S-wave information by generating common-focus-point gathers via application of the compressional-wave (P-wave) velocity on the converted seismic data. The second technique is to assume one common VP/ VS ratio to approximate two types of ratios, namely, the ratio of the average earth layer velocity and the ratio of the stacking velocity. The benefit is that we reduce two unknown ratios into one, so it can be easily scanned and picked in practice. The PS-wave images produced by this technology could be aligned with the PP-wave images such that both can be produced in the same coordinate system. The registration between the PP and PS images provides cross-validation of the migrated structures and a better estimation of underground rock and fluid properties. The S-wave velocity, computed from the picked optimal ratio, can be used not only for generating the PS-wave images, but also to ensure well registration between the converted-wave and P-wave images.

Bayesian linearized AVO inversion

Geophysics ◽  
2003 ◽  
Vol 68 (1) ◽  
pp. 185-198 ◽  
Author(s):  
Arild Buland ◽  
Henning Omre
Keyword(s):  
Wave Velocity ◽  
Posterior Distribution ◽  
Acoustic Impedance ◽  
Velocity Ratio ◽  
P Wave ◽  
Inversion Method ◽  
Inversion Algorithm ◽  
S Wave ◽  
Data Set ◽  
Avo Inversion

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.

AVO inversion and poroelasticity with P- and S-wave moduli

Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. N17-N24 ◽  
Author(s):  
Zhaoyun Zong ◽  
Xingyao Yin ◽  
Guochen Wu
Keyword(s):  
Seismic Data ◽  
Real Data ◽  
P Wave ◽  
Model Parameters ◽  
Initial Model ◽  
S Wave ◽  
Data Set ◽  
Statistical Probability ◽  
Avo Inversion ◽  
Poroelasticity Theory

The fluid term in the Biot-Gassmann equation plays an important role in reservoir fluid discrimination. The density term imbedded in the fluid term, however, is difficult to estimate because it is less sensitive to seismic amplitude variations. We combined poroelasticity theory, amplitude variation with offset (AVO) inversion, and identification of P- and S-wave moduli to present a stable and physically meaningful method to estimate the fluid term, with no need for density information from prestack seismic data. We used poroelasticity theory to express the fluid term as a function of P- and S-wave moduli. The use of P- and S-wave moduli made the derivation physically meaningful and natural. Then we derived an AVO approximation in terms of these moduli, which can then be directly inverted from seismic data. Furthermore, this practical and robust AVO-inversion technique was developed in a Bayesian framework. The objective was to obtain the maximum a posteriori solution for the P-wave modulus, S-wave modulus, and density. Gaussian and Cauchy distributions were used for the likelihood and a priori probability distributions, respectively. The introduction of a low-frequency constraint and statistical probability information to the objective function rendered the inversion more stable and less sensitive to the initial model. Tests on synthetic data showed that all the parameters can be estimated well when no noise is present and the estimated P- and S-wave moduli were still reasonable with moderate noise and rather smooth initial model parameters. A test on a real data set showed that the estimated fluid term was in good agreement with the results of drilling.

Inversion and interpretation of a 3D seismic data set from the Ouachita Mountains, Oklahoma

Geophysics ◽  
2009 ◽  
Vol 74 (2) ◽  
pp. B37-B45 ◽  
Author(s):  
Abuduwali Aibaidula ◽  
George McMechan
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
3D Visualization ◽  
Wave Impedance ◽  
P Wave ◽  
3D Seismic ◽  
S Wave ◽  
Ouachita Mountains ◽  
Data Set ◽  
3D Seismic Data

Acoustic impedance inversion (AI) and simultaneous angle-dependent inversion (SADI) of a 3D seismic data set characterize reservoirs of Mississippian Morrowan age in the triangle zone of the frontal Ouachita Mountains, Oklahoma. Acoustic impedance of the near-angle seismic data images the 3D spatial distributions of Wapanucka limestone and Cromwell sandstone. Lamé [Formula: see text] ([Formula: see text] and [Formula: see text]) and [Formula: see text] sections are derived from the P-wave and S-wave impedance ([Formula: see text] and [Formula: see text]) sections produced by the SADI. Lithology is identified from the gamma logs and [Formula: see text]. The [Formula: see text], [Formula: see text], and [Formula: see text] are interpreted in terms of a hydrocarbon distribution pattern. The [Formula: see text] is used to identify high [Formula: see text] regions that are consistent with high sand/shale ratio. The estimated impedances and derived Lamé parameter sections are consistent with the interpretation that parts of the Wapanucka limestone and Cromwell sandstone contain potential gas reservoirs in fault-bounded compartments. The Cromwell sandstone contains the main inferred reservoirs; the two largest of these are each [Formula: see text] in pore volume. The inversion results also explain the observed low production in previous wells because those did not sample the best compartments. We propose a single new well location that would penetrate both reservoirs; 3D visualization facilitates this recommendation.

Joint inversion of Rayleigh-wave dispersion and P-wave refraction data for laterally varying layered models

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. EN49-EN59 ◽  
Author(s):  
Daniele Boiero ◽  
Laura Valentina Socco
Keyword(s):  
Rayleigh Wave ◽  
Joint Inversion ◽  
A Priori ◽  
P Wave ◽  
Inversion Method ◽  
Model Parameters ◽  
Data Set ◽  
Wave Refraction ◽  
Velocity Models ◽  
Refraction Data

We implemented a joint inversion method to build P- and S-wave velocity models from Rayleigh-wave and P-wave refraction data, specifically designed to deal with laterally varying layered environments. A priori information available over the site and any physical law to link model parameters can be also incorporated. We tested and applied the algorithm behind the method. The results from a field data set revealed advantages with respect to individual surface-wave analysis (SWA) and body wave tomography (BWT). The algorithm imposed internal consistency for all the model parameters relaxing the required a priori assumptions (i.e., Poisson’s ratio level of confidence in SWA) and the inherent limitations of the two methods (i.e., velocity decreases for BWT).

Estimating elastic properties and attenuation factor from different frequency components of observed seismic data

2019 ◽  
Vol 220 (2) ◽  
pp. 794-805
Author(s):  
Huaizhen Chen
Keyword(s):  
Reflection Coefficient ◽  
Elastic Properties ◽  
Seismic Data ◽  
Wave Attenuation ◽  
Attenuation Factor ◽  
P Wave ◽  
S Wave ◽  
Data Set ◽  
Frequency Dependent ◽  
Frequency Components

SUMMARY Based on an attenuation model, we first express frequency-dependent P- and S-wave attenuation factors as a function of P-wave maximum attenuation factor, and then we re-express P- and S-wave velocities in anelastic media and derive frequency-dependent stiffness parameters in terms of P-wave maximum attenuation factor. Using the derived stiffness parameters, we propose frequency-dependent reflection coefficient in terms of P- and S-wave moduli at critical frequency and P-wave maximum attenuation factor for the case of an interface separating two attenuating media. Based on the derived reflection coefficient, we establish an approach to utilize different frequency components of observed seismic data to estimate elastic properties (P- and S-wave moduli and density) and attenuation factor, and following a Bayesian framework, we construct the objective function and an iterative method is employed to solve the inversion problem. Tests on synthetic data confirm that the proposed approach makes a stable and robust estimation of unknown parameters in the case of seismic data containing a moderate noise/error. Applying the proposed approach to a real data set illustrates that a reliable attenuation factor is obtained from observed seismic data, and the ability of distinguishing oil-bearing reservoirs is improved combining the estimated elastic properties and P-wave attenuation factor.

Sign in / Sign up

Export Citation Format

Share Document