Prestack migration by split‐step DSR

Geophysics ◽  
1996 ◽  
Vol 61 (5) ◽  
pp. 1412-1416 ◽  
Author(s):  
Alexander Mihai Popovici
Keyword(s):  
Phase Shift ◽  
Square Root ◽  
Variable Velocity ◽  
Step Method ◽  
Lateral Velocity ◽  
Data Set ◽  
Prestack Migration ◽  
Imaging Results ◽  
Velocity Variations ◽  
Migration Equation

The double‐square‐root (DSR) prestack migration equation, though defined for depth variable velocity, can be used to image media with strong velocity variations using a phase‐shift plus interpolation (PSPI) or split‐step correction. The split‐step method is based on applying a phase‐shift correction to the extrapolated wavefield—a correction that attempts to compensate for the lateral velocity variations. I show how to extend DSR prestack migration to lateral velocity media and exemplify the method by applying the new algorithm to the Marmousi data set. The split‐step DSR migration is very fast and offers excellent imaging results.

The Hartley transform in seismic imaging

Geophysics ◽  
2001 ◽  
Vol 66 (4) ◽  
pp. 1251-1257 ◽  
Author(s):  
Henning Kühl ◽  
Maurico D. Sacchi ◽  
Jürgen Fertig
Keyword(s):  
Fourier Transform ◽  
Phase Shift ◽  
Shift Operator ◽  
Three Dimensions ◽  
Lateral Velocity ◽  
Data Set ◽  
Prestack Migration ◽  
Hartley Transform ◽  
Alternative Means ◽  
The Fourier Transform

Phase‐shift migration techniques that attempt to account for lateral velocity variations make substantial use of the fast Fourier transform (FFT). Generally, the Hermitian symmetry of the complex‐valued Fourier transform causes computational redundancies in terms of the number of operations and memory requirements. In practice a combination of the FFT with the well‐known real‐to‐complex Fourier transform is often used to avoid such complications. As an alternative means to the Fourier transform, we introduce the inherently real‐valued, non‐symmetric Hartley transform into phase‐shift migration techniques. By this we automatically avoid the Hermitian symmetry resulting in an optimized algorithm that is comparable in efficiency to algorithms based on the real‐to‐complex FFT. We derive the phase‐shift operator in the Hartley domain for migration in two and three dimensions and formulate phase shift plus interpolation, split‐step migration, and split‐step double‐square‐root prestack migration in terms of the Hartley transform as examples. We test the Hartley phase‐shift operator for poststack and prestack migration using the SEG/EAGE salt model and the Marmousi data set, respectively.

Seismic depth imaging with the Gabor transform

Geophysics ◽  
2008 ◽  
Vol 73 (3) ◽  
pp. S91-S97 ◽  
Author(s):  
Yongwang Ma ◽  
Gary F. Margrave
Keyword(s):  
Fourier Transform ◽  
Phase Shift ◽  
Gabor Transform ◽  
Positioning Error ◽  
Lateral Velocity ◽  
Data Set ◽  
Depth Migration ◽  
Velocity Variations ◽  
Spatial Positioning ◽  
Wavefield Extrapolation

Wavefield extrapolation in depth, a vital component of wave-equation depth migration, is accomplished by repeatedly applying a mathematical operator that propagates the wavefield across a single depth step, thus creating a depth marching scheme. The phase-shift method of wavefield extrapolation is fast and stable; however, it can be cumbersome to adapt to lateral velocity variations. We address the extension of phase-shift extrapolation to lateral velocity variations by using a spatial Gabor transform instead of the normal Fourier transform. The Gabor transform, also known as the windowed Fourier transform, is applied to the lateral spatial coordinates as a windowed discrete Fourier transform where the entire set of windows is required to sum to unity. Within each window, a split-step Fourier phase shift is applied. The most novel element of our algorithm is an adaptive partitioning scheme that relates window width to lateral velocity gradient such that the estimated spatial positioning error is bounded below a threshold. The spatial positioning error is estimated by comparing the Gabor method to its mathematical limit, called the locally homogeneous approximation — a frequency-wavenumber-dependent phase shift that changes according to the local velocity at each position. The assumption of local homogeneity means this position-error estimate may not hold strictly for large scattering angles in strongly heterogeneous media. The performance of our algorithm is illustrated with imaging results from prestack depth migration of the Marmousi data set. With respect to a comparable space-frequency domain imaging method, the proposed method improves images while requiring roughly 50% more computing time.

