高速鐵路CP0基線解算中天頂對流層參數估計研究

任曉春，周東衛

(中鐵第一勘察設計院集團有限公司，西安 710043)

高速鐵路CP0基線解算中天頂對流層參數估計研究

任曉春，周東衛

(中鐵第一勘察設計院集團有限公司，西安 710043)

對流層延遲是GPS測量的重要誤差源之一，對模糊度解算及基線精度均有較大影響。高速鐵路CP0框架控制網基線解算中采用PWL分段線性法估計天頂對流層濕延遲參數的方法提高對流層折射改正精度，因此研究PWL分段線性法中參數估計的時間間隔對CP0基線解算的影響是十分必要的。介紹GPS對流層延遲的改正原理，通過設計不同解算方案采用工程測量數據對天頂對流層濕延遲參數估計的合理時間間隔進行研究。結果表明，每4～6 h估計一個天頂對流層濕延遲參數效果較好，能真實反映出對流層折射影響隨時間變化的趨勢，從而提高基線解的精度。

高速鐵路；GPS；對流層延遲；框架控制網；分段線性法；參數估計

1 概述

高速鐵路線路長、地區跨越幅度大且平面控制網沿線路呈帶狀布設，為了控制帶狀控制網的橫向擺動，并為平面控制測量提供統一的起算基準，實現勘察設計、施工建設和運營維護各階段控制網的“三網合一”，高速鐵路采用GPS精密定位測量技術，按一定間距布設建立了框架控制網(CP0)。CP0控制點布設間距為50～100 km，與國際IGS參考站或國家CGCS2000 A、B級GPS點進行聯測，采用GAMIT計算軟件進行中長基線解算，解算精度要求較高[1]。

對流層折射誤差不僅影響位置(特別是高程)精度，甚至會影響模糊度的解算，因此CP0長基線解算時必須顧及對流層折射誤差的改正精度，其改正一般通過選擇適宜的天頂對流層延遲模型及映射函數，并對天頂對流層濕延遲進行參數估計來實現的。考慮到天頂對流層濕延遲是引起對流層折射誤差的主要來源，且難于建立精確的模型[2]，GAMIT軟件中通過每隔一定時間間隔引入天頂對流層濕延遲參數進行估計的方法提高折射改正精度，即PWL分段線性法，因此，PWL中參數的估計時間間隔將對基線解算精度產生重要的影響。選取工程測量數據進行處理分析，研究了PWL分段線性法中參數估計時間間隔對CP0高精度數據處理的影響，得出了一些有益的結論。

2 對流層延遲改正原理

2.1 天頂對流層延遲及其影響

對流層延遲由干延遲與濕延遲兩部分組成。干延遲是由非水蒸氣部分的大氣延遲產生的，約占總延遲的90%，干延遲比較有規律，天頂方向可以1%的精度估算；濕延遲是由大氣中水蒸氣部分引起的，約占總延遲的10%，但濕延遲很復雜，影響因素較多，目前只能以10%～20%的精度估算[3]。雙差對流層延遲和天頂對流層延遲的關系可以用下面的公式來表示

?k?m

ωk

數據處理中可通過映射函數精確地估計出天頂對流層濕延遲(干延遲可通過模型計算得到)，然后求得信號傳播路徑上的對流層延遲，加以改正。在衛星高度角比較低時，對流層延遲表現出各向異性，即對流層延遲映射函數隨方位角的變化而不同，還需要估計對流層水平梯度的影響[6]。

研究表明[3]，在觀測誤差一定時，對流層延遲引起的基線誤差只與衛星與測站空間的幾何分布和測站的位置有關。當基線兩端的天頂對流層延遲差為2.25 m，GPS衛星的高度角為20°時，天頂對流層延遲引起的基線相對誤差約為1 ppm(1 ppm=10-6)，且基線越長，其影響越大。高速鐵路CP0控制點間距一般為50～100 km，由此引起的基線解算誤差為50～100 mm。對流層延遲引起基線的垂直分量誤差特別顯著，當高度角為20°時，每1 mm天頂對流層延遲差將引起基線垂直分量3 mm的誤差。由于測站溫度變化1 ℃或濕氣壓變化10-5hPa，將引起天頂延遲變化6～10 mm，結果可能引起(高度角為20°時)基線垂直分量18～30 mm的差異。因此，在高速鐵路CP0框架控制網基線解算過程中必須對天頂對流層延遲進行高精度改正。

2.2 天頂對流層延遲參數估計法

