Wave-equation migration velocity analysis by residual moveout fitting |

In this appendix I present the analytical development needed to compute all the terms in equation 16. Equations 14 and 15 provide the expression for computing the derivatives of the moveout parameters with respect to slowness as:

where the elements of the matrix are computed using either equation 9 or equation 11, and the elements of the matrix are given by

In this last expression . Given the moveout parametrization expressed in 7, and the previous expression simplifies into the following:

Furthermore, when , equation A-2 further simplifies into:

The derivative of the image vector with respect to slowness, are evaluated by applying the conventional wave-equation tomography operator that links perturbations in the slowness model to perturbations in the propagated wavefields by a first-order Born linearization of the wave equation.

Applying the chain rule to equation 1,
and taking into account the offset-to-angle transformation 2,
we can write

where the wavefields and are computed with the background slowness, and the wavefield derivatives and are computed by the conventional adjoint-state methodology that is at the basis of full waveform inversion and wave-equation tomography.

In more compact matrix notation the previous expression can be written as

where the matrices and are composed of the wavefields for every source and depth level, and properly shifted in space by the subsurface offset. For the computation of the gradient, we need to apply the adjoint operator that is:

Almomin and Tang (2010) present an equivalent, but different, derivation of an algorithm to compute the application of the operator , (or its adjoint) to a vector of slowness perturbations (or image perturbations).

Wave-equation migration velocity analysis by residual moveout fitting |

2010-11-26