基于擬合推估的中國(guó)區(qū)域電離層延遲精確建模方法

袁運(yùn)斌，李子申，王寧波，4，張寶成，霍星亮

基于擬合推估的中國(guó)區(qū)域電離層延遲精確建模方法

袁運(yùn)斌1，3，李子申1，2，王寧波1，3，4，張寶成1，5，霍星亮1，3

(1.中國(guó)科學(xué)院 測(cè)量與地球物理研究所，武漢 430077；2.中國(guó)科學(xué)院 光電研究院，北京 100094； 3.大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，武漢 430077；4.中國(guó)科學(xué)院大學(xué)，北京 100049；5.澳大利亞科廷大學(xué)，珀斯 6845)

高精度的格網(wǎng)電離層是廣域增強(qiáng)系統(tǒng)建設(shè)和應(yīng)用的重要信息之一，針對(duì)我國(guó)區(qū)域電離層復(fù)雜多變的特點(diǎn)，本文將電離層延遲變化分為趨勢(shì)項(xiàng)和隨機(jī)項(xiàng)兩部分，推導(dǎo)了利用擬合推估原理進(jìn)行中國(guó)區(qū)域電離層延遲建模的基本方法；在此基礎(chǔ)上，利用電離層活動(dòng)高年、中年和低年的觀測(cè)數(shù)據(jù)對(duì)方法的精度進(jìn)行了詳細(xì)評(píng)估。試驗(yàn)結(jié)果表明：基于擬合推估的電離層延遲建模方法的修正精度在電離層活動(dòng)高年約為0.4～0.7 m，在電離層活動(dòng)低年約為0.2 m；對(duì)于我國(guó)中緯度地區(qū)，基準(zhǔn)站間距可適當(dāng)放大至800 km左右，而在我國(guó)低緯度地區(qū)基準(zhǔn)站間距應(yīng)保持在350～450 km之間。

電離層延遲，全球?qū)Ш叫l(wèi)星系統(tǒng)，擬合推估，北斗廣域增強(qiáng)系統(tǒng)

0 引言

實(shí)時(shí)精密定位服務(wù)已成為全球衛(wèi)星導(dǎo)航系統(tǒng)發(fā)展與應(yīng)用的重要方向之一。廣域增強(qiáng)系統(tǒng)作為拓展衛(wèi)星導(dǎo)航技術(shù)應(yīng)用的重要基礎(chǔ)設(shè)施，可有效地提高大范圍區(qū)域內(nèi)導(dǎo)航定位的精度和可靠性[1-4]。廣域增強(qiáng)系統(tǒng)播發(fā)的電離層延遲修正信息，不僅可滿足系統(tǒng)服務(wù)范圍內(nèi)單頻用戶的電離層延遲修正，還可有效地輔助雙頻/多頻用戶精密定位的快速初始化[5-8]。

由于電離層活動(dòng)具有顯著的區(qū)域特性，并且與太陽(yáng)活動(dòng)水平密切相關(guān)，因此，廣域增強(qiáng)系統(tǒng)電離層延遲修正方法通常需要結(jié)合其服務(wù)區(qū)上空的電離層活動(dòng)特點(diǎn)而專門(mén)設(shè)計(jì)[5，9-11]。現(xiàn)有廣域電離層延遲修正方法可分為兩類：格網(wǎng)法和建模法[6-7，12-15]。格網(wǎng)法通過(guò)加權(quán)內(nèi)插技術(shù)，可有效描述局域電離層活動(dòng)特點(diǎn)，但是對(duì)于無(wú)觀測(cè)覆蓋的格網(wǎng)點(diǎn)，通常只能依賴于精度較低的電離層延遲經(jīng)驗(yàn)值；建模法盡管可提供區(qū)域內(nèi)任意處的電離層延遲信息，但難以有效顧及局部電離層變化特點(diǎn)。就我國(guó)廣域增強(qiáng)系統(tǒng)的服務(wù)區(qū)域而言，緯度跨越近50°，經(jīng)度跨越近70°，特別是部分低緯地區(qū)還受到電離層“赤道異?！庇绊懀煌暥葞?nèi)電離層變化差異顯著[16]，現(xiàn)有的電離層建模方法難以直接應(yīng)用于我國(guó)區(qū)域?qū)崿F(xiàn)高精度的電離層延遲修正。

