地震层析成像中LSQR算法的快速求解.pdfVIP

第33卷 第6期 物探化探计算技术 2011年11月 文章编号：1o01一l749(2011)O6—0632—04 地震层析成像中LSQR算法的快速求解 张 东 ，乔友峰 ，姜麟舜 ，杨 艳 ，秦前清 ，徐 凌2 (1．武汉大学，湖北 武汉 430072；2．中国石油勘探开发研究院 物探所，j匕京 100082) 摘 要：针对传统LSQR反演算法计算时间长、占用内存大的不足，这里提出了一种LSQR的快 速算法，通过线性三元组存储稀疏矩阵，实现相应的矩阵运算，使 LSQR算法的运算速度及 内存 占用都有较大地改善。 关键词：LSQR；地震层析成像；反演；稀疏矩阵 中图分类号：P631．4 文献标识码 ：A 程化为系数矩阵为方阵的方程，然后利用 Lanczos 0 前言 方法，求解方程的最 ／J~--乘解。LSQR算法是一种 迭代算法，具有计算量小的优点，并且能很容易地 层析成像 的迭代算法主要有：反投影技术 利用矩阵的稀疏性简化计算，适合求解大型稀疏问 BPT、代数重建技术 ART、联合迭代重建技术 题。对于方程组Ax=b的求解，LSQR算法的迭代 SIRT、共轭梯度法CG、最小平方正交分解法 LSQR 过程及具体的推导过程见参考文献[4]。 等-l】。】。BPT算法计算简单、快速 ，但分辨率低， ART和 SIRT算法 速度 快，但有 时不 收敛。 2 LSQR算法的特点分析 LSQR-3算法是 目前常用的方法，与其它迭代方法 相比，其在解奇异或病态问题时显示出更快的收 2．1 反演方程的系数矩阵 敛性及更好的结果。但LSQR算法 占用计算机内 对于反演方程组Ax=b，其中A为系数矩阵， 存较大，计算成本高，尤其是应用于三维地震层析 它的元素由每一条射线在每个网格中的长度构成。 成像中时，对内存的需求以及计算成本将会更高。 b为走时参数，它的每个元素对应一条射线在初始 如果对LSQR算法的求解进行一定的改进，使算法 模型中的走时和在理论模型中走时的差。通过研 在内存 占用上降低，同时提高运算速度，那么， 究发现，由于每一条射线只能经过某一部份网格， LSQR算法将会更好地应用于三维地震层析成像。 因此系数矩阵A是稀疏的。 作者在本文针对 LSQR存在的不足，提出了一 表 1(见下页)给出了算例模型矩阵A的元素 种快速的求解方法，使用线性三元组存储稀疏矩阵 统计，通过观察表 1可以发现，当最大炮检距越来 来实现矩阵运算，由于只涉及到非零元素，所以占 越大时，非零元素的比例趋于减小的趋势，而且这 用存贮空间少，运算速度快，特别适于求解系数为 一 比例都在 l％以内。由此可见，系数矩阵是极其 大型稀疏矩阵的方程组。 稀疏的。 2．2 反演方程的求解 1 LSQR算法原理 LSQR算法虽然是一种适合于大型稀疏矩阵 方程求解的算法，但经过实验我们发现，现有的 LSQR算法的思路是首先把任意稀疏矩阵方 LSQR算法，还是存在存储空间占用大，运算效率 收稿 日期 ：2011—05—25 改回 日期 ：2011一O8—30 6期 张东等：地震层析成像中I_SQR算法的快速求解 633 不高的问题。例如 MATLAB中的LSQR函数 ]， 该函数的原型为LSQR(A，B，TOL，MAXIT，X0)，其 中A为系数矩阵，B为