Predictive deconvolution

ADVERTISEMENT
From SEG Wiki
Revision as of 17:15, 7 October 2014 by Ageary (talk | contribs) (added link)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

The type 3 desired output, a time-advanced form of the input series, suggests a prediction process. Given the input x(t), we want to predict its value at some future time (t + α), where α is prediction lag. Wiener showed that the filter used to estimate x(t + α) can be computed by using a special form of the matrix equation (30) [1]. Since the desired output x(t + α) is the time-advanced version of the input x(t), we need to specialize the right side of equation (30) for the prediction problem.

Consider a five-point input time series x(t): (x0, x1, x2, x3, x4), and set α = 2. The autocorrelation of the input series is computed in Table 2-23, and the crosscorrelation between the desired output x(t + 2) and the input x(t) is computed in Table 2-24. Compare the results in Tables 2-23 and 2-24, and note that gi = ri+α for α = 2 and i = 0, 1, 2, 3, 4.

Equation (30), for this special case, is rewritten as follows:


(30)


(34)

The prediction filter coefficients a(t) : (a0, a1, a2, a3, a4) can be computed from equation (34) and applied to the input series x(t) : (x0, x1, x2, x3, x4) to compute the actual output y(t) : (y0, y1, y2, y3, y4) (Table 2-25). We want to predict the time-advanced form of the input; hence, the actual output is an estimate of the series x(t + α) : (x2, x3, x4), where α = 2. The prediction error series e(t) = x(t + α) − y(t) : (e2, e3, e4, e5, e6) is given in Table 2-26.

The results in Table 2-26 suggest that the error series can be obtained more directly by convolving the input series x(t) : (x0, x1, x2, x3, x4) with a filter with coefficients (1, 0, −a0, −a1, −a2, −a3, −a4) (Table 2-27). The results for (e2, e3, e4, e5, e6) are identical (Tables 2-26 and 2-27). Since the series (a0, a1, a2, a3, a4) is called the prediction filter, it is natural to call the series (1, 0, −a0, −a1, −a2, −a3, −a4) the prediction error filter. When applied to the input series, this filter yields the error series in the prediction process (Table 2-27).

Table 2-23. Autocorrelation lags of input series x(t) : (x0, x1, x2, x3, x4).
r1 = x0x1 + x1x2 + x2x3 + x3x4
r2 = x0x2 + x1x3 + x2x4
r3 = x0x3 + x1x4
r4 = x0x4
r5 = 0
r6 = 0
Table 2-24. Crosscorrelation of desired output x(t + α): (x2, x3, x4), α = 2, with input x(t) : (x0, x1, x2, x3, x4)
g0 = x0x2 + x1x3 + x2x4
g1 = x0x3 + x1x4
g2 = x0x4
g3 = 0
g4 = 0
Table 2-25. Convolution of prediction filter a(t): (a0, a1, a2, a3, a4) with input series x(t): (x0, x1, x2, x3, x4) to compute actual output y(t) : (y0, y1, y2, y3, y4).
y0 = a0x0
y1 = a1x0 + a0x1
y2 = a2x0 + a1x1 + a0x2
y3 = a3x0 + a2x1 + a1x2 + a0x3
y4 = a4x0 + a3x1 + a2x2 + a1x3 + a0x4
Table 2-26. The error series e(t) = x(t + α) − y(t) : (e2, e3, e4, e5, e6), α = 2. For y(t), see Table 2-25.
e2 = x2a0x0
e3 = x3a1x0a0x1
e4 = x4a2x0a1x1a0x2
e5 = 0 − a3x0a2x1a1x2a0x3
e6 = 0 − a4x0a3x1a2x2a1x3a0x4
Table 2-27. Convolution of prediction error filter coefficients (1, 0, −a0, −a1, −a2, −a3, −a4) with input series x(t) : (x0, x1, x2, x3, x4).
e0 = x0
e1 = x1
e2 = x2a0x0
e3 = x3a1x0a0x1
e4 = x4a2x0a1x1a0x2
e5 = 0 − a3x0a2x1a1x2a0x3
e6 = 0 − a4x0a3x1a2x2a1x3a0x4

Why place so much emphasis on the error series? Consider the prediction process as it relates to a seismic trace. From the past values of a time series up to time t, a future value can be predicted at time t + α, where α is the prediction lag. A seismic trace often has a predictable component (multiples) with a periodic rate of occurrence. According to assumption 6, anything else, such as primary reflections, is unpredictable.

Some may claim that reflections are predictable as well; this may be the case if deposition is cyclic. However, this type of deposition is not often encountered. While the prediction filter yields the predictable component (the multiples) of a seismic trace, the remaining unpredictable part, the error series, is essentially the reflection series.

Equation (34) can be generalized for the case of an n-long prediction filter and an α-long prediction lag.


(35)

Note that design of the prediction filters requires only autocorrelation of the input series.

There are two approaches to predictive deconvolution:

  1. The prediction filter (a0, a1, a2, …, an−1) may be designed using equation (35) and applied on input series as described in Figure 2.3-10.
  2. Alternatively, the prediction error filter (1, 0, 0, …, 0, −a0, −a1, −a2, …, −an−1) can be designed and convolved with the input series as described in Figure 2.3-11.

Now consider the special case of unit prediction lag, α = 1. For n = 5, equation (35) takes the following form:


(36)

By augmenting the right side to the left side, we obtain:


(37)

Add one row and move the negative sign to the column matrix that represents the filter coefficients to get:


(38)

where L = r0r1a0r2a1r3a2r4a3r5a4. Note that there are six unknowns, (a0, a1, a2, a4, a5, L), and six equations. Solution of these equations yields the unit-delay prediction error filter (1, −a0, −a1, −a2, −a4, −a5), and the quantity L — the error in the filtering process (Section B.5). We can rewrite equation (38) as follows:


(39)

where b0 = 1, bi = −ai, and i = 1, 2, 3, 4, 5. This equation has a familiar structure. In fact, except for the scale factor L, it has the same form as equation (31), which yields the coefficients for the least-squares zero-delay inverse filter. This inverse filter is therefore the same as the prediction error filter with unit prediction lag, except for a scale factor. Hence, spiking deconvolution actually is a special case of predictive deconvolution with unit prediction lag.


(31)

We now know that predictive deconvolution is a general process that encompasses spiking deconvolution. In general, the following statement can be made: Given an input wavelet of length (n + α), the prediction error filter contracts it to an α-long wavelet, where α is the prediction lag [2]. When α = 1, the procedure is called spiking deconvolution.

Figure 2.3-12 interrelates the various filters discussed in this chapter and indicates the kind of process they imply. From Figure 2.3-12, note that Wiener filters can be used to solve a wide range of problems. In particular, predictive deconvolution is an integral part of seismic data processing that is aimed at compressing the seismic wavelet, thereby increasing temporal resolution. In the limit, it can be used to spike the seismic wavelet and obtain an estimate for reflectivity.

References

  1. Robinson and Treitel, 1980, Robinson, E. A. and Treitel, S., 1980, Geophysical signal analysis: Prentice-Hall Book Co.
  2. Peacock and Treitel, 1969, Peacock, K. L. and Treitel, S., 1969, Predictive deconvolution — theory and practice: Geophysics, 34, 155–169.

See also

External links

find literature about
Predictive deconvolution
SEG button search.png Datapages button.png GeoScienceWorld button.png OnePetro button.png Schlumberger button.png Google button.png AGI button.png