周建營(yíng),陳國(guó)恒,陳梅森

(廣東省國(guó)土資源測(cè)繪院,廣東 廣州 510500)

0 引 言

隨著衛(wèi)星導(dǎo)航定位技術(shù)的迅速發(fā)展,我國(guó)在測(cè)繪地理信息產(chǎn)業(yè)的工程領(lǐng)域已進(jìn)入集全球定位系統(tǒng)(GPS)、格洛納斯衛(wèi)星導(dǎo)航系統(tǒng)(GLONASS)、北斗衛(wèi)星導(dǎo)航系統(tǒng)(BDS)和伽利略衛(wèi)星導(dǎo)航系統(tǒng)(Galileo)的全球衛(wèi)星導(dǎo)航系統(tǒng)(GNSS)時(shí)代,其中BDS的建設(shè)發(fā)展起到了重要作用. 在2012年,我國(guó)完成了北斗二號(hào)系統(tǒng)(BDS-2)5顆地球靜止軌道(GEO)衛(wèi)星、5顆傾斜地球同步軌道(IGSO)衛(wèi)星、4顆中圓地球軌道(MEO)衛(wèi)星的發(fā)射組網(wǎng),這種星座在低緯度地區(qū)及林區(qū)、城市交接區(qū)、山川峽谷區(qū)性能突出[1]. 按照目前BDS“三步走”的發(fā)展戰(zhàn)略,2020年將建成北斗三號(hào)系統(tǒng)(BDS-3),且在衛(wèi)星星座、信號(hào)設(shè)計(jì)等方面相較于BDS-2均有了很大改進(jìn)[2],這將為測(cè)繪生產(chǎn)工程控制網(wǎng)的建立提供多一種衛(wèi)星系統(tǒng)的選擇,而不需僅依靠GPS單星系統(tǒng). 部分文獻(xiàn)利用香港連續(xù)運(yùn)行參考站(CORS)或Multi-GNSS Experiment(MGEX)項(xiàng)目的觀測(cè)數(shù)據(jù)進(jìn)行BDS解算分析[3-5],一定程度上為BDS工程控制網(wǎng)數(shù)據(jù)處理提供參考. 而對(duì)于實(shí)際測(cè)繪生產(chǎn)工程控制網(wǎng),各測(cè)站的地理環(huán)境復(fù)雜,觀測(cè)條件相對(duì)較差. 文獻(xiàn)[6]以新疆地區(qū)某水利工程GNSS框架網(wǎng)7個(gè)控制點(diǎn)為例,比較分析了BDS和GPS的數(shù)據(jù)處理結(jié)果,試驗(yàn)說(shuō)明BDS可滿足該網(wǎng)的建設(shè)需要,但綜合比較后還是以GPS作為最終應(yīng)用成果;文獻(xiàn)[7]對(duì)高速鐵路基礎(chǔ)平面控制網(wǎng)(CPI)進(jìn)行BDS/GPS單獨(dú)系統(tǒng)的解算,結(jié)果表明BDS各項(xiàng)指標(biāo)的精度可用于高速鐵路測(cè)量控制網(wǎng)的建立.

本文以廣東省域內(nèi)某工程控制網(wǎng)為例,通過(guò)BDS/GPS單星系統(tǒng)的數(shù)據(jù)處理,分別從觀測(cè)數(shù)據(jù)質(zhì)量情況、基線解算與平差精度、坐標(biāo)結(jié)果差異性等指標(biāo),分析比較BDS/GPS單星系統(tǒng)下的建立情況,驗(yàn)證了BDS單星系統(tǒng)建立工程控制網(wǎng)的可行性.

1 工程控制網(wǎng)情況

本文選取的控制網(wǎng)位于廣東沿海地區(qū),該網(wǎng)按GNSS-B級(jí)的級(jí)別要求建立[8]. 工程區(qū)域布設(shè)了26個(gè)控制點(diǎn),以邊連式的觀測(cè)方式分成兩個(gè)同步環(huán),觀測(cè)時(shí)間為2019年第208-210天和第214-216天,接收機(jī)同步采集BDS/GPS數(shù)據(jù),同時(shí)聯(lián)合工程區(qū)域周邊19個(gè)GDCORS站一并進(jìn)行數(shù)據(jù)處理. 控制網(wǎng)點(diǎn)位布設(shè)如圖1所示.

