任超超，李靈巧，2，蔡成林，楊輝華，2

(1.桂林電子科技大學(xué) 電子工程與自動(dòng)化學(xué)院，廣西 桂林　541004；2.北京郵電大學(xué) 自動(dòng)化學(xué)院，北京　100876；3.桂林電子科技大學(xué) 廣西精密導(dǎo)航技術(shù)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室，廣西 桂林　541004)

?

星歷誤差的低次曲面模型改正方法

任超超1，李靈巧1，2，蔡成林3，楊輝華1，2

(1.桂林電子科技大學(xué) 電子工程與自動(dòng)化學(xué)院，廣西 桂林541004；2.北京郵電大學(xué) 自動(dòng)化學(xué)院，北京100876；3.桂林電子科技大學(xué) 廣西精密導(dǎo)航技術(shù)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室，廣西 桂林541004)

針對廣播星歷誤差限制了廣域差分定位精度的問題，提出了星歷誤差的低次曲面模型改正方法。通過6個(gè)均勻分布的差分站與5顆由GEO、MEO和IGSO衛(wèi)星組成的導(dǎo)航星座對該方法進(jìn)行仿真驗(yàn)證，結(jié)果表明低次曲面模型法能將星歷誤差改正到較高的精度。同時(shí)，針對6個(gè)差分站改正效果均勻性較差的問題，提出了11個(gè)差分站優(yōu)化布局方案，使得低次曲面模型法的改正效果更加均勻，定位精度達(dá)到1 m以內(nèi)。該方法對廣域差分定位精度有良好的改正效果，而且在GEO、MEO和IGSO衛(wèi)星導(dǎo)航星座中具有通用性，對我國衛(wèi)星導(dǎo)航星基增強(qiáng)系統(tǒng)建設(shè)有一定的參考價(jià)值。

星歷誤差；廣域差分；曲面模型；優(yōu)化布局

0　引言

廣播星歷是由全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system，GNSS)的地面控制中心確定和提供的，經(jīng)衛(wèi)星向全球所有用戶公開播發(fā)的一種預(yù)報(bào)星歷，其精度較差[1]。目前全球定位系統(tǒng)(global positioning system，GPS)的廣播星歷誤差在3～5 m，廣播星歷誤差限制了GPS的定位精度[2]。文獻(xiàn)[3]提出了星鐘和星歷誤差分離的廣域差分方法，將星鐘和星歷分離解算，有效地降低了精度衰減因子(dilution of precision，DOP)值，提高了差分用戶的定位精度。本文在這種方法的基礎(chǔ)上提出一種星歷誤差的低次曲面模型改正方法：在空間上將主站(與衛(wèi)星時(shí)間嚴(yán)格同步)觀測的偽距誤差分配給星鐘作為星鐘誤差，通過一站式改正方法[3]消除；然后將外場多個(gè)副站(與衛(wèi)星時(shí)間不嚴(yán)格同步)觀測的偽距誤差扣除星鐘誤差后的剩余量分配給星歷作為星歷誤差，在一定范圍的服務(wù)區(qū)內(nèi)建立曲面模型；并將曲面模型的幾個(gè)參數(shù)廣播給用戶，用戶根據(jù)自己所處的經(jīng)緯度，利用曲面模型參數(shù)計(jì)算出所在位置的剩余星歷誤差并直接改正。

1　星歷誤差及其空間分布

在衛(wèi)星導(dǎo)航星基增強(qiáng)系統(tǒng)中，各差分站時(shí)間保持完全同步非常困難。當(dāng)前采用衛(wèi)星雙向時(shí)間頻率傳遞技術(shù)[4]，各站之間仍會(huì)有0.3～1 ns的時(shí)間同步誤差；如果采用理想情況的矢量差分原理解算，站鐘時(shí)間同步誤差會(huì)吸收為星歷誤差，受逆定位DOP的誤差放大[5]影響，星歷改正數(shù)的解算精度并不高，目前GPS的星歷改正精度只能達(dá)到2.5 m[6]。

為了獲得更高的差分解算精度，文獻(xiàn)[7]提出一種星鐘和星歷誤差相對改正方法，這種改正并不是將某顆衛(wèi)星的星歷誤差真正完全分離，而是將相對于基準(zhǔn)星的星歷誤差分配給這顆星，它實(shí)質(zhì)上是同一副站對不同衛(wèi)星的徑向星歷誤差的相對量。外場差分站的偽距誤差經(jīng)過電離層、對流層、偽碼噪聲平滑等誤差預(yù)處理[8-9]，再經(jīng)過星鐘誤差一站式改正后，可認(rèn)為是由星歷誤差和站鐘偏差2部分組成的，其中站鐘偏差可用星間單差法[10]消除，最后剩余量為星歷誤差。

