反射地震走时层析成像中的大型稀疏矩阵压缩存储和求解

详细信息   

• 责任者：成谷^c中山大学^d1975^e博士 ; 张宝金^c广州海洋地质调查局
• 刊名：地球物理学进展
• 出版年：2008
• 卷期：23^b3总89
• 关键词：反射波法 ; 层析技术20^a地震勘探 ; 走时残差 ; 成象 ; 射线追踪 ; 正演问题 ; 反演问题 ; 数据处理
• 页码：674～680
• 索书号：P260.6/208-1
• 图表：^c表4参18插图
• 分类号：2100
• 文摘期号：2008/11

摘要

反射地震走时层析成像是一种精度较高的速度求取方法，最终可归结为线性方程组的求解。方程组具有很大的维数，常规解法需要很大的存储量和计算量。该文考虑到当投影函数取为走时残差，图象函数取为慢度残差时，灵敏度矩阵中的元素表示射线经过网格的长度的特殊物理意义，采用行索引的压缩存储方式，在射线追踪正演模拟过程中直接压缩存储灵敏度矩阵，在层析反演过程中利用压缩后的矩阵进行求解，大大降低了存储量和计算量。