# American Institute of Mathematical Sciences

doi: 10.3934/ipi.2020052

## Two-dimensional seismic data reconstruction using patch tensor completion

 1 School of Mathematics and Physics, China University of Geosciences, Wuhan 430074, China 2 School of Computer, Central China Normal University, Wuhan 430079, China

* Corresponding author: Lihua Fu

Received  October 2019 Revised  May 2020 Published  August 2020

Seismic data are often undersampled owing to physical or financial limitations. However, complete and regularly sampled data are becoming increasingly critical in seismic processing. In this paper, we present an efficient two-dimensional (2D) seismic data reconstruction method that works on texture-based patches. It performs completion on a patch tensor, which folds texture-based patches into a tensor. Reconstruction is performed by reducing the rank using tensor completion algorithms. This approach differs from past methods, which proceed by unfolding matrices into columns and then applying common matrix completion approaches to deal with 2D seismic data reconstruction. Here, we first re-arrange the seismic data matrix into a third-order patch tensor, by stacking texture-based patches that are divided from seismic data. Then, the seismic data reconstruction problem is formulated into a low-rank tensor completion problem. This formulation avoids destroying the spatial structure, and better extracts the underlying useful information. The proposed method is efficient and gives an improved performance compared with traditional approaches. The effectiveness of our patch tensor-based framework is validated using two classical tensor completion algorithms, low-rank tensor completion (LRTC), and the parallel matrix factorization algorithm (TMac), on both synthetic and field data experiments.

Citation: Qun Liu, Lihua Fu, Meng Zhang, Wanjuan Zhang. Two-dimensional seismic data reconstruction using patch tensor completion. Inverse Problems & Imaging, doi: 10.3934/ipi.2020052
##### References:

show all references

