A parallel computation method for scattered seismic waves using Nystrom discretization and FFT fast convolution
Xu YangYang,Sun JianGuo,Shang YaoDa
DOI: https://doi.org/10.6038/cjg2021O0391
2021-01-01
Abstract:The storage of coefficient matrix and solution of linear equations are the main factors of limiting the application of the Lippermann-Schwinger (L-S) integral equation to solving the seismic scattering problem. Because the integral part of the L-S equation is a spatial convolution, which can lead to a full rank matrix after discretization, thus forming a large or very large algebraic equation set. Therefore, when L-S equation is used to solve the seismic wave scattering problem, the scattering series method is generally used instead of the numerical method. However, the computational accuracy and convergence of the scattering series method strongly depend on the intensity of velocity perturbation, and one possible way to overcome such dependence is to improve the existing numerical methods or to establish new numerical solutions. Under the guidance of this idea, the L-S equation is first rewritten to obtain an integral equation equivalent to the original L-S equation (equivalent L-S equation). Then, the equivalent L-S equation is normalized point by point, and is discretized by the Nystrom method. The spatial convolution is calculated by Fast Fourier Transform(FFT). The coefficient matrix generated by the normalized equivalent L-S equation discretized by the Nystrom method is a Toeplitz matrix, and its Toeplitz property can be used to reduce the storage space. FFT is adopted to transform the convolution into a product, and the integral kernel only needs to be calculated once. Furthermore, in order to save computing time, MPI+OpenMP parallel mode combining process level and thread level are designed. Numerical experiments show that compared with the traditional integral equation numerical algorithm, the calculation scheme based on the equivalent L-S equation, Nystrom discretization and FFT fast convolution can greatly reduce the storage and CPU time required, and improve the calculation efficiency with sufficient accuracy.