Identifying, removing, and imaging P‐S conversions at salt‐sediment interfaces

Geophysics ◽  
2003 ◽  
Vol 68 (3) ◽  
pp. 1052-1059 ◽  
Author(s):  
Richard S. Lu ◽  
Dennis E. Willen ◽  
Ian A. Watson
Keyword(s):  
P Wave ◽  
Time Interval ◽  
S Wave ◽  
Converted Wave ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Time Migration ◽  
Wave Imaging ◽  
Wave Image ◽  
Converted Waves

The large velocity contrast between salt and the surrounding sediments generates strong conversions between P‐ and S‐wave energy. The resulting converted events can be noise on P‐wave migrated images and should be identified and removed to facilitate interpretation. On the other hand, they can also be used to image a salt body and its adjacent sediments when the P‐wave image is inadequate. The converted waves with smaller reflection and transmission angles and much larger critical angles generate substantially different illumination than does the P‐wave. In areas where time migration is valid, the ratio between salt thickness in time and the time interval between the P‐wave and the converted‐wave salt base on a time‐migrated image is about 2.6 or 1.3, depending upon whether the seismic wave propagates along one or both of the downgoing and upcoming raypaths in salt as the S‐wave, respectively. These ratios can be used together with forward seismic modeling and 2D prestack depth migration to identify the converted‐wave base‐of‐salt (BOS) events in time and depth and to correctly interpret the subsalt sediments. It is possible to mute converted‐wave events from prestack traces according to their computed arrival times. Prestack depth migration of the muted data extends the updip continuation of subsalt sedimentary beds, and improves the salt–sediment terminations in the P‐wave image. Prestack and poststack depth‐migrated examples illustrate that the P‐wave and the three modes of converted waves preferentially image different parts of the base of salt. In some areas, the P‐wave BOS can be very weak, obscured by noise, or completely absent. Converted‐wave imaging complements P‐wave imaging in delineating the BOS for velocity model building.

Examining the processing differences between P and P-S waves in a Rocky Mountain Foothills model

2014 ◽  
Vol 54 (2) ◽  
pp. 504
Author(s):  
Sanjeev Rajput ◽  
Michael Ring
Keyword(s):  
Rocky Mountain ◽  
P Wave ◽  
S Wave ◽  
Converted Wave ◽  
Wave Source ◽  
Depth Migration ◽  
S Waves ◽  
Full Waveform ◽  
Fluid Saturation ◽  
Converted Waves

For the past two decades, most of the shear-wave (S-wave) or converted wave (P-S) acquisitions were performed with P-wave source by making the use of downgoing P-waves converting to upgoing S-waves at the mode conversion boundaries. The processing of converted waves requires studying asymmetric reflection at the conversion point, difference in geometries and conditions of source and receiver, and the partitioning of energy into orthogonally polarised components. Interpretation of P-S sections incorporates the identification of P-S waves, full waveform modeling, correlation with P-wave sections and depth migration. The main applications of P-S wave imaging are to obtain a measure of subsurface S-wave properties relating to rock type and fluid saturation (in addition to the P-wave values), imaging through gas clouds and shale diapers, and imaging interfaces with low P-wave contrast but significant S-wave changes. This study examines the major differences in processing of P and P-S wave surveys and the feasibility of identifying converted mode reflections by P-wave sources in anisotropic media. Two-dimensional synthetic seismograms for a realistic rocky mountain foothills model were studied. A Kirchhoff-based technique that includes anisotropic velocities is used for depth migration of converted waves. The results from depth imaging show that P-S section help in distinguishing amplitude associated with hydrocarbons from those caused by localised stratigraphic changes. In addition, the full waveform elastic modeling is useful in finding an appropriate balance between capturing high-quality P-wave data and P-S data challenges in a survey.

Examining the processing differences between P and P-S waves in a Rocky Mountain Foothills model

2014 ◽  
Vol 54 (2) ◽  
pp. 536
Author(s):  
Sanjeev Rajput ◽  
Michael Ring
Keyword(s):  
Rocky Mountain ◽  
P Wave ◽  
S Wave ◽  
Converted Wave ◽  
Wave Source ◽  
Depth Migration ◽  
S Waves ◽  
Full Waveform ◽  
Fluid Saturation ◽  
Converted Waves

For the past two decades, most of the shear-wave (S-wave) or converted wave (P-S) acquisitions were performed with P-wave source by making the use of downgoing P-waves converting to upgoing S-waves at the mode conversion boundaries. The processing of converted waves requires studying asymmetric reflection at the conversion point, difference in geometries and conditions of source and receiver, and the partitioning of energy into orthogonally polarised components. Interpretation of P-S sections incorporates the identification of P-S waves, full waveform modeling, correlation with P-wave sections and depth migration. The main applications of P-S wave imaging are to obtain a measure of subsurface S-wave properties relating to rock type and fluid saturation (in addition to the P-wave values), imaging through gas clouds and shale diapers, and imaging interfaces with low P-wave contrast but significant S-wave changes. This study examines the major differences in processing of P and P-S wave surveys and the feasibility of identifying converted mode reflections by P-wave sources in anisotropic media. Two-dimensional synthetic seismograms for a realistic rocky mountain foothills model were studied. A Kirchhoff-based technique that includes anisotropic velocities is used for depth migration of converted waves. The results from depth imaging show that P-S section help in distinguishing amplitude associated with hydrocarbons from those caused by localised stratigraphic changes. In addition, the full waveform elastic modeling is useful in finding an appropriate balance between capturing high-quality P-wave data and P-S data challenges in a survey.