國(guó)內(nèi)外學(xué)者針對(duì)不同區(qū)域的電離層活動(dòng)特點(diǎn)，提出了多種針對(duì)廣域增強(qiáng)系統(tǒng)電離層延遲修正方法[9，11，17-23]；特別是，初步論證了基于Kriging插值方法建立廣域增強(qiáng)系統(tǒng)電離層格網(wǎng)的思路，在一定意義上將加權(quán)內(nèi)插法與整體建模法進(jìn)行了初步的融合[24-27]。但是，該方法在實(shí)際實(shí)施中僅利用簡(jiǎn)單的平面函數(shù)擬合區(qū)域電離層變化趨勢(shì)項(xiàng)，并固定了描述電離層變化相關(guān)性的協(xié)方差函數(shù)，使得其僅能夠滿足中緯度地區(qū)格網(wǎng)電離層延遲修正的需求。針對(duì)此，中國(guó)科學(xué)院測(cè)量與地球物理研究所進(jìn)一步發(fā)展了基于擬合推估理論的廣域增強(qiáng)系統(tǒng)電離層延遲修正方法。該方法利用球諧函數(shù)描述我國(guó)區(qū)域電離層變化趨勢(shì)項(xiàng)，通過(guò)建立自適應(yīng)的協(xié)方差函數(shù)描述區(qū)域電離層變化的隨機(jī)項(xiàng)，精確估計(jì)格網(wǎng)點(diǎn)電離層延遲及其精度信息[5，10，28]。本文正是在總結(jié)上述相關(guān)研究成果的基礎(chǔ)上，重點(diǎn)對(duì)該方法的基本原理及初步驗(yàn)證結(jié)果進(jìn)行介紹。

1 模型方法

1.1 電離層延遲模型

傳統(tǒng)電離層延遲建模中通常僅考慮電離層延遲整體變化趨勢(shì)，為了有效地顧及電離層延遲變化的局部特點(diǎn)，本文將電離層延遲分為趨勢(shì)項(xiàng)與隨機(jī)項(xiàng)兩部分，如式(1)所示。

(1)

在電離層薄層假設(shè)下，假設(shè)某電離層格網(wǎng)點(diǎn)(ionosphericgridpoint，IGP)周邊有效范圍內(nèi)的交叉點(diǎn)(ionosphericpiercepoint，IPP)個(gè)數(shù)為N[21]，則可得到交叉點(diǎn)電離層延遲的觀測(cè)方程及格網(wǎng)點(diǎn)電離層延遲的內(nèi)插方程，如式(2)所示

(2)

式(2)中，L1表示交叉點(diǎn)實(shí)測(cè)電離層延遲組成的列向量；V1表示其觀測(cè)誤差向量，對(duì)應(yīng)的協(xié)方差矩陣為DV1；A1表示觀測(cè)方程的設(shè)計(jì)矩陣，與交叉點(diǎn)分布有關(guān)；X表示電離層延遲趨勢(shì)項(xiàng)模型待估系數(shù)組成的列向量；Y1表示交叉點(diǎn)電離層隨機(jī)項(xiàng)組成的列向量，對(duì)應(yīng)協(xié)方差矩陣為DY1，假設(shè)V1與Y1是不相關(guān)的；L2表示格網(wǎng)點(diǎn)電離層延遲估值，A2表示內(nèi)插方程的設(shè)計(jì)矩陣，與A1形式相同；Y2表示格網(wǎng)點(diǎn)處的電離層隨機(jī)項(xiàng)。

第一步：將觀測(cè)方程改寫(xiě)成如式(3)所示的形式。

L1A+V1A=A1X

(3)

式(3)中：

(4)

式(4)即為參數(shù)X的觀測(cè)方程，L1A即可看作是A1X的觀測(cè)量，V1A表示其觀測(cè)誤差，在不考慮V1與Y1的相關(guān)性時(shí)，其對(duì)應(yīng)的協(xié)方差矩陣如式(5)所示。

DV1A=DV1+DY1

(5)

因此，可直接采用經(jīng)典最小二乘法則，得到參數(shù)X的最優(yōu)解如式(6)所示。

(6)

第二步：同樣，可根據(jù)經(jīng)典最小二乘準(zhǔn)則得到未知參數(shù)Y1B，即Y1的估值，如式(7)所示。

(7)

