Fast model-based migration velocity analysis and reflector shape estimation

Geophysics ◽  
2005 ◽  
Vol 70 (2) ◽  
pp. U9-U17 ◽  
Author(s):  
Weihong Fei ◽  
George A. McMechan

Migration velocity analysis can be made more efficient by preselecting the traces that contribute to a series of common-reflection-point (CRP) gathers and migrating only those traces. The data traces that contribute to a CRP for one reflection point on one layer are defined in a two-step procedure. First, poststack parsimonious Kirchhoff depth migration of zero-offset (or stacked) traces defines approximate reflector positions and orientations. Then, ray tracing from the reflection points for nonzero reflection angles defines the source and receiver locations of the data traces in the CRP gather. These traces are then prestack depth migrated, and the interval velocity model adjustment is obtained by fitting the velocity that maximizes the stack amplitude over the predicted (nonhyperbolic) moveout. A small number (2–3) of iterations converge to a 2D model of layer shape and interval velocity. Further efficiency is obtained by implementing layer stripping. The computation time is greatly reduced by combining parsimonious migration with processing only the salient portions of the whole seismic data set. The algorithm can handle lateral velocity variation within each layer as well as constant velocity. The computation time of the proposed algorithm is of the same order as that of the standard rms velocity scan method, but it does not have the inherent assumptions of the velocity scan method and is faster than current iterative prestack depth migration velocity analysis methods for typical field data.

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. S161-S167 ◽  
Author(s):  
Weihong Fei ◽  
George A. McMechan

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.


Geophysics ◽  
2005 ◽  
Vol 70 (3) ◽  
pp. U19-U27 ◽  
Author(s):  
Paul C. Sava ◽  
Biondo Biondi ◽  
John Etgen

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.


Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. U21-U28 ◽  
Author(s):  
Weihong Fei ◽  
George A. McMechan

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.


Geophysics ◽  
2003 ◽  
Vol 68 (2) ◽  
pp. 634-640 ◽  
Author(s):  
Paul C. Sava

Prestack Stolt residual migration can be applied to seismic images that are depth migrated using wavefield extrapolation techniques. This method has significant advantages over other methods that estimate interval velocity functions for depth migration. It is more accurate than methods that are based on focusing the stack of migrated images by a residual NMO operation, so it provides a more accurate estimate of the correct migration velocities. Also, it is conceptually simpler and easier to implement than traveltime tomography methods. Although the theory is developed assuming constant velocity, the method can be used for depth‐migrated images produced with smoothly varying velocity models, since the residually migrated images depend only on the ratio of the reference and updated velocities. This method closely resembles Stolt‐stretch techniques, so it inherits the Stolt method's speed and convenience. The main applications of this method are in migration velocity analysis (MVA), where it can be used to investigate the effects of gross velocity changes on the migrated image, and as a tool for residual image enhancement used by more sophisticated MVA methods, e.g., wave‐equation migration velocity analysis.


Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. U87-U96 ◽  
Author(s):  
Mamoru Takanashi ◽  
Ilya Tsvankin

One of the most serious problems in anisotropic velocity analysis is the trade-off between anisotropy and lateral heterogeneity, especially if velocity varies on a scale smaller than the maximum offset. We have developed a P-wave MVA (migration velocity analysis) algorithm for transversely isotropic (TI) models that include layers with small-scale lateral heterogeneity. Each layer is described by constant Thomsen parameters [Formula: see text] and [Formula: see text] and the symmetry-direction velocity [Formula: see text] that varies as a quadratic function of the distance along the layer boundaries. For tilted TI media (TTI), the symmetry axis is taken orthogonal to the reflectors. We analyzed the influence of lateral heterogeneity on image gathers obtained after prestack depth migration and found that quadratic lateral velocity variation in the overburden can significantly distort the moveout of the target reflection. Consequently, medium parameters beneath the heterogeneous layer(s) are estimated with substantial error, even when borehole information (e.g., check shots or sonic logs) is available. Because residual moveout in the image gathers is highly sensitive to lateral heterogeneity in the overburden, our algorithm simultaneously inverts for the interval parameters of all layers. Synthetic tests for models with a gently dipping overburden demonstrate that if the vertical profile of the symmetry-direction velocity [Formula: see text] is known at one location, the algorithm can reconstruct the other relevant parameters of TI models. The proposed approach helps increase the robustness of anisotropic velocity model-building and enhance image quality in the presence of small-scale lateral heterogeneity in the overburden.


Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 392-398 ◽  
Author(s):  
W.-J. Wu ◽  
L. Lines ◽  
A. Burton ◽  
H.-X. Lu ◽  
J. Zhu ◽  
...  

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.


1991 ◽  
Author(s):  
M. Turhan Taner ◽  
Richard W. Postma ◽  
Lee Lu ◽  
Edip Baysal

Sign in / Sign up

Export Citation Format

Share Document