趙喜江, 楊新宇, 霍大亮

(1.黑龍江科技大學(xué) 發(fā)展規(guī)劃與高等教育研究處, 哈爾濱 150022; 2.黑龍江科技大學(xué) 礦業(yè)工程學(xué)院, 哈爾濱 150022)

0 引 言

GNSS 技術(shù)是地面采集單點(diǎn)空間坐標(biāo)最為高效、快速的方式[1],但是通過 GNSS所獲取到的高程值并不能直接應(yīng)用于實(shí)際礦山測(cè)量工作,原因是GNSS 技術(shù)獲取的點(diǎn)的高程是依賴于參考橢球 WGS-84 坐標(biāo)系的大地高程,而我國礦山測(cè)量工作中,所采用的高程是以似大地水準(zhǔn)面為基準(zhǔn)的正常高或正高,在大地高與正常高之間存在差值,此差值即為高程異常,因此,利用RTK技術(shù)直接獲取高程時(shí),精度無法滿足礦區(qū)地表1∶500、1∶1 000比例尺測(cè)圖和開采沉陷監(jiān)測(cè)精度要求。隨著GNSS衛(wèi)星定位技術(shù)、現(xiàn)代通信技術(shù)、空間信息技術(shù)以及網(wǎng)絡(luò)技術(shù)的快速發(fā)展,許多大型礦區(qū)都在陸續(xù)建設(shè)覆蓋本礦區(qū)的CORS系統(tǒng),但是礦區(qū)建設(shè)的CORS系統(tǒng)時(shí),站點(diǎn)的初始參數(shù)由設(shè)備提供商的技術(shù)支持人員進(jìn)行簡(jiǎn)易設(shè)置,一般沒有充分考慮局部區(qū)域地理位置和地形地勢(shì)的影響,高程異常值精度一般為分米級(jí),而大型礦區(qū)地表變形監(jiān)測(cè)和工程測(cè)量的精度要求厘米級(jí),導(dǎo)致礦區(qū)CORS-RTK測(cè)量高程精度無法滿足測(cè)繪生產(chǎn)要求,制約了RTK技術(shù)在礦區(qū)地表大比例尺測(cè)圖和開采沉陷監(jiān)測(cè)工作中的應(yīng)用。因此,對(duì)于高程變化較大的大型礦區(qū),通過似大地水準(zhǔn)面精化[2]得到礦區(qū)高程異常值的最或是值來調(diào)整CORS系統(tǒng)的參數(shù)設(shè)置,對(duì)于提升RTK測(cè)量高程精度具有重大的工程意義。

表1 靜態(tài)約束平差基線最弱邊

1 研究工作與技術(shù)路線

1.1 測(cè)區(qū)概況

鶴崗市處于三江平原向小興安嶺山地過渡的明顯上升地段,地勢(shì)西北高東南低,可分為低山丘陵、漫崗、平原、溝谷及漫灘四種地貌類型,地勢(shì)由西北向東南傾斜,西北部為山區(qū),最高山峰小白山海拔1 022 m,境內(nèi)海拔200 m以上的山峰有34座。有梧桐河、鶴立河、阿凌達(dá)河、嘉蔭河等18條主要河流從西北流向東南。市區(qū)南部為丘陵漫崗,東南部為三江平原的邊緣地帶,平坦開闊,平均海拔80 m左右。鶴崗礦區(qū)位于鶴崗市東部,地理位置為東126°36′44″～127°31′46″,北緯45°40′20″～45°11′38″之間,從北端的平頂山到南端新華北山直線距離約50 km,從東端梧桐河?xùn)|到西端的新水源約30 km。礦區(qū)西北部和東南以丘陵地形為主,南部以河漫灘為主,礦區(qū)高程變化較大,地形地貌錯(cuò)綜復(fù)雜。測(cè)區(qū)道路較窄,時(shí)有大量運(yùn)輸煤炭的大型貨車通行,通視條件一般,通信信號(hào)覆蓋較好。

1.2 工作內(nèi)容

