反射地震走时层析成像中的大型稀疏矩阵压缩存储和求解
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
反射地震走时层析成像是一种精度较高的速度求取方法,最终可归结为线性方程组的求解.方程组具有很大的维数,常规解法需要很大的存储量和计算量.本文考虑到当投影函数取为走时残差,图象函数取为慢度残差时,灵敏度矩阵中的元素表示射线经过网格的长度的特殊物理意义,采用行索引的压缩存储方式,在射线追踪正演模拟过程中直接压缩存储灵敏度矩阵,在层析反演过程中利用压缩后的矩阵进行求解,大大降低了存储量和计算量.
Traveltime tomography of reflection seismic data is one of those methods which have higher accuracy to obtain velocity and finally comes down to the solution of the large sparse equation system.But very large storage amount and calculation amount are demanded when solving the equation matrix in regular way,due to the high dimension and the iterative solution style,which strongly restricts the practical application of tomography method.Thus the research on the methods reducing calculation amount and storage amount is very important and significative.In this paper,the methods of compression storage and solution of large sparse matrix are used in tomography inversion.The line-indexed compression storege is chosen,due to that the element of sensitivity matrix has special physical meaning that the value of the element represents the length of ray in each velocity cell,when considering traveltime residual as projection function and slowness residual as image function in Radon transformtion.Thus each line in the sensitivity matrix can be directly stored in a compression style in ray-tracing modeling,and the large amout of storage can be retrenched since a ray only passes very a few cells.Furthermore the tomography inversion can be easily implemented with compression-stored sensitivity matrix because only the multiplication operation between matrix(or matrix transpose) and vector is demanded in solution of tomography equation system,and the operation can be easily achieved.And research shows that more stroage amount and calculation amount can be retrenched if the sensitivity matrix has larger sparseness.
引文
[1]Wang Y H,Pratt R G.Sensitivities of seismic traveltimes andamplitudes in reflection tomography[J].Geophys.J.Int.,1997,131:618~642.
    [2]Murphy G E,Gray S H.Munual seismic reflection tomo-graphy[J].Geophysics,1999,64(5):1546~1552.
    [3]Kosloff D,Sudman Y.Uncertainty in determining interval ve-locities from surface reflection seismics[J].Geophysics,2002,67(3):952~963.
    [4]Williamson P R.A guide to the limits of resolution imposed byscattering in ray tomography[J].Geophysis,1991,56(2):202~207.
    [5]成谷,马在田,耿建华,曹景忠.反射地震走时层析成像发展回顾[J].勘探地球物理进展,2002,25(3):6~12.Cheng G,Ma Z T,Geng J H,Cao J Z.A Review on thegrowth of seismic tomography[J].Progress in Exploration Ge-ophysics,2002,25(3):6~12.
    [6]马在田,曹景忠,王加林,等编著.计算地球物理学概论[M].上海:同济大学出版社,1997.1-10,254~284,294~308.Ma Z T,Cao J Z,Wang J L,et al.Generality of Calculationgeophysics[M].Shanghai:Press of tongji university,1997,1-10,254~284,294~308.
    [7]杨文采,李幼铭编著.应用反射地震走时层析成像[M].北京:地质出版社,1993.142~160.Yang W Ci,Li Y M,et al.Application of reflection seismictraveltime tomography[M].Beijing:Press of Geology,1993,142~160.
    [8]Tieman H J.Investigating the velocity-depth ambiguity of re-flection traveltimes[J].Geophysics,1994,59(11):1763~1773.
    [9]Rathor B S.Short note:Velocity-depth ambiguity in the dip-ping reflector case[J].Geophysica,1997,62(5):1583~1585.
    [10]Meng Z,Bleistein N.On velocity/depth ambiguity in 3-D mi-gration velocity analysis[J].Geophysics,2001,66(1):256~260.
    [11]Lines L.Short note:Ambiguity in analysis of velocity anddepth[J].Geophysics,1993,58(4):596~597.
    [12]Bishop T N,Bube K P,Cutler R T,et al.Tomographic de-termation of velocity and depth in laterally varying media[J].Geophysis,1985,50(6):903~923.
    [13]Stork C.Singular value decomposition of the velocity-reflec-tor depth tradeoff,part 2:High-resolution analysis of a ge-neric model[J].Geophysics,1992,57(7):933~943.
    [14]Langan R T,Lerche I,Cutler R T.Tracing of rays throughheterogeneous media:An accurate and efficient procedure[J].Geophysics,1985,50(9):1456~1465.
    [15]成谷,张宝金.三角网格参数化在反射地震走时层析成像中的应用[J].中山大学学报(自然科学版),2006,45(5):128~132.Cheng G,Zhang B J.Application of Triangle Cell Parameter-ization in Travel-time Tomography of Reflection Seismic Data[J].ACTA Scientiarum Naturalium Universitatis Sunyatseni,2006,45(5):128~132.
    [16]成谷,马在田,张宝金.反射地震走时层析成像中存在的主要问题及应对策略[J].地球物理学进展,2003,18(3):512~518.Cheng G,Ma Z T,Zhang B J.Primary problems and accordingstrategies in seismic tomography[J].Progress in Geophysics,2003,18(3):512~518.
    [17]成谷,张宝金,马在田.浅谈反射地震走时层析中的正则化[J].地球物理学进展,2007,22(6):1715~1721.Cheng G,Zhang B J,Ma Z T.Regularization in traveltimetomography for reflection seismic data[J].Progress in Geo-physics,2007,22(6):1715~1721.
    [18]William H P,Saul A T,William T V,Brian P F著,傅组芸,赵梅娜,丁岩石等译,C数值算法(第二版)[M].北京:电子工业出版社,2004.50~65.William H P,Saul A T,William T V,Brian P F,Fu Z Y,Zhao M N,Ding Y S,et al.C Numerical Algorithm[M].Bei-jing:Press of Electronics Industry,2004,50~65.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心