一种基于三维扫描地质雷达的目标体坐标估计方法技术

本发明专利技术提供了一种基于三维扫描地质雷达的目标体坐标估计方法，该方法对地质雷达数据进行一维滤波，然后再通过信号极大值约束得到水平边界角和俯仰边界角，最后通过三维空间坐标转换得到地下目标体的三维边界坐标和中心坐标。本申请提出的一种基于三维扫描地质雷达的目标体坐标估计方法，能够准确获取地下目标体的三维边界坐标及中心坐标，有利于三维空间的目标体成像和解释。的目标体成像和解释。的目标体成像和解释。

【技术实现步骤摘要】
一种基于三维扫描地质雷达的目标体坐标估计方法


[0001]本专利技术属于地球物理探测领域，具体为一种基于三维扫描地质雷达的目标体坐标估计方法，该方法以三维超前探测地质雷达数据为基础，依据地质雷达探测区域在三维空间的球体分布特征，估计对应地下目标体的空间坐标。

技术介绍

[0002]在煤矿开采过程中，容易引发地质构造变化，隐伏的地质构造异常(如断层、陷落柱、褶曲等)极易导致矿井开采事故。主要的开采事故包括顶板破坏、煤和瓦斯突出、矿井涌水等。地质雷达作为有效的物探手段，适用于探测断层、陷落柱、褶曲等地质构造异常。三维空间扫描地质雷达通过旋转地质雷达天线的角度，用于工作面前方三维空间探测，也可用于侧帮、回采区等区域的定向探测，实现了地质构造的快速识别和三维成像。
[0003]三维空间扫描地质雷达向工作面发送脉冲形式的高频宽带电磁波，一部分经发射天线直接到达接收天线形成耦合波；经发射天线发射后遇空气与路面分界面后未进入地下而直接反射回来被接收天线接收，形成直达波，耦合波和直达波常常被合称直耦波，常常用来判断时间零点，作为地下目标体深度的判断参考依据，另一部分电磁波进入工作面前方，在介质中传播，当遇到存在电性差异的目标体(如断层、陷落柱、褶曲等)或不同介质界面时，电磁波便产生反射，反射电磁波返回到工作面时由接收天线所接收。接收信号经过采集、处理、存储，就成为包含有工作面前方介质信息的地质雷达数据。
[0004]三维空间扫描地质雷达通过旋转地质雷达天线的角度，用于工作面前方三维空间探测。测线在三维空间的形态为弧形，该弧形以天线位置为圆心。因此，地质雷达探测区域应呈现球体分布特征。通过本专利技术一种基于三维扫描地质雷达的目标体坐标估计方法可以估计对应地下目标体的空间三维坐标，更清晰准确的描述探测区域内的目标体范围及分布特征。

技术实现思路

