Reverse time migration angle gathers using Poynting vector and pseudospectral method

Author(s):  
Kwangjin Yoon
2013 ◽  
Author(s):  
Edvaldo S. Araujo ◽  
Reynam C. Pestana ◽  
Adriano W. G. dos Santos

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. S359-S376 ◽  
Author(s):  
Chen Tang ◽  
George A. McMechan

Because receiver wavefields reconstructed from observed data are not as stable as synthetic source wavefields, the source-propagation vector and the reflector normal have often been used to calculate angle-domain common-image gathers (ADCIGs) from reverse time migration. However, the existing data flows have three main limitations: (1) Calculating the propagation direction only at the wavefields with maximum amplitudes ignores multiarrivals; using the crosscorrelation imaging condition at each time step can include the multiarrivals but will result in backscattering artifacts. (2) Neither amplitude picking nor Poynting-vector calculations are accurate for overlapping wavefields. (3) Calculating the reflector normal in space is not accurate for a structurally complicated reflection image, and calculating it in the wavenumber ([Formula: see text]) domain may give Fourier truncation artifacts. We address these three limitations in an improved data flow with two steps: During imaging, we use a multidirectional Poynting vector (MPV) to calculate the propagation vectors of the source wavefield at each time step and output intermediate source-angle-domain CIGs (SACIGs). After imaging, we use an antitruncation-artifact Fourier transform (ATFT) to convert SACIGs to ADCIGs in the [Formula: see text]-domain. To achieve the new flow, another three innovative aspects are included. In the first step, we develop an angle-tapering scheme to remove the Fourier truncation artifacts during the wave decomposition (of MPV) while preserving the amplitudes, and we use a wavefield decomposition plus angle-filter imaging condition to remove the backscattering artifacts in the SACIGs. In the second step, we compare two algorithms to remove the Fourier truncation artifacts that are caused by the plane-wave assumption. One uses an antileakage FT (ALFT) in local windows; the other uses an antitruncation-artifact FT, which relaxes the plane-wave assumption and thus can be done for the global space. The second algorithm is preferred. Numerical tests indicate that this new flow (source-side MPV plus ATFT) gives high-quality ADCIGs.


2016 ◽  
Vol 24 (02) ◽  
pp. 1650005 ◽  
Author(s):  
Jiangang Xie ◽  
Zichao Guo ◽  
Hai Liu ◽  
Qing Huo Liu

We propose a pre-stack reverse time migration (RTM) seismic imaging method using the pseudospectral time-domain (PSTD) algorithm. Traditional pseudospectral method uses the fast Fourier transform (FFT) algorithm to calculate the spatial derivatives, but is limited by the wraparound effect due to the periodicity assumed in the FFT. The PSTD algorithm combines the pseudospectral method with a perfectly matched layer (PML) for acoustic waves. PML is a highly effective absorbing boundary condition that can eliminate the wraparound effect. It enables a wide application of the pseudospectral method to complex models. RTM based on the PSTD algorithm has advantages in the computational efficiency compared to traditional methods such as the second-order and high order finite difference time-domain (FDTD) methods. In this work, we implement the PSTD algorithm for acoustic wave equation based RTM. By applying the PSTD-RTM method to various seismic models and comparing it with RTM based on the eighth-order FDTD method, we find that PSTD-RTM method has better performance and saves more than 50% memory. The method is suitable for parallel computation, and has been accelerated by general purpose graphics processing unit.


Geophysics ◽  
2009 ◽  
Vol 74 (3) ◽  
pp. S57-S66 ◽  
Author(s):  
J. C. Costa ◽  
F. A. Silva Neto ◽  
M. R. Alcântara ◽  
J. Schleicher ◽  
A. Novais

The quality of seismic images obtained by reverse time migration (RTM) strongly depends on the imaging condition. We propose a new imaging condition that is motivated by stationary phase analysis of the classical crosscorrelation imaging condition. Its implementation requires the Poynting vector of the source and receiver wavefields at the imaging point. An obliquity correction is added to compensate for the reflector dip effect on amplitudes of RTM. Numerical experiments show that using an imaging condition with obliquity compensation improves reverse time migration by reducing backscattering artifacts and improving illumination compensation.


2013 ◽  
Author(s):  
Edvaldo S. Araujo ◽  
Reynam C. Pestana ◽  
Adriano W. G. dos Santos

Geophysics ◽  
2014 ◽  
Vol 79 (5) ◽  
pp. S163-S172 ◽  
Author(s):  
Edvaldo S. Araujo ◽  
Reynam C. Pestana ◽  
Adriano W. G. dos Santos

2006 ◽  
Vol 37 (1) ◽  
pp. 102-107 ◽  
Author(s):  
Kwangjin Yoon ◽  
Kurt J. Marfurt

Sign in / Sign up

Export Citation Format

Share Document