至此，可以得到未知參量X及隨機(jī)參量Y1的最優(yōu)估值。基于各交叉點(diǎn)實(shí)測(cè)電離層隨機(jī)項(xiàng)即可估計(jì)得到格網(wǎng)點(diǎn)電離層隨機(jī)項(xiàng)，如式(8)所示。

(8)

式(8)中，

(9)

1.2 電離層延遲協(xié)方差函數(shù)

合理的協(xié)方差函數(shù)是格網(wǎng)點(diǎn)電離層延遲估計(jì)的關(guān)鍵，通常推薦采用經(jīng)驗(yàn)的固定協(xié)方差函數(shù)[24，26]。但是，由于電離層變化受到多種因素的綜合控制，統(tǒng)一固定的協(xié)方差函數(shù)難以準(zhǔn)確地描述任意處任意時(shí)刻電離層變化的相關(guān)性。為此，本文嘗試采用球狀分段函數(shù)建立適合于中國(guó)區(qū)域的實(shí)時(shí)電離層協(xié)方差函數(shù)[32]，如式(10)～式(11)所示。

(10)

(11)

式(11)中，c0+c表示拱高值；c0表示基臺(tái)值；R表示變程值；利用統(tǒng)計(jì)得到的經(jīng)驗(yàn)變異函數(shù)值，基于最小二乘可解算c0，c，R三個(gè)待估參數(shù)。

如圖1所示，給出了基于上述函數(shù)擬合得到的2002年與2008年第76天中國(guó)低緯地區(qū)電離層變異函數(shù)及協(xié)方差函數(shù)?？梢钥吹剑们驙罘侄魏瘮?shù)可實(shí)時(shí)建立中國(guó)區(qū)域電離層協(xié)方差函數(shù)，反映不同交叉點(diǎn)處電離層延遲隨機(jī)項(xiàng)之間的相關(guān)性。有關(guān)實(shí)驗(yàn)數(shù)據(jù)詳見(jiàn)第2節(jié)。

圖1 基于球狀分段函數(shù)得到的2002年(A)與2008年(B)第76天06：15(UTC)中國(guó)低緯地區(qū)電離層變異函數(shù)及協(xié)方差函數(shù)

1.3 用戶電離層延遲估計(jì)方法

用戶接收到系統(tǒng)播發(fā)的格網(wǎng)點(diǎn)電離層延遲信息及實(shí)時(shí)協(xié)方差函數(shù)，即可計(jì)算衛(wèi)星視線方向上的電離層延遲改正值。如圖2所示，假設(shè)用戶某可視衛(wèi)星交叉點(diǎn)為P，其周圍4個(gè)格網(wǎng)點(diǎn)分別是A、B、C、D，則用戶電離層延遲改正值及其精度的計(jì)算方法，如式(12)所示。

圖2 廣域增強(qiáng)系統(tǒng)用戶端電離層延遲改正信息計(jì)算示意圖

(12)

(13)

式(13)中，C(dMN)表示基于協(xié)方差函數(shù)計(jì)算得到的點(diǎn)M與N之間的協(xié)方差。

(14)

式(14)中，IA、IB、IC與ID分別表示四個(gè)格網(wǎng)點(diǎn)A、B、C與D處的電離層延遲信息。

(15)

式(15)中，σA、σB、σC與σD分別表示格網(wǎng)點(diǎn)A、B、C與D處電離層延遲精度指標(biāo)。

2 試驗(yàn)驗(yàn)證

2.1 實(shí)驗(yàn)數(shù)據(jù)概況

為了有效地驗(yàn)證本文所介紹的中國(guó)區(qū)域電離層延遲建模方法的有效性，選擇2002年第76～82天、2008年第76～82天、2012年第76 ～82天共計(jì)三周實(shí)測(cè)GPS數(shù)據(jù)對(duì)本文方法的精度與可靠性進(jìn)行分析。2002年、2008年、2012年依次屬于電離層活動(dòng)高、低、中年，所選時(shí)段處于當(dāng)年春分附近。圖3和圖4分別給出了上述三個(gè)實(shí)驗(yàn)時(shí)段所選用的基準(zhǔn)站與監(jiān)測(cè)站分布。由于2002年與2008年中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)基準(zhǔn)站數(shù)量有限，共選擇了21個(gè)基準(zhǔn)站，基準(zhǔn)站在低緯地區(qū)分布較為稀疏；2012年共選擇29個(gè)基準(zhǔn)站，所增加的基準(zhǔn)站大部分分布于我國(guó)低緯度與邊境地區(qū)；監(jiān)測(cè)站共選擇8個(gè)，其中，中緯度與低緯度地區(qū)各4個(gè)，并且在中緯度和低緯度地區(qū)分別有一個(gè)監(jiān)測(cè)站(ULAB與TWTF)位于中國(guó)大陸境外。

