徐龍威　劉　暉　劉玉潔　舒　寶　張　明　錢　闖

1　武漢大學(xué)衛(wèi)星導(dǎo)航定位技術(shù)研究中心，武漢市珞喻路129號(hào)，430079 2　河南省遙感測(cè)繪院，鄭州市黃河路8號(hào)，450003

?

一種顧及GNSS系統(tǒng)間偏差的偽距單點(diǎn)定位方法

徐龍威1劉暉1劉玉潔2舒寶1張明1錢闖1

1武漢大學(xué)衛(wèi)星導(dǎo)航定位技術(shù)研究中心，武漢市珞喻路129號(hào)，430079 2河南省遙感測(cè)繪院，鄭州市黃河路8號(hào)，450003

對(duì)系統(tǒng)間偏差的成分和影響因素進(jìn)行分析，發(fā)現(xiàn)各項(xiàng)改正后殘余誤差對(duì)系統(tǒng)間偏差估值的影響可歸結(jié)為與衛(wèi)星編號(hào)和廣播星歷更新相關(guān)。提出一種新的顧及參與解算衛(wèi)星構(gòu)成的系統(tǒng)間偏差估計(jì)方案，新方案能夠在僅有4顆多模GNSS衛(wèi)星可見時(shí)提供可靠的位置服務(wù)。

多模GNSS；系統(tǒng)間偏差；衛(wèi)星編號(hào)；定位

對(duì)于常規(guī)導(dǎo)航用戶，觀測(cè)條件較差導(dǎo)致大量觀測(cè)誤差被系統(tǒng)間偏差(inter system bias, ISB)吸收[1-2]，ISB的穩(wěn)定性變?nèi)酢ontenbruck等[3]在分析CONGO網(wǎng)絡(luò)測(cè)站性能時(shí)，首次提出GPS/Galileo測(cè)距碼ISB，并從理論上證明ISB與接收機(jī)和測(cè)距碼的類型有關(guān)。Torre等[2]利用實(shí)測(cè)數(shù)據(jù)分別基于精密軌道和廣播星歷分析ISB的變化規(guī)律，證明了ISB與接收機(jī)類型相關(guān)的特性。文獻(xiàn)[4-5]利用ISB短時(shí)間內(nèi)變化不大的特性，將在可見衛(wèi)星充足歷元求得的ISB作為已知參數(shù)，引入可見衛(wèi)星不足歷元進(jìn)行系統(tǒng)間的偏差補(bǔ)償，減少觀測(cè)模型中未知參數(shù)的個(gè)數(shù)，在損失一定精度的條件下實(shí)現(xiàn)了GPS/GLONASS雙系統(tǒng)僅4顆星的定位解算。由于不同衛(wèi)星觀測(cè)誤差的差異明顯且會(huì)被ISB吸收，所以不同衛(wèi)星參與解算會(huì)得到不同的ISB估值。因此，在進(jìn)行系統(tǒng)間偏差補(bǔ)償時(shí)，顧及當(dāng)前歷元可見衛(wèi)星的觀測(cè)誤差，能夠減少Cai等[4]提到的精度損失。

1　系統(tǒng)間偏差構(gòu)成和影響因素

多模GNSS時(shí)空基準(zhǔn)和信號(hào)硬件延遲等因素的差異都會(huì)導(dǎo)致系統(tǒng)間偏差，可利用偽距觀測(cè)方程對(duì)系統(tǒng)間偏差的構(gòu)成和影響因素進(jìn)行推導(dǎo):

(1)

時(shí)空轉(zhuǎn)換后，常規(guī)偽距單點(diǎn)定位誤差模型為：

V=HX+L

(2)

