Common-image gathers in the incident phase-angle domain from reverse time migration in 2D elastic VTI media

Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. S197-S206 ◽  
Author(s):  
Qunshan Zhang ◽  
George A. McMechan
Keyword(s):  
Phase Angle ◽  
Incident Angle ◽  
Image Point ◽  
Common Source ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Wave Phase ◽  
Time Migration ◽  
Imaging Time ◽  
Vti Media

Reverse time migration (RTM) was implemented with a modified crosscorrelation imaging condition for data from 2D elastic vertically transversely isotropy (VTI) media. The computation cost was reduced because scalar qP- and qS-wavefield separations are performed in VTI media, for the source and receiver wavefields only at the RTM imaging time, to calculate the migrated qP and qS images. Angle-domain common-image gathers (CIGs) were extracted from qPqP and qPqS common-source RTM images. The local incident angle was produced as the difference between the qP-wave phase angle, obtained directly from the source wavefield polarization, and the normal to the reflector, calculated as the instantaneous wavenumber direction via a directional Hilbert transform of the stacked image. Angle-domain CIGs were extracted by reordering the prestack-migrated images by local incident phase angle, source by source. Vector decomposition of the source qP-wavefield was required to calculate the qP-wave phase polarization direction for each image point at its imaging time. RTM and angle-domain CIG extraction were successfully implemented and illustrated with a synthetic 2D elastic VTI example.

Plane-wave least-squares reverse time migration in complex VTI media

2016 ◽  
Author(s):  
Jinqiang Huang ◽  
Daojun Si ◽  
Zhenchun Li ◽  
Jianping Huang
Keyword(s):  
Plane Wave ◽  
Least Squares ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Vti Media

Direct vector-field method to obtain angle-domain common-image gathers from isotropic acoustic and elastic reverse time migration

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB135-WB149 ◽  
Author(s):  
Qunshan Zhang ◽  
George A. McMechan
Keyword(s):  
Direct Method ◽  
Incident Angle ◽  
Field Method ◽  
Propagation Direction ◽  
P Wave ◽  
Common Source ◽  
Spatial Gradient ◽  
Reverse Time ◽  
S Wave ◽  
Time Migration

We have developed an alternative (new) method to produce common-image gathers in the incident-angle domain by calculating wavenumbers directly from the P-wave polarization rather than using the dominant wavenumber as the normal to the source wavefront. In isotropic acoustic media, the wave propagation direction can be directly calculated as the spatial gradient direction of the acoustic wavefield, which is parallel to the wavenumber direction (the normal to the wavefront). Instantaneous wavenumber, obtained via a novel Hilbert transform approach, is used to calculate the local normal to the reflectors in the migrated image. The local incident angle is produced as the difference between the propagation direction and the normal to the reflector. By reordering the migrated images (over all common-source gathers) with incident angle, common-image gathers are produced in the incident-angle domain. Instantaneous wavenumber takes the place of the normal to the reflector in the migrated image. P- and S-wave separations allow both PP and PS common-image gathers to be calculated in the angle domain. Unlike the space-shift image condition for calculating the common-image gather in angle domain, we use the crosscorrelation image condition, which is substantially more efficient. This is a direct method, and is less dependent on the data quality than the space-shift method. The concepts were successfully implemented and tested with 2D synthetic acoustic and elastic examples, including a complicated (Marmousi2) model that illustrates effects of multipathing in angle-domain common-image gathers.

3D angle decomposition for elastic reverse time migration

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. S377-S389
Author(s):  
Yuting Duan ◽  
Paul Sava
Keyword(s):  
Field Data ◽  
Time Lag ◽  
Image Point ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Converted Wave ◽  
The Third ◽  
Time Migration ◽  
Straight Lines ◽  
2D And 3D

We have developed three approaches for 3D angle decomposition using elastic reverse time migration. The first approach uses time- and space-lag common-image point gathers computed from elastic wavefields. This method facilitates computing angle gathers at sparse and possibly irregularly distributed points in the image. The second approach transforms extended time-lag images to the angle domain using slant stacks along 4D surfaces, instead of using slant stacks along 2D straight lines. The third approach transforms space-lag common-image gathers to the angle domain. The three proposed methods solve a system of equations that handles dipping reflectors, and they yield angle gathers that are more accurate compared with those obtained via alternative existing methods. We have developed our methods using 2D and 3D synthetic and field data examples and found that they provide accurate opening and azimuth angles and they can handle steeply dipping reflectors and converted wave modes.

