nybox
(a) The
standard representation of the seismic acquisition grid. (b) Viable
alternative to conventional orthogonal drawings. (c) Fourier
transform of (a) with Nyquist sampling limits included. (d) Fourier
transform of (b) with Nyquist sampling limits included.
Figure 1 |

At the limit of this argument, we contend that it is much easier to de-couple completely the origins of these axis and plot them parallel to each other. Having performed a Fourier transform across space of both the source and receiver axes, we present them in Figure 2.

Figure 2

Viewing these three distinct axes separately aids in the
interpretation of this entire argument. Unfortunately, there is an
historic tendency when analyzing the acquisition grid coordinates to
include midpoint-offset, *mh*, axes as diagonal axes to those shown in
Figure 1. We will avoid the use of midpoint while
casting this presentation largely in the terms of shot-profile
migration, as well as explain later the development of our *x*-axis
during the imaging condition. Further, when they are superposed, an
incorrect stretching is implied. We will briefly consider the mid-point
axis, in order to highlight the danger of this practice.

The mapping transformation of energy from one coordinate frame to the other has been historically defined as:

(1) | ||

(2) |

(3) |

The cross-correlation imaging condition with subsurface offset Rickett and Sava (2001)

(4) |

The Appendix develops in rigorous detail the wavenumber limits
acceptable in the image to eliminate completely alias contamination. The
analysis of the problem centers around the effects of the
migration process on the data *grid*, without needing to consider the
values of the data on each grid node. We thus draw
an analogy to the body of work available from crystallography,
where structure can be analyzed mathematically without
need to know what atom resides at any particular location.
Thankfully, the regular Cartesian grid on which we normally acquire
and process seismic data is a simple rectilinear crystal,
though of several more dimensions than seen under a microscope.

The reference crystal we will consider will be the archetypal seismic grid where sources and receivers occupy all locations and share the same spacing increment. The suspicious or simply inquisitive reader can now turn to the appendix to work through the details of the following result. The maximum allowable wavenumbers, , to avoid artifacts due to migration operator aliasing is

(5) |

We consider three approaches to remove the aliasing
problems associated with the acquisition and subsampling situations
mentioned above during shot-profile migration. First, wavenumbers
from the source and receiver wavefields at each depth level are
band-limited to prevent the entry of aliased duplications
into the image during the imaging condition. This does not require
eliminating these components from the propagating wavefields, as we
can save appropriate portions of the wavefields in temporary buffers
for imaging condition evaluation. Second, a band-limited source
function, with a wavenumber spectrum limited to the cutoff frequencies
imposed by the resampled shot axis, is propagated throughout the
migration process. This effectively zeros energy in the aliased band
during the convolution in the imaging condition. No additional
computational overhead is required for the latter alternative, though
anti-aliasing by band-limited imaging requires two additional Fourier
transforms for a split-step Fourier migration strategy. It should be
noted, however, that both of these approaches will remove energy across
both *k*_{x} and *k*_{h} axes.

A third alternative is to restrict the wavenumbers of the subsurface
offset axis *k*_{h} during imaging. Casting the imaging condition in
terms of it Fourier dual can allow similar mitigation options.
Because *k*_{s}-*k*_{r}=*k*_{h}, we can select (*k*_{s},*k*_{r}) combinations during
the imaging condition that do not exceed our prescribed bandwidth.
The multiplication of the source and receiver wavefields shown above
takes the form of a convolution in the Fourier domain, which can be
utilized to insert our anti-aliasing criteria. Lastly,
decimating the receiver wavefield to match the shot increment, will be
discussed in more detail with reference to shot-geophone style migration.

5/23/2004