Reverse time migration

Geophysics ◽  
1983 ◽  
Vol 48 (11) ◽  
pp. 1514-1524 ◽  
Author(s):  
Edip Baysal ◽  
Dan D. Kosloff ◽  
John W. C. Sherwood
Keyword(s):  
Wave Amplitude ◽  
Synthetic Data ◽  
Data Sets ◽  
Field Observations ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Migration Method ◽  
Zero Offset ◽  
Time Extrapolation

Migration of stacked or zero‐offset sections is based on deriving the wave amplitude in space from wave field observations at the surface. Conventionally this calculation has been carried out through a depth extrapolation. We examine the alternative of carrying out the migration through a reverse time extrapolation. This approach may offer improvements over existing migration methods, especially in cases of steeply dipping structures with strong velocity contrasts. This migration method is tested using appropriate synthetic data sets.

An exact adjoint operation pair in time extrapolation and its application in least-squares reverse-time migration

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. H27-H33 ◽  
Author(s):  
Jun Ji
Keyword(s):  
Least Squares ◽  
Adjoint Operator ◽  
Synthetic Data ◽  
Forward Modeling ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Product Test ◽  
Least Squares Migration ◽  
Time Extrapolation

To reduce the migration artifacts arising from incomplete data or inaccurate operators instead of migrating data with the adjoint of the forward-modeling operator, a least-squares migration often is considered. Least-squares migration requires a forward-modeling operator and its adjoint. In a derivation of the mathematically correct adjoint operator to a given forward-time-extrapolation modeling operator, the exact adjoint of the derived operator is obtained by formulating an explicit matrix equation for the forward operation and transposing it. The programs that implement the exact adjoint operator pair are verified by the dot-product test. The derived exact adjoint operator turns out to differ from the conventional reverse-time-migration (RTM) operator, an implementation of wavefield extrapolation backward in time. Examples with synthetic data show that migration using the exact adjoint operator gives similar results for a conventional RTM operator and that least-squares RTM is quite successful in reducing most migration artifacts. The least-squares solution using the exact adjoint pair produces a model that fits the data better than one using a conventional RTM operator pair.

Fast 3D target-oriented reverse-time datuming

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCA141-WCA151 ◽  
Author(s):  
Shuqian Dong ◽  
Yi Luo ◽  
Xiang Xiao ◽  
Sergio Chávez-Pérez ◽  
Gerard T. Schuster
Keyword(s):  
Synthetic Data ◽  
Data Sets ◽  
Target Area ◽  
Complex Structures ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Salt Domes ◽  
Time Migration ◽  
Subsalt Imaging

Imaging of subsalt sediments is a challenge for traditional migration methods such as Kirchhoff and one-way wave-equation migration. Consequently, the more accurate two-way method of reverse-time migration (RTM) is preferred for subsalt imaging, but its use can be limited by high computation cost. To overcome this problem, a 3D target-oriented reverse-time datuming (RTD) method is presented, which can generate redatumed data economically in target areas beneath complex structures such as salt domes. The redatumed data in the target area then can be migrated inexpensively using a traditional migration method. If the target area is much smaller than the acquisition area, computation costs are reduced significantly by the use of a novel bottom-up strategy to calculate the extrapolated Green’s functions. Target-oriented RTD is tested on 2D and 3D SEG/EAGE synthetic data sets and a 3D field data set from the Gulf of Mexico. Results show that target-oriented RTD combined with standard migration can image sediments beneath complex structures accurately with much less calculation effort than full volume RTM. The requirement is that the area over the target zone is smaller than that of the acquisition survey.

Reverse time migration in tilted transversely isotropic media

2020 ◽  
Vol 38 (2) ◽  
Author(s):  
Razec Cezar Sampaio Pinto da Silva Torres ◽  
Leandro Di Bartolo
Keyword(s):  
Seismic Anisotropy ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Computer Clusters ◽  
Time Migration ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Tti Media

ABSTRACT. Reverse time migration (RTM) is one of the most powerful methods used to generate images of the subsurface. The RTM was proposed in the early 1980s, but only recently it has been routinely used in exploratory projects involving complex geology – Brazilian pre-salt, for example. Because the method uses the two-way wave equation, RTM is able to correctly image any kind of geological environment (simple or complex), including those with anisotropy. On the other hand, RTM is computationally expensive and requires the use of computer clusters. This paper proposes to investigate the influence of anisotropy on seismic imaging through the application of RTM for tilted transversely isotropic (TTI) media in pre-stack synthetic data. This work presents in detail how to implement RTM for TTI media, addressing the main issues and specific details, e.g., the computational resources required. A couple of simple models results are presented, including the application to a BP TTI 2007 benchmark model.Keywords: finite differences, wave numerical modeling, seismic anisotropy. Migração reversa no tempo em meios transversalmente isotrópicos inclinadosRESUMO. A migração reversa no tempo (RTM) é um dos mais poderosos métodos utilizados para gerar imagens da subsuperfície. A RTM foi proposta no início da década de 80, mas apenas recentemente tem sido rotineiramente utilizada em projetos exploratórios envolvendo geologia complexa, em especial no pré-sal brasileiro. Por ser um método que utiliza a equação completa da onda, qualquer configuração do meio geológico pode ser corretamente tratada, em especial na presença de anisotropia. Por outro lado, a RTM é dispendiosa computacionalmente e requer o uso de clusters de computadores por parte da indústria. Este artigo apresenta em detalhes uma implementação da RTM para meios transversalmente isotrópicos inclinados (TTI), abordando as principais dificuldades na sua implementação, além dos recursos computacionais exigidos. O algoritmo desenvolvido é aplicado a casos simples e a um benchmark padrão, conhecido como BP TTI 2007.Palavras-chave: diferenças finitas, modelagem numérica de ondas, anisotropia sísmica.

