GNSS廣播電離層全球與區域模型分析

楊禮朋 楊德宏

摘要：

廣播電離層精度監測評估是GNSS導航信息的重要組成部分。對GNSS單頻用戶而言，只有利用模型才能對電離層時延進行改正。GPS和Galileo已實現全球服務，而BDS目前還只能提供區域服務。由于3種模型所采用的電離層時延改正模型各異，所以為了客觀評估幾大導航系統的模型性能，也為了后續對模型和參數的進一步優化，需要選擇可靠的基準，對播發的電離層參數進行實時評估。考慮到雙頻改正模型的高精度特性，以雙頻改正模型得到的值作為基準，在全球和區域范圍選取測站數據，對電離層延遲改正模型從改正精度和改正率兩方面進行評估。研究結果可為廣播電離層監測評估以及單頻導航定位等應用提供參考。

關鍵詞：

電離層延遲；克羅布歇模型；全球定位系統；北斗導航系統

DOIDOI：10.11907/rjdk.172774

中圖分類號：TP301

文獻標識碼：A文章編號文章編號：16727800（2018）003005904

英文摘要Abstract：Broadcast ionosphere accuracy monitoring and evaluation is an important part of GNSS navigation information. For GNSS singlefrequency users， only the model is used to correct the ionospheric delay. GPS and Galileo have achieved global services， and BDS is still a regional service. In order to objectively evaluate the model performance of several navigation systems， and to optimize the model or parameters， we need to select a reliable benchmark for the proposed ionospheric parameters. Realtime evaluation. Considering the highprecision characteristics of the dualfrequency correction model， this paper evaluates the performance of the model from the global and regional selection stations by using the values obtained by the dualfrequency correction model. The ionospheric delay correction is corrected from the correction accuracy and the correction rate. The research results can provide reference for the application of radio and ionospheric monitoring and evaluation and single frequency navigation and positioning applications.

英文關鍵詞Key Words：ionospheric delay； KLOBUCHAR； GPS； BDS

0引言

電離層是地表大氣層被許多太陽光射線電離的部分，距地表60～1 000km。在電離層中，許多中性分子被電離，形成正離子和自由電子。這些粒子對導航電波的傳播將產生很大影響，在天頂方向可達十幾米[1]。因此，需要采用一定方法消除或減弱電離層延遲的影響。目前常用的電離層延遲改正方法有：經驗模型、雙頻改正模型和實測數據模型[2]。具體介紹如下：①經驗模型有本特（bent）模型、國際參考電離層（IRI）模型、KLOBUCHAR模型；②雙頻改正模型即利用全球衛星導航系統（Global Navigation Satellite System，GNSS）雙頻數據直接計算電離層延遲或組成無電離層延遲的組合觀測量（其電離層延遲改正效果一般可達99%），但對于大多數單頻接收機用戶而言并不適用；③實測數據模型是利用GNSS雙頻數據建立的滿足實時用戶短期預報需求的模型[3]，常用的有IGS的電離層格網模型和CODE的球諧函數模型[4]。對于單頻用戶而言，廣播電離層延遲改正模型是削弱電離層延遲的主要方法，我國的BDS（BEIDOU Navigation Satellite System，BDS）和美國的全球定位系統（Global Positioning System，GPS）在導航電文中都是提供KLOBUCHAR8參數的單層電離層改正模型（后面簡稱K8），除此之外還有伽利略衛星導航系統（Galileo Navigation Satellite System GAL）發播NeQuick2模型參數、歐洲定軌中心（Center for Orbit Determination in Europe CODE）的全球電離層模型等[56]。

1GNSS廣播電離層模型

1.1GPS KLOBUCHAR

克羅布歇模型是一種被單頻用戶接收機廣泛采用的模型，它是本特模型的演化，該模型假設將整個電離層壓縮在一個單層上，距地表350km[7]。夜晚的延遲為常數，白天則是余弦函數取正的部分[2]。此方式基本反映了電離層時延的周天特性。同時，大量研究資料表明，采用此模型將每天最大的時延時刻確定在14時左右。

在天頂方向調制的L1載波（f=1 575.42MHz）上測距碼的電離層改正時延Tg可表示為：

Tg=5*10-9+Acos2πp（t-14h）（1）

振幅A和周期P分別為：

A=∑3i=0αi（φm）i（2）

P=∑3i=0βi（φm）i（3）

式中，αi和βi從廣播星歷中獲取。

GPSK8模型是基于地磁坐標系的模型，通過映射函數將天頂方向的總電子含量VTEC投影至傳播方向[8]。模型中所涉及的8個參數是利用全球GNSS雙頻觀測數據解算得到的，該參數每天更新一組，通過廣播星歷播發給用戶。該模型的優點在于公式簡單、計算效率高。

