Seismic constant‐velocity remigration

Geophysics ◽  
1997 ◽  
Vol 62 (2) ◽  
pp. 589-597 ◽  
Author(s):  
Jörg Schleicher ◽  
Peter Hubral ◽  
German Höcht ◽  
Frank Liptow
Keyword(s):  
Constant Velocity ◽  
Wave Equations ◽  
Single Point ◽  
Velocity Model ◽  
Migration Velocity ◽  
Body Waves ◽  
Reflection Point ◽  
Ray Trajectories ◽  
Zero Offset ◽  
Migration Velocities

When a seismic common midpoint (CMP) stack or zero‐offset (ZO) section is depth or time migrated with different (constant) migration velocities, different reflector images of the subsurface are obtained. If the migration velocity is changed continuously, the (kinematically) migrated image of a single point on the reflector, constructed for one particular seismic ZO reflection signal, moves along a circle at depth, which we call the Thales circle. It degenerates to a vertical line for a nondipping event. For all other dips, the dislocation as a function of migration velocity depends on the reflector dip. In particular for reflectors with dips larger than 45°, the reflection point moves upward for increasing velocity. The corresponding curves in a Time‐migrated section are parabolas. These formulas will provide the seismic interpreter with a better understanding of where a reflector image might move when the velocity model is changed. Moreover, in that case, the reflector image as a whole behaves to some extent like an ensemble of body waves, which we therefore call remigration image waves. In the same way as physical waves propagate as a function of time, these image waves propagate as a function of migration velocity. Different migrated images can thus be considered as snapshots of image waves at different instants of migration velocity. By some simple plane‐wave considerations, image‐wave equations can be derived that describe the propagation of image waves as a function of the migration velocity. The Thales circles and parabolas then turn out to be the characteristics or ray trajectories for these image‐wave equations.

Transformation to zero offset for mode‐converted waves

Geophysics ◽  
1992 ◽  
Vol 57 (3) ◽  
pp. 474-477 ◽  
Author(s):  
Mohammed Alfaraj ◽  
Ken Larner
Keyword(s):  
Seismic Data ◽  
Constant Velocity ◽  
Reflection Point ◽  
P Waves ◽  
S Waves ◽  
Nmo Correction ◽  
Converted Waves ◽  
Zero Offset ◽  
Correct Data

The transformation to zero offset (TZO) of prestack seismic data for a constant‐velocity medium is well understood and is readily implemented when dealing with either P‐waves or S‐waves. TZO is achieved by inserting a dip moveout (DMO) process to correct data for the influence of dip, either before or after normal moveout (NMO) correction (Hale, 1984; Forel and Gardner, 1988). The TZO process transforms prestack seismic data in such a way that common‐midpoint (CMP) gathers are closer to being common reflection point gathers after the transformation.

Poststack migration of P-SV seismic data

Geophysics ◽  
1993 ◽  
Vol 58 (8) ◽  
pp. 1127-1135 ◽  
Author(s):  
Mark P. Harrison ◽  
Robert R. Stewart
Keyword(s):  
Phase Shift ◽  
Inhomogeneous Medium ◽  
Seismic Data ◽  
Migration Velocity ◽  
Velocity Function ◽  
Zero Offset ◽  
Migration Velocities

The exploding‐reflector model is only satisfied for zero‐offset P-SV data when [Formula: see text] is depth‐invariant. P-SV diffractions in a vertically‐inhomogeneous medium are approximately hyperbolic, and an expression for their migration velocity is derived. The resulting migration velocities are 6–11 percent less than the corresponding P-SV rms velocities. Migration of DMO‐corrected synthetic P-SV stacked data, using a conventional phase‐shift algorithm and the derived migration velocities, is found to adequately collapse diffractions, whereas migration using the rms velocity function gives significant overmigration.

scholarly journals Shot-gather time migration of planar reflectors without velocity model

Geophysics ◽  
2011 ◽  
Vol 76 (2) ◽  
pp. S93-S101 ◽  
Author(s):  
Andrej Bóna
Keyword(s):  
Synthetic Data ◽  
Velocity Model ◽  
Migration Velocity ◽  
Time Migration ◽  
Second Derivatives ◽  
Migration Method ◽  
Prestack Time Migration ◽  
Zero Offset ◽  
2D And 3D ◽  
Homogeneous Media

