基于三角域快速行进法的地震波走时计算
|
摘要
地震波走时计算是地震资料解释处理技术的重要组成部分,本文根据复杂地层构造中速度场分布的特点,设计了一种采用快速行进法基于三角网格的走时计算方法,针对计算效率优化和快速行进法在三角域上的计算格式进行了重点的研究,并根据地层限定条件对速度场进行网格剖分,在三角网格上用快速行进算法计算各点走时。与基于矩形网格的差分方法相比,该方法不需要对速度场边界进行任何平滑,无须通过细分网格来提高计算精度;可根据不同地质构造的复杂度进行变网格大小的剖分,网格剖分数目相对较少。最后通过计算实例进行了验证。
Travel-time computation is an important technique of seismic data process.According to the traits of complex velocity model,we present a travel-time computation using triangulated fast marching method.First,we get a triangulated mesh using 2d constrained triangulation.Then the triangulated fast marching method is introduced to compute travel-time.Compared to other finite difference wavefront computation method based on rectangle cell parameteration,this method based on triangle cell parameteration doesn't need any flatness on the interface of velocity model to improve the precision.It presents the complex velocity model with less cells through control the size of triangle on different part.The results of computation on Marmousi model demonstrate correctness of the proposed methods.
Travel-time computation is an important technique of seismic data process.According to the traits of complex velocity model,we present a travel-time computation using triangulated fast marching method.First,we get a triangulated mesh using 2d constrained triangulation.Then the triangulated fast marching method is introduced to compute travel-time.Compared to other finite difference wavefront computation method based on rectangle cell parameteration,this method based on triangle cell parameteration doesn't need any flatness on the interface of velocity model to improve the precision.It presents the complex velocity model with less cells through control the size of triangle on different part.The results of computation on Marmousi model demonstrate correctness of the proposed methods.
引文
[1]渥·伊尔马滋.地震资料分析[M].北京:石油工业出版社,2006.
[2]Vidale J.FInite-differendce calcaulation of traveltimes[J].Bull Seism Soc Am,1988,78(6):2062-2076
[3]Trier J Van,Symes W W.Upwind finite-difference calcu-lation of traveltimes[J].Geophysics,1991,56(6):812-821
[4]Asakawa E,Kawanaka T.Seismic ray tracing using linear traveltime interpolation[J].Geophysical Prospecting,1993,41(1):99-111
[5]王华忠,方正茂,匡斌等.任意复杂介质中的动态规划法地震波三维走时计算[J].地球物理学报,2001.44(增刊):179-189.
[6]高亮,李幼铭,陈旭容等.地震射线辛几何算法初探[J].地球物理学报,2000,43(3):402-410.
[7]J.A.Sethian.Level Set Methods and Fast Marching Methods[M].2nd edition.USA:Cambridge University Press,1999.
[8]Sethian J A,Popovici A M.3-D traveltime computa-tion using the fast marching method[J].Geophysics,1999,64(2):516-523.
[9]徐永安,谭建荣,杨钦,陈其明.二维任意域约束Delaunay三角化的实现[J].工程图学学报.1999,54(1):51-55.
[10]杨钦,刘瑞刚,孟宪海,张俊安.二维复杂限定条件的Delaunay三角化算法[J],计算机辅助设计与图形学学报,2007,19(2):145-150.
[11]成谷.地震反射走时层析理论与应用研究[D].上海:同济大学海洋与地球科学学院,2004.
[2]Vidale J.FInite-differendce calcaulation of traveltimes[J].Bull Seism Soc Am,1988,78(6):2062-2076
[3]Trier J Van,Symes W W.Upwind finite-difference calcu-lation of traveltimes[J].Geophysics,1991,56(6):812-821
[4]Asakawa E,Kawanaka T.Seismic ray tracing using linear traveltime interpolation[J].Geophysical Prospecting,1993,41(1):99-111
[5]王华忠,方正茂,匡斌等.任意复杂介质中的动态规划法地震波三维走时计算[J].地球物理学报,2001.44(增刊):179-189.
[6]高亮,李幼铭,陈旭容等.地震射线辛几何算法初探[J].地球物理学报,2000,43(3):402-410.
[7]J.A.Sethian.Level Set Methods and Fast Marching Methods[M].2nd edition.USA:Cambridge University Press,1999.
[8]Sethian J A,Popovici A M.3-D traveltime computa-tion using the fast marching method[J].Geophysics,1999,64(2):516-523.
[9]徐永安,谭建荣,杨钦,陈其明.二维任意域约束Delaunay三角化的实现[J].工程图学学报.1999,54(1):51-55.
[10]杨钦,刘瑞刚,孟宪海,张俊安.二维复杂限定条件的Delaunay三角化算法[J],计算机辅助设计与图形学学报,2007,19(2):145-150.
[11]成谷.地震反射走时层析理论与应用研究[D].上海:同济大学海洋与地球科学学院,2004.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心