A depth migration method based on the full‐wave reverse‐time calculation and local one‐way propagation

Author(s):  
Xiao‐Bi Xie ◽  
Ru‐Shan Wu
Geophysics ◽  
1994 ◽  
Vol 59 (4) ◽  
pp. 597-609 ◽  
Author(s):  
Wen‐Fong Chang ◽  
George A. McMechan

By combining and extending previous algorithms for 2-D prestack elastic migration and 3-D prestack acoustic migration, a full 3-D elastic prestack depth migration algorithm is developed. Reverse‐time extrapolation of the recorded data is by 3-D elastic finite differences; computation of the image time for each point in the 3-D volume is by 3-D acoustic finite differences. The algorithm operates on three‐component, vector‐wavefield common‐source data and produces three‐component vector reflectivity distributions. Converted P‐to‐S reflections are automatically imaged with the primary P‐wave reflections. There are no dip restrictions as the full wave equation is used. The algorithm is illustrated by application to synthetic data from three models; a flat reflector, a dipping truncated wedge overlying a flat reflector, and the classical French double dome and fault model.


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 ◽  
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 ◽  
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.


2021 ◽  
Vol 11 (7) ◽  
pp. 3010
Author(s):  
Hao Liu ◽  
Xuewei Liu

The lack of an initial condition is one of the major challenges in full-wave-equation depth extrapolation. This initial condition is the vertical partial derivative of the surface wavefield and cannot be provided by the conventional seismic acquisition system. The traditional solution is to use the wavefield value of the surface to calculate the vertical partial derivative by assuming that the surface velocity is constant. However, for seismic exploration on land, the surface velocity is often not uniform. To solve this problem, we propose a new method for calculating the vertical partial derivative from the surface wavefield without making any assumptions about the surface conditions. Based on the calculated derivative, we implemented a depth-extrapolation-based full-wave-equation migration from topography using the direct downward continuation. We tested the imaging performance of our proposed method with several experiments. The results of the Marmousi model experiment show that our proposed method is superior to the conventional reverse time migration (RTM) algorithm in terms of imaging accuracy and amplitude-preserving performance at medium and deep depths. In the Canadian Foothills model experiment, we proved that our method can still accurately image complex structures and maintain amplitude under topographic scenario.


Geophysics ◽  
1983 ◽  
Vol 48 (11) ◽  
pp. 1514-1524 ◽  
Author(s):  
Edip Baysal ◽  
Dan D. Kosloff ◽  
John W. C. Sherwood

Migration of stacked or zero‐offset sections is based on deriving the wave amplitude in space from wave field observations at the surface. Conventionally this calculation has been carried out through a depth extrapolation. We examine the alternative of carrying out the migration through a reverse time extrapolation. This approach may offer improvements over existing migration methods, especially in cases of steeply dipping structures with strong velocity contrasts. This migration method is tested using appropriate synthetic data sets.


2003 ◽  
Vol 46 (6) ◽  
pp. 1176-1185 ◽  
Author(s):  
Shengchang CHEN ◽  
Jingzhong CAO ◽  
Zaitian MA

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. Q27-Q37
Author(s):  
Yang Shen ◽  
Jie Zhang

Refraction methods are often applied to model and image near-surface velocity structures. However, near-surface imaging is very challenging, and no single method can resolve all of the land seismic problems across the world. In addition, deep interfaces are difficult to image from land reflection data due to the associated low signal-to-noise ratio. Following previous research, we have developed a refraction wavefield migration method for imaging shallow and deep interfaces via interferometry. Our method includes two steps: converting refractions into virtual reflection gathers and then applying a prestack depth migration method to produce interface images from the virtual reflection gathers. With a regular recording offset of approximately 3 km, this approach produces an image of a shallow interface within the top 1 km. If the recording offset is very long, the refractions may follow a deep path, and the result may reveal a deep interface. We determine several factors that affect the imaging results using synthetics. We also apply the novel method to one data set with regular recording offsets and another with far offsets; both cases produce sharp images, which are further verified by conventional reflection imaging. This method can be applied as a promising imaging tool when handling practical cases involving data with excessively weak or missing reflections but available refractions.


Sign in / Sign up

Export Citation Format

Share Document