圖1 控制網(wǎng)布設(shè)示意圖

為了客觀評(píng)價(jià)BDS/GPS單星系統(tǒng)的數(shù)據(jù)質(zhì)量情況、基線解算與平差精度、坐標(biāo)結(jié)果差異性等各項(xiàng)指標(biāo),本次實(shí)例處理流程均保持一致. 針對(duì)BDS/GPS多星觀測(cè)數(shù)據(jù)的rinex3.0格式采用G-Nut/Anubis和RTKLIB進(jìn)行數(shù)據(jù)質(zhì)量情況分析[9-10];廣播星歷和精密星歷均采用混合星歷產(chǎn)品;數(shù)據(jù)處理采用GAMIT/GLOBK 10.7,基線解算sittbl.表中JMGT站設(shè)置為0.05,0.05,0.05,歷元間隔為30 s,衛(wèi)星截至高度角為10°,采用對(duì)流層延遲估計(jì),海潮模型采用FES2004,其它參數(shù)均為默認(rèn)模式[11]. 平差計(jì)算以JMGT、FKGT、ZQGT、YJGT、GTCH、ZHGT共6個(gè)GDCORS站作為CGCS2000坐標(biāo)起算基準(zhǔn).

2 數(shù)據(jù)分析與比較

2.1 數(shù)據(jù)質(zhì)量情況

首先以其中一個(gè)測(cè)站B109為例,輸出全天各時(shí)段的可觀測(cè)衛(wèi)星情況如圖2所示.

圖2 可觀測(cè)衛(wèi)星情況

從圖2來(lái)看,測(cè)站能接收到的GPS衛(wèi)星(代碼為G)為32顆,BDS衛(wèi)星(代碼為C)為30顆,其中BDS同時(shí)接收了BDS-2和BDS-3衛(wèi)星,可觀測(cè)衛(wèi)星數(shù)BDS與GPS相差不大.

再對(duì)所有觀測(cè)站點(diǎn)進(jìn)行數(shù)據(jù)質(zhì)量檢查,分別從衛(wèi)星截至高度角>10°時(shí)的預(yù)期觀測(cè)值數(shù)量(Exp>10)、衛(wèi)星截至高度角>10°時(shí)的實(shí)際觀測(cè)值數(shù)量(Hav>10)、數(shù)據(jù)利用率(Ratio)、周跳比情況(o/slps、CSR)、多路徑誤差值(mp1、mp2、mp6、mp7)等各項(xiàng)指標(biāo)進(jìn)行統(tǒng)計(jì). 考慮各測(cè)站均為連續(xù)觀測(cè),此處僅列出26個(gè)控制點(diǎn)測(cè)站其中一天的質(zhì)量分析情況,如圖3所示.

(a)觀測(cè)值數(shù)量與數(shù)據(jù)利用率

(b)周跳比情況

(c)多路徑誤差值圖3 質(zhì)量檢查情況

由圖3(a)可知,各測(cè)站的數(shù)據(jù)可利用率均大于80%,符合相關(guān)要求. 其中80%～85%的有7個(gè)測(cè)站,最低值為B056,數(shù)據(jù)可利用率80.9%;85%～90%的有15個(gè)測(cè)站;大于90%的有4個(gè)測(cè)站,最高值為B099,數(shù)據(jù)可利用率91.16%. 周跳比情況以CSR值來(lái)評(píng)定[12],由圖3(b)可知,各測(cè)站均小于10,符合相關(guān)要求. 其中相對(duì)較大的有B056(CSR=7.19)、B106(CSR=9.80)、B109(CSR=9.35). 多路徑誤差值(mp)一般認(rèn)為小于50 cm即代表測(cè)站觀測(cè)條件良好,若mp值越小,說(shuō)明多路徑效應(yīng)影響越小[13]. 本次工程控制網(wǎng)BDS對(duì)應(yīng)的多路徑頻率索引分別為mp2、mp6和mp7,其中測(cè)站B072、B113可接收3個(gè)頻率,其余測(cè)站可接收2個(gè),GPS對(duì)應(yīng)的多路徑頻率索引分別為mp1和mp2. 由圖3(c)可知,測(cè)站中BDS的多路徑效應(yīng)影響總體小于GPS,所有測(cè)站的BDS多路徑誤差值均小于50 cm,符合相關(guān)要求,有22個(gè)測(cè)站的GPS多路徑誤差值均小于50 cm,有4個(gè)測(cè)站的GPS多路徑誤差值大于50 cm,分別為B056、B106、B109、B111.