Standard migration techniques require a velocity model. A new and fast prestack time migration method is presented that does not require a velocity model as an input. The only input is a shot gather, unlike other velocity-independent migrations that also require input of data in other gathers. The output of the presented migration is a time-migrated image and the migration velocity model. The method uses the first and second derivatives of the traveltimes with respect to the location of the receiver. These attributes are estimated by computing the gradient of the amplitude in a shot gather. The assumptions of the approach are a laterally slowly changing velocity and reflectors with small curvatures; the dip of the reflector can be arbitrary. The migration velocity corresponds to the root mean square (rms) velocity for laterally homogeneous media for near offsets. The migration expressions for 2D and 3D cases are derived from a simple geometrical construction considering the image of the source. The strengths and weaknesses of the methods are demonstrated on synthetic data. At last, the applicability of the method is discussed by interpreting the migration velocity in terms of the Taylor expansion of the traveltime around the zero offset.

Migration velocity analysis with the Kirchhoff integral

Geophysics ◽  
1991 ◽  
Vol 56 (3) ◽  
pp. 365-370 ◽  
Author(s):  
Y. C. Kim ◽  
R. Gonzalez
Keyword(s):  
Maximum Amplitude ◽  
Migration Velocity ◽  
Downward Continuation ◽  
Velocity Analysis ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Prestack Time Migration ◽  
Zero Offset ◽  
Kirchhoff Integral ◽  
Migration Velocities

To obtain accurate migration velocities, we must estimate the velocity at migrated depth points. Wavefront focusing analysis with downward continuation yields the rms velocity at migrated depth points; however, the large amount of computation required for downward continuation limits use of this approach for routine processing. The purpose of this paper is to present an implementation of the Kirchhoff integral which makes the wavefront focusing analysis practical for time‐migration velocity analysis. Downward continuation focuses the wavefront to the zero offset at the depth controlled by the velocity used for the continuation. The migration velocity is then determined from the depth where the focused wavefront attains the maximum amplitude. The flexibility of the Kirchhoff integral allows us to compute only the zero‐offset trace at each depth point and lets us avoid most of the computation for the downward continuation of unstacked data. Furthermore, since the velocity is obtained from the location where the focused wavefront shows the maximum amplitude, prestack time migration with the velocity from this technique produces the maximum amplitude for the subsurface reflector.

Ground-penetrating radar data diffraction focusing without a velocity model

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. H13-H24
Author(s):  
Nikos Economou ◽  
Antonis Vafidis ◽  
Maksim Bano ◽  
Hamdan Hamdan ◽  
Jose Ortega-Ramirez
Keyword(s):  
Ground Penetrating Radar ◽  
Constant Velocity ◽  
Electromagnetic Waves ◽  
Ground Truth ◽  
Velocity Model ◽  
Migration Velocity ◽  
Dominant Frequency ◽  
Limited Information ◽  
Velocity Range ◽  
Ground Penetrating

Ground-penetrating Radar (GPR) sections commonly suffer from strong scattered energy and weak reflectors with distorted lateral continuity. This is mainly due to the gradual variation of moisture with depth, dense lateral sampling of common-offset GPR traces (which are considered as zero-offset data), along with the small wavelength of the electromagnetic waves that is comparable to the size of the shallow subsurface dielectric heterogeneities. Focusing of the diffractions requires efficient migration that, in the presence of highly heterogeneous subsurface formations, can be improved by a detailed migration velocity model. Such a velocity model is difficult to develop because the common-offset antenna array is mostly used for its reduced time and cost in the data acquisition and processing stages. In such cases, migration processes are based on limited information from velocity analysis of clear diffractions, cores, or other ground truth knowledge, often leading to insufficient imaging. We have developed a methodology to obtain GPR sections with focused diffractions that is based on multipath summation, using weighted stacking (summation) of constant-velocity migrated sections over a predefined velocity range. The success of this method depends on the assignment of an appropriate weight, for each constant-velocity migrated section to contribute to the final stack, and the optimal width of the velocity range used. Additionally, we develop a postmultipath summation processing step, which consists of time-varying spectral whitening, to deal with the decrease of the dominant frequency due to attenuation effects and the additional degraded resolution expected by the constant migration summed images. This imaging strategy leads to GPR sections with sufficiently focused diffractions, enhancing the lateral and the temporal resolution, without the need to explicitly build a migration velocity model.

