Angle‐domain common‐image gathers in reverse‐time migration by combining the Poynting vector with local‐wavefield decomposition

2019 ◽  
Vol 67 (8) ◽  
pp. 2035-2060 ◽  
Author(s):  
Wu Chengliang ◽  
Wang Huazhong ◽  
Zhou Yang ◽  
Hu Jiangtao
Keyword(s):  
Poynting Vector ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Wavefield Decomposition

scholarly journals Angle-Weighted Reverse Time Migration With Wavefield Decomposition Based on the Optical Flow Vector

2021 ◽  
Vol 9 ◽  
Author(s):  
Chuang Xie ◽  
Peng Song ◽  
Xishuang Li ◽  
Jun Tan ◽  
Shaowen Wang ◽  
...  
Keyword(s):  
Optical Flow ◽  
Poynting Vector ◽  
Large Angle ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Reflection Angle ◽  
Flow Vector ◽  
Imaging Point ◽  
Wavefield Decomposition

Reverse time migration (RTM) is based on the two-way wave equation, so its imaging results obtained by conventional zero-lag cross-correlation imaging conditions contain a lot of low-wavenumber noises. So far, the wavefield decomposition method based on the Poynting vector has been developed to suppress these noises; however, this method also has some problems, such as unstable calculation of the Poynting vector, low accuracy of wavefield decomposition, and poor effect of large-angle migration artifacts suppression. This article introduces the optical flow vector method to RTM to realize high-precision wavefield decomposition for both the source and receiver wavefields and obtains four directions of wavefields: up-, down-, left-, and right-going. Then, the cross-correlation imaging sections of one-way propagation components of forward- and back-propagated wavefields are optimized and stacked. On this basis, the reflection angle of each imaging point is calculated based on the optical flow vector, and an attenuation factor related to the reflection angle is introduced as the weight to generate the optimal stack images. The tests of theoretical model and field marine seismic data illustrate that compared with the conventional RTM with wavefield decomposition based on the Poynting vector, the angle-weighted RTM with wavefield decomposition based on the optical flow vector proposed in this article can achieve wavefield decomposition for both the source and receiver wavefields and calculate the reflection angle of each imaging point more accurately and stably. Moreover, the proposed method adopts angle weighting processing, which can further eliminate large-angle migration artifacts and effectively improve the imaging accuracy of RTM.

Reducing artifacts of elastic reverse time migration by the deprimary technique

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. S569-S577 ◽  
Author(s):  
Yang Zhao ◽  
Houzhu Zhang ◽  
Jidong Yang ◽  
Tong Fei
Keyword(s):  
Complex Function ◽  
Noise Removal ◽  
Pure Mode ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Imaging Condition ◽  
Velocity Models ◽  
Time Migration ◽  
Wavefield Decomposition

Using the two-way elastic-wave equation, elastic reverse time migration (ERTM) is superior to acoustic RTM because ERTM can handle mode conversions and S-wave propagations in complex realistic subsurface. However, ERTM results may not only contain classical backscattering noises, but they may also suffer from false images associated with primary P- and S-wave reflections along their nonphysical paths. These false images are produced by specific wave paths in migration velocity models in the presence of sharp interfaces or strong velocity contrasts. We have addressed these issues explicitly by introducing a primary noise removal strategy into ERTM, in which the up- and downgoing waves are efficiently separated from the pure-mode vector P- and S-wavefields during source- and receiver-side wavefield extrapolation. Specifically, we investigate a new method of vector wavefield decomposition, which allows us to produce the same phases and amplitudes for the separated P- and S-wavefields as those of the input elastic wavefields. A complex function involved with the Hilbert transform is used in up- and downgoing wavefield decomposition. Our approach is cost effective and avoids the large storage of wavefield snapshots that is required by the conventional wavefield separation technique. A modified dot-product imaging condition is proposed to produce multicomponent PP-, PS-, SP-, and SS-images. We apply our imaging condition to two synthetic models, and we demonstrate the improvement on the image quality of ERTM.

Symplectic scheme and the Poynting vector in the reverse time migration

2013 ◽  
Author(s):  
Edvaldo S. Araujo ◽  
Reynam C. Pestana ◽  
Adriano W. G. dos Santos
Keyword(s):  
Poynting Vector ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Symplectic Scheme

Combining multidirectional source vector with antitruncation-artifact Fourier transform to calculate angle gathers from reverse time migration in two steps

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. S359-S376 ◽  
Author(s):  
Chen Tang ◽  
George A. McMechan
Keyword(s):  
Fourier Transform ◽  
Plane Wave ◽  
Poynting Vector ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Step ◽  
Imaging Condition ◽  
Time Migration ◽  
Reflection Image ◽  
Wave Decomposition

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.

scholarly journals A causal imaging condition for reverse time migration using the Discrete Hilbert transform and its efficient implementation on GPU

2019 ◽  
Vol 16 (5) ◽  
pp. 894-912
Author(s):  
Feipeng Li ◽  
Jinghuai Gao ◽  
Zhaoqi Gao ◽  
Xiudi Jiang ◽  
Wenbo Sun
Keyword(s):  
Image Quality ◽  
Hilbert Transform ◽  
Real Data ◽  
Processing Unit ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Discrete Hilbert Transform ◽  
Imaging Condition ◽  
Time Migration ◽  
Wavefield Decomposition

Abstract Reverse time migration (RTM) has shown a significant advantage over other imaging algorithms for imaging complex subsurface structures. However, low-wavenumber noise severely contaminates the image, which is one of the main issues in the RTM algorithm. To attenuate the undesired low-wavenumber noise, the causal imaging condition based on wavefield decomposition has been proposed. First, wavefield decompositions are performed to separate the wavefields as up-going and down-going wave components, respectively. Then, to preserve causality, it constructs images by correlating wave components that propagate in different directions. We build a causal imaging condition in this paper. Not only does it consider the up/down wavefield decomposition, but it also applies the decomposition on the horizontal direction to enhance the image quality especially for steeply dipping structures. The wavefield decomposition is conventionally achieved by the frequency-wavenumber (F-K) transform that is very computationally intensive compared with the wave propagation process of the RTM algorithm. To improve the efficiency of the algorithm, we propose a fast implementation to perform wavefield separation using the discrete Hilbert transform via the Graphics Processing Unit. Numerical tests on both the synthetic models and a real data example demonstrate the effectiveness of the proposed method and the efficiency of the optimized implementation scheme. This new imaging condition shows its ability to produce high image quality when applied to both the RTM stack image and also the angle domain common image gathers. The comparison of the total elapsed time for different methods verifies the efficiency of the optimized algorithm.

Sign in / Sign up

Export Citation Format

Share Document