# Gap deconvolution of a mixed-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 | 11 |

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

ISBN | 9781560801481 |

Store | SEG Online Store |

Much work has been done on the problem of deconvolution in the case of a mixed-delay (or nonminimum-delay) wavelet. See, for example, Eisner and Hampson (1990)^{[1]}; Ford (1978)^{[2]}; Lazear (1984)^{[3]}; Lines and Ulrych (1977)^{[4]}; Porsani and Ursin (1998)^{[5]}; Sacchi and Ulrych (2000)^{[6]}; Sengbush et al. (1987)^{[7]}; Tygel et al. (1991)^{[8]}; Ulrych et al. (1995)^{[9]}; and Yilmaz (1987)^{[10]}. Here, we give another approach.

Suppose various signature-deconvolution procedures are used to remove the source signature, absorption, and instrument response from a seismic trace. However, the process is not perfect, so a short causal nonminimum-delay residual wavelet *u* remains in their place. Thus, the resulting trace is , where *m* is the minimum-delay multiple and is the white reflectivity. Suppose the multiple consists only of the water-layer reverberation given above — that is, . For example, suppose the reverberation is

**(**)

The seismic wavelet *w* is defined as the convolution of the residual wavelet and the minimum-delay reverberation train. Thus, the seismic trace is . If the residual wavelet is the maximum-delay wavelet , then the seismic wavelet is the mixed-delay wavelet

**(**)

By inspection, the prediction operator (for a prediction distance of two) for the seismic wavelet is . We can verify that this indeed is the prediction operator by carrying out the convolution

**(**)

which is the causal part of the advanced seismic wavelet. Thus, the prediction-error operator is which, in fact, is the dereverberation operator. We carry out the convolution

**(**)

which is the residual wavelet *u*. If we deconvolve the seismic trace with this prediction-error operator, we obtain

**(**)

Thus, the deconvolved trace is the reflectivity smoothed by the residual wavelet.

By wavelet processing, this residual wavelet can be shaped into a desirable interpreter wavelet. We note that the deconvolution operator *f* is equal to the dereverberation . In effect, the nonminimum-delay residual wavelet falls thought the cracks in the gap-deconvolution process. This example illustrates that in the case of a nonminimum-delay seismic wavelet, gap deconvolution can be used judiciously in the correct circumstances to remove the reverberations.

Let us now give another example of gap deconvolution. We shall illustrate the method of predictive deconvolution for the case of a water layer, which produces a reverberation train. Suppose the water depth is 100 ft (30 m). Let the water-surface reflection coefficient be unity, and let the water-bottom reflection coefficient be 0.5. Let the velocity of the water be 5000 ft (830 m)/s. Then the two-way normal-incidence traveltime in the water layer is 2(100)/5000 = 0.040 s, which is 40 ms. Assume that the sampling increment is 4 ms, so that the reverberation cycle time in discrete units is 40/4 or 10.

Figure 9a shows the components of the convolutional model . Figure 9b depicts the prediction-error operator (PEO) *f* for a prediction distance equal to the cycle time, the deconvolved trace , the theoretical dereverberation operator *a*, and the smoothed reflectivity . How well the least-squares deconvolution worked can be determined as follows. First compare the prediction-error operator (PEO) with the dereverberation operator. The closer they are, the better the deconvolution is. Second, compare the deconvolved trace with the smoothed reflectivity. Again, the closer they are, the better the deconvolution is. Although there are differences in the details, the agreement is quite good.

Deconvolution removes the reverberation but leaves the nonminimum-phase residual wavelet in place (if the length of the residual wavelet is less than the reverberation cycle time). The important point is that the method of gap deconvolution works even if the residual wavelet is not minimum delay but is instead mixed delay. It is necessary to use the correct value of the prediction distance (namely, the reverberation cycle time) in the design of the prediction-error operator.

How does one determine the prediction distance? One way is to compute the autocorrelation of the seismic trace and see if it exhibits periodicity. If it does, the cycle time can be set equal to that period. A set of prediction distances clustering around such an estimate of cycle time should be used to determine empirically the “best-looking” dereverberation operator. Within statistical accuracy, such a method of predictive deconvolution can be described as follows: The prediction-error filter operating on the trace wipes out the reverberations so that the reverberations do not appear in the deconvolved trace. The residual wavelet (because its length is less than the cycle time) falls through the cracks of the prediction-error operator. Hence, the source wavelet does appear in the deconvolved trace. This deconvolved trace is equal to the reflectivity smoothed by the residual wavelet.

## References

- ↑ Eisner, E., and G. Hampson, 1990, Decomposition into minimum and maximum phase components: Geophysics,
**55**, 897-901. Discussion and reply: 1993, Geophysics,**58**, 1207-1213. - ↑ Ford, W. T., 1978, Optimum mixed delay spiking filters: Geophysics,
**43**, 125-132. - ↑ Lazear, G. D., 1984, An examination of the exponential decay method of mixed-phase wavelet estimation: Geophysics,
**49**, 2094-2099. - ↑ Lines, L. R., and T. J. Ulrych, 1977, The old and the new in seismic deconvolution and wavelet estimation: Geophysical Prospecting,
**25**, 512-540. - ↑ Porsani, M. J., and B. Ursin, 1998, Mixed-phase deconvolution: Geophysics,
**63**, 637-647. - ↑ Sacchi, M. D., and T. J. Ulrych, 2000, Nonminimum-phase wavelet estimation using higher order statistics: The Leading Edge,
**19**, no. 1, p. 80-83. - ↑ Sengbush, R. L., M. Hato, and H. Chang, 1987, Optimal compensation of time variance and non-minimum phase in Wiener-Robinson deconvolution,
*in*M. Bernabini, P. Carrion, G. Jacovitti, F. Rocca, S. Treitel, and M. Worthington, eds., Deconvolution and inversion: Blackwell Scientific Publications, 90-114. - ↑ Tygel, M., H. Huck, and P. Hubral, 1991, Mixed-delay wavelet deconvolution of the point source seismogram: Geophysics,
**56**, 1405-1411. - ↑ Ulrych, T. J., D. R. Velis, and M. D. Sacchi, 1995, Wavelet estimation revisited: The Leading Edge,
**14**, 1139-1143. - ↑ Yilmaz, Ö., 1987, Seismic data processing: SEG Investigations in Geophysics No. 2.

## Continue reading

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

Water reverberations | Prewhitening |

Previous chapter | Next chapter |

Deconvolution | Attributes |

## Also in this chapter

- Prediction-error filters
- Water reverberations
- Prewhitening
- Prediction distance
- Model-driven predictive deconvolution
- Convolutional model in the frequency domain
- Time-variant spectral whitening
- Model-based deconvolution
- Surface-consistent deconvolution
- Interactive earth-digital processing
- Appendix K: Exercises