Salt interpretation enabled by reverse-time migration

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE211-VE216 ◽  
Author(s):  
Jacobus Buur ◽  
Thomas Kühnel

Many production targets in greenfield exploration are found in salt provinces, which have highly complex structures as a result of salt formation over geologic time. Difficult geologic settings, steep dips, and other wave-propagation effects make reverse-time migration (RTM) the migration method of choice, rather than Kirchhoff migration or other (by definition approximate) one-way equation methods. Imaging of the subsurface using any depth-migration algorithm can be done successfully only when the quality of the prior velocity model is sufficient. The (velocity) model-building loop is an iterative procedure for improving the velocity model. This is done by obtaining certain measurements (residual moveout) on image gathers generated during the migration procedure; those measurements then are input into tomographic updating. Commonly RTM is applied around salt bodies, where building the velocity model fails essentially because tomography is ray-trace based. Our idea is to apply RTM directly inside the model-building loop but to do so without using the image gathers. Although the process is costly, we migrate the full frequency content of the data to create a high-quality stack. This enhances the interpretation of top and bottom salt significantly and enables us to include the resulting salt geometry in the velocity model properly. We demonstrate our idea on a 2D West Africa seismic line. After several model-building iterations, the result is a dramatically improved velocity model. With such a good model as input, the final RTM confirms the geometry of the salt bodies and basically the salt interpretation, and yields a compelling image of the subsurface.

Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. S81-S93 ◽  
Author(s):  
Mikhail M. Popov ◽  
Nikolay M. Semtchenok ◽  
Peter M. Popov ◽  
Arie R. Verdel

Seismic depth migration aims to produce an image of seismic reflection interfaces. Ray methods are suitable for subsurface target-oriented imaging and are less costly compared to two-way wave-equation-based migration, but break down in cases when a complex velocity structure gives rise to the appearance of caustics. Ray methods also have difficulties in correctly handling the different branches of the wavefront that result from wave propagation through a caustic. On the other hand, migration methods based on the two-way wave equation, referred to as reverse-time migration, are known to be capable of dealing with these problems. However, they are very expensive, especially in the 3D case. It can be prohibitive if many iterations are needed, such as for velocity-model building. Our method relies on the calculation of the Green functions for the classical wave equation by per-forming a summation of Gaussian beams for the direct and back-propagated wavefields. The subsurface image is obtained by cal-culating the coherence between the direct and backpropagated wavefields. To a large extent, our method combines the advantages of the high computational speed of ray-based migration with the high accuracy of reverse-time wave-equation migration because it can overcome problems with caustics, handle all arrivals, yield good images of steep flanks, and is readily extendible to target-oriented implementation. We have demonstrated the quality of our method with several state-of-the-art benchmark subsurface models, which have velocity variations up to a high degree of complexity. Our algorithm is especially suited for efficient imaging of selected subsurface subdomains, which is a large advantage particularly for 3D imaging and velocity-model refinement applications such as subsalt velocity-model improvement. Because our method is also capable of providing highly accurate migration results in structurally complex subsurface settings, we have also included the concept of true-amplitude imaging in our migration technique.


Geophysics ◽  
2020 ◽  
pp. 1-45
Author(s):  
German Garabito ◽  
Paul L. Stoffa ◽  
Yuri S. F. Bezerra ◽  
João L. Caldeira

The application of the reverse time migration (RTM) in land seismic data is still a great challenge due to its low quality, low signal-to-noise ratio, irregular spatial sampling, acquisition gaps, missing traces, etc. Therefore, prior to the application of this kind of depth migration, the input pre-stack data must be conveniently preconditioned, that is, it must be interpolated, regularized, and enhanced. There are several methods for seismic data preconditioning, but for 2D real land data, the regularization of pre-stack data based on common reflection surface (CRS) stack method provides high quality enhanced preconditioned data, which is suitable for pre-stack depth migration and velocity model building. This work demonstrates the potential of RTM combined with CRS-based pre-stack data regularization, applied to real land seismic data with low quality and irregularly sparse spatial sampled, from geologically complex areas with the presence of diabase sills and steep dip reflections. Usually, determining the wavelet of the seismic source from land data is a challenge, because of this, RTM migration is often applied using artificial sources (e.g. Ricker). In this work, from the power spectrum of the pre-stacked data, we determine the wavelet of the seismic source to apply the RTM to real land data. We present applications of the pre-stack data preconditioning based on CRS stack and of the RTM in 2D land data of Tacutu and Parnaiba Basins, Brazil. Comparisons with the standard Kirchhoff depth migration reveals that the RTM improves the quality and resolution of the migrated images.


Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB27-WB39 ◽  
Author(s):  
Zheng-Zheng Zhou ◽  
Michael Howard ◽  
Cheryl Mifflin