True-amplitude one-way wave equation migration in the mixed domain

Geophysics ◽  
2010 ◽  
Vol 75 (5) ◽  
pp. S199-S209 ◽  
Author(s):  
Flor A. Vivas ◽  
Reynam C. Pestana
Keyword(s):  
Wave Equation ◽  
Phase Shift ◽  
Numerical Experiments ◽  
Wave Equations ◽  
Laplacian Operator ◽  
Single Shot ◽  
Lateral Velocity ◽  
Data Set ◽  
Classical Formulation ◽  
Imaging Tool

One-way wave equation migration is a powerful imaging tool for locating accurately reflectors in complex geologic structures; however, the classical formulation of one-way wave equations does not provide accurate amplitudes for the reflectors. When dynamic information is required after migration, such as studies for amplitude variation with angle or when the correct amplitudes of the reflectors in the zero-offset images are needed, some modifications to the one-way wave equations are required. The new equations, which are called “true-amplitude one-way wave equations,” provide amplitudes that are equivalent to those provided by the leading order of the ray-theoretical approximation through the modification of the transverse Laplacian operator with dependence of lateral velocity variations, the introduction of a new term associated with the amplitudes, and the modification of the source representation. In a smoothly varying vertical medium,the extrapolation of the wavefields with the true-amplitude one-way wave equations simplifies to the product of two separable and commutative factors: one associated with the phase and equal to the phase-shift migration conventional and the other associated with the amplitude. To take advantage of this true-amplitude phase-shift migration, we developed the extension of conventional migration algorithms in a mixed domain, such as phase shift plus interpolation, split step, and Fourier finite difference. Two-dimensional numerical experiments that used a single-shot data set showed that the proposed mixed-domain true-amplitude algorithms combined with a deconvolution-type imaging condition recover the amplitudes of the reflectors better than conventional mixed-domain algorithms. Numerical experiments with multiple-shot Marmousi data showed improvement in the amplitudes of the deepest structures and preservation of higher frequency content in the migrated images.

Seismic depth migration with the Dirac equation

Geophysics ◽  
1999 ◽  
Vol 64 (3) ◽  
pp. 925-933 ◽  
Author(s):  
Ketil Hokstad ◽  
Rune Mittet
Keyword(s):  
Dirac Equation ◽  
Taylor Series ◽  
Rapid Expansion ◽  
Difference Operators ◽  
Square Root ◽  
Exact Linearization ◽  
Lateral Velocity ◽  
Depth Migration ◽  
Spatial Derivatives ◽  
Velocity Variations

We demonstrate the applicability of the Dirac equation in seismic wavefield extrapolation by presenting a new explicit one‐way prestack depth migration scheme. The method is in principle accurate up to 90° from the vertical, and it tolerates lateral velocity variations. This is achieved by performing the extrapolation step of migration with the Dirac equation, implemented in the space‐frequency domain. The Dirac equation is an exact linearization of the square‐root wave equation and is equivalent to keeping infinitely many terms in a Taylor series or continued‐fraction expansion of the square‐root operator. An important property of the new method is that the local velocity and the spatial derivatives decouple in separate terms within the extrapolation operator. Therefore, we do not need to precompute and store large tables of convolutional extrapolator coefficients depending on velocity. The main drawback of the explicit scheme is that evanescent energy must be removed at each depth step to obtain numerical stability. We have tested two numerical implementations of the migration scheme. In the first implementation, we perform depth stepping using the Taylor series approximation and compute spatial derivatives with high‐order finite difference operators. In the second implementation, we perform depth stepping with the Rapid expansion method and numerical differentiation with the pseudospectral method. The imaging condition is a generalization of Claerbout’s U / D principle. For both implementations, the impulse response is accurate up to 80° from the vertical. Using synthetic data from a simple fault model, we test the depth migration scheme in the presence of lateral velocity variations. The results show that the proposed migration scheme images dipping reflectors and the fault plane in the correct positions.

