Recursive integral time extrapolation of elastic waves using low-rank symbol approximation

Next: Energy-norm imaging condition Up: Theory Previous: Heterogeneous media

## Low-rank approximation

So far, we have laid out our basic theory of recursive integral time extrapolation of elastic waves. In mildly heterogeneous media, the Christoffel matrix is SPD. In strongly heterogeneous media, the Christoffel matrix becomes complex-valued and non-Hermitian. However, in both cases, the eigenvalues and eigenvectors of the Christoffel matrix become dependent on both spatial location and propagation direction, in other words, they are functions of both space and wavenumber . If these operators are implemented straightforwardly, one is faced with the daunting task of computing and storing the complete eigenvalue decomposition of the Christoffel matrix using all the combinations of and , leading to computational and memory complexity, where refers to the total number of mesh points in D. To perform wave extrapolation in the form of integral operators, one would have to multiply matrices with vectors in dimension of , leading to a computational complexity of . This is simply infeasible for practical applications.

In this work, to efficiently apply the derived Fourier Integral Operators (FIOs), we proposed to apply the low-rank decomposition (Fomel et al., 2013) on the mixed-domain wave extrapolation matrices. Take the wave extrapolation operator in equation 27, , as an example. We propose to apply low-rank approximation on each individual element of its expansion. For instance, the element, which operates on the x-component of the input vector wavefield and outputs to the x-component of the output vector wavefield, can be approximated as (Fomel et al., 2013):

 (29)

where and are sampled representative columns and rows from the original matrix , and are the numerical ranks of matrix , and the matrix is obtained from minimizing

 (30)

Similarly, and can be approximated as:
 (31) (32)

The computation of then becomes:
 (33)

The computation of and components can be carried out in a similar fashion. The computational cost of applying each FIO reduces to a complexity of , where is the complexity of one forward or inverse Fast Fourier Transform (FFT), and is the numerical rank of the low-rank approximation, which is for homogeneous media and for heterogeneous media.

 Recursive integral time extrapolation of elastic waves using low-rank symbol approximation

Next: Energy-norm imaging condition Up: Theory Previous: Heterogeneous media

2018-11-16