陳中新，朱麗強(qiáng)

(1.蘇州工業(yè)園區(qū)格網(wǎng)信息科技有限公司，江蘇蘇州 215021； 2.JSCORS蘇州分中心，江蘇蘇州 215021)

區(qū)域氣象特征對流層延遲模型的建設(shè)研究

陳中新1，2，朱麗強(qiáng)1，2

(1.蘇州工業(yè)園區(qū)格網(wǎng)信息科技有限公司，江蘇蘇州 215021； 2.JSCORS蘇州分中心，江蘇蘇州 215021)

在分析GPS隨機(jī)軟件中兩種常用的對流層延遲模型的基礎(chǔ)上，利用蘇州地區(qū)CORS網(wǎng)絡(luò)數(shù)據(jù)計(jì)算了模型與高度、緯度之間的關(guān)系，提出了在網(wǎng)絡(luò)RTK模式下建立區(qū)域氣象特征的對流層延遲模型的必要性。

對流層延遲；區(qū)域模型；網(wǎng)絡(luò)RTK

1 引 言

對于城市測繪而言，無論是靜態(tài)測量，還是實(shí)時(shí)動(dòng)態(tài)測量，利用站間距比較近而空間誤差相關(guān)性強(qiáng)的假設(shè)，通過觀測值之間求差的方法，可以有效削弱GPS信號傳播過程中的各種誤差。但是由于對流層直接與地面接觸并從地面得到能量，其溫度隨高度的上升而降低、大氣層質(zhì)量的99%集中于此，以及隨站間距的增加對流層延遲路徑上物理性質(zhì)不相似越明顯的特性，使得GPS信號延遲的80%來源于對流層[1]，因此，在與IGS站聯(lián)測、高精度GPS基線處理以及氣候復(fù)雜變化時(shí)的實(shí)時(shí)動(dòng)態(tài)，就不能通過簡單求差減弱對流層延遲的影響，必須采用合適的對流層模型進(jìn)行改正。

2 對流層延遲模型

對流層存在著各種各樣的復(fù)雜的物理過程，從而不斷改變氣壓、溫度、濕度沿高度和水平方向的分布，使得大氣折射率隨時(shí)間和空間而變化[2]。盡管可以采用探空儀探測精確的大氣剖面，但因其測量的時(shí)間長，而大氣結(jié)構(gòu)又隨時(shí)變化等因素，實(shí)踐中還是采用經(jīng)驗(yàn)?zāi)Ｐ透恼?、參?shù)估計(jì)等方法以減小對流層延遲影響。

測量型GPS的隨機(jī)軟件中對流層處理選項(xiàng)中的模型，應(yīng)該是具有代表性的常見模型。其中Trimble公司的TGO選項(xiàng)包括:Hopfield，Goad and Goodman，Saastamoinen，Black模型等；Topcon公司的Pinnicale選項(xiàng)包括:Goad and Goodman，分段常數(shù)法等；Leica公司的LGO選項(xiàng)包括:Hopfield，Saastamoinen，參數(shù)估計(jì)等。

2.1 Hopfield模型

Hopfield把對流層的折射對觀測值的影響分為干分量和濕分量兩部分，其中干分量主要取決于大氣的溫度和壓力，濕分量主要取決于信號傳播路徑上的大氣濕度和高度。引入映射函數(shù)，其模型方程為:

其中T為測站的絕對溫度，以K為單位；P為測站氣壓，以毫巴為單位；e為測站水氣壓，以毫巴為單位；E為衛(wèi)星的高度角，以度為單位；△trop(E)為對流層折射率改正數(shù)，以米為單位。

2.2 Saastamoinen模型

Saastamoinen模型把大氣分成三部分:從地面到10 km左右高度為對流層，假設(shè)其氣溫遞減率為6.5℃/km；從對流層頂?shù)?70 km左右高度為平流層，假設(shè)其氣溫為常數(shù)；70 km以外為電離層。天頂方向干濕延遲分量計(jì)算的簡約模型為:

