徐梅,嚴超,喬方,蘇迪

(安徽理工大學 測繪學院,安徽 淮南 232001)

0 引 言

北斗衛(wèi)星導航系統(tǒng)(BDS)是中國著眼于國家安全和經(jīng)濟社會發(fā)展需要,自主建設、獨立運行的衛(wèi)星導航系統(tǒng),是為全球用戶提供全天候、全天時、高精度的定位、導航和授時服務的國家重要空間基礎設施[1]。BDS的空間星座由對地靜止軌道(GEO)衛(wèi)星、傾斜地球同步軌道(IGSO)衛(wèi)星和中高軌道(MEO)衛(wèi)星三類衛(wèi)星組成,且全部發(fā)射三頻信號[2],與全球定位系統(tǒng)(GPS)、GLONASS和GALILEO相比,BDS增加了軌道高度較高、運動角速度較慢的GEO和IGSO衛(wèi)星。截止2017年8月,北斗二代衛(wèi)星在軌情況如表1所示。

由于GEO和IGSO衛(wèi)星的高軌道特性,其可見性要遠遠優(yōu)于MEO衛(wèi)星,能大幅度增加觀測時段內(nèi)的衛(wèi)星數(shù),對增強區(qū)域?qū)Ш叫阅芊浅Ｓ欣鸞3]。文獻[4]和[5]提出了附有約束的BDS單頻單歷元基線解算;文獻[6]基于雙頻偽距和寬巷載波相位雙差觀測方程,實現(xiàn)BDS雙頻單歷元基線解算;文獻[7]和[11]按照固定難度從小到大依次固定超寬巷、寬巷和窄巷或基頻模糊度,實現(xiàn)了BDS三頻單歷元基線解算。上述文獻對于BDS單頻、雙頻和三頻基線解算做了較為深入的研究,但這些研究僅限于從整體上分析北斗混合星座對于模糊度固定和定位精度的影響,并未深入研究IGSO衛(wèi)星的區(qū)域定位性能。對此,本文在已有的研究基礎上,采用內(nèi)蒙古地區(qū)和安徽地區(qū)的兩條短基線,比較分析IGSO衛(wèi)星對BDS三頻單歷元基線解算的影響。

表1 北斗二代衛(wèi)星在軌情況統(tǒng)計

聯(lián)系人: 徐梅E-mail:xm1019aust@163.com

1 函數(shù)模型和隨機模型

1.1 BDS三頻單歷元基線解算函數(shù)模型

三頻情況下,頻率的多樣性可以提高無幾何模糊度解算方法的可靠性[10]。BDS三頻觀測時,可以對載波觀測值進行組合得到超寬巷(EWL)組合,根據(jù)文獻[10]和文獻[11],選擇(0,-1,1)和(1,4,-5)兩個超寬巷組合,EWL載波觀測值與3個偽距觀測值PB1,PB2以及PB3組成的基于幾何模型的雙差觀測方程可表示為

(1)

采用上述方法,可以成功固定(0,-1,1)和(1,4,-5)這兩個WEL的ΔN(0,-1,1)和ΔN(1,4,-5)參數(shù),將其看成高精度的距離觀測值。由于短基線雙差大氣延遲誤差較小,可以采用無電離層觀測值消除電離層殘差,因此可采用下式來求解BDS的基礎模糊度:

(2)

利用固定的BDS雙差基礎模糊度,求解參數(shù)X,即:

VB1=AX+(LB1+λ1ΔN1)

(3)

采用經(jīng)典最小二乘估計得到參數(shù)X:

X=-(ATPA)-1(ATPL)

(4)

式中:P為權(quán)重;L=LB1+λ1ΔN1.

1.2 隨機模型

BDS不同軌道衛(wèi)星具有不同的性質(zhì),且不同頻點信號質(zhì)量之間具有差異,這是BDS單歷元基線解算中不可忽略的問題?；谛l(wèi)星高度角的隨機模型[12]和基于信噪比的隨機模型[13]在一定程度上能較為準確地反映觀測值的觀測質(zhì)量,而因BDS不同軌道衛(wèi)星高度不同的特點,也需考慮衛(wèi)地距的影響。綜上所述,所以可采用高度角、衛(wèi)地距與信噪比組合模型進行定權(quán)[14]:

