Efficient solution of large-scale matrix of acoustic wave equations in 3D frequency domain
Changcheng Li1 and Xiaofei Chen?2,1,3
1. Department of Earth and Space Sciences, Southern University of Science and Technology, Shenzhen 518055, China.
2. Shenzhen Key Labotory of Deep Offshore Oil and Gas Exploration Technology, Southern University of Science and Technology, Shenzhen 518055, China.
3. Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 511458, China.
Abstract In 3D frequency domain seismic forward and inversion calculation, the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation effi ciency. The frequency domain fi nite-diff erence forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation, and then the frequency domain wave fi eld value is obtained by solving the matrix equation. In this study, the predecessor’s optimized five-point method is extended to a 3D seven-point finite-difference scheme, and then a perfectly matched layer absorbing boundary condition (PML) is added to establish the corresponding matrix equation. In order to solve the complex matrix, we transform it to the equivalent real number domain to expand the solvable range of the matrix, and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods, and then use conjugate gradient algorithm to solve the problem. Previous studies have shown that in the conjugate gradient algorithm, the product of the matrix and the vector is the main factor that aff ects the calculation effi ciency. Therefore, this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm, thereby reducing the amount of calculation and storage.
Fund: This work is supported by the National Natural Science Foundation of China (Project U1901602&41790465), Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (GML2019ZD0203), Shenzhen Key Laboratory of Deep Off shore Oil and Gas Exploration Technology (Grant No.ZDSYS20190902093007855), Shenzhen Science and Technology Program (Grant No. KQTD20170810111725321), and the leading talents of Guangdong province program (Grant No.2016LJ06N652).
About author: LI Changcheng is currently a PhD student in Earth Science at the South University of Science and Technology, China. At present, its main work is focused on the research of software and hardware related to the detection of different scale medium structures in the interior of planets in seismology and Marsquake. E-mail: 1020168681@qq.com
Cite this article:
LI Chang-Cheng-1,CHEN Xiao-Fei-2,1 et al. Efficient solution of large-scale matrix of acoustic wave equations in 3D frequency domain[J]. APPLIED GEOPHYSICS, 2021, 18(3): 299-316.
2021, Vol. 18 
