付 青

(武漢大學(xué)遙感信息工程學(xué)院,湖北 武漢 430079)

0 引 言

精密單點定位(PPP)是利用精密衛(wèi)星軌道產(chǎn)品和高精度的誤差模型消除GNSS接收機觀測誤差,并采用最小二乘平差或濾波推估方法,實現(xiàn)單臺GNSS接收機的高精度定位技術(shù)[1]。PPP技術(shù)具有全天候、無需參考站和高精度定位的優(yōu)點,一直為科研項目和民用工程提供高可靠、高精度的保障,譬如準(zhǔn)實時GPS的氣象研究、地殼形變研究等[2-3]。但其必須有一個較長時間的收斂過程才能得到高精度穩(wěn)定解,正是該特點阻礙它在很多領(lǐng)域的進一步發(fā)展。

隨著全球GNSS系統(tǒng)的不斷完善,特別是北斗衛(wèi)星導(dǎo)航系統(tǒng)于2012年12月27日正式服務(wù)亞太地區(qū),研究GNSS組合定位再次成為很多學(xué)者的研究熱點。截止目前為止,BDS的在軌工作衛(wèi)星由5 GEO衛(wèi)星+5 IGSO衛(wèi)星+4 MEO衛(wèi)星組成?，F(xiàn)在對BDS關(guān)注熱點集中在單系統(tǒng)PPP定位精度、BDS與GPS組合定位精度和BDS與GPS、GLONASS三星組合定位精度等方面。YANG[4]等分析了BDS公眾用戶的可視衛(wèi)星數(shù)及其幾何精度值(DOP)。施闖[5]等利用武漢大學(xué)GNSS中心提供的精密星歷和精密鐘差產(chǎn)品解算得到BDS靜態(tài)PPP達到cm級,動態(tài)RTK可以達到5～10 cm.李偉[6]等分析了GPS單系統(tǒng)、BDS單系統(tǒng)、GPS與BDS組合系統(tǒng)三種定位模式,使用三個站兩天的觀測數(shù)據(jù)實驗,得出GPS與BDS組合定位系統(tǒng)的收斂時間相比單系統(tǒng)明顯減少,而在單系統(tǒng)本身衛(wèi)星星座分布良好時,組合定位系統(tǒng)定位精度提高不大。這些實驗都在單BDS系統(tǒng)或者BDS與GPS雙星組合系統(tǒng)進行,對于BDS與GPS、GLONASS三星組合系統(tǒng)研究非常少。本文旨在分析三星組合系統(tǒng)的收斂時間變化以及其定位精度。由于多星組合系統(tǒng)必須在一定的時間和坐標(biāo)系統(tǒng)進行,為此,本文統(tǒng)一了三星組合定位的時間系統(tǒng)和空間系統(tǒng),利用武漢九峰站全天觀測數(shù)據(jù)進行實驗,驗證三星組合定位能有效縮短PPP收斂時間以及提高在單GNSS系統(tǒng)衛(wèi)星數(shù)較少或者衛(wèi)星星座分布較差時的定位精度。

1 時空基準(zhǔn)統(tǒng)一

1.1 時間基準(zhǔn)統(tǒng)一

GPS的時間系統(tǒng)屬于原子時系統(tǒng),它的秒長和原子時秒長相同,但它和國際原子時存在著不同的原點,并且在任一時刻它們二者之間都存在一個常量偏差(19 s)。

北斗系統(tǒng)的時間基準(zhǔn)為北斗時(BDT)。BDT采用國際單位制(SI)秒為基本單位連續(xù)累計,不閏秒,起始歷元為2006年1月1日協(xié)調(diào)世界時(UTC)00時00分00秒,采用周和周內(nèi)秒計數(shù)。BDT與UTC 之間的閏秒信息在導(dǎo)航電文中傳播[7]。

GLONASS時間系統(tǒng)(GLONASST)屬于UTC時間系統(tǒng),它是基于GLONASS同步中心(CS)時間產(chǎn)生的,由于GLONASS控制部分本身存在的特性,使得GLONASST與俄羅斯維持的協(xié)調(diào)世界時UTC(SU)存在的整數(shù)差為3 h,此外它們之間還存在有1 ms 以內(nèi)的系統(tǒng)誤差τr[8].