BDS目前采用的廣播電離層模型同GPS一樣，也是8參數[910]。兩種模型不同之處在于GPS是日固地磁坐標系，而BDS是地理坐標系，參數每兩個小時播發一組，現階段北斗廣播電離層已實現區域服務，但BDS的全球廣播電離層模型尚未公布。

1.2Galileo NeQuick2模型

NeQuick2模型是一種描述電離層電子密度時空變化的半經驗電模型[11]，其模型特點是把電離層分層，通過模型計算出給定時間和位置的相應電子密度。隨著伽利略系統的逐漸完善，目前伽利略系統中所用的模型是第二版模型，簡稱NeQuick2。新、舊模型之間的差距在于改變了底部與頂部電離層解析式以及相關變量的計算方法。

伽利略模型參數的輸入信息分為兩部分，一部分是地磁緯度文件和系數文件，另一部分是地理位置、時間信息等數據[12]。伽利略分布的全球監測站通過傾斜電子含量來優化NeQuick2中的有效電離參數AZ，該模型參數每天更新一次。具體算法詳見參考文獻[1]、[1315]。伽利略全球分布站如圖1所示。

1.3評估策略

本文主要評估了包括BDS Klobuchar（以下簡稱BDSK8）、GPS Klobuchar（以下簡稱BDSK8）和GALILEIO NEQUICK模型（本文采用了最新的NEQUICK2模型進行評估）在內的3種廣播電離層模型，評估的基本思路如下：

（1）利用武漢大學IGS數據中心下載所需的測站觀測數據以及brdc和brdm兩種星歷文件；利用伯爾尼軟件中提供的網址下載CODE分析中心的DCB產品，包括P1P2*_ALL.DCB和P1C1*.DCB兩個文件。其中P1P2*_ALL.DCB既包含了衛星的DCB，也包含了部分測站的DCB信息，為提取實測數據作準備。

（2）利用IONCFG程序在指定路徑上輸入測站的觀測數據、brdm文件和DCB產品文件，提取得到實測電離層信息，其結果文件如例（bjfs3140.16i）所示；經過程序轉換得到評估電離層模型的基礎數據，其結果文件如例（bjfs314.16）所示。

（3）利用編寫的Fortran程序提取3種模型的廣播電離層參數。對于GPSK8而言，利用Fortran程序提取brdc星歷文件中頭文件的8個參數即可；對于BDSK8而言，同樣只需提取brdm星歷文件中頭文件的8個參數；對于GAL而言，可以利用MGEX測站數據信息，在測站文件中提取頭文件的3個參數即為GAL模型參數。

（4）將第三步提取的電離層信息放入指定文件目錄下，利用實測格網，程序將第二步得到的基礎數據進行解算，得到測站的廣播電離層改正率和改正精度，并進行對比分析和評估。

2算例與分析

本文從全球選取了33個IGS站，其中南半球有11個測站（UNSA、SANT、MGUE、RIO2、ZAMB、WIND、HRAO、SUTV、MOBS、KARR、TOW2、SAND）。北半球23個站（SHAO、TWTF、URUM、ULAB、YAKT、NOVM、IISC、ADIS、BSHM、ZECK、METS、RABT、BOGT、UNBJ、SCH2、THU3、FLIN、YELL、QUIN、FAIV、BADG、USUD、TWTF），選取2017年年積日32～59天4周的測站數據進行評估分析。利用GMT軟件與測站信息畫出所選測站位置和穿刺點信息，如圖2所示。

為了統計廣播電離層模型精度的評估結果，利用修正比例和改正精度評估3種廣播電離層延遲改正模型的性能。

PER=∑ni=11-VTECmod-VTECrefVTECref（4）

RMS=2∑ni=1（VTECmod-VTECref）2n （5）

式中，VTECmod是利用程序解算得到的值，而VTECref是基準值。

2.13種模型評估分析

利用基準VTEC和廣播電離層模型VTEC進行統計。由于BDSK8是區域服務，因此選取國內的3個IGS測站進行分析，3種模型精度對比如圖3～圖8所示。

為了直觀反映中國及其周邊區域3種模型的改正比例和改正精度，本文選取范圍為20°N-60°N，70°E-135°E的區域。此區域有7個IGS測站（SHAO、URUM、TWTF、BADG、USUD、ULAB、NOVM）。統計測站每天的改正率，然后取其平均值得到各系統在中國及其周邊區域的改正率，如表1、表2所示。

從3種廣播電離層延遲改正模型在中國區域的平均改正率P以及RMS，可以看出：

