【技术实现步骤摘要】
一种混合二维地震波走时计算方法
本专利技术涉及地震波走时计算领域,特别是一种混合二维地震波走时计算方法。
技术介绍
《工程地球物理学报》2009年第3期公开了张双杰等“快速推进法计算精度分析及改进”,介绍了两种提高快速推进法计算精度的方法,一是采用高阶差分格式进行计算;二是对震源点附近网格节点进行局部精细处理,在震源附近采用高阶差分格式进行计算,在剩余区域采用低阶差分格式进行计算。并且通过这两种方法对均匀介质模型进行地震波走时计算,实验结果得到了比较好的效果。《世界地质》2016年第3期公开了王乾龙等“基于完全三叉树的快速推进法地震波走时计算”,介绍了一种改进的快速推进地震波走时计算方法,将完全三叉树排序方法引入到地震波走时计算中,减少了窄带扩展中寻找最小走时值的时间,提高了整个算法的计算效率。并且通过基于完全三叉树的快速推进地震波走时计算方法对层状模型、Marmousi模型、Sigsbee2a模型进行了计算,实验结果得到了比较好的效果。通过以上例子可以看出,现有快速推进地震波走时计算方法在一定程度上能够提升计算精度,但实现过程复杂,提升的计算精度也有限。
技术实现思路
本专利技术所要解决的技术问题是提供一种混合二维地震波走时计算方法,通过灵活运用窄带技术将波前构建地震波走时计算方法和快速推进地震波走时计算方法无缝连接起来,即通过在震源附近小范围内使用计算精度较高的波前构建法,在剩余区域使用快速推进法计算走时,在改善原有快速推进法计算精度的同时,仍保留了其高效的特点。为解决上述技术问题,本专利技术采用的技术方案是:一种混合二维地震波走时计算方法,包括以下步骤: ...
【技术保护点】
一种混合二维地震波走时计算方法,其特征在于,包括以下步骤:步骤1:读入相关参数文件、速度模型,其中,所述参数文件包含速度模型的网格点数、网格间距和震源位置;步骤2:从炮点沿着不同方向发射射线并计算中心射线信息;步骤3:使用波前构建法计算射线覆盖范围内网格节点走时;步骤4:对所有网格点进行分类,将它们分为接受点、窄带点或远离点,即如果一个网格点的地震波走时已经计算过,并且它周围点的走时全部计算过,则这个点为接受点;如果一个网格点的地震波走时计算过,而它周围点的走时并非全部计算过,则这个点为窄带点;如果一个点的走时没有被计算过,则这个点为远离点;步骤5:将所有窄带点移入窄带内,构建初始窄带;步骤6:扩展窄带,直至窄带为空;在这一过程中,网格点走时是通过迎风差分法求解二维程函方程获得的,所述二维程函方程为:|▽τ|=s其中,τ为地震波走时,s为慢度,▽为梯度符号,上述公式中梯度项的迎风差分表达式为:
【技术特征摘要】
1.一种混合二维地震波走时计算方法,其特征在于,包括以下步骤:步骤1:读入相关参数文件、速度模型,其中,所述参数文件包含速度模型的网格点数、网格间距和震源位置;步骤2:从炮点沿着不同方向发射射线并计算中心射线信息;步骤3:使用波前构建法计算射线覆盖范围内网格节点走时;步骤4:对所有网格点进行分类,将它们分为接受点、窄带点或远离点,即如果一个网格点的地震波走时已经计算过,并且它周围点的走时全部计算过,则这个点为接受点;如果一个网格点的地震波走时计算过,而它周围点的走时并非全部计算过,则这个点为窄带点;如果一个点的走时没有被计算过,则这个点为远离点;步骤5:将所有窄带点移入窄带内,构建初始窄带;步骤6:扩展窄带,直至窄带为空;在这一过程中,网格点走时是通过迎风差分法求解二维程函方程获得的,所述二维程函方程为:|▽τ|=s其中,τ为地震波走时,s为慢度,▽为梯度符号,上述公式中梯度项的迎风差分表达式为:...
【专利技术属性】
技术研发人员:孙辉,孙章庆,张志厚,张平,孟繁昌,邹时贵,夏龙龙,邓瑞,李梦,黎佳宾,刘丁毅,钟雨田,张杰,
申请(专利权)人:西南交通大学,
类型:发明
国别省市:四川,51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。