式中，δx、δy、δz為接收機(jī)真實(shí)坐標(biāo)與近似坐標(biāo)之差，dtsys_GR、dtsys_GC、dtsys_GE為以GPS為基準(zhǔn)與其他GNSS系統(tǒng)的接收機(jī)鐘差參數(shù)之差，即系統(tǒng)間碼偏差。此時(shí)多種觀測(cè)誤差會(huì)導(dǎo)致ISB的估值包含時(shí)間系統(tǒng)間偏差、軌道誤差、衛(wèi)星鐘差、對(duì)流層延遲、多路徑、接收機(jī)硬件延遲以及觀測(cè)噪聲之差等。Montenbruck等[6]對(duì)GNSS廣播星歷偏差進(jìn)行分析，廣播星歷SSIRE在m級(jí)，同一時(shí)刻不同衛(wèi)星廣播星歷偏差不同，且GLONASS衛(wèi)星廣播星歷誤差受星歷更新影響具有明顯周期性。大氣延遲和多路徑效應(yīng)主要與衛(wèi)星信號(hào)傳播路徑相關(guān)，可以理解為與衛(wèi)星編號(hào)相關(guān)。另外，GLONASS衛(wèi)星頻率不同，導(dǎo)致不同衛(wèi)星的信號(hào)接收機(jī)硬件碼延遲不同[7]，也可被當(dāng)作與衛(wèi)星編號(hào)相關(guān)。衛(wèi)星觀測(cè)值的各項(xiàng)觀測(cè)誤差都會(huì)影響系統(tǒng)間偏差估計(jì)，這些誤差都可以視為與衛(wèi)星編號(hào)相關(guān)且短時(shí)間內(nèi)變化不大。因此，在估計(jì)ISB時(shí)顧及衛(wèi)星編號(hào)能夠明顯提高估值對(duì)當(dāng)前歷元的適應(yīng)性。

ISB的分析基于偽距單點(diǎn)定位數(shù)學(xué)模型，對(duì)流層延遲采用Saastamoinen模型改正，電離層采用消電離層組合(GPSL1/L2，GLONASSL1/L2，BDSB1/B2，GalileoE1/E5a)，實(shí)驗(yàn)分析均基于GPST。下文中，G代表GPS，R代表GLONASS，C代表BDS，E代表Galileo。

2　系統(tǒng)間偏差特性分析

2.1GPS作為參考系統(tǒng)合理性分析

本文以GPS為基準(zhǔn)系統(tǒng)。不同的GPS衛(wèi)星參與解算會(huì)對(duì)ISB估值造成影響，因此，本文利用不同數(shù)目的GPS衛(wèi)星參與解算，分析其對(duì)ISB估計(jì)的影響。如圖1所示(以G-RISB為例)，不同數(shù)目的GPS衛(wèi)星參與解算，ISB估值差異在5ns以內(nèi)。隨著參與解算的GPS衛(wèi)星數(shù)目的減少，觀測(cè)模型的穩(wěn)健性減弱，ISB穩(wěn)定性也隨之變差。參與解算的GPS衛(wèi)星為隨機(jī)選取，其他系統(tǒng)可見衛(wèi)星全部參與解算?？傊煌w數(shù)的GPS衛(wèi)星參與解算對(duì)ISB估值的影響遠(yuǎn)小于ISB估值的量級(jí)。

圖1　GPS衛(wèi)星對(duì)ISB的影響Fig.1　The influence of GPS satellite to inter system bias

2.2系統(tǒng)間偏差特性分析

選取6個(gè)不同類型接收機(jī)的MEGX跟蹤站(表1)2013-11-25的觀測(cè)數(shù)據(jù)，截止高度角設(shè)置為10°，采樣間隔30s。

從圖2～4可以看出，ISB與接收機(jī)類型相關(guān)。對(duì)于GLONASS(圖2)，ISB比較穩(wěn)定，24h變化區(qū)間在15ns以內(nèi)，與時(shí)間相關(guān)性很強(qiáng)，相鄰兩個(gè)歷元差異很小。但I(xiàn)SB經(jīng)常出現(xiàn)跳躍，且發(fā)生時(shí)刻均為GLONASS廣播星歷更新時(shí)刻，可知廣播星歷更新是ISB的重要影響因素。

