Reverse time migration with an exact two way illumination compensation

Geophysics ◽  
2021 ◽  
pp. 1-67
Author(s):  
Yuzhu Liu ◽  
Weigang Liu ◽  
Zheng Wu ◽  
Jizhong Yang
Keyword(s):  
Side Effects ◽  
Oil And Gas ◽  
Computational Cost ◽  
Receiver Side ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Subsurface Structures ◽  
Time Migration ◽  
Hessian Operator

Reverse time migration (RTM) has been widely used for imaging complex subsurface structures in oil and gas exploration. However, because only the adjoint of the forward Born modeling operator is applied to the seismic data in RTM, the output migration profile is biased in terms of the amplitude. To help partially balance the amplitude performance, the RTM image can be preconditioned with the inverse of the diagonal of the Hessian operator. Yet, existing preconditioning methods do not correctly consider the receiver-side effects, assuming that the receiver coverage is infinite or the velocity model is constant. We therefore provide a comparative study aiming to give a clearer understanding on the importance of incorporating the receiver-side effects by developing a frequency-domain scattering-integral reverse time migration (SI-RTM). In the proposed SI-RTM, the diagonal of the Hessian operator is explicitly computed in its exact formulation, and the source-side wavefield and receiver-side Green’s functions are obtained by solving the two-way wave equation. The computational cost is relatively affordable when compared with the more expensive least-squares RTM. In the comparative counterpart, the diagonal of the Hessian operator is approximated by the source-side illumination. We perform two synthetic numerical examples using an overthrust model and a complex reservoir model; the final migration images were significantly improved when the receiver-side effects were accurately considered. A third application of SI-RTM on one field data set acquired from the East China Sea further demonstrates the importance of incorporating the receiver-side effects in normalizing the RTM image. Findings of this study are expected to provide a theoretical basis for improving the ability of RTM imaging of subsurface structures, thereby critically advancing the application of geophysical techniques for imaging complex environments.

scholarly journals Fast Least-Squares Reverse Time Migration of OBN Down-Going Multiples

2021 ◽  
Vol 9 ◽  
Author(s):  
Yanbao Zhang ◽  
Yike Liu ◽  
Jia Yi ◽  
Xuejian Liu
Keyword(s):  
Least Squares ◽  
Oil And Gas ◽  
Computational Cost ◽  
Ocean Bottom ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Subsurface Structures ◽  
Time Migration ◽  
Resource Exploration ◽  
Imaging Algorithms

Nowadays the ocean bottom node (OBN) acquisition is widely used in oil and gas resource exploration and seismic monitoring. Conventional imaging algorithms of OBN data mainly focus on the processing of up-going primaries and down-going first-order multiples. Up-going multiples and higher-order down-going multiples are generally regarded as noise and should be eliminated or ignored in conventional migration methods. However, multiples carry abundant information about subsurface structures where primaries cannot achieve. To take full advantage of multiples, we propose a migration method using OBN down-going all-order multiples. And then the least-squares optimization algorithm is used to suppress crosstalks. Finally, a phase-encoding-based migration algorithm is developed to cut down the computational cost by blending several common receiver gathers together using random time delays and polarity reversals. Numerical experiments on the complex Marmousi model illustrate that the developed approach can enlarge the imaging area evidently, reduce the computational cost effectively, and enhance the image quality by suppressing crosstalks and improving the resolution.

Application of RTM 3D angle gathers to wide-azimuth data subsalt imaging

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB175-WB182 ◽  
Author(s):  
Yan Huang ◽  
Bing Bai ◽  
Haiyong Quan ◽  
Tony Huang ◽  
Sheng Xu ◽  
...  
Keyword(s):  
Model Updating ◽  
Velocity Model ◽  
Azimuth Angle ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Additional Information ◽  
Time Migration ◽  
Reflection Angle ◽  
Subsalt Imaging

The availability of wide-azimuth data and the use of reverse time migration (RTM) have dramatically increased the capabilities of imaging complex subsalt geology. With these improvements, the current obstacle for creating accurate subsalt images now lies in the velocity model. One of the challenges is to generate common image gathers that take full advantage of the additional information provided by wide-azimuth data and the additional accuracy provided by RTM for velocity model updating. A solution is to generate 3D angle domain common image gathers from RTM, which are indexed by subsurface reflection angle and subsurface azimuth angle. We apply these 3D angle gathers to subsalt tomography with the result that there were improvements in velocity updating with a wide-azimuth data set in the Gulf of Mexico.

