Tau migration and velocity analysis: Theory and synthetic examples

Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1331-1339 ◽  
Author(s):  
Tariq Alkhalifah
Keyword(s):  
Time Domain ◽  
Synthetic Data ◽  
Velocity Model ◽  
Migration Velocity ◽  
Surface Velocity ◽  
Velocity Analysis ◽  
Near Surface ◽  
Prestack Migration ◽  
Interval Velocity ◽  
Migration Velocity Analysis

Prestack migration velocity analysis in the time domain reduces the velocity‐depth ambiguity usually hampering the performance of prestack depth‐migration velocity analysis. In prestack τ migration velocity analysis, we keep the interval velocity model and the output images in vertical time. This allows us to avoid placing reflectors at erroneous depths during the velocity analysis process and, thus, avoid slowing down its convergence to the true velocity model. Using a 1D velocity update scheme, the prestack τ migration velocity analysis performed well on synthetic data from a model with a complex near‐surface velocity. Accurate velocity information and images were obtained using this time‐domain method. Problems occurred only in resolving a thin layer where the low resolution and fold of the synthetic data made it practically impossible to estimate velocity accurately in this layer. This 1D approach also provided us reasonable results for synthetic data from the Marmousi model. Despite the complexity of this model, the τ domain implementation of the prestack migration velocity analysis converged to a generally reasonable result, which includes properly imaging the elusive top‐of‐the‐reservoir layer.

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.

Seismic imaging and velocity analysis for an Alberta Foothills seismic survey

Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 721-732 ◽  
Author(s):  
Lanlan Yan ◽  
Larry R. Lines
Keyword(s):  
Seismic Imaging ◽  
Migration Velocity ◽  
Surface Velocity ◽  
Seismic Survey ◽  
Velocity Analysis ◽  
Prestack Depth Migration ◽  
Near Surface ◽  
Depth Migration ◽  
Migration Velocity Analysis ◽  
Imaging Results

Seismic imaging of complex structures from the western Canadian Foothills can be achieved by applying the closely coupled processes of velocity analysis and depth migration. For the purposes of defining these structures in the Shaw Basing area of western Alberta, we performed a series of tests on both synthetic and real data to find optimum imaging procedures for handling large topographic relief, near‐surface velocity variations, and the complex structural geology of steeply dipping formations. To better understand the seismic processing problems, we constructed a typical foothills geological model that included thrust faults and duplex structures, computed the model responses, and then compared the performance of different migration algorithms, including the explicit finite difference (f-x) and Kirchhoff integral methods. When the correct velocity was used in the migration tests, the f-x method was the most effective in migration from topography. In cases where the velocity model was not assumed known, we determined a macrovelocity model by performing migration/velocity analysis by using smiles and frowns in common image gathers and by using depth‐focusing analysis. In applying depth imaging to the seismic survey from the Shaw Basing area, we found that imaging problems were caused partly by near‐surface velocity problems, which were not anticipated in the modeling study. Several comparisons of different migration approaches for these data indicated that prestack depth migration from topography provided the best imaging results when near‐surface velocity information was incorporated. Through iterative and interpretive migration/velocity analysis, we built a macrovelocity model for the final prestack depth migration.

The finite-frequency sensitivity kernel for migration residual moveout and its applications in migration velocity analysis

Geophysics ◽  
2008 ◽  
Vol 73 (6) ◽  
pp. S241-S249 ◽  
Author(s):  
Xiao-Bi Xie ◽  
Hui Yang
Keyword(s):  
Wave Equation ◽  
Synthetic Data ◽  
Velocity Model ◽  
Migration Velocity ◽  
Data Sets ◽  
Velocity Analysis ◽  
Depth Migration ◽  
Sensitivity Kernel ◽  
Migration Velocity Analysis ◽  
Tomography Method

We have derived a broadband sensitivity kernel that relates the residual moveout (RMO) in prestack depth migration (PSDM) to velocity perturbations in the migration-velocity model. We have compared the kernel with the RMO directly measured from the migration image. The consistency between the sensitivity kernel and the measured sensitivity map validates the theory and the numerical implementation. Based on this broadband sensitivity kernel, we propose a new tomography method for migration-velocity analysis and updating — specifically, for the shot-record PSDM and shot-index common-image gather. As a result, time-consuming angle-domain analysis is not required. We use a fast one-way propagator and multiple forward scattering and single backscattering approximations to calculate the sensitivity kernel. Using synthetic data sets, we can successfully invert velocity perturbations from the migration RMO. This wave-equation-based method naturally incorporates the wave phenomena and is best teamed with the wave-equation migration method for velocity analysis. In addition, the new method maintains the simplicity of the ray-based velocity analysis method, with the more accurate sensitivity kernels replacing the rays.

CRP-based seismic migration velocity analysis

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. U21-U28 ◽  
Author(s):  
Weihong Fei ◽  
George A. McMechan
Keyword(s):  
Velocity Model ◽  
Migration Velocity ◽  
Data Sets ◽  
Velocity Analysis ◽  
Lateral Velocity ◽  
Depth Migration ◽  
Interval Velocity ◽  
Model Parameterization ◽  
Migration Velocity Analysis ◽  
Local Interval

A new migration velocity analysis is developed by combining the speed of parsimonious prestack depth migration with velocity adjustments estimated within and across common-reflection-point (CRP) gathers. The proposed approach is much more efficient than conventional tomographic velocity analysis because only the traces that contribute to a series of CRP gathers are depth migrated at each iteration. The local interval-velocity adjustments for each CRP are obtained by maximizing the stack amplitude over the predicted (nonhyperbolic) moveout in each CRP gather; this does not involve retracing rays. At every iteration, the velocity in each pixel is updated by averaging over all the predicted velocity updates. Finally, CRP positions and orientations are updated by parsimonious migration, and rays are retraced to define new CRP gathers for the next iteration; this ensures internal consistency between the updated velocity model and the CRP gather. Because the algorithm has a gridded-model parameterization, no explicit representation or fitting of reflectors is involved. Strong lateral-velocity variations, such as those found at salt flanks, can be handled. Application to synthetic and field data sets show that the proposed algorithm works effectively and efficiently.

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.

Time-migration velocity analysis by image-wave propagation of common-image gathers

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE161-VE171 ◽  
Author(s):  
J. Schleicher ◽  
J. C. Costa ◽  
A. Novais
Keyword(s):  
Wave Propagation ◽  
Seismic Event ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Reference Velocity ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Velocity Image ◽  
The Common

Image-wave propagation or velocity continuation describes the variation of the migrated position of a seismic event as a function of migration velocity. Image-wave propagation in the common-image gather (CIG) domain can be combined with residual-moveout analysis for iterative migration velocity analysis (MVA). Velocity continuation of CIGs leads to a detection of those velocities in which events flatten. Although image-wave continuation is based on the assumption of a constant migration velocity, the procedure can be applied in inhomogeneous media. For this purpose, CIGs obtained by migration with an inhomogeneous macrovelocity model are continued starting from a constant reference velocity. The interpretation of continued CIGs, as if they were obtained from residual migrations, leads to a correction formula that translates residual flattening velocities into absolute time-migration velocities. In this way, the migration velocity model can be improved iteratively until a satisfactory result is reached. With a numerical example, we found that MVA with iterative image continuation applied exclusively to selected CIGs can construct a reasonable migration velocity model from scratch, without the need to build an initial model from a previous conventional normal-moveout/dip-moveout velocity analysis.

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.

Sign in / Sign up

Export Citation Format

Share Document