三维旅行时场B样条插值射线追踪方法
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
传统的三维射线追踪算法由于受到计算区域网格离散化的限制,会导致追踪结果的射线路径出现不应有的弯折和相邻射线合并现象,进而影响追踪精度。本文根据地震初至波路径垂直于波前面的原理,提出基于B样条插值的旅行时场梯度反向射线追踪方法。该方法通过三次B样条插值得到空间连续的旅行时场,并计算旅行时场梯度;然后从检波点开始,沿着旅行时场梯度的负方向进行追踪,即可得到从震源到检波点的最小旅行时路径。由于通过样条插值能获得连续的旅行时场,故追踪步长不受网格尺寸的限制,从而克服计算区域离散化所带来的问题。应用具体算例对该方法的效果进行了验证,结果表明,该方法在计算精度和效率上都有明显改善。应用于地震初至波层析成像的正演计算时,能够使收敛误差更小。
Grid dissection of computational areas in conventional 3-D seismic ray tracing algorithms always lead to some zigzags and unexpected mergence in ray paths which may bring lower precision to calculation.To solve these problems,a new ray tracing scheme based on B-spline interpolation in backward step is presented in this paper.We calculate traveltime at any point in the area by cubic B-spline interpolation,and then get the field of traveltime gradient.According to the Fermat's principle,the ray path is always perpendicular to the wavefront which is the same of the negative traveltime gradient.First arrival ray path then can be traced from receiver point to source point along the negative gradient step by step.Since the traveltime field reconstructed by cubic B-spline interpolation is continuous and smooth,the tracing step-size has no limitation of the grid size,and the tracing error introduced by area discretization can be reduced obviously.The results of two experimental models revealed the improvement in precision and efficiency of the new algorithm.And less convergence error can also be reached when the new algorithm is applied to forward modeling for seismic traveltime tomography.
引文
[1]李宇,杨德义,邓辉等.初至旅行时层析反演近地表模型精度分析.石油地球物理勘探,2011,46(2):221~225 Li Yu,Yang Deyi,Deng Hui et al.Analysis on in-version accuracy of near-surface model inversed by first break traveltime tomography.OGP,2011,46(2):221~225
    [2]和锐,杨建思,张翼.地震层析成像方法综述.CT理论与应用研究,2007,16(1):35~48 He Rui,Yang Jiansi,Zhang Yi.A review on the technology of seismic tomography.Computerized Tomography Theory and Applications,2007,16(1):35~48
    [3]张钋,刘洪,李幼铭.射线追踪方法的发展现状.地球物理学进展,2000,15(1):36~45 Zhang Po,Liu Hong,Li Youming.The situation and progress of ray tracing method research.Progress in Geophysics,2000,15(1):36~45
    [4]Vidale J E.Finite-difference calculation of traveltimes.Bull Seism Soc Am,1988,78(6),2062~2076
    [5]Vidale J E.Finite-difference calculation of traveltimes in three dimensions.Geophysics,1990,55(5):521~526
    [6]李文杰,魏修成,宁俊瑞等.一种改进的利用程函方程计算旅行时的方法.石油地球物理勘探,2008,43(5):589~594 Li Wenjie,Wei Xiucheng,Ning Junrui et al.A method using improved eikonal equation to compute travel time.OGP,2008,43(5):589~594
    [7]单联瑜,刘连升.三维高阶快速步进波前重构旅行时计算精度分析.石油地球物理勘探,2007,42(4):413~417 Shan Lianyu,Liu Liansheng.Analysis on computa-tion precision of wavefront-reconstructing traveltime for3-D high-order rapid stepping algorithm.OGP,2007,42(4):413~417
    [8]马德堂,朱光明,范廷恩.二维TTI介质中初至波旅行时的搜索算法.石油地球物理勘探,2011,46(5):710~714 Ma Detang,Zhu Guangming,Fan Tingen.A search algorithm for first break traveltime in2DTTI medi-um.OGP,2011,46(5):710~714
    [9]Asakawa E,Kanawa T.Seismic ray tracing using lin-ear traveltime interpolation.Geophysical Prospec-ting,1993,41(1):99~111
    [10]Moser T J.Shortest path calculation of seismic rays.Geophysics,1991,56(1):59~67
    [11]卞爱飞,於文辉.三维最短路径法射线追踪及改进.天然气工业,2006,26(5):43~45 Bian Aifei,Yu Wenhui.3-D shortest path method of ray tracing and its improvement.Natural Gas Indus-try,2006,26(5):43~45
    [12]聂建新,杨慧珠.地震波旅行时二次/线性联合插值法.清华大学学报(自然科学版),2003,43(11):1495~1498 Nie Jianxin,Yang Huizhu.Quadratic/linear travel time interpolation of seismic ray-tracing.Journal of Tsinghua University(Science and Technology Edi-tion),2003,43(11):1495~1498
    [13]黄翼坚,朱光明,敦敏.直达波旅行时线性插值算法.石油地球物理勘探,2010,45(2):225~229 Huang Yijian,Zhu Guangming,Dun Min.Direct wave traveltime linear interpolation algorithm.OGP,2010,45(2):225~229
    [14]张赛民,周竹生,陈灵君等.对旅行时进行抛物型插值的地震射线追踪方法.地球物理学进展,2007,22(1):43~48 Zhang Saimin,Zhou Zhusheng,Chen Lingjun et al.Seismic ray-tracing method of applying parabolic in-terpolation to travel-time.Progress in Geophysics,2007,22(1):43~48
    [15]张东,谢宝莲,杨艳等.一种改进的线性旅行时插值射线追踪算法.地球物理学报,2009,52(1):200~205 Zhang Dong,Xie Baolian,Yang Yan et al.A ray tracing method based on improved linear traveltime interpolation.Chinese Journal of Geophysics,2009,52(1):200~205
    [16]张东,童海杰,杨艳等.基于波前扩展的线性旅行时插值射线追踪算法.武汉大学学报(理学版),2010,56(6):649~654 Zhang Dong,Tong Haijie,Yang Yan et al.Linear traveltime interpolation algorithm based on expanding wavefronts.Journal of Wuhan University(Natural Science Edition),2010,56(6):649~654
    [17]梅胜全,邓飞,钟本善等.基于改进的双线性旅行时插值的三维射线追踪.物探化探计算技术,2010,32(2):152~157 Mei Shengquan,Deng Fei,Zhong Benshan et al.The 3Dray tracing method base on the improved bilinear traveltime interpolation.Computing Techniques for Geophysical and Geochemical Exploration,2010,32(2):152~157
    [18]Zhang Jianzhong,Huang Yueqin,Song Linping et al.Fast and accurate3-D ray tracing using bilinear trav-eltime interpolation and the wave front group mar-ching.Geophysical Journal International,2011,184(3):1327~1340
    [19]张东,傅相如,杨艳等.基于LTI和网格界面剖分的三维地震射线追踪算法.地球物理学报,2009,52(9):2370~2376 Zhang Dong,Fu Xiangru,Yang Yan et al.3-D seis-mic ray tracing algorithm based on LTI and partition of grid interface.Chinese Journal of Geophysics,2009,52(9):2370~2376
    [20]刘根洪,何仁杰.微分几何与计算几何.四川成都:四川大学出版社,2000

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