(▲：基準(zhǔn)站；●：監(jiān)測(cè)站)圖3 2002年與2008年實(shí)驗(yàn)所選基準(zhǔn)站與監(jiān)測(cè)站分布示意圖

(▲：基準(zhǔn)站；●：監(jiān)測(cè)站)圖4 2012年實(shí)驗(yàn)所選基準(zhǔn)站與監(jiān)測(cè)站分布示意圖

2.2 電離層延遲建模精度分析

(16)

圖5給出了實(shí)驗(yàn)期間上述各監(jiān)測(cè)站修正后殘余誤差的概率分布圖，其中，第一行與第二行的三幅圖分別表示中緯度和低緯度地區(qū)的監(jiān)測(cè)站，分布曲線對(duì)應(yīng)的1σ在圖例中各監(jiān)測(cè)站后的括弧中給出，單位為m。總體上看，在本實(shí)驗(yàn)基準(zhǔn)站布設(shè)下，我國(guó)及周邊中緯度地區(qū)，電離層格網(wǎng)精度在電離層活動(dòng)高年優(yōu)于0.4m，在電離層活動(dòng)低年優(yōu)于0.2m；我國(guó)低緯度地區(qū)，電離層格網(wǎng)精度在電離層活動(dòng)高年約為0.4～0.7m，在電離層活動(dòng)低年優(yōu)于0.2m。

圖5 2002/2008/2012年實(shí)驗(yàn)期間不同監(jiān)測(cè)站格網(wǎng)電離層修正殘余誤差統(tǒng)計(jì)分布

對(duì)比圖5中不同監(jiān)測(cè)站的統(tǒng)計(jì)結(jié)果，在中緯度地區(qū)，2002年電離層格網(wǎng)在境外ULAB監(jiān)測(cè)站的修正精度為0.34m，相對(duì)于境內(nèi)其他三個(gè)監(jiān)測(cè)站上降低約20%，這主要是由于電離層TEC模型外推誤差造成的；BJFS監(jiān)測(cè)站附近布設(shè)有基準(zhǔn)站BJSH(站間距約為80km)，從而使得BJFS站的修正精度最高(0.20m)；ZHNZ監(jiān)測(cè)站次之(0.26m)，URUM監(jiān)測(cè)站最差(0.28m)，其主要仍是由基準(zhǔn)站布設(shè)不均勻造成的，在ZHNZ監(jiān)測(cè)站西部布設(shè)有基準(zhǔn)站XIAA(站間距約為400km)，而距離URUM監(jiān)測(cè)站最近的基準(zhǔn)站為WUSH(站間距約800km)。2008年電離層格網(wǎng)在BJFS與ZHNZ監(jiān)測(cè)站上的精度基本相當(dāng)(0.10m)，而在URUM與ULAB監(jiān)測(cè)站上的精度分別為0.15m與0.18m；2012年電離層格網(wǎng)在BJFS、URUM、ZHNZ監(jiān)測(cè)站上的精度基本相當(dāng)，約為0.17m，在URUM監(jiān)測(cè)站上精度(0.16m)的提高主要得益于附近增加布設(shè)的3個(gè)基準(zhǔn)站(XJAL、XJWQ、XJTZ)，從而使得監(jiān)測(cè)站與基準(zhǔn)站的站間距縮短至450km左右。

在低緯度地區(qū)，2002年電離層格網(wǎng)在WUHN監(jiān)測(cè)站附近修正精度最高(0.39m)，XIAG監(jiān)測(cè)站附近次之(0.56m)，TWTF監(jiān)測(cè)站附近最低(0.68m)；WUHN監(jiān)測(cè)站附近精度相對(duì)較高主要是由于其附近布設(shè)WHJF基準(zhǔn)站(站間距約為20km)；XIAG監(jiān)測(cè)站附近布設(shè)有基準(zhǔn)站KMIN(站間距約為280km)；TWTF監(jiān)測(cè)站附近電離層格網(wǎng)主要依靠XIAM基準(zhǔn)站(站間距約350km)上電離層觀測(cè)值進(jìn)行外推，再加上低緯地區(qū)電離層TEC變化梯度較大，導(dǎo)致其精度出現(xiàn)明顯衰減；比較而言，在電離層活動(dòng)低年，電離層格網(wǎng)在TWTF監(jiān)測(cè)站附近的精度未出現(xiàn)明顯衰減。

