# American Institute of Mathematical Sciences

December  2020, 14(6): 985-1000. 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, 2020, 14 (6) : 985-1000. 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] Meng Ding, Ting-Zhu Huang, Xi-Le Zhao, Michael K. Ng, Tian-Hui Ma. Tensor train rank minimization with nonlocal self-similarity for tensor completion. Inverse Problems & Imaging, 2021, 15 (3) : 475-498. doi: 10.3934/ipi.2021001 [2] Liangliang Ma. Stability of hydrostatic equilibrium to the 2D fractional Boussinesq equations. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021068 [3] José Raúl Quintero, Juan Carlos Muñoz Grajales. On the existence and computation of periodic travelling waves for a 2D water wave model. Communications on Pure & Applied Analysis, 2018, 17 (2) : 557-578. doi: 10.3934/cpaa.2018030 [4] Weihua Jiang, Xun Cao, Chuncheng Wang. Turing instability and pattern formations for reaction-diffusion systems on 2D bounded domain. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021085 [5] Peng Tong, Xiaogang Ma. Design of differentiated warranty coverage that considers usage rate and service option of consumers under 2D warranty policy. Journal of Industrial & Management Optimization, 2021, 17 (4) : 1577-1591. doi: 10.3934/jimo.2020035 [6] Xin-Guang Yang, Rong-Nian Wang, Xingjie Yan, Alain Miranville. Dynamics of the 2D Navier-Stokes equations with sublinear operators in Lipschitz-like domains. Discrete & Continuous Dynamical Systems, 2021, 41 (7) : 3343-3366. doi: 10.3934/dcds.2020408 [7] Habib Ammari, Josselin Garnier, Vincent Jugnon. Detection, reconstruction, and characterization algorithms from noisy data in multistatic wave imaging. Discrete & Continuous Dynamical Systems - S, 2015, 8 (3) : 389-417. doi: 10.3934/dcdss.2015.8.389 [8] Mingchao Zhao, You-Wei Wen, Michael Ng, Hongwei Li. A nonlocal low rank model for poisson noise removal. Inverse Problems & Imaging, 2021, 15 (3) : 519-537. doi: 10.3934/ipi.2021003 [9] Abdeslem Hafid Bentbib, Smahane El-Halouy, El Mostafa Sadek. Extended Krylov subspace methods for solving Sylvester and Stein tensor equations. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021026 [10] Xiaofei Liu, Yong Wang. Weakening convergence conditions of a potential reduction method for tensor complementarity problems. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021080 [11] Tao Wu, Yu Lei, Jiao Shi, Maoguo Gong. An evolutionary multiobjective method for low-rank and sparse matrix decomposition. Big Data & Information Analytics, 2017, 2 (1) : 23-37. doi: 10.3934/bdia.2017006 [12] Marion Darbas, Jérémy Heleine, Stephanie Lohrengel. Numerical resolution by the quasi-reversibility method of a data completion problem for Maxwell's equations. Inverse Problems & Imaging, 2020, 14 (6) : 1107-1133. doi: 10.3934/ipi.2020056 [13] Mirela Kohr, Sergey E. Mikhailov, Wolfgang L. Wendland. Dirichlet and transmission problems for anisotropic stokes and Navier-Stokes systems with L∞ tensor coefficient under relaxed ellipticity condition. Discrete & Continuous Dynamical Systems, 2021  doi: 10.3934/dcds.2021042 [14] Tobias Breiten, Sergey Dolgov, Martin Stoll. Solving differential Riccati equations: A nonlinear space-time method using tensor trains. Numerical Algebra, Control & Optimization, 2021, 11 (3) : 407-429. doi: 10.3934/naco.2020034 [15] Manfred Einsiedler, Elon Lindenstrauss. On measures invariant under diagonalizable actions: the Rank-One case and the general Low-Entropy method. Journal of Modern Dynamics, 2008, 2 (1) : 83-128. doi: 10.3934/jmd.2008.2.83 [16] Hirofumi Notsu, Masato Kimura. Symmetry and positive definiteness of the tensor-valued spring constant derived from P1-FEM for the equations of linear elasticity. Networks & Heterogeneous Media, 2014, 9 (4) : 617-634. doi: 10.3934/nhm.2014.9.617 [17] Yumi Yahagi. Construction of unique mild solution and continuity of solution for the small initial data to 1-D Keller-Segel system. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021099 [18] Brian Ryals, Robert J. Sacker. Bifurcation in the almost periodic $2$D Ricker map. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021089 [19] Lin Yang, Yejuan Wang, Tomás Caraballo. Regularity of global attractors and exponential attractors for $2$D quasi-geostrophic equations with fractional dissipation. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021093 [20] Yangrong Li, Fengling Wang, Shuang Yang. Part-convergent cocycles and semi-convergent attractors of stochastic 2D-Ginzburg-Landau delay equations toward zero-memory. Discrete & Continuous Dynamical Systems - B, 2021, 26 (7) : 3643-3665. doi: 10.3934/dcdsb.2020250

2019 Impact Factor: 1.373