劉 鈍,甄衛(wèi)民,馮 健,鄧忠新,馬寶田

(中國(guó)電波傳播研究所,山東青島266107)

0 引 言

電離層閃爍是影響衛(wèi)星導(dǎo)航系統(tǒng)定位性能的重要因素之一。電離層閃爍對(duì)衛(wèi)星導(dǎo)航系統(tǒng)的影響包括三個(gè)方面:對(duì)接收機(jī)內(nèi)部環(huán)路的影響,對(duì)接收機(jī)定位性能的影響,對(duì)衛(wèi)星導(dǎo)航系統(tǒng)總體性能的影響。電離層閃爍可以造成接收機(jī)接收到衛(wèi)星信號(hào)的信噪比下降,從而影響接收機(jī)的環(huán)路跟蹤精度,引起載波相位測(cè)量和偽碼測(cè)量的精度降低,閃爍嚴(yán)重時(shí)會(huì)引起接收機(jī)環(huán)路的失鎖,尤其是載波跟蹤環(huán)路的失鎖。偽距測(cè)量誤差的增大,以及衛(wèi)星信號(hào)失鎖引起的用戶定位可用DOP的增大,將影響用戶的定位精度。此外,與電離層閃爍相關(guān)的電離層不均勻性的存在將引起背景電離層延遲的增大,進(jìn)而影響單頻用戶的定位精度。接收機(jī)定位誤差的增大將影響系統(tǒng)定位性能的實(shí)現(xiàn),并進(jìn)一步影響系統(tǒng)的完好性、連續(xù)性和可用性的實(shí)現(xiàn)。

因此,電離層閃爍的影響是衛(wèi)星導(dǎo)航系統(tǒng)應(yīng)用領(lǐng)域關(guān)注的一個(gè)重要問題。電離層閃爍的影響可通過對(duì)閃爍期間用戶接收機(jī)定位性能的變化進(jìn)行分析獲得[1]。由于衛(wèi)星導(dǎo)航系統(tǒng)的復(fù)雜性,以及電離層閃爍隨時(shí)間、地點(diǎn)的變化特性,必須通過仿真方法開展進(jìn)一步的電離層閃爍對(duì)系統(tǒng)總體性能的影響分析[2]。我國(guó)南方低緯地區(qū)是電離層閃爍的高發(fā)區(qū),因此,針對(duì)我國(guó)區(qū)域開展電離層閃爍對(duì)衛(wèi)星導(dǎo)航系統(tǒng)性能的影響研究具有重要意義。

在電離層閃爍對(duì)衛(wèi)星導(dǎo)航系統(tǒng)影響的仿真中,建立閃爍情況下的接收機(jī)環(huán)路誤差模型,仿真分析接收機(jī)環(huán)路測(cè)量誤差隨電離層閃爍參量的變化;建立用戶偽距測(cè)量誤差模型,利用用戶接收機(jī)定位算法獲得用戶的定位誤差變化;利用電離層閃爍模型,給出電離層閃爍參量(如閃爍指數(shù)S4)隨時(shí)間、季節(jié)、地理位置的變化分布,仿真分析衛(wèi)星導(dǎo)航系統(tǒng)性能受電離層閃爍影響的分布特性,并作為閃爍對(duì)系統(tǒng)影響分析的基礎(chǔ)。

1 電離層閃爍對(duì)接收機(jī)的影響仿真

1.1 無(wú)電離層閃爍情況下接收機(jī)跟蹤環(huán)路的誤差模型

衛(wèi)星導(dǎo)航接收機(jī)接收信號(hào)的載噪比C/N0可由如下模型獲得

式中:C為接收機(jī)天線端接收到的最小信號(hào)強(qiáng)度,該值可通過衛(wèi)星導(dǎo)航系統(tǒng)相應(yīng)的ICD[3]文件獲得;B為實(shí)際測(cè)量值與ICD文件中規(guī)定的最小信號(hào)強(qiáng)度之間的差值;對(duì)于GPS系統(tǒng)一般取為3 dB;Ga為天線增益;L為接收機(jī)接收通道處理中的信號(hào)損耗;N0為噪聲功率;I為干擾功率;仿真中不考慮干擾的影響,即I=0.

式中:Bn為PLL環(huán)路或DLL環(huán)路的帶寬;η接收機(jī)預(yù)檢測(cè)積分時(shí)間;d為相關(guān)器間距。c/n0=100.1C/N0.

