Assuming seismic data in a suitable domain is low rank while missing traces or noises increase the rank of the data matrix,the rank⁃reduced methods have been applied successfully for seismic interpolation and denoising.These rank⁃reduced methods mainly include Cadzow reconstruction that uses eigen decomposition of the Hankel matrix in the f⁃x(frequency⁃spatial)domain,and nuclear⁃norm minimization(NNM)based on rigorous optimization theory on matrix completion(MC).In this paper,a low patch⁃rank MC is proposed with a random⁃overlapped texture⁃patch mapping for interpolation of regularly missing traces in a three⁃dimensional(3D)seismic volume.The random overlap plays a simple but important role to make the low⁃rank method effective for aliased data.It shifts the regular column missing of data matrix to random point missing in the mapped matrix,where the missing data increase the rank thus the classic low⁃rank MC theory works.Unlike the Hankel matrix based rank⁃reduced method,the proposed method does not assume a superposition of linear events,but assumes the data have repeated texture patterns.Such data lead to a low⁃rank matrix after the proposed texture⁃patch mapping.Thus the methods can interpolate the waveforms with varying dips in space.A fast low⁃rank factorization method and an orthogonal rank⁃one matrix pursuit method are applied to solve the presented interpolation model.The former avoids the singular value decomposition(SVD)computation and the latter only needs to compute the large singular values during iterations.The two fast algorithms are suitable for large⁃scale data.Simple averaging realizations of several results from different random⁃overlapped texture⁃patch mappings can further increase the reconstructed signal⁃to⁃noise ratio(SNR).Examples on synthetic data and field data are provided to show successful performance of the presented method.
展开▼
