丁樂樂,王珍,潘宇明,李蕾

(天津勘察院,天津 300191)

0 引 言

北斗地基增強系統(tǒng)是國家統(tǒng)一規(guī)劃建設的以北斗衛(wèi)星導航系統(tǒng)(BDS)為主,兼容其他GNSS衛(wèi)星導航系統(tǒng)的地基增強系統(tǒng)。自2012年我國BDS實現(xiàn)亞太區(qū)域定位后,已有多個省市完成了北斗增強系統(tǒng)的升級。目前,基準站升級策略是在原參考站基礎上,共用一個天線,增加一臺接收機,并且增加單獨的數(shù)據(jù)處理系統(tǒng),兩套數(shù)據(jù)處理系統(tǒng)的觀測數(shù)據(jù)及服務均未進行融合處理。

參考站間模糊度固定是北斗地基增強系統(tǒng)服務的核心技術[1]。只有快速、準確的固定參考站間整周模糊度,才能實時地為用戶提供高精度的誤差改正信息。國內外相關學者針對參考站間模糊度固定方法做過許多研究,并且取得了諸多成果[2-5]。主要方法可分為兩個方面:1)利用偽距、相位觀測數(shù)據(jù)進行差分或線性組合整周模糊度,最為經(jīng)典的算法為MW寬巷+窄巷組合確定整周模糊度;2)利用實數(shù)模糊度及其方差協(xié)方差陣進行矩陣變換,代表性的是最小二乘降相關平差法(LAMBDA)方法。

北斗地基增強系統(tǒng)形成了“一天線+雙接收機”的觀測模式，每條基線上存兩個零基線約束條件,由于零基線雙差過程信號傳播路徑相同,消除了大部分誤差,很容易實現(xiàn)模糊度固定?；诖?本文在已有模糊度固定方法的基礎上,增加零基線約束條件,進行相應的矩陣變換,以期提高模糊度固定的速度和效率。

1 GNSS雙差原始觀測方程建立

GNSS雙差觀測方程形成時,星間差在各自系統(tǒng)內求差。雙差觀測方程可表示為

λ

λ·ΔεΦ,

(1)

聯(lián)系人: 丁樂樂 E-mail: dinglelecumt@126.com

一般選取參考星原則為選取衛(wèi)星高度角最大的衛(wèi)星作為參考星。某時刻同步觀測到m顆GPS衛(wèi)星和n顆BDS衛(wèi)星,且分別選擇1顆衛(wèi)星高度角最大的GPS和BDS為各自系統(tǒng)的參考星,形成m+n-2個雙差觀測方程。

2 模糊度解算參數(shù)估計方法

參數(shù)估計通常采用擴展Kalman濾波模型。參考站間基線可認為已知,則參數(shù)估計主要包括模糊度、對流層延遲。通過使用擴展Kalman濾波,狀態(tài)向量x和它的方差協(xié)方差陣P可以與每個歷元tk觀測向量yk通過下式估計:

(2)

式中:zk、H、R分別表示量測向量、系數(shù)矩陣和觀測誤差的方差協(xié)方差陣。在擴展Kalman濾波中,假定系統(tǒng)模型為線性的,狀態(tài)方程時間更新和它的方差協(xié)方差陣表示如下:

(3)

其中:A為狀態(tài)轉移矩陣;Q為系統(tǒng)噪聲的方差協(xié)方差。

基于GNSS雙差觀測模型的卡爾曼濾波狀態(tài)方程的狀態(tài)向量x如下:

式中:P0初值為零矩陣; 狀態(tài)噪聲陣Qk的對角陣元素分別為各狀態(tài)量的方差。由于參考站間基線已知,位置參數(shù)的先驗過程噪聲為0.1 m,雙差模糊度參數(shù)的過程噪聲為20周。大于40 km的基線,狀態(tài)向量中包含對流層估計參數(shù),對流層延遲過程噪聲為0.3 m.電離層延遲誤差可通過無電離層模型消除,參數(shù)估計中可不估計雙差電離層延遲。

觀測向量Zk用雙差載波和偽距觀測值組成:

Zk= (λ1·

Hk為雙差觀測方程的線性化系數(shù):

假定觀測r顆衛(wèi)星。則系數(shù)矩陣中的各項都含有r-1行。E為(r-1)×(r-1)的單位陣。O代表零矩陣。Rk為觀測值噪聲矩陣,這里假定各顆衛(wèi)星的觀測值之間相互獨立,并且設定偽距和載波的標準差分別為σP和σφ.

Kalman濾波獲得載波相位模糊度實數(shù)值及其方差協(xié)方差矩陣后,即可采用整數(shù)LAMBDA解算整周模糊度。

3 GNSS參考站間模糊度固定方法

考慮參考站間距在20～70 km,GNSS參考站間模糊度固定策略采用“無線電離層+寬巷模糊度組合”。在中長基線解算時,無電離層組合觀測值,消除電離層影響,同時,利用組合系數(shù)保持模糊度的整周特性,形成無電離層組合模糊度,構建觀測方程可表示為

(4)

式中: LC代表無電離層組合觀測值; L1、L2為載波原始距離觀測值;N1、N2為載波原始相位觀測值;f1、f2代表GNSS觀測值信號的兩個頻率,不同衛(wèi)星系統(tǒng)頻率不同;λ1、λ2代表GNSS觀測值信號的兩個波長;ρ表示衛(wèi)星到接收機間幾何距離;由于雙差觀測值消除了接收機鐘差和衛(wèi)星鐘差,電離層延遲誤差也已經(jīng)消除。對流層誤差不能完全消除但能夠利用模型映射函數(shù)對對流層延遲殘差進行估計。

由于寬巷雙差整周模糊度具有較長波長,其模糊度相對容易固定,其固定不受基線長短的影響,對低高度角的衛(wèi)星也不敏感。

(5)

4 基于零基線約束的參考站間模糊度矩陣變換

1) 構建線性變換系數(shù)進行變換。

(6)

(7)

式中:QZZ為變換后做差模糊度參數(shù);Nzz的方差協(xié)方差陣;QAC·ZZ為NAC與Nzz之間的協(xié)方差陣。

式(2)變換后的模糊度參數(shù)表達為

NAC-NAD=-NCD,

NAC-NBC=NAB,

NAC-NBD=NAB-NCD.

(8)

由于A/B、C/D接收機為同一天線兩臺接收機的模式,在A與B、C與D之間進行雙差時衛(wèi)星至接收機的電離層延遲、對流層延遲均相同,因此,能夠實現(xiàn)A/B、C/D模糊度參數(shù)的單歷元固定,相應的NAB-NCD能夠實現(xiàn)固定。模糊度A/C及其方差協(xié)方差更新方程為

(9)

利用模糊度之間的相關性,修正了NAC模糊度以及NAC模糊度的方差協(xié)方差;利用LAMBDA方法進行模糊度固定。

5 數(shù)據(jù)處理及結果分析

試驗選取天津市北斗地基增強系統(tǒng)中的兩個參考站KC04和KC09站作為參考試驗站,實例數(shù)據(jù)采用2016年5月1日的天寶接收機和南方接收機1 s數(shù)據(jù)進行試驗,接收機觀測數(shù)據(jù)分別命名為NF04、TB04、NF09、TB09,參考站間距離為60.2 km.讀取以上兩個參考站四個接收機的觀測值數(shù)據(jù),對四個接收機的觀測數(shù)據(jù)進行周跳探測等數(shù)據(jù)預處理。數(shù)據(jù)處理流程如圖1所示。

圖1 零基線約束的參考站間模糊度計算流程圖

按照論文第2、3節(jié),形成無電離層組合觀測方程,利用Kalman濾波模型進行參數(shù)估計,估計參數(shù)為無電離層實數(shù)模糊度以及對流層延遲參數(shù)。通過第3節(jié)給定的系數(shù)將無電離層組合模糊度分解為寬巷模糊度和窄巷模糊度;寬巷模糊度對高度角的影響不敏感且易固定,一般通過多歷元取均值即可固定。

圖2示出了基于觀測值直接組合求解的雙差浮點寬巷模糊度。從圖中可知,只要通過數(shù)個歷元的平滑即可得到寬巷周周模糊度。通過取整可得到此寬巷模糊度值為2 619 635周。寬巷模糊度一般在幾個歷元就可以正確固定。

設置觀測值的偽距與相位權重比為1∶100,卡爾曼濾波中可以設置坐標狀態(tài)參數(shù)的精度為0.05,模糊度參數(shù)的精度設置為20.由于增加了多余觀測值,同時包含有零基線約束,收斂過程要快得多。從圖3中可知,模糊度正確固定后的整周數(shù)為942 774周,零基線多余觀測下的參考站間窄巷實數(shù)模糊度收斂速度要快得多,收斂結果可以看出,零基線約束條件下的窄巷實數(shù)模糊度更接近固定后的整周模糊度。

圖2 MW組合雙差模糊度實數(shù)解

圖3 卡爾曼濾波某衛(wèi)星的窄巷實數(shù)模糊度

對流層延遲不能通過觀測值組合消除,論文處理方式為模型改正+參數(shù)估計的方法消除對流延遲的影響。對流層延遲利用天頂延遲和映射函數(shù)模型共同表達,對流層延遲參數(shù)按照逐歷元進行估計,分為干延遲和濕延遲。逐歷元的參數(shù)估計結果如圖4所示。