表1　MEGX站接收機(jī)類型以及兼容系統(tǒng)

圖2　不同類型接收機(jī)獲得觀測(cè)值的G-R ISB估值Fig.2　G-R ISB estimates for different types of receiver observations

從圖3看出，所有測(cè)站的ISB估值相鄰歷元間差異很小，CUT0和MAL2的ISB穩(wěn)定性較強(qiáng)。這兩個(gè)站ISB穩(wěn)定性好的主要原因是其可見BDS衛(wèi)星大于6顆，ISB估值為多顆衛(wèi)星觀測(cè)值的加權(quán)平均，衛(wèi)星數(shù)越多穩(wěn)定性越好。BRUX僅能觀測(cè)到3顆左右BDS衛(wèi)星，導(dǎo)致其ISB穩(wěn)定性較差。

圖3　不同類型接收機(jī)獲得觀測(cè)值的 G-C ISB估值Fig.3　G-C ISB estimates for different types of receiver observations

如圖4，目前Galileo在軌衛(wèi)星較少，導(dǎo)致G-E ISB穩(wěn)定性略差，不同接收機(jī)間的ISB差異也十分明顯。

圖4　不同類型接收機(jī)獲得觀測(cè)值的G-E ISB估值Fig.4　G-E ISB estimates for different types of receiver

2.3單顆衛(wèi)星對(duì)ISB影響分析

對(duì)單顆衛(wèi)星ISB的分析見圖5～7。單顆衛(wèi)星的ISB估計(jì)方法與系統(tǒng)間ISB估計(jì)方法類似。以1顆GLONASS衛(wèi)星的ISB為例，所有GPS衛(wèi)星+1顆GLONASS衛(wèi)星解算獲得該顆GLONASS衛(wèi)星的ISB。

由圖5可知，GLONASS衛(wèi)星ISB變化與時(shí)間相關(guān)，具有明顯的周期性，且與廣播星歷更新時(shí)間相符。與下文其他GNSS系統(tǒng)相比，不同衛(wèi)星對(duì)應(yīng)的ISB差異更明顯，R14與R08對(duì)應(yīng)的ISB相差20 ns左右。這主要是由于GLONASS采用頻分多址技術(shù)，不同的衛(wèi)星信號(hào)頻率不同，導(dǎo)致接收機(jī)硬件延遲差異。

圖5　不同GLONASS衛(wèi)星G-R ISB(GOP6站)Fig.5　G-R ISB for different GLONASS satellites(GOP6 station)

圖6為不同BDS衛(wèi)星ISB變化規(guī)律，較差穩(wěn)定且無(wú)明顯周期性，但與時(shí)間的相關(guān)性明顯，星歷歷元估值差異在3 ns以內(nèi)。當(dāng)衛(wèi)星高度角較低時(shí)，ISB與衛(wèi)星高度角的相關(guān)性非常明顯，這主要是由于BDS碼偏差的作用。

圖6　不同BDS衛(wèi)星G-C ISB (CUT0站)Fig.6　G-C ISB for different BDS satellites(CUT0 station)

如圖7所示，不同的Galileo衛(wèi)星ISB相差不大，單顆衛(wèi)星ISB在固定區(qū)間內(nèi)隨機(jī)分布，穩(wěn)定性優(yōu)于GLONASS和BDS，變化范圍一般在10 ns以內(nèi)。觀測(cè)時(shí)段GPST 00:00～02:00之間，ISB出現(xiàn)劇烈跳動(dòng)是由E20衛(wèi)星觀測(cè)噪聲造成。單顆Galileo衛(wèi)星ISB未出現(xiàn)類似GLONASS的周期性波動(dòng)，也未出現(xiàn)類似BDS的與高度角相關(guān)變化。但目前Galileo衛(wèi)星過(guò)少，廣播星歷播發(fā)穩(wěn)定性不足，許多時(shí)段的廣播星歷衛(wèi)星信號(hào)群延差改正參數(shù)(BGD)均不可用。

