- One possibility is to flatten events in angle-domain common-image gathers (ADCIG) by using residual moveout.
- Another possibility is to use residual migration to flatten ADCIGs and, at the same time, focus diffractions which can be observed in common-offset sections.

In principle, both focusing in space (along the midpoint axis) and focusing in offset are velocity indicators, and they should be used together to achieve the highest accuracy. In this paper, however, we emphasize migration velocity analysis using only focusing of diffractions along the spatial axes.

In our current implementation, we use prestack Stolt residual migration Sava (2003); Stolt (1996) as the image enhancement operator (). This residual migration operator applied to the background image creates new images (), functions of a scalar parameter (), which represents the ratio of a new slowness model relative to the background one:

(8) |

(9) |

We address this challenge by using *linearized image perturbations* .
We run residual migration for a large number of parameters
and pick at every location the value where the image
is best focused. Then we estimate at every
point the gradient of the image relative to the parameter
and construct the image perturbations using the following relation:

(10) |

The main benefit of constructing image perturbations with equation (10) is that we avoid the danger of subtracting images that are out of phase. In fact, we do not subtract images at all, but we simply construct the image perturbation that corresponds to a particular map of residual migration parameters (). In this way, we honor the information from residual migration, but we are safe relative to the limits of the first-order Born approximation.

Figures 1 and 2 illustrate the migration velocity analysis methodology using residual migration and linearized image perturbations.

Figure 1

Figure 1 shows three simple models with diffractors and
reflectors with a constant velocity *v*=2000 m/s.
We use these three models to illustrate different situations:
an isolated diffractor at location *x*=2000 m and depth *z*=900 m,
(Figure 1, left),
the same diffractor flanked by other diffractors at *z*=1100 m
(Figure 1, middle),
and finally the same diffractor next to a short reflector at *z*=1100 m
(Figure 1, right).

We migrate each synthetic datum with an incorrect velocity, *v*=1800 m/s,
and then run residual migration with various velocity ratios
from to .From top to bottom, each row corresponds to a different velocity ratio
as follows: 0.7, 0.8, 0.9, 1.0, 1.1.
For all residual migration examples, we have eliminated the vertical
shift induced by the different velocities, such that only the
diffraction component of residual migration is left. Thus, we can better
compare focusing of various events without being distracted by their
vertical movement.

The images at are the best focused images. Since both the backgrounds and the perturbations are constant, the images focus at a single ratio parameter. The ratio difference between the original images at and the best focused images at is .In general, the images focus at different ratios at different locations; therefore is a spatially variable function.

Using the background images and the measured , we compute the linearized image perturbations (Figure 2, top), and using the WEMVA operator we compute the corresponding slowness perturbations after 15 linear iterations (Figure 2, bottom). The image perturbations closely resemble the background image (Figure 1, fourth row from top), with a phase shift and appropriate scaling with the measured .

For all models in Figures 1 and
2, we measure focusing on a single event
(the main diffractor at *x*=2000 m and depth *z*=900 m), but assign
the computed to other elements of the image in the
vicinity of this diffractor. The rationale for doing so is that we can
assume that all elements at a particular location are influenced
by roughly the same part of the model. Therefore, not only is
a priori separation of the diffractors from the reflectors not
required, but the additional elements present in the image perturbation
add robustness to the inversion.

DIFLAmva
Migration velocity analysis for the three simple
synthetic models in Figure 1.
The top row depicts image perturbations, and
the bottom row depicts slowness perturbations
obtained after 15 linear iterations.
Figure 2 |

5/23/2004