孫勇，韓昱

（安徽省地球物理地球化學(xué)勘查技術(shù)院，安徽合肥 230022）

0 引言

在國(guó)家地質(zhì)找礦的大背景下，全球定位系統(tǒng)（Global Positioning System，GPS）作為一門相對(duì)于傳統(tǒng)測(cè)量而言的新興技術(shù)，在新世紀(jì)得到了快速發(fā)展，GPS 技術(shù)作業(yè)效率高、精度高的顯著優(yōu)點(diǎn)，在地質(zhì)找礦方面也得到迅速的普及［1］。本文就GPS 技術(shù)在物探測(cè)量方面的應(yīng)用，進(jìn)行系統(tǒng)闡述，對(duì)其優(yōu)缺點(diǎn)進(jìn)行總結(jié)，為今后技術(shù)應(yīng)用提出一點(diǎn)個(gè)人見解。

1 評(píng)價(jià)區(qū)自然地理概況

“安徽省懷遠(yuǎn)-蒙城一帶金、鐵多金屬礦產(chǎn)調(diào)查評(píng)價(jià)”項(xiàng)目為安徽省2016 年度第二批公益性地質(zhì)項(xiàng)目（皖國(guó)土資〔2016〕71 號(hào)），編號(hào)：2016-g-2-13。項(xiàng)目調(diào)查評(píng)價(jià)區(qū)（以下簡(jiǎn)稱“評(píng)價(jià)區(qū)”）主要位于安徽省懷遠(yuǎn)縣、蒙城縣一帶，東至固鎮(zhèn)縣城-鳳陽(yáng)縣劉府一線，西至利辛縣江集鎮(zhèn)-馬店鎮(zhèn)一線，南至利辛縣闞疃鎮(zhèn)-懷遠(yuǎn)縣馬城鎮(zhèn)一線，北至蒙城縣小澗鎮(zhèn)-固鎮(zhèn)縣城一線。包括闞疃集、尚塘集、樂土鎮(zhèn)、常興鎮(zhèn)、龍亢鎮(zhèn)、唐集鎮(zhèn)、古城等27幅1∶50 000圖幅，除查區(qū)東部懷遠(yuǎn)荊山、涂山等地有少量露頭，其余均為覆蓋區(qū)，總計(jì)面積約為6 924 km2。

區(qū)域內(nèi)地形開闊、地勢(shì)平坦，除東部荊山、涂山夾淮河而對(duì)峙外，其余均為平原，在平原地貌單元中，由于河流變遷、交互沉積和歷次黃河南泛覆蓋及人工開河筑壩等影響，局部地面不平整，具“大平、小不平”的特點(diǎn)。據(jù)此又分為湖地、崗坡地、灣地三種小類型。區(qū)內(nèi)分布有北淝河、渦河、阜蒙新河、芡河、茨淮新河、渦河等河流，其中渦河、芡河在懷遠(yuǎn)流入淮河。區(qū)內(nèi)地形等級(jí)據(jù)原地質(zhì)礦產(chǎn)部地球物理地球化學(xué)勘查局編制的《全國(guó)1∶50 000地形圖地形等級(jí)集》，綜合地形等級(jí)為2.5級(jí)。

評(píng)價(jià)區(qū)地處平原，區(qū)內(nèi)交通便利，S307 省道連接懷遠(yuǎn)、蒙城兩縣，S305 省道連接蒙城、利辛兩縣，G36寧洛高速、S203 省道、S308 省道等均從區(qū)內(nèi)穿過，徐（州）—淮（北）—阜（陽(yáng)）鐵路、徐（州）—宿（州）—蚌（埠）—合（肥）高鐵分別從調(diào)查區(qū)西北及東北角通過。村村通工程使得水泥路已覆蓋區(qū)內(nèi)所有行政村。

2 作業(yè)技術(shù)安排

