空間緯度對于差分系統(tǒng)電離層格網(wǎng)性能的影響＊

馮 煒 陳向東 吳 星 葉修松

(1)北京環(huán)球信息應(yīng)用開發(fā)中心,北京 100094 2)北京特種工程設(shè)計研究院,北京 100028 3)海軍海洋測繪研究所,天津 300061 4)西安衛(wèi)星測控中心,西安 710043)

空間緯度對于差分系統(tǒng)電離層格網(wǎng)性能的影響＊

馮 煒1)陳向東1)吳 星2,3)葉修松4)

(1)北京環(huán)球信息應(yīng)用開發(fā)中心,北京 100094 2)北京特種工程設(shè)計研究院,北京 100028 3)海軍海洋測繪研究所,天津 300061 4)西安衛(wèi)星測控中心,西安 710043)

基于廣域精密實時定位系統(tǒng)試算的電離層格網(wǎng)改正數(shù)據(jù),使用兩種加權(quán)方法對參考站上空信號穿刺點的電離層延遲進行了估算,并利用估算出的穿刺點電離層格網(wǎng)計算延遲和實際測量的穿刺點電離層延遲進行了比較。結(jié)果表明對于不同位置的測站,中高緯度格網(wǎng)改正較好,低緯度較差,這與電離層的特性有關(guān),并不是加權(quán)內(nèi)插校正方案造成的。

差分系統(tǒng);廣域精密實時定位系統(tǒng);電離層;格網(wǎng);加權(quán)內(nèi)插

1 引言

確定電磁波在穿越電離層時產(chǎn)生的電離層延遲是提高衛(wèi)星導(dǎo)航系統(tǒng)精度的重要問題。隨著衛(wèi)星導(dǎo)航技術(shù)的不斷發(fā)展和人們對精度要求的不斷提高,很多學(xué)者相繼提出了各種電離層改正模型和方法。為了使眾多的單頻用戶獲得精度更高、更可靠的導(dǎo)航服務(wù),廣域差分系統(tǒng) (WAAS)應(yīng)運而生。它通過差分手段向服務(wù)區(qū)內(nèi)的用戶提供各類差分校正信息以及完好性狀態(tài)[1,2]。

差分系統(tǒng)用戶的電離層延遲校正問題是目前廣域差分研究中的熱點問題之一,本文將基于廣域?qū)崟r精密定位與示范系統(tǒng)試算的電離層格網(wǎng)數(shù)據(jù),分析電離層格網(wǎng)校正算法在不同的空間緯度上的性能差別。具體方法是：結(jié)合廣域?qū)崟r精密定位與示范系統(tǒng)電離層估計軟件試算出的格網(wǎng)數(shù)據(jù),同時使用Junkins加權(quán)法與雙線性模型兩種加權(quán)方法對不同緯度的兩個基準(zhǔn)站進行格網(wǎng)校正。利用雙頻實測電離層延遲對比格網(wǎng)內(nèi)插電離層延遲的策略,對格網(wǎng)校正問題中緯度的影響進行詳細的分析。

2 穿刺點實測電離層延遲的獲取

雙頻能夠有效且精確地實時解算穿刺點的電離層延遲,本文利用雙頻實測電離層延遲對比格網(wǎng)內(nèi)插電離層延遲,需要實時解算參考站穿刺點上的電離層延遲。

2.1 雙頻平滑偽距算法[1]

由雙頻觀測值有[3]：

由式 (1)、(2)聯(lián)立可得：

其中：

2.2 電離層投影函數(shù)

TEC是方向的函數(shù),其數(shù)值將隨信號傳播路徑的高度角及方位角的變化而變化。其中方位角影響較小,可以忽略其帶來的變化。而在計算穿刺點TEC的變化時,采用天頂方向的總電子含量來討論,并記為VTEC。VTEC是高度角為 90°時的 TEC,某一時刻某一地點的VTEC值是唯一的,它是 TEC中的最小值。

因此需要引入電離層單層模型 SLM的概念,根據(jù) Klobuchar模型[4],單層距地面的高度取 H=450 km。

三角函數(shù)型 SLM投影函數(shù)為：

z′為穿刺點處的天頂距,使用改進投影函數(shù),按照高度角的不同用一個分段權(quán)函數(shù)與原函數(shù)相乘[5]。

P為分段權(quán),E為衛(wèi)星高度角。計算穿刺點處的天頂距 z′的公式如下：

R為測站處地球半徑,R0為地球平均半徑,H為單層距地面的高度,z為接收機處的衛(wèi)星天頂距。

由以上可得：

使用電離層單層模型 SLM可實現(xiàn)傳播路徑上的電離層延遲到單層模型垂直方向上延遲之間的轉(zhuǎn)換,改進 SLM投影函數(shù)可以更為精確地描述投影特征,這樣可以方便地使用實際測量數(shù)據(jù)確定穿刺點上的電離層延遲。

3 差分系統(tǒng)的格網(wǎng)校正