[0005]本专利技术的目的是以三维超前探测地质雷达数据为基础，依据地质雷达探测区域在三维空间的球体分布特征，估计对应地下目标体的空间坐标。
[0006]本方法按照以下步骤获取地下目标体的三维边界坐标及中心坐标：
[0007]步骤(A1)：对地质雷达数据D(α
i
,β
j
,k)，其中i＝0,1,...,N
‑
1；j＝0,1,...,N
‑
1；k＝1,2,...,H，α
i
为水平角，β
j
为俯仰角，k为采样点序号，H为采样点总数，首先进行一维滤波，对于某一水平角α
i
和俯仰角β
j
下的一维雷达数据D(α
i
,β
j
,k)，经过滤波函数F()滤波后得到
[0008]D
f
(α
i
,β
j
,k)＝F(D(α
i
,β
j
,k))
[0009]其中，D(α
i
,β
j
,k)经过滤波后的数据为D
f
(α
i
,β
j
,k)；
[0010]步骤(A2)：对于滤波后的雷达数据D
f
(α
i
,β
j
,k)，不妨令采样点序号k＝1，信号极大
值为max(D
f
(α
i
,β
j
,1))，求解
[0011][0012]其中，λ为常数，得到对应的i＝U1,j＝E1，其中U1∈{0,1,...,N
‑
1},E1∈{0,1,...,N
‑
1}，水平边界角为和俯仰边界角为依次求出不同采样点序号下的水平边界角和俯仰边界角，得到步骤(A3)：对于通过三维空间坐标转换，其三维坐标可以表示为其中，v为地下介质中的电磁波平均波速，Δt为单位采样点对应的探测时窗；
[0013]步骤(A4)：依次求出的三维坐标，维坐标，将这H个三维坐标连接起来，得到地下目标体的三维边界坐标；
[0014]步骤(A5)：三维边界坐标的均值为地下目标体的中心坐标。
[0015]本专利技术有以下优点：
[0016]1、通过本专利技术提出的一种基于三维扫描地质雷达的目标体坐标估计方法，能够准确获取地下目标体的三维边界坐标及中心坐标，有利于三维空间的目标体成像和解释。
附图说明
[0017]图1本专利技术技术方法流程图
[0018]图2本专利技术地质雷达数据三维空间成像
具体实施方式
[0019]本专利技术的目的是以三维超前探测地质雷达数据为基础，依据地质雷达探测区域在三维空间的球体分布特征，估计对应地下目标体的空间坐标。
[0020]本专利技术一种基于三维扫描地质雷达的目标体坐标估计方法分为两类情形，在第一种情形下，经过一维滤波，探测区域内有且只有1个目标体，此时具体步骤如下：
[0021](1)对于滤波后的地质雷达数据D
f
(α
i
,β
j
,k)，不妨令采样点序号k＝1，信号极大值为max(D
f
(α
i
,β
j
,1))，求解
[0022][0023]其中，λ为常数，得到对应的i＝U1,j＝E1，其中U1∈{0,1,...,N
‑
1},E1∈{0,1,...,N
‑
1}，水平边界角为和俯仰边界角为依次求出不同采样点序号下的水平边界角和俯仰边界角，得到
[0024]在第二种情形下，经过一维滤波，探测区域内有2个及以上目标体，此时具体步骤如下：
[0025](1)对于滤波后的地质雷达数据D
f
(α
i
,β
j
,k)，不妨令采样点序号k＝1，信号极大值为max(D
f
(α
i
,β
j
,1))，求解
[0026][0027]其中，λ为常数，得到对应的i＝U1,j＝E1，其中U1∈{0,1,...,N
‑
1},E1∈{0,1,...,N
‑
1}，水平边界角为和俯仰边界角为依次求出不同采样点序号下的水平边界角和俯仰边界角，得到
[0028](2)将上一步骤中水平边界角和俯仰边界角范围内的数据置为0，得到新的雷达数据D
′
f
(α
i
,β
j
,k)，不妨令采样点序号k＝1，信号极大值为max(D
′
f
(α
i
,β
j
,1))，求解
[0029][0030]其中，λ为常数，得到对应的i＝U1′
,j＝E1′
，其中U1′
∈{0,1,...,N
‑
1},E1′
∈{0,1,..本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于三维扫描地质雷达的目标体坐标估计方法，该方法以三维超前探测地质雷达数据为基础，依据地质雷达探测区域在三维空间的球体分布特征，估计对应地下目标体的空间坐标，具体步骤如下：步骤(A1)：对地质雷达数据D(α
i
,β
j
,k)，其中j＝0,1,...,N
‑
1；k＝1,2,...,H，α
i
为水平角，β
j
为俯仰角，k为采样点序号，H为采样点总数，首先进行一维滤波，对于某一水平角α
i
和俯仰角β
j
下的一维雷达数据D(α
i
,β
j
,k)，经过滤波函数F()滤波后得到D
f
(α
i
,β
j
,k)＝F(D(α
i
,β
j
,k))其中，D(α
i
,β
j
...

【专利技术属性】
技术研发人员：乔旭，杨峰，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：