針對(duì)本次作業(yè)的精度要求，遵循“由整體到局部，由高級(jí)到低級(jí)，先控制后碎部”的原則［2］，根據(jù)測(cè)區(qū)實(shí)際情況，優(yōu)先進(jìn)行GPS 控制測(cè)量，求取7 個(gè)參數(shù)后，再進(jìn)行RTK 碎部放樣。本次在馬店、劉圩、支子湖3 個(gè)區(qū)塊部署了面積性重力測(cè)量，工作比例尺為1∶25 000，網(wǎng)度為250 m×50 m。面積性物探工作設(shè)計(jì)重力點(diǎn)位中誤差≤±10 m，高程中誤差≤±0.2 m，采用RTK 代入?yún)?shù)后，進(jìn)行測(cè)點(diǎn)放樣。在馬店、劉圩、立倉(cāng)三個(gè)區(qū)塊部署了面積性地面高精度磁法測(cè)量，工作比例尺為1∶10 000，網(wǎng)度為100 m×40 m。采用手持GPS 代入?yún)?shù)后，完成測(cè)點(diǎn)放樣。設(shè)計(jì)要求手持GPS點(diǎn)位中誤差≤±10 m，高程中誤差≤±10 m；GPS 測(cè)點(diǎn)定位精度檢查采用同精度檢查，檢查點(diǎn)數(shù)量不少于總點(diǎn)數(shù)的3%，三維坐標(biāo)較差不大于測(cè)點(diǎn)中誤差的2倍。主要測(cè)量參數(shù)為：根據(jù)測(cè)區(qū)地形圖比例尺采用6°分帶，坐標(biāo)系采用1954 年北京坐標(biāo)系統(tǒng)，高程采用1985 國(guó)家高程基準(zhǔn)，投影采用高斯正形投影，根據(jù)測(cè)區(qū)經(jīng)緯度中央子午線選為117°，GPS 坐標(biāo)系統(tǒng)轉(zhuǎn)換參數(shù)分別是RTK 7 個(gè)參數(shù)和手持導(dǎo)航型GPS5個(gè)參數(shù)等[3]。

3 GPS控制測(cè)量

GPS 控制測(cè)量主要為地面物探測(cè)點(diǎn)放樣及高程測(cè)量提供物控點(diǎn)?？紤]到篩選的重點(diǎn)區(qū)塊相隔較遠(yuǎn)，在實(shí)際生產(chǎn)中，各區(qū)塊獨(dú)立組網(wǎng)施測(cè)。

3.1 馬店區(qū)塊

生產(chǎn)儀器為天寶r4-3 型GPS 雙頻接收機(jī)，通過了符合國(guó)家計(jì)量檢定要求的質(zhì)檢單位的檢定。

收集C 級(jí)控制點(diǎn)4 個(gè)，加密控制點(diǎn)1 個(gè)。加密控制點(diǎn)采用混凝土標(biāo)石不銹鋼十字絲標(biāo)志，GPS 控制點(diǎn)名稱按所在地名拼音首字母命名。

控制網(wǎng)布置采用圖形強(qiáng)度較好和作業(yè)效率較高的邊連式GPS 網(wǎng)，作業(yè)方法為GPS 靜態(tài)定位測(cè)量。GPS 接收天線安置對(duì)中誤差不大于3 mm（用鋼尺在互為120°的三個(gè)方向上，測(cè)量后取平均值），衛(wèi)星信號(hào)接收同步觀測(cè)時(shí)間不小于60 min。觀測(cè)衛(wèi)星高度截止角大于15°，歷元間隔5 s，野外觀測(cè)數(shù)據(jù)GPS 自動(dòng)記錄。

共測(cè)GPS基準(zhǔn)站5個(gè)，基線12條。網(wǎng)中平均基線長(zhǎng)度為12.033 km。基線解為雙差固定解。重復(fù)基線最大誤差ds=6 mm，小于規(guī)范允許值17 mm［4］。

GPS控制網(wǎng)二維約束平差控制網(wǎng)使用4個(gè)已知點(diǎn)強(qiáng)制約束解算；高程采用4個(gè)已知高程點(diǎn)參數(shù)擬合。

GPS 靜態(tài)測(cè)量和RTK 測(cè)量使用儀器為天寶r4-3 型GPS 雙頻接收機(jī)，儀器標(biāo)稱精度：靜態(tài)平面±3 mm+0.5 mm/km，靜態(tài)高程±5 mm+0.5 mm/km；動(dòng)態(tài)平面±0.8 cm+1 mm/km，動(dòng)態(tài)高程±1.5 cm+1 mm/km。其擬合解算結(jié)果見表1。

表1 馬店區(qū)塊約束平差、高程擬合解算結(jié)果Table 1.Results of constrained adjustment and elevation fitting for the Madian block

由表1 可知，水平殘差小于規(guī)范要求的20 mm,垂直殘差小于規(guī)范要求的40 mm。

3.2 支子湖區(qū)塊