因此，為了保證格網(wǎng)電離層具有一定的修正精度，對(duì)于我國(guó)中緯度地區(qū)，基準(zhǔn)站間距可適當(dāng)放大至800km左右，而在我國(guó)低緯度地區(qū)基準(zhǔn)站間距應(yīng)保持在350～450km之間。

3 結(jié)束語(yǔ)

針對(duì)我國(guó)區(qū)域電離層TEC變化特點(diǎn)，基于擬合推估理論初步建立了我國(guó)廣域增強(qiáng)系統(tǒng)格網(wǎng)電離層建模的方法。該方法在利用球諧函數(shù)優(yōu)良的數(shù)學(xué)特性描述電離層變化趨勢(shì)項(xiàng)的基礎(chǔ)上，通過(guò)實(shí)時(shí)建立反映中國(guó)區(qū)域電離層相關(guān)性的協(xié)方差函數(shù)，精確估計(jì)了格網(wǎng)點(diǎn)處電離層延遲及其精度信息，有效顧及了局部地區(qū)電離層變化特點(diǎn)；在顯著提高局部地區(qū)格網(wǎng)點(diǎn)電離層延遲信息估計(jì)精度的同時(shí)，保證了區(qū)域電離層延遲信息格網(wǎng)的整體精度與可靠性。

基于不同電離層活動(dòng)水平下(2002年、2008年與2012年)的實(shí)測(cè)數(shù)據(jù)，分別對(duì)本文所提出格網(wǎng)電離層建模方法在中國(guó)高、中、低緯度地區(qū)的修正精度進(jìn)行了分析。初步結(jié)果顯示，基于擬合推估的電離層延遲建模方法的修正精度在電離層活動(dòng)高年約為0.4～0.7m，在電離層活動(dòng)低年優(yōu)于0.2m。本文所提出的電離層建模方法將為我國(guó)北斗廣域增強(qiáng)系統(tǒng)建立中電離層延遲修正提供重要參考。

[1]LIXing-xing，GEMao-rong，ZHANGHong-ping，etal.TheGFZReal-timeGNSSPrecisePositioningServiceSystemandItsAdaptionforCOMPASS[J].AdvancesinSpaceResearch，2012，51(6)：1008-1018.

[2]JUANJM，SANZJ，HERNáNDEZ-PAJARESM，etal.WideAreaRTK：ASatelliteNavigationSystemBasedonPreciseReal-timeIonosphericModelling[J].RadioScience，2012，47(RS2016)：1-14.

[3]WALTERT，ENGEP，REDDANP.ModernizingWAAS[C]//Proceedingsofthe17thInternationalTechnicalMeetingoftheSatelliteDivisionofTheInstituteofNavigation(IONGNSS2004).LongBeach，CA：InstituteofNavigation，2004：1683-1690.

[4]AGARWALS.WAASAugmentationforCatIIandCatIIIOperations[R].Bombay：DepartmentofAerospaceEngineering，IndianInstituteofTechnology，2009.

[5] 李子申.GNSS/Compass電離層時(shí)延修正及TEC監(jiān)測(cè)理論與方法研究[D].武漢：中國(guó)科學(xué)院測(cè)量與地球物理研究所，2012.

[6]HERNNDEZ-PAJARESM，JUANJ，SANZJ，etal.TheIonosphere：Effects，GPSModelingandtheBenefitsforSpaceGeodeticTechniques[J].JournalofGeodesy，2011，85(12)：887-907.

[7]WARNANTR，F(xiàn)OELSCHEU，AQUINOM，etal.MitigationofIonosphericEffectsonGNSS[J].AnnalsofGeophysics，2010，52(3-4)：373-390.

[8]DINGWen-wu,OUJi-kun,LIZi-shen,etal.AnInstantaneousRe-initializationMethodofRealTimeKinematicbyAddingIonosphericDelayConstraints[J].ChineseJournalofGeophysics，2014，57(3)：280-292.

