許 艷 楊元喜 許國昌

1)長安大學(xué)地質(zhì)工程與測繪學(xué)院，西安 710054

2)中國衛(wèi)星導(dǎo)航定位與應(yīng)用管理中心，北京 100088

3)地理信息工程國家重點實驗室，西安 710054

4)德國地球科學(xué)研究中心，波茨坦14473

1 引言

對流層延遲是影響GPS 定位精度的關(guān)鍵因素［1］。對流層延遲的處理方法有對流層模型改正法，引入描述對流層影響的附加待估參數(shù)法［2－4］等。在高精度GPS 測量中，由于受模型誤差、氣象元素測定誤差，尤其是氣象元素代表性誤差的影響，直接建立對流層延遲模型難以取得理想效果。如果將模型求得的值作為近似值，通過較嚴(yán)格的平差計算來估計其精確值是一種較好的方法［5］。

南極地區(qū)GNSS 定位具有與內(nèi)陸地區(qū)不同的特性。文獻(xiàn)［6］對武漢站和南極地區(qū)測站進(jìn)行對比，發(fā)現(xiàn)與中低緯度測站相比，極地定位具有更多的可視衛(wèi)星，但多數(shù)為低高度角。而相關(guān)研究表明對流層延遲對低高度角衛(wèi)星觀測的影響較大［7，8］。所以在極地GNSS 定位中需要尋求更有效的改正對流層延遲的方法以提高定位精度。本文分析了利用標(biāo)準(zhǔn)氣象元素與實際氣象數(shù)據(jù)計算的對流層改正數(shù)的差別，并進(jìn)行GPS 定位。結(jié)果表明，使用實際氣象數(shù)據(jù)模型，在保持同樣定位精度的基礎(chǔ)上可降低觀測數(shù)據(jù)的截止高度角，即可以利用更多的低高度角的衛(wèi)星觀測信息，提高觀測冗余，提高導(dǎo)航定位的精度。

2 對流層延遲改正模型

對流層延遲可表示為:

式中，δ 為對流層總延遲，δdz為天頂方向?qū)α鲗痈煞至垦舆t，F(xiàn)d為相應(yīng)的對流層干分量投影函數(shù)，δwz為天頂方向?qū)α鲗訚穹至垦舆t，F(xiàn)w為相應(yīng)的對流層濕分量投影函數(shù)［3］。

目前，在精密單點定位中常用的對流層模型有Saastamoinen、Hopfield，而投影函數(shù)模型有Neill、VMF1、GMF 等［3］。本文中對流層延遲改正使用Saastamoinen 模型和Neill 投影函數(shù)模型。

Saastamoinen 模型為［5］:

式中，z 為衛(wèi)星的天頂角，P是測站大氣壓，e是水汽壓，T是絕對溫度，δ 為對流層延遲總量。cosz表示投影函數(shù)，可以將其替換為其他的投影函數(shù)，本文中將其替換為Neill 投影函數(shù)。B 為附加的改正項，是測站高程的列表函數(shù)。

在對流層延遲改正模型中要用到測站上的氣象元素，實踐中一般都使用標(biāo)準(zhǔn)氣象元素，即依據(jù)海平面上一組標(biāo)準(zhǔn)的氣象元素值和測站高程H 求得測站上的氣象元素值。

3 對流層延遲改正模型對南極定位精度的影響

3.1 改進(jìn)的對流層延遲模型

實驗數(shù)據(jù)為2011年2月22日丹麥國家空間中心(DNSC)在南極地區(qū)航空重力觀測中的地面靜態(tài)參考站rot1 的數(shù)據(jù)。數(shù)據(jù)采樣率為1 秒，觀測時間為08:40:38—16:53:27，約8 小時12 分鐘。

對南極靜態(tài)參考站rot1 進(jìn)行PPP 處理時，分別使用標(biāo)準(zhǔn)氣象元素法和rot1 站當(dāng)天的實際氣象數(shù)據(jù)進(jìn)行對流層延遲改正模型的計算(取截止高度角為15°)，得到其對流層延遲改正值分別為δstandard和δreal，以及兩種改正值之差δreal－δstandard如表1。

表1 標(biāo)準(zhǔn)氣象元素和實際氣象數(shù)據(jù)計算對流層延遲改正值及兩者結(jié)果之差(單位:m)Tab.1 Tropospheric delays calculated by standard meteorological data and by real meteorological data(unit:m)

從表1 可以看出，分別使用標(biāo)準(zhǔn)氣象元素和實際氣象數(shù)據(jù)計算得到的對流層延遲存在明顯差別，最大可達(dá)0.25 m。通過比較，可以得到兩者間的比例因子α 為:

用α 對原對流層改正模型δ 進(jìn)行處理，得到改進(jìn)后的對流層改正模型δ'為:

3.2 兩種對流層延遲模型對南極定位精度的影響比較

3.2.1 兩種模型僅進(jìn)行對流層模型改正

