Tomographic velocity analysis and wave-equation depth migration in an overthrust terrain: A case study from the Tuha Basin, China

2018 ◽  
Vol 6 (1) ◽  
pp. T1-T13
Author(s):  
Bin Lyu ◽  
Qin Su ◽  
Kurt J. Marfurt
Keyword(s):  
Wave Equation ◽  
Data Quality ◽  
Model Building ◽  
Velocity Model ◽  
Lateral Velocity ◽  
Near Surface ◽  
Depth Migration ◽  
Time Migration ◽  
Head Waves ◽  
Velocity Variations

Although the structures associated with overthrust terrains form important targets in many basins, accurately imaging remains challenging. Steep dips and strong lateral velocity variations associated with these complex structures require prestack depth migration instead of simpler time migration. The associated rough topography, coupled with older, more indurated, and thus high-velocity rocks near or outcropping at the surface often lead to seismic data that suffer from severe statics problems, strong head waves, and backscattered energy from the shallow section, giving rise to a low signal-to-noise ratio that increases the difficulties in building an accurate velocity model for subsequent depth migration. We applied a multidomain cascaded noise attenuation workflow to suppress much of the linear noise. Strong lateral velocity variations occur not only at depth but near the surface as well, distorting the reflections and degrading all deeper images. Conventional elevation corrections followed by refraction statics methods fail in these areas due to poor data quality and the absence of a continuous refracting surface. Although a seismically derived tomographic solution provides an improved image, constraining the solution to the near-surface depth-domain interval velocities measured along the surface outcrop data provides further improvement. Although a one-way wave-equation migration algorithm accounts for the strong lateral velocity variations and complicated structures at depth, modifying the algorithm to account for lateral variation in illumination caused by the irregular topography significantly improves the image, preserving the subsurface amplitude variations. We believe that our step-by-step workflow of addressing the data quality, velocity model building, and seismic imaging developed for the Tuha Basin of China can be applied to other overthrust plays in other parts of the world.

Depth migration by the Gaussian beam summation method

Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. S81-S93 ◽  
Author(s):  
Mikhail M. Popov ◽  
Nikolay M. Semtchenok ◽  
Peter M. Popov ◽  
Arie R. Verdel
Keyword(s):  
Wave Equation ◽  
Model Building ◽  
Velocity Structure ◽  
Velocity Model ◽  
Summation Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Depth Migration ◽  
Time Migration ◽  
Ray Methods

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.

Overthrust imaging with tomo‐datuming: A case study

Geophysics ◽  
1998 ◽  
Vol 63 (1) ◽  
pp. 25-38 ◽  
Author(s):  
Xianhuai Zhu ◽  
Burke G. Angstman ◽  
David P. Sixta
Keyword(s):  
Velocity Model ◽  
Surface Velocity ◽  
Time Shift ◽  
Lateral Velocity ◽  
Prestack Depth Migration ◽  
Near Surface ◽  
Depth Migration ◽  
Static Corrections ◽  
Velocity Variations ◽  
Kirchhoff Method

Through the use of iterative turning‐ray tomography followed by wave‐equation datuming (or tomo‐datuming) and prestack depth migration, we generate accurate prestack images of seismic data in overthrust areas containing both highly variable near‐surface velocities and rough topography. In tomo‐datuming, we downward continue shot records from the topography to a horizontal datum using velocities estimated from tomography. Turning‐ray tomography often provides a more accurate near‐surface velocity model than that from refraction statics. The main advantage of tomo‐datuming over tomo‐statics (tomography plus static corrections) or refraction statics is that instead of applying a vertical time‐shift to the data, tomo‐datuming propagates the recorded wavefield to the new datum. We find that tomo‐datuming better reconstructs diffractions and reflections, subsequently providing better images after migration. In the datuming process, we use a recursive finite‐difference (FD) scheme to extrapolate wavefield without applying the imaging condition, such that lateral velocity variations can be handled properly and approximations in traveltime calculations associated with the raypath distortions near the surface for migration are avoided. We follow the downward continuation step with a conventional Kirchhoff prestack depth migration. This results in better images than those migrated from the topography using the conventional Kirchhoff method with traveltime calculation in the complicated near surface. Since FD datuming is only applied to the shallow part of the section, its cost is much less than the whole volume FD migration. This is attractive because (1) prestack depth migration usually is used iteratively to build a velocity model, so both efficiency and accuracy are important factors to be considered; and (2) tomo‐datuming can improve the signal‐to‐noise (S/N) ratio of prestack gathers, leading to more accurate migration velocity analysis and better images after depth migration. Case studies with synthetic and field data examples show that tomo‐datuming is especially helpful when strong lateral velocity variations are present below the topography.