前已述及，對流層延遲模型可以大幅修正折射誤差，但由于受到模型誤差、氣象元素誤差以及實際大氣狀態和大氣模型間差異的影響，模型的改正精度是有限的。目前針對上述問題的解決方法主要是采用參數估計法。參數估計是指將利用各種經典模型求得的對流層延遲改正視為一種初始的近似值，然后將天頂對流層延遲的剩余誤差作為待定參數，與待定點坐標以及其他的待求量一起估計求解，通過平差計算來估計其精確值。常用的參數估計方法有隨機過程法和分段線性法兩種。

由于天頂方向的濕分量延遲服從一階高斯—馬爾科夫隨機過程[2]，因此隨機過程法采用隨機游走模型來描述天頂對流層濕延遲隨時間的變化規律，實際計算中可通過建立狀態方程和觀測方程一起構成卡爾曼濾波方程動態估計天頂對流層濕延遲分量[2]。盡管該方法是最理想的對流層延遲估計方法，但要求每個歷元都要有一個待估參數，估計參數太多，適用于動態實時解算，無法滿足高速鐵路CP0框架控制網高精度靜態解算的要求。在該法基礎上提出分段線性法，該方法將整個時段分為若干個子區間，每個區間各引入一個天頂對流層濕延遲參數進行估計。如在觀測時段內，每隔一定的歷元間隔k取一個狀態未知數，用步長k·Δt的離散隨機過程來表示對流層延遲隨時間的隨機變化。設相鄰兩個待估參數為X(i)和X(i+k)，則狀態方程為

分段線性法適用于時段長、天氣變化不太規則的場合，但引入的參數個數較多。GAMIT 軟件中通常采用PWL分段線性法，該方法將整個觀測時段分為若干個子區間，每個區間各引入一個參數，在每個區間內，各測站天頂方向的對流層濕延遲折射隨時間呈線性變化。假設某一子區間內i到i+k之間的歷元j滿足線性變化規律，則有

將X(i)代入相應的觀測方程，就可以得到以X(i)為未知數的法方程。同時

式中，k為某一子區間中的歷元個數；t為子區間時間間隔；I為數據采樣率。

通過選取適當大小的參數k，以獲得適當數量的天頂對流層濕延遲待估參數X(t)、X(i+k)，從而可以使用最小二乘方法估計出天頂對流層濕延遲分量參數。例如對于24 h單日的觀測數據，若每隔2 h估計1個天頂延遲因子，則PWL方法為每個測站生成13 個天頂對流層濕延遲參數，在平差過程中與其他未知參數一并進行解算。

3 天頂對流層濕延遲參數選取方法對比研究

3.1 試驗方案

選擇國內某高速鐵路在勘測設計階段建立的CP0框架控制網2007年2月8日(DOY39)和10日(DOY41)兩天的觀測數據進行計算分析。控制網共布設6個控制點KY01、SP02、CC03、DH04、LL05和HE06，點位間距為67～99 km，共聯測2個IGS參考站點BJFS和SUWN。

為了比較不同天頂對流層濕延遲參數選取方法對基線解算結果的影響，選用GAMIT 10.40基線解算軟件分析數據，在其他參數設置均相同時，改變參數估計時間間隔，設計如下6種試驗方案。

方案一：不進行天頂對流層濕延遲參數估計；

方案二：每24 h估計1個天頂對流層濕延遲參數，共估計2個參數；

方案三：每8 h估計1個天頂對流層濕延遲參數，共估計4個參數；

方案四：每6 h估計1個天頂對流層濕延遲參數，共估計5個參數；

方案五：每4 h估計1個天頂對流層濕延遲參數，共估計7個參數；

方案六：每2 h估計1個天頂對流層濕延遲參數，共估計13個參數。

GAMIT 10.40軟件中基線解算其他相關參數設置見表1。

表1 GAMIT參數設置

3.2 基線精度評定方法

(1)標準化均方根誤差NRMS

GAMIT解算結果中的標準化均方根誤差NRMS(Normalized Root Mean Square)可以用來衡量單時段解出的基線值偏離其加權平均值的程度，是從歷元的模糊度解算中得出的殘差，是衡量GAMIT解算結果的一個重要指標。根據國內外GPS數據處理經驗，其值一般應小于0.3，若NRMS值太大，則說明處理過程中周跳可能未得到完全修復。其計算公式如下[7]

(2)基線分量改正量

一般來說，基線分量的改正量不能大于其站點坐標約束量的2倍。否則，需要查看Q文件或autcln.sum文件是否有大量觀測數據被autcln模塊運行時所刪除。根據規范[8]要求，同時段觀測值的數據剔除率宜小于10%。

(3)基線重復率