Elastic reverse‐time migration

Geophysics ◽  
1987 ◽  
Vol 52 (10) ◽  
pp. 1365-1375 ◽  
Author(s):  
Wen‐Fong Chang ◽  
George A. McMechan
Keyword(s):  
Finite Difference ◽  
Vertical Component ◽  
Common Source ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Imaging Condition ◽  
Time Migration ◽  
Two Component ◽  
Synthetic Test ◽  
Data Extrapolation

Elastic, prestack, reverse‐time, finite‐difference migration of two‐component seismic surface data requires data extrapolation and application of an imaging condition. Data extrapolation involves synchronous driving of the vertical‐component and horizontal‐component finite‐difference meshes with the time reverse of the recorded vertical and horizontal traces, respectively. Extrapolation uses the coupled elastic wave equation for variable velocity solved with a second‐order, explicit finite‐difference scheme. The imaging condition at any point in the grid is the one‐way traveltime from the source to that point. Elastic migrations of both synthetic test data and real two‐component common‐source gathers produce simpler images than acoustic migrations because of the coalescing of double reflections (compressional waves and shear waves) into single loci.

3-D prestack migration in anisotropic media

Geophysics ◽  
1993 ◽  
Vol 58 (1) ◽  
pp. 79-90 ◽  
Author(s):  
Zhengxin Dong ◽  
George A. McMechan
Keyword(s):  
A Priori ◽  
Three Dimensional ◽  
Synthetic Data ◽  
Anisotropic Media ◽  
P Wave ◽  
Common Source ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Prestack Migration ◽  
Time Migration

A three‐dimensional (3-D) prestack reverse‐time migration algorithm for common‐source P‐wave data from anisotropic media is developed and illustrated by application to synthetic data. Both extrapolation of the data and computation of the excitation‐time imaging condition are implemented using a second‐order finite‐ difference solution of the 3-D anisotropic scalar‐wave equation. Poorly focused, distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation; well focused, clear images are obtained using anisotropic extrapolation. A priori estimation of the 3-D anisotropic velocity distribution is required. Zones of anomalous, directionally dependent reflectivity associated with anisotropic fracture zones are detectable in both the 3-D common‐ source data and the corresponding migrated images.

Pseudo-acoustic reverse-time migration in VTI media and its application to fracture-vuggy imaging

2014 ◽  
Author(s):  
Liu Wenqing ◽  
Wang Yuchao ◽  
Yong Xueshan ◽  
Wang Yanxiang ◽  
Shao Xichun
Keyword(s):  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Vti Media

Angle- and azimuth-domain common-image gathers by reverse time migration for 3D elastic isotropic media

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. S151-S167
Author(s):  
Zabihollah Khaksar ◽  
George A. McMechan
Keyword(s):  
Incident Angle ◽  
Migration Velocity ◽  
Spatial Gradient ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Disk Storage ◽  
Time Migration ◽  
Isotropic Media ◽  
The Hilbert Transform ◽  
Elastic Particle

A 2D algorithm for angle-domain common-image gather (CIG) calculation is extended and modified to produce 3D elastic angle and azimuth CIGs. The elastic seismic data are propagated with the elastic particle displacement wave equation, and then the PP-reflected and PS-converted waves are separated by divergence and curl calculations during application of the excitation-time imaging condition. The incident angles and azimuths are calculated using source propagation directions and the reflector normals. The source propagation direction vector is computed as the spatial gradient of the incident 3C P-wavefield. The vector normal to the reflector is calculated using the Hilbert transform. Ordering the migrated images with respect to incident angles for a fixed azimuth bin, or with respect to azimuths for a fixed incident angle bin, creates angle- or azimuth-domain CIGs, respectively. Sorting the azimuth gathers by the incident angle bins causes a shift to a greater depth for too-high migration velocity and to a smaller depth for too-low migration velocity. For the sorted incident angle gathers, the velocity-dependent depth moveout is within the angle gathers and across the azimuth gathers. This method is compared with three other 3D CIG algorithms with respect to the number of calculations and their disk storage and RAM requirements; it is three to six orders of magnitude faster and requires two to three orders of magnitude less disk space. The method is successfully tested with data for a modified part of the SEG/EAGE overthrust model.

Sign in / Sign up

Export Citation Format

Share Document