Multidimensional de-aliased Cadzow reconstruction of seismic records

ABSTRACTWe tested a strategy for beyond-alias interpolation of seismic data using Cadzow reconstruction. The strategy enables Cadzow reconstruction to be used for interpolation of regularly sampled seismic records. First, in the frequency-space (f-x) domain, we generated a Hankel matrix from the spatial samples of the low frequencies. To perform interpolation at a given frequency, the spatial samples were interlaced with zero samples and another Hankel matrix was generated from the zero-interlaced data. Next, the rank-reduced eigen-decomposition of the Hankel matrix at low frequencies was used for beyond-alias preconditioning of the Hankel matrix at a given frequency. Finally, antidiagonal averaging of the conditioned Hankel matrix produced the final interpolated data. In addition, the multidimensional extension of the proposed algorithm was explained. The proposed method provides a unifying thread between reduced-rank Cadzow reconstruction and beyond alias f-x prediction error interpolation. Synthetic an...

[1]  S. Spitz.  Seismic trace interpolation in the F-X domain , 1991 .

[2]  Mauricio D. Sacchi,et al.  Multistep autoregressive reconstruction of seismic records , 2007 .

[3]  Jon F. Claerbout,et al.  An algorithm for interpolation in the pyramid domain , 2010 .

[4]  Mauricio D. Sacchi,et al.  Multifrequency singular spectrum analysis , 2009 .

[5]  S. Trickett.  F-xy eigenimage noise suppression , 2003 .

[6]  Lynn Burroughs,et al.  Rank-Reduction-Based Trace Interpolation , 2010 .

[7]  Kristopher A. Innanen,et al.  Seismic data interpolation using a fast generalized Fourier transform , 2011 .

[8]  Per-Gunnar Martinsson,et al.  Randomized algorithms for the low-rank approximation of matrices , 2007, Proceedings of the National Academy of Sciences.

[9]  Mauricio D. Sacchi,et al.  On sampling functions and Fourier reconstruction methods , 2010 .

[10]  Milton J. Porsani,et al.  Seismic trace interpolation using half-step prediction filters , 1999 .

[11]  F. Herrmann,et al.  Simply denoise: Wavefield reconstruction via jittered undersampling , 2008 .

[12]  M. Sacchi,et al.  Simultaneous seismic data denoising and reconstruction via multichannel singular spectrum analysis , 2011 .

[13]  Mauricio D. Sacchi,et al.  A Fast Rank Reduction Method For the Reconstruction of 5D Seismic Volumes , 2011 .

[14]  Mauricio D. Sacchi,et al.  f-x adaptive seismic-trace interpolation , 2009 .

[15]  Necati Gulunay,et al.  Seismic trace interpolation in the Fourier transform domain , 2003 .

[16]  R. Abma,et al.  3D interpolation of irregular data with a POCS algorithm , 2006 .

[17]  Mauricio D. Sacchi,et al.  Minimum weighted norm interpolation of seismic records , 2004 .

[18]  Mauricio D. Sacchi,et al.  Beyond alias hierarchical scale curvelet interpolation of regularly and irregularly sampled seismic data , 2010 .