于 耕,李大武，陳志強

(1.沈陽航空航天大學 a.民用航空學院，b.電子信息工程學院,沈陽 110136;2.中國國際貨運航空有限公司 飛行技術(shù)部，北京 100621)

完好性指的是衛(wèi)星導航系統(tǒng)由于自身或外界的影響停止導航服務亦或保護級異常發(fā)出告警信號的情況下,系統(tǒng)能夠自我檢測與報警的能力。 它是衡量導航系統(tǒng)優(yōu)劣的重要指標,對于用戶的生命財產(chǎn)安全有著重大影響。

當前運用廣泛且成熟的完好性算法主要有接收機自主完好性監(jiān)測(Receiver Autonomous Integrity Montior,RAIM)及地面完好性通道(Ground Integrity Channel,GIC),前者屬于用戶終端的性能要求,后者屬于系統(tǒng)頂層設計的范疇[1]。

完好性的好壞主要由報警閾值、告警時間以及完好概率等指標來反應,而無論采用何種完好性算法,針對不同星座數(shù)目計算出的結(jié)果都會不盡相同[2-3]。完好性相關(guān)參數(shù)的計算是當前完好性研究的熱點。某種程度上,導航系統(tǒng)的完好性水平是用戶最為關(guān)注的指標。航空導航的最大特點不僅要達到較高的精度要求,而且更要能達到一定的安全性指標,從而最大滿足用戶的需求。

通過查閱資料了解到當今研究現(xiàn)狀如下:基于北斗星座的完好性保護級算法已有人研究[4],主要做法是綜合利用中國北斗星座以及美國的GPS星座對保護級算法及其結(jié)果進行研究分析;此外也有高校對雙星座下陸基增強系統(tǒng)的完好性算法進行了探究[5];但在北斗星座上應用星基增強系統(tǒng)所做的研究還很少,因此這將是本文研究的重點。

鑒于我國的星基增強系統(tǒng)還處在建設初期,由于目前條件的限制,很多功能尚處在研究模擬階段,以我國自主發(fā)展的北斗衛(wèi)星導航系統(tǒng)[6-7]為基礎,驗證了星基增強系統(tǒng)在中國部分地區(qū)保護級結(jié)果運用的可行性。文中重點探究了桂林地區(qū)不同抽樣時刻的保護級水平以及調(diào)取的6個參考站在同一時刻的保護級水平對比,分析了勉縣地區(qū)一段時間內(nèi)保護級變化水平狀況。利用中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡的10個GNSS參考站的觀測數(shù)據(jù),仿真分析了完好性重要參數(shù)水平保護級(HPL)與垂直保護級(VPL)。

1 完好性保護級算法原理

1.1 星基增強系統(tǒng)及完好性生成概述

北斗星基增強系統(tǒng)的構(gòu)成如圖1所示。

北斗星基增強系統(tǒng)由四個部分組成,分別為空間部分、地面端、用戶端和數(shù)據(jù)鏈路?？臻g部分指的是北斗系統(tǒng)的空間星座,包括地球同步衛(wèi)星(GEO)、中圓軌道衛(wèi)星(MEO)和傾斜地球同步軌道衛(wèi)星(IGSO)。其中所有衛(wèi)星均可向用戶發(fā)送基本導航信息,而GEO衛(wèi)星在發(fā)送基本導航信息的同時還可向用戶發(fā)送北斗系統(tǒng)的廣域差分信息、格網(wǎng)電離層信息和完好性信息。地面部分包括監(jiān)測站、中心處理站和上傳注入站。參考站使用雙頻接收機接收fB1、fB2頻率碼和載波相位觀測量的同時還采集氣象參數(shù),并經(jīng)過預處理后送入中心處理站[8]。中心處理站將收集到的數(shù)據(jù)進行計算得到誤差修正信息和完好性信息,并通過上傳注入站發(fā)送給地球同步通信衛(wèi)星。用戶端的主要任務是同時接收衛(wèi)星的導航信息和GEO廣播的差分及完好性信息,并利用這些信息進行定位解算、完好性分析等。數(shù)據(jù)鏈路包括參考站與中心處理站之間的通信鏈路和衛(wèi)星與用戶之間的通信鏈路。其中參考站與中心處理站由公用數(shù)據(jù)傳輸網(wǎng)、電話網(wǎng)或?qū)Ｓ猛ㄐ欧绞竭B接;衛(wèi)星與用戶之間采用衛(wèi)星通信方式,廣播的信息需按一定格式進行編碼。

