Improved AVA imaging for vector-based elastic reverse time migration using angle-domain dot-product imaging condition

Reverse time migration (RTM) is achieved by a forward and reverse time propagation of source and receiver wavefields respectively, followed by an imaging condition. The quality of the image from any RTM implementation depends directly on the method’s ability to separate the true reflection data from the backscattered correlation noise between the source and receiver wavefields. In this paper we demonstrate the application of an inverse scattering imaging condition that significantly reduces this backscattered correlation noise and results in a much higher quality subsurface image than is achieved by typical RTM imaging conditions.

Vector-based excitation amplitude imaging condition for elastic RTM

The earth [Formula: see text] filtering causes poor illumination of the subsurface. Compensating for the attenuated amplitude and distorted phase is crucial during elastic reverse time migration (ERTM) to improve the imaging quality. Conventional [Formula: see text]-compensated ERTM ([Formula: see text]-ERTM) methods tend to boost the attenuated energy to inverse the [Formula: see text] effects. These methods usually suffer from severe numerical instability because of the unlimited amplification of the high-frequency noise. Low-pass filtering is generally used to stabilize the process, however, at the expense of precision. We have developed a stable compensation approach in this paper. Based on the decoupled fractional Laplacians viscoelastic wave equation, two compensation operators are obtained by extrapolating wavefield in the dispersion-only and viscoelastic media. Because no explicit amplification is included, these two operators are absolutely stable for implementation. To improve the division morbidity for calculating the compensation operators, we use the excitation amplitude criterion and embed the operators into a vector-based [Formula: see text]-compensated excitation amplitude imaging condition. With the derived imaging condition, we could compensate for the absorption accurately without needing to concern the stability issue. The [Formula: see text]-ERTM results for noise-free data are carried out over a simple layered model and a more realistic Marmousi model with an attenuating area to verify the feasibility of the proposed approach. The migration results for noisy data from the Marmousi model further prove that the proposed method performs better fidelity, adaptability, and antinoise performance compared with conventional compensation method with low-pass filtering.

Elastic reverse time migration (ERTM), which is capable of using multicomponent seismic data, provides not only an improvement of the P-P image compared to the one derived from acoustic RTM, but also more opportunities to understand the subsurface through converted wave images (P-S, S-P, and S-S images). However, the polarity reversals in P-S and S-P images and cross-talk noises generated in S-P and S-S images degrade the converted wave images of ERTM. To overcome these problems, we derive a new P-S converted wave imaging condition for 2D ERTM based on wavefield separation techniques. The proposed imaging condition, called converted-wave guided (CWG) imaging condition, incorporates an extra term that represents the sign and wavelength of S-waves converted from source wavefields into the zero-lag cross-correlation imaging condition for P-S imaging. The extra term compensates for the polarity reversal of separated S-waves from receiver wavefields because the converted S-waves from source wavefields also have the change in polarity. In addition, since this CWG imaging condition produces images where P- and S-waves from source wavefields and S-waves from receiver wavefields coincide, image resolution is enhanced without generating spurious events. Our approach is motivated by the specific feature of ERTM that generates converted waves at the reflection points (conventional imaging points) when proper elastic models are used. Through a numerical experiment with a simple elastic model, we demonstrate that the proposed CWG imaging condition successfully corrects the polarity reversal and provides higher image resolution. We also test our migration algorithm on a synthetic ocean bottom cable (OBC) dataset created using the Marmousi-II model. The P-S image obtained from CWG imaging condition shows continuous events and improved image resolution.

It is important to apply certain intrinsic properties (e.g., the wave propagation directions and particle motion information) in reconstructed wavefields in order to solve the image problems exhibited in elastic reverse time migration (RTM). This paper presents a procedure to decompose the source and receiver wavefields into local plane waves as well as to separate them into pure P and S modes. We generate the partial PP and PS images by crosscorrelation of these plane waves along different directions, and then formulate an imaging condition as a product of an angle-domain operator and the partial images. The new angle-domain imaging condition substantially reduces the artifacts in the PP image and produces the PS image with correct polarizations. However, the imaging procedure involves intensive computations. A numerical algorithm, which can greatly reduce the calculation, is proposed to improve the efficiency.