scholarly journals Stable wide‐angle Fourier finite‐difference downward extrapolation of 3‐D wavefields

Geophysics ◽  
2002 ◽  
Vol 67 (3) ◽  
pp. 872-882 ◽  
Author(s):  
Biondo Biondi
Keyword(s):  
Finite Difference ◽  
Continuation Method ◽  
Phase Error ◽  
Azimuthal Anisotropy ◽  
Laplacian Operator ◽  
Lateral Velocity ◽  
Data Set ◽  
Frequency Dependent ◽  
Reference Velocity ◽  
Velocity Variations

I present an unconditionally stable, implicit finite‐difference operator that corrects the constant‐velocity phase‐shift operator for lateral velocity variations. The method is based on the Fourier finite‐difference (FFD) method. Contrary to previous results, my correction operator is stable even when the medium velocity has sharp discontinuities, and the reference velocity is higher than the medium velocity. The stability of the new correction enables the definition of a new downward‐continuation method based on the interpolation of two wavefields: the first wavefield is obtained by applying the FFD correction starting from a reference velocity lower than the medium velocity; the second wavefield is obtained by applying the FFD correction starting from a reference velocity higher than the medium velocity. The proposed Fourier finite‐difference plus interpolation (FFDPI) method combines the advantages of the FFD technique with the advantages of interpolation. A simple and economical procedure for defining frequency‐dependent interpolation weight is presented. When the interpolation step is performed using these frequency‐dependent interpolation weights, it significantly reduces the residual phase error after interpolation, the frequency dispersion caused by the discretization of the Laplacian operator, and the azimuthal anisotropy caused by splitting. Tests on zero‐offset data from the SEG‐EAGE salt data set show that the FFDPI method improves the imaging of a fault reflection with respect to a similar interpolation scheme that uses a split‐step correction for adapting to lateral velocity variations.

Prestack partial migration

Geophysics ◽  
1980 ◽  
Vol 45 (12) ◽  
pp. 1753-1779 ◽  
Author(s):  
Özdoan Yilmaz ◽  
Jon F. Claerbout
Keyword(s):  
Seismic Data ◽  
Partial Migration ◽  
Velocity Variation ◽  
Square Root ◽  
Seismic Data Processing ◽  
Lateral Velocity ◽  
Data Set ◽  
Significant Term ◽  
Separable Approximation ◽  
Root Operator

Conventional seismic data processing can be improved by modifying wide‐offset data so that dipping events stack coherently. A procedure to achieve this improvement is proposed here, which is basically a “partial” migration of common offset sections prior to stack. It has an advantage over full migration before stack in that, in the case of the latter, the final product is a migrated section. However, the prestack partial migration provides the interpreter with a high‐quality common midpoint (CMP) stacked section which can be subsequently migrated. The theory of prestack partial migration is based on the double square‐root equation, which describes seismic imaging with many shots and receivers. The double square‐root operator in midpoint‐offset space can be separated approximately into two terms, one involving only migration effects and the other involving only moveout correction. This separation provides an analysis of conventional processing. Estimation of errors in the separation yields the equation for prestack partial migration. Extension of the theory for separable approximation to incorporate lateral velocity variation yields a significant term proportional to the product of the first powers of offset, dip, and lateral velocity gradient. This term was used to obtain a rough estimate of lateral velocity variation from a field data set.

scholarly journals Compensating finite‐difference errors in 3-D migration and modeling

Geophysics ◽  
1991 ◽  
Vol 56 (10) ◽  
pp. 1650-1660 ◽  
Author(s):  
Zhiming Li
Keyword(s):  
Wave Equation ◽  
Phase Shift ◽  
Finite Difference ◽  
Interpolation Method ◽  
Splitting Method ◽  
Three Dimensional ◽  
Lateral Velocity ◽  
Frequency Space ◽  
Difference Grid ◽  
Imaging Results