然而星鐘誤差一站式改正法只能保證主控站視向星歷誤差消除，對于非主控站，星歷誤差會(huì)隨著主控站和用戶之間距離的增加而增加。星歷誤差最終會(huì)反映到用戶的測量結(jié)果上；因此可以將星歷誤差投影到站星矢徑上進(jìn)行分析，用星歷誤差的用戶等效距離誤差(user equivalent range error，UERE)來表示。本文進(jìn)一步分析了4顆衛(wèi)星星歷誤差的UERE空間分布規(guī)律，如圖1所示。

圖1　星歷誤差的UERE空間分布

衛(wèi)星1、4、5為中地球軌道(medium Earth orbit，MEO)衛(wèi)星，衛(wèi)星2為地球靜止軌道(geosynchronous Earth orbit，GEO)衛(wèi)星，衛(wèi)星3為傾斜地球同步軌道(inclined geosynchronous satellite orbit，IGSO)衛(wèi)星；其中衛(wèi)星1作為基準(zhǔn)星不改正，主控站為西安站(110°E，34°N)。從圖1中可以看出：1)在74°E～138°E與9°N～49°N服務(wù)區(qū)內(nèi)對不同的衛(wèi)星，與主控站處于同一緯度但經(jīng)度不同的用戶，其徑向星歷誤差的UERE在1.5 m以內(nèi)，與主控站處于同一經(jīng)度但緯度不同的用戶，其徑向星歷誤差的UERE在1 m以內(nèi)；2)主控站附近星歷誤差最小，距離主控站越遠(yuǎn)星歷誤差越大，服務(wù)區(qū)內(nèi)星歷誤差的大小與用戶所處的經(jīng)緯度有關(guān)；3)在74°E～138°E與9°N～49°N服務(wù)區(qū)內(nèi)，徑向星歷誤差在空間上呈現(xiàn)一個(gè)平滑的曲面形狀。

2　低次曲面模型法原理

從圖1給出的星歷誤差的空間分布可以看出，星歷誤差在空間上能用低次曲面近似表征?；谶@一點(diǎn)，本文提出利用低次曲面模型改正星歷誤差：利用多個(gè)差分站的星歷誤差，在一定的服務(wù)區(qū)內(nèi)建立低次曲面模型；并將曲面模型的幾個(gè)參數(shù)廣播給用戶，用戶根據(jù)自己所處的經(jīng)緯度，利用曲面模型參數(shù)計(jì)算出所在位置的星歷誤差并直接改正。

低次曲面模型可表征為：

(1)

(2)

(3)

a4B2+a5L3+a6B3；

(4)

a4LB+a5B2+a6L3+a7B3；

(5)

a5B2+a6L3+a7BL2+a8LB2+a9B3。

(6)

式(1)至(6)中：V表示星歷誤差；L和B分別表示經(jīng)度和緯度；ai(i=0，1，2，…,9)為第i個(gè)系數(shù)陣，迭代解為

(7)

式中Δρ表示一站式相對改正消除星鐘誤差后的剩余星歷誤差，其模型對應(yīng)的迭代矩陣H有6種形式，式(8)以模型3為例表示為

(8)

式中：H的行數(shù)等于差分站數(shù)量；Ln和Bn分別表示差分站n的經(jīng)度和緯度。

根據(jù)模型參數(shù)計(jì)算用戶所在位置的星歷誤差為

(9)

式中Hu表示用戶所在位置對應(yīng)的迭代矩陣，有6種形式，式(10)以模型3為例表示為

(10)

式中Lu和Bu分別表示用戶所在位置的經(jīng)度和緯度。

改正后的殘余星歷誤差

ε=Δρ-V。

(11)

3　實(shí)驗(yàn)與結(jié)果分析

3.1曲面模型改正效果