Multicomponent prestack depth migration using the elastic screen method

Geophysics ◽  
2005 ◽  
Vol 70 (1) ◽  
pp. S30-S37 ◽  
Author(s):  
Xiao-Bi Xie ◽  
Ru-Shan Wu
Keyword(s):  
Synthetic Data ◽  
Data Sets ◽  
Converted Wave ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
S Waves ◽  
Angle Correction ◽  
Scattering Effects ◽  
Wave Imaging ◽  
Migration Method

A 3D multicomponent prestack depth-migration method is presented. An elastic-screen propagator based on one-way wave propagation with a wide-angle correction is used to extrapolate both source and receiver wavefields. The elastic-screen propagator neglects backscattered waves but can handle forward multiple-scattering effects, such as focusing/defocusing, diffraction, interference, and conversions between P- and S-waves. Vector-imaging conditions are used to generate a P-P image and a P-S converted-wave image. The application of the multicomponent elastic propagator and vector-imaging condition preserves more information carried by the elastic waves. It also solves the polarization problem of converted-wave imaging. Partial images from different sources with correct polarizations can be stacked to generate a final image. Numerical examples using 2D synthetic data sets are presented to show the feasibility of this method.

Elastic least-squares reverse time migration with hybrid l1/l2 misfit function

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. S271-S291 ◽  
Author(s):  
Bingluo Gu ◽  
Zhenchun Li ◽  
Peng Yang ◽  
Wencai Xu ◽  
Jianguang Han
Keyword(s):  
Least Squares ◽  
Seismic Data ◽  
Synthetic Data ◽  
Adjoint Equations ◽  
P Wave ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Time Migration ◽  
Wave Image

We have developed the theory and synthetic tests of elastic least-squares reverse time migration (ELSRTM). In this method, a least-squares reverse time migration algorithm is used to image multicomponent seismic data based on the first-order elastic velocity-stress wave equation, in which the linearized elastic modeling equations are used for forward modeling and its adjoint equations are derived based on the adjoint-state method for back propagating the data residuals. Also, we have developed another ELSRTM scheme based on the wavefield separation technique, in which the P-wave image is obtained using P-wave forward and adjoint wavefields and the S-wave image is obtained using P-wave forward and S-wave adjoint wavefields. In this way, the crosstalk artifacts can be minimized to a significant extent. In general, seismic data inevitably contain noise. We apply the hybrid [Formula: see text] misfit function to the ELSRTM algorithm to improve the robustness of our ELSRTM to noise. Numerical tests on synthetic data reveal that our ELSRTM, when compared with elastic reverse time migration, can produce images with higher spatial resolution, more-balanced amplitudes, and fewer artifacts. Moreover, the hybrid [Formula: see text] misfit function makes the ELSRTM more robust than the [Formula: see text] misfit function in the presence of noise.

scholarly journals Three-component amplitude versus offset analysis

1989 ◽  
Vol 20 (2) ◽  
pp. 257
Author(s):  
D.R. Miles ◽  
G. Gassaway ◽  
L. Bennett ◽  
R. Brown
Keyword(s):  
Seismic Data ◽  
P Wave ◽  
S Wave ◽  
Converted Wave ◽  
P Waves ◽  
S Waves ◽  
Avo Analysis ◽  
Avo Inversion ◽  
Amplitude Versus Offset ◽  
Converted Waves

Three-component (3-C) amplitude versus offset (AVO) inversion is the AVO analysis of the three major energies in the seismic data, P-waves, S-waves and converted waves. For each type of energy the reflection coefficients at the boundary are a function of the contrast across the boundary in velocity, density and Poisson's ratio, and of the angle of incidence of the incoming wave. 3-C AVO analysis exploits these relationships to analyse the AVO changes in the P, S, and converted waves. 3-C AVO analysis is generally done on P, S, and converted wave data collected from a single source on 3-C geophones. Since most seismic sources generate both P and S-waves, it follows that most 3-C seismic data may be used in 3-C AVO inversion. Processing of the P-wave, S-wave and converted wave gathers is nearly the same as for single-component P-wave gathers. In split-spread shooting, the P-wave and S-wave energy on the radial component is one polarity on the forward shot and the opposite polarity on the back shot. Therefore to use both sides of the shot, the back shot must be rotated 180 degrees before it can be stacked with the forward shot. The amplitude of the returning energy is a function of all three components, not just the vertical or radial, so all three components must be stacked for P-waves, then for S-waves, and finally for converted waves. After the gathers are processed, reflectors are picked and the amplitudes are corrected for free-surface effects, spherical divergence and the shot and geophone array geometries. Next the P and S-wave interval velocities are calculated from the P and S-wave moveouts. Then the amplitude response of the P and S-wave reflections are analysed to give Poisson's ratio. The two solutions are then compared and adjusted until they match each other and the data. Three-component AVO inversion not only yields information about the lithologies and pore-fluids at a specific location; it also provides the interpreter with good correlations between the P-waves and the S-waves, and between the P and converted waves, thus greatly expanding the value of 3-C seismic data.

