# Nonminimum-delay wavelet

Series | Geophysical References Series |
---|---|

Title | Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing |

Author | Enders A. Robinson and Sven Treitel |

Chapter | 9 |

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

ISBN | 9781560801481 |

Store | SEG Online Store |

What is the structure of a nonminimum-delay wavelet? Let the signature be the real-value, finite-length, nonminimum-delay wavelet . Let us call anything that has a flat magnitude spectrum white and anything that has a curved magnitude spectrum colored. The canonical representation is the key to signature deconvolution. The *canonical representation* states that a given nonminimum-delay wavelet *s* is equal to the convolution of (1) its minimum-delay counterpart *b*, which is a wavelet with the same color as the given nonminimum-delay wavelet *s* but with a minimum-phase spectrum, and (2) the all-pass wavelet *p*, which has a white (i.e., flat) magnitude spectrum and which carries the extra phase. The canonical representation of the signature is given by

**(**)

The signature *s* is known. The problem is to find its two components.

Let be the Fourier transform of the signature. The energy spectrum is defined as

**(**)

If we let , we can rewrite this as

**(**)

The autocorrelation is defined as

**(**)

The autocorrelation is symmetric; that is, The energy spectrum can be written as

**(**)

Thus, if we are given the signature *s*, we can compute its energy spectrum as well as its autocorrelation.

We now want to consider the *inverse problem*: Given the energy spectrum (or, equivalently, the autocorrelation), we want to find the wavelet that yields this energy spectrum. This inverse problem, as it stands, is not unique in that many wavelets have the given amplitude spectrum Only one of these wavelets is minimum delay. This exceptional wavelet can be determined because it is possible to determine the minimum-phase spectrum from knowledge of its amplitude spectrum . Thus, we have enough information to compute the minimum-delay counterpart of the nonminimum-delay signature. This is how we do it:

We write the energy spectrum as the *Z*-transform

**(**)

We see that is a polynomial of degree 2*N*. Because the autocorrelation is symmetric, it follows that if and only if . Thus, *Z* is a root of this polynomial if and only if is a root. Moreover, because the polynomial has real coefficients, it follows that for every complex root, the corresponding complex-conjugate root must occur. For example, suppose the autocorrelation is (2, 5, 2), where the center point 5 is at time index 0. The *Z*-transform of the autocorrelation is . The polynomial can be factored as . We see that -2 is a root and that is also a root.

Let us compute the roots of the polynomial. Any root of modulus one gives rise to an equal-delay component. Such situations are treated in Robinson (1967b)^{[1]}. Here, we assume that there are no roots of modulus one. There are 2*N* roots in all. Of these roots, *N* of them, say, , will have modulus greater than one. The remaining *N* roots will have modulus less than one. Let us form the polynomial

**(**)

where the constant is determined by requiring that the wavelet *b* have the same energy as the signature *s*. In other words, the constant is determined by the requirement that

**(**)

The wavelet *b* is thus the desired minimum-delay counterpart of the signature *s*.

The inverse of the minimum-delay counterpart *b* can be obtained by carrying out the polynomial division

**(**)

Let us now describe a computing method to obtain the coefficients . Write equation **40** as

**(**)

If we let , we can rearrange terms to give

**(**)

For this equation to hold, we see that

**(**)

Right away, we find the value . We shall put this value in the equation for and solve the equation for . Then we put these values in the equation for and solve the equation for . We continue in this way to solve the equations to obtain the inverse .

Next, we use the expression

**(**)

to obtain an expression for the *all-pass* filter *p*. Once this all-pass operator has been found, its inverse follows without further computation by simply reversing the order of its coefficients. That is, the inverse of *p* is simply its reverse (with respect to the time index 0):

**(**)

Thus, we have obtained the components of the canonical representation of the signature, and the components of the signature inverse, namely

**(**)

## References

- ↑ Robinson, E. A., 1967b, Predictive decomposition of time series with application to seismic exploration: Geophysics,
**32**, 418-484.

## Continue reading

Previous section | Next section |
---|---|

Convolutional model | Signature deconvolution |

Previous chapter | Next chapter |

Synthetics | Deconvolution |

## Also in this chapter

- Wavelets
- The shaping filter
- Spiking filter
- White convolutional model
- Wavelet processing
- All-pass filter
- Convolutional model
- Signature deconvolution
- Vibroseis
- Dual-sensor wavelet estimation
- Deconvolution: Einstein or predictive?
- Summary
- Appendix I: Exercises