圖7　不同Galileo衛(wèi)星G-E ISB(MAL2站)Fig.7　G-E ISB for different Galileo satellites (MAL2 station)

3　ISB在定位中的應(yīng)用

3.1可見衛(wèi)星不足條件下多模組合單點(diǎn)定位

多模組合單點(diǎn)定位通常采用增加鐘差參數(shù)的方法處理，最少需要3+n顆(n為參與解算的GNSS系統(tǒng)種類)可見衛(wèi)星才能實(shí)現(xiàn)定位解算。實(shí)時(shí)低精度導(dǎo)航用戶的定位環(huán)境比較復(fù)雜，經(jīng)常出現(xiàn)衛(wèi)星信號(hào)被遮擋，導(dǎo)致只有少數(shù)衛(wèi)星可見。此時(shí)，通常利用ISB短時(shí)間內(nèi)的穩(wěn)定性，將基于觀測(cè)值充足歷元求取的ISB作為先驗(yàn)值，引入觀測(cè)模型以減少未知參數(shù)個(gè)數(shù)，實(shí)現(xiàn)衛(wèi)星數(shù)不足條件下的單點(diǎn)定位。然而不同的衛(wèi)星參與解算會(huì)得到不同的ISB估值，此時(shí)的ISB先驗(yàn)值不適用于當(dāng)前歷元?；诒疚奶岢龅腎SB估計(jì)方法，利用一個(gè)存儲(chǔ)下來(lái)的歷史歷元觀測(cè)值，分別求取當(dāng)前觀測(cè)值不足歷元可見衛(wèi)星與GPS系統(tǒng)間的ISB，并用于當(dāng)前歷元定位。

3.2實(shí)驗(yàn)分析

利用6個(gè)MGEX跟蹤站(表1)2013-11-25的觀測(cè)數(shù)據(jù)，設(shè)計(jì)3種實(shí)驗(yàn)方案，并根據(jù)方案要求隨機(jī)選取相應(yīng)顆數(shù)的衛(wèi)星觀測(cè)值。

方案1：隨機(jī)選取4顆GPS衛(wèi)星進(jìn)行單點(diǎn)定位；方案2：隨機(jī)選取來(lái)自多系統(tǒng)的4顆衛(wèi)星，引入不考慮衛(wèi)星號(hào)的ISB進(jìn)行定位；方案3：隨機(jī)選取來(lái)自多系統(tǒng)的4顆衛(wèi)星，引入顧及衛(wèi)星號(hào)的ISB進(jìn)行定位。

方案1中，隨機(jī)選取參與解算的GPS衛(wèi)星。方案2和方案3中，CUT0、MAL2和BRUX站為在當(dāng)前歷元可見衛(wèi)星中隨機(jī)選取1顆GPS、1顆GLONASS、1顆BDS和1顆Galileo衛(wèi)星。GOP6、CONZ和WTZZ站不能接收到BDS測(cè)距信號(hào)，于是隨機(jī)選取2顆GPS、1顆GLONASS和1顆Galileo衛(wèi)星。設(shè)置存儲(chǔ)歷元每5 min更新一次，截止高度角10°，歷元間隔30 s， PDOP閾值40。方案2根據(jù)歷史歷元中所有觀測(cè)值求取的ISB進(jìn)行改正。方案3為求取當(dāng)前歷元可見衛(wèi)星對(duì)應(yīng)的ISB，并改正當(dāng)前歷元相應(yīng)的觀測(cè)值，在求取ISB時(shí)，歷史歷元中所有GPS觀測(cè)值均參與解算。

表2為上述3種方案24 h定位結(jié)果?？梢钥闯觯桨?定位精度最高。

表2　定位偏差RMS統(tǒng)計(jì)

4　結(jié)　語(yǔ)

