吉林大学学报(地球科学版)

• 地球探测与信息技术 • 上一篇    下一篇

旅行时线性插值射线追踪提高计算精度和效率的改进方法

叶佩1 ,李庆春2   

  1. 1.中国国土资源航空物探遥感中心,北京100083;
    2.长安大学地质工程与测绘学院,西安710054
  • 收稿日期:2012-05-10 出版日期:2013-01-26 发布日期:2013-01-26
  • 作者简介:叶佩(1986-),男,助理工程师,主要从事综合物探方法处理与解释研究工作,E-mail:yepei886@126.com
  • 基金资助:

    国家“863”计划项目A类课题(2007AA06Z103)

Improvements of Linear Traveltime Interpolation Ray Tracing for the Accuracy and Efficiency

Ye Pei1,Li Qingchun2   

  1. 1.China Aero Geophysical Survey and Remote Sensing Center for Land and Resources, Beijing100083, China;
    2 College of Geology Engineering and Geomatics, Chang’an University, Xi’an710054, China
  • Received:2012-05-10 Online:2013-01-26 Published:2013-01-26

摘要:

旅行时线性插值(LTI)射线追踪算法是基于线性假设的,在向前处理过程中仅用按行或按列扫描的方法来计算节点旅行时没有考虑逆向传播射线,导致其计算精度与网格剖分大小有关,在处理复杂介质时会使得追踪出来的射线路径不一定满足最短旅行时。因此,笔者提出了两点改进措施:在向前处理时需采用全方位循环的方法来计算节点最小旅行时;在网格边界加入次生节点。模型试算结果表明:采用全方位循环的LTI法考虑了来自各个方向的射线,可提高其对复杂模型的适应能力;在节点间距相同的情况下,网格边界插入次生节点的LTI法较传统的LTI法计算精度至少可以提高一个数量级,同时,计算速度也更快;随着节点间距剖分的越精细,计算耗时下降也越明显,计算速度较传统的方法可提高n~10n倍。

关键词: 射线追踪, 旅行时线性插值, 全方位循环, 次生节点

Abstract:

The linear traveltime interpolation(LTI) method is based on linear assumption so that its precision depends on grid division, and the rays transported from the reverse direction are not considered when we recalculate traveltime row by row or column by column,and the computation is unstable in the complex models. Here, we put forward two improvements for the LTI method. First, calculating the minimum travel-time of all grid nodes by all-round cycle strategy during ray tracing. Secondly, inserting secondary nodes in the original grid boundary. Model tests show that the all-round cycle strategy of LTI considers the races from all direction, so it’s more stable for complex models.In the case of same grid division, comparing with the traditional LTI, the LTI method with secondary inserted nodes has higher precision that be improved at least an order of magnitude.At the same time, the calculation speed of our method is faster and time consuming decreases more obviously when the grids become smaller.The calculation speed of our method can be n-10n times than the traditional one’s.

Key words: ray tracing, linear traveltime interpolation, all-around cycle, secondary node

中图分类号: 

  • P631.4
[1] 孙章庆, 汪登科, 韩复兴. 复杂海底各种地震波的射线追踪与运动学特征[J]. 吉林大学学报(地球科学版), 2019, 49(4): 1169-1181.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 刘建峰,迟效国,周燕,王铁夫,金巍,周建波,董春艳,黎广荣. 小兴安岭东北部金林岩体全岩-角闪石Rb-Sr年龄[J]. J4, 2005, 35(06): 690 -0693 .
[2] 黄冠星, 孙继朝, 张英, 刘景涛, 张玉玺, 荆继红. 珠江三角洲污灌区地下水重金属含量及其相互关系[J]. J4, 2011, 41(1): 228 -234 .
[3] 肖长来,梁秀娟,崔建铭,兰盈盈,张君,李书兰,梁瑞奇,郑策. 确定含水层参数的全程曲线拟合法[J]. J4, 2005, 35(06): 751 -0755 .
[4] 郭振华,王璞珺,印长海,黄玉龙. 松辽盆地北部火山岩岩相与测井相关系研究[J]. J4, 2006, 36(02): 207 -0214 .
[5] 孙永河,付晓飞,吕延防,付广,阎冬. 地震泵抽吸作用与油气运聚成藏物理模拟[J]. J4, 2007, 37(1): 98 -0104 .
[6] 谢忠雷,陈卓,孙文田,尹波. 不同茶园茶叶氟含量及土壤氟的形态分布[J]. J4, 2008, 38(2): 293 -0298 .
[7] 贾军涛,王璞珺,邵 锐,程日辉,张 斌,侯景涛,李金龙,边伟华. 松辽盆地东南缘营城组地层序列的划分与区域对比[J]. J4, 2007, 37(6): 1110 -1123 .
[8] 康立明,任战利. 多参数定量研究流动单元的方法--以鄂尔多斯盆地W93井区为例[J]. J4, 2008, 38(5): 749 -0756 .
[9] 薛永超, 程林松. 白豹油田长8油藏成岩储集相[J]. J4, 2011, 41(2): 365 -371 .
[10] 姜纪沂, 张宇东, 谷洪彪, 左兰丽. 基于灰色关联熵的地下水环境演化模式判别模型研究[J]. J4, 2009, 39(6): 1111 -1116 .