電離層延遲隨時間、地點的改變而變化。為各個較小的區(qū)域分別提供實時電離層改正是提高精度的有效途徑,電離層格網(wǎng)改正技術(shù)正是這樣一種分別對各個小區(qū)域提供近乎實時電離層校正的方法。廣域差分系統(tǒng)一般使用此種方法向單頻用戶廣播電離層格網(wǎng)數(shù)據(jù)[1]。

電離層格網(wǎng)是基于一種人為規(guī)定的球面網(wǎng)格。該球面的中心與地心重合,在假想球面上也定義了相應(yīng)的經(jīng)線和緯線,電離層網(wǎng)格點就分布在該假想球面上。網(wǎng)格點的間隔為 5°,如果廣域差分系統(tǒng)能實時地提供各網(wǎng)格點的垂直電離層延遲改正值,用戶就可以利用網(wǎng)格內(nèi)插法獲得非常精確的電離層延遲改正[7]。

3.1 電離層格網(wǎng)校正的兩種加權(quán)算法

用戶接收到廣域增強系統(tǒng)廣播的網(wǎng)格點電離層延遲后,可采用內(nèi)插法計算用戶穿透點的電離層延遲,即利用穿透點所在網(wǎng)格頂點的校正數(shù)據(jù)進行加權(quán)計算。

用戶穿透點的垂直電離層延遲 IU計算公式為：

其中,K為用于內(nèi)插的格點個數(shù),一般為 4(圖 1)。但當(dāng) 4個格點中的某一個不可用時,如剩余的 3個點包圍了用戶穿透點,則用這 3個點計算。

當(dāng) K為 4,權(quán)可取為：

式 (11)x=Δ λ/(λ2-λ1),y=Δφ/(φ2-φ1),W權(quán)的計算可用

及 Junkins加權(quán)模型

計算。式中 Junkins加權(quán)模型是 FAA的推薦模型[1,2]。

圖1 4點內(nèi)插示意圖Fig.1 Sketch of four-point interpolation

4 結(jié)果分析

鑒于中國地殼運動 GPS監(jiān)測網(wǎng)的公開性較好、分布也較為合理,是驗證程序使用比較廣泛的數(shù)據(jù)源,所以本文采用的 GPS數(shù)據(jù)流是利用中國地殼運動 GPS監(jiān)測網(wǎng) 2008年 1月8日的 22個參考站數(shù)據(jù)模擬出的實時數(shù)據(jù)流 (采樣率為 30 s)。格網(wǎng)數(shù)據(jù)是利用上述測站數(shù)據(jù)通過廣域?qū)崟r定位系統(tǒng)電離層估計軟件計算出的格網(wǎng)改正數(shù)據(jù)。所用的站點如下：BJFS、GUAO、SHAO、WUSH、BJSN、HLAR、SU IY、XI AG、CHUN、J I XN、TA I N、XI AM、DLHA、LHAS、TASH、XN I N、DXI N、LUZH、T WTF、GUAN、Q I ON、WHJF。以這 22個基準(zhǔn)站的模擬實時數(shù)據(jù)計算了我國上空的 5°×5°電離層參考面上的每個網(wǎng)格結(jié)點的垂直電離層延遲。圖 2反應(yīng)了BJFS基準(zhǔn)站上空 20分鐘時間內(nèi) 25號衛(wèi)星的穿刺點情況。

通過圖 2可以看出該衛(wèi)星在這段時間內(nèi)的穿刺點是在 (30°,120°)、(30°,125°)、(35°,120°)(35°, 125°)4個格網(wǎng)點內(nèi)。圖 3顯示了該段時間內(nèi)廣域?qū)崟r精密定位系統(tǒng)計算的該 4個格網(wǎng)點上的垂直電離層延遲值變化情況。

圖2 站點坐標(biāo)示意圖Fig.2 Sketch of site positions

圖3 4個格網(wǎng)點上的電離層垂直延遲Fig.3 Ionospheric vertical delays at the four grid points

利用穿刺點周圍 4格網(wǎng)點數(shù)據(jù)使用兩種加權(quán)方法計算穿刺點電離層垂直延遲并與利用第二節(jié)方法計算出的穿刺點雙頻實測電離層垂直延遲做差,可反應(yīng)出兩種算法的內(nèi)插校正精度。圖 4為兩種算法得出的 2008年 1月 8日 2：00到 2008年 1月 8日2：45的垂直延遲與雙頻實測值的差值圖。

從圖 4可以清晰地反應(yīng)出兩種加權(quán)方法的差異,由于單位 TECU對測距的影響很小,所以這兩種方法在性能上區(qū)別并不是很大,但是 FAA推薦模型的穩(wěn)定性和相關(guān)性似乎更好。這個結(jié)論可以更好的幫助我們觀察緯度對于格網(wǎng)有怎樣的影響。

圖5是與上圖相同數(shù)據(jù)時間內(nèi),兩種加權(quán)方法分別在 BJFS和 YONG參考站與實測值的差值情況。