1.2 電離層閃爍情況下接收機(jī)跟蹤環(huán)路的誤差模型

電離層閃爍對(duì)接收機(jī)環(huán)路影響的仿真中,建立接收機(jī)端閃爍信號(hào)的模型,將接收機(jī)環(huán)路作為一個(gè)系統(tǒng)并建立相應(yīng)的通道模型,利用隨機(jī)信號(hào)通過線性系統(tǒng)的分析方法獲得環(huán)路跟蹤誤差與閃爍參量、接收機(jī)參量間的關(guān)系,并作為仿真分析的基礎(chǔ)。

文獻(xiàn)[2]中給出了典型GPS接收機(jī)結(jié)構(gòu)中,電離層閃爍對(duì)接收機(jī)環(huán)路跟蹤誤差的影響,這里直接利用文中給出的環(huán)路誤差公式,并忽略掉接收機(jī)晶振噪聲引起的環(huán)路跟蹤誤差,具體公式推導(dǎo)參考相關(guān)文獻(xiàn)[2]。

閃爍情況下,接收機(jī)載波跟蹤環(huán)路(PLL)和偽碼跟蹤環(huán)路(DLL)的誤差可以表示為

式中:T為相位閃爍譜中1 Hz處的譜強(qiáng)度;p為相位閃爍譜的譜指數(shù);S4為電離層閃爍指數(shù);k為接收機(jī)跟蹤環(huán)路階數(shù),仿真中取為3;fn為環(huán)路自然頻率;Bn為PLL環(huán)路或DLL環(huán)路的帶寬;η為接收機(jī)預(yù)檢測(cè)積分時(shí)間。

進(jìn)一步的可以由偽碼跟蹤環(huán)路DLL的誤差στ獲得以m為單位的偽碼測(cè)量誤差,

WCA=293.0523 m,為GPS CA碼的碼元長(zhǎng)度。

2 電離層閃爍對(duì)衛(wèi)星導(dǎo)航系統(tǒng)的影響仿真

2.1 用戶測(cè)量偽距的精度評(píng)估

仿真中的偽距測(cè)量精度模型參考文獻(xiàn)[4]與文獻(xiàn)[5],并采用以下的模型估計(jì)獲得

F為電離層傾斜轉(zhuǎn)換因子

式中:RE,h分別為地球半徑和電離層球殼模型的高度(一般取為350 km);E為衛(wèi)星觀測(cè)仰角。

2.2 用戶定位結(jié)果的精度評(píng)估

利用大地測(cè)量學(xué)公式可以獲得衛(wèi)星導(dǎo)航用戶的定位結(jié)果及其精度的評(píng)估[6]。

用戶定位測(cè)量的誤差方程為

用戶定位解算一般采用最小二乘方法獲得,解得的定位結(jié)果,及定位解的精度估計(jì)為

式中:L為觀測(cè)向量;V為殘差向量;X為待求解的未知數(shù),包括用戶三維位置和接收機(jī)鐘差估計(jì);G為狀態(tài)矩陣,與用戶位置和衛(wèi)星位置有關(guān)。仿真中,由于用戶位置已知,且衛(wèi)星位置可由星歷計(jì)算獲得,因此,狀態(tài)矩陣G可以獲得。W為加權(quán)矩陣,可利用(7)式中對(duì)每顆觀測(cè)衛(wèi)星的偽距測(cè)量精度估計(jì)獲得。用戶定位解算方法及精度評(píng)定可參考相關(guān)衛(wèi)星大地測(cè)量學(xué)資料[6]。

3 仿真結(jié)果及分析

3.1 對(duì)接收機(jī)影響的仿真結(jié)果及分析

利用1.2節(jié)中建立的模型,可以獲得在不同閃爍強(qiáng)度下的接收機(jī)跟蹤誤差。其中,譜強(qiáng)度T和譜指數(shù)P通過對(duì)實(shí)際測(cè)量的電離層閃爍數(shù)據(jù)進(jìn)行分析獲得。k為接收機(jī)跟蹤環(huán)路階數(shù),仿真中取為3。fn為相應(yīng)的環(huán)路自然頻率,取為 1.91 Hz。PLL環(huán)路和DLL環(huán)路帶寬分別取10 Hz和0.1 Hz作為典型值。對(duì)于GPS接收機(jī)而言,預(yù)檢測(cè)積分時(shí)間一般為0.02 s(GPS導(dǎo)航電文的數(shù)據(jù)比特率為50 bps),而WAAS接收機(jī)一般為0.002 s(WAAS播發(fā)信息的數(shù)據(jù)比特率為500 bps)。圖1和圖2給出了典型的GPS接收機(jī)和WAAS接收機(jī)分別在閃爍指數(shù)S4為0,0.6,0.705情況下仿真獲得的接收機(jī)環(huán)路跟蹤誤差。

