Least-squares prediction and smoothing
|Series||Geophysical References Series|
|Title||Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing|
|Author||Enders A. Robinson and Sven Treitel|
|Store||SEG Online Store|
An important operation is signal extrapolation or, in other words, prediction. Smoothing is yet another important signal operation. Often, the desired quantity is observable only after it has been corrupted or altered by mixture — perhaps additively — with other signals. Such a corrupting signal commonly is called noise, even for cases in which the signal has a physically legitimate presence. For example, reverberations often are regarded as noise, but reverberations are a physical necessity. Such signal-generated noise contains valuable information in its own right, and the trick is to use this information to an advantage whenever possible.
Often it is important to ascertain, in the least-squares sense, what the data would have been like without contamination by noise. Smoothing might be the complete problem to be addressed; alternatively, it might be combined with a prediction problem, which means that we wish to know what the uncontaminated signal will do in the future. Whereas the smoothing problem is clearly distinguishable from the prediction problem, mixed problems involving elements of both are greatly important. Indeed, good smoothing performance usually depends on introduction of a sufficient delay. If the delay is negative, filter performance suffers. On the other hand, the filter becomes a smoothing predictor, which often is a useful tool. Seismic processing must be innovative. Most geophysical data-processing methods are mathematical and physical hybrids and are based on a particular geophysical model.
Before an actual digital filter can be designed, it is necessary to build a model. The time-honored least-squares model requires three signals: the input signal , the actual output signal , and the desired output signal . The filter is a finite-impulse-response (FIR) filter. Let such a filter k have a finite number of coefficients, given by
This filter is said to be of length N. The output is the convolution of the filter with the input; that is, . More specifically, the output value at time n is given by the discrete convolution formula
The input signal and the desired output signal are known. The problem is to find the filter by the least-squares criterion. At time instant n, the value of the output is and the value of the desired output is . The difference between the desired output and the actual output is the error . The quantity
is the error energy. The filter f that minimizes the value of the error energy is the optimum filter in the least-squares sense. Setting its partial derivatives with respect to each of the filter coefficients equal to zero minimizes the error energy. The result of that minimization is the set of a system of N linear simultaneous equations known as the normal equations. In their matrix form, the normal equations are
which were derived in detail in Chapter 9.
The known quantities in this system are the autocorrelation of the input signal and the positive-lag values of crosscorrelation of the desired output signal with the input signal. The unknown quantities are the values of the filter coefficients. The autocorrelation of the input is
whereas the crosscorrelation between the desired output and the input is
The solution of the normal equations (equations 4) yields the filter coefficients. The square matrix in the normal equations has the autocorrelation coefficients arranged in Toeplitz form. That is, all the autocorrelation coefficients along any diagonal (the main diagonal or any subdiagonal) are the same. Because of this Toeplitz structure, the well-known Levinson algorithm can be used in the solution of the normal equations (Robinson and Treitel, 2000, p. 163-169).
Let us now derive the least-squares prediction filter, and let denote this prediction filter. Such a filter uses the input signal’s past values to predict that signal’s future values. The input is the signal , and the desired output is the time-advanced version of the input. Let the prediction distance be given by the positive integer . The input to the filter is the input at present time n. The desired output is the input at the future time . The prediction filter is designed so that the output at the present time n is an optimum estimate of the future value . If this estimate is denoted by , the filter’s action can be represented by the convolution
Equation 7 says that the prediction operator acts on an input up to time n and estimates its value at some future time . Normal equations 4 now can be used.
The crosscorrelation between the desired output and the input is given by the positive lag coefficients of the crosscorrelation between the time-advanced trace and the trace
As equation 8 shows, the crosscorrelation between the desired output and the input is equal to the autocorrelation of the input for lags greater than or equal to . Thus, the normal equations for the prediction filter can be written as
Their solution yields the coefficients of the optimum prediction filter we seek.
- ↑ Robinson, E. A., and S. Treitel, 2000, Geophysical signal analysis: SEG.
|Previous section||Next section|
|Model used for deconvolution||The prediction-error filter|
|Previous chapter||Next chapter|
|Wavelet Processing||Fine Points|
Also in this chapter
- Model used for deconvolution
- The prediction-error filter
- Spiking deconvolution
- Gap deconvolution
- Tail shaping and head shaping
- Seismic deconvolution
- Piecemeal convolutional model
- Time-varying convolutional model
- Random-reflection-coefficient model
- Implementing deconvolution
- Canonical representation
- Appendix J: Exercises