| Abstract: | On seismograms recorded at sea bubble pulse oscillations can present a serious problem to an interpreter. We propose a new approach, based on generalized linear inverse theory, to the solution of the debubbling problem. Under the usual assumption that a seismogram can be modelled as the convolution of the earth's impulse response and a source wavelet we show that estimation of either the wavelet or the impulse response can be formulated as a generalized linear inverse problem. This parametric approach involves solution of a system of equations by minimizing the error vector (ΔX = Xobs– Xcal) in a least squares sense. One of the most significant results is that the method enables us to control the accuracy of the solution so that it is consistent with the observational errors and/or known noise levels. The complete debubbling procedure can be described in four steps: (1) apply minimum entropy deconvolution to the observed data to obtain a deconvolved spike trace, a first approximation to the earth's response function; (2) use this trace and the observed data as input for the generalized linear inverse procedure to compute an estimated basic bubble pulse wavelet; (3) use the results of steps 1 and 2 to construct the compound source signature consisting of the primary pulse plus appropriate bubble oscillations; and (4) use the compound source signature and the observed data as input for the generalized linear inverse method to determine the estimated earth impulse response—a debubbled, deconvolved seismogram. We illustrate the applicability of the new approach with a set of synthetic seismic traces and with a set of field seismograms. A disadvantage of the procedure is that it is computationally expensive. Thus it may be more appropriate to apply the technique in cases where standard analysis techniques do not give acceptable results. In such cases the inherent advantages of the method may be exploited to provide better quality seismograms. |
