一种用于地质雷达层析探测的直射线追踪算法
【专利摘要】本发明涉及一种用于地质雷达层析探测的直射线追踪算法，属于地球物理反演算法领域，该方法包括如下步骤：建立平面直角相对坐标系；对剖面进行网格划分与编号；确定各网格的左右标点坐标；确定发射点与接收点的坐标；计算射线参数；计算射线在各网格内的长度；检验计算结果；将计算结果进行稀疏存储。本发明根据射线的左右截距与斜率三个参数追踪射线走向，只需计算射线经过的网格，相比传统算法遍历所有网格，大大降低了对内存和时间的消耗，提高了算法的效率与精度，最后将计算结果以稀疏格式存储，可节省存储空间并直观显示射线的分布状况。
【专利说明】一种用于地质雷达层析探测的直射线追踪算法
【技术领域】
[0001]本发明属于地球物理反演算法领域，涉及一种用于地质雷达层析探测的直射线追踪算法。
【背景技术】
[0002]地质雷达技术是一种快速、无损、高效的地球物理探测技术，近年来在矿井的应用逐渐增多，如煤矿含水层、陷落柱、裂隙带和瓦斯突出等灾害隐患的探测。目前地质雷达主要采用反射回波技术，探测深度局限于30-50米。由于现代化矿井回采区巷道之间的距离往往大于200米，使得地质雷达在探测采煤工作面内灾害源时无法达到所需要的探测深度。
[0003]层析探测技术是地质雷达领域一个新的方向，通过采集直达波的走时或能量进行反演计算，得到探测剖面的地质结构信息，探测距离可达到反射回波技术的2倍以上。反演计算的实质是求解一个矩阵方程组，其系数矩阵为射线在剖面离散化后的网格模型中各个网格内的长度。因此，射线追踪结果的精度与反演计算的精度密切相关。传统的射线追踪算法是针对每一条射线，均遍历整个网格模型，由于网格数目众多，并且每条射线经过的网格个数占整个模型的网格个数的比例很小，射线在大多数网格内的长度均为0，导致算法不仅存储量大、计算耗时长，并且效率很低。
【发明内容】
[0004]本发明所要解决的技术问题是针对现有射线追踪算法存储量大、计算耗时长，并且效率很低的缺点而提供一种用于地质雷达层析探测的直射线追踪算法。
[0005]本发明解决其技术问题采用的技术方案是:一种用于地质雷达层析探测的直射线追踪算法,具体包括:
51、建立平面直角相对坐标系:将地质雷达探测到的剖面全部置于坐标系的第一象限，其中层析探测方向平行于Z轴，测线方向平行于_7轴；坐标轴原点位于剖面顶点；层析探测垂直距离为4，剖面的测线长度为Ζ,，则剖面的四个端点的坐标分别为(0，O), {Lx,O), (O, L〉, {Lx, L》；
52、对剖面进行网格划分与编号:在剖面内绘制边长为D的正方形网格，以坐标轴原点处的网格为起始，编号为I ;进一步的以^轴方向为行，轴方向为列，逐行依次编号，从而将剖面划分为#个网格；
53、确定各网格的左右标点坐标:将网格4个顶点中ymJt)的顶点确定为左标点，Cr-，YmJ的顶点确定为右标点，第i个网格的左标点的坐标记为Ori, 7λ.)，第i个网格的右标点的坐标记为irx” ryt) , i=l, 2，— ,N ；
54、确定发射点与接收点的坐标:根据发射点、接收点与探测剖面的相对位置，确定发射点、接收点在坐标系中的位置，发射点坐标记为7>,)，接收点坐标记为Ofer, Re)；进一步的定位射线的发射点所在的网格编号/^与接收点所在的网格编号;S5、计算射线参数-Ji1= {Rey - Trr) / {Rex - Trx)，k2 =
【权利要求】
1.一种用于地质雷达层析探测的直射线追踪算法，其特征在于，所述方法包含如下步骤: 建立平面直角相对坐标系:将地质雷达探测到的剖面全部置于坐标系的第一象限，其中层析探测方向平行于Z轴，测线方向平行于_7轴；坐标轴原点位于剖面顶点；层析探测垂直距离为4，剖面的测线长度为Z,，则剖面的四个端点的坐标分别为(O，O), {Lx, O), (O,y)，(Ζχ，; 51、对剖面进行网格划分与编号:在剖面内绘制边长为D的正方形网格，以坐标轴原点处的网格为起始，编号为I ;进一步的以Z轴方向为行，轴方向为列，逐行依次编号，从而将剖面划分为#个网格； 52、确定各网格的左右标点坐标:将网格4个顶点中、xmin,ymin)的顶点确定为左标点，Cr-，y,J的顶点确定为右标点，第i个网格的左标点的坐标记为Ori, 7λ.)，第i个网格的右标点的坐标记为irx” ryt) , i=l, 2，— ,N ； 53、确定发射点与接收点的坐标:根据发射点、接收点与探测剖面的相对位置，确定发射点、接收点在坐标系中的位置，发射点坐标记为7>,)，接收点坐标记为Ofer, Re)；进一步的定位射线的发射点所在的网格编号/^与接收点所在的网格编号;
54、计算射线参数-Ji1= {Rey - Try) / {Rex - Trx), k2 =
2.如权利要求1所述的方法，其特征在于，所述步骤SI中的剖面为二维矩形剖面。
3.如权利要求1所述的方法，其特征在于，所述步骤S4中的发射点的横坐标均为0，接收点的横坐标均为4。
4.如权利要求1所述的方法，其特征在于，所述步骤S4中定位发射点所在的网格编号Tl1的方法为:
5.如权利要求1所述的方法，其特征在于，所述步骤S4中定位接收点所在的网格编号n2的方法为:
6.如权利要求1所述的方法，其特征在于，所述步骤S6中计算之=O的射线在各网格内的长度的方法为:若/?,≤/?，射线在网格编号在[/?,，n2]范围的网格内的长度均为D ;若射线在网格编号在[/?, Jij]范围的网格内的长度均为D ;其他网格内的长度均为O。
7.如权利要求1所述的方法，其特征在于，所述步骤S6中计算之^O的射线在各网格内的长度的方法为:.561、计算射线与网格左边界所在直线的交点的纵坐标力，根据公式:4=?-?计算左截距 .562、计算射线与网格/7,右边界所在直线的交点的纵坐标&，根据公式:计算右截距4; .563、根据如下公式计算射线在网格内的长度:
【文档编号】G01V7/02GK104020508SQ201410265945
【公开日】2014年9月3日 申请日期:2014年6月16日 优先权日:2014年6月16日
【发明者】杨峰, 杜翠, 郑晶, 彭猛, 张媛媛, 刘童, 陈爱梅 申请人:中国矿业大学(北京)