(5)

式中:S為縮放因子,S=

SNRB1、SNRB2、SNRB3分別為B1、B2、B3頻點的信噪比值；ρ表示衛(wèi)地距。

2 實驗分析

本文選擇了位于中國安徽省淮南和內(nèi)蒙古自治區(qū)鄂爾多斯2個地區(qū)所觀測的BDS三頻短基線數(shù)據(jù)(<2 km)。其中,安徽淮南地區(qū)數(shù)據(jù)采樣開始日期為GPS時2017年7月4日4時16分,取3600個歷元數(shù)據(jù),采樣間隔1 s,基準站為JZJK,流動站為CS01;內(nèi)蒙古自治區(qū)鄂爾多斯地區(qū)數(shù)據(jù)采樣開始日期為GPS時2016年10月30日零時47分,取3600個歷元數(shù)據(jù),采樣間隔1 s,基準站為JZTJ,流動站為CTJH。衛(wèi)星截止高度角均設置為10°。

圖1示出了觀測時段內(nèi)IGSO衛(wèi)星信噪比隨高度角的變化關(guān)系。表明觀測時段內(nèi)可見5顆IGSO衛(wèi)星(兩條基線均沒觀測到C13號衛(wèi)星),但是可見性有所不同,這是因為5顆IGSO衛(wèi)星均勻分布在3個傾斜同步軌道面上,其運動軌跡以赤道為對稱軸的“8”字形,南北跨度較大,且每顆IGSO衛(wèi)星覆蓋范圍限定于某一固定區(qū)域;IGSO衛(wèi)星B1、B2、B3頻點信噪比與高度角之間有一定的趨勢性關(guān)系,即:衛(wèi)星高度角越大,信噪比越大(CJTH站C08因為高度角變化較小,大約1°,所以規(guī)律性不明顯)。

圖1 流動站IGSO衛(wèi)星信噪比隨高度角的變化

根據(jù)IGSO衛(wèi)星可見性,逐次增加一顆IGSO衛(wèi)星,具體方案如下:

方案一:排除觀測時段內(nèi)所有的IGSO衛(wèi)星,僅利用GEO和MEO衛(wèi)星進行定位,參考衛(wèi)星選用C03;

方案二:在方案一的基礎上,兩測區(qū)均增加C06;

方案三:在方案二的基礎上,安徽測區(qū)增加C07,內(nèi)蒙古測區(qū)增加C08;

方案四:在方案三的基礎上,兩測區(qū)均增加C09;

方案五:在方案四的基礎上,安徽測區(qū)增加C10,內(nèi)蒙古測區(qū)增加C07;

方案六:在方案五的基礎上,安徽測區(qū)增加C08,內(nèi)蒙古測區(qū)增加C10.

2.1 模糊度解算結(jié)果分析

對上述兩條基線采用BDS三頻單歷元基線解算的數(shù)學模型進行單歷元模糊度解算,解算結(jié)果的有效性檢驗可以用成功率進行衡量[15-16]。安徽測區(qū)基線以淮南礦業(yè)(集團)有限責任公司朱集東礦GNSS網(wǎng)數(shù)據(jù)處理結(jié)果作為參考值,內(nèi)蒙古測區(qū)基線以淮南礦業(yè)(集團)有限責任公司唐家會煤礦GNSS網(wǎng)數(shù)據(jù)處理結(jié)果作為參考值。分別將各歷元固定解與參考值進行比較,相同則認為是正確的。

成功率=正確固定歷元數(shù)/有效歷元數(shù)

兩條基線采用不同方案模糊度解算的成功率如表2所示。

表2 不同方案下模糊度解算成功率統(tǒng)計

從表2可知:

