Hankel Low-Rank Approximation for Seismic Noise Attenuation

The low-rankness property of the Hankel matrix formulated from the clean seismic data corresponding to a few number of linear events has been successively leveraged in many low-rank (LR) approximation methods for seismic data denoising. The common scheme in these rank-reduction methods is to compute the best LR approximation of the formulated Hankel matrix and then obtain the denoised data from the LR matrix. However, without utilizing the Hankel structure when computing the LR approximation, if we rearrange the denoised data into a Hankel matrix, it is in general not exactly LR as expected. In this paper, we propose a Hankel LR (HLR) approximation method to simultaneously exploit both the Hankel structure and the LR property underlying the clean seismic data. The formulated HLR approximation problem is solved by an alternating-minimization-based algorithm. We provide rigorously convergence analysis of the proposed algorithm. The superior performance of the proposed HLR approximation method is demonstrated on both synthetic and field seismic data.