One‐pass three‐dimensional (3-D) depth migration potentially offers more accurate imaging results than does conventional two‐pass 3-D migration for variable velocity media. Conventional one‐pass 3-D migration, using the method of finite‐difference inline and crossline splitting, however, creates large errors in the image of complex structures. These errors are due to paraxial wave‐equation approximation of the one‐way wave equation, inline‐crossline splitting, and finite‐difference grid dispersion. To compensate for these errors, and still preserve the efficiency of the conventional finite‐difference splitting method, a phase‐correction operator is derived by minimizing the difference between the ideal 3-D migration (or modeling) and the actual, conventional 3-D migration (or modeling). For frequency‐space 3-D finite‐difference migration and modeling, the compensation operator is implemented using either the phase‐shift, or phase‐shift‐plus‐interpolation method, depending on the extent of lateral velocity variations. The compensation operator increases the accuracy of handling steep dips, suppresses the inline and crossline splitting error, and reduces finite‐difference grid dispersions.

Velocity analysis by iterative profile migration

Geophysics ◽  
1989 ◽  
Vol 54 (6) ◽  
pp. 718-729 ◽  
Author(s):  
Kamal Al‐Yahya
Keyword(s):  
Initial Velocity ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Lateral Velocity ◽  
Prestack Migration ◽  
Subsequent Step ◽  
One Step ◽  
Velocity Variations ◽  
And Migration

In conventional seismic processing, velocity analysis is performed by using the normal moveout (NMO) equation which is based on the assumption of flat, horizontal reflectors. Imaging by migration (either before or after stack) is done normally in a subsequent step using these velocities. In this paper, velocity analysis and imaging are combined in one step, and migration itself is used as a velocity indicator. Because, unlike NMO, migration can be formulated for any velocity function, migration‐based velocity analysis methods are capable of handling arbitrary structures, i.e., those with lateral velocity variations. In the proposed scheme, each shot gather (profile) is migrated with an initial depth‐velocity model. Profile migration is implemented in the (x, ω) domain, but the actual implementation of profile migration is not critical, as long as it is not done in a spatial‐wavenumber domain, which would preclude handling of lateral velocity variations. After migration with an initial velocity model, the velocity error is estimated, and the initial velocity model is updated; the process is repeated until convergence is achieved. The velocity analysis is based on the principle that after prestack migration with the correct velocity model, an image in a common‐receiver gather (CRG) is aligned horizontally regardless of structure. The deviation from horizontal alignment is therefore a measure of the error in velocity. If the migration velocity is lower than the velocity of the medium, events curve upward, whereas if the migration velocity is higher than the velocity of the medium, events curve downward.

The migration of common midpoint slant stacks

Geophysics ◽  
1984 ◽  
Vol 49 (3) ◽  
pp. 237-249 ◽  
Author(s):  
Richard Ottolini ◽  
Jon F. Claerbout
Keyword(s):  
Seismic Data ◽  
Lateral Velocity ◽  
Velocity Inversion ◽  
Interval Velocity ◽  
Reflection Seismic ◽  
Straightforward Method ◽  
Velocity Surface ◽  
Velocity Variations ◽  
Migration Equation ◽  
Ray Parameter

Reflection seismic data can be imaged by migrating common midpoint slant stacks. The basic method is to assemble slant stack sections from the slant stack of each common midpoint gather at the same ray parameter. Earlier investigators have described migration methods for slant stacked shot profiles or common receiver gathers instead of common midpoint gathers. However, common midpoint slant stacks enjoy the practical advantages of midpoint coordinates. In addition, the migration equation makes no approximation for steep dips, wide offsets, or vertical velocity variations. A theoretical disadvantage is that there is no exact treatment of lateral velocity variations. Slant stack migration is a method of “migration before stack.” It solves the dip selectivity problem of conventional stacking, particularly when horizontal reflectors intersect steep dipping reflectors. The correct handling of all dips also improves lateral resolution in the image. Slant stack migration provides a straightforward method of measuring interval velocity after migration has improved the seismic data. The kinematics (traveltime treatment) of slant stack migration is also accurate for postcritical reflections and refractions. These events transform into a p-τ surface with the additional dimension of midpoint. The slant stack migration equation converts the p-τ surface into a depth‐midpoint velocity surface. As with migration in general, the effects of dip are automatically accounted for during velocity inversion.

Sign in / Sign up

Export Citation Format

Share Document