公式中符號的意義與Hopfield模型一致，f(B，H)為測站所在緯度和高程的函數(shù)。

映射函數(shù)可采用Marini連分式映射函數(shù)，干濕分量映射函數(shù)相同，均為:

干分量映射函數(shù)系數(shù)定義為:

濕分量映射函數(shù)系數(shù)定義為:

3 模型應(yīng)用的驗(yàn)證與分析

3.1 高精度靜態(tài)后處理

城市高精度測量控制網(wǎng)都是采用靜態(tài)測量技術(shù)，為了確保成果的質(zhì)量，一般比規(guī)范要求的觀測時(shí)間增加50%甚至100%，并嚴(yán)格控制觀測過程質(zhì)量，但是往往忽略或不重視基線處理過程中需要的測站氣象元素的測定。因此，在沒有實(shí)測氣象元素的情況下，利用加拿大DIPOP軟件提供的溫度、氣壓與高度的關(guān)系[3]計(jì)算模型的延遲量:

其中T0＝288.15，以度為單位；P0＝1 013.25，以毫巴為單位；e0＝11.691，以毫巴為單位；h為海拔高，以米為單位。因此，對于緯度30°的地區(qū)，由于測站高差、不同衛(wèi)星高度角引起的對流層延遲影響，如表1所示。

緯度30°地區(qū)對流層模型延遲誤差(單位/m) 表1

由表1可知，在同一地區(qū)(緯度變化不大)，海拔較低時(shí)(例如:高程約20 m)在觀測高度角不低于15°的情況下，兩種典型延遲模型間的誤差變化量均在1 cm內(nèi)。在海拔上升到約 500 m時(shí)，兩種典型延遲模型間的誤差變化量在 50 cm左右。而當(dāng)海拔為6 000 m時(shí)，模型誤差間的差值達(dá)到幾米，甚至十幾米。從模型延遲對高度的梯度分析，Hopfield模型的梯度與高度無關(guān)，是溫度的函數(shù)，而Saastamoinen模型與高度有關(guān)，隨著高度的變化比較明顯。隨著測站高程的增大，兩者之差明顯變大。Hopefield模型只按照單層計(jì)算，而Saastamoinen模型將溫度作為常數(shù)分兩層計(jì)算，并對高度的變化反應(yīng)靈敏，因此，其精度要比Hopfield模型改正效果要好。

因此，在不同地區(qū)進(jìn)行GPS網(wǎng)數(shù)據(jù)后處理的時(shí)候，要甄別使用合適的對流層模型。近年來，在大型、高精度的GPS控制網(wǎng)數(shù)據(jù)處理經(jīng)驗(yàn)，一般采用Saastamoinen模型來估計(jì)對流層延遲。為了進(jìn)一步驗(yàn)證模型間的延遲差異，對蘇州地區(qū)平均邊長約120 km的大型控制網(wǎng)，采用隨機(jī)軟件設(shè)置不同大氣模型進(jìn)行基線的處理，其結(jié)果如表2所示?？梢?，在蘇州地區(qū)使用不同的對流層模型對結(jié)果基本沒有影響，成果和表1的理論分析相一致。

不同對流層模型對基線解算的影響 表2

3.2 網(wǎng)絡(luò)RTK實(shí)時(shí)定位

在常規(guī)RTK測量中，由于基站與流動(dòng)站之間的距離一般控制在10 km范圍內(nèi)，相對于GPS信號的傳播路徑來說是個(gè)微小量，因此，在不使用模型的情況下，直接求雙差就能極大削弱對流層延遲。但在網(wǎng)絡(luò)RTK中，基站的間距達(dá)到60 km～100 km，只有當(dāng)流動(dòng)站與基準(zhǔn)站的距離較近、高差較小的時(shí)候，利用差分可以有效減小對流層延遲誤差的影響，由于對流層相關(guān)性減弱，在大多數(shù)時(shí)候通過雙差不能簡單地消除對流層的影響。

