程國祥

（安徽理工大學空間信息與測繪工程學院，安徽 淮南 232001）

1 GPS/BDS 實時精密單點定位性能的研究背景

精密單點定位 （Precise Point Position，PPP）技術自提出以來一直是衛(wèi)星導航定位高精度應用的前沿和熱點[1]。目前，通過國際的全球?qū)Ш叫l(wèi)星系統(tǒng) （Global Navigation Satellite System，GNSS）服務機構(gòu)或分析中心提供的精密軌道及鐘差產(chǎn)品，在全球范圍內(nèi)使用單臺接收機就能完成PPP，實現(xiàn)厘米級至米級的定位精度[2-3]。

PPP 所依靠的精密產(chǎn)品往往具有一定的時延，不能滿足PPP 的實時應用。為了滿足對實時PPP（Real Time PPP，RTPPP）的需求，國際的GNSS服務機構(gòu)于2013 年開始提供實時精密產(chǎn)品。目前，多家服務機構(gòu)已提供多系統(tǒng)實時服務，已有一些學者對各服務機構(gòu)提供的實時服務進行了研究，其中大部分是基于 GPS/GLONASS 或 GPS/GLONASS/Galileo 的多系統(tǒng)組合RTPPP，但是其側(cè)重點各有不同[4-6]。

但是，在北斗衛(wèi)星導航系統(tǒng)（BeiDou navigation satellite System，BDS）RTPPP 服務方面，目前僅有法國國家空間研究中心 （Centre National d'études Spatiales，CNES）的實時精密產(chǎn)品CLK93 提供BDS實時狀態(tài)空間表示 （State Space Representation，SSR）改正信息，為BDS 提供實時產(chǎn)品支持，且僅支持BDS-2 衛(wèi)星，包括5 顆GEO 衛(wèi)星、6 顆IGSO衛(wèi)星、3 顆MEO 衛(wèi)星。對于BDS 單系統(tǒng)RTPPP 研究已經(jīng)初步開展，且證明了BDS 對多系統(tǒng)組合PPP性能提升做出的貢獻[7-9]。

為了進一步對 BDS RTPPP 及 GPS/BDS 組合RTPPP 性能進行分析，本文基于GPS 與BDS 進行組合RTPPP 研究，以期為BDS RTPPP 推廣提供參考。本文采用CNES CLK93 播發(fā)的BDS 實時軌道、鐘差產(chǎn)品，對 5 個多 GNSS 實驗 （Multi-GNSS Experiment，MGEX） 測站 2019 年 DOY 098 數(shù)據(jù)，在BDS，GPS，BDS/GPS 3 種模式下進行RTPPP 解算，并從收斂時間和定位精度兩方面對其性能進行分析。

2 實時精密單點定位模型

2.1 觀測方程組

在進行PPP 時，GNSS 基本觀測方程組為

式中：P 和Φ 分別為偽距和載波的觀測值；ρ 為衛(wèi)星和接收機之間的集合距離；c 為光速；r 和s 分別為接收機和衛(wèi)星；dtr和dts分別為接收機和衛(wèi)星鐘差；dr和ds分別為接收機和衛(wèi)星與頻率相關的碼偏差；T 為對流層延遲；ε（P）和ε（Φ）分別為位居和載波觀測值的測量噪聲及未模型化誤差；λ 為對應頻率的波長；br和bs分別為接收機和衛(wèi)星端的相位偏差；N 為模糊度。對于雙頻PPP，目前主要組合方式有消電離層 （Ionosphere Free，IF）組合方式、UofC模型組合方式和非組合模型組合方式。選取IF組合方式消除一階電離層的影響[10-14]，表達式方程組為

式中：PIF和ΦIF分別為IF 組合后的偽距和載波的觀測值；dr，IF和分別為IF 組合后的接收機和衛(wèi)星端的碼偏差；ε（PIF）和 ε（ΦIF）為 IF 組合的測量噪聲及未模型化誤差；λIF為IF 組合后的頻率波長；br，IF和為IF 組合后的接收機和衛(wèi)星的相位偏差；NIF為IF 組合后的模糊度；其余參數(shù)同式 （1）。

