阻尼LSQR算法在多波地震反演中的应用
|
摘要
多波VTI介质AVA储层参数联合反演的关键是解Jacobi矩阵。本文从Lanczos迭代出发,结合双对角化和QR分解,推导出了阻尼Lanczos算法流程。通过比较SVD算法和LSQR算法,显示了阻尼Lanczos算法在效率和精度方面的优越性。该算法应用于多波地震勘探资料反演中,获得的目标层储层参数:纵波速度、横波速度和密度,证明了该算法的有效性。
The key of the joint inversion of reservoir parameters for the Multi-wave VTI media AVA is the solution of the Jacobi matrix.This paper based on the Lanczos iteration,deduced Lanczos algorithm flow damping by combining double diagonalization and QR decomposition.Compared with SVD algorithm and LSQR algorithm,it was showed that the damped Lanczos algorithm was superior in its efficiency and accuracy.Then,we applied this algorithm to the Multi-wave seismic data inversion to prove the effectiveness of the algorithm.
The key of the joint inversion of reservoir parameters for the Multi-wave VTI media AVA is the solution of the Jacobi matrix.This paper based on the Lanczos iteration,deduced Lanczos algorithm flow damping by combining double diagonalization and QR decomposition.Compared with SVD algorithm and LSQR algorithm,it was showed that the damped Lanczos algorithm was superior in its efficiency and accuracy.Then,we applied this algorithm to the Multi-wave seismic data inversion to prove the effectiveness of the algorithm.
引文
[1]徐树芳.矩阵计算的理论与方法[M].北京:北京大学出版社,1995:319～328.
[2]张赛民.跨孔地震层析成像研究[D].长沙:中南大学,2003.
[3]裘渔洋.线性约束矩阵最小二乘问题:理论与算法[D].杭州:浙江大学,2007.
[4]裴顺平.中国大陆上地幔顶部体波速度层析成像[D].北京:中国地震局地球物理所,2002.
[2]张赛民.跨孔地震层析成像研究[D].长沙:中南大学,2003.
[3]裘渔洋.线性约束矩阵最小二乘问题:理论与算法[D].杭州:浙江大学,2007.
[4]裴顺平.中国大陆上地幔顶部体波速度层析成像[D].北京:中国地震局地球物理所,2002.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心