requency domain wave equation forward modeling using gaussian elimination with static pivoting

• Published in: Applied geophysics Vol. 8; no. 1; pp. 60 - 68
• Main Author: Song Jian-Yong Zheng Xiao-Dong Zhang Yan Xu Ji-Xiang Qin Zhen Song Xue-Juan
• Format: Journal Article
• Language: English
• Published: 2011
• Subjects: Marmousi模型; 分布式并行计算; 正演模拟; 波动方程; 计算效率; 静电消除
• ISSN: 1672-7975; 1993-0658

Frequency domain wave equation forward modeling is a problem of solving large scale linear sparse systems which is often subject to the limits of computational efficiency and memory storage. Conventional Gaussian elimination cannot resolve the parallel computation of huge data. Therefore, we use the Gaussian elimination with static pivoting （GESP） method for sparse matrix decomposition and multi-source finite-difference modeling. The GESP method does not only improve the computational efficiency but also benefit the distributed parallel computation of matrix decomposition within a single frequency point. We test the proposed method using the classic Marmousi model. Both the single-frequency wave field and time domain seismic section show that the proposed method improves the simulation accuracy and computational efficiency and saves and makes full use of memory. This method can lay the basis for waveform inversion.