The vertical image projection method for migration stretch compensation

2021 ◽  
pp. 1-55
Author(s):  
Arash JafarGandomi
Keyword(s):  
Real Data ◽  
Amplitude Variation ◽  
Seismic Migration ◽  
Amplitude Variation With Offset ◽  
Amplitude Inversion ◽  
Image Projection ◽  
And Migration ◽  
Avo Attributes ◽  
Conditioning Approach ◽  
The Impact

True amplitude inversion is often carried out without taking into account migration distortions to the wavelet. Seismic migration leaves a dip-dependent effect on the wavelet that can cause significant inaccuracies in the inverted impedances obtained from conventional inversion approaches based on 1D vertical convolutional modelling. Neglecting this effect causes misleading inversion results and leakage of dipping noise and migration artifacts from higher frequency bands to the lower frequencies. I have observed that despite dip-dependency of this effect, low-dip and flat events may also suffer if they are contaminated with cross-cutting noise, steep migration artifacts, and smiles. In this paper I propose an efficient, effective and reversible data pre-conditioning approach that accounts for dip-dependency of the wavelet and is applied to migrated images prior to inversion. My proposed method consists of integrating data with respect to the total wavenumber followed by the differentiation with respect to the vertical wavenumber. This process is equivalent to applying a deterministic dip-consistent pre-conditioning that projects the data from the total wavenumber to the vertical wavenumber axis. This preconditioning can be applied to both pre- and post-stack data as well as to amplitude variation with offset (AVO) attributes such as intercept and gradient before inversion. The vertical image projection methodology that I propose here reduces the impact of migration artifacts such as cross-cutting noise and migration smiles and improves inverted impedances in both synthetic and real data examples. In particular I show that neglecting the proposed pre-conditioning leads to anomalously higher impedance values along the steeply dipping structures.

The impact of reservoir scale on amplitude variation with offset

2016 ◽  
Vol 65 (3) ◽  
pp. 736-746 ◽  
Author(s):  
Chao Xu ◽  
Jianxin Wei ◽  
Bangrang Di
Keyword(s):  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
The Impact

scholarly journals Stochastic reservoir characterization using prestack seismic data

Geophysics ◽  
2004 ◽  
Vol 69 (4) ◽  
pp. 978-993 ◽  
Author(s):  
Jo Eidsvik ◽  
Per Avseth ◽  
Henning Omre ◽  
Tapan Mukerji ◽  
Gary Mavko
Keyword(s):  
Reservoir Characterization ◽  
A Priori ◽  
Rock Physics ◽  
Real Data ◽  
Spatial Coupling ◽  
The North ◽  
Amplitude Versus Offset ◽  
Likelihood Model ◽  
Avo Attributes ◽  
The Impact

Reservoir characterization must be based on information from various sources. Well observations, seismic reflection times, and seismic amplitude versus offset (AVO) attributes are integrated in this study to predict the distribution of the reservoir variables, i.e., facies and fluid filling. The prediction problem is cast in a Bayesian setting. The a priori model includes spatial coupling through Markov random field assumptions and intervariable dependencies through nonlinear relations based on rock physics theory, including Gassmann's relation. The likelihood model relating observations to reservoir variables (including lithology facies and pore fluids) is based on approximations to Zoeppritz equations. The model assumptions are summarized in a Bayesian network illustrating the dependencies between the reservoir variables. The posterior model for the reservoir variables conditioned on the available observations is defined by the a priori and likelihood models. This posterior model is not analytically tractable but can be explored by Markov chain Monte Carlo (MCMC) sampling. Realizations of reservoir variables from the posterior model are used to predict the facies and fluid‐filling distribution in the reservoir. A maximum a posteriori (MAP) criterion is used in this study to predict facies and pore‐fluid distributions. The realizations are also used to present probability maps for the favorable (sand, oil) occurrence in the reservoir. Finally, the impact of seismic AVO attributes—AVO gradient, in particular—is studied. The approach is demonstrated on real data from a turbidite sedimentary system in the North Sea. AVO attributes on the interface between reservoir and cap rock are extracted from 3D seismic AVO data. The AVO gradient is shown to be valuable in reducing the ambiguity between facies and fluids in the prediction.