1) 方案一中只存在GEO衛(wèi)星的情況下(兩測區(qū)在測段內(nèi)均未接收到MEO衛(wèi)星信號),仍能達定位要求,查看共視數(shù)據(jù)發(fā)現(xiàn),在兩測區(qū)的觀測時段內(nèi)均可接收到5顆GEO衛(wèi)星,充分體現(xiàn)了GEO衛(wèi)星作為地球靜止軌道衛(wèi)星對于中國區(qū)域的增強效果,但也說明只有GEO衛(wèi)星參與解算時,BDS三頻模糊度解算成功率較低,大約為30%;

2) 與方案一相比,方案二雖然只是增加了1顆IGSO衛(wèi)星,但是成功率至少提升了246.5%,成功率高達90%以上,體現(xiàn)了IGSO衛(wèi)星對于BDS三頻模糊度解算成功率具有重要意義;

3) 由方案二、三、四、五、六可知,安徽測區(qū)增加3顆IGSO衛(wèi)星模糊度解算成功率最高,內(nèi)蒙古測區(qū)增加2顆IGSO衛(wèi)星模糊度解算成功率最高,結(jié)合圖1,這主要是兩測區(qū)可見性最好的2～3顆IGSO衛(wèi)星均參與計算,沒有衛(wèi)星升降對其解算結(jié)果造成影響;當增加4～5顆IGSO衛(wèi)星時,模糊度解算成功率相對于2～3顆有所下降,這可能是引入了觀測質(zhì)量較差的IGSO衛(wèi)星。

2.2 定位精度分析

本文在進行精度分析時,僅從定位結(jié)果角度出發(fā),判斷模糊度正確固定的歷元的解算結(jié)果,與參考值作差,得到兩條基線北(N)、東(E)、高程(U)三方向的定位結(jié)果,如圖2、圖3所示,并利用內(nèi)符合精度(STD)進行精度分析,如表3所示。

圖2 安徽測區(qū)基線在N、E、U方向上的定位偏差 (a)方案一； (b)方案二； (c)方案三； (d)方案四；(e)方案五； (f)方案六

圖3 內(nèi)蒙古測區(qū)基線在N、E、U方向上的定位偏差 (a)方案一； (b)方案二； (c)方案三； (d)方案四；(e)方案五； (f)方案六

表3 內(nèi)符合精度STD統(tǒng)計 單位：mm

基線模式方案一二三四五六安徽測區(qū)N54.05.05.04.44.64.9E3.23.03.63.34.06.1U39.611.712.713.326.831.2內(nèi)蒙古測區(qū)N51.024.15.05.15.15.1E2.72.42.22.12.32.4U45.524.09.513.013.813.5

結(jié)合圖2、圖3和表3可知:

1) 方案一中,N方向定位精度遠遠低于E方向,這是因為參與方案一解算的衛(wèi)星為5顆GEO衛(wèi)星,5顆GEO衛(wèi)星分別定點于東經(jīng)58.75°、80°、110.5°、140°和160°,而我國位于北半球;

2) 與方案一相比,方案二、三、四逐步增加一顆IGSO衛(wèi)星,N、U方向定位誤差相對于方案一均明顯減少;在方案五、六中,增加兩顆可見性較差的衛(wèi)星時,N、E、U方向定位誤差相對于方案三、四均有所降低,這是因為方案五、六中增加的兩顆衛(wèi)星高度角較低,且信噪比明顯低于其他衛(wèi)星;

3) 由方案二、三、四可知,分別增加1～3顆可見性較好的IGSO衛(wèi)星時,N、E、U方向定位精度達到最佳;方案二、三、四、五、六中,U方向定位偏差基本偏向于同一側(cè),這可能是GEO衛(wèi)星多路徑效應引起的;解算結(jié)果N、E方向優(yōu)于U方向,E方向優(yōu)于N方向,誤差分布與目前BDS服務性能一致。

3 結(jié)束語

