反射地震技术在成矿地质背景与深部矿产勘查中的应用_现状与前景.docVIP

反射地震走时层析成像中的大型稀疏矩阵压缩存 反射地震走时层析成像中的大型 稀疏矩阵压缩存储和求解 成 谷1 , 张宝金2 ( 1 . 中山大学地球科学系 ,广州 510275 ; 2 . 广州海洋地质调查局 ,广州 510760) 摘 要 反射地震走时层析成像是一种精度较高的速度求取方法 ,最终可归结为线性方程组的求解 . 方程组具有很 大的维数 ,常规解法需要很大的存储量和计算量 . 本文考虑到当投影函数取为走时残差 ,图象函数取为慢度残差时 , 灵敏度矩阵中的元素表示射线经过网格的长度的特殊物理意义 ,采用行索引的压缩存储方式 ,在射线追踪正演模拟 过程中直接压缩存储灵敏度矩阵 ,在层析反演过程中利用压缩后的矩阵进行求解 ,大大降低了存储量和计算量 . 关键词 反射地震 ,走时层析 ,稀疏矩阵 ,压缩存储 ,压缩求解 中图分类号 P631 文献标识码 A 文章编号 100422903 (2008) 0320674207 Compression storage and sol ution of large sparse matrix in traveltime tomogra phy of ref lect ion seismic data C H EN G Gu1 , Z H A N G Bao2ji n2 ( 1 . De p a rt ment o f Ea rt h S ciences , S u n Y at2S en U ni ve rsi t y , Gu an g z hou 510275 , Chi na ; 2 . Guan g z hou M a ri ne Geol o g ical S u rve y , Gu an g z hou 510760 , Chi na) Abstract Traveltime to mo grap hy of ref lectio n sei smic data i s o ne of t ho se met ho ds w hich have higher accuracy to o b2 t ain velocit y a nd finally co me s do w n to t he sol utio n of t he lar ge spa r se equatio n system. But ver y la r ge sto rage amo unt a nd calculatio n a mo unt a re demanded w hen solving t he equatio n mat rix i n regular way , due to t he high dimensio n a nd t he iterative sol utio n st yle , w hich st ro ngly rest rict s t he p ractical applicatio n of to mo grap hy met ho d. Thus t he re2 sea rch o n t he met ho ds reducing calculatio n a mo unt a nd sto rage a mo unt i s ver y impo rta nt a nd significative. In t hi s pa2 per , t he met ho ds of co mp re ssio n sto rage a nd sol utio n of lar ge spa r se mat rix a re used in to mo grap hy i nver sio n. The line2i ndexed co mp ressio n sto rege i s cho sen , due to t hat t he element of sensitivit y mat rix ha s sp ecial p hysical mea ning t hat t he val ue of t he element rep resent s t he lengt h of ray in each velocit y cell , w hen co nsidering t raveltime re sidual a s p rojectio n f unctio n a nd slo w ne ss residual a s image f unctio n in Rado n t ra nsfo r mtio n. Thus ea