Traveltime and relative geometric spreading approximation in elastic orthorhombic medium

Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. C153-C162 ◽  
Author(s):  
Shibo Xu ◽  
Alexey Stovas ◽  
Hitoshi Mikada
Keyword(s):  
Seismic Data ◽  
Real Data ◽  
P Wave ◽  
Form Expression ◽  
Amplitude Variation ◽  
Rational Form ◽  
Amplitude Variation With Offset ◽  
Practical Applications ◽  
Numerical Tests ◽  
Geometric Spreading

Wavefield properties such as traveltime and relative geometric spreading (traveltime derivatives) are highly essential in seismic data processing and can be used in stacking, time-domain migration, and amplitude variation with offset analysis. Due to the complexity of an elastic orthorhombic (ORT) medium, analysis of these properties becomes reasonably difficult, where accurate explicit-form approximations are highly recommended. We have defined the shifted hyperbola form, Taylor series (TS), and the rational form (RF) approximations for P-wave traveltime and relative geometric spreading in an elastic ORT model. Because the parametric form expression for the P-wave vertical slowness in the derivation is too complicated, TS (expansion in offset) is applied to facilitate the derivation of approximate coefficients. The same approximation forms computed in the acoustic ORT model also are derived for comparison. In the numerical tests, three ORT models with parameters obtained from real data are used to test the accuracy of each approximation. The numerical examples yield results in which, apart from the error along the y-axis in ORT model 2 for the relative geometric spreading, the RF approximations all are very accurate for all of the tested models in practical applications.

The Shuey-Radon transform

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. V197-V206 ◽  
Author(s):  
Ali Gholami ◽  
Milad Farshad
Keyword(s):  
Radon Transform ◽  
Seismic Data ◽  
Matching Pursuit ◽  
Real Data ◽  
Transform Domain ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Linear System Of Equations ◽  
Matching Pursuit Algorithm ◽  
Zero Offset

The traditional hyperbolic Radon transform (RT) decomposes seismic data into a sum of constant amplitude basis functions. This limits the performance of the transform when dealing with real data in which the reflection amplitudes include the amplitude variation with offset (AVO) variations. We adopted the Shuey-Radon transform as a combination of the RT and Shuey’s approximation of reflectivity to accurately model reflections including AVO effects. The new transform splits the seismic gather into three Radon panels: The first models the reflections at zero offset, and the other two panels add capability to model the AVO gradient and curvature. There are two main advantages of the Shuey-Radon transform over similar algorithms, which are based on a polynomial expansion of the AVO response. (1) It is able to model reflections more accurately. This leads to more focused coefficients in the transform domain and hence provides more accurate processing results. (2) Unlike polynomial-based approaches, the coefficients of the Shuey-Radon transform are directly connected to the classic AVO parameters (intercept, gradient, and curvature). Therefore, the resulting coefficients can further be used for interpretation purposes. The solution of the new transform is defined via an underdetermined linear system of equations. It is formulated as a sparsity-promoting optimization, and it is solved efficiently using an orthogonal matching pursuit algorithm. Applications to different numerical experiments indicate that the Shuey-Radon transform outperforms the polynomial and conventional RTs.

Reducing 3-D acquisition footprint for 3-D DMO and 3-D prestack migration

Geophysics ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1177-1183 ◽  
Author(s):  
Anat Canning ◽  
Gerald H. F. Gardner
Keyword(s):  
Spatial Distribution ◽  
Velocity Analysis ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Prestack Migration ◽  
Avo Analysis ◽  
And Migration

