预条件共轭梯度法在地震数据重建方法中的应用 Application of the preconditioned conjugate gradient method to reconstruction of seismic data.pdf

第56卷第4期 地 球 物 理 学 报 V01．56，No．4 CHlNESE 0FGEOPHYSICS Apr．，2013 2013年4月 JOURNAL cj920130426． HuoZ ofthe co methodtoreconstructionofseismicdata． J F． nj gradient Z，XiongD，ZhangApplicationpreconditionedugate C^i以PsP-厂．GPo户^ys．(inChinese)，2013，56(4)：1321一1330，doi：10．6038／cj920130426． 预条件共轭梯度法在地震数据重建方法中的应用 霍志周1’2，熊 登3，张剑锋1 1中国科学院地质与地球物理研究所，北京 100029 2中国科学院大学，北京100049 3中国石油东方地球物理公司物探技术研究中心，河北涿州 072751 摘 要 基于最小平方的Fourier地震数据重建方法最终转化为求解一个线性方程组，其系数矩阵是Toeplitz矩 阵，可以用共轭梯度法求解该线性方程组．共轭梯度法的迭代次数受系数矩阵病态程度的影响，地震数据的非规则 采样程度越高，所形成的系数矩阵病态程度越高，就越难收敛和得到合理的计算结果．本文研究了基于Toeplitz矩 阵的不同预条件的构造方法，以及对共轭梯度法收敛性的影响．通过预条件的使用，加快了共轭梯度法的迭代速 度，改进了共轭梯度算法的收敛性，提高了计算的效率．数值算例和实际地震数据重建试验证明J，预条件共轭梯度 法对计算效率有很大的提高． 关键词 预条件，共轭梯度法，地震数据重建，Fourier方法 doi：10．6038／cj920130426中图分类号P631 收稿日期2012一06-06，2012一09_28收修定稿 ofthe method ApplicationpreconditionedcoIlj ugategradient toreconstructionofseismicdata HU0 Zhi—Zhoul”，XIONG Den93，ZHANGJian—Fen91 lJn盯n“￡P 100029，C^inn o，(：kozogyⅡndGPo户^ys掂5，(冼inPsPAfddPmyo，SciFncPs，BP玎i”g 2 100049，(1^inⅡ U”i御r5i￡yo，(1^2nPjPAf日dPmyo，SciPcPs。BP0ing 3 RPsP日rf^CPn￡rP，BGP Z^“02^o“072751，Chi月Ⅱ GPo声^y5ic口ZTPf^niq“e hf．o，CⅣPC．HP6Pi AbstractSelsmlc datareconstructionbasedonthe Fouriermethodis least—squares ultimately