Deconvolución Gap de un ondícula de retardo mixto

From SEG Wiki
Jump to navigation Jump to search
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.
This page is a translated version of the page Gap deconvolution of a mixed-delay wavelet and the translation is 50% complete.
Other languages:
Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing
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
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 9.  (a) The components of the convolutional model of a reverberating trace. (b) The prediction-error operator, the deconvolved trace, the dereverberation operator, and the smoothed reflectivity for the trace shown in Figure 7a.

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.


  1. 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.
  2. Ford, W. T., 1978, Optimum mixed delay spiking filters: Geophysics, 43, 125-132.
  3. Lazear, G. D., 1984, An examination of the exponential decay method of mixed-phase wavelet estimation: Geophysics, 49, 2094-2099.
  4. Lines, L. R., and T. J. Ulrych, 1977, The old and the new in seismic deconvolution and wavelet estimation: Geophysical Prospecting, 25, 512-540.
  5. Porsani, M. J., and B. Ursin, 1998, Mixed-phase deconvolution: Geophysics, 63, 637-647.
  6. 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.
  7. 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.
  8. Tygel, M., H. Huck, and P. Hubral, 1991, Mixed-delay wavelet deconvolution of the point source seismogram: Geophysics, 56, 1405-1411.
  9. Ulrych, T. J., D. R. Velis, and M. D. Sacchi, 1995, Wavelet estimation revisited: The Leading Edge, 14, 1139-1143.
  10. Yilmaz, Ö., 1987, Seismic data processing: SEG Investigations in Geophysics No. 2.

Sigue leyendo

Sección previa Siguiente sección
Reverberaciones en agua Preblanqueo
Capítulo previo Siguiente capítulo
Deconvolución Atributos

Tabla de contenido

También en este capítulo

Vínculos externos

find literature about
Gap deconvolution of a mixed-delay wavelet/es
SEG button search.png Datapages button.png GeoScienceWorld button.png OnePetro button.png Schlumberger button.png Google button.png AGI button.png