使用儀器及施測(cè)方法同上。收集C 級(jí)控制點(diǎn)4個(gè)，加密控制點(diǎn)2 個(gè)。共測(cè)GPS 基準(zhǔn)站6 個(gè)，基線18條。網(wǎng)中平均基線長(zhǎng)度為12.194 km?；€解為雙差固定解。重復(fù)基線最大誤差ds=15 mm，小于規(guī)范允許值28 mm。GPS控制網(wǎng)二維約束平差控制網(wǎng)使用4個(gè)已知點(diǎn)強(qiáng)制約束解算；高程采用4 個(gè)已知高程點(diǎn)參數(shù)擬合，其擬合解算結(jié)果見表2。

表2 支子湖區(qū)塊約束平差、高程擬合解算結(jié)果Table 2.Results of constrained adjustment and elevation fitting for the Zhizihu block

由表2可知，水平殘差小于規(guī)范要求的20 mm，垂直殘差小于規(guī)范要求的40 mm。

3.3 劉圩與立倉(cāng)區(qū)塊

使用儀器為廣州市三鼎光電儀器有限公司產(chǎn)T66型衛(wèi)星定位接收機(jī)5 臺(tái)，全部?jī)x器均經(jīng)過法定計(jì)量檢測(cè)檢定合格，在作業(yè)期間，所有參與作業(yè)的儀器均在有效檢定期內(nèi)。

由于立倉(cāng)工區(qū)范圍已布設(shè)有E 級(jí)GPS 控制網(wǎng)（原西賈莊-羅集工區(qū)），且劉圩工區(qū)與立倉(cāng)工區(qū)距離較近，故本次在原西賈莊-羅集工區(qū)E 級(jí)GPS 控制網(wǎng)基礎(chǔ)上，對(duì)原控制網(wǎng)向東擴(kuò)展并重新進(jìn)行平差以取得適合本次作業(yè)的參數(shù)，網(wǎng)中包括擬作為平面起算的C 級(jí)GPS 控制點(diǎn)6 個(gè)，原E 級(jí)GPS 控制點(diǎn)13 個(gè)，新設(shè)GPS控制點(diǎn)4 個(gè)（K100、K200、K300、K400）。新增點(diǎn)均選擇在視野開闊、交通便利、地基相對(duì)穩(wěn)固的位置，點(diǎn)位周圍無大功率電力及無線信號(hào)干擾且無大面積水域，所有新增點(diǎn)均埋設(shè)了普通混凝土標(biāo)石。

作業(yè)時(shí)，以圖形強(qiáng)度較好的大地四邊形作為主要基本圖形，使用4 臺(tái)套GPS 接收機(jī)進(jìn)行同步觀測(cè)；以單三角形作為輔助基本圖形，使用3 臺(tái)套GPS 接收機(jī)進(jìn)行同步觀測(cè)。觀測(cè)過程中，高度截止角均大于15°，同步觀測(cè)衛(wèi)星數(shù)均大于4 顆，同步觀測(cè)時(shí)間均大于40 min。GPS 網(wǎng)平差計(jì)算采用南方測(cè)繪儀器有限公司的GNSS 靜態(tài)解算軟件（版本號(hào)4.05.150324）。本次GPS 靜態(tài)測(cè)網(wǎng)全網(wǎng)共有23 個(gè)測(cè)站點(diǎn)，按照規(guī)范中觀測(cè)時(shí)段數(shù)≥1.6的要求，應(yīng)觀測(cè)37個(gè)測(cè)站數(shù)據(jù)，實(shí)際參與組網(wǎng)平差的測(cè)站數(shù)據(jù)共51個(gè)，觀測(cè)數(shù)量滿足規(guī)范要求。

全網(wǎng)最大邊長(zhǎng)為35 350 m，最小邊長(zhǎng)為3 426 m，平均邊長(zhǎng)為10 515 m。按照儀器精度計(jì)算的整網(wǎng)基線測(cè)量標(biāo)準(zhǔn)差σ＝10.6 mm,在此基礎(chǔ)上放大一倍的高程網(wǎng)基線測(cè)量標(biāo)準(zhǔn)差σ＝21.2 mm。標(biāo)準(zhǔn)差計(jì)算公式如下：

式中：σ為標(biāo)準(zhǔn)差，mm；a為固定誤差，mm；b為比例誤差系數(shù)，10-6mm；d為基線平均邊長(zhǎng)，km。

