Coupled equations for reverse time migration in transversely isotropic media

Geophysics ◽  
2010 ◽  
Vol 75 (1) ◽  
pp. S11-S22 ◽  
Author(s):  
Paul J. Fowler ◽  
Xiang Du ◽  
Robin P. Fletcher
Keyword(s):  
Dispersion Relation ◽  
Fourth Order ◽  
Transversely Isotropic ◽  
Shear Velocity ◽  
Second Order ◽  
Vertical Shear ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Coupled Equations ◽  
Time Migration

Reverse time migration (RTM) images reflectors by using time-extrapolation modeling codes to synthesize source and receiver wavefields in the subsurface. Asymptotic analysis of wave propagation in transversely isotropic (TI) media yields a dispersion relation describing coupled P- and SV-wave modes. This dispersion relation can be converted into a fourth-order scalar partial differential equation (PDE). Increased computational efficiency can be achieved using equivalent coupled second-order PDEs. Analysis of the corresponding dispersion relations as matrix eigenvalue systems allows one to characterize all possible coupled linear second-order systems equivalent to a given linear fourth-order PDE and to determine which ones yield optimally efficient finite-difference implementations. Setting the shear velocity along the axis of symmetry to zero yields a simpler approximate TI wave equation that is more efficient to implement. This simpler approximation, however, can become unstable for some plausible combinations of anisotropic parameters. The same eigensystem analysis can be applied using finite vertical shear velocity to obtain solutions that avoid these instability problems.

Reverse time migration in tilted transversely isotropic (TTI) media

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCA179-WCA187 ◽  
Author(s):  
Robin P. Fletcher ◽  
Xiang Du ◽  
Paul J. Fowler
Keyword(s):  
Wave Equation ◽  
Numerical Solutions ◽  
Transversely Isotropic ◽  
Shear Velocity ◽  
Vertical Axis ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Axis Of Symmetry ◽  
Tti Media

Reverse time migration (RTM) exhibits great advantages over other imaging methods because it is based on computing numerical solutions to a two-way wave equation. It does not suffer from dip limitation like one-way downward continuation techniques do, thus enabling overturned reflections to be imaged. As well as correctly handling multipathing, RTM has the potential to image internal multiples when the boundaries responsible for generating the multiples are present in the model. In isotropic media, one can use a scalar acoustic wave equation for RTM of pressure data. In anisotropic media, P- and SV-waves are coupled together so, formally, elastic wave equations must be used for RTM. A new wave equation for P-waves is proposed in tilted transversely isotropic (TTI) media that can be solved as part of an acoustic anisotropic RTM algorithm, using standard explicit finite differencing. If the shear velocity along the axis of symmetry is set to zero, stable numerical solutions can be computed for media with a vertical axis of symmetry and [Formula: see text] not less than [Formula: see text]. In TTI media with rapid variations in the direction of the axis of symmetry, setting the shear velocity along the axis of symmetry to zero can cause numerical solutions to become unstable. A solution to this problem is proposed that involves using a small amount of nonzero shear velocity. The amount of shear velocity added is chosen to remove triplications from the SV wavefront and to minimize the anisotropic term of the SV reflection coefficient. We show modeling and high-quality RTM results in complex TTI media using this equation.

Analysis of higher‐order, finite‐difference schemes in 3-D reverse‐time migration

Geophysics ◽  
1996 ◽  
Vol 61 (3) ◽  
pp. 845-856 ◽  
Author(s):  
Wen‐Jing Wu ◽  
Larry R. Lines ◽  
Han‐Xing Lu
Keyword(s):  
Finite Difference ◽  
Fourth Order ◽  
Higher Order ◽  
Second Order ◽  
Difference Schemes ◽  
Finite Difference Schemes ◽  
Scalar Wave ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration

The design and usefulness of practical algorithms for 3-D reverse‐time depth migration is examined, while demonstrating data applications of the algorithm. We evaluate quantitatively the accuracy of the finite‐difference operator from second‐order to eighth‐order for the scalar wave equation by comparing numerical and analytical solutions. The results clearly show the advantage of using higher‐order, finite‐difference schemes, especially from second‐order to fourth‐order for space derivatives. Hence, a finite‐difference method with the accuracy of fourth‐order in space and second‐order in time is applied to 3-D full scalar wave equations in reverse‐time migration. Considerable savings in CPU and memory are obtained by using larger horizontal grid spacings than the vertical spacings. Therefore, we derive dispersion and stability conditions for such unequal grid spacing. The 3-D reverse‐time migration of data from the Hibernia Field shows an image improvement over 2-D migrations, despite the fact that much of the structure was considered to be approximately 2-D. Stable, nondispersive, finite‐difference methods of higher order can allow for tractable and efficient 3-D reverse‐time migration solutions.

Decoupled equations for reverse time migration in tilted transversely isotropic media