2.2 參數(shù)估計及誤差處理策略

通過IF 組合消除電離層延遲一階項的影響后，RTPPP 待估參數(shù)包括測站坐標、接收機鐘差、天頂對流層延遲以及模糊度參數(shù)，使用卡爾曼（Kalman）濾波進行參數(shù)估計。隨機模型選取高度角模型，采用以衛(wèi)星高度角定權的隨機模型時，可以根據(jù)正弦函數(shù)模型建立隨機模型，表達式方程組為

BDS RTPPP 處理策略與GPS 類似，由文獻[15-16]可知，GEO 軌道衛(wèi)星實時產(chǎn)品精度較差，因此在進行高度角定權時對GEO 衛(wèi)星觀測值進行降權。取截止高度角為7°。RTPPP 的參數(shù)估計及誤差處理策略見表1。

表1 RTPPP 的參數(shù)估計及誤差處理策略

3 GPS/BDS 實時精密單點定位實驗及其分析

3.1 實驗說明

為了分析 GPS/BDS RTPPP 精度，選取 5 個MGEX測站 （分別為DAEJ，HKSL，HKWS，JFNG，PERT 測站）的2019 年DOY 98 觀測時段，通過BKG NTRIP Caster 接收CLK93 實時數(shù)據(jù)流對廣播星歷改正，獲取并存儲實時軌道及鐘差產(chǎn)品 （https://igs.bkg.bund.de/ntrip/orbits），實時軌道和鐘差數(shù)據(jù)采用實時SSR 改正信息。按照表1 的數(shù)據(jù)處理策略[17]，分別對實驗數(shù)據(jù)以1 d 為單位進行解算，將國際 GNSS 服務 （International GNSS Service，IGS）周解坐標作為參考值，統(tǒng)計不同測站不同時段PPP的收斂時間及定位誤差，從均方根誤差 （Root Mean Square Error，RMSE）角度對定位精度進行分析。

3.2 實驗結(jié)果分析

為了直觀呈現(xiàn)3 種方案下N，E，U 方向收斂趨勢，選取HKSL 測站定位誤差序列，見第58 頁圖1。從圖1 可以看出，3 種方案收斂趨勢較為明顯，GPS/BDS 組合及GPS 方案均能在很短時間內(nèi)收斂至穩(wěn)定值，但是BDS 系統(tǒng)RTPPP 序列存在較長時間未收斂的情況。圖1 給出了3 種方案下的可視衛(wèi)星數(shù)，由此可以看出由于GPS/BDS 雙系統(tǒng)有更多的可視衛(wèi)星數(shù)及更好的空間幾何結(jié)構(gòu)，因此表現(xiàn)性能優(yōu)于單系統(tǒng)解，且可明顯看出單BDS 表現(xiàn)性能較差，體現(xiàn)在波動較大、收斂較慢。

3.2.1 收斂時間分析

第 58 頁圖 2 為 GPS，BDS 和 GPS/BDS 3 種方案下5 個測站的RTPPP 收斂時間。對于收斂時間定義如下：從開始歷元起，某一歷元偏差小于10 cm 并且之后連續(xù)20 歷元的平均偏差也滿足此要求，則從開始歷元到該歷元的觀測時間為收斂時間。從圖2 可以看出，單BDS RTPPP 收斂時間較長，收斂到分米級均需要70 min 以上；單GPS 收斂時間優(yōu)于單BDS RTPPP，最長需要60.5 min，最短需要23 min；GPS/BDS 組合收斂時間最短，均能在20 min 左右收斂到10 cm。從平均數(shù)值角度評價，GPS，BDS 和GPS/BDS 3 種方案的平均收斂時間分別為38.8 min，94.7 min 和21.6 min，GPS/BDS雙系統(tǒng)RTPPP 收斂時間相對于單GPS 和單BDS 分別提高了44%，77%。

圖1 HKSL 測站的收斂序列及各系統(tǒng)可視衛(wèi)星數(shù)

圖2 3 種方案下各個測站的RTPPP 收斂時間

3.2.2 定位精度分析