本次研究工作在鶴崗礦區(qū)2021年建成的14個(gè)C級(jí)GNSS控制點(diǎn)(圖1)和3個(gè)CORS站等已有控制網(wǎng)觀測(cè)成果基礎(chǔ)上,依據(jù)鶴崗市區(qū)2020年7月建立的新鶴II-33二等水準(zhǔn)點(diǎn),參考富力煤礦辦公樓、興安南橋頭、峻德礦辦公樓、二跨橋、公鐵交叉口、水泥廠、梧桐河水文站、新華北山等原有礦區(qū)三等水準(zhǔn)高程控制點(diǎn),通過對(duì)礦區(qū)14個(gè)C級(jí)GNSS控制網(wǎng)點(diǎn)中的356ps、新水源、梧桐河?xùn)|、鳥山、東線場(chǎng)、新華30隊(duì)等六個(gè)外圍點(diǎn)和礦務(wù)局、峻德救護(hù)隊(duì)、鳥山等3個(gè)CORS站進(jìn)行三等水準(zhǔn)聯(lián)測(cè),再利用C級(jí)GNSS控制點(diǎn)中平頂山、陰日山和新華北山三個(gè)已知三等水準(zhǔn)高程值,進(jìn)行礦區(qū)高程平差結(jié)算和高程擬合,完成礦區(qū)似大地水準(zhǔn)面精化,得到全部14個(gè)控制點(diǎn)三等水準(zhǔn)高程和3個(gè)CORS站的高程異常值,進(jìn)而調(diào)整CORS站高程異常參數(shù)和RTK參數(shù)設(shè)置。為了檢驗(yàn)精度,在礦區(qū)主要礦井工業(yè)廣場(chǎng)附近布設(shè)10個(gè)近礦點(diǎn),平面用E級(jí)GNSS控制點(diǎn)施測(cè),高程用三等水準(zhǔn)施測(cè)。通過RTK實(shí)際10個(gè)近礦點(diǎn)高程來進(jìn)行精度驗(yàn)證。

圖1 鶴崗礦區(qū)C級(jí)GNSS控制點(diǎn)和CORS站點(diǎn)示意

1.3 技術(shù)路線

利用前期項(xiàng)目實(shí)測(cè)礦區(qū)首級(jí)GNSS控制網(wǎng)及CORS站坐標(biāo)和大地高程,采用三等水準(zhǔn)測(cè)量的方法聯(lián)測(cè)了礦區(qū)外圍控制點(diǎn)高程,完成了礦區(qū)高程平差計(jì)算,通過平面模型法、二次曲面模型法、多面函數(shù)擬合法高程擬合比較,優(yōu)選了結(jié)合EGM2008 模型的組合法,運(yùn)用“移去-恢復(fù)”法精化礦區(qū)似大地水準(zhǔn)面模型[3-5],得到鶴崗礦區(qū)全部14個(gè)首級(jí)GNSS控制點(diǎn)及3個(gè)CORS站處的高程及高程異常值,通過調(diào)整CORS高程異常參數(shù)設(shè)置,在礦區(qū)主要礦井的近礦點(diǎn)上通過實(shí)際高程測(cè)量進(jìn)行精度驗(yàn)證。

2 平面及高程控制測(cè)量

2.1 GNSS近礦點(diǎn)選點(diǎn)遵循的原則

(1)在10°～15°高度角以上不能有成片障礙物,點(diǎn)的周圍200 m范圍內(nèi)不能有強(qiáng)電磁波干擾源,如大功率無線電發(fā)射設(shè)施、高壓輸電線等。

(2)為避免或減少多路徑效應(yīng)的發(fā)生,近礦點(diǎn)應(yīng)遠(yuǎn)離對(duì)電磁波信號(hào)反射強(qiáng)烈的地形、地物,如高層建筑、成片水域等。

(3)為保證后期驗(yàn)證實(shí)驗(yàn)質(zhì)量,近礦點(diǎn)選在覆蓋整個(gè)鶴崗礦區(qū)的鳥山、益新、富力、興安、竣德等主力生產(chǎn)礦井工業(yè)廣場(chǎng)附近,保證交通便利且點(diǎn)位應(yīng)易于保存處。