為了驗(yàn)證上述6種低次曲面模型的改正效果，本文利用6個(gè)均勻分布的差分站(西安、長春、上海、烏魯木齊、三亞、拉薩)與5顆由GEO、MEO和IGSO衛(wèi)星組成的混合星座進(jìn)行實(shí)驗(yàn)。衛(wèi)星1、4、5選用GPS中的MEO衛(wèi)星(PRN7、PRN8、PRN10)，衛(wèi)星2、3分別選用北斗衛(wèi)星導(dǎo)航系統(tǒng)中的GEO衛(wèi)星(PRN3)和IGSO衛(wèi)星(PRN9)；這5顆衛(wèi)星在北京時(shí)間04∶00—05∶00對這6個(gè)差分站都是可視的，其中衛(wèi)星1作為基準(zhǔn)星不改正。實(shí)驗(yàn)數(shù)據(jù)中衛(wèi)星1、4、5的真實(shí)星歷用國際GNSS服務(wù)組織發(fā)布的精密星歷代替，衛(wèi)星2、3的真實(shí)星歷用武漢大學(xué)GNSS研究中心發(fā)布的精密星歷代替，廣播星歷以及偽距觀測值來自中科院國家授時(shí)中心的雙頻接收機(jī)，采樣頻率為5min，通過內(nèi)插讓真實(shí)星歷和廣播星歷均為1s一次。在74°E～138°E與9°N～49°N格網(wǎng)范圍內(nèi)，偽距誤差經(jīng)過電離層誤差改正、對流層誤差改正、偽碼噪聲平滑、星鐘誤差一站式改正、站鐘偏差消除后，剩余量為星歷誤差，通過曲面模型法在經(jīng)度和緯度方向每隔4°計(jì)算一個(gè)改正后的殘余星歷誤差ε，用ε的最大值Max、最小值Min和平均數(shù)E表征模型的精度，用ε的標(biāo)準(zhǔn)差(來表征模型改正效果的均勻性。實(shí)驗(yàn)結(jié)果如表1所示。

從表1可以看出：1)對GEO、MEO和IGSO這3種混合星座，曲面模型改正后的殘余星歷誤差的均值都小于0.1 m，誤差的標(biāo)準(zhǔn)差都小于0.1 m，這表明低次曲面模型法達(dá)到了較高的精度；2)模型3改正后GEO、MEO和IGSO這3種混合星座衛(wèi)星的殘余星歷誤差的最大值、最小值、平均數(shù)和標(biāo)準(zhǔn)差的絕對值都是最小的，這表明模型3的改正效果最好，其中，殘余星歷誤差小于0.15 m，誤差的均值小于0.06 m，誤差的標(biāo)準(zhǔn)差小于0.06 m。圖2給出模型3改正后的殘余星歷誤差空間分布，可以直觀地看出模型3的改正效果。

表1　曲面模型改正后的殘余星歷誤差　m

圖2　模型3改正后的殘余星歷誤差空間分布

3.2差分站優(yōu)化布局

利用6個(gè)均勻分布的差分站建立的低次曲面模型3雖然能將星歷誤差的均值改正到0.06 m；但是誤差的標(biāo)準(zhǔn)差大于均值，殘余星歷誤差的最大值達(dá)到0.15 m，說明改正效果不夠均勻：因此需要考慮差分站優(yōu)化布局[8-10]。

從圖2可以看出，由6個(gè)均勻分布的差分站建立的低次曲面模型改正后的視向殘余星歷誤差較大的區(qū)域大約在(100°E～120°E，20°N～40°N)和(60°E～80°E，10°N～20°N)之間；因此考慮在包頭、喀什、長沙、昆明和東沙群島各增加1個(gè)差分站(對這5顆衛(wèi)星可視)，由11個(gè)差分站優(yōu)化布局建立低次曲面模型，選擇改正效果最好的模型3進(jìn)行仿真實(shí)驗(yàn)。實(shí)驗(yàn)結(jié)果如表2所示。

表2　差分站優(yōu)化布局后的曲面模型改正結(jié)果　m

從表2可以看出：對GEO、MEO和IGSO這3種混合星座，11個(gè)差分站優(yōu)化布局建立的低次曲面模型能將星歷誤差改正到0.09 m以內(nèi)，誤差的均值在0.05 m以內(nèi)，誤差的標(biāo)準(zhǔn)差小于0.02 m。圖3給出差分站優(yōu)化布局后的曲面模型改正效果，可以直觀地看出曲面模型的改正精度得到很大提升，改正效果非常均勻。

圖3　差分站優(yōu)化布局后的曲面模型改正效果

圖4用靶圖給出11站差分前后的主站(西安站)、副站(以長沙站為代表)和2個(gè)網(wǎng)內(nèi)驗(yàn)證用戶(武漢和成都)的定位偏差。靶心表示測繪位置；與靶心的距離表示定位偏差?？梢钥闯觯翰罘智暗亩ㄎ黄钤?0～30 m之間；經(jīng)過電離層、對流層、偽碼噪聲平滑、星鐘誤差一站式改正、星歷誤差的低次曲面模型改正后的定位偏差在1 m以內(nèi)。

圖4　定位偏差

4　結(jié)束語

通過建立低次曲面模型與差分站優(yōu)化布局，可以得出如下結(jié)論：

1)經(jīng)過星鐘誤差一站式改正后，星歷誤差能夠用本文提出的低次曲面模型法改正，其中二次曲面模型的改正效果最好。

2)低次曲面模型法的改正效果與地面差分站布局有很強(qiáng)的相關(guān)性，本文提出的11個(gè)差分站優(yōu)化布局方案建立的曲面模型可以將中國陸地和近海區(qū)域的徑向星歷誤差改正到0.09 m以內(nèi)，誤差的均值在0.05 m以內(nèi)，誤差的標(biāo)準(zhǔn)差小于0.02 m，改正效果非常均勻，差分定位偏差修正到1 m以內(nèi)。

