Comparison of methods for extracting ADCIGs from RTM

Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. S89-S103 ◽  
Author(s):  
Hu Jin ◽  
George A. McMechan ◽  
Huimin Guan

Methods for extracting angle-domain common-image gathers (ADCIGs) during 2D reverse-time migration fall into three main categories; direction-vector-based methods, local-plane-wave decomposition methods, and local-shift imaging condition methods. The direction-vector-based methods, which use either amplitude gradients or phase gradients, cannot handle overlapping events because of an assumption of one propagation direction per imaging point per imaging time; however, the ADCIGs from the direction-vector-based methods have the highest angle resolution. A new direction-vector-based method using instantaneous phase gradients in space and time gives the same propagation directions and ADCIGs as those obtained by the Poynting vector or polarization vector based methods, where amplitudes are large. Angles calculated by the phase gradients have larger uncertainties at smaller amplitudes, but they do not significantly degrade the ADCIGs because they contribute only small amplitudes. The local-plane-wave decomposition and local-shift imaging condition methods, implemented either by a Fourier transform or by a slant stack transform, can handle overlapping events, and produce very similar angle gathers. ADCIGs from both methods depend on the local window size in which the transforms are done. In small local windows, both methods produce ADCIGs with low noise, but also with low angle resolution; in large windows, they have high angle resolution, but contain smeared artifacts.

Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. S69-S77 ◽  
Author(s):  
Xiongwen Wang ◽  
Jianliang Qian ◽  
Huazhong Wang

Because angle-domain common-image gathers (ADCIGs) from reverse time migration (RTM) are capable of obtaining the correct illumination of a subsurface geologic structure, they provide more reliable information for velocity model building, amplitude-variation versus angle analysis, and attribute interpretation. The approaches for generating ADCIGs mainly consist of two types: (1) indirect approaches that convert extended image gathers into ADCIGs and (2) direct approaches that first obtain propagating angles of wavefronts and then map the imaging result to the angle domain. In practice, however, generation of ADCIGs usually incurs high computational cost, poor resolution, and other drawbacks. To generate efficient ADCIGs using RTM methods, we have introduced a novel approach to obtain polarization vectors — directions of particle motion — from the Cauchy wavefield (CWF) and an efficient localized plane-wave decomposition algorithm to implement the angle-domain imaging condition. The CWF is a wavefield constructed from the Cauchy condition of the wave equation at any given time, and it only contains negative frequencies of the original wavefield so that the polarization vector is obtained from the local CWF in the wavenumber domain. With polarization vectors at our disposal, we have further developed an efficient localized plane-wave decomposition algorithm to implement the angle-domain imaging condition. Numerical examples have indicated that the new approach is able to handle complex wave phenomenon and has advantages in illuminating subsurface structure.


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.


Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. S1-S9 ◽  
Author(s):  
Jiangtao Hu ◽  
Huazhong Wang ◽  
Xiongwen Wang

Angle-domain common imaging gathers (ADCIGs) are important input data for migration velocity analysis and amplitude variation with angle analysis. Compared with Kirchhoff migration and one-way wave equation migration, reverse time migration (RTM) is the most accurate imaging method in complex areas, such as the subsalt area. We have developed a method to generate ADCIGs from RTM using analytic wavefield propagation and decomposition. To estimate the wave-propagation direction and angle by spatial Fourier transform during the time domain wave extrapolation, we have developed an analytic wavefield extrapolation method. Then, we decomposed the extrapolated source and receiver wavefields into their local angle components (i.e., local plane-wave components) and applied the angle-domain imaging condition to form ADCIGs. Because the angle-domain imaging condition is a convolution imaging condition about the source and receiver propagation angles, it is costly. To increase the efficiency of the angle-domain imaging condition, we have developed a local plane-wave decomposition method using matching pursuit. Numerical examples of synthetic and real data found that this method could generate high-quality ADCIGs. And these examples also found that the computational cost of this approach was related to the complexity of the source and receiver wavefields.


