Residual migration can also be used to create an image perturbation. In its simplest form, we can build it as a difference between an improved image () and the reference image ()
(9) |
If we define , we can also write the discrete version of the image perturbation as
(10) |
(11) |
(12) |
The Equations (7) and (12) are very similar, which comes at no surprise since they effectively represent the same thing: the perturbation of the image given a perturbation of the slowness field, or equivalently, a perturbation of the (ratio) field. We will use Equation (12) to create the image perturbation, which we will then backproject in the slowness space using the adjoint of the WEMVA equation (7).
Equation (12) offers the possibility to build the image perturbation directly. We achieve this by computing three elements: the derivative of the image with respect to the depth wavenumber, and two weighting functions, one for the derivative of the depth wavenumber with respect to the velocity ratio parameter (), and the other one for the magnitude of the perturbation from the reference to the improved image.
Firstly, the image derivative in the Fourier domain, , is straightforward to compute in the space domain as
(13) |
Secondly, we can obtain the weighting representing the derivative of the depth wavenumber with respect to the velocity ratio parameter, , starting from the double square root (DSR) equation written for prestack Stolt residual migration Sava (2000):
(14) |
The derivative of with respect to is
(15) |
(16) |
Finally, can be picked from the set of residually migrated images at various values of the parameter Sava (2000). The main criterion that should be used is the flatness of the angle-domain image gathers, although in principle other derived parameters, such as stack power or semblance, can be used as well.