(4)相鄰點(diǎn)間最小距離應(yīng)為平均距離的1/2～1/3;相鄰點(diǎn)間最大距離應(yīng)為平均距離的2～3倍;特殊情況下個(gè)別點(diǎn)的間距允許超出精度和邊長(zhǎng)表中的規(guī)定。

2.2 GNSS近礦點(diǎn)布設(shè)

為了檢驗(yàn)研究成果,以礦務(wù)局、峻德救護(hù)大隊(duì)2個(gè)CORS站為起始點(diǎn),在礦區(qū)5個(gè)主力生產(chǎn)礦井的工業(yè)廣場(chǎng)附近布設(shè)鳥山1、2,益新1、2,富力1、2,興安1、2,峻德5、6等10個(gè)E級(jí)GNSS近礦平面和高程控制點(diǎn),作為實(shí)驗(yàn)高程精度驗(yàn)證點(diǎn),如圖2所示。

圖2 鶴崗礦區(qū)主力礦井E級(jí)GNSS近礦點(diǎn)分布

2.3 GNSS近礦控制點(diǎn)平面控制測(cè)量

2.3.1 外業(yè)施測(cè)

使用南方極點(diǎn)RTK慣導(dǎo)版C對(duì)近礦點(diǎn)進(jìn)行平面控制測(cè)量,各靜態(tài)點(diǎn)觀測(cè)不少于2個(gè)時(shí)段,每個(gè)時(shí)段觀測(cè)45 min,與鶴崗市內(nèi)CORS站共同進(jìn)行解算。

2.3.2 內(nèi)業(yè)數(shù)據(jù)處理及精度

解算軟件采用南方地理數(shù)據(jù)處理平臺(tái)軟件SGO(簡(jiǎn)稱SGO)。對(duì)于相鄰點(diǎn)距離較近的短邊保證進(jìn)行同步觀測(cè)。靜態(tài)觀測(cè)及解算精度符合CJJ 8—010《城市測(cè)量規(guī)范》、CJJ/T73—010《衛(wèi)星定位城市測(cè)量技術(shù)規(guī)范》要求。

GNSS近礦控制點(diǎn)的內(nèi)業(yè)平差計(jì)算過程如圖3所示,近礦控制點(diǎn)平面精度,其中靜態(tài)約束平差基線最弱邊(表1)和平面最弱點(diǎn)(表2),基線和點(diǎn)位誤差滿足規(guī)程規(guī)定。

圖3 靜態(tài)約束平差基線最弱邊

表2 靜態(tài)約束平差平面最弱點(diǎn)

2.4 礦區(qū)高程控制測(cè)量

2.4.1 水準(zhǔn)測(cè)量外業(yè)

根據(jù)研究工作的精度需要,課題組采用三等水準(zhǔn)建立了覆蓋礦區(qū)的高程控制網(wǎng)。在礦區(qū)已有高程控制點(diǎn)現(xiàn)場(chǎng)踏勘和現(xiàn)狀調(diào)查基礎(chǔ)上,確定了采用原礦區(qū)11個(gè)C級(jí)GNSS控制點(diǎn)和本次選定的10個(gè)E級(jí)GNSS近礦點(diǎn)作為平面控制點(diǎn)的同時(shí)兼做礦區(qū)高程控制點(diǎn),采用CGCS2000坐標(biāo)系和1985國家高程基準(zhǔn)。為了保證精度,堅(jiān)持每測(cè)站觀測(cè)程序“后前前后”“單程雙站(變儀器高)”觀測(cè)。三絲能讀數(shù),視距小于100 m,前后視距差和累計(jì)差不超限,累計(jì)實(shí)測(cè)三等水準(zhǔn)路線約80 km。

2.4.2 水準(zhǔn)測(cè)量精度檢查

選擇地形起伏較大的最長(zhǎng)閉合環(huán)鶴立河-興安南橋-公鐵交叉口-鶴立河的閉合水準(zhǔn)路線(圖4)來驗(yàn)證水準(zhǔn)施測(cè)精度,閉合差及限差見表3。

圖4 閉合水準(zhǔn)路線

表3 閉合線路閉合差及限差統(tǒng)計(jì)

3 似大地水準(zhǔn)面精化及效果驗(yàn)證