基線重復率是衡量GPS基線解算結果質量的重要指標之一，可以衡量多時段基線解的解算質量。其計算公式如下[7]

進一步以單時段解的基線重復性為觀測值，用線性擬合方法求出基線重復率的常數部分a以及與邊長成比例的系數部分b

(4)基線網平差精度

基線網三維約束平差后基線向量的殘差和點位坐標精度也可作為評判基線解算結果好壞的一個參考指標。基線網三維約束平差通過固定IGS參考站在某一參考框架參考歷元下的三維地心坐標，對上述6種方案解算的基線結果分別平差，求出各條基線的殘差、各測站點的坐標及其精度指標。基線向量殘差和CP0點位中誤差越小，說明基線解算精度越高。

3.3 試驗結果分析

(1)NRMS

圖1是不同試驗方案DOY39和DOY41兩個單天解算結果的NRMS值，“none”代表不估計天頂對流層濕延遲參數。圖中NRMS值均小于0.3，說明基線解算質量合格。從圖中可以看出，不估計和估計天頂對流層濕延遲參數對于NRMS值的影響不大，且不同方案兩個單天解的解算結果沒有表現出明顯的規律性，因此僅僅從NRMS值來判定何種方案更好是遠遠不夠的。

圖1 不同方案解算結果的NRMS值

(2)基線分量改正量

圖2是不同試驗方案DOY39和DOY41兩個單天解算結果的XYZ基線分量最大改正量。計算方案中對IGS參考站點BJFS和SUWN進行強約束，設置值分別為0.03、0.03、0.05 m，其余點為非約束點。圖中基線分量最大改正量為0.035 m，均小于其站點坐標約束量的2倍，說明基線解算質量合格。從圖中可以看出，不估計和估計天頂對流層濕延遲參數對基線分量改正量的影響不大，不同方案兩個單天解的解算結果沒有表現出明顯的規律性，因此也無法從基線分量改正量來判定何種方案更好。

圖2 不同方案基線分量最大改正量

(3)基線重復率

表2為不同試驗方案解算結果站心地平坐標系下的基線向量和邊長的重復率。基線重復率是根據基線解算結果采用公式(7)擬合得到的，由常數部分和系數部分組成，常數部分單位為 mm，系數部分基線長度L的單位為m。

表2 不同方案解算結果的基線重復率

為了更直觀地反映基線重復率誤差的大小，對50 km長度基線的重復率所引起的絕對誤差進行了計算，計算結果見圖3。

圖3 不同方案50 km長度基線重復率引起的絕對誤差

從表2和圖3可以看出。

①天頂對流層濕延遲參數估計個數的不同，對基線N、E分量的影響約為2 mm，對U分量的影響約為9 mm，對基線解算產生的誤差主要體現在U分量上。

②當對流層處于相對平靜狀態時，估計和不估計天頂對流層濕延遲基線重復率變化不顯著。

③估計不同個數的天頂對流層濕延遲參數基線解算的效果比較接近，但是仍然有區別，估計7個天頂對流層濕延遲參數的L基線重復率是最小的，其余方案的L基線重復率基本相當。

(4)基線網平差精度

基線網三維約束平差固定IGS參考站BJFS(參考框架ITRF2000，參考歷元1997.0)的三維地心坐標，對上述6種方案解算的基線結果分別平差，求出各測站點的三維坐標及精度指標，IGS參考站SUWN與其相應框架相應歷元下的已知坐標XYZ之差均不超過10 cm。三維約束平差計算得到各基線三維向量的殘差和點位絕對精度分別如圖4、圖5所示。

圖4 不同方案約束平差基線向量殘差

圖5 不同方案約束平差點位絕對精度

從圖4、圖5可以看出，單天解算時估計7個天頂對流層濕延遲參數時基線向量殘差最小、點位坐標精度最好，估計5個參數的效果次之，估計13個參數的基線向量殘差最大，不估計天頂對流層濕延遲參數的點位坐標精度最差。事實上，估計5～7個天頂對流層濕延遲參數，其結果不論是基線重復率還是平差精度指標都是非常接近的。這是因為天頂對流層濕延遲參數估計個數太少，不足以描述觀測時段內測站上空的中性大氣折射變化情況，特別是存在較大濕度導致對流層變化劇烈的情況；參數估計的太多又增加了平差解算的參數個數，降低了方程解算精度。因此，對于CP0框架控制網來說，每4～6 h估計1個天頂對流層濕延遲參數時基線解算效果較好。

4 結語

本文選用工程測量數據對高速鐵路CP0框架控制網數據處理過程中天頂對流層濕延遲參數估計及選取方法進行了研究，研究結論如下。

