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

Abstract

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 effi ciency 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 fi nite-difference modeling. The GESP method does not only improve the computational effi ciency but also benefi t 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 effi ciency and saves and makes full use of memory. This method can lay the basis for waveform inversion.