On the application of the fast Kalman algorithm to adaptive deconvolution of seismic data

Abstract

The application of a recently proposed fast implementation of the recursive least squares algorithm, called the Fast Kalman Algorithm (FKA) to adaptive deconvolution of seismic data is discussed. The newly proposed algorithm does not require the storage and updating of a matrix to calculate the filter gain, and hence is computationally very efficient. Furthermore, it has an interesting structure yielding both the forward and backward prediction residuals of the seismic trace and thus permits the estimation of a "smoothed residual" without any additional computations. The paper also compares the new algorithm with the conventional Kalman algorithm (CKA) proposed earlier [3] for seismic deconvolution. Results of experiments on simulated as well as real data show that while the FKA is a little more sensitive to the choice of some initial parameters which need to be selected carefully, it can yield comparable performance with greatly reduced computational effort.