scholarly journals Wave-equation time migration

Geophysics ◽  
2020 ◽  
pp. 1-58
Author(s):  
Sergey Fomel ◽  
Harpreet Kaur
Keyword(s):  
Wave Equation ◽  
Time Domain ◽  
Field Data ◽  
Model Building ◽  
Lateral Velocity ◽  
Depth Migration ◽  
Time Migration ◽  
Expensive Model ◽  
Approximate Equations ◽  
Ray Coordinates

Time migration, as opposed to depth migration, suffers from two well-known shortcomings: (1)approximate equations are used for computing Green’s functions inside the imaging operator; (2) in case of lateral velocity variations, the transformation between the image ray coordinates andthe Cartesian coordinates is undefined in places where the image rays cross. We show that thefirst limitation can be removed entirely by formulating time migration through wave propagationin image-ray coordinates. The proposed approach constructs a time-migrated image without relyingon any kind of traveltime approximation by formulating an appropriate geometrically accurateacoustic wave equation in the time-migration domain. The advantage of this approach is that thepropagation velocity in image-ray coordinates does not require expensive model building and canbe approximated by quantities that are estimated in conventional time-domain processing. Synthetic and field data examples demonstrate the effectiveness of the proposed approach and show that theproposed imaging workflow leads to a significant uplift in terms of image quality and can bridge thegap between time and depth migrations. The image obtained by the proposed algorithm is correctlyfocused and mapped to depth coordinates it is comparable to the image obtained by depth migration.

3-D prestack Kirchhoff depth migration: From prototype to production in a massively parallel processor environment

Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 546-556 ◽  
Author(s):  
Herman Chang ◽  
John P. VanDyke ◽  
Marcelo Solano ◽  
George A. McMechan ◽  
Duryodhan Epili
Keyword(s):  
Model Building ◽  
Velocity Model ◽  
Production Scale ◽  
Data Set ◽  
Depth Migration ◽  
Prestack Migration ◽  
Prototype Development ◽  
Time Migration ◽  
Intel Paragon ◽  
Kirchhoff Depth Migration

Portable, production‐scale 3-D prestack Kirchhoff depth migration software capable of full‐volume imaging has been successfully implemented and applied to a six‐million trace (46.9 Gbyte) marine data set from a salt/subsalt play in the Gulf of Mexico. Velocity model building and updates use an image‐driven strategy and were performed in a Sun Sparc environment. Images obtained by 3-D prestack migration after three velocity iterations are substantially better focused and reveal drilling targets that were not visible in images obtained from conventional 3-D poststack time migration. Amplitudes are well preserved, so anomalies associated with known reservoirs conform to the petrophysical predictions. Prototype development was on an 8-node Intel iPSC860 computer; the production version was run on an 1824-node Intel Paragon computer. The code has been successfully ported to CRAY (T3D) and Unix workstation (PVM) environments.

Response by Arthur E. Barnes to the reply by the authors

Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1947-1947 ◽  
Author(s):  
Arthur E. Barnes
Keyword(s):  
Lateral Velocity ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Time Migration ◽  
Pulse Distortion ◽  
Velocity Variations ◽  
Zero Offset ◽  
Complex Geology

I appreciate the thoughtful and thorough response given by Tygel et al. They point out that even for a single dipping reflector imaged by a single non‐zero offset raypath, pulse distortion caused by “standard processing” (NM0 correction‐CMP sort‐stack‐time migration) and pulse distortion caused by prestack depth migration are not really the same, because the reflecting point is mispositioned in standard processing. Within a CMP gather, this mispositioning increases with offset, giving rise to “CMP smear.” CMP smear degrades the stack, introducing additional pulse distortion. Where i‐t is significant, and where lateral velocity variations or reflection curvature are large, such as for complex geology, the pulse distortion of standard processing can differ greatly from that of prestack depth migration.