2.2 基線解算與平差精度

根據(jù)上述的數(shù)據(jù)質(zhì)量情況,對(duì)該26個(gè)測(cè)站與周邊19個(gè)GDCORS站進(jìn)行BDS/GPS單星系統(tǒng)數(shù)據(jù)處理. 限于篇幅,本文列出了數(shù)據(jù)質(zhì)量較差的B109測(cè)站,以及數(shù)據(jù)質(zhì)量較好的B099測(cè)站較長(zhǎng)和較短基線解算統(tǒng)計(jì)情況,同時(shí)輸出兩個(gè)測(cè)站的天空及相位殘差圖,如表1和圖4所示.

表1 基線解算統(tǒng)計(jì) m

(a)B109-BDS

(b)B109-GPS

(c)B099-BDS

(d)B099-GPS圖4 天空及相位殘差圖

從表1可知,數(shù)據(jù)質(zhì)量較差的B109測(cè)站在N、E值上的BDS/GPS差值均為毫米級(jí),而在U的差值上均為厘米級(jí),其中較長(zhǎng)基線B109-FKGT的U差值為4.09 cm,較短基線B109-B111的U差值為3.04 cm. 數(shù)據(jù)質(zhì)量較好的B099測(cè)站在N、E、U值上的BDS/GPS無(wú)論是較長(zhǎng)基線或較短基線,其差值均為毫米級(jí),差異性不大. 由圖4可知,B109測(cè)站天空?qǐng)D中的正負(fù)殘差相較于B099測(cè)站更大,表明B109測(cè)站受多路徑或水汽影響較大. 從各時(shí)段的衛(wèi)星運(yùn)動(dòng)情況來(lái)比較,BDS的衛(wèi)星運(yùn)動(dòng)軌跡明顯少于GPS,與可觀測(cè)衛(wèi)星情況不符,從解算結(jié)果文件來(lái)看,原因是GAMIT/GLOBK10.7目前只解算了BDS-2衛(wèi)星. 從相位殘差來(lái)比較,與GPS相比,BDS的相位殘差紅線波動(dòng)較大,點(diǎn)位分布較發(fā)散. B099測(cè)站的相位殘差紅線波動(dòng)小于B109測(cè)站,點(diǎn)位分布更集中,這也證明了B099測(cè)站的數(shù)據(jù)質(zhì)量是優(yōu)于B109測(cè)站的.

針對(duì)BDS/GPS的數(shù)據(jù)質(zhì)量差異,對(duì)第209天所有測(cè)站的N、E、U基線分量進(jìn)行精度比較,如圖5所示.

圖5 各基線分量精度比較

由圖5可知,各測(cè)站基線解算精度比較中,BDS在N、E、U方向的精度均低于GPS,其中BDS的N方向平均精度為4.6 mm,E方向平均精度為6.2 mm,U方向平均精度為32.9 mm;GPS的N方向平均精度為3.7 mm,E方向平均精度為4.4 mm,U方向平均精度為17.2 mm. BDS基線解算精度整體低于GPS.

再以JMGT、FKGT、ZQGT、YJGT、GTCH、ZHGT共6個(gè)GDCORS站作為CGCS2000坐標(biāo)起算基準(zhǔn),合并6天的基線解算結(jié)果進(jìn)行平差計(jì)算,比較N、E、U方向的平差精度,如圖6所示.

圖6 平差精度比較

由圖6可知,平差計(jì)算中BDS的N方向平均精度為1.74 mm,略低于GPS的N方向平均精度1.39 mm;BDS的E方向平均精度為2.28 mm,略低于GPS的E方向平均精度1.61 mm;BDS的U方向平均精度12.09 mm,低于GPS的U方向平均精度6.25 mm. 整體來(lái)看,BDS在N、E、U方向平差精度均低于GPS,其中U方向相對(duì)更顯著.

