Local vertical seismic profiling (VSP) elastic reverse-time migration and migration resolution: Salt-flank imaging with transmitted P-to-S waves

Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. S35-S49 ◽  
Author(s):  
Xiang Xiao ◽  
W. Scott Leaney
Keyword(s):  
Velocity Model ◽  
Staggered Grid ◽  
Image Point ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
S Waves ◽  
Seismic Profiling ◽  
Vertical Seismic Profiling ◽  
Time Migration

To avoid the defocusing effects of propagating waves through salt and overburden with an inaccurate overburden velocity model, we introduce a vertical seismic profiling (VSP) local elastic reverse-time-migration (RTM) method for salt-flank imaging by transmitted P-to-S waves. This method back-projects the transmitted PS waves using a local velocity model around the well until they are in phase with the back-projected PP waves at the salt boundaries. The merits of this method are that it does not require the complex overburden and salt-body velocities and it automatically accounts for source-side statics. In addition, the method accounts for kinematic and dynamic effects, including anisotropy, absorption, and all other unknown rock effects outside of this lo-cal subsalt velocity model. Numerical tests on an elastic salt model and offset 2D VSP data in the Gulf of Mexico, using a finite-difference time-domain staggered-grid RTM scheme, partly demonstrate the effectiveness of this method over interferometry PS-PP transmission migration and local acoustic RTM. Our method separates elastic wavefields to vector P- and S-wave velocity components at the trial image point and achieves better resolution than local acoustic RTM and interferometric transmission migration. The analytical formulas of migration resolution for local acoustic and elastic RTM show that the migration illumination is limited by data frequency and receiver aperture, and the spatial resolution is lower than standard poststack and prestack migration. This new method can image salt flanks as well as subsalt reflectors.

SP- and SS-imaging for 3D elastic reverse time migration

Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. A1-A6 ◽  
Author(s):  
Xufei Gong ◽  
Qizhen Du ◽  
Qiang Zhao
Keyword(s):  
Three Dimensional ◽  
Scalar Wave ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Sh Waves ◽  
S Waves ◽  
Time Migration ◽  
Underlying Principle ◽  
Imaging Conditions

Three-dimensional elastic reverse time migration has been confronted with the problem of generating scalar images with vector S-waves. The underlying principle for solving this problem is to convert the vector S-waves into scalars. Previous methods were mainly focused on PS-imaging, but they usually cannot work properly on SP- and SS-cases. The complexity of SP- and SS-imaging arises from the fact that the incident S-wave has unpredictable relationship with the raypath plane. We have suggested that S-wave should be treated separately as SV- and SH-waves, which keep predictable relationships with the raypath plane. First, the elastic wavefield is separated into P- and S-waves using the Helmholtz decomposition. Then, we evaluate the normal direction of the raypath plane at each imaging grid. Next, we separate the vector S-wave obtained with curl operator into SH- and SV-waves, both of which are scalars. Finally, correlation imaging conditions are implemented to those scalar wave modes to produce scalar SV-P, SV-SV, and SH-SH images.

Elastic wavefield separation based on the Helmholtz decomposition

Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. S173-S183 ◽  
Author(s):  
Hejun Zhu
Keyword(s):  
Vector Fields ◽  
Helmholtz Decomposition ◽  
Elastic Media ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
S Waves ◽  
Poisson Solver ◽  
Time Migration ◽  
Imaging Conditions

Divergence and curl operators used for the decomposition of P- and S-wave modes in elastic reverse time migration (RTM) change the amplitudes, units, and phases of extrapolated wavefields. I separate the P- and S-waves in elastic media based on the Helmholtz decomposition. The decomposed wavefields based on this approach have the same amplitudes, units, and phases as the extrapolated wavefields. To avoid expensive multidimensional integrals in the Helmholtz decomposition, I introduce a fast Poisson solver to efficiently solve the vector Poisson’s equation. This fast algorithm allows us to reduce computational complexity from [Formula: see text] to [Formula: see text], where [Formula: see text] is the total number of grid points. Because the decomposed P- and S-waves are vector fields, I use vector imaging conditions to construct PP-, PS-, SS-, and SP-images. Several 2D numerical examples demonstrate that this approach allows us to accurately and efficiently decompose P- and S-waves in elastic media. In addition, elastic RTM images based on the vector imaging conditions have better quality and avoid polarity reversal in comparison with images based on the divergence and curl separation or direct component-by-component crosscorrelation.

