摘要
To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algorithm.Our procedure is:(1) Subdivide a triangle grid based on the Delaunay triangular subdivision criterion and the relationships of the points,lines,and the surfaces in the subdividing area.(2) Define the topology relationships and related concepts of triangular unit ray tracing.(3) The source point and wave arrival points at any time compose the propagating plane wave and the minimum traveltime and secondary source positions are calculated during the plane wave propagation.We adopt the hyperbolic approximation global algorithm for secondary source retrieving.(4) By minimum traveltime ray tracing,collect the path from receiver to source points with the neighborhood point's traveltime and the direction of the secondary source.Numerical simulation examples are given to test the algorithm.The results show that the triangular net ray tracing method demonstrates model subdivision flexibility,precise velocity discontinuity interfaces,and accurate computations.
针对矩形网射线追踪存在的模型剖分灵活性差、速度界面描述精度差等问题,研究了复杂结构三角网最小走时射线追踪全局算法。(1)根据剖分区域点、线、面的结构关系,遵循Delaunay三角剖分的优化准则进行三角网格剖分;(2)定义三角单元射线追踪的拓扑关系;(3)波源点及某一时刻波到达的每一个节点点构成波行面,在波行面扩展过程中计算节点的最小走时和次级源位置,实际次级源检索采用双曲线近似算法;(4)利用各节点走时和次级源方向信息,通过最小走时搜索,拾取从接收点到源点的射线路径。数值模拟结果表明,三角网射线追踪方法模型剖分时灵活性强、速度间断面的描述精度高,追踪结果准确。
基金
the National Natural Science Foundation of China(Grant Nos.50974081,50774051, 51034003)
the Foundation for the Author of National Excellent Doctoral Dissertation of PR China(Grant No.200958)
the Program for Changjiang Scholars and Innovative Research Team in University(Grant No.IRT0843)
the Open Research Fund Program of the State Key Laboratory of Coal Resources and Safe Mining(China University of Mining and Technology, Beijing)(No.2007-07)