地震走时层析成像方法中大型稀疏矩阵方程解算方法比较.docVIP

地震走时层析成像方法中大型稀疏矩阵方程解算方法比较 X 杨海波 (成都理工大学信息工程学院,四川成都　610059 　　摘　要:在地震走时层析成像方法中,对大型稀疏矩阵方程解算方的稳定性和计算效率的研究是一个重要问题。本文列举了五种用于解决此类问题的常用算法:SVD 、SIRT 、Jacobi-ICCG 、SSOR-ICCD 和LSQR,并在求解一般问题和病态问题的试验中,分析各方法的误差、计算效率和稳定性;进而使用各算法对同一地质模型进行速度成像,并对结果进行比较,给出针对各种方法的客观评价。 关键词:地震走时层析成像;稳定性;计算效率;误差分析 　　中图分类号:P 631.4+43　　文献标识码:A 　　文章编号:1006—7981(201010—0051—021　问题的提出 地震走时层析成像是根据地震波走时观测数据求解地球内部结构的方法。在这项研究中,根据对问题的地质了解假定一个初始模型速度模型,用射线追踪方法得到射线的理论走时,解决模型正问题,再借助于模型参数扰动,构成初始模型的理论值和观测值之间的线性方程,反演速度结构,修改初始模型重复上述过程,直到获得满意的结果。一般而言,地震层析成像反演问题可以归结为解不适定方程组: Ax=B (1通过构造正规方程:A T Ax =A T b (2解出超定方程的最小二乘解,就可以得到这个断面中各个单元的慢度,从而推测出断面中不同物质的