在模糊度域內進行線性變換;利用A/C-A/D、A/C-B/C、A/C-B/D線性組合的模糊度單元歷元可實現(xiàn)固定的特性,對變換后的模糊度利用LAMBDA方法固定;利用模糊度之間的相關性,用固定的零基線模糊度修正A/C窄巷模糊度實數(shù)解及其方差協(xié)方差陣;利用LAMBDA方法對A/C窄巷模糊度進行固定。

圖4 對流層天頂延遲估計結果

圖5 零基線約束模糊度固定與原始模糊度固定Ratio值對比

模糊度固定時的Ratio值對比如圖5所示。其中,N1表示分離的模糊度固定時的Raito值,ZN零基線模糊度固定時的Raito值,N1′表示加零基線約束后模糊度固定時的Raito值。從圖5中可以看出,零基線參考站間模糊度固定結果較原始的模糊度固定速度和效率有極大的提高。通過統(tǒng)計可以看出,零基線模糊度固定率為100%,原始N1模糊度固定率為95.06%,零基線約束后的模糊度固定率為99.91%.從前500歷元看,原始模糊度固定在300歷元左右實現(xiàn)固定,固定時間約為5 min.零基線約束后模糊度在第1個歷元就已經(jīng)固定。從初始化速度考慮,零基線參考站間模糊度可以在極短時間里固定。

為進一步零基線約束條件下模糊度固定的比率,選擇近6個小時,采樣率為1 s的觀測數(shù)據(jù)進行初始化。每60個歷元分別采用兩種方法固定載波相位模糊度,Ratio閾值設置為2.初始化固定情況如表1所示.由表1可知,零基線約束模糊度固定方法極大地提高了模糊度解算的成功率。

表1 零基線約束模糊度固定統(tǒng)計表

6 結束語

論文論述了GNSS零基線參考站間模糊度快速固定方法,該方法基于“一天線+雙接收機”設計進行參考站間模糊度固定,并利用了零基線參考站的觀測數(shù)據(jù)進行模糊度解算,增加了多余觀測,將一個天線下的兩臺接收機觀測數(shù)據(jù)融合處理。同時,利用零基線作為模糊度固定的約束條件,用于解算參考站間模糊度。與傳統(tǒng)參考站間模糊度固定方法比,加快了模糊度收斂速度,極大地提高了模糊度固定的速度和成功率。

[1] 李征航,張小紅.衛(wèi)星導航定位新技術及高精度數(shù)據(jù)處理方法[M].武漢:武漢大學出版社,2009.

[2] 陽仁貴,袁運斌,歐吉坤.相位實時差分技術應用于飛行器交會對接研究[J]. 中國科學:物理學力學天文學,2015,40(5):651-657.

[3] 祝會忠,徐愛功,高猛等.BDS網(wǎng)絡RTK中距離參考站整周模糊度單歷元解算方法[J].測繪學報,2016,45(1):50-57.

[4] 鄧建,王慶,潘樹國.網(wǎng)絡RTK參考站間低仰角衛(wèi)星模糊度快速解算方法[J].儀器儀表學報,2010,31(6):1201-1206.

[5] 韓保民,歐吉坤.一種附約束的單頻單歷元GPS雙差相位解算方法[J]. 測繪學報,2002, 31 (4):300-302.

[6] 高星偉,李毓麟,葛茂榮.GPS/GLONASS相位差分的數(shù)據(jù)處理方法[J].測繪科學,2004,29(2)22-24.

[7] 周樂韜,黃丁發(fā).一種參考站間雙差模糊度快速解算策略[J].大地測量與地球動力學,2006,26 (4):34-40.

[8] 羅孝文,歐吉坤.中長基線GPS網(wǎng)絡RTK模糊度快速解算的一種新方法[J].武漢大學學報(信息科學版), 2007, 32(2):156-159.

[9] 唐衛(wèi)明,劉經(jīng)南，施闖，等.三步法確定網(wǎng)絡RTK基準站雙差模糊度[J].武漢大學學報(信息科學版), 2007,32 (4):305-308.

[10]李博峰,沈云中,周澤波.中長基線三頻GNSS模糊度的快速算法[J]. 測繪學報,2009, 38 (4):296-301.

[11]YANG, Y X, HE H B, XU G C. A new adaptively robust filtering for kinematic geodetic positioning [J]. Journal of Geodesy, 2001，73(2-3): 109-116.

[12]TEUNISSEN P J G, VERHAGEN S.The GNSS ambiguity ratio-test revisited: A better way of using it[J].Survey Review,2009，41(312):138-151.