由于基站具有全天候、連續(xù)觀測的特點(diǎn)，系統(tǒng)有足夠多的數(shù)據(jù)來精確求解對流層的天頂延遲參數(shù)，從而模擬出區(qū)域范圍內(nèi)的大氣模型，系統(tǒng)控制中心使用的是根據(jù)觀測數(shù)據(jù)不斷修正的高精度的對流層延遲模型。但是對于流動(dòng)站而言，需要通過一兩個(gè)歷元的觀測得到高精度的定位，實(shí)際上當(dāng)流動(dòng)站與基準(zhǔn)站的高差超過6 000 m時(shí)，如果不修正對流層延遲，由對流層延遲引起的天頂方向延遲偏差達(dá)5 m[4]。

對流層延遲的殘差對于流動(dòng)站與CORS網(wǎng)絡(luò)而言，主要是兩者之間的天頂方向相對延遲差，對高程的影響特別顯著。它所引起的高程誤差[5]:

其中△dtrop為天頂方向殘余對流層延遲，Z為天頂距，△h為引起的高程誤差，E為衛(wèi)星截止角。根據(jù)上式和表1計(jì)算在緯度30°的地區(qū)，當(dāng)基站和流動(dòng)站高差不大的時(shí)候，不同模型延遲量的較差對高程的影響，其結(jié)果如表3所示；而當(dāng)基站和流動(dòng)站高差特別大的時(shí)候，同一模型延遲量的較差對高程的影響，其結(jié)果如表4所示。

Hopfield和Saastamoinen對流層延遲量之差對高程的影響(單位/m) 表3

Hopfield模型計(jì)算對流層延遲量對高程的影響(單位/m) 表4

由表3可知，在低海拔地區(qū)，流動(dòng)站與參考站的高差不大的情況下，即使采用不同的對流層延遲模型所引起的高程誤差也在厘米級；而在高海拔地區(qū)當(dāng)采用不同模型時(shí)所引起的高程誤差可能達(dá)到幾十米。由表4可知，當(dāng)流動(dòng)站與參考站之間的高差達(dá)到 6 000 m，即使采用同一對流層延遲模型改正，所引起的高程誤差也要有十幾米。

因此，即使在網(wǎng)絡(luò)RTK中的短基線解算也必須謹(jǐn)慎處理對流層延遲誤差，從而提高實(shí)時(shí)RTK的精度和成果的可靠性，原先的一些經(jīng)典模型改正，其靜力學(xué)延遲部分的改正精度可以達(dá)到厘米級，而濕延遲部分的殘余影響還比較大，一般仍會(huì)有數(shù)厘米的殘差。所以，非常有必要在多部門聯(lián)建的省CORS系統(tǒng)的基礎(chǔ)上，通過對對流層延遲隨時(shí)間的變化隨機(jī)模型[6]，自主修正高程偏差的距離相關(guān)對流層網(wǎng)絡(luò)內(nèi)插模型[7]以及用隨機(jī)游走方法來估計(jì)對流層延遲殘差[8]等方法的實(shí)踐和優(yōu)化，得到符合區(qū)域大氣條件的對流層延遲模型，從而提高實(shí)時(shí)定位精度和成果的可靠性。

5 結(jié) 論

盡管連續(xù)參考站系統(tǒng)已被國內(nèi)公認(rèn)為城市基礎(chǔ)設(shè)施之一，網(wǎng)絡(luò)RTK技術(shù)能實(shí)時(shí)獲取高精度的點(diǎn)位坐標(biāo)而具有非常廣闊的前景，但是對于這一新興的GPS差分技術(shù)，有必要進(jìn)行系統(tǒng)而嚴(yán)密的研究。從目前的一些測試和實(shí)踐表明，一方面對流層延遲誤差仍是該技術(shù)的主要誤差源之一；另一方面由于基準(zhǔn)站和流動(dòng)站獲取的數(shù)據(jù)量以及精度要求不同，因而對流層改正模型也應(yīng)有差異。

