何倩倩，張 瑞，劉 楊

（1.武漢大學 衛(wèi)星導航定位工程技術(shù)研究中心，武漢 430079；2.武漢大學 測繪學院，武漢 430079）

1 引言

自選擇可用性 （selective availability，SA）政策取消之后，電離層延遲成為全球衛(wèi)星導航系統(tǒng)（global navigation satellite system，GNSS）測量中影響最嚴重的誤差源之一［1］。伽利略衛(wèi)星導航系統(tǒng) （Galileo navigation satellite system，Galileo）廣播星歷擬采用NeQuick模型，全球衛(wèi)星定位系統(tǒng) （global positioning system，GPS）廣播星歷采用Klobuchar模型，這些模型都是在全球范圍內(nèi)大尺度長周期平均的經(jīng)驗模型，電離層延遲修正精度不理想。而利用GNSS雙頻觀測值建立電離層延遲模型，能夠精確地估計電離層延遲，是改正電離層延遲誤差的有效方法。如今，國際上利用地基GPS數(shù)據(jù)建立電離層延遲模型已經(jīng)取得了重大進展。國際GNSS服務組織 （international GNSS service，IGS）自1998年成立以來，電離層工作組于2003年正式提供電離層垂直總電子含量（vertical total electron content，VTEC）服務；目前，美國噴氣推進實驗室 （jet propulsion laboratory，JPL）、歐洲歐洲定軌中心 （center for orbit determination in Europe，CODE）、歐洲航天局（European space agency，ESA）及西班牙加泰羅尼亞理工大學 （technical university of catalonia，UPC）等4家分析中心正在從事全球電離層延遲日常數(shù)據(jù)處理和產(chǎn)品發(fā)布。國內(nèi)，武漢大學、上海天文臺、中國科學院測量與地球物理研究所等單位也開展了相關(guān)研究。海量地基GPS雙頻觀測為電離層延遲建模研究提供了高分辨率時空覆蓋的數(shù)據(jù)源［2-3］。盡管如此，穿刺點的數(shù)量及空間分布、觀測精度影響著建模精度。聯(lián)合使用不同的導航系統(tǒng)進行電離層的監(jiān)測與預報研究，增加了測站可視衛(wèi)星的數(shù)量，并極大改善了穿刺點的幾何分布，更加有利于區(qū)域電離層模型的建立。因此，利用GNSS多頻觀測數(shù)據(jù)進行電離層的監(jiān)測與預報研究是當前的熱點問題。目前，GNSS系統(tǒng)尚未全面建成，要同時獲得GPS、格洛納斯衛(wèi)星導 航 系 統(tǒng) （global navigation satellite system，GLONASS）、Galileo、北斗衛(wèi)星導航系統(tǒng) （BeiD－ou navigation satellite system，BDS）等多系統(tǒng)的雙頻觀測值，只能通過仿真實現(xiàn)。本文基于Analytical Graphics公司開發(fā)的衛(wèi)星工具軟件包 （satellite tool kit，STK），利用軟件編程技術(shù)實現(xiàn)了GNSS多頻觀測數(shù)據(jù)的仿真，采用GNSS多頻觀測數(shù)據(jù)來構(gòu)建中國區(qū)電離層球諧函數(shù)模型，并對該區(qū)域模型結(jié)果進行了分析 。研究結(jié)果表明：利用GNSS多衛(wèi)星系統(tǒng)頻觀測值建立中國區(qū)域電離層球諧函數(shù)模型的精度要優(yōu)于單一系統(tǒng)。

2 GNSS多頻觀測數(shù)據(jù)

本文基于Analytical Graphics公司開發(fā)的衛(wèi)星工具軟件包STK，對GPS、GLONASS、Galileo及BDS等導航系統(tǒng)進行了星座仿真，并根據(jù)GNSS定位原理和數(shù)學模型，利用軟件編程技術(shù)實現(xiàn)了GNSS多頻觀測數(shù)據(jù)的仿真。相比于使用硬件模擬實現(xiàn)多模導航系統(tǒng)的仿真，本文使用的方法則具有技術(shù)簡單、成本較低及研究周期短的優(yōu)勢。

2.1 GNSS衛(wèi)星星座仿真

進行星座仿真時，采用的GPS、GLONASS、Galileo以及BDS基本性能及參數(shù)如表1所示。據(jù)此所得衛(wèi)星坐標文件則為2012年1月26日全天24h，以15min為采樣間隔，基于 WGS－84參考框架和GPS時時間系統(tǒng)的全部衛(wèi)星坐標信息，并將其轉(zhuǎn)換為標準SP3格式。