Migration is an important step in seismic data processing for oil and gas exploration. The accuracy of migration directly affects the accuracy of subsequent oil and gas reservoir characterization. Reverse-time migration is one of the most accurate migration methods at present. Multi-wave and multicomponent seismic data contain more P- and S-wave information. Making full use of multi-wave and multicomponent seismic data can offer more information about underground structure and lithology, as well as improve the accuracy of seismic exploration. Elastic reverse-time migration (ERTM) has no dip restriction and can be applied to image multi-wave and multicomponent seismic data in complex structural areas and some special lithology structures. However, the surface topography of complex regions has an influence on wavefield and seriously degrades the quality of ERTM’s migration results. We developed a new ERTM method to migrate multi-wave and multicomponent seismic data in the region with complex surface topography. We first fill the layers between the highest and lowest undulating surface with near-surface elastic parameters in a complex topography model to obtain a new model with a horizontal surface. This allows the finite difference (FD) method based on the regular rectangular grid to be used to numerically solve elastic wave equations in the model with complex topography. The decoupled wave equations are used to generate source P- and S-waves and receiver P- and S-waves to reduce crosstalk artefacts in ERTM. A topography-related filter is further used to remove the influence of surface topography on migration results. The scalar imaging condition is also applied to generate PP and PS migration images. Some numerical examples with different complex topographies demonstrate that our proposed ERTM method can remove the influence of complex topography on ERTM’s images and effectively generate high-quality ERTM images.

Polarity reversal is a well-known problem in elastic reverse time migration, and it is closely related to the imaging conditions. The dot product of source and receiver wavefields is a stable and efficient way to construct scalar imaging conditions for decomposed elastic vector wavefields. However, for PP images, the dot product introduces an angle-dependent factor that will change the polarity of image amplitudes at large opening angles, and it is also contaminated by low-wavenumber artifacts when sharp contrasts exist in the velocity model. Those two problems can be suppressed by muting the reflections with large opening angles at the expense of losing useful information. We have developed an elastic inverse-scattering imaging condition that can retain the initial polarity of the image amplitude and significantly reduce the low-wavenumber noise. For PS images, much attention is paid to the polarity-reversal problem at the normal incidence, and the dot-product-based imaging condition successfully avoids this kind of polarity reversal. There is another polarity-reversal problem arising from the sign change of the PS reflection coefficient at the Brewster angle. However, this sign change is often neglected in the construction of a stacked PS image, which will lead to reversed or distorted phases after stacking. We suggested using the S-wave impedance kernel used in elastic full-waveform inversion but only in the PS mode as an alternative to the dot-product imaging condition to alleviate this kind of polarity-reversal problem. In addition to dot-product-based imaging conditions, we analytically compare divergence- and curl-based imaging conditions and the elastic energy norm-based imaging condition with the presented imaging conditions to identify their advantages and weaknesses. Two numerical examples on a two-layer model and the SEAM 2D model are used to illustrate the effectiveness and advantages of the presented imaging conditions in suppressing low-wavenumber noise and correcting the polarity-reversal problem.

Non-artifact vector P- and S-wave separation for elastic reverse time migration

Summary Borehole dipole shear reflection imaging has attracted widespread attention in the evaluation of heterogeneous reservoirs in recent years due to its good performance in reflection amplitude, detecting depth and directivity. However, the conventional migration method based on oneway acoustic equation has been proved to be improper for the wave excited by dipole source in this paper. An elastic wave reverse time migration (RTM) is developed for more exactly imaging of dipole reflected wave to find the location of effective reservoir. The problem of storage in RTM is solved by a boundary storage strategy improved for borehole environment in our research, which has been proved to be very effective by the simulated data. The processing of 2D simulated data of formation models with cave, crack and formation interface indicate the good performance of this migration method. The migration results of the field data with this method agree well with the borehole electric imaging logging data. Furthermore, our researches reveal that the effects of the unknown formation parameters, shielding interface of former reflector, would make the interpretation of the migration image difficult.

The conventional reverse time migration utilizes regularly sampled computational grids to simulate wave propagation. Selecting the appropriate grid sampling is important for computational accuracy and efficiency. In general, the uniform-size grid cannot represent the complexity of the geology well. The grid may appear sparse in the low-velocity zone, especially in shallow depths where dispersion may occur. Conversely, it may appear excessively dense in the high-velocity zone, such as at greater depths or within a salt body, which results in higher computational memory and time consumption. To overcome these issues, we developed an efficient and accurate adaptive variable grid discretization method that automatically selects the vertical grid size based on the velocity, depth, and dominant frequency of the wavelet in elastic medium. Then we reformulated the elastic equations based on the adaptive variable grid by introducing a mapping relationship. To test the effectiveness, accuracy, and efficiency of the equation, we implemented it to both the forward propagation and migration of elastic wavefield. Synthetic numerical examples demonstrate that our proposed method can achieve elastic wavefield separation and no significant dispersion phenomenon. The multi-component imaging accuracy of reverse time migration is nearly equivalent to the traditional method, while significantly improving computational efficiency and saving storage space.

Multicomponent joint elastic reverse time migration is one of important parts in multicomponent seismic data processing. We analyse the implementation of elastic reverse time migration and the polarity reversal in PS and SP images. We present a method using Poynting vector to correct the polarity reversal in elastic reverse time migration. The polarity reversal correction using Poynting vector can be done with the imaging conditions at every time step of the elastic RTM. And the whole correction procedure can be carried out in common shot domain rather than angle domain, which is low cost and convenient. Examples show this method can accurately locate the normal incidence point and correct the polarity reversal for both simple and complex models.