Use of RTM full 3D subsurface angle gathers for subsalt velocity update and image optimization: Case study at Shenzi field

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB27-WB39 ◽  
Author(s):  
Zheng-Zheng Zhou ◽  
Michael Howard ◽  
Cheryl Mifflin
Keyword(s):  
Model Building ◽  
Transverse Isotropy ◽  
Velocity Model ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Edge Diffraction ◽  
Time Migration ◽  
Image Optimization ◽  
Free Generation

Various reverse time migration (RTM) angle gather generation techniques have been developed to address poor subsalt data quality and multiarrival induced problems in gathers from Kirchhoff migration. But these techniques introduce new problems, such as inaccuracies in 2D subsurface angle gathers and edge diffraction artifacts in 3D subsurface angle gathers. The unique rich-azimuth data set acquired over the Shenzi field in the Gulf of Mexico enabled the generally artifact-free generation of 3D subsurface angle gathers. Using this data set, we carried out suprasalt tomography and salt model building steps and then produced 3D angle gathers to update the subsalt velocity. We used tilted transverse isotropy RTM with extended image condition to generate full 3D subsurface offset domain common image gathers, which were subsequently converted to 3D angle gathers. The angle gathers were substacked along the subsurface azimuth axis into azimuth sectors. Residual moveout analysis was carried out, and ray-based tomography was used to update velocities. The updated velocity model resulted in improved imaging of the subsalt section. We also applied residual moveout and selective stacking to 3D angle gathers from the final migration to produce an optimized stack image.

First-order particle velocity equations of decoupled P- and S-wavefields and their application in elastic reverse time migration

Geophysics ◽  
2021 ◽  
pp. 1-78
Author(s):  
Zhiyuan Li ◽  
Youshan Liu ◽  
Guanghe Liang ◽  
Guoqiang Xue ◽  
Runjie Wang
Keyword(s):  
Particle Velocity ◽  
Computational Cost ◽  
Staggered Grid ◽  
Additional Stress ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Computational Accuracy ◽  
First Order ◽  
Time Migration ◽  
Stress Variable

The separation of P- and S-wavefields is considered to be an effective approach for eliminating wave-mode cross-talk in elastic reverse-time migration. At present, the Helmholtz decomposition method is widely used for isotropic media. However, it tends to change the amplitudes and phases of the separated wavefields compared with the original wavefields. Other methods used to obtain pure P- and S-wavefields include the application of the elastic wave equations of the decoupled wavefields. To achieve a high computational accuracy, staggered-grid finite-difference (FD) schemes are usually used to numerically solve the equations by introducing an additional stress variable. However, the computational cost of this method is high because a conventional hybrid wavefield (P- and S-wavefields are mixed together) simulation must be created before the P- and S-wavefields can be calculated. We developed the first-order particle velocity equations to reduce the computational cost. The equations can describe four types of particle velocity wavefields: the vector P-wavefield, the scalar P-wavefield, the vector S-wavefield, and the vector S-wavefield rotated in the direction of the curl factor. Without introducing the stress variable, only the four types of particle velocity variables are used to construct the staggered-grid FD schemes, so the computational cost is reduced. We also present an algorithm to calculate the P and S propagation vectors using the four particle velocities, which is simpler than the Poynting vector. Finally, we applied the velocity equations and propagation vectors to elastic reverse-time migration and angle-domain common-image gather computations. These numerical examples illustrate the efficiency of the proposed methods.

From acoustic to elastic inverse extended Born modeling: a first insight in the marine environment

Geophysics ◽  
2021 ◽  
pp. 1-73
Author(s):  
Milad Farshad ◽  
Hervé Chauris
Keyword(s):  
Least Squares ◽  
Marine Environment ◽  
Computational Cost ◽  
Physical Parameters ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Time Migration ◽  
Physical Density ◽  
Acoustic Approximation

Elastic least-squares reverse time migration is the state-of-the-art linear imaging technique to retrieve high-resolution quantitative subsurface images. A successful application requires many migration/modeling cycles. To accelerate the convergence rate, various pseudoinverse Born operators have been proposed, providing quantitative results within a single iteration, while having roughly the same computational cost as reverse time migration. However, these are based on the acoustic approximation, leading to possible inaccurate amplitude predictions as well as the ignorance of S-wave effects. To solve this problem, we extend the pseudoinverse Born operator from acoustic to elastic media to account for the elastic amplitudes of PP reflections and provide an estimate of physical density, P- and S-wave impedance models. We restrict the extension to marine environment, with the recording of pressure waves at the receiver positions. Firstly, we replace the acoustic Green's functions by their elastic version, without modifying the structure of the original pseudoinverse Born operator. We then apply a Radon transform to the results of the first step to calculate the angle-dependent response. Finally, we simultaneously invert for the physical parameters using a weighted least-squares method. Through numerical experiments, we first illustrate the consequences of acoustic approximation on elastic data, leading to inaccurate parameter inversion as well as to artificial reflector inclusion. Then we demonstrate that our method can simultaneously invert for elastic parameters in the presence of complex uncorrelated structures, inaccurate background models, and Gaussian noisy data.

