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.

A nonuniform mesh method in the Floquet parameter domain for wave scattering by periodic surfaces

Sound Field Separating on Arbitrary Surfaces Enclosing a Sound Scatterer Based on Combined Integral Equations.

An adaptive finite element method for high-frequency scattering problems with smoothly varying coefficients

Wave scattering with time-periodic coefficients: Energy estimates and harmonic formulations

A Fast Hybrid Galerkin Method for High-Frequency Acoustic Scattering

A Bloch transform based numerical method for the rough surface scattering problems

A Spectral Boundary Element Method for Scattering Problems

On Reducing High Computational Cost with Periodic Wavelets in Solving Two-Dimensional Acoustic Radiation and Scattering

Two New Schemes for Numerical Modeling of Acoustic Scattering

Time-harmonic Acoustic Scattering from a Non-Locally Perturbed Trapezoidal Surface

A spectrally accurate direct solution technique for frequency-domain scattering problems with variable media

Trapped acoustic waves and raindrops: high-order accurate integral equation method for localized excitation of a periodic staircase

An Adaptive Finite Element Method for the Wave Scattering with Transparent Boundary Condition

Acoustic wave diffraction by a quadrant of sound-soft scatterers

A boundary integral formulation for the prediction of acoustic scattering from periodic structures

A fast direct solver for quasi-periodic scattering problems

An FFT-based algorithm for efficient computation of Green's functions for the Helmholtz and Maxwell's equations in periodic domains

Simulating time-harmonic acoustic wave effects induced by periodic holes/inclusions on surfaces

A novel integral equation for scattering by locally rough surfaces and application to the inverse problem

A parallel computation method for scattered seismic waves using Nystrom discretization and FFT fast convolution