從收斂后RMSE 角度對定位精度進行分析，RMSE 為RTPPP 解算結(jié)果與IGS 周解文件坐標之間的差值。圖3 為GPS，BDS 和GPS/BDS 方案下5 個測站的N，E，U 方向收斂后定位精度，即3 種方案下5 個測站的單天解定位偏差結(jié)果。由圖3 可以看出，3 種方案下，E，N 方向定位精度優(yōu)于U 方向，單GPS 定位精度較好，除DAEJ 測站U 方向偏差在3 cm 左右，其余收斂后偏差均在2 cm 以內(nèi)；同樣地，GPS/BDS 組合定位精度與GPS 定位精度相當，均優(yōu)于2 cm；BDS 定位精度最差，N 方向定位精度優(yōu)于5 cm，E 方向定位精度在5 cm 左右，除PERT 站U 方向定位精度為15 cm，其余測站U 方向定位精度均優(yōu)于5 cm。

圖3 3 種方案下各個測站的N，E，U 方向收斂后定位精度

表2 為3 種方案下N，E，U 方向收斂后的平均定位精度，即各方案下N，E，U 方向定位結(jié)果平均值。由表2 中各個方向的平均定位精度可以看出，GPS/BDS 在N，E，U 方向收斂后定位精度分別為0.94 cm，1.63 cm，2.25 cm，相對于單BDS 的2.73 cm，5.91 cm，7.16 cm，分別提高了 66%，72%，69%，改善效果明顯。單BDS 定位精度偏低的可能原因是目前BDS 全球跟蹤站數(shù)目較少，地面控制網(wǎng)不能很好地覆蓋全球，導致實時軌道及鐘差改正信息解算精度低于GPS。

表2 3 種方案下N，E，U 方向收斂后的平均定位精度 （m）

3.2.3 ZTD 精度分析

通過列舉DAEJ 及HKSL 測站的對流層天頂延遲 （Zenith Tropospheric Delay，ZTD）序列 （見第59 頁圖4）可以看出，GPS/BDS 所求ZTD 序列均處于2.5 m 左右，而BDS 所求ZTD 精度稍差于GPS，在初期因為可視衛(wèi)星數(shù)較少，需要較長時間收斂到穩(wěn)定值，并且BDS 所求ZTD 序列存在波動情況。在使用PPP 進行對流層反演過程中，需對BDS 進行特別關注。

圖4 DAEJ，HKSL 測站的PPP 解算ZTD 序列

4 結(jié)論

本文采用CNES 的多系統(tǒng)實時軌道和鐘差信息，實現(xiàn)了GPS，BDS 以及GPS/BDS 組合RTPPP，并對定位精度進行分析。結(jié)果表明，GPS/BDS 組合能夠較好地提高單BDS 收斂時間，單GPS 系統(tǒng)收斂效率與 GPS/BDS 組合相當，GPS/BDS 雙系統(tǒng)RTPPP 平均收斂時間為21.6 min，相對于單BDS 的94.7 min、單GPS 的38.8 min，分別提高了 77%，44%。定位精度上，GPS/BDS 組合對單GPS 改善效果不明顯，且兩種方案下N，E，U 方向的定位精度均優(yōu)于3 cm，單BDS 定位精度較差，N 方向定位精度優(yōu)于5 cm，E 方向定位精度在5 cm 左右，U 方向除個別測站外定位精度也在5 cm 左右；GPS/BDS 在N，E，U 方向收斂后的定位精度分別為0.94 cm，1.63 cm，2.25 cm，相對于單 BDS 的2.73 cm，5.91 cm，7.16 cm，分別提高了 66%，72%，69%，改善效果明顯。

由此可知，在RTPPP 方面，GPS 的加入使得BDS RTPPP 效果改善較大，一部分原因在于目前CNES 等機構(gòu)提供的公開RTPPP 服務僅局限于BDS-2 衛(wèi)星，隨著BDS-3 組網(wǎng)完成及BDS-3 系統(tǒng)完善的導航定位和授時 （Positioning, Navigation and Timing，PNT）服務體系，BDS RTPPP 的性能有望進一步提高。