表1 GNSS系統(tǒng)仿真參數(shù)

2.2 GNSS觀測值仿真

GNSS實際觀測中，誤差來源主要包括電離層延遲、對流層延遲、衛(wèi)星星歷誤差、衛(wèi)星鐘差、接收機鐘差及噪聲及多路徑效應等，可以使用國際全球?qū)Ш叫l(wèi)星系統(tǒng)服務 （international GNSS service，IGS）提供的精密星歷、電離層、對流層、接收機鐘差等產(chǎn)品仿真各種誤差影響，也可以選擇合適的誤差模型直接仿真各種誤差影響。表2顯示了本文在進行GNSS觀測值仿真時所用的幾種主要的誤差模型及其精度。

表2 仿真誤差模型及其精度

根據(jù)衛(wèi)星導航系統(tǒng)的測量原理，偽距和載波相位定位的觀測方程分別為

式中，ρi為第i顆衛(wèi)星至接收機的距離，（Xi，Yi，Zi） 為 衛(wèi) 星 坐 標，（XR，YR，ZR） 為 接 收 機坐標，c為真空中光速，ΔtR為接收機鐘差，Δti為第i顆衛(wèi)星的鐘差，為第i顆衛(wèi)星信號的對流層延遲距離改正，為第i顆衛(wèi)星信號的電離層延遲距離改正，φi為第i顆衛(wèi)星的載波相位觀測值，λ為載波的波長，N為整周模糊度。衛(wèi)星坐標由仿真而得的SP3精密星歷，以30s間隔的Neville內(nèi)插法獲得，插值階數(shù)為10，接收機坐標則取自中國區(qū)域內(nèi)盡可能覆蓋全國的27個觀測站 （測站分布如圖1所示），接收機鐘差由高斯隨機模型模擬，衛(wèi)星鐘差設為0。

圖1 仿真實驗測站分布

2.3 GNSS仿真觀測值精度分析

對仿真觀測數(shù)據(jù)的驗證有很多方式，但對Galileo及BDS等系統(tǒng)來說，目前尚無法采集到足夠的可用的實際觀測數(shù)據(jù)。本文仿真的數(shù)據(jù)是為后續(xù)區(qū)域電離層模型的建立提供可用的輸入數(shù)據(jù)，所以只需要檢驗仿真觀測數(shù)據(jù)的可用性是否具有足夠的精度即可。本文選取部分站作為仿真參考站，用單點定位軟件處理仿真觀測數(shù)據(jù)，并與其實際坐標進行比較。

本文模擬了2012年1月26日1－23h中國區(qū)域內(nèi)27個GNSS觀測站的觀測數(shù)據(jù)，采樣率為30s。以哈爾濱 （HRBN）為參考站，分析仿真觀測數(shù)據(jù)的可用性，結(jié)果如表3所示。

表3 HRBN站解算坐標與參考站實際坐標對比

考慮到仿真中誤差處理還不夠全面，衛(wèi)星鐘差為零，對于接收機噪聲采用的處理方法也比較簡單，沒有考慮多路徑誤差，所以該仿真觀測數(shù)據(jù)可以滿足后續(xù)的研究需要。

3 利用GNSS雙頻觀測值建立中國區(qū)域電離層延遲模型

電離層分布在離地面約50－1 000km的高度范圍內(nèi)，在GNSS電離層研究領(lǐng)域常常引入單層假設［4］，在這一假設下，總電子含量成為與位置、時間有關(guān)的數(shù)學量，從而可以通過數(shù)學模型，如多項式模型［5］、球諧函數(shù)模型等進行擬合。在提供大尺度總電子含量 （total electron content，TEC）產(chǎn)品方面，CODE現(xiàn)已研制出全球TEC現(xiàn)報系統(tǒng)并實現(xiàn)了TEC現(xiàn)報結(jié)果的網(wǎng)上準實時發(fā)布，向用戶提供在日固坐標系和電離層單層模型下全球范圍的各時段的球諧函數(shù)展開式的系數(shù)，用戶則根據(jù)這些系數(shù)求得自己所需范圍TEC值［6-11］。本文將利用GNSS雙頻觀測值來建立中國區(qū)域電離層低階球諧函數(shù)模型。

球諧函數(shù)擬合法，是利用一系列基于電離層單層模型和日固坐標系下的球諧函數(shù)，來擬合各段時空范圍內(nèi)的TEC分布情況。在獲取了絕對電離層延遲信息之后，根據(jù)衛(wèi)星和測站坐標獲取可穿刺點的位置信息，并將其轉(zhuǎn)換到地磁日固坐標系下，據(jù)此則可以建立二維電離層模型的觀測方程［12］