Numeric implementation of wave-equation migration velocity analysis operators

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE145-VE159 ◽  
Author(s):  
Paul Sava ◽  
Ioan Vlad
Keyword(s):  
Wave Equation ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Estimation ◽  
Velocity Analysis ◽  
Velocity Models ◽  
Image Optimization ◽  
Migration Velocity Analysis ◽  
Zero Offset ◽  
Wavefield Extrapolation

Wave-equation migration velocity analysis (MVA) is a technique similar to wave-equation tomography because it is designed to update velocity models using information derived from full seismic wavefields. On the other hand, wave-equation MVA is similar to conventional, traveltime-based MVA because it derives the information used for model updates from properties of migrated images, e.g., focusing and moveout. The main motivation for using wave-equation MVA is derived from its consistency with the corresponding wave-equation migration, which makes this technique robust and capable of handling multipathing characterizing media with large and sharp velocity contrasts. The wave-equation MVA operators are constructed using linearizations of conventional wavefield extrapolation operators, assuming small perturbations relative to the background velocity model. Similar to typical wavefield extrapolation operators, the wave-equation MVA operators can be implemented in the mixed space-wavenumber domain using approximations of differentorders of accuracy. As for wave-equation migration, wave-equation MVA can be formulated in different imaging frameworks, depending on the type of data used and image optimization criteria. Examples of imaging frameworks correspond to zero-offset migration (designed for imaging based on focusing properties of the image), survey-sinking migration (designed for imaging based on moveout analysis using narrow-azimuth data), and shot-record migration (also designed for imaging based on moveout analysis, but using wide-azimuth data). The wave-equation MVA operators formulated for the various imaging frameworks are similar because they share elements derived from linearizations of the single square-root equation. Such operators represent the core of iterative velocity estimation based on diffraction focusing or semblance analysis, and their applicability in practice requires efficient and accurate implementation. This tutorial concentrates strictly on the numeric implementation of those operators and not on their use for iterative migration velocity analysis.

scholarly journals Zero-offset sections with a deblurring filter in the time domain

Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. S239-S249
Author(s):  
Shihang Feng ◽  
Oz Yilmaz ◽  
Yuqing Chen ◽  
Gerard T. Schuster
Keyword(s):  
Time Domain ◽  
Constant Velocity ◽  
Field Data ◽  
Velocity Model ◽  
Velocity Models ◽  
Time Migration ◽  
Prestack Time Migration ◽  
Zero Offset ◽  
Image Volume ◽  
The Time Domain

The conventional common-midpoint stack is not equivalent to the zero-offset section due to the existence of velocity uncertainty. To obtain a zero-offset reflection section that preserves most reflections and diffractions, we have developed a velocity-independent workflow for reconstructing a high-quality zero-offset reflection section from prestack data with a deblurring filter. This workflow constructs a migration image volume by prestack time migration using a series of constant-velocity models. A deblurring filter for each constant-velocity model is applied to each time-migration image to get a deblurred image volume. To preserve all events in the image volume, each deblurred image panel is demigrated and then summed over the velocity axis. Compared with the workflow without a deblurring filter, the composite zero-offset reflection section has higher resolution and fewer migration artifacts. We evaluate applications of our method to synthetic and field data to validate its effectiveness.

3D common-reflection-point-based seismic migration velocity analysis

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. S161-S167 ◽  
Author(s):  
Weihong Fei ◽  
George A. McMechan
Keyword(s):  
Synthetic Data ◽  
Computation Time ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Reflection Point ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Migration Velocity Analysis ◽  
Spatial Coordinates

Three-dimensional prestack depth migration and depth residual picking in common-image gathers (CIGs) are the most time-consuming parts of 3D migration velocity analysis. Most migration-based velocity analysis algorithms need spatial coordinates of reflection points and CIG depth residuals at different offsets (or angles) to provide updated velocity information. We propose a new algorithm that can analyze 3D velocity quickly and accurately. Spatial coordinates and orientations of reflection points are provided by a 3D prestack parsimonious depth migration; the migration involves only the time samples picked from the salient reflection events on one 3D common-offset volume. Ray tracing from the reflection points to the surface provides a common-reflection-point (CRP) gather for each reflection point. Predicted (nonhyperbolic) moveouts for local velocity perturbations, based on maximizing the stacked amplitude, give the estimated velocity updates for each CRP gather. Then the velocity update for each voxel in the velocity model is obtained by averaging over all predicted velocity updates for that voxel. Prior model constraints may be used to stabilize velocity updating. Compared with other migration velocity analyses, the traveltime picking is limited to only one common-offset volume (and needs to be done only once); there is no need for intensive 3D prestack depth migration. Hence, the computation time is orders of magnitude less than other migration-based velocity analyses. A 3D synthetic data test shows the algorithm works effectively and efficiently.