為了驗證比較改進(jìn)后的對流層模型δ'對定位精度的影響，采用精密單點定位技術(shù)對rot1 站數(shù)據(jù)進(jìn)行處理時分別使用原對流層模型δ 和改進(jìn)后的對流層模型δ'。在精密單點定位中，一般是將精密星歷和精密衛(wèi)星鐘差代入觀測方程，固定衛(wèi)星軌道和消去衛(wèi)星鐘差項，并且采用雙頻觀測值消除電離層影響，故計算的數(shù)學(xué)模型為［2］:

式中，c 為真空中光速，Δt 為接收機鐘差，δ 為對流層延遲影響，α 為比例因子，εΦ為多路徑、觀測噪聲等未模型化的誤差影響，Φ 為消除了電離層影響的組合觀測值，υΦ為觀測誤差，λ 為波長，ρ 為信號發(fā)射時刻的衛(wèi)星位置到信號接收時刻接收機位置之間的幾何距離，N 為消除了電離層影響的組合觀測值的整周未知數(shù)。

對誤差方程進(jìn)行線性化得:

其中P 為觀測權(quán)矩陣，且

整理可得

當(dāng)分別使用對流層模型δ 和δ'計算時，其對流層模型誤差δ－δ'對定位結(jié)果的影響為

在截止高度角取10°和15°時分別使用原對流層模型δ 和改進(jìn)后的對流層模型的計算結(jié)果如表2。

分析表2 可以得到:

1)若采用原對流層模型δ 進(jìn)行對流層延遲改正，與截止高度角取10°的情況相比，截止高度角取15°時其定位精度在X、Y、Z 方向上分別提高了3.49%、8.76%和1.10%。

2)在截止高度角取15°的情況下，當(dāng)采用改進(jìn)后的對流層模型δ'進(jìn)行對流層改正時，其定位精度與使用原對流層模型δ 結(jié)果相比，在X、Y、Z 方向上分別提高了約2.81%、3.03%和2.65%。在截止高度角取10°的情況下，當(dāng)采用改進(jìn)后的對流層模型δ'進(jìn)行對流層改正時，其定位精度與使用原對流層模型δ 結(jié)果相比，在X、Y、Z 方向上分別提高了約10.08%、10.60%和10.06%。在截止高度角取5°的情況下，采用改進(jìn)后的模型δ'進(jìn)行對流層延遲改正，其定位精度在X、Y、Z 方向上為0.024 6、0.020 6、0.070 9 m，與截止高度角取15°時采用原模型δ 進(jìn)行對流層延遲改正的定位精度相當(dāng)。

根據(jù)分析可得采用對流層改進(jìn)模型δ'可以使南極定位在取相同截止高度角的情況下提高定位精度;而在保持同等定位精度的條件下，截止高度角可以從常規(guī)的15°降低到5°。

表2 兩種對流層改正模型的定位精度(單位:m)Tab.2 Positioning precision calculated by two tropospheric corrected models(unit:m)

3.2.2 兩種模型進(jìn)行對流層模型改正和參數(shù)估計

利用模型進(jìn)行對流層改正后，干分量的部分可以改正90%以上，而濕分量部分的殘余影響還比較大［9－12］。因此，在精密定位中，我們一般利用參數(shù)估計方法將對流層延遲當(dāng)作待定參數(shù)進(jìn)行估計求解，而將改正模型所求得的對流層延遲改正視為一種初始近似值。

分別使用對流層模型δ 和改進(jìn)后的模型δ'進(jìn)行對流層延遲模型改正，并將該模型改正值當(dāng)作先驗值，再將對流層延遲當(dāng)作待定參數(shù)進(jìn)行序貫最小二乘［1，13］求解。其處理模型如下:

設(shè)觀測方程為

相應(yīng)的誤差方程為

式中，L 為n ×1 階觀測向量，A 為n ×m 階設(shè)計矩陣，Δ 為n×1 階誤差向量，V 為改正數(shù)向量，X 為包含對流層參數(shù)在內(nèi)的m×1 階具有先驗統(tǒng)計特性的未知數(shù)向量，為估值向量。假設(shè)參數(shù)驗前期望，協(xié)方差矩陣QΔ)，且觀測值獨立，L 與獨立。依據(jù)經(jīng)典最小二乘(LS)貝葉斯估計原理可得

或

驗后協(xié)方差陣為

對流層模型誤差對參數(shù)貝葉斯估式的影響函數(shù)為

分別將兩種模型改正值當(dāng)作先驗值進(jìn)行參數(shù)估計后得到兩者的定位結(jié)果X、Y、Z 坐標(biāo)差值結(jié)果如表3。

從表3 可以看出，分別將采用模型δ 和改進(jìn)后的模型δ'計算得到的對流層延遲值作為先驗值，再

表3 兩種模型對對流層延遲參數(shù)聯(lián)合求解定位結(jié)果(單位:m)Tab.3 Comparison of positioning results between models and in co-determination of the tropospheric parameter(unit:m)