式中，β為穿刺點的地磁緯度，s為穿刺點的日固經(jīng)度，為正則化勒讓德級數(shù)，anm、bnm為待估的系數(shù)，STEC（β，s）為從 GNSS觀測值 （此處特指P1P2型接收機，其他類型觀測值需劃歸到統(tǒng)一基準下）中獲取的絕對電離層延遲信息，M為對應的投影函數(shù)，DCBs為衛(wèi)星的P2－P1間硬件延遲偏差，DCBr為測站接收機的P2－P1硬件延遲偏差。IGS提供的產(chǎn)品一般是P1－P2硬件延遲偏差，所以與式 （3）觀測方程獲得的結(jié)果符號相反。從式（3）中可以看出，測站與衛(wèi)星的硬件延遲在沒有其他已知條件的情況下是無法分離的，要分離測站和衛(wèi)星的硬件延遲，通常需要附加一定的約束條件。因為所用的GNSS觀測值在仿真是并沒有加入測站和衛(wèi)星的硬件延遲，這里暫時不必考慮硬件延遲的分離。

4 多衛(wèi)星系統(tǒng)GNSS電離層建模分析

4.1 穿刺點分析

利用GNSS雙頻觀測值建模時，各穿刺點處上的VTEC作為觀測值，穿刺點的數(shù)量以及分布是影響建模精度的主要因素之一，穿刺點分布密集的區(qū)域擬合精度高，而穿刺點分布稀疏的區(qū)域擬合精度低。相比于單一的GPS、Galileo或者是BDS，多GNSS導航系統(tǒng)具有明顯的優(yōu)勢：它集合了各個導航定位系統(tǒng)，增加了可視衛(wèi)星的數(shù)量，改善了星座分布及其精度衰減因子 （dilution of precision，DOP）值，增加了穿刺點的數(shù)量，并改善了穿刺點的分布。圖2所示為2012年1月26日1時第一個歷元內(nèi)，中國區(qū)域的穿刺點分布情況。

4.2 多衛(wèi)星系統(tǒng)GNSS電離層建模精度分析

本文分別利用BDS、GPS及Galileo導航系統(tǒng)，以及GPS＋GLONASS＋BDS組合導航系統(tǒng)的雙頻觀測值，使用4階球諧函數(shù)模型，模擬了2012年1月26日中國區(qū)域內(nèi)的電離層狀況，輸出了在GPST時間系統(tǒng)下，2012年1月26日3－4h，7－8h，11－12h，15－16h，19－20h等5個時段的電離層延遲量，并與GNSS觀測值中的實際電離層延遲量進行對比。

圖2 中國區(qū)域IPP分布圖

從高緯度地區(qū)至低緯度地區(qū)，選取了包頭站（BAOT）、武漢站 （WUHN）及廣州站 （GUAN）三個典型的GNSS觀測站，采取以下策略進行對比分析與驗證：1）利用建模得到的1個多系統(tǒng)模型及3個單系統(tǒng)模型計算出的包頭站 （BAOT）、武漢站（WUHN）、廣州站 （GUAN）全天的VTEC信息投影到斜路徑方向，求得斜路徑電離層延遲信息；2）以GNSS觀測值誤差仿真時所加入的電離層延遲作為真值，視其為實際的電離層延遲量，減掉上述步驟求出的模型斜路徑絕對電離層延遲，將其差值的絕對值與實際延遲量相比。模型誤差在實際延遲量中所占的百分比越低，則模型精度越高；反之，則模型精度越低。

為了比較多/單導航系統(tǒng)觀測值建模的平均水平，使用了電離層模型誤差與實際延遲量的平均百分比。圖3至圖5則直觀地顯示和對比了位于高緯到低緯的測站，使用多/單導航系統(tǒng)觀測值建立區(qū)域電離層模型的平均水平。