因此，建立精確的區(qū)域?qū)α鲗友舆t改正模型越來越必要和迫切，這就需要我們創(chuàng)新性地研究對流層延遲對GPS測量的影響，通過各種模型的比較分析，找出修正對流層延遲較為有效的方法，建立起區(qū)域氣象特征對流層延遲模型，從而可以提高網(wǎng)絡(luò)RTK成果可靠性、改善高精度基線處理的數(shù)據(jù)歸算，以及優(yōu)化基于CORS系統(tǒng)的高精度GPS數(shù)據(jù)仿真。而近年來由測繪局、氣象局、地震局等部門聯(lián)合建設(shè)的GPS連續(xù)觀測網(wǎng)絡(luò)，比如江蘇CORS[9]，能提供的觀測數(shù)據(jù)類型越來越全面，為研究反演計(jì)算GPS信號在對流層中傳播的時(shí)延，建立區(qū)域氣象特征的對流層延遲模型提供了有利的條件。

[1]霍夫曼－韋倫霍夫，利希特內(nèi)格爾，瓦斯勒.程鵬飛，蔡艷輝等譯.全球衛(wèi)星導(dǎo)航系統(tǒng)[M].北京:測繪出版社，2009

[2]謝益溪，J拉菲涅特，JP蒙等.電波傳播—超短波、微波、毫米波[M].北京:電子工業(yè)出版社，1990

[3]曲偉菁，朱文耀，宋淑麗等.三種對流層延遲改正模型精度評估[J].天文學(xué)報(bào)，2008(1):113～122

[4]何海波.高精度GPS動(dòng)態(tài)測量及質(zhì)量控制[D].博士學(xué)位論文.中國人民解放軍信息工程大學(xué)，2002

[5]U.Hugentobler，S.Schaer，P.Fridez.Bernese GPS Software (Version 4.2)[M].Stronomical Institute University ofBerne.2001

[6]邵占英，葛茂榮，劉經(jīng)南.GPS定位中對流層折射率隨機(jī)模型的研究[J].地殼形變與地震，1996，16(2):1～7

[7]熊永良，黃丁發(fā)，丁曉利等.虛擬參考站技術(shù)中對流層誤差建模方法研究[J].測繪學(xué)報(bào)，2006，35(2):118～121，132

[8]劉經(jīng)南，葉世榕.GPS非差相位精密單點(diǎn)定位技術(shù)探討[J].武漢大學(xué)學(xué)報(bào):信息科學(xué)版，2002，27(3):234～240

[9]安艷輝，史照良.全球衛(wèi)星導(dǎo)航定位連續(xù)運(yùn)行參考站網(wǎng)的技術(shù)探討——以JSCORS徐州試點(diǎn)工程為例[J].現(xiàn)代測繪，2006，29(6):3～6

Research on the Modeling of Tropospheric Delay Based on a Region Atmosphere Features

Chen ZhongXin1，2，Zhu LiQiang1，2
(1.Suzhou Industrial Park Geone Information Technology CO.，Ltd.Suzhou，215021，China；2.JSCORS Suzhou Center，Suzhou 215021，China)

This paper bases on analysing the two familiar tropospheric delay model，using the SZCORS date to prove the relationship between altitude and latitude.Then It indicates the necessary to build the modeling of tropospheric delay based on a region atmosphere features.

tropospheric delay；region model；network RTK

1672－8262(2010)04－93－04

P228

A

2009—11—23

陳中新(1974—)，男，高級工程師，研究方向:為GNSS連續(xù)運(yùn)行參考站系統(tǒng)的建設(shè)、應(yīng)用及相關(guān)領(lǐng)域的研發(fā)。

現(xiàn)代工程測量國家測繪局重點(diǎn)實(shí)驗(yàn)室經(jīng)費(fèi)資助項(xiàng)目(TJES0907)