Effective AMO implementation in the log-stretch, frequency-wavenumber domain |

For the purpose of this discussion we define stretching of a
single-dimension space as any transformation from one space to
another that has the following property: at least an arbitrarily chosen sequence
of two consecutive, equal in length, intervals in the input space is
transformed into a sequence of two consecutive, equal in length,
intervals in the output space. Stretching an x-space to a y-space will be denoted as

where is a positive real number whose value does not matter for the purpose of this discussion. As it can be seen in Fig. 2, if we keep the same sampling rate ( ), aliasing can occur when doing the reverse transformation, from x to y. In order to avoid aliasing, we need to compute , the largest accceptable sampling rate in the y domain. This can sometimes lead to a larger number of samples in the domain, and thus to larger computational expense. This can be limited to some extent if the signal in the -space has been bandpassed, as is often the case with seismic data, with the largest frequency present in the data () smaller than the Nyquist frequency given by the sampling rate (). Thus, we can replace in our calculations with

which will result in a larger than that computed using , the sampling rate in the space.

strali
Illustration of how aliasing can occur
while stretching: if the same sampling rate is used for the -space
(lower plot) as for the -space (upper plot), serious aliasing will
occur when transforming back to -space. This will not happen if the
sampling rate in the -space is smaller than or equal to
Figure 2. |
---|

pystrali
Illustration of how aliasing can occur
while stretching: if the same sampling rate is used for the -space
(lower plot) as for the -space (upper plot), serious aliasing will
occur when transforming back to -space. This will not happen if the
sampling rate in the -space is smaller than or equal to
Figure 3. |
---|

In order to compute
, we will consider two points in the space, as seen in Fig. 2, such as

The largest sampling rate in the -space that will not result in aliasing is , the minimum possible value of . Suppose there is a value that minimizes . Then,

In particular, in the case of log-stretch, given by equation (1), if plays the role of from the equation above, then

will be minimum when is as large as possible, thus minimizing the expression under the logarithm. How large can get? Since the length of the seismic trace is limited to a value ,

because is the equivalent of from eq. (7) and Fig. 2. Thus, we get

Effective AMO implementation in the log-stretch, frequency-wavenumber domain |

2013-03-03