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

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.

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 ◽  
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 ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1200-1209 ◽  
Author(s):  
Jinming Zhu ◽  
Larry Lines ◽  
Sam Gray

Reliable seismic depth migrations require an accurate input velocity model. Inaccurate velocity estimates will distort point diffractors into smiles or frowns on a depth section. For both poststack and prestack migrated sections, high velocities cause deep smiles while low velocities cause shallow frowns on migrated gathers. However, for prestack images in the offset domain, high velocities cause deep frowns while low velocities cause shallow smiles. If the velocity is correct, there will be no variation in the depth migration as a function of offset and no smiles or frowns in the offset domain. We explain migration responses both mathematically and graphically and thereby provide the basis for depth migration velocity analysis.


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 ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1331-1339 ◽  
Author(s):  
Tariq Alkhalifah

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.


Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. U1-U8 ◽  
Author(s):  
Bingbing Sun ◽  
Tariq Alkhalifah

Macro-velocity model building is important for subsequent prestack depth migration and full-waveform inversion. Wave-equation migration velocity analysis uses the band-limited waveform to invert for velocity. Normally, inversion would be implemented by focusing the subsurface offset common-image gathers. We reexamine this concept with a different perspective: In the subsurface offset domain, using extended Born modeling, the recorded data can be considered as invariant with respect to the perturbation of the position of the virtual sources and velocity at the same time. A linear system connecting the perturbation of the position of those virtual sources and velocity is derived and solved subsequently by the conjugate gradient method. In theory, the perturbation of the position of the virtual sources is given by the Rytov approximation. Thus, compared with the Born approximation, it relaxes the dependency on amplitude and makes the proposed method more applicable for real data. We determined the effectiveness of the approach by applying the proposed method on isotropic and anisotropic vertical transverse isotropic synthetic data. A real data set example verifies the robustness of the proposed method.


Sign in / Sign up

Export Citation Format

Share Document