（1）在中國區域，BDS廣播電離層模型精度高于 GPS、Galileo，改正率平均值為73%，RMS優于6個TECU，并且表現出良好的一致性。

（2）GPSK8模型是3種模型中精度最低的，平均改正率低于65%。

2.2全球模型南北半球評估分析

目前GPS和GAL已實現了全球服務，為了對比兩種模型在南北半球的精度，選取2016年11月30天的數據進行分析。按照白天和夜晚分別統計，并以測站的改正精度TECU值的大小作為評估標準來衡量。

伽利略模型在南半球解算時由于UNSA測站改正精度較低，其整體改正精度大于17個TECU，所以在統計時予以剔除。

GPS模型在南半球解算時UNSA和RIO2兩個測站改正精度大于17個TECU，在統計時予以剔除。

由表5、表6可知，在北半球區域，伽利略模型的精度高于GPS的模型精度。根據兩種模型在南北半球的對比，在北半球伽利略的模型精度要優于GPS，在南半球GPS要優于伽利略。同時，對同一種模型而言，北半球的改正精度優于南半球的精度。

3結語

本文主要開展了GNSS電離層延遲改正模型的評估工作，包括在中國區域和全球范圍內的評估。首先基于雙頻電離層延遲改正模型精確提取電離層信息，然后以此為基準，開展廣播電離層延遲改正模型的評估工作，主要結論如下：

（1）在中國及其周邊區域，BDSK8模型的改正率優于GPSK8和NeQuick2模型，其改正精度好于6個TECU。

（2）對比兩種模型在南北半球的精度可以看出，無論是GPS還是GAL，北半球的改正精度均優于南半球。

（3）在全球模型評估中，在北半球區域，GAL模型比GPSK8模型體現出了更好的精度優勢，GAL的整體改正精度無論是在白天還是夜晚都比GPSK8小很多。

（4）在全球模型評估中，在南半球區域，GPSK8模型在白天的精度與GAL相當，改正精度在1個TECU范圍內。整體而言，GPSK8的精度稍高于GAL，但兩者并未表現出顯著的精度優勢。

本文給出了雙頻改正模型評估廣播電離層延遲改正的基本思路及相關分析結果，對廣播電離層延遲改正模型的監測評估及單頻導航定位應用都具有一定參考價值。

參考文獻參考文獻：

[1]張強，趙齊樂.北斗衛星導航系統Klobuchar模型精度評估[C].第四屆中國衛星導航學術年會論文集，2014.

[2]李征航，張小紅.衛星導航定位新技術及高精度數據處理方法[M].武漢：武漢大學出版社，2009.

[3]魏傳軍.基于地基GNSS觀測數據的電離層延遲改正研究[D].西安：長安大學，2014.

[4]耿長江.利用地基GNSS數據實時監測電離層延遲理論與方法研究[D].武漢：武漢大學，2011.

[5]YUAN Y， J OU.Ionospheric Eclipse Factor Method （IEFM） for determining the ionospheric delay using GPS data[J]. Progress In Natural Science， 14（9）：800804.

[6]YUAN Y， C C TSCHERNING， P KNUDSEN，et al.The ionospheric eclipse factor method （IEFM） and its application to determining the ionospheric delay for GPS[J]. Journal of Geodesy， 82（1）：18.

[7]章紅平.基于地基GPS的中國區域電離層監測與延遲改正研究[D].上海：中國科學院上海天文臺，2006.

[8]李征航，黃勁松.GPS測量與數據處理[M].武漢：武漢大學出版社.2010.

[9]李子申.GNSS/COMPASS電離層時延修正及TEC監測理論與方法研究[D].武漢：中國科學院測量與地球物理研究所，2012.

[10]張強，趙齊樂.北斗衛星導航系統電離層模型精度的研究[C].第四屆北斗衛星導航學術年會，2013.

[11]謝杰，張博，侯博，等.Galileo衛星導航系統廣播 NEQUICK模型及其參數擬合[J].空間科學學報，2012，32（6）：881886.

[12]WU XIAOLI， HU XIAOGONG， WANG GANG.Evaluation of compass ionospheric model in GNSS positioning[J].Advances in Space Research，2013，51（6）：959968.

[13]李施佳，蔡昌盛，戴吾蛟，等.NEQICK電離層改正模型的性能評估[J].測繪工程，2015，24（1）：4043.

[14]李懷建，羅曉武，郭大軍.NEQUICK模型算法研究[J].第三屆中國衛星導航學術年會，2012.

[15]薛軍琛.中國地區電離層延遲函數模型建立與精度估計[D].青島：山東科技大學，2010.

責任編輯（責任編輯：黃健）