全網(wǎng)組網(wǎng)基線共75 條，解算模式選擇雙差固定解，解算結(jié)果均滿足要求。重復(fù)基線共14 條，基線長(zhǎng)度較差最大值為13.5 mm，對(duì)應(yīng)重復(fù)基線長(zhǎng)度較差允許值為42.9 mm。

按照3 層深度搜索模式得到三邊閉合環(huán)總數(shù)122個(gè)，其中同步環(huán)49 個(gè)，異步環(huán)73 個(gè)，環(huán)內(nèi)坐標(biāo)分量閉合差最大為8.6 mm，對(duì)應(yīng)環(huán)內(nèi)坐標(biāo)分量允許閉合差為16.7 mm；環(huán)閉合差最大為19.9 mm，對(duì)應(yīng)環(huán)允許閉合差為28.9 mm。

無約束平差中，基線分量改正數(shù)最大值為34.9 mm，對(duì)應(yīng)基線分量改正數(shù)允許值為163.1 mm。約束平差中，基線分量改正數(shù)最大值為43.5 mm，對(duì)應(yīng)基線分量改正數(shù)允許值為163 mm。平面最弱點(diǎn)為WK13 點(diǎn)，點(diǎn)位精度為±5.49 mm，高程擬合采用二次曲面擬合方式，內(nèi)符合精度中誤差為±59.1 mm；最弱點(diǎn)為K400，高程擬合精度為±9.4 mm。以上所有結(jié)果均滿足規(guī)范相應(yīng)要求。

GPS 靜態(tài)測(cè)量和RTK 測(cè)量使用儀器為南方T66型衛(wèi)星定位接收機(jī)GPS雙頻接收機(jī)，儀器標(biāo)稱精度：靜態(tài)平面為±3 mm+1 mm/km，靜態(tài)高程為±5 mm+1 mm/km；動(dòng)態(tài)平面為±1cm+1mm/km，動(dòng)態(tài)高程為±2cm+1mm/km。

4 碎部測(cè)量

控制測(cè)量完成后，馬店、劉圩、支子湖三個(gè)區(qū)塊所需7 個(gè)參數(shù)分別都帶入RTK 中，每個(gè)區(qū)塊選取基站點(diǎn)后，用RTK手簿自帶軟件，做控制點(diǎn)測(cè)量，從最近已知點(diǎn)引點(diǎn)至基站位置，分別從2個(gè)已知點(diǎn)引點(diǎn)至基站，選取得到的2 個(gè)基站坐標(biāo)平均數(shù)，作為最終基站坐標(biāo)。最終坐標(biāo)合格率都為100%。針對(duì)1∶25 000 面積性重力測(cè)量高程中誤差≤±0.2 m的精度要求，采用RTK單基站模式進(jìn)行設(shè)計(jì)點(diǎn)位的點(diǎn)放樣。為方便測(cè)量后質(zhì)檢，采用竹簽定點(diǎn)，懸掛小紙旗或紅布條作為臨時(shí)標(biāo)記。

GPS 測(cè)量基準(zhǔn)站的架設(shè)有專人看守，負(fù)責(zé)維護(hù)基準(zhǔn)站，保證儀器狀態(tài)良好，確保各流動(dòng)站全天正常工作。流動(dòng)站RTK 測(cè)點(diǎn)觀測(cè)衛(wèi)星數(shù)不少于4 顆，PDOP值小于6，高度截止角大于15°，移動(dòng)站與基準(zhǔn)站距離一般不大于10 km。移動(dòng)站和基站數(shù)據(jù)鏈采用中國(guó)移動(dòng)的GPRS網(wǎng)絡(luò)模式。

工作中若遇密集建筑物、大片水域或其他不利地形等情況，一般在1/2點(diǎn)距范圍內(nèi)可作適當(dāng)位移，超出此范圍作棄點(diǎn)處理。

點(diǎn)位、高程、重力測(cè)點(diǎn)質(zhì)量檢查采用“一同三不同”的方法隨野外生產(chǎn)進(jìn)度同步展開，質(zhì)量檢查點(diǎn)在時(shí)間上和空間上分布基本均勻。各項(xiàng)指標(biāo)均優(yōu)于設(shè)計(jì)精度。各工區(qū)統(tǒng)計(jì)結(jié)果見表3。

表3 重力測(cè)點(diǎn)質(zhì)量檢查統(tǒng)計(jì)Table 3.Statistics of gravity survey point quality inspection