圖1 電離層閃爍對(duì)用戶接收機(jī)PLL環(huán)路的跟蹤誤差影響

從圖中可以看到,隨載噪比的降低和閃爍影響的增大,接收機(jī)環(huán)路跟蹤誤差增大。由于WAAS系統(tǒng)播發(fā)信息的數(shù)據(jù)比特率高于GPS系統(tǒng),WAAS接收機(jī)的預(yù)檢測(cè)積分時(shí)間一般小于GPS接收機(jī)的預(yù)檢測(cè)積分時(shí)間,因此,WAAS接收機(jī)更容易受到電離層閃爍的影響。衛(wèi)星導(dǎo)航系統(tǒng)中應(yīng)合理設(shè)計(jì)下行鏈路的信息量和信息速率,以減小電離層閃爍可能造成的下行信號(hào)的中斷。

圖2 電離層閃爍對(duì)用戶接收機(jī)DLL環(huán)路的跟蹤誤差影響

圖中給出了信號(hào)載噪比大于30 dB-Hz時(shí)接收機(jī)的環(huán)路跟蹤誤差。當(dāng)信號(hào)的載噪比進(jìn)一步降低時(shí),環(huán)路跟蹤誤差將進(jìn)一步增大。但在較低載噪比時(shí),利用上述模型仿真給出環(huán)路跟蹤誤差的意義不大,此時(shí)接收機(jī)一般處于失鎖的臨界狀態(tài),其輸出的偽距具有很大的不確定性。閃爍期間實(shí)際測(cè)量數(shù)據(jù)的分析也表明,此時(shí)單頻用戶定位誤差可以達(dá)到幾百米甚至上千米的量級(jí),定位結(jié)果已不可信[1]。

3.2 對(duì)中國(guó)區(qū)域用戶定位精度的仿真結(jié)果及分析

電離層閃爍對(duì)中國(guó)區(qū)域用戶的定位影響仿真分析中,選定經(jīng)度60°E～140°E,緯度10°N～55°N范圍內(nèi)的區(qū)域,按1°×1°的網(wǎng)格進(jìn)行劃分,每個(gè)網(wǎng)格點(diǎn)作為用戶已知位置,分析電離層閃爍情況下用戶的定位誤差分布,重點(diǎn)分析了單頻定位用戶的受影響情況。仿真中,首先假定沒有電離層閃爍發(fā)生時(shí),該區(qū)域內(nèi)用戶的定位誤差分布,然后考慮存在電離層閃爍時(shí),區(qū)域內(nèi)用戶的定位誤差分布。通過有無(wú)電離層閃爍情況下用戶定位誤差分布的比較分析,說(shuō)明電離層閃爍對(duì)衛(wèi)星導(dǎo)航系統(tǒng)定位精度的影響。

仿真分析中利用GISM電離層閃爍模型給出電離層閃爍指數(shù)S4在上述區(qū)域內(nèi)的分布,仿真的時(shí)間取為2003年10月12日的14:00 UT(對(duì)應(yīng)于北京時(shí)間22:00),模型輸入的F10.7參數(shù)取為110,影響頻率為1575.42 MHz(GPS L1信號(hào)頻率)。圖3為GISM模型給出的中國(guó)區(qū)域電離層閃爍指數(shù)S4的分布。

用戶可視衛(wèi)星是否受電離層閃爍影響,以及電離層閃爍影響的程度,由用戶至衛(wèi)星視線路徑與電離層的穿刺點(diǎn)處的電離層閃爍情況確定,這里假定采用電離層球殼模型,球殼高度取為350 km。由衛(wèi)星信號(hào)受到的電離層閃爍影響進(jìn)一步獲得用戶接收機(jī)對(duì)該衛(wèi)星的偽距跟蹤誤差估計(jì)。

圖3 電離層閃爍指數(shù)S4的區(qū)域分布

