Migration Velocity Analysis for TI Media in the Presence of Quadratic Lateral Velocity Variation

Author(s):  
M. Takanashi
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.


Author(s):  
J. Panizzardi ◽  
C. Andreoletti ◽  
N. Bienati ◽  
L. Casasanta

Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 164-175 ◽  
Author(s):  
Hans J. Tieman

The most common migration velocity analysis algorithms, iterative profile migration, focusing analysis, and stack power analysis are based on restrictive subsurface assumptions that may cause the methods to break down in the presence of lateral velocity variations. Typically, the subsurface is assumed to have either a constant velocity, or at most, a depth variable velocity. Recent innovations include the use of traveltime inversion philosophy to invert migration measurements of the curvature as a function of offset exhibited by an event following migration. Traveltime inversion makes few assumptions regarding the subsurface, but is a rather unstable process. Thus, an important question is, “Under what conditions do the traditional methods break down and make the use of traveltime inversion methods mandatory?” Marrying the common migration analysis with tomography results in a set of equations that, while useful for generating updates from migration measurements, are too complex for answering the above question. However, by restricting the subsurface to low relief structures and assuming small angle wave propagation, these updating equations can be approximated by forms that are identical to the traditional updating equations, except for a factor that is dependent upon the magnitude, position, and spatial wavelength of potential lateral velocity variations. These simpler equations indicate that even relatively long wavelength anomalies can cause updates to point in the wrong direction and iterative procedures to diverge, and under certain conditions, cause the accuracy of these updates to decrease dramatically.


Geophysics ◽  
1997 ◽  
Vol 62 (4) ◽  
pp. 1238-1249 ◽  
Author(s):  
Zhenyue Liu

Prestack depth migration provides a powerful tool for velocity analysis in complex media. Both prominent approaches to velocity analysis—depth‐focusing analysis and residual‐curvature analysis, rely on approximate formulas to update velocity. Generally, these formulas are derived under the assumptions of horizontal reflector, lateral velocity homogeneity, or small offset. Therefore, the conventional methods for updating velocity lack accuracy and computational efficiency when velocity has large, lateral variations. Here, based on ray theory, I find the analytic representation for the derivative of imaged depths with respect to migration velocity. This derivative function characterizes a general relationship between residual moveout and residual velocity. Using the derivative function and the perturbation method, I derive a new formula to update velocity from residual moveout. In the derivation, I impose no limitation on offset, dip, or velocity distribution. Consequently, I revise the residual‐curvature‐analysis method for velocity estimation in the postmigrated domain. Furthermore, my formula provides sensitivity and error estimation for migration‐based velocity analysis, which is helpful in quantifying the reliability of the estimated velocity. The theory and methodology in this paper have been tested on synthetic data (including the Marmousi data).


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 ◽  
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 ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. U73-U85 ◽  
Author(s):  
Saleh M. Al-Saleh ◽  
Jianwu Jiao

We introduce an integrated wave-equation technique for migration velocity analysis (MVA) that consists of three steps: (1) forming the extended data, (2) approximating the correct transmitted wavefield, and (3) using wavefield tomography to update the velocity model. In the first step, the crosscorrelation imaging condition is relaxed to produce other nonzero-lag common image gathers (CIG) that, combined, form a common image cube (CIC). Slicing the CIC at different crosscorrelation lags forms a series of CIGs. Flattened events will occur in the CIGs at a lag other than the zero-lag when an incorrect velocity model is used in the migration. In the second step, for each event on the CIG, we pick the focusing depth and crosscorrelation lag at which it is flattest. We then model a Green’s function by seeding a source at the focusing depth using one-way wave equation modeling, then shift the modeled wavefield with the focusing crosscorrelation lag. This process is repeated for the other primary events at different lateral and vertical positions. The result is a set of modeled data whose wavefield approximates the wavefield that would have been generated if the correct velocity model had been used to simulate these gathers. We then apply wavefield tomography on these data-driven modeled data to update the velocity model. Our inversion scheme is based on wave-equation traveltime tomography that can update the velocity model in the presence of large velocity errors and a complex environment. Tests on synthetic and real 2D seismic data confirm the method’s effectiveness in building velocity models in complex structural areas that have large lateral velocity variations.


Sign in / Sign up

Export Citation Format

Share Document