Using compact representations of Green’s functions we derive a common framework for reverse time migration (RTM) and Kirchhoff migration. These compact Green’s functions (CGFs) are 3D volumes containing traveltimes and amplitudes for the N most representative events in the upcoming/downgoing decomposed 4D wavefields originating from a point source. Within this framework, we implement an RTM algorithm using a multivalued excitation time/amplitude imaging condition. This new approach produces four complementary imaging volumes (different combinations of source and receiver decomposed wavefields) and angle/azimuth gathers with computational effort less than 15% greater than that of plain (one image and no gathers) RTM algorithms. The advantages of separating the image volume into four complementary volumes are well established in the literature (low-frequency noise separation and turning wave imaging); however, its use has been limited by the computational cost. Despite using two source propagations to decompose the source wavefield, we reduce the computations to less than 20% of a single source propagation by performing finite-difference propagation with half the frequency limit used in the receiver wavefield propagation. The combination of CGF and an excitation time/amplitude imaging condition allows receiver wavefield decomposition with only one wavefield propagation. Our RTM algorithm constructs angle/azimuth gathers using a postmigration computation of the source and receiver wavefield’s propagation directions. To compute the propagation directions after migration, we use a new concept: the cumulative wavefield volumes, which are 3D, imaging-condition-guided compressions, of the 4D source and receiver wavefields. We also use CGF to implement a Kirchhoff migration algorithm that produces four complementary image volumes with RTM-like quality. Furthermore, we present synthetic and field data examples to clarify the new concepts and illustrate the results obtained using these methods.

Five alternative algorithms were evaluated to circumvent the excessive storage requirement imposed by saving source wavefield snapshots used for the crosscorrelation image condition in 2D prestack elastic reverse time migration. We compared the algorithms on the basis of their ability, either to accurately reconstruct (not save) the source wavefield or to use an alternate image condition so that neither saving nor reconstruction of full wavefields was involved. The comparisons were facilitated by using the same (velocity-stress) extrapolator in all the algorithms, and running them all on the same hardware. We assumed that there was enough memory in a node to do an extrapolation, and that all input data were stored on disk rather than residing in random-access memory. This should provide a fair and balanced comparison. Reconstruction of the source wavefield from boundary and/or initial values reduced the required storage to a very small fraction of that needed to store source wavefield snapshots for conventional crosscorrelation, at the cost of adding an additional source extrapolation. Reverse time checkpointing avoided recursive forward recomputation. Two nonreconstructive imaging conditions do not require full snapshot storage or an additional extrapolation. Time-binning the imaging criteria removed the need for image time searching or sorting. Numerical examples using elastic data from the Marmousi2 model showed that the quality of the elastic prestack PP and PS images produced by the cost-optimized alternative algorithms were (virtually) identical to the higher cost images produced by traditional crosscorrelation.

Elastic wave extrapolation in the time domain is significant for an elastic wave equation-based processing. To improve the simulation reliability and accuracy of decoupled elastic P- and S- waves, we propose the staggered-grid lowrank finite-difference method based on the elastic wave decomposition. For elastic wave propagation, a lowrank finite-difference method based on the staggered grid is derived to improve the accuracy. Regarding the application of the decoupled elastic wave equation, we derive the finite-difference scheme coefficients which are dependent on velocity. Based on the elastic wave decomposition and plane wave theories, we formulate the elastic wave-extrapolation operators, which contain trigonometric adjustment factors. Accordingly, by applying the lowrank method to approximating the operators, the finite- difference scheme is designed to discretize the decoupled wave equation. The derivation processing implies the combination of elastic wave-mode decomposition and extrapolation. The proposed method enables elastic P- and S-waves to extrapolate in the time-space domain separately and produces accurate P-and S-wave components simultaneously. Dispersion analysis suggests that our proposed method is reliable and accurate in a wide range of wavenumber. Numerical simulation tests on a simple model and the Marmousi2 model validate the accuracy and effectiveness of the method, showing its ability in handling complex structures. Although the operators are accurate only when the medium is homogeneous, they are of high accuracy when the velocity gradient is quite small and are applicable when the velocity gradient is large. The subsequent results of reverse time migration for the Marmousi2 model also suggest that the proposed method is enough to serve as an extrapolator in elastic reverse time migration.

Studies have been conducted to demonstrate the impact of correct TTI parameterization on acoustic RTM (reverse time migration) imaging results. In this study, we extend the study to investigate the TTI effect on the elastic reverse time migration. We expect to see improvement in terms of better kinematics positioning as well as better focusing by including proper TTI parameterization.