Tensor tree decomposition as a rank‐reduction method for pre‐stack interpolation

Abstract

AbstractTensors have been proposed to represent pre‐stack seismic data, particularly for seismic data denoising and reconstruction. They naturally permit describing multi‐dimensional seismic signals that depend on time (or frequency) and source and receiver coordinates. A tensor representation aims to preserve the information embedded in the multi‐linear array in a reduced space. Such a representation is part of many algorithms for seismic data reconstruction via tensor completion methodologies. We investigate and apply one particular tensor tree representation to seismic data reconstruction. The Tensor Tree decomposition methodology permits decomposing a high‐order tensor into third‐order tensors. The technique relies on the truncated singular value decomposition, which branches the tensors into low‐dimensional tensors. As a benefit, the tensor tree allows us to reorganize the tensor into a matrix that demands the least singular values to have an optimal low‐rank approximation. We have developed an algorithm that uses the tensor tree for data reconstruction in an iterative optimization scheme and directly compared it to the parallel matrix factorization. At first, we demonstrated the proposed methodology results with five‐dimensional synthetic shot data and then moved forward with five‐dimensional field data, where we analysed it both pre and post‐stack. The tensor tree performs well in reconstructing both synthetic and field data with high fidelity, at the same level as the well‐established parallel matrix factorization.