Hybrid migration: A cost‐effective 3-D depth‐imaging technique

Geophysics ◽  
1997 ◽  
Vol 62 (2) ◽  
pp. 568-576 ◽  
Author(s):  
Young C. Kim ◽  
Worth B. Hurt, ◽  
Louis J. Maher ◽  
Patrick J. Starich
Keyword(s):  
Model Building ◽  
Cost Effective ◽  
Velocity Model ◽  
Depth Image ◽  
Hybrid Technique ◽  
Prestack Depth Migration ◽  
Depth Imaging ◽  
Depth Migration ◽  
Time Migration ◽  
Prestack Time Migration

The transformation of surface seismic data into a subsurface image can be separated into two components—focusing and positioning. Focusing is associated with ensuring the data from different offsets are contributing constructively to the same event. Positioning involves the transformation of the focused events into a depth image consistent with a given velocity model. In prestack depth migration, both of these operations are achieved simultaneously; however, for 3-D data, the cost is significant. Prestack time migration is much more economical and focuses events well even in the presence of moderate velocity variations, but suffers from mispositioning problems. Hybrid migration is a cost‐effective depth‐imaging approach that uses prestack time migration for focusing; inverse migration for the removal of positioning errors; and poststack depth migration for proper positioning. When lateral velocity changes are moderate, the hybrid technique can generate a depth image that is consistent with a velocity field. For very complex structures that require prestack depth migration, the results of the hybrid technique can be used to create a starting velocity model, thereby reducing the number of iterations for velocity model building.

Salt interpretation enabled by reverse-time migration

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE211-VE216 ◽  
Author(s):  
Jacobus Buur ◽  
Thomas Kühnel
Keyword(s):  
Model Building ◽  
Velocity Model ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Seismic Line ◽  
Depth Migration ◽  
Time Migration ◽  
Velocity Model Building ◽  
Migration Method

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.

3D angle gathers from reverse time migration

Geophysics ◽  
2011 ◽  
Vol 76 (2) ◽  
pp. S77-S92 ◽  
Author(s):  
Sheng Xu ◽  
Yu Zhang ◽  
Bing Tang
Keyword(s):  
Wave Equation ◽  
Model Building ◽  
Velocity Model ◽  
Azimuth Angle ◽  
Three Dimensions ◽  
Practical Implementation ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Reflection Angle

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.

Datum correction by wave‐equation extrapolation

Geophysics ◽  
1988 ◽  
Vol 53 (10) ◽  
pp. 1311-1322 ◽  
Author(s):  
V. Shtivelman ◽  
A. Canning
Keyword(s):  
Wave Equation ◽  
Synthetic Data ◽  
Real Data ◽  
Static Solution ◽  
Velocity Model ◽  
Lateral Velocity ◽  
Static Correction ◽  
Slow Velocity ◽  
Velocity Variations ◽  
Seismic Sections

Seismic sections are usually datum corrected by static shifting. For small differences in elevation and slow velocity variations between the input datum and the output datum, static shifting is a sufficiently accurate datum correction procedure. However, for significant differences in elevations and a more complicated velocity model, the accuracy of the static solution may prove to be insufficient; and a more exact method should be used. In this paper, we study the limitations of the static method of datum correction and develop simple and effective extrapolation schemes based on the wave equation, schemes which lead to more accurate datum correction. The distortions of seismic events caused by static correction are illustrated by a number of simple examples. To reduce the distortions, we propose a number of extrapolation schemes based on the asymptories of the Kirchhoff integral solution of the 2-D scalar wave equation. Application of the extrapolation algorithms to synthetic data shows that they provide accurate datum corrections even for a nonplanar input datum and vertical and lateral velocity variations. The algorithms have been successfully applied to real data.

Seismic wave-equation demigration/migration

Geophysics ◽  
2010 ◽  
Vol 75 (3) ◽  
pp. S111-S119 ◽  
Author(s):  
Hervé Chauris ◽  
Mondher Benjemaa
Keyword(s):  
Wave Equation ◽  
Source Term ◽  
Model Building ◽  
Synthetic Data ◽  
Velocity Model ◽  
Single Scattering ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Velocity Models ◽  
Time Migration

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.

Sign in / Sign up

Export Citation Format

Share Document