1. 引言
CORS-RTK技术给现代矿山测量带来技术变革 [1] [2] [3] [4] 。CORS-RTK技术较之传统测量技术可以提高测量精度、减少工作量、提高工作效率、节约资金等优势,已在矿山测量中获得广泛应用 [3] - [9] 。但,CORS-RTK测量高程成果的精度显著低于平面测量精度,特别是正常高测量精度远低于传统几何水准测量精度,大大地限制了其在煤矿岩移沉降监测上的应用 [2] [10] - [15] 。
基于此,本文拟基于矿区大似地水准面和区间高程异常变化率构建CORS-RTK高程修正模型,并进行实例验证,以期拓展CORS-RTK的应用范围,获得其在岩移沉降监测方面的广泛应用。
2. 修正模型
CORS-RTK测量可获得监测点平面位置及高程信息,其中高程为以大地水准面为参考面的大地高经由仪器内置模型拟合出监测点高程异常进而转换为点位正常高,大地高与正常高关系见式(1) [11] [12] [13] [14] 。但由于CORS-RTK内置转换模型通常为大范围区域似大地水准面模型,有时无法精确反应小区域内的高程异常变化趋势,故需对CORS-RTK测量中高程异常变化项进行改正。
(1)
2.1. 高程异常变化率修正模型
地面任意两控制点a、b之间的高程异常变化率为:
(2)
式中,
为
、
方向的高程异常变化率,单位为
;
为
、
点间的平距,单位为
;
为
、
点间的大地高差,单位为
;
为
、
点间的正常高高差,单位为
。利用控制点计算两点间的高程异常变化率,为了更准确,可以是多点间高程异常变化率的平均值。
根据式(2),a、b两点间的高程异常差为:
(3)
其中,
。
若在高程异常变化方向(ab)布设岩移观测点,经RTK测量后可得到观测点j与a或者b的平距
或
,由(3)式可得观测点与已知点的高程异常差
或
。则CORS-RTK的高程异常改正为:
(4)
式中,
为
点改正后的高程异常值,
或
为
或
点实测高程异常值,
或
为
点高程异常修正值。结合(1)式,则CORS-RTK测量的高程修正模型为:
(5)
2.2. 二次曲面拟合修正模型
二次曲面高程异常值拟合模型,即:
(6)
式中,
为拟合系数,
为
点的高程异常值,至少需要6个已知控制点高程异常值。
通过矿区及临近区域已知点GPS数据和相应水准数据对(6)式进行解算,得到区域高程异常模型。在岩移观测点j上进行CORS-RTK观测,获得观测点的平面位置信息。通过训练好的二次曲面拟合模型解算出观测点的高程异常值
,以
作为j点RTK的高程异常修正值即
。j点正常高有:
(7)
3. 算例及分析
3.1. 测区概况及已有数据
煤矿位于济宁市某县,矿区所在位置地面起伏较为平坦,多为农田。为了研究地表岩移规律,在二采区地表建立地表岩移观测站,沿工作面走向方向布设一条观测线N,设置了62个观测点,在观测线N的中点且垂直与N的方向上布设观测线E,设置了58个观测点。数据采集采用济南CORS基站,应用RTK技术进行平面和高程测量。开采面沉降区域外共有11个C级GPS已知点且均进行了三等水准联测。测区已知点分布、二采区地表范围(矩形区域)及观测线见图1。
![](//html.hanspub.org/file/4-2840182x53_hanspub.png)
Figure 1. Second mining area range and distribution of known point
图1. 二采区范围及已知点分布
3.2. 修正计算
为验证模型修正准确度与两种模型的修正优劣,在考虑测区已有数据的前提下进行以下方案设计:利用测区已知数据分别计算出二采区的高程异常拟合模型,选取二采区地表观测线N方向3个岩移观测点(A01、A02、A03)与E方向3个岩移观测点(S01、S02、S03)进行修正计算,计算结果与6个岩移观测点的四等水准测量结果对比。
1) 高程异常变化率修正计算
依据图1,在二采区地表沉降范围外N、E观测线两端布设四个GPS控制点与测区已知点进行C级GPS静态测量并进行三等水准联测,经网平差计算后作为二采区岩移观测的控制点。二采区控制点分布图见图2。
利用式(1)、(2)计算观测线方向高程异常变化率,计算结果见表1。
根据表1和式(2)~(5),可得6个CORS-RTK岩移观测点的高程异常修正值,及修正后的正常高,计算结果见表2。
2) 二次曲面拟合修正计算
利用测区已知数据、式(1)及最小二乘法可得测区二次曲面高程异常拟合模型
(8)
![](//html.hanspub.org/file/4-2840182x55_hanspub.png)
Figure 2. Distribution of control points in the second mining area
图2. 二采区控制点分布
![](Images/Table_Tmp.jpg)
Table 1. Calculation results of elevation abnormal change rate
表1. 高程异常变化率计算结果
利用式(8)及6个RTK岩移观测点观测数据可得各点高程异常修正值,及修正后的正常高,计算结果见表3。
3) 精度分析
据表2、表3,可得二采区地表观测线N、E方向上的正常高残差中误差以及与其相对应的四等水准测量限差。计算结果见表4。
3.3. 结果分析
据表2,高程异常变化率修正模型最大修正残差为4 mm,最小为2 mm,残差中误差为3.4 mm。据表3,二次曲面修正模型最大修正残差为4 mm,最小为1 mm,残差中误差为2.9 mm。对比表1、表2,两种修正模型均能达到良好的修正效果。据表3,二采区地面沉降区在N、E方向上经过两种模型修正后的CORS-RTK点正常高均能达到四等水准测量的要求,且二次曲面修正模型总体上优于高程异常变化率修正模型。
![](Images/Table_Tmp.jpg)
Table 2. Correction results of CORS-RTK rock-movement monitoring points
表2. CORS-RTK岩移监测点修正结果
![](Images/Table_Tmp.jpg)
Table 3. Correction results of CORS-RTK rock-movement monitoring points
表3. CORS-RTK岩移监测点修正结果
4. 结束语
CORS-RTK技术在矿山测量中具有明显技术优势,鉴于其在煤矿岩移沉降监测方面的限制,构建了利用CORS-RTK进行较高精度正常高测量修正模型。利用实例进行了验证,结果表明:在合理限制和外业观测条件良好的情况下,利用修正模型获得了测点正常高精度的明显提高,可以达到四等几何水准测量的精度。该模型可用于煤矿岩移沉降监测,具有很好的推广应用价值。
NOTES
*通讯作者。