Multidirectional-vector-based elastic reverse time migration and angle-domain common-image gathers with approximate wavefield decomposition of P- and S-waves

Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. S57-S79 ◽  
Author(s):  
Chen Tang ◽  
George A. McMechan
Keyword(s):  
Wave Mode ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
S Waves ◽  
Time Migration ◽  
Mode Separation ◽  
Reflection Image ◽  
Wavefield Decomposition ◽  
Imaging Conditions

Elastic reverse time migration (E-RTM) has limitations when the migration velocities contain strong contrasts. First, the traditional scheme of P/S-wave mode separation is based on Helmholtz’s equations, which ignore the conversion between P- and S-waves at the current separation time. Thus, it contains an implicit assumption of the constant shear modulus and requires smoothing the heterogeneous model to approximately satisfy a locally constant condition. Second, the vector-based imaging condition needs to use the reflection-image normal, and it also cannot give the correct polarity of the PP image in all possible conditions. Third, the angle-domain common-image gathers (ADCIGs) calculated using the Poynting vectors (PVs) do not consider the wave interferences that happen at each reflector. Therefore, smooth models are often used for E-RTM. We relax this condition by proposing an improved data flow that involves three new contributions. The first contribution is an improved system of P/S-wave mode separation that considers the converted wave generated at the current time, and thus it does not require the constant-shear-modulus assumption. The second contribution is the new elastic imaging conditions based on multidirectional vectors; they can give the correct image polarity in all possible conditions without knowledge of the reflection-image normal. The third contribution is two methods to calculate multidirectional propagation vectors (PRVs) for RTM images and ADCIGs: One is the elastic multidirectional PV, and the other uses the sign of wavenumber-over-frequency ([Formula: see text]) ratio obtained from an amplitude-preserved approximate-propagation-angle-based wavefield decomposition to convert the particle velocities into multidirectional PRVs. The robustness of the improved data flow is determined by several 2D numerical examples. Extension of the schemes into 3D and amplitude-preserved imaging conditions is also possible.

Elastic reverse time migration of marine walkaway vertical seismic profiling data

Geophysics ◽  
1998 ◽  
Vol 63 (5) ◽  
pp. 1685-1695 ◽  
Author(s):  
Ketil Hokstad ◽  
Rune Mittet ◽  
Martin Landrø
Keyword(s):  
Vertical Component ◽  
Horizontal Component ◽  
Reciprocal Relationship ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Seismic Profiling ◽  
Vertical Seismic Profiling ◽  
Depth Images ◽  
Time Migration ◽  
Walkaway Vsp

Walkaway vertical seismic profiling (VSP) acquisition with three‐component geophones allows for direct measurement of compressional as well as shear energy. This makes full elastic reverse time migration an attractive alternative for imaging data. We present results from elastic reverse time migration of a marine walkaway VSP acquired offshore Norway. The reverse time migration scheme is based on a high‐order finite‐difference solution to the two‐way elastic wave equation. Depth images of the subsurface are constructed by correlation of forward‐ and back‐propagated elastic wavefields. In the walkaway VSP configuration, the number of shots is much larger than the number of geophone levels. Using processing methods operating in the shot/receiver domain, it is advantageous to use the reciprocal relationship between the walkaway VSP and the reverse VSP configurations. We do this by imaging each component of each geophone level as a reverse VSP common shot gather. The final images are constructed by stacking partial images from each level. The depth images obtained from the vertical components reveal the major characteristics of the geological structure below geophone depth. A graben in the base Cretaceous unconformity and a faulted coal layer can be identified. The horizontal components are more difficult to image. Compared to the vertical components, the horizontal component images are more corrupted by migration artifacts. This is because the horizontal component images are more sensitive to aperture effects and to the shear‐wave velocity macromodel. When converted to two‐way time, the migration results tie well with the surface seismic section. Comparison of fully elastic and acoustic reverse time migration shows that the vertical component is dominantly PP-reflected events, whereas the horizontal components get important contributions from PS-converted energy. The horizontal components also provide higher resolution because of the shorter wavelength of the shear waves.