Migration velocity analysis from locally coherent events in 2‐D laterally heterogeneous media, Part I: Theoretical aspects

Geophysics ◽  
2002 ◽  
Vol 67 (4) ◽  
pp. 1202-1212 ◽  
Author(s):  
Hervé Chauris ◽  
Mark S. Noble ◽  
Gilles Lambaré ◽  
Pascal Podvin
Keyword(s):  
Cost Function ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Estimation ◽  
Velocity Analysis ◽  
Seismic Reflection Data ◽  
Velocity Models ◽  
Migration Velocity Analysis ◽  
Paraxial Ray ◽  
The Cost

We present a new method based on migration velocity analysis (MVA) to estimate 2‐D velocity models from seismic reflection data with no assumption on reflector geometry or the background velocity field. Classical approaches using picking on common image gathers (CIGs) must consider continuous events over the whole panel. This interpretive step may be difficult—particularly for applications on real data sets. We propose to overcome the limiting factor by considering locally coherent events. A locally coherent event can be defined whenever the imaged reflectivity locally shows lateral coherency at some location in the image cube. In the prestack depth‐migrated volume obtained for an a priori velocity model, locally coherent events are picked automatically, without interpretation, and are characterized by their positions and slopes (tangent to the event). Even a single locally coherent event has information on the unknown velocity model, carried by the value of the slope measured in the CIG. The velocity is estimated by minimizing these slopes. We first introduce the cost function and explain its physical meaning. The theoretical developments lead to two equivalent expressions of the cost function: one formulated in the depth‐migrated domain on locally coherent events in CIGs and the other in the time domain. We thus establish direct links between different methods devoted to velocity estimation: migration velocity analysis using locally coherent events and slope tomography. We finally explain how to compute the gradient of the cost function using paraxial ray tracing to update the velocity model. Our method provides smooth, inverted velocity models consistent with Kirchhoff‐type migration schemes and requires neither the introduction of interfaces nor the interpretation of continuous events. As for most automatic velocity analysis methods, careful preprocessing must be applied to remove coherent noise such as multiples.

Migration velocity analysis based on focusing diffractions generated using the CRS wavefront attributes: Application to real land data

Geophysics ◽  
2021 ◽  
pp. 1-50
Author(s):  
German Garabito ◽  
José Silas dos Santos Silva ◽  
Williams Lima
Keyword(s):  
Time Domain ◽  
Real Data ◽  
Normal Incidence ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Velocity Models ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
The Time Domain

In land seismic data processing, the prestack time migration (PSTM) image remains the standard imaging output, but a reliable migrated image of the subsurface depends on the accuracy of the migration velocity model. We have adopted two new algorithms for time-domain migration velocity analysis based on wavefield attributes of the common-reflection-surface (CRS) stack method. These attributes, extracted from multicoverage data, were successfully applied to build the velocity model in the depth domain through tomographic inversion of the normal-incidence-point (NIP) wave. However, there is no practical and reliable method for determining an accurate and geologically consistent time-migration velocity model from these CRS attributes. We introduce an interactive method to determine the migration velocity model in the time domain based on the application of NIP wave attributes and the CRS stacking operator for diffractions, to generate synthetic diffractions on the reflection events of the zero-offset (ZO) CRS stacked section. In the ZO data with diffractions, the poststack time migration (post-STM) is applied with a set of constant velocities, and the migration velocities are then selected through a focusing analysis of the simulated diffractions. We also introduce an algorithm to automatically calculate the migration velocity model from the CRS attributes picked for the main reflection events in the ZO data. We determine the precision of our diffraction focusing velocity analysis and the automatic velocity calculation algorithms using two synthetic models. We also applied them to real 2D land data with low quality and low fold to estimate the time-domain migration velocity model. The velocity models obtained through our methods were validated by applying them in the Kirchhoff PSTM of real data, in which the velocity model from the diffraction focusing analysis provided significant improvements in the quality of the migrated image compared to the legacy image and to the migrated image obtained using the automatically calculated velocity model.

Sign in / Sign up

Export Citation Format

Share Document