三類衛(wèi)星的時間系統(tǒng)詳細描述如表1所示。

表1 三類衛(wèi)星時間系統(tǒng)的比較

1.2 坐標(biāo)基準(zhǔn)統(tǒng)一

GPS坐標(biāo)系統(tǒng)采用的是WGS-84坐標(biāo)系。該坐標(biāo)系是地心空間直角坐標(biāo)系,它的原點為地球質(zhì)心,Z軸指向國際時間服務(wù)機構(gòu)(BIH)1984.0定義的協(xié)議地球極(CTP)方向,X軸指向BIH 1984.0的零子午面和CTP赤道的交點,Y軸與Z軸、X軸垂直構(gòu)成右手坐標(biāo)系。BDS坐標(biāo)系統(tǒng)采用的是我國2000大地坐標(biāo)系統(tǒng)(CGCS2000),它的定義和國際地球參考系統(tǒng)(ITRS)相一致[7]。GLONASS坐標(biāo)系統(tǒng)以前采用的蘇聯(lián)1985地心坐標(biāo)系(SGS85),1993年以后改為使用PZ90.02[8-9]坐標(biāo)系,該坐標(biāo)系是由俄羅斯進行地面網(wǎng)與空間網(wǎng)聯(lián)合平差后所建立的。WGS-84、CGCS2000和PZ90.02坐標(biāo)系的參考橢球定義的基本參數(shù)如表2所示。

表2 GPS、BDS與GLONASS坐標(biāo)系統(tǒng)參數(shù)

GPS與BDS兩個坐標(biāo)系的主要不同在于參考橢球的扁率f有微小的差異,而這種差異在當(dāng)前的測量精度水平中可以忽略,對于在坐標(biāo)系定義上的比較,可以認為在同一歷元下CGCS2000和WGS-84在坐標(biāo)系的實現(xiàn)精度范圍內(nèi),兩者的坐標(biāo)是一致的[9]。

當(dāng)前,大多數(shù)學(xué)者認為使用俄羅斯MCC計算得到的七參數(shù)作為二個坐標(biāo)系轉(zhuǎn)換參數(shù)的精度最高,該七參數(shù)是通過全球激光跟蹤測軌數(shù)據(jù)計算得到,PZ90.02與WGS84之間的坐標(biāo)轉(zhuǎn)換公式[10]為

(1)

2 多星組合定位模型

2.1 原始觀測量數(shù)學(xué)模型

在多衛(wèi)星組合定位中,由GNSS接收機獲取的偽距和載波相位觀測量可以給出二者的觀測方程[11]。其中,g表示GPS衛(wèi)星,l表示GLONASS衛(wèi)星,c表示BDS衛(wèi)星。

(2)

式中:下標(biāo)j(=1,2)為衛(wèi)星頻率號;P為偽距觀測量;Φ為載波相位觀測量;dtr、dts分別為接收機與衛(wèi)星鐘差;c為光速;dtrop為對流層延遲量;dion為電離層延遲量;λ為載波波長;N為整周模糊度;ePj、εΦj分別為偽距與載波相位的觀測噪聲；ρ為信號發(fā)射時刻的衛(wèi)星與信號接收時刻的GNSS接收機天線相位中心的幾何距離。

2.2 PPP非差組合模型

通常情況,采用無電離層組合模型消除偽距與載波相位中的電離層延遲影響,其表達式為

(3)

2.3 擴展卡爾曼濾波

由于在組合定位中,觀測方程是非線性的,故本文采用擴展卡爾曼濾波(EKF)[12]。其表達式為

(4)

當(dāng)觀測方程經(jīng)過線性化后,狀態(tài)量與其相關(guān)系數(shù)矩陣的時間更新為

(5)

3 實驗及其結(jié)果分析

