探究不同緯度對北斗電離層延遲的影響

，

(1.沈陽航空航天大學 民用航空學院，沈陽 110136；2.沈陽航空航天大學 電子信息工程學院，沈陽 110136)

0 引言

電離層是地球大氣的重要電離區域，容易使穿過其中的無線電導航信號發生信號衰減、信號彎曲、信號折射等色散效應，這種現象在衛星導航中容易引起電離層延遲，從而會給導航定位精度帶來誤差。為了減小電離層延遲引起的誤差，提高導航定位精度，需要對電離層延遲改正誤差進行建模。目前國內外專家學者提出了很多電離層改正的方法，但均存在一定的不確定性[1-2]。因此，針對星基增強系統(SBAS，Satellite-Based Augmentation System)中的電離層延遲誤差，有兩種解決方法：針對單頻用戶采用Klobuchar模型計算電離層延遲，而雙頻用戶則采用格網電離層模型進行誤差校正。采用上述兩種方法，用戶能夠獲得更為精確的電離層誤差修正。

Klobuchar模型是用于廣域差分(WADGPS)中用戶測站進行電離層誤差校正的理想方法[3]，它能夠直觀簡潔地反映出電離層的周期變化特性[4]。并且為了清晰的體現電離層的變化特性，文中加入隨時間變化的三角余弦函數，這樣可以更準確的計算出不同時間段的電離層延遲誤差。太陽活動對電離層的影響是全球性的[5]，是引起電離層變化的主要因素，其中相對太陽黑子數是影響太陽活動強度的因素之一。因此，考慮這一因素，對Klobuchar模型的周期引入相對太陽黑子數這一變量。

格網電離層模型把電離層描述為地球表面上空約350km處的殼層,在此殼層上建立矩形或其他形式的格網[6]。隨著空間需求的不斷提高、大規模電離層格網的建立，更加需要優化格網電離層模型來提高導航定位精度，這對我國的太空發展計劃有重要的戰略意義。格網電離層模型中，參考面需根據一定的經度和緯度間隔網格化，網格結點播發SBAS中的電離層延遲誤差。主控站按照各個參考站實時監測的穿透點的垂直電離層延遲值，通過內插方法獲得準確的網格結點的垂直電離層延遲誤差。其原理與基于參考站電離層改正法根據各參考站數據測得流動站的誤差延遲項相似。在基于參考站電離層改正方法的幾種內插模型中，KRG內插法精度最高[7]。因此，采用反距離加權法和KRG內插法對兩個不同緯度的參考站進行格網電離層延遲校正，并與雙頻實測數據做對比，分析影響格網電離層模型修正延遲性能的因素。

1 Klobuchar模型的改進

導致模型誤差主要體現在，固定了夜間垂直延遲常數、余弦函數的振幅、周期以及相位值[8]。電離層的活動強度與太陽活動密切相關，固定夜間延遲值及相位會分別引入20%-30%和3%-10%的模型誤差[8]。因此，本文對Klobuchar模型作出如下改進:1)僅憑電文參數α算出白天三角余弦曲線的振幅是存在誤差的，因此需加入振幅改正系數對振幅改正；2)由于太陽活動的平均特性會導致電離層發生許多變化，二者密切相關。因此，需考慮相對太陽黑子數這一參數對模型周期的影響；3)夜間時延值為固定常數5 ns，這將會增加一定的模型誤差[9]。據相關資料得知，夜間VETC值可近似看作隨時間的變化的三角余弦函數。為此，需加入余弦函數來改進夜間時延值，使其具有隨時間變化的特性，即：

振幅：A2→A2(1+K1)

優化模型的數學表達式為：

(1)

式中，t為衛星與接收機傳輸的信號過程中穿過電離層處的地方時。A2為白天電離層延遲余弦函數的振幅，由an計算得到；T’為余弦曲線的周期，用bn系數求得;GPS衛星廣播星歷中導航電文發播天頂電離層改正的參數的an和bn值[10]；N為相對太陽黑子數，數據由太陽預報活動中心提供；K1、K2為振幅改正系數:

(2)

n的值通過顯著性檢驗自動選取[11]，取n=2。迭代次數取10。通過非線性迭代方法求解K1、K2收斂值,了解認識迭代函數對收斂性的影響，知道當迭代函數滿足什么條件時，振幅改正系數收斂。為了保障模型各參數本身的特性,引入監測網采集觀測量,可以實現連續解算,不斷的更新電離層延遲參數，從而使改進Klobuchar模型電離層延遲值接近實測電離層延遲值，即：