進(jìn)行對流層延遲參數(shù)估計求解，得到兩者的定位結(jié)果X、Y、Z 坐標(biāo)差值很小，其中X 方向上差值約1 mm，Z 方向上差值約4 mm，Y 方向上基本一致。將兩者定位結(jié)果取平均值作為標(biāo)準(zhǔn)值，分別與使用模型δ 和改進(jìn)后的模型δ'進(jìn)行對流層延遲模型改正后的定位結(jié)果見表4。

從表4 可以看出，相對于使用模型δ 進(jìn)行對流層改正，使用改進(jìn)后的模型δ'改正對流層延遲，能使定位結(jié)果更精確，更接近標(biāo)準(zhǔn)值。

表4 δ 模型和δ'模型改正定位結(jié)果與標(biāo)準(zhǔn)值的比較(單位:m)Tab.4 Comparison of positioning results δ and δ' standard values calculated by modelsand(unit:m)

4 結(jié)論

1)取相等截止高度角的情況下，使用實際氣象數(shù)據(jù)的對流層延遲改正模型可以提高定位精度。

2)當(dāng)截止高度角從15°降低為5°時，通過使用實際氣象數(shù)據(jù)的對流層延遲改正模型，兩者可以保持基本同等的定位精度，這樣就可以利用更多的低高度角的衛(wèi)星觀測，提高觀測冗余度。

3)相對于使用標(biāo)準(zhǔn)氣象元素進(jìn)行對流層模型改正，使用實際氣象數(shù)據(jù)改進(jìn)的對流層模型進(jìn)行對流層延遲改正，能使定位結(jié)果更精確，更接近參考值。因此，對于高精度要求的用戶，應(yīng)該采用實際氣象數(shù)據(jù)進(jìn)行對流層延遲的改正。

1 戴吾蛟，等.區(qū)域精密對流層延遲建模［J］.武漢大學(xué)學(xué)報(信息科學(xué)版)，2011，36(4):392－396.(Dai Wujiao，et al.Model costructing for regional precise tropospheric delay［J］.Geomatics and Information Science of Wuhan University，2011，36(4):392－396)

2 葉世榕，張雙成，劉經(jīng)南.精密單點定位方法估計對流層延遲精度分析［J］.武漢大學(xué)學(xué)報(信息科學(xué)版)，2008，33(8):788－791.(Ye Shirong，Zhang Shuangcheng and Liu Jingnan.Analysis of the precision for tropospheric delay estimation by precise point positioning［J］.Geomatics and Information Science of Wuhan University，2008，33(8):788－791)

3 聶建亮.GPS 精密單點定位算法及故障研究［D］.長安大學(xué)，2010.(Nie Jianliang.Study on the algorithm and hitch in precise point positioning［D］.Changan University，2010)

4 葛茂榮，劉經(jīng)南.GPS 定位中對流層折射估計研究［J］.測繪學(xué)報，1996，25(4):285－291.(Ge Maorong and Liu Jingnan.Study on the tropospheric refraction(delay)estimation in global positioning system［J］.Acta Geodaetica et Cartographica Sinica，1996，25(4):285－291)

5 李征航，黃勁松.GPS 測量與數(shù)據(jù)處理［M］.武漢:武漢大學(xué)出 版 社，2010.(Li Zhenghang and Huang Jingsong.GPS surveying and data processing［M］.Wuhan:Wuhan University Press，2010)

6 Xu Yan，Yang Yuanxi and Xu Guochang.Precise determination of GNSS trajectory in the Antarctic airborne kinematic positioning［A］.Lecture Notes in Electrical Engineering——China Satellite Navigation Conference 2012 Proceedings［C］.2012，159:95－105.

7 佘娣，謝劭峰，彭家頔.衛(wèi)星高度角變化對大氣折射的影響［J］.全球定位系統(tǒng)，2011，2:25－28.(She Di，Xie Shaofeng and Peng Jiadi.The influence of satellite elevation angle changing on atmospheric refraction［J］.GNSS World of China，2011，2:25－28)

8 任超，等.低高度角衛(wèi)星信號對提高對流層估計精度的影響分析［J］.大地測量與地球動力學(xué)，2011，(6):124－127.(Ren Chao，et al.Effects of sign of the low elevation angle GPS satellite on estimation of tropospheric delay［J］.Journal of Geodesy and Geodynamics，2011，(6):124－127)

9 Brunner F K and Welsch W M.Effect of the troposphere on GPS measurements［J］.GPS World，1993，4(1):42－51.

10 Collins J P and Langley R B.Estimating the residual tropospheric delay for airborne differential GPS positioning［R］.ION GPS 1997，Kansas City，1997.

11 Tsujii T，et al.Estimation of residual tropospheric delay for high-altitude vehicles:towards precise positioning of a stratosphere airship［R］.ION GPS 2000，Utah，2000.

12 楊元喜.自適應(yīng)動態(tài)導(dǎo)航定位［M］.北京:測繪出版社，2006.(Yang Yuanxi.Adaptive navigation and kinematic positioning［M］.Beijing:Surveying and Mapping Press，2006)

13 Yang Yuanxi.Robust bayesian estimation［J］.Bulletin Geodesique，1991，65(03):145－150.