A two‐pass approximation to 3-D prestack migration

Geophysics ◽  
1996 ◽  
Vol 61 (2) ◽  
pp. 409-421 ◽  
Author(s):  
Anat Canning ◽  
Gerald H. F. Gardner
Keyword(s):  
Computation Time ◽  
Three Dimensions ◽  
Velocity Analysis ◽  
Depth Migration ◽  
Prestack Migration ◽  
Time Migration ◽  
Data Volume ◽  
Land Survey ◽  
Zero Offset ◽  
The One

A two‐pass approximation to 3-D Kirchhoff migration simplifies the migration procedure by reducing it to a succession of 2-D operations. This approach has proven very successful in the zero‐offset case. A two‐pass approximation to 3-D migration is described here for the prestack case. Compared to the one‐pass approach, the scheme presented here provides significant reduction in computation time and a relatively simple data manipulation scheme. The two‐pass method was designed using velocity independent prestack time migration (DMO‐PSI) applied in the crossline direction, followed by conventional prestack depth migration in the inline direction. Velocity analysis, an important part of prestack migration, is also included in the two‐pass scheme. It is carried out as a 2-D procedure after 3-D effects are removed from the data volume. The procedure presented here is a practical full volume 3-D prestack migration. One of its main benefits is a realistic and efficient iterative velocity analysis procedure in three dimensions. The algorithm was designed in the frequency domain and the computational scheme was optimized by processing individual frequency slices independently. Irregular trace distribution, a feature that characterizes most 3-D seismic surveys, is implicitly accounted for within the two‐pass algorithm. A numerical example tests the performance of the two‐pass 3-D prestack migration program in the presence of a vertical velocity gradient. A 3-D land survey from a fold and thrust belt region was used to demonstrate the algorithm in a complex geological setting. The results were compared with images from other 2-D and 3-D migration schemes and show improved resolution and higher signal content.

scholarly journals Recursive prestack depth migration using CFP gathers

Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. S273-S283 ◽  
Author(s):  
Jan Thorbecke ◽  
A. J. Berkhout
Keyword(s):  
Seismic Data ◽  
Signal To Noise Ratio ◽  
Velocity Analysis ◽  
Signal To Noise ◽  
Depth Migration ◽  
Prestack Migration ◽  
Data Volume ◽  
The Common ◽  
Complex Geology ◽  
Focus Point

The common-focus-point technology (CFP) describes prestack migration by focusing in two steps: emission and detection. The output of the first focusing step represents a CFP gather. This gather defines a shot record that represents the subsurface response resulting from a focused source wavefield. We propose applying the recursive shot-record, depth-migration algorithm to the CFP gathers of a seismic data volume and refer to this process as CFP-gather migration. In the situation of complex geology and/or low signal-to-noise ratio, CFP-based image gathers are easier to interpret for nonalignment than the conventional image gathers. This makes the CFP-based image gathers better suited for velocity analysis. This important property is illustrated by examples on the Marmousi model.

Prestack depth migration of an Alberta Foothills data set—The Husky experience

Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 392-398 ◽  
Author(s):  
W.-J. Wu ◽  
L. Lines ◽  
A. Burton ◽  
H.-X. Lu ◽  
J. Zhu ◽  
...  
Keyword(s):  
Velocity Analysis ◽  
Reverse Time ◽  
Prestack Depth Migration ◽  
Data Set ◽  
Depth Migration ◽  
Geological Interpretation ◽  
Prestack Migration ◽  
Velocity Models ◽  
Depth Images ◽  
Migration Velocity Analysis

We produce depth images for an Alberta Foothills line by iteratively using a number of migration and velocity analysis techniques. In imaging steeply dipping layers of a foothills data set, it is apparent that thrust belt geology can violate the conventional assumptions of elevation datum corrections and common midpoint (CMP) stacking. To circumvent these problems, we use migration from topography in which we perform prestack depth migration on the data using correct source and receiver elevations. Migration from topography produces enhanced images of steep shallow reflectors when compared to conventional processing. In addition to migration from topography, we couple prestack depth migration with the continuous adjustment of velocity depth models. A number of criteria are used in doing this. These criteria require that our velocity estimates produce a focused image and that migrated depths in common image gathers be independent of source‐receiver offset. Velocity models are estimated by a series of iterative and interpretive steps involving prestack migration velocity analysis and structural interpretation. Overlays of velocity models on depth migrations should generally show consistency between velocity boundaries and reflection depths. Our preferred seismic depth section has been produced by using prestack reverse‐time depth migration coupled with careful geological interpretation.

Reply by the authors to Arthur E. Barnes

Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1944-1946
Author(s):  
M. Tygel ◽  
J. Schleicher ◽  
P. Hubral
Keyword(s):  
Destructive Interference ◽  
Depth Migration ◽  
Time Migration ◽  
Pulse Distortion ◽  
Nmo Correction ◽  
Important Processing ◽  
Offset Time ◽  
Zero Offset ◽  
Seismic Reflections ◽  
Processing Steps