A wavefield-separation-based elastic least-squares reverse time migration

Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. S279-S297 ◽  
Author(s):  
Bingluo Gu ◽  
Zhenchun Li ◽  
Jianguang Han
Keyword(s):  
Least Squares ◽  
Synthetic Data ◽  
P Wave ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
S Waves ◽  
Numerical Tests ◽  
Time Migration

Elastic least-squares reverse time migration (ELSRTM) has the potential to provide improved subsurface reflectivity estimation. Compared with elastic RTM (ERTM), ELSRTM can produce images with higher spatial resolution, more balanced amplitudes, and fewer artifacts. However, the crosstalk between P- and S-waves can significantly degrade the imaging quality of ELSRTM. We have developed an ELSRTM method to suppress the crosstalk artifacts. This method includes three crucial points. The first is that the forward and backward wavefields are extrapolated based on the separated elastic velocity-stress equation of P- and S-waves. The second is that the separated vector P- and S-wave residuals are migrated to form reflectivity images of Lamé constants [Formula: see text] and [Formula: see text] independently. The third is that the reflectivity images of [Formula: see text] and [Formula: see text] are obtained by the vector P-wave wavefields achieved in the backward extrapolation of the separated vector P-wave residuals and the vector S-wave wavefields achieved in the backward extrapolation of the separated vector S-wave residuals, respectively. Numerical tests with synthetic data demonstrate that our ELSRTM method can produce images free of crosstalk artifacts. Compared with ELSRTM based on the coupled wavefields, our ELSRTM method has better convergence and higher accuracy.

Vector-wave-based elastic reverse time migration of ocean-bottom 4C seismic data

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. S333-S343 ◽  
Author(s):  
Pengfei Yu ◽  
Jianhua Geng ◽  
Jiqiang Ma
Keyword(s):  
Seismic Data ◽  
Ocean Bottom ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Vector Wave ◽  
S Waves ◽  
Time Migration ◽  
Multicomponent Seismic ◽  
Wavefield Extrapolation

The acoustic-elastic coupled equation (AECE) has several advantages when compared with conventional scalar-wave-based elastic reverse time migration (ERTM) methods used to image ocean-bottom multicomponent seismic data. In particular, vector-wave-based ERTM requires vectorial P- and S-waves on the source and receiver sides, but these cannot be directly obtained from wavefield extrapolation using AECE. Therefore, we have developed a P- and S-wave vector decomposition (VD) approach within AECE; this approach enables the deduction of a novel VD-based AECE, from which vectorial P- and S-waves can be obtained directly via wavefield extrapolation. We are also able to derive a new formulation suitable for vector-wave-based ERTM of ocean-bottom multicomponent seismic data that can generate a phase-preserved PS-image. Three synthetic examples illustrate the validity and effectiveness of our new method.

Application of RTM 3D angle gathers to wide-azimuth data subsalt imaging

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB175-WB182 ◽  
Author(s):  
Yan Huang ◽  
Bing Bai ◽  
Haiyong Quan ◽  
Tony Huang ◽  
Sheng Xu ◽  
...  
Keyword(s):  
Model Updating ◽  
Velocity Model ◽  
Azimuth Angle ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Additional Information ◽  
Time Migration ◽  
Reflection Angle ◽  
Subsalt Imaging

The availability of wide-azimuth data and the use of reverse time migration (RTM) have dramatically increased the capabilities of imaging complex subsalt geology. With these improvements, the current obstacle for creating accurate subsalt images now lies in the velocity model. One of the challenges is to generate common image gathers that take full advantage of the additional information provided by wide-azimuth data and the additional accuracy provided by RTM for velocity model updating. A solution is to generate 3D angle domain common image gathers from RTM, which are indexed by subsurface reflection angle and subsurface azimuth angle. We apply these 3D angle gathers to subsalt tomography with the result that there were improvements in velocity updating with a wide-azimuth data set in the Gulf of Mexico.

Sign in / Sign up

Export Citation Format

Share Document