圖3至圖5的結(jié)果顯示：相比于利用單系統(tǒng)數(shù)據(jù)建立電離層模型，GNSS多系統(tǒng)建模精度有一定提高，且位于中緯度地區(qū)的武漢站和低緯地區(qū)的廣州站優(yōu)于緯度相對較高的包頭站。初步分析原因如下：GNSS系統(tǒng)中，GPS有32顆衛(wèi)星，北斗系統(tǒng)為30顆星，衛(wèi)星數(shù)量較多，則提供的穿刺點數(shù)量就較多，利于提高球諧函數(shù)模型的精度；BDS因采用了3顆地球靜止軌道衛(wèi)星 （geostationary earth orbits，GEO）、3顆傾斜地球同步軌道衛(wèi)星（inclined geo－synchronous orbits，IGSO）及分布在3個軌道面內(nèi)的24顆、傾角為55°的中圓軌道衛(wèi)星 （medium earth orbit，MEO），使得IGSO能充分利用GEO的優(yōu)點，同時克服了其高緯度區(qū)始終是低仰角的問題，故使用BDS觀測值建立的中國區(qū)域電離層球諧函數(shù)模型精度較好，且很多時段其建模精度要優(yōu)于GPS觀測值的建模精度；GLONASS只有24顆星，衛(wèi)星數(shù)量少，使得電離層建模精度不高。另外，觀測值仿真僅采用了Klobuchar模型仿真了中國區(qū)域的電離層變化，然而，Klobuchar模型結(jié)構(gòu)簡單，無法精確反映中國區(qū)域電離層變化的復雜性，從而使得仿真得到的電離層變化相對較為平穩(wěn)，在該情況下，電離層建模的精度主要取決于所采用的電離層觀測值數(shù)據(jù)及其分布；而數(shù)據(jù)仿真中低緯度和中緯度地區(qū)的電離層穿刺點分布較多且相對較為均勻，這個非常有利于電離層建模，而在高緯度地區(qū)，特別是包頭附近，電離層穿刺點相對較少，且分布也不均勻，電離層建模精度勢必會降低。

圖3 包頭站區(qū)域模型誤差占實際延遲量的平均百分比

圖4 武漢站區(qū)域模型誤差占實際延遲量的平均百分比

圖5 廣州站區(qū)域模型誤差占實際延遲量的平均百分比

5 結(jié)束語

本文基于STK軟件模擬了BDS、GPS、Galileo以及GLONASS衛(wèi)星系統(tǒng)衛(wèi)星導航系統(tǒng)，并采用全國分布的27個測站，分析了單系統(tǒng)、多系統(tǒng)組合對建立中國區(qū)域電離層低階球諧函數(shù)模型的影響。仿真結(jié)果表明：對于中國區(qū)域，利用GNSS多衛(wèi)星系統(tǒng)觀測值建立中國區(qū)域電離層球諧函數(shù)模型的平均精度要優(yōu)于單一系統(tǒng)，我國BDS單系統(tǒng)觀測值的建模精度與GNSS多系統(tǒng)建模精度相近，模型誤差最大時僅有實際延遲量的14%，二者均優(yōu)于GPS單系統(tǒng)觀測值的建模精度，其模型誤差最大可達20%，而GLONASS單系統(tǒng)觀測值建模精度最差，模型誤差最大時達實際延遲量的35%。

［1］LIU Zhi－zhao.Ionosphere Tomographic Modeling and Applications Using Global Positioning System （GPS）Measurements［D］.Calgary：University of Calgary，2004.

［2］章紅平，韓文慧，黃玲，等.地基GNSS全球電離層延遲建模［J］.武漢大學學報：信息科學，2012，37（10）：1186－1189.

［3］邱蕾，羅和平，王澤民.GPS網(wǎng)絡RTK流動站的對流層內(nèi)插改正分析［J］.測繪工程，2011，20（5）：70－73.

［4］SCHAER S.Mapping and Predicting the Earth’s Ionosphere Using the Global Positioning System［D］.Bern：The University of Bern，1999.

［5］章紅平，施闖，唐衛(wèi)明.地基GPS區(qū)域電離層多項式模型與硬件延遲統(tǒng)一解算分析［J］.武漢大學學報：信息科學版，2008.33（8）：805－808.

［6］袁運斌.基于GPS的電離層監(jiān)測及延遲改正理論與方法的研究［D］.北京：中國科學院研究生院，2002.

［7］李衛(wèi)國.基于GPS觀測值構(gòu)建區(qū)域電離層延遲模型的研究［D］.合肥：合肥工業(yè)大學，2009.

［8］章紅平.基于地基GPS的中國區(qū)域電離層監(jiān)測與延遲改正研究［D］.上海：中國科學院上海天文臺，2006.

［9］韓玲.區(qū)域GPS電離層TEC監(jiān)測、建模和應用［D］.上海：中國科學院上海天文臺，2006.

［10］黃逸丹.區(qū)域電離層GPS監(jiān)測及應用研究［D］.上海：中國科學院上海天文臺，2007.

［11］何玉晶.GPS電離層延遲改正及其擾動監(jiān)測的分析研究［D］.鄭州：解放軍信息工程大學，2006.

［12］耿長江.利用地基GNSS數(shù)據(jù)實時監(jiān)測電離層延遲理論與方法研究［D］.武漢：武漢大學，2011.