We highly appreciate the useful remarks of Dr. Barnes relating our work to well‐known practical seismic processing effects. This is of particular interest as normal‐moveout (NMO) correction and post‐stack time migration are still two very important processing steps. Most exploration geophysicists know about the significance of pulse distortions known as “NM0 stretch” and “frequency shifting due to zero‐offset time migration.” As a result of the discussion of Dr. Barnes, it should now be possible to better appreciate the importance of our very general formulas (27) describing the pulse distortion of seismic reflections from an arbitrarily curved subsurface reflector when subjected to a prestack depth migration in 3‐D laterally inhomogeneous media. This discussion thus relates in particular to such important questions as how to correctly sample signals in the time or depth domain in order to avoid spatial aliasing, or how to stack seismic data without loss of information due to destructive interference of wavelets of different lengths.

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.

Wave-equation migration velocity analysis by focusing diffractions and reflections

Geophysics ◽  
2005 ◽  
Vol 70 (3) ◽  
pp. U19-U27 ◽  
Author(s):  
Paul C. Sava ◽  
Biondo Biondi ◽  
John Etgen
Keyword(s):  
Wave Equation ◽  
Migration Velocity ◽  
Velocity Estimation ◽  
Velocity Analysis ◽  
Depth Migration ◽  
Interval Velocity ◽  
Migration Velocity Analysis ◽  
Inversion Procedure ◽  
Velocity Information ◽  
Zero Offset

We propose a method for estimating interval velocity using the kinematic information in defocused diffractions and reflections. We extract velocity information from defocused migrated events by analyzing their residual focusing in physical space (depth and midpoint) using prestack residual migration. The results of this residual-focusing analysis are fed to a linearized inversion procedure that produces interval velocity updates. Our inversion procedure uses a wavefield-continuation operator linking perturbations of interval velocities to perturbations of migrated images, based on the principles of wave-equation migration velocity analysis introduced in recent years. We measure the accuracy of the migration velocity using a diffraction-focusing criterion instead of the criterion of flatness of migrated common-image gathers that is commonly used in migration velocity analysis. This new criterion enables us to extract velocity information from events that would be challenging to use with conventional velocity analysis methods; thus, our method is a powerful complement to those conventional techniques. We demonstrate the effectiveness of the proposed methodology using two examples. In the first example, we estimate interval velocity above a rugose salt top interface by using only the information contained in defocused diffracted and reflected events present in zero-offset data. By comparing the results of full prestack depth migration before and after the velocity updating, we confirm that our analysis of the diffracted events improves the velocity model. In the second example, we estimate the migration velocity function for a 2D, zero-offset, ground-penetrating radar data set. Depth migration after the velocity estimation improves the continuity of reflectors while focusing the diffracted energy.

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.

Fourier prestack migration by equivalent wavenumber

Geophysics ◽  
1999 ◽  
Vol 64 (1) ◽  
pp. 197-207 ◽  
Author(s):  
Gary F. Margrave ◽  
John C. Bancroft ◽  
Hugh D. Geiger
Keyword(s):  
Fourier Transform ◽  
Inverse Fourier Transform ◽  
Square Root ◽  
Constant Frequency ◽  
Prestack Migration ◽  
Change Of Variables ◽  
Time Migration ◽  
Data Volume ◽  
And Migration ◽  
Amplitude Preservation

Fourier prestack migration is reformulated through a change of variables, from offset wavenumber to a new equivalent wavenumber, which makes the migration phase shift independent of horizontal wavenumber. After the change of variables, the inverse Fourier transform over horizontal wavenumber can be performed to create unmigrated, but fully horizontally positioned, gathers at each output location. A complete prestack migration then results by imaging each gather independently with a poststack migration algorithm. This equivalent wavenumber migration (EWM) is the Fourier analog of the space‐time domain method of equivalent offset migration (EOM). The latter is a Kirchhoff time‐migration technique which forms common scatterpoint (CSP) gathers for each migrated trace and then images those gathers with a Kirchhoff summation. These CSP gathers are formed by trace mappings at constant time, and migration velocity analysis is easily done after the gathers are formed. Both EWM and EOM are motivated by the algebraic combination of a double square root equation into a single square root. This result defines equivalent wavenumber or offset. EWM is shown to be an exact reformulation of prestack f-k migration. The EWM theory provides explicit Fourier integrals for the formation of CSP gathers from the prestack data volume and the imaging of those gathers to form the final prestack migrated result. The CSP gathers are given by a Fourier mapping, at constant frequency, of the unmigrated spectrum followed by an inverse Fourier transform. The mapping requires angle‐dependent weighting factors for full amplitude preservation. The imaging expression (for each CSP gather) is formally identical to poststack migration with the result retained only at zero equivalent offset. Through a numerical simulation, the impulse responses of EOM and EWM are shown to be kinematically identical. Amplitude scale factors, which are exact in the constant velocity EWM theory, are implemented approximately in variable velocity EOM.