對(duì)平面擬合法、二次曲面擬合法、多面函數(shù)擬合法和結(jié)合EGM2008模型的組合法進(jìn)行對(duì)比試驗(yàn),確定上述方法在礦區(qū)內(nèi)的適用性和擬合精度,優(yōu)選出結(jié)合EGM2008模型的“移去-恢復(fù)”法來建立礦區(qū)似大地水準(zhǔn)面,求出高程異常值最或是值的方法。

3.1 結(jié)合EGM2008模型的“移去-恢復(fù)”

原理是基于近年由NGA(美國國家地理空間情報(bào)局)釋放EGM2008全球超高階地球重力場(chǎng)模型,采用該模型以及GPS/水準(zhǔn)數(shù)據(jù)采用“移去-恢復(fù)”法[6-7],利用礦區(qū)C級(jí)GNSS平面控制和三等水準(zhǔn)測(cè)量得到的坐標(biāo)和高程測(cè)量成果,獲得高精度的區(qū)域似大地水準(zhǔn)面。具體步驟如下。

(1)移去。選擇鶴崗礦區(qū)范圍內(nèi)具有 GPS 數(shù)據(jù)與水準(zhǔn)數(shù)據(jù)的17個(gè)重合點(diǎn),利用重合點(diǎn)的正常高和大地高求出各點(diǎn)對(duì)應(yīng)的高程異常值,利用EGM2008 模型能夠得到各點(diǎn)的重力高程異常,據(jù)此可以計(jì)算出對(duì)應(yīng)點(diǎn)的剩余高程異常量。

(2)擬合剩余高程異常量。針對(duì)在對(duì)重力場(chǎng)模型的利用之下所獲得重力異常,通常而言我們會(huì)將其視為無誤差值。在此的擬合數(shù)據(jù)為這17個(gè)點(diǎn)的剩余高程異常,采用擬合函數(shù)解算后,即可得到相應(yīng)的擬合參數(shù),將待求點(diǎn)的平面坐標(biāo)通過確定參數(shù)后的模型計(jì)算后,剩余高程異常值就能夠通過內(nèi)插形式所獲得。

(3)恢復(fù)。利用重力場(chǎng)模型求出重力高程異常,與內(nèi)插得到的剩余高程異常相加,就能夠得到各待求點(diǎn)的高程異常,據(jù)此可以求解出待求點(diǎn)的正常高。

3.2 鶴崗礦區(qū)高程異常等值線

通過結(jié)合EGM2008模型的“移去-恢復(fù)”,計(jì)算得到的17個(gè)點(diǎn)的高程異常值,進(jìn)而生成鶴崗礦區(qū)高程異常等值線圖(圖5)。

圖5 鶴崗礦區(qū)高程異常等值線

3.3 成果應(yīng)用及精度驗(yàn)證

3.3.1 高程異常改正

文中以鶴崗礦區(qū)LEICA GNSS Spider CORS網(wǎng)絡(luò)參考站高程異常值修正過程為例,說明礦區(qū)高程異常成果在 CORS系統(tǒng)的應(yīng)用。

①創(chuàng)建一個(gè)新的遠(yuǎn)程站點(diǎn)服務(wù)器連接。啟動(dòng)軟件,從“菜單欄”中點(diǎn)擊“服務(wù)器”,選擇“新建站點(diǎn)服務(wù)器”,輸入服務(wù)器名稱、服務(wù)器密碼,輸入服務(wù)器的PC名稱或者TCP/IP地址。

②配置通信設(shè)備和參數(shù)。打開創(chuàng)建的服務(wù)器,設(shè)置傳感器類型和傳感器連接方式,并且通過瀏覽日志欄檢查是否已經(jīng)建立了傳感器通訊。

③設(shè)置常規(guī)站點(diǎn)參數(shù)。輸入站點(diǎn)名稱,并通過下拉列表為該站點(diǎn)接收機(jī)上記錄的數(shù)據(jù)文件選擇下載間隔時(shí)間。

④配置接收機(jī)坐標(biāo)等信息。輸入本CORS站點(diǎn)精確的經(jīng)緯度、大地高和高程異常值,(高程異常值即為通過擬合計(jì)算求出的本CORS站點(diǎn)處的高程異常值),配置界面如圖6所示。