The acquisition patterns of 3-D surveys often have a significant effect on the results of dip moveout (DMO) or prestack migration. When the spatial distribution of input traces is irregular, results from DMO and migration are contaminated by artifacts. In many cases, the footprint of the acquisition patterns can be seen on the migrated section and may result in incorrect interpretation. This phenomena also has a very significant effect on the feasibility of conducting amplitude variation with offset (AVO) analysis after 3-D prestack migration or after 3-D DMO, and also may affect velocity analysis. We propose a simple enhancement to migration and DMO programs that acts to minimize acquisition artifacts.

Iterative multiparameter waveform inversion of precritical reflection data using prestack time Kirchhoff approximation

Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. R15-R27 ◽  
Author(s):  
Hassan Khaniani ◽  
John C. Bancroft ◽  
Eric von Lunen
Keyword(s):  
Waveform Inversion ◽  
Computational Cost ◽  
Real Data ◽  
Forward Modeling ◽  
Kirchhoff Approximation ◽  
Amplitude Variation With Offset ◽  
Reflection Data ◽  
Elastic Wave Scattering ◽  
Iterative Inversion ◽  
And Migration

We have studied elastic wave scattering and iterative inversion in the context of the Kirchhoff approximation. The approach is more consistent with the weak-contrast reflectivity functions of Zoeppritz equations as compared to the Born approximation. To reduce the computational cost associated with inversion, we demonstrated the use of amplitude-variation-with-offset (AVO) analysis, prestack time migrations (PSTMs), and the corresponding forward modeling in an iterative scheme. Forward modeling and migration/inversion operators are based on the double-square-root (DSR) equations of PSTM and linearized reflectivity functions. All operators involved in the inversion, including the background model for DSR and AVO, are defined in P-to-P traveltime and are updated at each iteration. Our method is practical for real data applications because all operators of the inversion are known to be applicable for standard methods. We have evaluated the inversion on synthetic and real data using the waveform characteristics of P-to-P and P-to-S data.

Seismic amplitude variation with offset inversion in a vertical transverse isotropy medium

2019 ◽  
Vol 7 (3) ◽  
pp. T581-T593 ◽  
Author(s):  
Mark Sams ◽  
Annushia Annamalai ◽  
Jeremy Gallop
Keyword(s):  
Transverse Isotropy ◽  
Low Frequency ◽  
Seismic Inversion ◽  
Well Log ◽  
Amplitude Variation ◽  
Frequency Model ◽  
Isotropic Solution ◽  
Amplitude Variation With Offset ◽  
Vertical Transverse Isotropy ◽  
The Impact

Vertical transverse isotropy (VTI) will affect seismic inversion, but it is not possible to solve for the full set of anisotropic elastic parameters from amplitude variation with offset inversion because there exists an isotropic solution to every VTI problem. We can easily approximate the pseudoisotropic properties that result from the isotropic solution to the anisotropic problem for well-log data. We can then use these well-log properties to provide a low-frequency model for inversion and/or a framework for interpreting either absolute or relative inversion results. This, however, requires prior knowledge of the anisotropic properties, which are often unavailable or poorly constrained. If we ignore anisotropy and assume that the amplitude variations caused by VTI are going to be accounted for by effective wavelets, the inversion results would be in error: The impact of anisotropy is not merely a case of linear scaling of seismic amplitudes for any particular angle range. Ignoring VTI does not affect the prediction of acoustic impedance, but it does affect predictions of [Formula: see text] and density. For realistic values of anisotropy, these errors can be significant, such as predicting oil instead of brine. If the anisotropy of the rocks is known, then we can invert for the true vertical elastic properties using the known anisotropy coefficients through a facies-based inversion. This can produce a more accurate result than solving for pseudoelastic properties, and it can take advantage of the sometimes increased separation of isotropic and anisotropic rocks in the pseudoisotropic elastic domain. Because the effect of anisotropy will vary depending on the strength of the anisotropy and the distribution of the rocks, we strongly recommend forward modeling for each case prior to inversion to understand the potential impact on the study objectives.

Nonlinear amplitude-variation-with-offset inversion for Lamé parameters using a direct inversion method

