Theory of differential offset continuation |

Dip moveout represents a particular case of offset continuation for the output offset equal to zero. In this section, I consider the DMO case separately in order to compare the solutions of equation (1) with the Fourier-domain DMO operators, which have been the standard for DMO processing since Hale's outstanding work (Hale, 1983,1984).

Equation (64) transforms to the time-wavenumber domain
with the help of integral tables:

Setting the output offset to zero, we obtain the following DMO-like
integral operators in the - domain:

the wavenumber corresponds to the midpoint axis , and is the zeroth-order Bessel function. The Fourier transform of (83) and (84) with respect to the time variable reduces to known integrals (Gradshtein and Ryzhik, 1994) and creates explicit DMO-type operators in the frequency-wavenumber domain, as follows:

where

It is interesting to note that the first term of the continuation to zero offset (85) coincides exactly with the imaginary part of Hale's DMO operator (Hale, 1984). However, unlike Hale's, operator (82) is causal, which means that its impulse response does not continue to negative times. The non-causality of Hale's DMO and related issues are discussed in more detail by Stovas and Fomel (1996).

Though Hale's DMO is known to provide correct reconstruction of the geometry of zero-offset reflections, it does not account properly for the amplitude changes (Black et al., 1993). The preceding section of this paper shows that the additional contribution to the amplitude is contained in the second term of the OC operator (64), which transforms to the second term in the DMO operator (82). Note that this term vanishes at the input offset equal to zero, which represents the case of the inverse DMO operator.

Considering the inverse DMO operator as the continuation from zero
offset to a non-zero offset, we can obtain its representation in the
- domain from equations (77-79) as

Comparing operator (90) with Ronen's version of inverse DMO
(Ronen, 1987), one can see that if Hale's DMO is denoted
by
, then Ronen's inverse DMO is
, while the amplitude-preserving inverse (90)
is
. Here is the derivative
operator
, and
stands for the adjoint operator defined by the dot-product test

In high-frequency asymptotics, the difference between the amplitudes of the two inverses is simply the Jacobian term , asymptotically equal to . This difference corresponds exactly to the difference between Black's definition of amplitude preservation (Black et al., 1993) and the definition used in Born DMO (Bleistein, 1990; Liner, 1991), as discussed above. While operator (90) preserves amplitudes in the Born DMO sense, Ronen's inverse satisfies Black's amplitude preservation criteria. This means Ronen's operator implies that the ``geometric spreading'' correction (multiplication by time) has been performed on the data prior to DMO.

To construct a one-term DMO operator, thus avoiding the estimation of the offset derivative in (72), let us consider the problem of inverting the inverse DMO operator (90). One of the possible approaches to this problem is the least-squares iterative inversion, as proposed by Ronen (1987). This requires constructing the adjoint operator, which is Hale's DMO (or its analog) in the case of Ronen's method. The iterative least-squares approach can account for irregularities in the data geometry (Ronen, 1994; Ronen et al., 1991) and boundary effects, but it is computationally expensive because of the multiple application of the operators. An alternative approach is the asymptotic inversion, which can be viewed as a special case of preconditioning the adjoint operator (Chemingui and Biondi, 1996; Liner and Cohen, 1988). The goal of the asymptotic inverse is to reconstruct the geometry and the amplitudes of the reflection events in the high-frequency asymptotic limit.

According to Beylkin's theory of asymptotic inversion, also known as
the *generalized Radon transform* (Beylkin, 1985), two operators of the
form

constitute a pair of asymptotically inverse operators ( matching in the high-frequency asymptotics) if

where is the ``Beylkin determinant''

With respect to the high-frequency asymptotic representation, we can
recast (90) in the equivalent form by moving the time
derivative under the integral sign:

where

The amplitude factor (100) corresponds exactly to that of
Born DMO (Bleistein, 1990) in full accordance with the conclusions of the
asymptotic analysis of the offset-continuation amplitudes. An
analogous result can be obtained with the different definition of
amplitude preservation proposed by Black et al. (1993). In
the time-and-space domain, the operator asymptotically analogous to
(99) is found by applying either the stationary phase
technique (Liner, 1990; Black et al., 1993) or Goldin's method of
discontinuities (Goldin, 1990,1988), which is the
time-and-space analog of Beylkin's asymptotic inverse theory
(Stovas and Fomel, 1996). The time-and-space asymptotic DMO operator takes the
form

Theory of differential offset continuation |

2014-03-26