Delay in general
Series | Geophysical References Series |
---|---|
Title | Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing |
Author | Enders A. Robinson and Sven Treitel |
Chapter | 7 |
DOI | http://dx.doi.org/10.1190/1.9781560801610 |
ISBN | 9781560801481 |
Store | SEG Online Store |
The concepts of minimum delay and maximum delay have been defined for two-length wavelets. For example, take the minimum-delay wavelet (7, 3), in which the coefficient 7 occurs at time index 0 and the coefficient 3 occurs at time index 1. The reverse is the maximum-delay wavelet (3, 7), in which the coefficient 3 occurs at time index 0 and the coefficient 7 occurs at time index 1. Now let us define these concepts for wavelets of greater length. To do so, take a group of dipoles, each with the minimum-delay wavelet on the left (Table 3).
The convolution of the minimum-delay wavelets in the group (i.e., the convolution of the wavelets on the left) yields a minimum-delay wavelet
( )
We shall now define maximum delay in general. The convolution of the maximum-delay wavelets in the group of dipoles (i.e., the convolution of the wavelets on the right [Table 3]) yields a maximum-delay wavelet. Thus, (0.25, 0.5, 1, 2) is a maximum-delay wavelet. We see that the maximum-delay wavelet is the reverse of the minimum-delay wavelet (2, 1, 0.5, 0.25), as we would expect. The reverse of a finite-length minimum-delay wavelet is thus a maximum-delay wavelet.
Let us now define mixed delay in general. The convolution of one wavelet from each dipole in the group, chosen so that we have a mixture of minimum-delay and maximum-delay wavelets, gives a mixed-delay wavelet. Thus,
Dipoles | Minimum delay | Maximum delay |
---|---|---|
Dipole 1 | (2, 1) | (1, 2) |
Dipole 2 | (i, 0.5) | (0.5, –i) |
Dipole 3 | (–i, 0.5) | (0.5, i) |
( )
is a mixed-delay wavelet.
The inverse of the Z-transform of the two-term equal-delay wavelet is
( )
Thus, the inverse is the discrete form of the complex exponential for positive time. Similarly, the inverse of the Z-transform of the two-term equal-delay wavelet is
( )
Thus, the Z-transform of the causal sine (with parameter ) is
( )
where = 0,1,2,3,... is the unit step, and which is
( )
We note that the denominator is the Z-transform of the convolution of the equal-delay wavelet with the equal-delay wavelet . This convolution gives the wavelet , whose inverse is seen to have the Z-transform
( )
However, this expression is equal to Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle {1/}\left(Z{\rm \ sin\ }\theta \right) } times the Z-transform of the causal sine.
Thus, we divide the causal sine by and advance the result by one time unit to obtain the inverse of the wavelet Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \left({1\ }{,\ }-{2\ \ cos\ }\theta {,\ 1}\right) } . In other words, this inverse is the sinusoid . For example, let us take the case when the frequency is f = 12.5 Hz. Then the parameter is Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \theta =2\pi f\Delta t=2\pi \left({12.5}\right)\left({0.004}\right)=\pi {/10}) } . The one-sided sinusoid for this frequency is shown in Figure 12.
For the damping factor Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \alpha } where , let us form the wavelet
( )
Its inverse is a damped one-sided sinusoid. Figure 13 shows a damped 12.5-Hz sinusoid for .
How many wavelets can be generated from N dipoles? If every two-length wavelet in a group of N dipoles is distinct, there are possible (N + 1)-length wavelets that can be generated by taking one two-length wavelet from each dipole. Of course, if some of the two-length wavelets are not distinct, some of the (N + 1)-length wavelets will not be distinct either. In any case, all the (N + 1)-length wavelets generated in this way from one group of dipoles will have the same autocorrelation. These (N + 1)-length wavelets, so generated, are said to belong to a suite of wavelets all with the same autocorrelation.
Now let us show that the minimum-delay concept is unique. Of the possible (N + 1)-length wavelets in the suite, one is a minimum-delay wavelet, another is a maximum-delay wavelet, and the others are mixed-delay wavelets. However, these -length wavelets are not the only wavelets in the suite because other wavelets exist that have greater length but the same autocorrelation. Thus, those other wavelets are also members of the suite. If we consider all the members of the suite, there is only one minimum-delay wavelet, which is the (N + 1)-length minimum-delay wavelet that we have found. Thus, the minimum-delay concept is unique.
In contrast, the maximum-delay concept is a relative one, because a wavelet is a maximum-delay wavelet with respect to its length. In our example, the wavelet (0.25, 0.5, 1, 2) is the maximum-delay four-length wavelet of the suite. By delaying this wavelet by one time unit, we obtain the wavelet (0, 0.25, 0.5, 1, 2), which is the maximum-delay five-length wavelet of the suite. Similarly, (0, 0, 0.25, 0.5, 1, 2) is the maximum-delay six-length wavelet of the suite, and so forth. All the other wavelets of the suite - that is, the wavelets other than the minimum-delay wavelet and the maximum-delay wavelets - are mixed-delay wavelets.
The Z-transform of a filter’s impulse response is called the filter’s transfer function. A stable causal time-invariant linear filter is said to be (strictly) a minimum-delay filter (or a minimum-phase filter) if the zeros (or roots) of its transfer function lie outside the unit circle.
The transfer function of a causal FIR filter is the polynomial
( )
This polynomial is called a minimum-delay polynomial if the filter is minimum delay. It is called a maximum-delay polynomial if the filter is maximum delay. It is called a mixed-delay polynomial if the filter is mixed delay.
The inverse of a minimum-delay filter A(Z) is , which also is a minimum-delay filter. An important type of stable causal IIR (infinite-impulse-response) filter is the recursive filter, which has its transfer function given by the rational function with the provision that the denominator polynomial A(Z) is a minimum-delay polynomial (i.e., all the roots of the polynomial lie outside of the unit circle). This provision ensures that the recursive filter is stable (Treitel and Robinson, 1964^{[1]}). The locations of the zeros of the numerator polynomial B(Z) do not need to be specified for the recursive filter to be stable. However, if it is further specified that the numerator polynomial B(Z) is a minimum-delay polynomial, the recursive filter is a minimum-delay filter.
References
- ↑ Treitel, S., and E. A. Robinson, 1964, The stability of digital filters: IEEE Transactions on Geoscience Electronics, GE-2, 6–18.
Continue reading
Previous section | Next section |
---|---|
Illustrations of spectra | Energy |
Previous chapter | Next chapter |
Frequency | Synthetics |
Also in this chapter
- Wavelets
- Fourier transform
- Z-transform
- Delay: Minimum, mixed, and maximum
- Two-length wavelets
- Illustrations of spectra
- Energy
- Autocorrelation
- Canonical representation
- Zero-phase wavelets
- Symmetric wavelets
- Ricker wavelet
- Appendix G: Exercises