圖1 BD-SBAS系統(tǒng)構(gòu)成

SBAS完好性的計算流程描述如下:

(1)從觀測衛(wèi)星那里可以獲知航空電子輸入模擬器的飛行器動態(tài)以及星座信息,結(jié)合數(shù)據(jù)庫中已有的掩飾矩陣構(gòu)成數(shù)據(jù)基礎;

(2)從觀測衛(wèi)星那里可以得出快速及長期誤差模型與機載接收機誤差模型,而由觀測衛(wèi)星與飛行器動態(tài),又可以得到電離層與對流層誤差模型;

(3)以上模型構(gòu)成總的誤差模型,由總誤差模型結(jié)合觀測衛(wèi)星信息就可以得到相關(guān)計算時用到的矩陣,分別是觀測矩陣、加權(quán)矩陣和投影矩陣;

(4)利用總誤差模型與矩陣信息計算出垂直與水平保護級。

SBAS系統(tǒng)的導航電文信息包括本衛(wèi)星基本導航信息、全部衛(wèi)星歷書、與其他系統(tǒng)時間同步信息、北斗系統(tǒng)完好性及差分信息和格網(wǎng)點電離層信息[9-10]。但基于中國北斗的星基增強系統(tǒng)由于尚處于建設初期,地球同步衛(wèi)星(GEO)還沒能播發(fā)增強導航信息的相關(guān)參數(shù),所以不能針對保護級和完好性進行實際驗證。因此只能依據(jù)當前條件,對系統(tǒng)所要求的性能進行理論研究與模擬分析,從而在一定程度上驗證了SBAS系統(tǒng)在中國的可行性。

1.2 完好性保護級的計算

1.2.1 用戶差分偽距誤差UDRE

用戶差分偽距誤差UDRE是對各個監(jiān)測站點接收機所接收到的數(shù)據(jù)進行衛(wèi)星天線相位中心改正、接收機鐘差改正、格網(wǎng)電離層延遲改正、對流層改正、相對論效應改正、地球自轉(zhuǎn)改正、固體潮和海潮等潮汐改正后,可得到改正后的偽距Rm,同時可根據(jù)廣播星歷和已知的監(jiān)測站坐標得到衛(wèi)星與監(jiān)測站之間的距離Rc,則兩者之差dR即是UDRE。

dR=Rm-Rc

(1)

用戶差分偽距誤差指的是:給星基增強系統(tǒng)服務區(qū)域內(nèi)的用戶提供各項改正誤差的綜合參數(shù),且在一定置信度內(nèi)保障用戶安全。由監(jiān)測站得到的觀測量通過地球自轉(zhuǎn)改正、廣域差分格網(wǎng)電離層延遲改正、對流層延遲改正等各項改正后,可得到監(jiān)測站與各個衛(wèi)星之間的偽距,該偽距稱為觀測偽距,用Rm表示。由監(jiān)測站自身的已知坐標及通過廣播星歷得到的衛(wèi)星三維坐標可知監(jiān)測站與各個衛(wèi)星之間的偽距稱為計算偽距,用Rc表示。

對各個監(jiān)測站的觀測量進行數(shù)據(jù)處理,可得到同一時刻不同監(jiān)測站與相同衛(wèi)星的dR。通過對dR在一定置信度下進行分析,可得到該時刻相應衛(wèi)星的UDRE值,即

(2)

(3)

用戶通過計算UDRE值可得到當前偽距的誤差限值,并根據(jù)UDRE值對解算的位置信息進行完好性判斷[11]。

