# Integral DMO correction

Series | Investigations in Geophysics |
---|---|

Author | Öz Yilmaz |

DOI | http://dx.doi.org/10.1190/1.9781560801580 |

ISBN | ISBN 978-1-56080-094-1 |

Store | SEG Online Store |

In migration principles, we reviewed the migration process based on Kirchhoff summation. Dip-moveout correction also can be formulated as an integration process ^{[1]}^{[2]}^{[3]}. The integral DMO correction is particularly the preferred method for data with irregular spatial sampling and 3-D data with large variations in source-receiver azimuths (processing of 3-D seismic data).

The traveltime trajectory associated with the dip-moveout correction operator is given by the following equation (Section E.3):

**(**)

Equation (**21**) describes an ellipse with the following properties (Figure 5.1-12):

- Semi-major axis in midpoint
*y*_{0}direction:*a*=*h*. - Semi-minor axis in time
*τ*_{0}direction:*b*=*t*._{n}

The ellipse of equation (**21**) in the *y*_{0} − *τ*_{0} plane describes the *impulse response* of a dip-moveout operator applied to nonzero-offset data with offset 2*h*. In Figure 5.1-12, the coordinate of midpoint *M* is *y _{n}* and

*y*

_{0}coordinate is referenced to

*y*. Note also that the maximum lateral extent of the ellipse — the aperture of the DMO operator, is equal to offset 2

_{n}*h*. Figures 5.1-6b and 5.1-11 show the DMO ellipses associated with the impulse responses of the frequency-wavenumber and log-stretch DMO operators.

Analogous to the semicircle superposition technique for migration (migration principles), DMO correction can be viewed as mapping of an amplitude *A*_{0} at time *t _{n}* on a normal-moveout-corrected trace at midpoint

*y*to an amplitude

_{n}*A*

_{1}at time

*τ*

_{0}on a trace at midpoint

*y*

_{0}. The vertical excursion Δ

*t*and the horizontal excursion Δ

_{DMO}*y*denoted in Figure 5.1-12 are given by equations (

_{DMO}**9a**,

**9b**). If trace spacing in the midpoint direction is Δ

*y*, then the lateral excursion is Δ

*y*/Δ

_{DMO}*y*traces, the maximum excursion being

*h*/Δ

*y*traces.

**(**)

**(**)

While the kinematics of the DMO correction operator is given by equation (**21**), the question remains as to the amplitude and the phase of the operator. Although the method was first described by Deregowski and Rocca ^{[1]}, a formal derivation of integral DMO correction with amplitude preserving characteristics is given by Black ^{[3]}. Mathematical treatise of the problem is quite involved, and we shall only refer to the results of the analysis.

Rewrite equation (**21**) explicitly in terms of the normal-moveout-corrected time variable *t _{n}* as

**(**)

where

**(**)

and *y*_{0} ≤ *h*. Given the output sample time *τ*_{0} on the DMO ellipse, equation (**22**) gives the input sample time *t _{n}* (Figure 5.1-12). The output sample value

*P*(

_{out}*y*

_{0},

*τ*

_{0};

*h*) is computed by summing over the input sample values

*P*(

_{in}*y*,

_{n}*t*;

_{n}*h*) over the DMO operator aperture

**(**)

where Δ*y* is the trace spacing in midpoint direction.

**Figure 5.1-6**Intermediate results from DMO processing the nonzero-offset synthetic data derived from the depth model in Figure 5.1-3: (a) common-offset sections with offset range from 50 to 1550 m and an increment of 300 m sorted from the NMO-corrected gathers as in Figure 5.1-4c; (b) impulse responses of the DMO operators applied to the common-offset gathers; (c) common-offset sections as in (a) after DMO correction; (d) CMP gathers sorted from the common-offset sections as in (c) at midpoint locations from 32 to 63 as denoted in Figure 5.1-3 with an increment of 3.

Equation (**23**) is adapted from Black ^{[3]} and is similar to equation (4-5) which describes the Kirchhoff summation. For the 2-D application of DMO correction, the *ρ*(*t _{n}*) filter has an amplitude spectrum of the form with

*ω*being the temporal frequency associated with the input time variable

_{n}*t*, and a phase spectrum equal to

_{n}*π*/4.

In the integral implementation of DMO correction by Deregowski and Rocca ^{[1]}, and Deregowski ^{[2]}, the term *α* is set to unity in equation (**24**). Moreover, in the Liner ^{[4]} and Bleistein ^{[5]} implementation of integral DMO correction, the term 2*α*^{2} − 1 in equation (**24**) is replaced with *α*^{2}(2*α*^{2} − 1). Nevertheless, within the context of a conventional processing sequence which includes geometric spreading correction prior to DMO correction, the amplitude scaling (2*α*^{2} − 1)/h described here preserves relative amplitudes. Other DMO amplitude scaling strategies include those suggested by Sorin and Ronen ^{[6]} and Gardner and Forel ^{[7]}.

In practical implementations of the integral DMO correction, a user-defined aperture commonly is imposed on the operator to avoid aliasing along the steep flanks of the DMO ellipse, especially at late times. When the DMO ellipse is truncated before it reaches its fullest lateral extent, amplitude distribution along the elliptic trajectory is adjusted accordingly, such that the amplitude is tapered to zero at the truncation point on the ellipse. Figure 5.1-13 shows the impulse responses of an integral DMO operator based on equation (**24**) for 1000-m, 2000-m, and 3000-m offsets. Compare with the impulse responses in Figure 5.1-11 and note the truncation at steep flanks of the DMO ellipses. The problem of spatial aliasing due to undersampling and the adverse effect of irregular sampling on DMO correction are particularly relevant for 3-D data; as such, these issues will be dealt with in processing of 3-D seismic data.

## References

- ↑
^{1.0}^{1.1}^{1.2}Deregowski and Rocca (1981), Deregowski, S. M. and Rocca, F., 1981, Geometrical optics and wave theory for constant-offset sections in layered media: Geophys. Prosp., 29, 374–406. - ↑
^{2.0}^{2.1}Deregowski, 1987, Deregowski, S. M., 1987, An integral implementation of dip-moveout: Geophys. Trans. of Geophys. Inst. of Hungary, 33, 11–22. - ↑
^{3.0}^{3.1}^{3.2}Black et al. (1993), Black, J., Schleicher, K. L. and Zhang, L., 1993, True-amplitude imaging and dip moveout: Geophysics, 58, 47–66. - ↑ Liner, 1990, Liner, C. L., 1990, General theory and comparative anatomy of dip-moveout: Geophysics, 55, 595–607.
- ↑ Bleistein (1990), Bleistein, N., 1990, Born DMO revisited: 60th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1366–1369.
- ↑ Sorin and Ronen (1989), Sorin, V. and Ronen, J., 1989, Ray-geometrical analysis of dip moveout amplitude distribution: Geophysics, 54, 1333–1335.
- ↑ Gardner and Forel (1990), Gardner, G. F. H. and Forel, D., 1990, Amplitude preservation equations for DMO: Geophysics, 55, 485–487.

## See also

- Prestack partial migration
- Frequency-wavenumber DMO correction
- Log-stretch DMO correction
- Velocity errors
- Variable velocity
- Turning-wave migration