Various reverse time migration (RTM) angle gather generation techniques have been developed to address poor subsalt data quality and multiarrival induced problems in gathers from Kirchhoff migration. But these techniques introduce new problems, such as inaccuracies in 2D subsurface angle gathers and edge diffraction artifacts in 3D subsurface angle gathers. The unique rich-azimuth data set acquired over the Shenzi field in the Gulf of Mexico enabled the generally artifact-free generation of 3D subsurface angle gathers. Using this data set, we carried out suprasalt tomography and salt model building steps and then produced 3D angle gathers to update the subsalt velocity. We used tilted transverse isotropy RTM with extended image condition to generate full 3D subsurface offset domain common image gathers, which were subsequently converted to 3D angle gathers. The angle gathers were substacked along the subsurface azimuth axis into azimuth sectors. Residual moveout analysis was carried out, and ray-based tomography was used to update velocities. The updated velocity model resulted in improved imaging of the subsalt section. We also applied residual moveout and selective stacking to 3D angle gathers from the final migration to produce an optimized stack image.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. S469-S475 ◽  
Author(s):  
Carlos Alberto da Costa Filho ◽  
Andrew Curtis

The objective of prestack depth migration is to position reflectors at their correct subsurface locations. However, migration methods often also generate artifacts along with physical reflectors, which hamper interpretation. These spurious reflectors often appear at different spatial locations in the image depending on which migration method is used. Therefore, we have devised a postimaging filter that combines two imaging conditions to preserve their similarities and to attenuate their differences. The imaging filter is based on combining the two constituent images and their envelopes that were obtained from the complex vertical traces of the images. We have used the method to combine two images resulting from different migration schemes, which produce dissimilar artifacts: a conventional migration method (equivalent to reverse time migration) and a deconvolution-based imaging method. We show how this combination may be exploited to attenuate migration artifacts in a final image. A synthetic model containing a syncline and stochastically generated small-scale heterogeneities in the velocity and density distributions was used for the numerical example. We compared the images in detail at two locations where spurious events arose and also at a true reflector. We found that the combined imaging condition has significantly fewer artifacts than either constituent image individually.


Geophysics ◽  
2011 ◽  
Vol 76 (2) ◽  
pp. S77-S92 ◽  
Author(s):  
Sheng Xu ◽  
Yu Zhang ◽  
Bing Tang

Common-image gathers are an important output of prestack depth migration. They provide information needed for velocity model building and amplitude and phase information for subsurface attribute interpretation. Conventionally, common-image gathers are computed using Kirchhoff migration on common-offset/azimuth data volumes. When geologic structures are complex and strong contrasts exist in the velocity model, the complicated wave behaviors will create migration artifacts in the image gathers. As long as the gather output traces are indexed by any surface attribute, such as source location, receiver location, or surface plane-wave direction, they suffer from the migration artifacts caused by multiple raypaths. These problems have been addressed in a significant amount of work, resulting in common-image gathers computed in the reflection angle domain, whose traces are indexed by the subsurface reflection angle and/or the subsurface azimuth angle. Most of these efforts have concentrated on Kirchhoff and one-way wave-equation migration methods. For reverse time migration, subsurface angle gathers can be produced using the same approach as that used for one-way wave-equation migration. However, these approaches need to be revisited when producing high-quality subsurface angle gathers in three dimensions (reflection angle/azimuth angle), especially for wide-azimuth data. We have developed a method for obtaining 3D subsurface reflection angle/azimuth angle common-image gathers specifically for the amplitude-preserved reverse time migration. The method builds image gathers with a high-dimensional convolution of wavefields in the wavenumber domain. We have found a windowed antileakage Fourier transform method that leads to an efficient and practical implementation. This approach has generated high-resolution angle-domain gathers on synthetic 2.5D data and 3D wide-azimuth real data.