2017 ◽  
Vol 5 (3) ◽  
pp. SL57-SL67 ◽  
Author(s):  
Guangsen Cheng ◽  
Xingyao Yin ◽  
Zhaoyun Zong
Keyword(s):  
Nonlinear Equation ◽  
Real Data ◽  
Seismic Inversion ◽  
Inversion Method ◽  
Nonlinear Inversion ◽  
Reservoir Prediction ◽  
Amplitude Variation ◽  
Linear Inversion ◽  
Amplitude Variation With Offset ◽  
Avo Inversion

Prestack seismic inversion is widely used in fluid indication and reservoir prediction. Compared with linear inversion, nonlinear inversion is more precise and can be applied to high-contrast situations. The inversion results can be affected by the parameters’ sensitivity, so the parameterization of nonlinear equations is very significant. Considering the poor nonlinear amplitude-variation-with-offset (AVO) inversion results of impedance and velocity parameters, we adjust the parameters of the nonlinear equation, avoid the inaccuracy caused by parameters sensitivity and get the ideal nonlinear AVO inversion results of the Lamé parameters. The feasibility and stability of the nonlinear equation based on the Lamé parameters and method are verified by the model and the real data examples. The resolution and the lateral continuity of nonlinear inversion results are better compared with the linear inversion results.

A method to compensate for migration stretch to improve the resolution of amplitude variation with offset, S-impedance (ZS), and density (ρ)

2020 ◽  
Vol 8 (4) ◽  
pp. T687-T699
Author(s):  
Swetal Patel ◽  
Francis Oyebanji ◽  
Kurt J. Marfurt
Keyword(s):  
Matching Pursuit ◽  
Azimuthal Anisotropy ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Prestack Migration ◽  
Fort Worth Basin ◽  
Seismic Wavelet ◽  
Compensation Factor ◽  
Nmo Correction ◽  
And Migration

Because of their improved leverage against ground roll and multiples, as well as the ability to estimate azimuthal anisotropy, wide-azimuth 3D seismic surveys routinely now are acquired over most resource plays. For a relatively shallow target, most of these surveys can be considered to be long offset as well, containing incident angles up to 45°. Unfortunately, effective use of the far-offset data often is compromised by noise and normal moveout (NMO) (or, more accurately, prestack migration) stretch. The conventional NMO correction is well-known to decrease the frequency content and distort the seismic wavelet at far offsets, sometimes giving rise to tuning effects. Most quantitative interpreters work with prestack migrated gathers rather than unmigrated NMO-corrected gathers. However, prestack migration of flat reflectors suffers from the same limitation called migration stretch. Migration stretch leads to lower S-impedance ([Formula: see text]) and density ([Formula: see text]) resolution estimated from inversion, misclassification of amplitude variation with offset (AVO) types, and infidelity in amplitude variation with azimuth (AVAZ) inversion results. We have developed a matching pursuit algorithm commonly used in spectral decomposition to correct the migration stretch by scaling the stretched wavelets using a wavelet compensation factor. The method is based on hyperbolic moveout approximation. The corrected gathers show increased resolution and higher fidelity amplitudes at the far offsets leading to improvement in AVO classification. Correction for migration stretch rather than conventional “stretch-mute” corrections provides three advantages: (1) preservation of far angles required for accurate [Formula: see text] inversion, (2) improvement in the vertical resolution of [Formula: see text] and [Formula: see text] volumes, and (3) preservation of far angles that provide greater leverage against multiples. We apply our workflow to data acquired in the Fort Worth Basin and retain incident angles up to 42° at the Barnett Shale target. Comparing [Formula: see text], [Formula: see text], and [Formula: see text] of the original gather and migration stretch-compensated data, we find an insignificant improvement in [Formula: see text], but a moderate to significant improvement in resolution of [Formula: see text] and [Formula: see text]. The method is valid for reservoirs that exhibit a dip of no more than 2°. Consistent improvement is observed in resolving thick beds, but the method might introduce amplitude anomalies at far offsets for tuning beds.

Sign in / Sign up

Export Citation Format

Share Document