Imaging conditions for elastic reverse time migration

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. S95-S111 ◽  
Author(s):  
Wei Zhang ◽  
Ying Shi
Keyword(s):  
Decomposition Method ◽  
Inner Product ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Imaging Condition ◽  
Subsurface Structures ◽  
Time Migration ◽  
Imaging Conditions ◽  
Wave Modes

Elastic reverse time migration (RTM) has the ability to retrieve accurately migrated images of complex subsurface structures by imaging the multicomponent seismic data. However, the imaging condition applied in elastic RTM significantly influences the quality of the migrated images. We evaluated three kinds of imaging conditions in elastic RTM. The first kind of imaging condition involves the crosscorrelation between the Cartesian components of the particle-velocity wavefields to yield migrated images of subsurface structures. An alternative crosscorrelation imaging condition between the separated pure wave modes obtained by a Helmholtz-like decomposition method could produce reflectivity images with explicit physical meaning and fewer crosstalk artifacts. A drawback of this approach, though, was that the polarity reversal of the separated S-wave could cause destructive interference in the converted-wave image after stacking over multiple shots. Unlike the conventional decomposition method, the elastic wavefields can also be decomposed in the vector domain using the decoupled elastic wave equation, which preserves the amplitude and phase information of the original elastic wavefields. We have developed an inner-product imaging condition to match the vector-separated P- and S-wave modes to obtain scalar reflectivity images of the subsurface. Moreover, an auxiliary P-wave stress image can supplement the elastic imaging. Using synthetic examples with a layered model, the Marmousi 2 model, and a fault model, we determined that the inner-product imaging condition has prominent advantages over the other two imaging conditions and generates images with preserved amplitude and phase attributes.

Acoustic-elastic coupled equation for ocean bottom seismic data elastic reverse time migration

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. S333-S345 ◽  
Author(s):  
Pengfei Yu ◽  
Jianhua Geng ◽  
Xiaobo Li ◽  
Chenlong Wang
Keyword(s):  
Boundary Conditions ◽  
Particle Velocity ◽  
Ocean Bottom ◽  
Receiver Side ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
New Equation ◽  
Obs Data ◽  
Coupled Equation

Conventionally, multicomponent geophones used to record the elastic wavefields in the solid seabed are necessary for ocean bottom seismic (OBS) data elastic reverse time migration (RTM). Particle velocity components are usually injected directly as boundary conditions in the elastic-wave equation in the receiver-side wavefield extrapolation step, which causes artifacts in the resulting elastic images. We have deduced a first-order acoustic-elastic coupled equation (AECE) by substituting pressure fields into the elastic velocity-stress equation (EVSE). AECE has three advantages for OBS data over EVSE when performing elastic RTM. First, the new equation unifies wave propagation in acoustic and elastic media. Second, the new equation separates P-waves directly during wavefield propagation. Third, three approaches are identified when using the receiver-side multicomponent particle velocity records and pressure records in elastic RTM processing: (1) particle velocity components are set as boundary conditions in receiver-side vectorial extrapolation with the AECE, which is equal to the elastic RTM using the conventional EVSE; (2) the pressure component may also be used for receiver-side scalar extrapolation with the AECE, and with which we can accomplish PP and PS images using only the pressure records and suppress most of the artifacts in the PP image with vectorial extrapolation; and (3) ocean-bottom 4C data can be simultaneously used for elastic images with receiver-side tensorial extrapolation using the AECE. Thus, the AECE may be used for conventional elastic RTM, but it also offers the flexibility to obtain PP and PS images using only pressure records.

Surface-offset gathers from elastic reverse time migration and velocity analysis

Geophysics ◽  
2020 ◽  
Vol 85 (1) ◽  
pp. S47-S64
Author(s):  
Yang Zhao ◽  
Tao Liu ◽  
Xueyi Jia ◽  
Hongwei Liu ◽  
Zhiguang Xue ◽  
...  
Keyword(s):  
Velocity Ratio ◽  
Computational Cost ◽  
Cost Effective ◽  
Velocity Estimation ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Time Migration ◽  
The Common ◽  
Wave Decomposition