Prestack reverse time imaging in tunnels based on the decoupled nonconversion elastic-wave equation

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. S371-S383
Author(s):  
Yuxiao Ren ◽  
Zhichao Yang ◽  
Bin Liu ◽  
Xinji Xu ◽  
Yangkang Chen
Keyword(s):  
Seismic Data ◽  
Accurate Method ◽  
Wave Mode ◽  
P Wave ◽  
Tunnel Construction ◽  
Reverse Time ◽  
S Wave ◽  
Tunnel Face ◽  
Time Migration ◽  
Converted Waves

The safety and efficiency of tunnel construction depend on the knowledge of complex geologic conditions. Hence, an accurate forward-prospecting technique is required to detect unexpected inhomogeneous geologic structures ahead of tunnel construction. As an accurate method to image geologic heterogeneities, elastic reverse time migration (ERTM) is introduced to the field of tunnel forward prospecting. However, in the tunnel environment, ERTM images may suffer from the interference of crosstalk artifacts, which are caused by converted waves on tunnel surfaces. Therefore, considering the actual influence of the tunnel body, the decoupled nonconversion elastic equation was incorporated into traditional ERTM. This method prevents the generation of converted waves but ensures independent P- and S-wave propagation. In addition, wave-mode separation for raw seismic data is required in our approach. Synthetic examples based on the real geologic environment of tunnels show that our method produces satisfactory results for P-wave and S-wave imaging and the S-wave can produce a better imaging effect in tunnels. Finally, we apply our method to the seismic data obtained from a real highway tunnel construction site to demonstrate its performance in real-world applications. The results indicate that the migrated images can help to accurately constrain the geologic formations ahead of the tunnel face.

Vector-based elastic reverse time migration based on scalar imaging condition

Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. S111-S127 ◽  
Author(s):  
Qizhen Du ◽  
ChengFeng Guo ◽  
Qiang Zhao ◽  
Xufei Gong ◽  
Chengxiang Wang ◽  
...  
Keyword(s):  
Alternative Methods ◽  
Polarity Reversal ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Converted Wave ◽  
Imaging Condition ◽  
Time Migration ◽  
Wavefield Separation ◽  
Wave Modes

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.

Converted-Wave Prestack Depth Migration of North Sea Salt Domes

1999 ◽  
Author(s):  
X. Zhu ◽  
J. Langhammer ◽  
D. King ◽  
E. Madtson ◽  
H. K. Helgesen ◽  
...  
Keyword(s):  
North Sea ◽  
Sea Salt ◽  
Converted Wave ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Salt Domes

Common-Reflection-Point-Based Prestack Depth Migration for Imaging Lithosphere in Python: Application to the Dense Warramunga Array in Northern Australia

2020 ◽  
Vol 91 (5) ◽  
pp. 2890-2899 ◽  
Author(s):  
Weijia Sun ◽  
Brian L. N. Kennett
Keyword(s):  
Open Source ◽  
Receiver Function ◽  
Function Analysis ◽  
Vertical Component ◽  
P Wave ◽  
Northern Australia ◽  
Reflection Point ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Conversion Point

Abstract We exploit estimates of P-wave reflectivity from autocorrelation of transmitted teleseismic P arrivals and their coda in a common reflection point (CRP) migration technique. The approach employs the same portion of the vertical-component seismogram, as in standard Ps receiver function analysis. This CRP prestack depth migration approach has the potential to image lithospheric structures on scales as fine as 4 km or less. The P-wave autocorrelation process and migration are implemented in open-source software—the autocorrelogram calculation (ACC) package, which builds on the widely used the seismological Obspy toolbox. The ACC package is written in the open-source and free Python programming language (3.0 or newer) and has been extensively tested in an Anaconda Python environment. The package is simple and friendly to use and runs on all major operating systems (e.g., Windows, macOS, and Linux). We utilize Python multiprocessing parallelism to speed up the ACC on a personal computer system, or servers, with multiple cores and threads. The application of the ACC package is illustrated with application to the closely spaced Warramunga array in northern Australia. The results show how fine-scale structures in the lithospheric can be effectively imaged at relatively high frequencies. The Moho ties well with conventional H−κ receiver analysis and deeper structure inferred from stacked autocorrelograms for continuous data. CRP prestack depth migration provides an important complement to common conversion point receiver function stacks, since it is less affected by surface multiples at lithospheric depths.

Sign in / Sign up

Export Citation Format

Share Document