[9]ALTSHULERE，CORMIERD，GOH.ImprovementtotheWAASIonosphericAlgorithms[C]//Proceedingsofthe17thInternationalTechnicalMeetingoftheSatelliteDivisionofTheInstituteofNavigation(IONGNSS2002).Portland：InstituteofNavigation，2002：2256-2261.

[10]袁運(yùn)斌.基于GPS的電離層監(jiān)測(cè)及延遲改正理論與方法的研究[D].武漢：中國(guó)科學(xué)院測(cè)量與地球物理研究所，2002.

[11]HERNNDEZ-PAJARESM，ZORNOZAJ，SUBIRANAJ，etal.EGNOSTestBedIonosphericCorrectionsundertheOctoberandNovember2003Storms[J].IEEETransactionsonGeoscienceandRemoteSensing，2005，43(10)：2283-2293.

[12]LIZi-shen，YUANYun-bin，WANGNing-bo，etal.SHPTS：towardsaNewMethodforGeneratingPreciseGlobalIonosphericTECMapBasedonSphericalHarmonicandGeneralizedTrigonometricSeriesFunctions[J].JournalofGeodesy，2015，89(4)：331-345.

[13]YUANYun-bin，OUJi-kun.ModelsandMethodsforPreciseDeterminationofIonosphericDelaysUsingGPS[J].ProgressInNaturalScience，2007，17(2)：187-196.

[14]RATNAMDV，SARMAAD.ModelingofIndianIonosphereusingMMSEEstimatorforGAGANApplications[J].TheJournalofIndianGeophysicalUnion，2006，10(4)：303-312.

[15]YUANYun-bin，OUJi-kun.AGeneralizedTrigonometricSeriesFunctionModelforDeterminingIonosphericDelay[J].ProgressInNaturalScience，2004，14(11)：1010-1014.

[16]霍星亮.基于GNSS的電離層形態(tài)監(jiān)測(cè)與延遲模型研究[D].武漢：中國(guó)科學(xué)院測(cè)量與地球物理研究所，2008.

[17]SHUKLAAK，DASS，NAGORIN，etal.Two-shellIonosphericModelforIndianRegion：ANovelApproach[J].IEEETransactionsonGeoscienceandRemoteSensing，2009，47(8)：2407-2412.

[18]RATNAMDV，SARMAAD.ModelingofIndianIonosphereusingMMSEEstimatorforGAGANApplications[J].TheJournalofIndianGeophysicalUnion，2006，10(4)：303-312.

[19]ACHARYAR，NAGORIN，JAINN，etal.IonosphericStudiesfortheImplementationofGAGAN[J].IndianJournalofRadioandSpacePhysics，2007，36(5)：394-404.

[20]HERNNDEZ-PAJARESM，ZOMOZAJMJ，SUBIRANAJS，etal.FeasibilityofWide-areaSubdecimeterNavigationwithGALILEOandModernizedGPS[J].IEEETransactionsonGeoscienceandRemoteSensing，2003，41(9)：2128-2131.

[21]YUANYun-bin，OUJi-kun.DifferentialAreasforDifferentialStations(DADS)：ANewMethodofEstablishingGridIonosphericModel[J].ChineseScienceBulletin，2002，47(12)：1033-1036.

[22]王一舉，袁洪，萬(wàn)衛(wèi)星.觀測(cè)站稀疏地區(qū)的WAAS電離層時(shí)延網(wǎng)格修正算法[J].空間科學(xué)學(xué)報(bào).2003，33(2)：1-5.

[23]王一舉，周建華.改進(jìn)的WAAS電離層時(shí)延網(wǎng)格修正算法[J].空間科學(xué)學(xué)報(bào).2002，22(2)：129-135.

[24]BLANCHJ，WALTERT，ENGEP.IonosphericEstimationUsingExtendedKrigingforaLowLatitudeSBAS[C]//Proceedingsofthe17thInternationalTechnicalMeetingoftheSatelliteDivisionofTheInstituteofNavigation(IONGNSS2004).SanDiego，CA：TheInstituteofNavigation，2004：21-24.

[25]SPARKSL，BLANCHJ，PANDYAN.EstimatingIonosphericDelayUsingKriging：2.ImpactonSatellite‐BasedAugmentationSystemAvailability[J].RadioScience，2011，46(RS0D22)：1-10.