圖6 配置接收機(jī)信息

⑤配置完成并上傳。

⑥配置接收機(jī)記錄參數(shù)。設(shè)置采樣率、文件長(zhǎng)度和自動(dòng)刪除記錄文件的時(shí)間閾值,文件長(zhǎng)度為傳感器創(chuàng)建數(shù)據(jù)記錄的新文件的時(shí)間間隔。

⑦連接至流動(dòng)站。流動(dòng)站通過輸入TCP/IP端口與接收機(jī)進(jìn)行連接,連接成功后,即可接收相關(guān)的參數(shù)和數(shù)據(jù)。

⑧實(shí)測(cè)礦區(qū)地面點(diǎn)坐標(biāo)及高程。

根據(jù)“創(chuàng)建項(xiàng)目-設(shè)置參數(shù)-連接流動(dòng)站”的步驟?？梢酝ㄟ^該軟件中的“轉(zhuǎn)換”功能對(duì)相應(yīng)的參數(shù)、參考橢球、投影和大地水準(zhǔn)面模型/殘差進(jìn)行定義,并通過RTCM信息將它們波發(fā)至流動(dòng)站,流動(dòng)站不需要手動(dòng)設(shè)置任何轉(zhuǎn)換參數(shù)即可進(jìn)行坐標(biāo)轉(zhuǎn)換。

3.3.2 似大地水準(zhǔn)面精化效果檢驗(yàn)

為了驗(yàn)證結(jié)合EGM2008模型的“移去-恢復(fù)”法對(duì)礦區(qū)似大地水準(zhǔn)面精化的精度情況,課題組依據(jù)似大地水準(zhǔn)面精化成果對(duì)鶴崗礦區(qū)內(nèi)已經(jīng)建成的礦務(wù)局、峻德和鳥山3個(gè)CORS站點(diǎn)的高程異常值進(jìn)行重新賦值,并用RTK實(shí)測(cè)了分布在礦區(qū)內(nèi)5個(gè)煤礦的近礦點(diǎn)高程進(jìn)行檢驗(yàn),驗(yàn)證實(shí)驗(yàn)采用E級(jí)GNSS近礦點(diǎn)和CORS站點(diǎn)的實(shí)測(cè)三等水準(zhǔn)高程值與似大地水準(zhǔn)面精化高程擬合后RTK實(shí)測(cè)高程進(jìn)行對(duì)照檢驗(yàn)(表4)。

表4 實(shí)測(cè)三等水準(zhǔn)高程與似大地水準(zhǔn)面精化后RTK實(shí)測(cè)值對(duì)照統(tǒng)計(jì)

由實(shí)驗(yàn)結(jié)果可以看出,鶴崗礦區(qū)范圍內(nèi)依托CORS站的RTK采集高程點(diǎn)的誤差從不小于 20 cm提高到不大于2 cm。實(shí)驗(yàn)證明通過礦區(qū)似大地水準(zhǔn)面精化和高程擬合,可以顯著削弱高程異常對(duì)RTK高程測(cè)量帶來的影響。

4 結(jié) 論

(1)研究結(jié)果證明,基于礦區(qū)實(shí)測(cè)GNSS平面控制成果和三等水準(zhǔn)高程控制成果,結(jié)合EGM2008模

型,運(yùn)用“移去-恢復(fù)”法精化礦區(qū)范圍內(nèi)似大地水準(zhǔn)面模型的技術(shù)路線可行。

(2)諸如鶴崗礦區(qū)這類地表高程變化較復(fù)雜的大型礦區(qū),采用文中精化技術(shù)路線得到的礦區(qū)似大地水準(zhǔn)面模型和高程異常值,采用改進(jìn)后的CORS-RTK系統(tǒng)進(jìn)行高程測(cè)量誤差小于2 cm,可以極大地提高RTK高程測(cè)量的精度。

(3)改進(jìn)后的CORS-RTK系統(tǒng)可以較好的滿足礦區(qū)范圍內(nèi)大比例尺測(cè)圖和開采沉陷日常監(jiān)測(cè)的精度需要,極大地拓寬了RTK在礦區(qū)地面高程測(cè)量應(yīng)用場(chǎng)景。