(1)天頂對流層延遲對基線解算結果的影響與基線長度有關，基線越長，影響越大，CP0數據處理時必須予以精確改正。

(2)天頂對流層濕延遲參數估計個數的不同，對50 km長度基線N、E分量的影響約為2 mm，對U分量的影響約為9 mm，對基線解算產生的誤差主要體現在U分量上；當對流層處于相對平靜狀態時，估計和不估計天頂對流層濕延遲基線重復率變化不顯著。

(3)天頂對流層濕延遲對基線網平差精度影響較大，單天解估計7個天頂對流層濕延遲參數時基線向量殘差最小、點位坐標精度最好。事實上，估計5～7個天頂對流層濕延遲參數，其結果不論是基線重復率還是網平差精度指標都是非常接近的。因此，在使用GAMIT軟件進行CP0基線解算時，將天頂對流層濕延遲參數的估計間隔設置為4～6 h效果較好，能真實反映出對流層折射影響隨時間變化的趨勢，從而提高基線解的精度。

[1] 中華人民共和國鐵道部.TB 10601—2009 高速鐵路工程測量規范[S].北京:中國鐵道出版社， 2009．

[2] 李征航，黃勁松.GPS測量與數據處理[M ].2版.武漢:武漢大學出版社，2011．

[3] 歐吉坤.GPS測量的中性大氣折射改正的研究[J].測繪學報，1998，27(1):31-36．

[4] Niell A. E. Global mapping functions for the atmosphere delay at radio Wavelengths[J]. Journal of Geophysical Research， 1996，101(B2):3227-3246．

[5] Saastamoinen J. Contribution to the Theory of Atmospheric Refraction[J]. Bulletin Geodesique， 1973，107:13-34．

[6] Yoaz E， Peter M. Estimating horizontal gradients of tropospheric path delay with a single GPS receiver[J]. Journal of Geophysical Research， 1998，103(B3):5019-5035．

[7] 李征航，張小紅.衛星導航定位新技術及高精度數據處理方法[M].武漢:武漢大學出版社，2009．

[8] GB/T18314—2001 全球定位系統(GPS)測量規范[S].北京:中國標準出版社，2001．

[9] 周東衛.GNSS參考站網絡的對流層完備性監測技術研究[J].工程勘察，2012，40(10):65-70．

[10] 葛茂榮，劉經南.GPS定位中對流層折射估計研究[J].測繪學報，1996，25(4):255-291．

[11] 周樂韜.連續運行參考站網絡實時動態定位理論、算法和系統實現[D].成都:西南交通大學，2007．

[12] MIT. Department of Earth， Atmospheric， and Planetary Sciences. GPS Analysis at MIT (Release 10.4)， 2010．

[13] 余興勝.構建橋涵勘測技術新體系的研究[J].鐵道標準設計，2013(2):73-76.

[14] 田林亞，祖力比亞·阿布都熱西提等.高速鐵路測量中高斯平面坐標與斜軸墨卡托平面坐標的轉換[J].鐵道標準設計，2013(4):13-16.

Research on Estimation of Parameters of Zenith Troposphere with CP0 Baseline Resolution for High-speed Railway

Ren Xiaochun， Zhou Dongwei

(China Railway First Surveying and Design Institute Group Ltd.， Xi’an 710043, China)

The refraction of the tropospheric delay on radio signals， which is one of the important error sources in GPS surveying， has great influences on the accuracy of ambiguity and baseline resolution. The method to estimate zenith tropospheric wet delay parameters of PWL model， which may be applied to CP0 baseline resolution of high-speed railway， improves the accuracy of the tropospheric delay correction effectively. Thus， it is necessary to conduct a research on the time interval selection of zenith tropospheric wet delay parameters of PWL model. In this paper， the principle of GPS tropospheric delay correction is introduced， and then surveying data in different resolution schemes are analyzed to determine a reasonable time interval of zenith tropospheric wet delay parameters of PWL model. The results show that the estimation of zenith tropospheric wet delay parameter every 4 to 6 hours has better baseline resolution effect， which truly indicates the trend of changes of the tropospheric refraction effect with time， and， thereby， improves the accuracy of baseline resolution．

High-speed railway; GPS; Tropospheric delay; Frame control network; CP0; PWL; parameter estimation

2014-01-15；

：2014-02-12

任曉春(1962—)，男，教授級高級工程師，1982年畢業于西南交通大學鐵道航空勘察專業，工學碩士。

1004-2954(2014)11-0036-05

U238； P228

：A

10.13238/j.issn.1004-2954.2014.11.009