Application of the preconditioned conjugate gradient method to reconstruction of seismic data

作者：Huo Zhi-Zhou^*; Xiong Deng; Zhang Jian-Feng

来源：Acta Geophysica Sinica, 2013, 56(4): 1321-1330.

DOI：10.6038/cjg20130426

• Precondition
• Conjugate Gradient Method
• Seismic data reconstruction
• Fourier method

Seismic data reconstruction based on the least-squares Fourier method is ultimately transformed into solving a linear equation. The coefficient matrix is Toeplitz matrix. The conjugate gradient method can be used to solve the linear equations. Pathological extent of matrix affects iterations of the conjugate gradient method. The more irregular sampled seismic data is, the more pathological the matrix is, then it is more difficult to get convergence and reasonable results. We study different construction methods of preconditions based on Toeplitz matrix and the effects of convergence of the conjugate gradient method. Through the use of preconditions, we can speed up the iterative speed of the conjugate gradient method, improve the convergence of conjugate gradient method and the efficiency of computation. Numerical examples and real seismic data reconstruction experiment show that the preconditioned conjugate gradient method has greatly improved efficiency of calculation.

• 出版日期2013-4
• 单位中国科学院大学; 中国科学院地质与地球物理研究所兰州油气资源研究中心; 中国科学院地质与地球物理研究所