[26]SPARKSL，BLANCHJ，PANDYAN.EstimatingIonosphericDelayUsingKriging：1.Methodology[J].RadioScience，2011，46(RS0D21)：1-13.

[27]ORSR，HERNNDEZ-PAJARESM，JUANJ，etal.ImprovementofGlobalionosphericVTECMapsbyUsingKrigingInterpolationTechnique[J].JournalofAtmosphericandSolar-TerrestrialPhysics，2005，67(16)：1598-1609.

[28]LIZi-shen，YUANHong，YUANYun-binetal.ModelingtheIonosphericDelayBasedonCollocationforBDSWideAreaAugmentationSystem[C]//Proceedingsofthe27thInternationalTechnicalMeetingofTheSatelliteDivisionoftheInstituteofNavigation(IONGNSS+ 2014).Tampa，F(xiàn)lorida：TheInstituteofNavigation，2014：2457-2466.

[29]周江文.再論擬合推估[J].測(cè)繪學(xué)報(bào).2001，30(4)：283-285.

[30]周江文.擬合推估新解之一——兩步解法[J].測(cè)繪學(xué)報(bào).2002，31(3)：189-191.

[31]楊元喜，劉念.擬合推估兩步極小解法[J].測(cè)繪學(xué)報(bào).2002，31(3)：192-195.

[32]OLEAR.GeostatisticsforNaturalResourcesEvaluation[M].NewYork：OxfordUniversityPress，1999.

[33]張寶成.GNSS非差非組合精密單點(diǎn)定位的理論方法與應(yīng)用研究[D].武漢：中國(guó)科學(xué)院測(cè)量與地球物理研究所，2012.

Precise Modeling of Ionospheric Delay over China Region Based on Collocation

YUANYun-bin1，3，LIZi-shen1，2，WANGNing-bo1，3，4，ZHANGBao-cheng1，3，5，HUOXing-liang1，3

(1.Institute of Geodesy and Geophysics，Chinese Academy of Sciences，Wuhan 430077，China； 2.Academy of Opto-Electronics，Chinese Academy of Sciences，Beijing 100094，China； 3.State Key Laboratory of Geodesy and Earth Dynamics，Wuhan 430077，China； 4.University of Chinese Academy of Sciences，Beijing 100049，China； 5.Curtin University，Perth 6845，Australia)

Grid-based ionospheric information with high accuracy is crucial to the Wide Area Augmentation System (WAAS) establishment and operation.With the specialty and complexity of the ionospheric conditions over China region，an approach is proposed to model the Chinese regional ionosphere using collocation after splitting the ionospheric delays into a deterministic trend together with a stochastic component.The usefulness of the method is validated with GPS measurements collected during high-，medium- and low-solar activities over China.It follows from our analyses that，the accuracy of our method resides within 0.4-0.7 m during high-solar activity year，whereas it is reduced to approximately 0.2 m during low-solar activity year.Considering the mid-latitude area of China，the distance between reference stations can maintain about 800 km.For the low-latitude area，however，such a distance should be limited between 350 km and 450 km.

ionospheric delay，GNSS，collocation，BDS wide area augmentation system

袁運(yùn)斌，李子申，王寧波，等.基于擬合推估的中國(guó)區(qū)域電離層延遲精確建模方法[J].導(dǎo)航定位學(xué)報(bào),2015,3(3):49-55.(YUAN Yun-bin,LI Zi-shen,WANG Ning-bo,et al.Precise Modeling of Ionospheric Delay over China Region Based on Collocation[J].Journal of Navigation and Positioning,2015,3(3):49-55.)

10.16547/j.cnki.10-1096.20150310.

2015-05-18

973項(xiàng)目(2012CB25604)，國(guó)家自然科學(xué)基金(41231064、41304034、41204031、41321063)，中國(guó)科學(xué)院交叉創(chuàng)新團(tuán)隊(duì)(KZZD-EW-TZ-05)，北京市自然科學(xué)基金(4144094)，國(guó)際GNSS監(jiān)測(cè)評(píng)估系統(tǒng)(IGMAS)分析中心建設(shè)與運(yùn)行維護(hù)項(xiàng)目，大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金(SKLGED2014-3-1-E)。

袁運(yùn)斌(1971—)，男，江西余干人，研究員，主要從事衛(wèi)星導(dǎo)航電離層方面的研究。

P228

A

2095-4999(2015)-03-0049-07