BJFS參考站位于北京房山是中緯度地區(qū), YONG參考站位于我國西沙群島永興島,從圖 5和圖 4可以看出兩種加權(quán)的方法本身區(qū)別不大,但是中緯度地區(qū)使用格網(wǎng)計算的電離層精度明顯高,基本可以到達0.5TECU的量級,低緯度的用戶站電離層校正精度相對不高只能到 1.3 TECU左右;格網(wǎng)改正的性能和空間緯度密切相關(guān),但和格網(wǎng)內(nèi)插校正方法無關(guān)。

圖4 BJFS兩種加權(quán)算法與實測電離層垂直延遲差值Fig.4 Differences between ionospheric vertical delays calculated with the two weighted algorithm and measured atBJFS

圖5 BJFS和 YONG參考站計算值與實測值的差值Fig.5 Differences between ionospheric vertical delays calculated with the two weighted algorithm and measured atBJFS and YONG

5 結(jié)論

實驗結(jié)果表明,電離層格網(wǎng)算法是單頻用戶有效地快速獲得電離層改正參數(shù)的途徑。FAA推薦的加權(quán)內(nèi)插校正方法相關(guān)性更好,并且很穩(wěn)定。它和雙線性加權(quán)內(nèi)插校正方法給單頻用戶帶來的精度提升基本相當(dāng)。同時發(fā)現(xiàn)格網(wǎng)內(nèi)插校正的精度變化和加權(quán)內(nèi)插校正方法無關(guān),但和緯度變化相關(guān)。對于電離層水平梯度不大并且電子濃度較低的中高緯度地區(qū),格網(wǎng)算法本身的改正效果很好。但是對于電子變化更加復(fù)雜、電子濃度較高的低緯度地區(qū)改正精度相對比較差。

致謝 感謝武漢大學(xué)國家 GPS工程中心的支持。

1 袁運斌.基于 GPS的電離層監(jiān)測及延遲改正理論與方法的研究[D].中國科學(xué)院研究生院,2002.(Yuan Yunbin. Study on theories andmethodsof correcting ionospheric delay and monitoring ionosphere based on GPS[D].Graduate U-niversity of Chinese Academy of Sciences,2002)

2 U S Depart ment of Transportation and Federal Aviation Administration.Specification for W ide Are Augmentation System(WAAS)[S].USA：FAA,1999.

3 Steven C First and kamran Ghssemt.GPS the next generation[J].Proceeding of the IEEE,1999,86(3)：59.

4 LogiVitharsson,et al.Satellite autonomous integritymonitoring and its role in enhancing GPS user performance[S]. USA：Standford University and NavAstro CO.2004.

5 楊力.大氣對 GPS測量影響的理論與研究[D].解放軍信息工程大學(xué),2001.(YangLi.The theory and research of atmosphere affection to GPS surveying[D].Information EngineeringUniversity,2001)

6 袁運斌,歐吉坤.GPS觀測數(shù)據(jù)中的儀器偏差對確定電離層延遲的影響及處理方法 [J].測繪學(xué)報,1999,28 (2)：110-114.(Yuan Yunbin and Ou Jikun.The effects of instrumental bias in GPS observations on deter mining ionospheric delays and the methods of its calibration[J].Acta Geodaetica et Cartographica sinica,1999,28(2)：110-114)

7 袁運斌.利用 IGS的 GPS資料確定全球電離層 TEC的初步結(jié)果與分析[J].自然科學(xué)進展,2003,(8)：885-888. (Yuan Yunbin.The preliminary result and analysis of global ionosphere TEC using GPS data of IGS[J].Progress inNatural Science,2003,(8)：885-888)

INFLUENCE OF SPACE LATITUDE ON D IFFERENTIAL I ONOSPHERIC GRI D SYSTEM

FengWei1),Chen Xiangdong1),Wu Xing2,3)and Ye Xiusong4)

(1)B eijing Application and Developm ent Center of Round-the-W orld Infor m ation,B eijing 100094 2)B eijing Special Type Engineering Design and Reasearch Institute,B eijing 100028 3)N aval Institute of Hydrographic Surveying and Charting,Tianjin 300061 4)X i’an Satellite Control Center,X i’an 710043)

On the basis of the ionospheric grid-correction data which are tested by wide-area precise real-time positioning system.The ionospheric delay of the signal puncture point over the reference station is estimated by using two weighting methods.Futhermore the esti mated ionospheric grid delay of the signal puncture point has been compared with those tested in the actualmeasurement.For the stations located in different positions,in the middle or high-latitude the grid correction is better than that in low-latitude.This result reveals that it is concerning to the characteristics of the ionosphere,but caused byweighted interpolation correction scheme.

differential system;wide-area real-time precision positioning system;ionospheric;grid;weighted interpolation

1671-5942(2011)04-0135-04

2010-12-29

國家自然科學(xué)基金(40774031);全國優(yōu)秀博士學(xué)位論文作者專項基金(200344);國家 863計劃項目 (2007AA120603)

馮煒,男,1982年生,碩士,主要從事測量與導(dǎo)航工程研究.E-mail：pladsps@sina.com

P228.41

A