最佳路径射线追踪算法研究
|
摘要
射线追踪方法在地震波理论研究、地震波反演、层析成像中有重要作用,其中最佳路径射线追踪法现在应用较广泛。最佳路径射线追踪法是根据Fermat原理及图论中的最佳路径理论,先计算出所有未知节点的最小旅行时,然后从接收点开始,反向追踪炮点到接收点的射线路径。建立了均匀介质、含有高速体及低速体的介质、层状介质和连续介质的二维模型,用最佳路径法对它们进行了射线追踪。发现该方法在网格间距较大时,追踪出的射线呈"之"字型,计算走时比实际走时偏大。但是该方法能计算出直达波、折射波和回折波等初至波,适用于各种复杂的介质模型。
Ray tracing method is of great importance in seismic wave theory research seismic inversion and tomography,migration processes.Shortest path ray tracing algorithm is more widly used lately.The shortest path ray tracing method is based on Fermat's principle and the shortest path graph theory.First calculate the travel time of the unknown point,then starting from the receiver,pursue the path from the exciton to the receiver in reverse.In this article we establish the two-dimensional model of equality,medium with high speed or low speed,samdwich and continuity medium,using the shortest path method tracing the ray.In this article we find when the grid parameters is big,the result will be zigzag,and the calculated time is bigger than the real time.This method can calculate the nonstop wave、refracted wave and inflexion wave,and can be applied in any complex model.
Ray tracing method is of great importance in seismic wave theory research seismic inversion and tomography,migration processes.Shortest path ray tracing algorithm is more widly used lately.The shortest path ray tracing method is based on Fermat's principle and the shortest path graph theory.First calculate the travel time of the unknown point,then starting from the receiver,pursue the path from the exciton to the receiver in reverse.In this article we establish the two-dimensional model of equality,medium with high speed or low speed,samdwich and continuity medium,using the shortest path method tracing the ray.In this article we find when the grid parameters is big,the result will be zigzag,and the calculated time is bigger than the real time.This method can calculate the nonstop wave、refracted wave and inflexion wave,and can be applied in any complex model.
引文
[1]张永刚.地震波场数值模拟方法[J].石油物探,2003,42(2):143-148.
[2]张钋,刘洪,李幼铭.射线追踪方法的发展现状[J].地球物理学进展,2000,15(1):36-45.
[3]Vidale J E.Finite-differencecal culation of taveltrimes[J].Bull.,Seis.Sos.Am,1988,78(6):2062-2076.
[4]Vidale J E.Finite-difference calculation of taveltimes inthree dimension[J].Geophysics,1990,55(5):521-526.
[5]Van Trier J,Symes W W.Upwind finite-difference calcu-lation of traveltimes[J].Geophysics,1991,56(6):812-821.
[6]黄联捷,李幼铭,吴如山.用于图象重建的波前射线追踪[J].地球物理学报,1992,35(2):223-233.
[7]Sava P,Fomel S.Huygens wavefront tracing:A robust al-ternative to raytracing[J].68 Ann.Internat.Mtg.,Soc.Expl.Geophys.,Expanded Abstracts,1998:1961-1964.
[8]Symes W W.Aslowness matching finite difference methodfor traveltimes beyond transmission caustics[J].68 Ann.Internat.Mtg.,Soc.Expl.Geophys.,Expanded Abstracts,1998:1945-1948.
[9]高亮,李幼铭,陈旭荣,等.地震辛几何算法初探[J].地球物理学报,2002,45(3):402-410.
[10]陈景波,秦孟兆.射线追踪、辛几何算法与波场数值模拟[J].计算物理,2001,18(6):481-487.
[11]秦孟兆,陈景波.Moslov渐进理论与辛几何算法[J].地球物理学报,2000,43(4):522-533.
[12]Moser TJ.Shortest path calculation of seismic rays[J].Geophysics,1991,56(1):59-67.
[13]张建中,陈世军,余大祥.最短路径射线追踪方法及其改进[J].地球物理学进展,2000,15(1):146-150.
[14]许琨,吴律,王妙月.改进Moser法射线追踪[J].地球物理学进展,1998,13(4):61-66.
[15]张美根,程冰洁,李小凡,等.一种最短路径射线追的快速算法[J].地球物理学报,2006,49(5):1467-174.
[16]李志辉,刘争平.最短路径法射线追踪的MATLAB实现[J].工程地质计算机应用,2004,4:16-21.
[17]苏晓生.掌握MATLAB6.0及其工程应用[M].北京:科学出版社,2002.
[2]张钋,刘洪,李幼铭.射线追踪方法的发展现状[J].地球物理学进展,2000,15(1):36-45.
[3]Vidale J E.Finite-differencecal culation of taveltrimes[J].Bull.,Seis.Sos.Am,1988,78(6):2062-2076.
[4]Vidale J E.Finite-difference calculation of taveltimes inthree dimension[J].Geophysics,1990,55(5):521-526.
[5]Van Trier J,Symes W W.Upwind finite-difference calcu-lation of traveltimes[J].Geophysics,1991,56(6):812-821.
[6]黄联捷,李幼铭,吴如山.用于图象重建的波前射线追踪[J].地球物理学报,1992,35(2):223-233.
[7]Sava P,Fomel S.Huygens wavefront tracing:A robust al-ternative to raytracing[J].68 Ann.Internat.Mtg.,Soc.Expl.Geophys.,Expanded Abstracts,1998:1961-1964.
[8]Symes W W.Aslowness matching finite difference methodfor traveltimes beyond transmission caustics[J].68 Ann.Internat.Mtg.,Soc.Expl.Geophys.,Expanded Abstracts,1998:1945-1948.
[9]高亮,李幼铭,陈旭荣,等.地震辛几何算法初探[J].地球物理学报,2002,45(3):402-410.
[10]陈景波,秦孟兆.射线追踪、辛几何算法与波场数值模拟[J].计算物理,2001,18(6):481-487.
[11]秦孟兆,陈景波.Moslov渐进理论与辛几何算法[J].地球物理学报,2000,43(4):522-533.
[12]Moser TJ.Shortest path calculation of seismic rays[J].Geophysics,1991,56(1):59-67.
[13]张建中,陈世军,余大祥.最短路径射线追踪方法及其改进[J].地球物理学进展,2000,15(1):146-150.
[14]许琨,吴律,王妙月.改进Moser法射线追踪[J].地球物理学进展,1998,13(4):61-66.
[15]张美根,程冰洁,李小凡,等.一种最短路径射线追的快速算法[J].地球物理学报,2006,49(5):1467-174.
[16]李志辉,刘争平.最短路径法射线追踪的MATLAB实现[J].工程地质计算机应用,2004,4:16-21.
[17]苏晓生.掌握MATLAB6.0及其工程应用[M].北京:科学出版社,2002.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心