反射地震走时层析成像中的大型稀疏矩阵压缩存储和求解
Compression storage and solution of large sparse matrix in traveltime tomography of reflection seismic data作者机构:中山大学地球科学系广州510275 广州海洋地质调查局广州510760
出 版 物:《地球物理学进展》 (Progress in Geophysics)
年 卷 期:2008年第23卷第3期
页 面:674-680页
核心收录:
学科分类:081801[工学-矿产普查与勘探] 08[工学] 0818[工学-地质资源与地质工程]
基 金:广东省博士科研启动基金项目(32110-4203205)资助
摘 要:反射地震走时层析成像是一种精度较高的速度求取方法,最终可归结为线性方程组的求解.方程组具有很大的维数,常规解法需要很大的存储量和计算量.本文考虑到当投影函数取为走时残差,图象函数取为慢度残差时,灵敏度矩阵中的元素表示射线经过网格的长度的特殊物理意义,采用行索引的压缩存储方式,在射线追踪正演模拟过程中直接压缩存储灵敏度矩阵,在层析反演过程中利用压缩后的矩阵进行求解,大大降低了存储量和计算量.