由于在UDRE的計算過程中,各項改正及一些偶然誤差無法徹底消除,導致UDRE的估計出現(xiàn)部分偏差,無法和理論真值保持一致。此偏差同時會影響UDRE及相關(guān)導航性能的降低。但對于局部誤差可采取平滑偽距及限制高度角等方法進行改正[12],計算偽距差時可加入如下限制條件:

(1)保證觀測到相同衛(wèi)星的監(jiān)測站數(shù)目為一個以上,否則認為UDRE為“未被監(jiān)測”;

(2)如果觀測到的相同衛(wèi)星且高度角大于15°的監(jiān)測站數(shù)目為2個以上(含兩個),則只使用大于15°的觀測值;

(3)如果觀測到的相同衛(wèi)星且高度角大于15°的監(jiān)測站數(shù)目不足2個,則只使用觀測到的衛(wèi)星高度角最大的2個觀測值;

UDRE每秒計算一次,按更新率需求(如6秒)周期性廣播給用戶。

UDRE分成16等級,按預先設定的等級關(guān)系,計算出對應的等級值UDREI,發(fā)布給用戶參考。等級值中包含“未監(jiān)測”、“不可用”等信息,見表1。

表1 UDRE等級分列

1.2.2 SBAS保護等級公式

保護級的計算公式如式(4)所示。

(4)

VPLSBAS=KVdU

這里的幾個K都是固定值,其中KH,NPA代表非精密進情況下的取值,為6.18,此時系統(tǒng)根據(jù)水平導航來引導用戶進近;KH,PA指有水平/垂直方向?qū)Ш揭嗷蛴?無垂直向?qū)нM近時的取值,為6.0;KV值為5.33。dmajor代表隨機誤差橢圓半長軸,具體公式如式(5)所示。

(5)

(6)

依次表示為:

東軸線上超過真實誤差模型的分布方差;

北軸線上超過真實誤差模型的分布方差;

東軸和北軸線上模型分布的協(xié)方差;

縱軸上超過真實誤差模型的分布方差。

另外,seast,i代表東方向上第i顆衛(wèi)星偽距的定位誤差偏導數(shù);

snorth,i代表北方向上第i顆衛(wèi)星偽距的定位誤差偏導數(shù);

sU,i代表垂直方向上第i顆衛(wèi)星偽距的定位誤差偏導數(shù)。

(7)

投影矩陣S具體表示為

(GT·W·G)-1·GT·W

(8)

第i行幾何矩陣G:

Gi=[-cosElisinAzi-cosElicosAzi-sinEli1]

(9)

其中Eli、Azi表示衛(wèi)星i的仰角和方位角。

加權(quán)矩陣W表示為

(10)

1.3 總方差分析與計算

(11)

這里δUDRE指代用戶差分偽距誤差系數(shù),通常取1,σi,UDRE則表示用戶差分偽距方差。

(12)

(13)

第三個,機載接收機方差的值為

(14)

最后是對流層誤差方差的計算公式如式(15)所示。

σi,tropo=(σTVE·m(Eli))

σTVE=0.12 meters

(15)

公式中的m(Eli)代表對流層校正映射函數(shù)。

2 數(shù)據(jù)仿真與分析

基于以上保護級算法的分析,再結(jié)合采集到的中國大陸環(huán)境監(jiān)測網(wǎng)的監(jiān)測數(shù)據(jù),對中國區(qū)域內(nèi)10個監(jiān)測站附近接收到的衛(wèi)星數(shù)據(jù)進行仿真,計算得出相關(guān)地區(qū)某時刻與一段時間內(nèi)的保護級水平及其分析。圖2給出了10個監(jiān)測站在中國大陸的地理分布,如圖所示基本上可以覆蓋中國全境。

圖2 選取的監(jiān)測站分布圖

2.1 桂林參考站不同時刻保護級水平

數(shù)據(jù)采集時間為當日的北京時間1時到4時,采樣間隔為1 s。在桂林監(jiān)測站隨機抽取了5個時刻,對保護級水平進行了仿真,旨在探究保護級的隨機性結(jié)果是否符合實際應用。結(jié)果如圖3所示。