影響ISB的各項(xiàng)誤差源均可看作與衛(wèi)星編號(hào)相關(guān)，因此提出一種基于衛(wèi)星編號(hào)的ISB估計(jì)方案。以GPS為參考系統(tǒng)，利用不同類型接收機(jī)的觀測(cè)數(shù)據(jù)估計(jì)ISB，發(fā)現(xiàn)不同衛(wèi)星參與解算ISB估值存在明顯差異。最后，將顧及衛(wèi)星號(hào)和廣播星歷更新的ISB估計(jì)方案應(yīng)用于可見衛(wèi)星不足時(shí)的偽距單點(diǎn)定位，其精度明顯改善，且適用于多種類型接收機(jī)。

[1]徐龍威，劉暉，張明，等.不同截止高度角多模GNSS組合單點(diǎn)定位性能分析[J]. 大地測(cè)量與地球動(dòng)力學(xué),2015,35(6)：987-991(Xu Longwei, Liu Hui, Zhang Ming, et al. Performance Assessment of Multi-GNSS Single Point Positioning with Different Cut-off Elevation Angles[J]. Journal of Geodesy and Geodynamics,2015, 35(6):987-991)

[2]Torre A D, Caporali A. An Analysis of Intersystem Biases for Multi-GNSS Positioning[J]. GPS Solutions, 2015, 19(2): 297-307

[3]Montenbruck O, Hauschild A, Hessels U. Characterization of GPS/GIOVE Sensor Stations in the CONGO Network[J]. GPS Solutions, 2011, 15(3): 193-205

[4]Cai C S, Yang G. A Combined GPS/GLONASS Navigation Algorithm for Use with Limited Satellite Visibility [J].Journal of Navigation, 2009, 62(4): 671-685

[5]Angrisano A, Gaglione S, Gioia C. Performance Assessment of GPS/GLONASS Single Point Positioning in an Urban Environment[J]. Acta Geodaetica et Geophysica, 2013, 48(2): 149-161

[6]Montenbruck O, Steigenberger P, Hauschild A. Broadcast Versus Precise Ephemerides: A Multi-GNSS Perspective[J]. GPS Solutions, 2015, 19(2):321-333

[7]張雪. GNSS系統(tǒng)時(shí)間偏差監(jiān)測(cè)精度改善及評(píng)估[D].西安：中國(guó)科學(xué)院國(guó)家授時(shí)中心,2014(Zhang Xue. Research on the Accuracy Improvement and Assessment of GNSS System Time Offset Monitoring[D]. Xi’an：National Time Service Center, CAS, 2014)

Foundation support:National Key Research and Development Program of China, No.2016YFB0800405

About the first author:XU Longwei, PhD candidate, majors in multi-GNSS navigation and positioning, E-mail：xlw_ltu2012@163.com.

A Multi-GNSS Pseudorange Positioning Method with Inter System Bias

XULongwei1LIUHui1LIUYujie2SHUBao1ZHANGMing1QIANChuang1

1GNSS Research Center, Wuhan University, 129 Luoyu Road, Wuhan 430079, China 2Institute of Remote Sensing and Surveying and Mapping of Henan Province, 8 Huanghe Road, Zhengzhou 450003, China

In this paper, with the derivation of code observation equation, the structure and influence factors of inter system bias are analyzed. We discover that all factors for low precision navigation can be treated as related to satellite number. A new inter system bias parameter estimation solution is proposed. Finally, the new solution is applied to multi-GNSS positioning, providing reliable positioning information for real-time low precision users when only 4 multi-GNSS satellites are visible.

multi-GNSS; inter system bias; satellite number; positioning

2015-09-18

徐龍威，博士生，主要研究方向?yàn)槎嗄NSS導(dǎo)航定位，E-mail： xlw_ltu2012@163.com。

10.14075/j.jgg.2016.09.013

1671-5942(2016)09-0813-04

P228

A

項(xiàng)目來(lái)源：國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFB0800405)。