GPS衛(wèi)星的星歷誤差估計(jì)中,直接利用了美國(guó)GPS現(xiàn)代化計(jì)劃中精度工程的統(tǒng)計(jì)分析結(jié)果[7],這里取為1.7 m。電離層延遲誤差估計(jì)中,利用CODE提供的電離層圖作為電離層延遲的真實(shí)分布[8]。用戶采用全球單頻電離層模型對(duì)電離層延遲誤差進(jìn)行修正,其精度一般可以達(dá)到60%。σMul通過對(duì)監(jiān)測(cè)站環(huán)境下測(cè)量的偽距精度隨仰角的變化統(tǒng)計(jì)分析得到,一般取為0.5 m。由于對(duì)流層延遲殘差主要由濕延遲引起,而天頂方向濕延遲一般為0.3 m,仿真中取濕延遲的不確定度為30%,σTrop取為0.1 m.

圖4給出了沒有電離層閃爍發(fā)生情況下用戶的定位誤差分布,單位為m?？梢钥闯?沒有電離層閃爍發(fā)生時(shí),仿真區(qū)域內(nèi)的用戶定位誤差分布比較平滑,一般在15 m以內(nèi)。在該區(qū)域的西南部分,定位誤差略大,這是由于仿真時(shí)刻該區(qū)域的可視GPS衛(wèi)星數(shù)較周圍地區(qū)少(如圖5所示)。同時(shí),該區(qū)域的本地時(shí)為傍晚前后(18∶00 LT),電離層還未達(dá)到象東部區(qū)域的平靜狀態(tài)(22∶00 LT)。但該區(qū)域的用戶定位誤差最大也在16 m左右,并未出現(xiàn)異常值。

圖6為電離層閃爍發(fā)生情況下,用戶的定位誤差分布。與圖4相比,可以看出:

1)中緯地區(qū)用戶的定位精度沒有明顯變化,因?yàn)殡婋x層閃爍主要發(fā)生在低緯地區(qū),對(duì)中緯地區(qū)用戶一般沒有影響。

2)用戶定位精度發(fā)生降低的區(qū)域要大于電離層閃爍發(fā)生的區(qū)域。這時(shí)因?yàn)殡m然用戶沒有位于電離層閃爍影響區(qū),但其接收的信號(hào)穿越閃爍區(qū)域,并以較低仰角到達(dá)用戶接收機(jī),由于信號(hào)載噪比較低,信號(hào)受閃爍影響反而更明顯,因此將影響用戶定位精度。

3)低緯地區(qū)用戶定位誤差分布發(fā)生明顯變化,尤其是在東南區(qū)域,用戶定位誤差明顯增大,并存在一個(gè)明顯的定位誤差較大區(qū)域,定位誤差最大值近50 m。

產(chǎn)生該誤差較大區(qū)域的原因在于:

1)電離層閃爍引起用戶接收機(jī)測(cè)量誤差的增大;

從圖2可以看出,在強(qiáng)電離層閃爍情況下(S4＞0.6),電離層閃爍可以造成接收機(jī)偽距跟蹤測(cè)量誤差增大3倍以上,甚至更多。接收機(jī)偽距測(cè)量誤差的增大引起用戶定位誤差的增大。

2)電離層閃爍造成用戶接收機(jī)跟蹤衛(wèi)星的失鎖;

圖7給出了電離層閃爍情況下,用戶可視空間失鎖衛(wèi)星數(shù)的分布。由于接收機(jī)的載波跟蹤環(huán)路更容易受到閃爍的影響[9-10],因此,仿真中采用10°作為接收機(jī)失鎖的門限值,該值是通過對(duì)載波跟蹤環(huán)路失鎖平均時(shí)間進(jìn)行分析獲得的結(jié)論[2]?？梢钥闯?在強(qiáng)電離層閃爍發(fā)生的東南區(qū)域,用戶可視空間衛(wèi)星發(fā)生多顆衛(wèi)星失鎖的現(xiàn)象,最大時(shí)用戶失去4～5顆可視衛(wèi)星,嚴(yán)重影響了用戶定位的DOP值,引起用戶定位誤差的增大。

圖7 電離層閃爍情況下用戶接收機(jī)失鎖的衛(wèi)星數(shù)分布

需要說(shuō)明的是,上述仿真時(shí)間的選擇是基于已有的閃爍觀測(cè)資料而確定的[11]。已有觀測(cè)表明,2003年10月我國(guó)南方區(qū)域有經(jīng)常性的電離層閃爍發(fā)生?？梢灶A(yù)見,在太陽(yáng)活動(dòng)的最高時(shí)期(2012年),我國(guó)南方區(qū)域的電離層閃爍現(xiàn)象可能更強(qiáng)烈。因此,在即將到來(lái)的新的太陽(yáng)活動(dòng)高年中,我們將進(jìn)一步加強(qiáng)實(shí)際的電離層閃爍觀測(cè)工作,并在上述基礎(chǔ)上,結(jié)合實(shí)際測(cè)量開展更深入的分析工作。