Angle-domain common-image gathers (ADCIGs) from elastic reverse time migration (ERTM) are valuable tools for seismic elastic velocity estimation. Traditional ADCIGs are based on the concept of common-offset domains, but common-shot domain implementations are often favored for computational cost considerations. Surface-offset gathers (SOGs) built from common-offset migration may serve as an alternative to the common-shot ADCIGs. We have developed a theoretical kinematic framework between these two domains, and we determined that the common SOG gives an alternative measurement of kinematic correctness in the presence of incorrect velocity. Specifically, we exploit analytical expressions for the image misposition between these two domains, with respect to the traveltime perturbation caused by velocity errors. Four formulations of the PP and PS residual moveout functions are derived and provide insightful information of the velocity error, angle, and PS velocity ratio contained in ERTM gathers. The analytical solutions are validated with homogeneous examples with a series of varied parameters. We found that the SOGs may perform in the way of simplicity and linearity as an alternative to the common-shot migration. To make a full comparison with ADCIGs, we have developed a cost-effective workflow of ERTM SOGs. A fast vector P- and S-wave decomposition can be obtained via spatial gradients at selected time steps. A selected ERTM imaging condition is then modified in which the migration is done by offset groups between each source and receiver pair for each P- and S-wave decomposition. Two synthetic (marine and land) examples are used to demonstrate the feasibility of our methods.

Acoustic and elastic numerical wave simulations by recursive spatial derivative operators

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. T167-T174 ◽  
Author(s):  
Dan Kosloff ◽  
Reynam C. Pestana ◽  
Hillel Tal-Ezer
Keyword(s):  
Synthetic Data ◽  
Analytic Solutions ◽  
Numerical Dispersion ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Lamb’S Problem ◽  
Dynamic Elasticity ◽  
Time Migration ◽  
Recursive Operators

A new scheme for the calculation of spatial derivatives has been developed. The technique is based on recursive derivative operators that are generated by an [Formula: see text] fit in the spectral domain. The use of recursive operators enables us to extend acoustic and elastic wave simulations to shorter wavelengths. The method is applied to the numerical solution of the 2D acoustic wave equation and to the solution of the equations of 2D dynamic elasticity in an isotropic medium. An example of reverse-time migration of a synthetic data set shows that the numerical dispersion can be significantly reduced with respect to schemes that are based on finite differences. The method is tested for the solutions of the equations of dynamic elasticity by comparing numerical and analytic solutions to Lamb’s problem.

Mini-batch Least-Squares Reverse Time Migration In A Deep Learning Framework

Geophysics ◽  
2020 ◽  
pp. 1-61
Author(s):  
Janaki Vamaraju ◽  
Jeremy Vila ◽  
Mauricio Araya-Polo ◽  
Debanjan Datta ◽  
Mohamed Sidahmed ◽  
...  
Keyword(s):  
Deep Learning ◽  
Least Squares ◽  
Spatial Resolution ◽  
Data Science ◽  
Computational Cost ◽  
Computation Time ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Learning Framework ◽  
Time Migration

Migration techniques are an integral part of seismic imaging workflows. Least-squares reverse time migration (LSRTM) overcomes some of the shortcomings of conventional migration algorithms by compensating for illumination and removing sampling artifacts to increase spatial resolution. However, the computational cost associated with iterative LSRTM is high and convergence can be slow in complex media. We implement pre-stack LSRTM in a deep learning framework and adopt strategies from the data science domain to accelerate convergence. The proposed hybrid framework leverages the existing physics-based models and machine learning optimizers to achieve better and cheaper solutions. Using a time-domain formulation, we show that mini-batch gradients can reduce the computation cost by using a subset of total shots for each iteration. Mini-batch approach does not only reduce source cross-talk but also is less memory intensive. Combining mini-batch gradients with deep learning optimizers and loss functions can improve the efficiency of LSRTM. Deep learning optimizers such as the adaptive moment estimation are generally well suited for noisy and sparse data. We compare different optimizers and demonstrate their efficacy in mitigating migration artifacts. To accelerate the inversion, we adopt the regularised Huber loss function in conjunction. We apply these techniques to 2D Marmousi and 3D SEG/EAGE salt models and show improvements over conventional LSRTM baselines. The proposed approach achieves higher spatial resolution in less computation time measured by various qualitative and quantitative evaluation metrics.

Sign in / Sign up

Export Citation Format

Share Document