3)本文提出的方法在MEO、GEO和IGSO衛(wèi)星導(dǎo)航星座中具有通用性，對我國衛(wèi)星導(dǎo)航星基增強(qiáng)系統(tǒng)建設(shè)有一定的參考價(jià)值。

[1]焦海松，王紅芳，姚飛娟.衛(wèi)星星歷誤差對 GPS 定位精度的影響與分析[J].全球定位系統(tǒng)，2009，34(1)：24-28.

[2]周立. GPS 廣播星歷誤差影響診斷與預(yù)測模型研究[J].測繪科學(xué)，2005，30(4)：37-38，4.

[3]CAI Chenglin，LI Xiaohui，WU Haitao. An analysis of the wide area differential method of geostationary orbit satellites[J]. Science in China Series G：Physics，Mechanics and Astronomy，2009，52(2)：310-314.

[4]BAUNCH A. Time and frequency comparisons using radiofrequency signals from satellites[J]. Comptes Rendus Physique，2015，16(5)：471-479.

[5]王澤民，孟泱，伍岳，等. GPS、Galileo 及其組合系統(tǒng)導(dǎo)航定位的DOP值分析[J].武漢大學(xué)學(xué)報(bào)·信息科學(xué)版，2006，31(1)：9-11，38.

[6]曹新運(yùn)，王堅(jiān)，劉超. DOP值進(jìn)行 GNSS 相對定位預(yù)報(bào)的適用性研究[J].大地測量與地球動(dòng)力學(xué)，2013，33(1)：81-85.

[7]CAI Chenglin，LI Xiaohui，WU Haitao. A virtual reference satellite differential method on relative correction of satellite ephemeris errors[J].Science in China Series G：Physics，Mechanics and Astronomy，2010，53(12)：2300-2306.

[8]王甫紅.高精度星載GPS實(shí)時(shí)定軌卡爾曼濾波模型[J].武漢大學(xué)學(xué)報(bào)·信息科學(xué)版，2010，35(6)：653-656.

[9]易重海，陳永奇，朱建軍，等.一種基于 IGS 超快星歷的區(qū)域性實(shí)時(shí)精密單點(diǎn)定位方法[J].測繪學(xué)報(bào)，2011，40(2)：226-231.

[10]劉吉華，歐吉坤，孫保琪，等.基于星間單差法的GEO衛(wèi)星精密定軌[J].武漢大學(xué)學(xué)報(bào)·信息科學(xué)版，2011，36(1)：24-28.

Low-order surface model correction for ephemeris errors

REN Chaochao1，LI Lingqiao1，2，CAI Chenglin3，YANG Huihua1，2

(1.School of Electronic Engineering and Automation，Guilin University of Electronic Technology，Guilin，Guangxi 541004，China；2.School of Automation，Beijing University of Posts and Telecommunications，Beijing 100876，China；3.Guangxi Key Laboratory of Precision Navigation Technology and Application，Guilin University of Electronic Technology，Guilin，Guangxi 541004，China)

Aiming at the problem that the positioning accuracy of wide area differential system is limited by ephemeris errors，a method of low-order surface model correction for ephemeris errors was presented in this paper. Six evenly distributed differential stations and five satellites comprises GEO，MEO and IGSO were used in simulations. The results showed that this correction could achieve a high precision. Meanwhile，the positioning accuracy was corrected to one meter and the effect was more evenly by eleven optimized distribution differential stations. This method would have not only a high positioning accuracy but also versatility in a variety of GEO，MEO and IGSO satellite navigation constellation，and provide a reference for the construction of China’s satellite-based augmentation system of satellite navigation.

ephemeris errors；wide area differential；surface model；optimized distribution

2015-12-18

國家自然科學(xué)基金項(xiàng)目(61263028)；廣西自然科學(xué)基金重點(diǎn)項(xiàng)目(2012GXNSFDA053027，2014GXNSFBA118280)；中國科學(xué)院時(shí)間頻率基準(zhǔn)重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(Y000YR1S01)；桂林電子科技大學(xué)研究生教育創(chuàng)新計(jì)劃資助項(xiàng)目(GDYCSZ201478)。

任超超(1991—)，男，安徽蕪湖人，碩士研究生，研究方向?yàn)樾l(wèi)星導(dǎo)航與信息處理。

10.16547/j.cnki.10.-1096.20160301.

P228

A

2095-4999(2016)03-0001-06

引文格式：任超超，李靈巧，蔡成林，等.星歷誤差的低次曲面模型改正方法[J].導(dǎo)航定位學(xué)報(bào)，2016，4(3)：1-5,25.(REN Chaochao，LI Lingqiao，CAI Chenglin，et al.Low-order surface model correction for ephemeris errors[J].Journal of Navigation and Positioning，2016，4(3)：1-5,25.)