Least-squares reverse time migration via linearized waveform inversion using a Wasserstein metric

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. S411-S423
Author(s):  
Peng Yong ◽  
Jianping Huang ◽  
Zhenchun Li ◽  
Wenyuan Liao ◽  
Luping Qu
Keyword(s):  
Least Squares ◽  
Waveform Inversion ◽  
Synthetic Data ◽  
Gradient Algorithm ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Wasserstein Metric ◽  
Time Migration ◽  
Imaging Results ◽  
Primal Dual

Least-squares reverse time migration (LSRTM), an effective tool for imaging the structures of the earth from seismograms, can be characterized as a linearized waveform inversion problem. We have investigated the performance of three minimization functionals as the [Formula: see text] norm, the hybrid [Formula: see text] norm, and the Wasserstein metric ([Formula: see text] metric) for LSRTM. The [Formula: see text] metric used in this study is based on the dynamic formulation of transport problems, and a primal-dual hybrid gradient algorithm is introduced to efficiently compute the [Formula: see text] metric between two seismograms. One-dimensional signal analysis has demonstrated that the [Formula: see text] metric behaves like the [Formula: see text] norm for two amplitude-varied signals. Unlike the [Formula: see text] norm, the [Formula: see text] metric does not suffer from the differentiability issue for null residuals. Numerical examples of the application of three misfit functions to LSRTM on synthetic data have demonstrated that, compared to the [Formula: see text] norm, the hybrid [Formula: see text] norm and [Formula: see text] metric can accelerate LSRTM and are less sensitive to non-Gaussian noise. For the field data application, the [Formula: see text] metric produces the most reliable imaging results. The hybrid [Formula: see text] norm requires tedious trial-and-error tests for the judicious threshold parameter selection. Hence, the more automatic [Formula: see text] metric is recommended as a robust alternative to the customary [Formula: see text] norm for time-domain LSRTM.

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.

3-D prestack migration in anisotropic media

Geophysics ◽  
1993 ◽  
Vol 58 (1) ◽  
pp. 79-90 ◽  
Author(s):  
Zhengxin Dong ◽  
George A. McMechan
Keyword(s):  
A Priori ◽  
Three Dimensional ◽  
Synthetic Data ◽  
Anisotropic Media ◽  
P Wave ◽  
Common Source ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Prestack Migration ◽  
Time Migration

A three‐dimensional (3-D) prestack reverse‐time migration algorithm for common‐source P‐wave data from anisotropic media is developed and illustrated by application to synthetic data. Both extrapolation of the data and computation of the excitation‐time imaging condition are implemented using a second‐order finite‐ difference solution of the 3-D anisotropic scalar‐wave equation. Poorly focused, distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation; well focused, clear images are obtained using anisotropic extrapolation. A priori estimation of the 3-D anisotropic velocity distribution is required. Zones of anomalous, directionally dependent reflectivity associated with anisotropic fracture zones are detectable in both the 3-D common‐ source data and the corresponding migrated images.

Least-squares reverse time migration in TTI media using a pure qP-wave equation

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. S199-S216
Author(s):  
Xinru Mu ◽  
Jianping Huang ◽  
Jidong Yang ◽  
Xu Guo ◽  
Yundong Guo
Keyword(s):  
Wave Equation ◽  
Least Squares ◽  
Seismic Imaging ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
The Matrix ◽  
The Stability

Anisotropy is a common phenomenon in subsurface strata and should be considered in seismic imaging and inversion. Seismic imaging in a vertical transversely isotropic (VTI) medium does not take into account the effects of the tilt angles, which can lead to degraded migrated images in areas with strong anisotropy. To correct such waveform distortion, reduce related image artifacts, and improve migration resolution, a tilted transversely isotropic (TTI) least-squares reverse time migration (LSRTM) method is presented. In the LSRTM, a pure qP-wave equation is used and solved with the finite-difference method. We have analyzed the stability condition for the pure qP-wave equation using the matrix method, which is used to ensure the stability of wave propagation in the TTI medium. Based on this wave equation, we derive a corresponding demigration (Born modeling) and adjoint migration operators to implement TTI LSRTM. Numerical tests on the synthetic data show the advantages of TTI LSRTM over VTI RTM and VTI LSRTM when the recorded data contain strong effects caused by large tilt angles. Our numerical experiments illustrate that the sensitivity of the adopted TTI LSRTM to the migration velocity errors is much higher than that to the anisotropic parameters (including epsilon, delta, and tilted angle parameters), and its sensitivity to the epsilon model and tilt angle is higher than that to the delta model.

Sign in / Sign up

Export Citation Format

Share Document