IZ-I=0

(3)

I為監測站實測的電離層延遲值。

(4)

化為等價方程

(5)

非線性迭代方法求解K1、K2收斂值流程如圖1所示。

圖1 求解K1、K2收斂值流程圖

得到各參數值后，應根據傾斜因子轉換成穿透點垂直入射的電離層延遲誤差:

I(t)=FIZ

(6)

傾斜因子反映的是GPS信號分別垂直入射和傾斜入射到同一個電離層穿透點所引起的電離層延遲的比率。計算公式如下 :

(7)

其中：Re為地球半徑；H為電離層參考面的平均高度；E為GPS衛星相對監測站的仰角。

2 格網電離延遲模型

2.1 格網電離層延遲修正過程

1)將電離層參考面根據一定的經緯度間隔分割成一定數量的網格，網格的結點作為SBAS系統播發電離層延遲誤差的依附點。

2) 參考站根據雙頻接收機接收到的信息實時計算出其穿透點電離層延遲值和北斗衛星在參考面上穿透點的經緯度。

3)參考站的數據無線通信網絡傳送至主控站，計算出電離層參考面上每個網格結點的垂直電離層延遲值，并調制成導航數據經地球同步通訊衛星發送至區域內用戶。

4)用戶通過接收到的數據計算出穿透點的經緯度，然后根據網格結點電離層延遲數據借助相應算法得到電離層延遲值。

2.2 格網點垂直延遲值的計算

電離層格網點的延遲計算值對電離層延遲的結果準確度有一定的影響，因此如何由穿透點處的垂直電離層延遲計算格網點處垂直電離層延遲是 保證格網電離層模型精度的關鍵[12]。

由于監測站采用雙頻P碼接收機，電離層的延遲數據可通過比較f1、f2的P碼觀測數據實時獲得：

(8)

其中：I為電離層穿透點的傾斜延遲值，ρ1、ρ2分別為f1、f2頻率下的P碼偽距。

電離層穿透點的傾斜延遲根據傾斜因子F轉換為電離層穿透點的垂直延遲值：

(9)

2.2.1 反距離加權法

在SBAS中，主控站通常利用反距離加權法計算出參考面內網格結點的垂直電離層延遲值。幾何示意圖參見圖2。

圖2 格網點幾何示意圖

反距離加權算法的原理是：參考面涵蓋范圍內的每個網格結點，各個監測站監測的穿透點數據按照距離該網格結點的遠近賦予不同的加權值，即穿透點距離網格結點越近，加權值越大，其對應的性質就會越相似；距離越遠加權值越小，則它們的相似性越小；距離超出一定范圍，加權值為零[13]。具體數學表達式如下：

(10)

式中，R為所選穿透點的距離值，Hj、Lj為格網點的經緯度,Hi、Li為穿透點的經緯度[14]。

2.2.2 KRG內插模型

KGR內插法又稱空間自協方差最佳插值法。它考慮在空間分布的屬性，對某一空間位置有影響的范圍內的條件值進行內插[7]。數學模型為：

(11)

(12)

Cij=[C1jC2j…Cnj]T

(13)

Cv為兩穿透點間的相關函數矩陣，相關函數Cmn=lmax-lmn,lmax取穿透點間距離最長的基線長度。

en=[1 1 … 1]T,λ為拉格朗日乘數，用于限制條件數。表達式為：

(14)

各參數確定后，即可求得系數矩陣為：

(15)

因此，格網點延遲值為：

Dj=aDi

(16)

2.3 用戶穿透點電離層延遲的計算

雖然用戶能夠接收到GEO轉發的SBAS區域內網格結點的垂直電離層延遲數據，但對于用戶來說其電離層穿透點不可能時刻都在網格結點上，所以需利用廣播的電離層延遲數據對其各個穿透點進行內插。內插后的電離層穿透點的垂直延遲為：

(17)

其中：ΦPP、λPP分別為電離層穿透點的緯度和經度，τvi為格網點的電離層垂直延遲。

根據加權函數模型f(x,y)=xy得到具體的加權函數為：

W1=xPPyPP

W2=(1-xPP)yPP

W3=(1-xPP)(1-yPP)

W4=xPP(1-yPP)

(18)

最后，穿透點垂直延遲值乘以傾斜因子即為電離層校正延遲:

(19)

2.4 格網點電離層垂直改正誤差GIVE的計算