Geophysics ◽  
2010 ◽  
Vol 75 (3) ◽  
pp. S111-S119 ◽  
Author(s):  
Hervé Chauris ◽  
Mondher Benjemaa

Reverse-time migration is a well-known method based on a single-scattering approximation; it is designed to obtain seismic images in the case of a complex subsurface. It can, however, be a very time-consuming task because the number of computations is directly proportional to the number of processed sources. In the context of velocity model-building, iterative approaches require that one derives a series of migrated sections for different velocity models. We propose to replace the summation over sources by a summation over depth offsets or time delays defined in the subsurface. For that, we have developed a new relationship between two migrated sections obtained for two different velocity models. Starting from one of the two images, we obtain a second section correctly and efficiently. For each time delay, we compute a generalized source term by extending the concept of exploding reflector to nonzero offset. We obtain the final migrated section by solving the same wave equation in the perturbed model with the modified source term. Our work included testing the methodology on 2D synthetic data sets, particularly when the initial and perturbed velocity models differ greatly.


Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB169-WB174 ◽  
Author(s):  
Shuo Ji ◽  
Tony Huang ◽  
Kang Fu ◽  
Zhengxue Li

For deep-water Gulf of Mexico, accurate salt geometry is critical to subsalt imaging. This requires the definition of both external and internal salt geometries. In recent years, external salt geometry (i.e., boundaries between allochthonous salt and background sediment) has improved a great deal due to advances in acquisition, velocity model building, and migration algorithms. But when it comes to defining internal salt geometry (i.e., intrasalt inclusions or dirty salt), no efficient method has yet been developed. In common industry practices, intrasalt inclusions (and thus their velocity anomalies) are generally ignored during the model building stages. However, as external salt geometries reach higher levels of accuracy, it becomes more important to consider the once-ignored effects of dirty salt. We have developed a reflectivity-based approach for dirty salt velocity inversion. This method takes true-amplitude reverse time migration stack volumes as input, then estimates the dirty salt velocity based on reflectivity under a 1D assumption. Results from a 2D synthetic data set and a real 3D Wide Azimuth data set demonstrated that the reflectivity inversion scheme significantly improves the subsalt image for certain areas. In general, we believe that this method produces a better salt model than the traditional clean salt velocity approach.


Geophysics ◽  
1989 ◽  
Vol 54 (3) ◽  
pp. 309-318 ◽  
Author(s):  
Moshe Reshef ◽  
David Kessler

This work deals with the practical aspects of three‐dimensional (3-D) poststack depth migration. A method, based on depth extrapolation in the frequency domain, is used for the migration. This method is suitable for structures with arbitrary velocity variation, and the number of computations required can be directly related to the complexity of the given velocity function. We demonstrate the superior computational efficiency of this method for 3-D depth migration relative to the reverse‐time migration method. The computational algorithm used for the migration is designed for a multi‐processor machine (Cray-XMP/48) and takes advantage of advanced disk technologies to overcome the input/output (I/O) problem. The method is demonstrated with both synthetic and field data. The migration of a typical 3-D data volume can be accomplished in only a few hours.


2020 ◽  
Vol 17 (3) ◽  
pp. 525-538 ◽  
Author(s):  
Yimin Sun ◽  
Young Seo Kim ◽  
Shan Qu ◽  
Eric Verschuur

Abstract Joint migration inversion is a recently proposed technology, accommodating velocity model building and seismic migration in one integrated process. Different from the widely accepted full waveform inversion technology, it uses imaging parameters, i.e. velocities and reflectivities of the subsurface, to parameterize its solution space. The unique feature of this new technology is its explicit capability to exploit multiples in its inversion scheme, which are treated as noise by most current technologies. In this paper, we comprehensively evaluate the state-of-the-art joint migration inversion technology from various angles: we first benchmark its performance, on both velocity model building and seismic imaging, against that of the well-accepted workflow comprising full waveform inversion and reverse-time migration using a fully controlled 2D realistic synthetic dataset. Next, we demonstrate its application on a 2D field dataset. Last, we use another 2D synthetic dataset to clearly illustrate the challenges the current joint migration inversion technology is facing. With this paper, we transparently reveal the pros of cons of the current joint migration inversion, and we will also point out the imminent research directions joint migration inversion technology should focus on in the next phase for it to be more widely accepted by the geophysics community.


Sign in / Sign up

Export Citation Format

Share Document