Adaptive singular spectrum analysis for seismic denoising and interpolation

Abstract

We have developed an adaptive singular spectrum analysis (ASSA) method for seismic data denoising and interpolation purposes. Our algorithm iteratively updates the singular-value decomposition (SVD) of current spatial patches using the most recently added spatial sample. The method reduces the computational cost of classic singular spectrum analysis (SSA) by requiring QR decompositions on smaller matrices rather than the factorization of the entire Hankel matrix of the data. A comparison between results obtained by the ASSA and SSA methods, in which the SVD applies to all of the traces at once, proves that the ASSA method is a valid way to cope with spatially varying dips. In addition, a comparison of the ASSA method with the windowed SSA method indicates gains in efficiency and accuracy. Synthetic and real data examples illustrate the effectiveness of our method.