為了驗證網絡插值算法的可行性，需計算格網點電離層垂直改正誤差GIVE。GIVE對特定時段內的的每個電離層穿透點的延遲修正殘差統計，計算格網點四周穿透點誤差限值。計算過程[15]為：

2)采用雙頻法得到穿透點電離層延遲校正值DIPP(t)，將該穿刺點的電離層延遲改正計算值與估計值進行比較，即：

(20)

3)統計該時段內的穿透點eIPP(t)的誤差限值：

(21)

其中:κ(Pr)為置信水平是99.9%的置信度，取值為3.2905。

(22)

(23)

4)網點相鄰的的網格中穿刺點的GIVE值為：

(24)

(25)

3 數據處理與分析

選取武漢、桂林、勉縣、太原、鹽城、廈門各參考站的觀測數據，采集時間為2016年3月20日12:00-15:00，采樣間隔為1 s。分別用Klobuchar模型、改進的Klobuchar模型和雙頻校正模型計算電離層延遲改正誤差。以武漢站和桂林站為例，仿真結果如圖3和圖4所示。

圖3 武漢站3種模型電離層延遲對比

圖4 桂林站3種模型電離層延遲對比

可以看出，改進的Klobuchar模型修正電離層延遲誤差可達3～5 m，高于原模型修正精度1～2 m，具有更好的修正效果，且更加接近于雙頻實測值。

分別采用反距離加權法和KRG內插法計算格網點電離層延遲值，進而得到用戶電離層延遲誤差的精度，比較兩種內插模型的效果。

圖5和圖6分別表示太原站和桂林站利用2種內插法的格網模型進行定位時的改正效果。表1和表2分別列出了兩站當日取樣時段內修正電離層延遲誤差的精度情況。可以看出，在桂林站使用兩種格網內插模型校正的電離層延遲值明顯比太原站更接近于雙頻實測值。

太原站(37.71°,112.43°)屬于中緯度地區，桂林站(25.19°,110.31°)屬于低緯度地區，低緯度參考站的格網模型校正精度相對較高，說明格網模型的改正性能和空間緯度密切相關[16]。

圖5 太原站3種模型電離層延遲對比

圖6 桂林站3種模型電離層延遲對比

模型最大值/m最小值/m平均值/m雙頻校正11.74244.90089.0127反距離加權15.20397.507012.1416KRG14.765810.266011.8209

表2 桂林站3種模型的精度情況

根據參考站實時監測的穿透點數據，分別采用上節敘述兩種內插法計算GIVE值。GIVE值在通常低于10 m ，最佳保持在1. 5 m 以下[17]。

圖7和圖8為格網點(27. 5°N，110°E)在兩站監測時段內的GIVE的分布情況。可以看出，兩站使用KRG內插法格網模型的修正誤差均在正常范圍內，驗證了算法的可行性。太原地區KRG內插法格網模型GIVE值相對更小，說明在中緯度地區更能凸顯KRG內插法的優勢。表1也可以看出KRG內插法格網模型對電離層延遲值的改正結果更接近于雙頻校正值。由此可見，在使用格網模型校正效果不佳的中緯度地區，使用KRG內插法格網模型修正電離層延遲誤差的精度相對更準確。

圖7 太原站兩種模型GIVE對比

圖8 桂林站兩種模型GIVE對比

4 結論

本文對Klobuchar模型的振幅和夜間延遲值進行不斷改進，得到誤差數據補償定位，使定位更加精確。考慮到與太陽活動強度有關的相對太陽黑子數這一因素，對Klobuchar模型的周期進行修正，建立了改進的Klobuchar電離層延遲模型，并雙頻模型進行評估。此外，還研究了北斗星基增強系統中格網電離層延遲算法以及計算格網點垂直延遲值的方法。以參考站電離層改正方法為基礎，采用KEG內插法計算格網點垂直延遲值的計算方法進行改進，計算格網點GIVE值。結果表明：

1)改進的Klobuchar模型比原模型提高了修正效果，更接近雙頻實測電離層延遲值，從而提高定位精度。

2)空間緯度是影響格網電離層模型修正延遲性能的因素之一。對于低緯度地區，兩種格網內插算法校正精度基本相當，改正誤差GIVE值也均在正常范圍內。而中緯度地區格網模型改正精度相對較差，但使用KRG內插法改正電離層延遲誤差的精度略高于反距離加權法，而且GIVE值更可靠。因此，中緯度地區可采用KRG內插法的格網模型計算用戶電離層延遲誤差。