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-compensated reverse time migration in tilted transversely isotropic media

Geophysics ◽  
2020 ◽  
pp. 1-79
Author(s):  
Ali Fathalian ◽  
Daniel O. Trad ◽  
Kristopher A. Innanen
Keyword(s):  
Time Domain ◽  
Transversely Isotropic ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Modeling And Analysis ◽  
Seismic Amplitude ◽  
Time Migration ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
The Time Domain

Anisotropy and absorption are critical to the modeling and analysis of seismic amplitude,phase, and traveltime data. To neglect any of these phenomena, which are often bothoperating simultaneously, degrades the resolution and interpretability of migrated images.However, a full accounting of anisotropy and anelasticity is computationally complex andexpensive. One strategy for accommodating these aspects of wave propagation, while keepingcost and complexity under control, is to do so within an acoustic approximation. Weset up a procedure for solving the time-domain viscoacoustic wave equation for tilted transverselyisotropic (TTI) media, based on a standard linear solid model and, from this, developa viscoacoustic reverse time migration (Q-RTM) algorithm. In this approach, amplitudecompensation occurs within the migration process through a manipulation of attenuationand phase dispersion terms in the time domain differential equations. Specifically, theback-propagation operator is constructed by reversing the sign only of the amplitude lossoperators, but not the dispersion-related operators, a step made possible by reformulatingthe absorptive TTI equations such that the loss and dispersion operators appear separately.The scheme is tested on synthetic examples to examine the capacity of viscoacoustic RTM to correct for attenuation, and the overall stability of the procedure.

Scalar and vector imaging based on wave mode decoupling for elastic reverse time migration in isotropic and transversely isotropic media

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. S383-S398 ◽  
Author(s):  
Chenlong Wang ◽  
Jiubing Cheng ◽  
Børge Arntsen
Keyword(s):  
Computational Cost ◽  
Transversely Isotropic ◽  
Polarity Reversal ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Time Migration ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Wave Modes

Recording P- and S-wave modes acquires more information related to rock properties of the earth’s interior. Elastic migration, as a part of multicomponent seismic data processing, potentially offers a great improvement over conventional acoustic migration to create a spatial image of some medium properties. In the framework of elastic reverse time migration, we have developed new scalar and vector imaging conditions assisted by efficient polarization-based mode decoupling to avoid crosstalk among the different wave modes for isotropic and transversely isotropic media. For the scalar imaging, we corrected polarity reversal of zero-lag PS images using the local angular attributes on the fly of angle-domain imaging. For the vector imaging, we naturally used the polarization information in the decoupled single-mode vector fields to automatically avoid the polarity reversal and to estimate the local angular attributes for angle-domain imaging. Examples of increasing complexity in 2D and 3D cases found that the proposed approaches can be used to obtain a physically interpretable image and angle-domain common-image gather at an acceptable computational cost. Decoupling and imaging the 3D S-waves involves some complexity, which has not been addressed in the literature. For this reason, we also attempted at illustrating the physical contents of the two separated S-wave modes and their contribution to seismic full-wave imaging.

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.

Anisotropic migration velocity analysis using reverse-time migration

Geophysics ◽  
2014 ◽  
Vol 79 (1) ◽  
pp. R13-R25 ◽  
Author(s):  
Wiktor Waldemar Weibull ◽  
Børge Arntsen
Keyword(s):  
Wave Equation ◽  
Seismic Anisotropy ◽  
Migration Velocity ◽  
Image Point ◽  
Velocity Analysis ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Isotropic Media

Seismic anisotropy, if not accounted for, can cause significant mispositioning of the reflectors in depth-migrated images. Accounting for anisotropy in depth migration requires velocity analysis tools that can estimate the anisotropic background velocity field. We extended wave equation migration velocity analysis to deal with 2D tilted transverse isotropic media. The velocities were obtained automatically by nonlinear optimization of the focusing and stack power of common-image point gathers constructed using an extended imaging condition. We used the elastic two-way wave equation to reconstruct the wavefields needed for the image and gradient computations. This led to an anisotropic migration velocity analysis algorithm based on reverse-time migration. We illustrated the method with synthetic and field data examples based on marine surface seismic acquisition. The results showed that the method significantly improves the quality of the depth-migrated image. However, as is common in the case of velocity analysis using surface seismic data, the estimation of anisotropic parameters seems to be strongly nonunique.

Elastic least-squares reverse time migration in vertical transverse isotropic media

Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. S539-S553 ◽  
Author(s):  
Jidong Yang ◽  
Hejun Zhu ◽  
George McMechan ◽  
Houzhu Zhang ◽  
Yang Zhao
Keyword(s):  
Wave Equation ◽  
Least Squares ◽  
Transversely Isotropic ◽  
Image Resolution ◽  
Total Variation Regularization ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Isotropic Media ◽  
Band Limited

Using adjoint-based elastic reverse time migration, it is difficult to produce high-quality reflectivity images due to the limited acquisition apertures, band-limited source time function, and irregular subsurface illumination. Through iteratively computing the Hessian inverse, least-squares migration enables us to reduce the point-spread-function effects and improve the image resolution and amplitude fidelity. By incorporating anisotropy in the 2D elastic wave equation, we have developed an elastic least-squares reverse time migration (LSRTM) method for multicomponent data from the vertically transversely isotropic (VTI) media. Using the perturbed stiffness parameters [Formula: see text] and [Formula: see text] as PP and PS reflectivities, we linearize the elastic VTI wave equation and obtain a Born modeling (demigration) operator. Then, we use the Lagrange multiplier method to derive the corresponding adjoint wave equation and reflectivity kernels. With linearized forward modeling and adjoint migration operators, we solve a linear inverse problem to estimate the subsurface reflectivity models for [Formula: see text] and [Formula: see text]. To reduce the artifacts caused by data over-fitting, we introduce total-variation regularization into the reflectivity inversion, which promotes a sparse solution in terms of the model derivatives. To accelerate the convergence of LSRTM, we use source illumination to approximate the diagonal Hessian and use it as a preconditioner for the misfit gradient. Numerical examples help us determine that our elastic VTI LSRTM method can improve the spatial resolution and amplitude fidelity in comparison to adjoint migration.

Sign in / Sign up

Export Citation Format

Share Document