scholarly journals Alternative stable qP wave equations in TTI media with their applications for reverse time migration

2015 ◽  
Vol 12 (5) ◽  
pp. 734-744 ◽  
Author(s):  
Yang Zhou ◽  
Huazhong Wang ◽  
Wenqing Liu
Keyword(s):  
Wave Equations ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Tti Media

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.

Visco-acoustic least-squares reverse time migration in TTI media and application to OBN data

2019 ◽  
Author(s):  
Side Jin ◽  
Henning Kuehl ◽  
Michael Kiehn ◽  
Rene-Edward Plessix ◽  
Martina Wittmann-Hohlbein
Keyword(s):  
Least Squares ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Tti Media

Elastic reverse time migration using acoustic propagators

Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. S399-S408 ◽  
Author(s):  
Yunyue Elita Li ◽  
Yue Du ◽  
Jizhong Yang ◽  
Arthur Cheng ◽  
Xinding Fang
Keyword(s):  
Wave Equations ◽  
Mode Conversion ◽  
Computational Cost ◽  
Body Waves ◽  
Constant Density ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Isotropic Constant ◽  
Time Migration

Elastic wave imaging has been a significant challenge in the exploration industry due to the complexities in wave physics and numerical implementation. We have separated the governing equations for P- and S-wave propagation without the assumptions of homogeneous Lamé parameters to capture the mode conversion between the two body waves in an isotropic, constant-density medium. The resulting set of two coupled second-order equations for P- and S-potentials clearly demonstrates that mode conversion only occurs at the discontinuities of the shear modulus. Applying the Born approximation to the new equations, we derive the PP, PS, SP, and SS imaging conditions from the first gradients of waveform matching objective functions. The resulting images are consistent with the physical perturbations of the elastic parameters, and, hence, they are automatically free of the polarity reversal artifacts in the converted images. When implementing elastic reverse time migration (RTM), we find that scalar wave equations can be used to back propagate the recorded P-potential, as well as individual components in the vector field of the S-potential. Compared with conventional elastic RTM, the proposed elastic RTM implementation using acoustic propagators not only simplifies the imaging condition, it but also reduces the computational cost and the artifacts in the images. We have determined the accuracy of our method using 2D and 3D numerical examples.

True-amplitude linearized waveform inversion with the quasi-elastic wave equation

Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. R827-R844 ◽  
Author(s):  
Zongcai Feng ◽  
Gerard Schuster
Keyword(s):  
Wave Equation ◽  
Elastic Wave ◽  
Wave Equations ◽  
Waveform Inversion ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Elastic Wave Equation ◽  
P Waves ◽  
Time Migration

We present a quasi-elastic wave equation as a function of the pressure variable, which can accurately model PP reflections with elastic amplitude variation with offset effects under the first-order Born approximation. The kinematic part of the quasi-elastic wave equation accurately models the propagation of P waves, whereas the virtual-source part, which models the amplitudes of reflections, is a function of the perturbations of density and Lamé parameters [Formula: see text] and [Formula: see text]. The quasi-elastic wave equation generates a scattering radiation pattern that is exactly the same as that for the elastic wave equation, and only requires the solution of two acoustic wave equations for each shot gather. This means that the quasi-elastic wave equation can be used for true-amplitude linearized waveform inversion (also known as least-squares reverse time migration) of elastic PP reflections, in which the corresponding misfit gradients are with respect to the perturbations of density and the P- and S-wave impedances. The perturbations of elastic parameters are iteratively updated by minimizing the [Formula: see text]-norm of the difference between the recorded PP reflections and the predicted pressure data modeled from the quasi-elastic wave equation. Numerical tests on synthetic and field data indicate that true-amplitude linearized waveform inversion using the quasi-elastic wave equation can account for the elastic PP amplitudes and provide a robust estimate of the perturbations of P- and S-wave impedances and, in some cases, the density. In addition, true-amplitude linearized waveform inversion provides images with a wider bandwidth and fewer artifacts because the PP amplitudes are accurately explained. We also determine the 2D scalar quasi-elastic wave equation for P-SV reflections and the 3D vector equation for PS reflections.

High-order finite-difference approximations to solve pseudoacoustic equations in 3D VTI media

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. T225-T235 ◽  
Author(s):  
Leandro Di Bartolo ◽  
Leandro Lopes ◽  
Luis Juracy Rangel Lemos
Keyword(s):  
Finite Difference ◽  
Wave Equations ◽  
Computational Cost ◽  
Transversely Isotropic ◽  
High Order ◽  
Staggered Grid ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Grid Scheme

Pseudoacoustic algorithms are very fast in comparison with full elastic ones for vertical transversely isotropic (VTI) modeling, so they are suitable for many applications, especially reverse time migration. Finite differences using simple grids are commonly used to solve pseudoacoustic equations. We have developed and implemented general high-order 3D pseudoacoustic transversely isotropic formulations. The focus is the development of staggered-grid finite-difference algorithms, known for their superior numerical properties. The staggered-grid schemes based on first-order velocity-stress wave equations are developed in detail as well as schemes based on direct application to second-order stress equations. This last case uses the recently presented equivalent staggered-grid theory, resulting in a staggered-grid scheme that overcomes the problem of large memory requirement. Two examples are presented: a 3D simulation and a prestack reverse time migration application, and we perform a numerical analysis regarding computational cost and precision. The errors of the new schemes are smaller than the existing nonstaggered-grid schemes. In comparison with existing staggered-grid schemes, they require 25% less memory and only have slightly greater computational cost.

Sign in / Sign up

Export Citation Format

Share Document