上述仿真工作初步建立了電離層閃爍對(duì)衛(wèi)星導(dǎo)航系統(tǒng)影響的分析框架,進(jìn)一步的工作中,將考慮采用利用實(shí)測(cè)數(shù)據(jù)實(shí)現(xiàn)的電離層閃爍模型,以代替基于統(tǒng)計(jì)分析建立的GISM模型。

4 結(jié) 論

利用電離層閃爍模型,結(jié)合接收機(jī)模型和衛(wèi)星導(dǎo)航系統(tǒng)用戶定位模型,可以實(shí)現(xiàn)電離層閃爍對(duì)衛(wèi)星導(dǎo)航系統(tǒng)定位性能影響的仿真分析,研究評(píng)估電離層閃爍對(duì)用戶定位性能的影響,及受影響的范圍分布情況。

電離層閃爍造成導(dǎo)航接收機(jī)跟蹤環(huán)路的誤差增大,影響接收機(jī)的偽距測(cè)量精度;電離層閃爍嚴(yán)重時(shí),可以造成接收機(jī)的失鎖,引起用戶定位中DOP的增大。兩種因素共同影響用戶的定位精度。

我國(guó)低緯地區(qū)的電離層閃爍可以引起用戶較大的定位誤差,尤其是在太陽(yáng)活動(dòng)高年,電離層閃爍可能會(huì)引起一個(gè)區(qū)域內(nèi)用戶定位精度的嚴(yán)重降低,因而是對(duì)該區(qū)域內(nèi)用戶導(dǎo)航應(yīng)用的嚴(yán)重威脅。

在即將到來(lái)的新的太陽(yáng)活動(dòng)高年中,將進(jìn)一步加強(qiáng)實(shí)際的電離層閃爍觀測(cè)工作,考慮利用我國(guó)的電離層閃爍實(shí)測(cè)數(shù)據(jù),建立更真實(shí)的電離層閃爍參量分布模型,并結(jié)合未來(lái)衛(wèi)星導(dǎo)航系統(tǒng)新的信號(hào)體制,進(jìn)一步開展電離層閃爍對(duì)衛(wèi)星導(dǎo)航系統(tǒng)性能影響的分析評(píng)估工作。

[1] 劉 鈍,馮 健,鄧忠新,甄衛(wèi)民.電離層閃爍對(duì)全球?qū)Ш叫l(wèi)星系統(tǒng)(GNSS)的定位影響分析[J].全球定位系統(tǒng),2009,34(6):1-8.

[2] Robert S C,Bakry El-Arini M,Christopher J.Tzong-Yang Hsiao H.Modeling the effects of ionospheric scintillation on GPS/SBAS availability[C]//ION Annual Meeting,2000:563-576.

[3] ICD-GPS-200.Navstar GPS space segment/navigation user interfaces[R].Arinc Research Corporation,1993.

[4] RTCA Special Committee 159.Minimum operational performance standards for airborne equipment using global positioning system/wide area augmentation system,RTCA/DO-229 Change 3[R].Washington:RTCA,1997.

[5] Corrigan T M,Hartranft J F,Levy L J,Parker K E,et al.GPS risk assessment Study(final report)[R].The Johns Hopkins University.Applied Physics Laboratory,1999.

[6] 周忠謨.GPS衛(wèi)星測(cè)量原理與應(yīng)用[M].北京:測(cè)繪出版社,1992.

[7] 李 躍,邱致和.導(dǎo)航與定位[M].北京:國(guó)防工業(yè)出版社,2008.

[8] Stenfan S.M apping and predicting the earth's ionosphere using the global positioning system[D].University of Bern,1999.

[9] Elliott D K.GPS原理與應(yīng)用[M].邱致和,王萬(wàn)義譯,北京:電子工業(yè)出版社,2002.

[10] James Bao-Yen T.GPS軟件接收機(jī)基礎(chǔ)[M].陳軍,潘高峰,等譯.2版.北京:電子工業(yè)出版社,2007.

[11] 甄衛(wèi)民,馮 健,陳 麗,韓一平.多站多路徑GPS信號(hào)研究低緯電離層不均勻體[J].電波科學(xué)學(xué)報(bào),2007,22(1):138-142.