Prestack and poststack migration of crooked-line seismic reflection data: A case study from the South Portuguese Zone fold belt, southwestern Iberia

Geophysics ◽  
2007 ◽  
Vol 72 (2) ◽  
pp. B9-B18 ◽  
Author(s):  
C. Schmelzbach ◽  
C. Juhlin ◽  
R. Carbonell ◽  
J. F. Simancas
Keyword(s):  
Seismic Reflection ◽  
Fold Belt ◽  
The South ◽  
Prestack Migration ◽  
Time Migration ◽  
Straight Line ◽  
Prestack Time Migration ◽  
2D Imaging ◽  
The Common ◽  
Zero Offset

Crooked-line 2D seismic reflection survey geometries violate underlying assumptions of 2D imaging routines, affecting our ability to resolve the subsurface reliably. We compare three crooked-line imaging schemes involving prestack and poststack time migration using the 2D IBERSEIS deep seismic reflection profile running over the South Portuguese Zone thrust-and-fold belt to obtain crisp high-resolution images of the shallow crust. The crust is characterized by a complex subsurface geometry with conflicting dips of up to [Formula: see text]. In summary, the three schemes are (1) normal-moveout (NMO) corrections, dip-moveout (DMO) corrections, common-midpoint (CMP) stacking, CMP projection, and poststack time migration; (2) NMO corrections, DMO corrections, CMP projection, zero-offset time migration of the common-offset gathers, and CMP stacking; (3) CMP projection, prestack time migration in the common-offset domain, and CMP stacking. An essential element of all three schemes is a CMP projection routine, projecting the CMPs first binned along individual segments for preprocessing onto one straight line, which is parallel to the general dip direction of the subsurface structures. After CMP projection, the data satisfy the straight-line assumption of 2D imaging routines more closely. We observe that the prestack time-migration scheme yields comparable or more coherent synthetic and field-data images than the other two DMO-based schemes along the parts of the profile where the acquisition overall follows a straight line. However, the schemes involving DMO corrections are less plagued by migration artifacts than the prestack time-migration scheme along profile parts where the acquisition line is crooked. In particular, prominent migration artifacts on the prestack migrated synthetic data can be related to significant variations in source-receiver azimuths for which 2D prestack migration cannot account. Thus, the processing scheme including DMO corrections, CMP projection, and zero-offset migration of common-offset gathers offers a reliable and effective alternative to prestack migration for crooked-line 2D seismic reflection processing.

Generalized nonhyperbolic moveout approximation

Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. U9-U18 ◽  
Author(s):  
Sergey Fomel ◽  
Alexey Stovas
Keyword(s):  
Functional Form ◽  
Velocity Analysis ◽  
Numerical Examples ◽  
Time Migration ◽  
Analytical Approximations ◽  
Offset Time ◽  
Zero Offset ◽  
Two Parameters

Reflection moveout approximations are commonly used for velocity analysis, stacking, and time migration. A novel functional form for approximating the moveout of reflection traveltimes at large offsets is introduced. In comparison with the classic hyperbolic approximation, which uses only two parameters (zero-offset time and moveout velocity), this form involves five parameters that can be determined, in a known medium, from zero-offset computations and from tracing one nonzero-offset ray. It is called a generalized approximation because it reduces to some known three-parameter forms with a particular choice of coefficients. By testing the accuracy of the proposed approximation with analytical and numerical examples, the new approximation is shown to bring an improvement in accuracy of several orders of magnitude compared to known analytical approximations, which makes it as good as exact for many practical purposes.

True-amplitude prestack depth migration

Geophysics ◽  
2007 ◽  
Vol 72 (3) ◽  
pp. S155-S166 ◽  
Author(s):  
Feng Deng ◽  
George A. McMechan
Keyword(s):  
Velocity Ratio ◽  
Synthetic Data ◽  
Transmission Losses ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Prestack Migration ◽  
Time Migration ◽  
Geometric Spreading

Most current true-amplitude migrations correct only for geometric spreading. We present a new prestack depth-migration method that uses the framework of reverse-time migration to compensate for geometric spreading, intrinsic [Formula: see text] losses, and transmission losses. Geometric spreading is implicitly compensated by full two-way wave propagation. Intrinsic [Formula: see text] losses are handled by including a [Formula: see text]-dependent term in the wave equation. Transmission losses are compensated based on an estimation of angle-dependent reflectivity using a two-pass recursive reverse-time prestack migration. The image condition used is the ratio of receiver/source wavefield amplitudes. Two-dimensional tests using synthetic data for a dipping-layer model and a salt model show that loss-compensating prestack depth migration can produce reliable angle-dependent reflection coefficients at the target. The reflection coefficient curves are fitted to give least-squares estimates of the velocity ratio at the target. The main new result is a procedure for transmission compensation when extrapolating the receiver wavefield. There are still a number of limitations (e.g., we use only scalar extrapolation for illustration), but these limitations are now better defined.

Sign in / Sign up

Export Citation Format

Share Document