實驗數(shù)據(jù)來源于武漢九峰站2013年8月28日全天觀測數(shù)據(jù),該站接收機型號為TRIMBLE NETR9,它可以獲取當(dāng)前所有的在軌導(dǎo)航衛(wèi)星系統(tǒng)的觀測數(shù)據(jù),并且設(shè)置它的觀測時間為GPST,衛(wèi)星高度截止角為15°.該數(shù)據(jù)采樣間隔30 s,全天歷元總數(shù)為2 880.對全天實測數(shù)據(jù)按照如下6種方案進行實驗:①GPS單系統(tǒng);②BDS單系統(tǒng);③GLONASS單系統(tǒng);④GPS與BDS組合雙系統(tǒng);⑤GPS與GLONASS組合雙系統(tǒng);⑥GPS、BDS與GLONASS組合三系統(tǒng)。

圖1示出了六種方案中可視衛(wèi)星數(shù)，從圖1可以看出,全天時間中GLONASS可視衛(wèi)星數(shù)最少,GPS可視衛(wèi)星數(shù)最多,BDS處于二者之間。表3示出了6種方案中可視衛(wèi)星個數(shù)情況,從該表可以看出全天時間中GPS平均可視衛(wèi)星顆數(shù)7.7顆;BDS平均可視衛(wèi)星數(shù)7.1顆;GLONASS平均可視衛(wèi)星數(shù)6.0顆;三星組合定位的平均可視衛(wèi)星顆數(shù)最高,達到20.4顆。可視衛(wèi)星顆數(shù)越多,系統(tǒng)的定位精度亦更高。

圖1 可視衛(wèi)星數(shù)曲線

表3可視衛(wèi)星數(shù)

定位方案可視衛(wèi)星數(shù)(顆) 最大值最小值平均值 GPS1247.7 BDS867.1 GLONASS956.0 GPS/BDS19714.5 GPS/GLONASS191013.7 GPS/BDS/GLONASS261320.4

圖2 6種方案X、Y、Z三個坐標(biāo)偏差圖

圖2示出了6種方案X、Y、Z三個方向的偏差，從圖2可以看出,GPS單系統(tǒng)定位性能較好,收斂后全天X、Y、Z三個方向都處于10 cm以內(nèi),主要是因為GPS系統(tǒng)建設(shè)較早,系統(tǒng)穩(wěn)定,衛(wèi)星已實現(xiàn)全球覆蓋,星座空間分布較好。BDS單系統(tǒng)定位性能次之,收斂后全天X、Y、Z三個方向都處于20 cm以內(nèi),這是因為BDS系統(tǒng)建設(shè)正處于第二個階段,衛(wèi)星覆蓋亞太地區(qū),暫未實現(xiàn)全球覆蓋。GLONASS單系統(tǒng)收斂后全天X、Y、Z三個方向處于40 cm以內(nèi),其定位精度受GLONASS衛(wèi)星本身硬件延遲以及精密衛(wèi)星軌道精度影響較大。從圖2中還可以看出,GPS單系統(tǒng)收斂時間為180 min左右(采樣頻率30 s),BDS單系統(tǒng)收斂時間為200 min左右,GLONASS單系統(tǒng)收斂時間為240 min左右,而對于GPS與BDS組合定位的收斂時間為150 min左右,GPS,BDS和GLONASS三星組合定位的收斂時間只需要90 min左右,明顯縮短單系統(tǒng)的收斂時間。所以，對于GPS與BDS雙系統(tǒng)組合定位,能夠有效的提高BDS單系統(tǒng)定位性能,有效的減少BDS單系統(tǒng)的收斂時間,該組合系統(tǒng)的定位精度優(yōu)于5 cm,同樣對于GLONASS與GPS組合定位,有效地提高了單GLONASS系統(tǒng)的定位精度。對于三星組合定位系統(tǒng),與單系統(tǒng)比較,收斂時間最短,定位精度最高。

從表4可以看出,GPS單系統(tǒng)的X、Y、Z三個方向平均偏差分別為8.81 cm、6.10 cm和5.02 cm,BDS單系統(tǒng)的X、Y、Z三個方向平均偏差分別為8.90 cm、9.53 cm和5.39 cm,GLONASS單系統(tǒng)的X、Y、Z三個方向平均偏差分別為2.95 cm、18.62 cm和5.04 cm,三星組合系統(tǒng)的定位精度最高,X、Y、Z三個方向平均偏差分別為8.00 cm、1.40 cm和0.37 cm.三星組合定位精度的提高,主要是由于組合定位系統(tǒng)的空間可視衛(wèi)星明顯增加如表3所示,改善了單系統(tǒng)的衛(wèi)星空間分布,尤其是在單系統(tǒng)出現(xiàn)遮擋,可視衛(wèi)星較少,組合定位系統(tǒng)的精度要明顯高于單系統(tǒng)。

