Optimizing Finite-Difference Operator in Seismic Wave Numerical Modeling

The finite-difference method is widely used in seismic wave numerical simulation, imaging, and waveform inversion. In the finite-difference method, the finite difference operator is used to replace the differential operator approximately, which can be obtained by truncating the spatial convolution s...

Full description

Bibliographic Details
Main Authors: Hui Li, Yuan Fang, Zhiguo Huang, Mengyao Zhang, Qing Wei
Format: Article
Language:English
Published: MDPI AG 2022-04-01
Series:Algorithms
Subjects:
Online Access:https://www.mdpi.com/1999-4893/15/4/132
Description
Summary:The finite-difference method is widely used in seismic wave numerical simulation, imaging, and waveform inversion. In the finite-difference method, the finite difference operator is used to replace the differential operator approximately, which can be obtained by truncating the spatial convolution series. The properties of the truncated window function, such as the main and side lobes of the window function’s amplitude response, determine the accuracy of finite-difference, which subsequently affects the seismic imaging and inversion results significantly. Although numerical dispersion is inevitable in this process, it can be suppressed more effectively by using higher precision finite-difference operators. In this paper, we use the krill herd algorithm, in contrast with the standard PSO and CDPSO (a variant of PSO), to optimize the finite-difference operator. Numerical simulation results verify that the krill herd algorithm has good performance in improving the precision of the differential operator.
ISSN:1999-4893