Published as Geophysical Journal International, 222, 1824-1845, (2020)

Five-dimensional seismic data reconstruction using the optimally damped rank-reduction method

Yangkang Chen% latex2html id marker 3434
\setcounter{footnote}{1}\fnsymbol{footnote}% latex2html id marker 3434
ootnote}{1}<SPAN CLASS=% latex2html id marker 3435
\setcounter{{}{0}\fnsymbol{{}, Min Bai% latex2html id marker 3436
\setcounter{footnote}{1}\fnsymbol{footnote}, Zhe Guan% latex2html id marker 3437
\setcounter{footnote}{2}\fnsymbol{footnote}, Qingchen Zhang% latex2html id marker 3438
\setcounter{footnote}{1}\fnsymbol{footnote}, Mi Zhang% latex2html id marker 3439
\setcounter{footnote}{3}\fnsymbol{footnote}, and Hang Wang% latex2html id marker 3440
\setcounter{footnote}{1}\fnsymbol{footnote}% latex2html id marker 3440
ootnote}{1}<SPAN CLASS=% latex2html id marker 3441
% latex2html id marker 3442
\setcounter{footnote}{1}\fnsymbol{footnote}School of Earth Sciences
Zhejiang University
Hangzhou, Zhejiang Province, China, 310027
Applied Physics Program
Rice University
Houston, TX, USA
State Key Laboratory of Petroleum Resources and Prospecting
China University of Petroleum
Fuxue Road 18th
Beijing, China, 102200


It is difficult to separate additive random noise from spatially coherent signal using a rank-reduction method that is based on the truncated singular value decomposition (TSVD) operation. This problem is due to the mixture of the signal and the noise subspaces after the TSVD operation. This drawback can be partially conquered using a damped rank reduction method, where the singular values corresponding to effective signals are adjusted via a carefully designed damping operator. The damping operator works most powerfully in the case of a small rank and a small damping factor. However, for complicated seismic data, e.g., multi-channel reflection seismic data containing highly curved events, the rank should be large enough to preserve the details in the data, which makes the damped rank reduction method less effective. In this paper, we develop an optimal damping strategy for adjusting the singular values when a large rank parameter is selected so that the estimated signal can best approximate the exact signal. We first weight the singular values using optimally calculated weights. The weights are theoretically derived by solving an optimization problem that minimizes the Frobenius-norm difference between the approximated signal components and the exact signal components. The damping operator is then derived based on the initial weighting operator to further reduce the residual noise after the optimal weighting. The resulted optimally damped rank reduction method is nearly an adaptive method, i.e., insensitive to the rank parameter. We demonstrate the performance of the proposed method on a group of synthetic and real five-dimensional seismic data.