表4 統(tǒng)計6種方案X、Y、Z三個坐標(biāo)偏差

4 結(jié)束語

隨著我國北斗衛(wèi)星導(dǎo)航系統(tǒng)加入GNSS行列,研究BDS與其他系統(tǒng)之間的組合定位尤為重要,特別是在單系統(tǒng)受所處環(huán)境影響導(dǎo)致可視衛(wèi)星顆數(shù)較少時,組合系統(tǒng)能夠有效的提高單系統(tǒng)的定位精度。本文在這個背景下,對BDS、GPS與GLONASS三類衛(wèi)星的單系統(tǒng)定位、雙系統(tǒng)組合定位以及三星組合定位6種情況進行實驗分析,得到以下結(jié)論:

1)組合定位系統(tǒng)收斂時間相比單系統(tǒng)收斂時間有所減少,這對于改善單系統(tǒng)PPP需要長時間收斂時間具有重要意義。

2)組合定位系統(tǒng)能夠有效的提高單系統(tǒng)的可視衛(wèi)星顆數(shù),改善組合系統(tǒng)的空間衛(wèi)星分布,這也是組合定位系統(tǒng)能夠提高單系統(tǒng)定位精度的前提保證。特別是在環(huán)境惡劣情況下,單系統(tǒng)衛(wèi)星顆數(shù)較少時,組合定位系統(tǒng)能夠發(fā)揮明顯優(yōu)勢。同時,GNSS組合定位必將是未來定位技術(shù)的發(fā)展趨勢,本文為GNSS定位提供一定的參考意義。

[1]KOUBA J,HEROUS H. Precise point positioning using IGS orbit and clock products [J]. GPS Solutions,2011,5(2):12-28.

[2]GENDT G,DICK G,REIGBER C,etal. Near real time GPS water vapor monitoring for numerical weather prediction in Germany [J]. Meteorological Society Japan,2004,82(1B):361-370.

[3]CALAIS E,HAN J Y,DEMETS C,etal. Deformation of the North American plate interior from a decade of continuous GPS measurements [J]. Geophysical Reserach,2006,111,B06402,doi:10.1029/2005JB004253.

[4]YANG Yuanxi,LI Jinlong,XU Junyi,etal. Contribution of the compass satellite navigation system to global PNT users [J]. Chinese Science Bulletin,2011,56(26):2813-2819.

[5]SHI Chung,ZHAO Qile,LI Min,etal. Precise orbit determination of Beidou satellites with precise positioning [J]. Science China: Earth Sciences,2012,55(7):1079-1086.

[6]LI Wei,TEUNISSEN P J G,ZHANG Baocheng,etal. Precise point positioning using GPS and Compass observations [C]//Proceedings of the 4th China satellite navigation conference (CSNC),Wuhan,China,15-17 May,2013.

[7]中國衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室.北斗衛(wèi)星導(dǎo)航系統(tǒng)空間信號接口控制文件公開服務(wù)信號B1 (2.0版)[R]. 2013.

[8]Russian Institute of Space Device Engineering. GLONASS interface control document(ICD)(version5.1)[R]. 2008.

[9]黨亞民，成英燕，薛樹強.大地坐標(biāo)系統(tǒng)及其應(yīng)用[M].北京：測繪出版社,2010.

[10]李鶴峰，黨亞民，秘金鐘，等.BDS與GPS、GLONASS多模整合導(dǎo)航定位時空統(tǒng)一[J].大地測量與地球動力學(xué),2013,33(4):73-78.

[11]韓保民，歐吉坤.基于GPS非差觀測值進行精密單點定位研究[J].武漢大學(xué)學(xué)報·信息科學(xué)版,2003,28(4):409-412.

[12]ARTHUR G. Applied optimal estimation [M]. The M. I. T. Press,1974.

[13]IGS. JFNG site information form international GNSS service [EB/OL].ftp://igscb.jpl.nasa.gov/pub/station/general/sitelog_instr.txt,2013.