Geophysics ◽  
1989 ◽  
Vol 54 (10) ◽  
pp. 1339-1343 ◽  
Author(s):  
S. C. Singh ◽  
G. F. West ◽  
C. H. Chapman

The delay‐time (τ‐p) parameterization, which is also known as the plane‐wave decomposition (PWD) of seismic data, has several advantages over the more traditional time‐distance (t‐x) representation (Schultz and Claerbout, 1978). Plane‐wave seismograms in the (τ, p) domain can be used for obtaining subsurface elastic properties (P‐wave and S‐wave velocities and density as functions of depth) from inversion of the observed oblique‐incidence seismic data (e.g., Yagle and Levy, 1985; Carazzone, 1986; Carrion, 1986; Singh et al., 1989). Treitel et al. (1982) performed time migration of plane‐wave seismograms. Diebold and Stoffa (1981) used plane‐wave seismograms to derive a velocity‐depth function. Decomposing seismic data also allows more rapid modeling, since it is faster to compute synthetic seismograms in the (τ, p) than in the (t, x) domain. Unfortunately, the transformation of seismic data from the (t, x) to the (τ, p) domain may produce artifacts, such as those caused by discrete sampling, of the data in space.


Geophysics ◽  
1982 ◽  
Vol 47 (10) ◽  
pp. 1375-1401 ◽  
Author(s):  
Sven Treitel ◽  
P. R. Gutowski ◽  
D. E. Wagner

A point‐source seismic recording can be decomposed into a set of plane‐wave seismograms for arbitrary angles of incidence. Such plane‐wave seismograms possess an inherently simple structure that make them amenable to existing inversion methods such as predictive deconvolution. Implementation of plane‐wave decomposition (PWD) takes place in the frequency‐wavenumber domain under the assumption of radial symmetry. This version of PWD is equivalent to slant stacking if allowance is made for the customary use of linear recording arrays on the surface of a three‐dimensional medium. An imaging principle embodying both kinematic as well as dynamic characteristics allows us to perform time migration of the plane‐wave seismograms. The imaging procedure is implementable as a two‐dimensional filter whose independent variables are traveltime and angle of incidence.


Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. S111-S127 ◽  
Author(s):  
Qizhen Du ◽  
ChengFeng Guo ◽  
Qiang Zhao ◽  
Xufei Gong ◽  
Chengxiang Wang ◽  
...  

The scalar images (PP, PS, SP, and SS) of elastic reverse time migration (ERTM) can be generated by applying an imaging condition as crosscorrelation of pure wave modes. In conventional ERTM, Helmholtz decomposition is commonly applied in wavefield separation, which leads to a polarity reversal problem in converted-wave images because of the opposite polarity distributions of the S-wavefields. Polarity reversal of the converted-wave image will cause destructive interference when stacking over multiple shots. Besides, in the 3D case, the curl calculation generates a vector S-wave, which makes it impossible to produce scalar PS, SP, and SS images with the crosscorrelation imaging condition. We evaluate a vector-based ERTM (VB-ERTM) method to address these problems. In VB-ERTM, an amplitude-preserved wavefield separation method based on decoupled elastic wave equation is exploited to obtain the pure wave modes. The output separated wavefields are both vectorial. To obtain the scalar images, the scalar imaging condition in which the scalar product of two vector wavefields with source-normalized illumination is exploited to produce scalar images instead of correlating Cartesian components or magnitude of the vector P- and S-wave modes. Compared with alternative methods for correcting the polarity reversal of PS and SP images, our ERTM solution is more stable and simple. Besides these four scalar images, the VB-ERTM method generates another PP-mode image by using the auxiliary stress wavefields. Several 2D and 3D numerical examples are evaluated to demonstrate the potential of our ERTM method.


Sign in / Sign up

Export Citation Format

Share Document