與GPS相比,BDS的高軌道GEO衛(wèi)星和IGSO衛(wèi)星是其一大優(yōu)勢,除了其可見時段長、有利于增加用戶的可見衛(wèi)星數(shù)之外,高軌道衛(wèi)星信號更容易被接收機接收,同時能有效減少低仰角衛(wèi)星引入的多路徑誤差影響,尤其是在高樓建筑物阻擋的惡劣環(huán)境下,北斗系統(tǒng)特殊的星座構(gòu)型能有效保證定位結(jié)果的穩(wěn)定性。特別是IGSO衛(wèi)星不僅具有GEO衛(wèi)星的區(qū)域服務能力,還可以同時保證對高緯度地區(qū)的覆蓋。本文利用我國內(nèi)蒙古地區(qū)和安徽地區(qū)的實測數(shù)據(jù),分析了IGSO對BDS三頻單歷元基線解算的影響。由數(shù)據(jù)分析可得:1) IGSO衛(wèi)星B1、B2、B3頻點信噪比隨著高度角的增大而增大;2) 當只有GEO衛(wèi)星參與定位時,BDS三頻模糊度解算成功率低于30%,N、U方向定位誤差較大,能達到5 cm左右;3) 按可見性逐次增加IGSO衛(wèi)星參與定位時,增加1～3顆可見性較好的IGSO衛(wèi)星時,模糊度解算成功率最高,且定位精度最佳;4)解算結(jié)果N、E方向優(yōu)于U方向,E方向優(yōu)于N方向,誤差分布與目前BDS服務性能一致。

[1] 北斗衛(wèi)星導航系統(tǒng).北斗衛(wèi)星導航系統(tǒng)介紹[EB/OL].(2010-01)[2017-05-23]. http://www.beidou.gov.cn/xtjs.html.

[2] 楊元喜. 北斗衛(wèi)星導航系統(tǒng)的進展、貢獻與挑戰(zhàn)[J]. 測繪學報, 2010, 39(1):1-6.

[3] 吳甜甜,張云,劉永明, 等. 北斗/GPS組合定位方法[J]. 遙感學報,2014, 18(5):1087-1097.

[4] 王靜,趙興旺,劉超.基于正則化的GPS/BDS單頻單歷元模糊度固定[J].大地測量與地球動力學,2016,36(12):1083-1087.

[5] 張書畢,劉鑫,宋冰,等. 附有約束的BDS單頻單歷元改進型Par Lambda算法[J]. 中國礦業(yè)大學學報, 2017(1):201-208.

[6] 唐衛(wèi)明,鄧辰龍,高麗峰. 北斗單歷元基線解算算法研究及初步結(jié)果[J]. 武漢大學學報(信息科學版), 2013, 38(8):897-901.

[7] TANG W, DENG C, SHI C,etal. Triple-frequency carrier ambiguity resolution for Beidou navigation satellite system[J]. GPS Solutions, 2014, 18(3):335-344.

[8] 劉炎炎,葉世榕,江鵬,等. 基于北斗三頻的短基線單歷元模糊度固定[J]. 武漢大學學報(信息科學版), 2015, 40(2):209-213.

[9] 何俊,劉萬科,張小紅. 北斗短基線三頻實測數(shù)據(jù)單歷元模糊度固定[J]. 武漢大學學報(信息科學版), 2015, 40(3):361-365.

[10]謝建濤,郝金明,韓聰,等. BDS單歷元TCAR算法優(yōu)化研究[J]. 測繪科學技術(shù)學報, 2016, 33(1):6-10.

[11]呂偉才,高井祥,王堅,等. 北斗三頻約束的短基線模糊度單歷元算法[J]. 中國礦業(yè)大學學報, 2015, 44(6):1090-1096.

[12]RIZOS C. Principles and practice of GPS surveying[M].School of Geomatic Engineering, The University of New South Wales, Sney, 1997.

[13]LANGLEY R. GPS receiver system noise[J].GPS World, 1997(8):40-45.

[14]劉乾坤,隋立芬,肖國銳,等. 北斗系統(tǒng)差分碼偏差解算中一種新的定權(quán)方法[J]. 測繪科學技術(shù)學報, 2015, 32(5):473-478.

[15]劉經(jīng)南,鄧辰龍,唐衛(wèi)明. GNSS整周模糊度確認理論方法研究進展[J]. 武漢大學學報(信息科學版), 2014, 39(9):1009-1015.

[16]DENG C L, TANG W M. Reliable single-epoch ambiguity resolution for short baselines using GPS/BeiDou system[J]. GPS Solutions, 2014, 3(18):375-386.