圖3 桂林參考站各采樣點保護級

圖中左邊黑色條形為非精密進近時的水平保護級數(shù)值；中間灰色條形為精密進近在水平或垂直導航亦或有/無垂直引導定位器情況下的水平保護級；右邊白色條形為垂直保護級。

可以看出:

(1)兩種情況下的水平保護級的數(shù)值很相近,都在4～5米之間波動,最大、最小值分別為4.748米、4.191米和4.609米、4.069米;

(2)垂直保護級的數(shù)值更小,在2～3米之間波動,最大、最小值分別為2.506米、2.260米。

因為規(guī)定中說明水平報警限值不超過17米,垂直報警限值不超過5.3米即可滿足二、三類進近標準[13],因此仿真結(jié)果表明桂林地區(qū)的保護級水平完全能滿足飛機的進近要求,而且精度很高。

2.2 不同參考站的保護級對比

對各個參考站采集到的數(shù)據(jù)以及同一時刻不同地區(qū)的保護級水平進行了仿真計算與對比,旨在探究相同時刻不同地點接收不同衛(wèi)星播發(fā)的信息對保護級數(shù)值水平以及實際應用的影響。結(jié)果如圖4所示,其中橫坐標1～6分別代表武漢、桂林、勉縣、太原、鹽城、廈門等6地參考站。

圖4 各參考站保護級對比

圖中左邊黑色柱形代表在水平導航進近時的保護級水平；中間灰色柱形代表水平或垂直導航亦或在有/無垂直引導情況下的水平保護級數(shù)值；右邊白色柱形代表垂直保護級水平。

可以看出:

(1)各個參考站的水平保護級數(shù)值很相近,都在4米上下波動且不低于3米,也不高于5米,兩者的最大值、最小值分別為:4.712米、4.030米和4.575米、3.912米。

(2)垂直保護級都在2～3米之間波動,最大值2.438米,最小值2.240米。

各個參考站之間的保護級數(shù)值略有不同,這與參考站位置的選取以及計算過程中產(chǎn)生的誤差有關(guān)[14],整體上得到的數(shù)據(jù)還是很理想的,表明能夠滿足導航進近的高精度要求。

2.3 勉縣地區(qū)連續(xù)保護級水平

由于實際應用中保護級水平是連續(xù)變化的,因此本文選取勉縣地區(qū)監(jiān)測數(shù)據(jù)對其保護級水平的連續(xù)性結(jié)果進行了仿真探討，實驗結(jié)果如圖5所示。

圖中加點曲線代表非精密進近時的水平保護級數(shù)值；黑色曲線代表精密進近在水平或垂直導航亦或有/無垂直引導定位器情況下的水平保護級；淺灰色曲線為垂直保護級。

圖5 勉縣保護級水平

3 結(jié)論

本文以北斗導航系統(tǒng)為基礎,探究了星基增強系統(tǒng)在我國部分地區(qū)的完好性性能,且對主要參數(shù)——水平保護級與垂直保護級分別在離散與連續(xù)性方面進行了研究,并做了仿真計算。結(jié)果表明,星基增強系統(tǒng)適用于中國區(qū)域內(nèi)用戶,而且能保持一定的精度水平,從側(cè)面探究了星基增強系統(tǒng)的實際可行性。本文通過對桂林參考站隨機采樣點的仿真計算對比,以及6個參考站同一時刻的保護級參數(shù)對比,以及勉縣地區(qū)一段時間內(nèi)保護級變化水平分析,得出結(jié)論如下:

(1)BD-SBAS完好性保護級算法能夠很好的處理相關(guān)衛(wèi)星數(shù)據(jù),并最終給出理想的保護級數(shù)值,由此驗證了保護級算法在中國地區(qū)的可行性。

(2)以不同的視角探究了保護級的變化情況,一定程度上真實地反映了保護級的變化情況,對以后BD-SBAS實際應用具有一定的借鑒意義。