Geophysics ◽  
2012 ◽  
Vol 77 (2) ◽  
pp. T37-T45 ◽  
Author(s):  
Ge Zhan ◽  
Reynam C. Pestana ◽  
Paul L. Stoffa
Keyword(s):  
Broad Band ◽  
Transversely Isotropic ◽  
Rapid Expansion ◽  
P Wave ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Coupled Equations ◽  
Time Migration ◽  
Tti Media ◽  
Wave Modes

Conventional modeling and migration for tilted transversely isotropic (TTI) media may suffer from numerical instabilities and shear wave artifacts due to the coupling of the P-wave and SV-wave modes in the TTI coupled equations. Starting with the separated P- and SV-phase velocity expressions for vertical transversely isotropic (VTI) media, we extend these decoupled equations for modeling and reverse time migration (RTM) in acoustic TTI media. Compared with the TTI coupled equations published in the geophysical literature, the new TTI decoupled equations provide a more stable solution due to the complete separation of the P-wave and SV-wave modes. The pseudospectral method is the most convenient method to implement these equations due to the form of wavenumber expressions and has the added benefit of being highly accurate and thus avoiding numerical dispersion. The rapid expansion method (REM) in time is employed to produce a broad band numerically stable time evolution of the wavefields. Synthetic results validate the proposed TTI decoupled equations and show that modeling and RTM in TTI media with the decoupled equations remain numerically stable even for models with strong anisotropy and sharp contrasts.

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.

Q least-squares reverse time migration based on the first-order viscoacoustic quasidifferential equations

Geophysics ◽  
2021 ◽  
pp. 1-65
Author(s):  
Yingming Qu ◽  
Yixin Wang ◽  
Zhenchun Li ◽  
Chang Liu
Keyword(s):  
Wave Equation ◽  
Least Squares ◽  
Surface Topography ◽  
Density Perturbation ◽  
Second Order ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
First Order ◽  
Time Migration ◽  
Grid Scheme

Seismic wave attenuation caused by subsurface viscoelasticity reduces the quality of migration and the reliability of interpretation. A variety of Q-compensated migration methods have been developed based on the second-order viscoacoustic quasidifferential equations. However, these second-order wave-equation-based methods are difficult to handle with density perturbation and surface topography. In addition, the staggered grid scheme, which has an advantage over the collocated grid scheme because of its reduced numerical dispersion and enhanced stability, works in first-order wave-equation-based methods. We have developed a Q least-squares reverse time migration method based on the first-order viscoacoustic quasidifferential equations by deriving Q-compensated forward-propagated operators, Q-compensated adjoint operators, and Q-attenuated Born modeling operators. Besides, our method using curvilinear grids is available even when the attenuating medium has surface topography and can conduct Q-compensated migration with density perturbation. The results of numerical tests on two synthetic and a field data sets indicate that our method improves the imaging quality with iterations and produces better imaging results with clearer structures, higher signal-to-noise ratio, higher resolution, and more balanced amplitude by correcting the energy loss and phase distortion caused by Q attenuation. It also suppresses the scattering and diffracted noise caused by the surface topography.

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.

Separating quasi-P-wave in transversely isotropic media with a vertical symmetry axis by synthesized pressure applied to ocean-bottom seismic data elastic reverse time migration

Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. C295-C307 ◽  
Author(s):  
Pengfei Yu ◽  
Jianhua Geng ◽  
Chenlong Wang
Keyword(s):  
Transversely Isotropic ◽  
P Wave ◽  
Ocean Bottom ◽  
Receiver Side ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Isotropic Media ◽  
Wave Imaging ◽  
Obs Data

Quasi-P (qP)-wavefield separation is a crucial step for elastic P-wave imaging in anisotropic media. It is, however, notoriously challenging to quickly and accurately obtain separated qP-wavefields. Based on the concepts of the trace of the stress tensor and the pressure fields defined in isotropic media, we have developed a new method to rapidly separate the qP-wave in a transversely isotropic medium with a vertical symmetry axis (VTI) by synthesized pressure from ocean-bottom seismic (OBS) data as a preprocessing step for elastic reverse time migration (ERTM). Another key aspect of OBS data elastic wave imaging is receiver-side 4C records back extrapolation. Recent studies have revealed that receiver-side tensorial extrapolation in isotropic media with ocean-bottom 4C records can sufficiently suppress nonphysical waves produced during receiver-side reverse time wavefield extrapolation. Similarly, the receiver-side 4C records tensorial extrapolation was extended to ERTM in VTI media in our studies. Combining a separated qP-wave by synthesizing pressure and receiver-side wavefield reverse time tensorial extrapolation with the crosscorrelation imaging condition, we have developed a robust, fast, flexible, and elastic imaging quality improved method in VTI media for OBS data.

Sign in / Sign up

Export Citation Format

Share Document