全區(qū)統(tǒng)計(jì)測(cè)地平面中誤差為±0.057 m，優(yōu)于設(shè)計(jì)10 m 要求；高程中誤差為±0.025 m，優(yōu)于設(shè)計(jì)0.2 m 要求。由表3 可知，平面和高程中誤差均小于設(shè)計(jì)2 倍中誤差的要求，質(zhì)量檢查率小于設(shè)計(jì)3%。

RTK質(zhì)檢所采用公式［5］分別為：

（1）RTK定位平面中誤差εg的計(jì)算公式如下：

式中：εxi，εyi為第i點(diǎn)x和y方向原始觀測(cè)與檢查觀測(cè)之差；i為檢查點(diǎn)號(hào)，i=1，2，…，n。

（2）RTK定位高程中誤差εh的計(jì)算公式如下：

式中：δhi為第i點(diǎn)高程原始觀測(cè)與檢查觀測(cè)之差；n為檢查點(diǎn)數(shù)，i=1，2，…，n。

針對(duì)在馬店、劉圩、立倉(cāng)三個(gè)區(qū)塊部署的1∶10 000面積性地面高精度磁法測(cè)量，采用帶入?yún)?shù)后的手持GPS 進(jìn)行測(cè)點(diǎn)放樣［5］。首先計(jì)算手持GPS 測(cè)量所需參數(shù)，用WGS84坐標(biāo)系X、Y、Z、A、F值，減去我國(guó)坐標(biāo)系的對(duì)應(yīng)值，得出實(shí)現(xiàn)坐標(biāo)系轉(zhuǎn)換的5個(gè)參數(shù)，即DX、DY、DH、DA、DF，所需WGS84 坐標(biāo)和平差后的目標(biāo)坐標(biāo)由E級(jí)GPS控制測(cè)量得到。

參數(shù)驗(yàn)證：將求得的5 個(gè)參數(shù)輸入手持GPS 中，在應(yīng)用區(qū)域內(nèi)選擇2 個(gè)以上控制點(diǎn)進(jìn)行實(shí)測(cè)，所測(cè)坐標(biāo)與已知坐標(biāo)比對(duì)，如果最大誤差不超過15 m，平均誤差不超過10 m，則計(jì)算參數(shù)可以應(yīng)用。實(shí)際測(cè)量誤差最大為5 m，平均誤差值為4 m滿足以上條件。

手持GPS 質(zhì)量檢查也采用同精度“一同三不同”方式進(jìn)行，隨工作進(jìn)度隨機(jī)進(jìn)行，統(tǒng)計(jì)結(jié)果見表4。由表4 可知，三個(gè)工作區(qū)塊質(zhì)檢檢查率最小的馬店也大于設(shè)計(jì)3%的要求，平面和高程中誤差均小于設(shè)計(jì)精度。

表4 磁法測(cè)量質(zhì)量檢查統(tǒng)計(jì)Table 4.Statistics of magnetic survey quality inspection

5 結(jié)束語(yǔ)

本次測(cè)量工作根據(jù)設(shè)計(jì)和規(guī)范要求結(jié)合物探測(cè)量的特點(diǎn)，有針對(duì)性地安排了GPS控制測(cè)量和碎部測(cè)量，從精度上來說完全滿足設(shè)計(jì)要求，控制測(cè)量點(diǎn)的數(shù)量和密度滿足設(shè)計(jì)規(guī)范要求，分布完全覆蓋3 個(gè)測(cè)區(qū)，但是由于GPS數(shù)據(jù)具有偶然性和信號(hào)易受各種因素干擾，所以每次測(cè)量觀測(cè)3次取平均數(shù)為最終結(jié)果。RTK 單基準(zhǔn)站模式下，流動(dòng)站坐標(biāo)是由一個(gè)基準(zhǔn)站來確定的，因此可靠性較差，所以每次測(cè)量完成后進(jìn)行質(zhì)檢時(shí)，質(zhì)檢會(huì)選擇另一個(gè)基站（坐標(biāo)已知），通過質(zhì)檢測(cè)量數(shù)據(jù)表可知誤差符合規(guī)范要求。手持GPS 測(cè)量主要是注意避開影響信號(hào)的障礙物和大范圍水面等因素?？偟膩碚f，GPS 測(cè)量在工作中操作簡(jiǎn)單，效率高。相信通過不斷總結(jié)經(jīng)驗(yàn)和技術(shù)改進(jìn)，GPS 技術(shù)一定能更好地為物探測(cè)量提供服務(wù)。