2.3 坐標(biāo)結(jié)果差異性

為進(jìn)一步比較BDS/GPS單星工程控制網(wǎng)的數(shù)據(jù)處理差異,本文分別提取平差文件中26個(gè)待求測(cè)站的CGCS2000三維坐標(biāo)值,以BDS/GPS坐標(biāo)較差的絕對(duì)值來(lái)比較分析其結(jié)果的差異性,如圖7所示.

圖7 坐標(biāo)結(jié)果差異

由圖7可知,26個(gè)待求測(cè)站中BDS/GPS坐標(biāo)結(jié)果在X方向上平均較差為5.8 mm,較差小于10 mm的有22個(gè),10～20 mm的有3個(gè),大于20 mm的有1個(gè),最大較差為B111測(cè)站的23.42 mm;Y方向上平均較差為13.51 mm,較差小于1 mm的有11個(gè),10～20 mm的有10個(gè),大于20 mm的有5個(gè),最大較差為B111測(cè)站的46.45 mm;Z方向上平均較差為4.41 mm,較差小于10 mm的有23個(gè),10～20 mm的有3個(gè),沒(méi)有大于20 mm的測(cè)站,最大較差為B054的15.2 mm. 從結(jié)果差異性來(lái)看,BDS/GPS的坐標(biāo)結(jié)果在X、Z方向的差異相對(duì)更小,Y方向的差異相對(duì)更大,而且數(shù)據(jù)質(zhì)量較差的測(cè)站,其坐標(biāo)結(jié)果差異性也相對(duì)更大,但無(wú)論數(shù)據(jù)質(zhì)量如何,這種差異性均為毫米級(jí)或者厘米級(jí). 由于只處理了BDS-2衛(wèi)星,BDS-3暫無(wú)法一并進(jìn)行解算,導(dǎo)致BDS的數(shù)據(jù)處理衛(wèi)星數(shù)減少了一半,影響了數(shù)據(jù)處理精度和坐標(biāo)結(jié)果.

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

本文利用廣東省域內(nèi)某工程控制網(wǎng)的應(yīng)用實(shí)例,分析比較了BDS/GPS單星系統(tǒng)下的建立情況,得到了以下幾點(diǎn)結(jié)論:

1)測(cè)站的BDS/GPS數(shù)據(jù)質(zhì)量情況大部分良好,其中數(shù)據(jù)可利用率均大于80%,反映周跳比的CSR值均小于10,所有26個(gè)測(cè)站的BDS和22個(gè)測(cè)站的GPS多路徑誤差均小于50 cm,有4個(gè)測(cè)站的GPS多路徑誤差值大于50 cm,BDS的多路徑效應(yīng)影響要小于GPS,BDS可見(jiàn)衛(wèi)星數(shù)與GPS相差不大.

2)數(shù)據(jù)質(zhì)量較好的測(cè)站,BDS/GPS基線解算的各分量值均為毫米級(jí)的差異;質(zhì)量較差的測(cè)站,BDS/GPS基線解算的N、E分量值為毫米級(jí)的差異,U分量值為厘米級(jí)的差異. 整體的數(shù)據(jù)處理方面,BDS的基線解算與平差精度低于GPS,且在U方向上相對(duì)更為顯著.

3)數(shù)據(jù)質(zhì)量較好的測(cè)站,BDS/GPS得到的坐標(biāo)結(jié)果差異性平均只有毫米級(jí),完全可替代GPS作為工程應(yīng)用的最終解算成果;在數(shù)據(jù)質(zhì)量較差的測(cè)站,Z方向上平均為毫米級(jí)的差異,X、Y方向上平均為厘米級(jí)的差異,可滿足一般工程控制網(wǎng)的要求.

4)由于目前GAMIT/GLOBK 10.7只能處理BDS-2的數(shù)據(jù),導(dǎo)致BDS-3衛(wèi)星數(shù)據(jù)尚未能有效參與解算,隨著B(niǎo)DS-3的全球組網(wǎng)和相關(guān)高精度數(shù)據(jù)解算程序的更新,BDS的解算精度會(huì)進(jìn)一步提高.