##### References:
Illustration of texture-patch pre-transformation and patch tensor pre-transformation
Singular value plots for the mode-1, mode-2, and mode-3 unfolding matrices of patch tensor of the field data. (a) field data. (b) the singular value plots for the mode-1 unfolding matrix of the patch tensor of field data (blue) and field data with missing columns (red star). (c) the singular value plots for the mode-2 unfolding matrix of the patch tensor of field data (blue) and field data with missing columns (red star). (d) the singular value plots for the mode-3 unfolding matrix of the patch tensor of field data (blue) and field data with missing columns (red star)
Illustration of unfolding a 3D tensor into a matrix and inversely folding a matrix into a 3D tensor
Synthetic seismic data recovery. (a) Original seismic data. (b) Corrupted data with 50% randomly missing traces. (c) Recovered signal via APG. (d) Recovered signal via LMaFit. (e)Recovered signal via LRTC. (f) Recovered signal via TMac
Comparison of the 9th trace taken from original synthetic data and reconstructed data processed by different methods: (a) APG, (b) LMaFit, (c) LRTC, and (d) TMac
Comparison of f-k spectra of synthetic data. (a) Original data. (b) Corrupted data. (c) Recovered data by APG. (d) Recovered data using the LMaFit algorithm. (e) Recovered data using the LRTC algorithm. (f) Recovered data using the TMac algorithm
Reconstruction results for the APG, LMaFit, LRTC, and TMac algorithms. (a) Signal-to-noise-ratio (SNR) versus the sampling ratio. (b) Computational time versus the sampling ratio
Reconstruction of post-stack seismic data. (a) Original data. (b) Corrupted data with 50% randomly missing traces
Comparison of the 58th trace taken from the original post-stack data and reconstructed data and the differences between them. Panels (a), (c), (e), and (g) show this reconstructed single trace using the APG, LMaFit, LRTC, and TMac algorithms, respectively. Panels (b), (d), (f), and (h) show the differences between the original single trace and the reconstructed single traces when using the APG, LMaFit, LRTC, and TMac algorithms, respectively
Comparison of f-k spectra for post-stack seismic data. (a) Original data. (b) Corrupted data. (c) Data recovered using the APG algorithm. (d) Data recovered using the LMaFit algorithm. (e) Data recovered using the LRTC algorithm. (f) Data recovered using the TMac algorithm
Reconstruction results of the post-stack seismic data. (a) The SNR versus the sampling ratio for the APG, LMaFit, LRTC, and TMac algorithms. (b) Amplitude spectrum comparison of the 58th single trace. The spectrums of the original data, data from the LMaFit and data from TMac are displayed in black, red, and blue lines, respectively
The reconstructed SNR for post-stack seismic data with the increase of iteration numbers in different patch size by LRTC
 [1] Yangyang Xu, Ruru Hao, Wotao Yin, Zhixun Su. Parallel matrix factorization for low-rank tensor completion. Inverse Problems & Imaging, 2015, 9 (2) : 601-624. doi: 10.3934/ipi.2015.9.601 [2] Yi Yang, Jianwei Ma, Stanley Osher. Seismic data reconstruction via matrix completion. Inverse Problems & Imaging, 2013, 7 (4) : 1379-1392. doi: 10.3934/ipi.2013.7.1379 [3] Wei Wan, Weihong Guo, Jun Liu, Haiyang Huang. Non-local blind hyperspectral image super-resolution via 4d sparse tensor factorization and low-rank. Inverse Problems & Imaging, 2020, 14 (2) : 339-361. doi: 10.3934/ipi.2020015 [4] Juan Manuel Reyes, Alberto Ruiz. Reconstruction of the singularities of a potential from backscattering data in 2D and 3D. Inverse Problems & Imaging, 2012, 6 (2) : 321-355. doi: 10.3934/ipi.2012.6.321 [5] Ke Wei, Jian-Feng Cai, Tony F. Chan, Shingyu Leung. Guarantees of riemannian optimization for low rank matrix completion. Inverse Problems & Imaging, 2020, 14 (2) : 233-265. doi: 10.3934/ipi.2020011 [6] Joshua Hudson, Michael Jolly. Numerical efficacy study of data assimilation for the 2D magnetohydrodynamic equations. Journal of Computational Dynamics, 2019, 6 (1) : 131-145. doi: 10.3934/jcd.2019006 [7] Thomas März, Andreas Weinmann. Model-based reconstruction for magnetic particle imaging in 2D and 3D. Inverse Problems & Imaging, 2016, 10 (4) : 1087-1110. doi: 10.3934/ipi.2016033 [8] Zhouchen Lin. A review on low-rank models in data analysis. Big Data & Information Analytics, 2016, 1 (2&3) : 139-161. doi: 10.3934/bdia.2016001 [9] Tomáš Oberhuber, Tomáš Dytrych, Kristina D. Launey, Daniel Langr, Jerry P. Draayer. Transformation of a Nucleon-Nucleon potential operator into its SU(3) tensor form using GPUs. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020383 [10] Yanfei Wang, Dmitry Lukyanenko, Anatoly Yagola. Magnetic parameters inversion method with full tensor gradient data. Inverse Problems & Imaging, 2019, 13 (4) : 745-754. doi: 10.3934/ipi.2019034 [11] Magdalena Czubak, Nina Pikula. Low regularity well-posedness for the 2D Maxwell-Klein-Gordon equation in the Coulomb gauge. Communications on Pure & Applied Analysis, 2014, 13 (4) : 1669-1683. doi: 10.3934/cpaa.2014.13.1669 [12] Adam Larios, Yuan Pei. Approximate continuous data assimilation of the 2D Navier-Stokes equations via the Voigt-regularization with observable data. Evolution Equations & Control Theory, 2020, 9 (3) : 733-751. doi: 10.3934/eect.2020031 [13] Jan Boman, Vladimir Sharafutdinov. Stability estimates in tensor tomography. Inverse Problems & Imaging, 2018, 12 (5) : 1245-1262. doi: 10.3934/ipi.2018052 [14] Fan Wu. Conditional regularity for the 3D Navier-Stokes equations in terms of the middle eigenvalue of the strain tensor. Evolution Equations & Control Theory, 2020  doi: 10.3934/eect.2020078 [15] Gianluca Crippa, Elizaveta Semenova, Stefano Spirito. Strong continuity for the 2D Euler equations. Kinetic & Related Models, 2015, 8 (4) : 685-689. doi: 10.3934/krm.2015.8.685 [16] Bernd Kawohl, Guido Sweers. On a formula for sets of constant width in 2d. Communications on Pure & Applied Analysis, 2019, 18 (4) : 2117-2131. doi: 10.3934/cpaa.2019095 [17] Ka Kit Tung, Wendell Welch Orlando. On the differences between 2D and QG turbulence. Discrete & Continuous Dynamical Systems - B, 2003, 3 (2) : 145-162. doi: 10.3934/dcdsb.2003.3.145 [18] Julien Cividini. Pattern formation in 2D traffic flows. Discrete & Continuous Dynamical Systems - S, 2014, 7 (3) : 395-409. doi: 10.3934/dcdss.2014.7.395 [19] Géry de Saxcé, Claude Vallée. Structure of the space of 2D elasticity tensors. Discrete & Continuous Dynamical Systems - S, 2013, 6 (6) : 1525-1537. doi: 10.3934/dcdss.2013.6.1525 [20] Igor Kukavica, Amjad Tuffaha. On the 2D free boundary Euler equation. Evolution Equations & Control Theory, 2012, 1 (2) : 297-314. doi: 10.3934/eect.2012.1.297

2019 Impact Factor: 1.373

## Tools

Article outline

Figures and Tables