(中鐵第四勘察設計院集團有限公司，湖北武漢　430070)

鐵路勘測過程中，線路沿線的高程數(shù)據(jù)是確定線路縱斷面的基礎數(shù)據(jù)。這些高程數(shù)據(jù)要么是采用幾何水準測量方法測量的正常高，要么是采用GPS高程擬合方式擬合出的正常高。傳統(tǒng)的幾何水準測量雖然精度相對較高，但作業(yè)效率低，尤其在相對高差大的山區(qū)及荒漠無人區(qū)，不能很好地滿足大量高程數(shù)據(jù)快速采集的要求。近年來，GPS技術以其高精度、全天候、低成本、高效率的優(yōu)點被廣泛應用于鐵路勘測的各個環(huán)節(jié)。GPS測量的參考橢球為WGS84橢球，GPS直接測得的高程是相對于WGS84橢球面的大地高，與勘測過程中采用的正常高系統(tǒng)不一致，故將GPS測得的大地高用GPS高程擬合的方法轉(zhuǎn)換為工程中需要的正常高是一項必要工作。常規(guī)的GPS高程擬合方式需要聯(lián)測較為密集的高程點，且擬合精度較低。

1　EGM2008地球重力場模型

地球重力場模型是采用一類基本參數(shù)的集合，用逼近或擬合的方式來描述和表示地球重力場的模型。2008年4月，NGA(美國國家地理空間情報局National Geospatial-Intelligence Agency)發(fā)布了最新的全球重力場模型，命名為EGM2008[1]。EGM2008研究和總結(jié)了以往地球重力場模型理論和經(jīng)驗，采用了各種先進的基礎數(shù)據(jù)，其中包括：

①以PGM2007為基礎演變出來的PGM2007B為參考模型；

②利用了帶有ITG-GRACE03S位系數(shù)信息以及相應協(xié)方差信息的GRACE衛(wèi)星衛(wèi)星跟蹤數(shù)據(jù)；

③利用了衛(wèi)星測高數(shù)據(jù)(TOPEX)和全球5′×5′空間分辨率的地面重力異常數(shù)據(jù)；

④利用了高精度、廣域覆蓋的地面重力數(shù)據(jù)。

EGM2008重力場模型的球諧系數(shù)階次完全至2 159階，擴展到2 190階，相當于約9 km的模型空間分辨率，是目前最高分辨率和精度的重力場模型。EGM2008重力場模型提供了包括2 190階次的全球重力場模型、全球5′×5′網(wǎng)格重力異常、全球5′×5′、2.5′×2.5′網(wǎng)格大地水準面、全球5′×5′網(wǎng)格垂線偏差(ξ，η)的最終成果。

針對我國情況，國內(nèi)學者進行了深入的研究：EGM2008重力場模型提供的高程異常在我國大陸地區(qū)的總體精度為20 cm，中東部地區(qū)12 cm，北方地區(qū)達到9 cm，西部地區(qū)為24 cm[2]。應用EGM2008地球重力模型進行的大地水準面計算、重力位計算、不同尺度重力場分析計算等將更為精確。

2　GPS高程擬合

如圖1所示，H為大地高，Hg為正高，Hr為正常高，hg為大地水準面到地球橢球面的距離，ξ為似大地水準面到地球橢球面的高差，叫做高程異常。

圖1　高程異常

GPS測量的參考橢球體位WGS-84橢球，測得的大地高H是基于WGS-84橢球表面的大地高H84，我國鐵路勘測中采用的國家1985高程基準的85高程H85是基于似大地水準面的正常高Hr，這種WGS-84橢球面與似大地水準面之間的差距叫高程異常，可表示為

ξ=H-Hr=H84-H85

通常采用GPS高程擬合的方法獲取這些正常高，因此計算出精度可靠的高程異常十分重要。GPS高程轉(zhuǎn)換的基本原理是通過測得GPS點的大地高和已知該GPS點的正常高獲得該點的對應高程異常值。由于高程異常在小區(qū)域范圍內(nèi)的變化相對平緩，可以用一些擬合的方法對區(qū)域內(nèi)的高程異常值進行估計，然后利用已知點正常高計算出擬合系數(shù)，求出待定點的高程異常值，進而獲得該點的正常高。

3　基于EGM2008模型的GPS高程擬合方法

在以地球質(zhì)心為坐標原點、Z軸與地球自轉(zhuǎn)軸相重合坐標系中的球坐標中，以r表示向徑，θ表示緯度，λ表示經(jīng)度，由Bruns公式可以推導出地球表面上任意一點P的高程異常

根據(jù)物理大地測量學的理論，高程異?？杀硎緸?/p>

ξ=ξGM+ξΔG+ξT

其中ξT表示短波部分，可以通過求解地形改正得到；ξΔG為中波部分，可以通過求解重力異常的邊值得到；ξGM為長波項，可由重力場模型計算出[4]。采用基于EGM2008地球重力場模型中的重力場中長波信息可以求解ξGM項，而其中的中波、短波部分可通過GPS測量中若干已知正常高的GPS點逼近，擬合出其他點ξΔG、ξT兩項值。

基于EGM2008模型的GPS高程擬合方法在實際工程應用中具體步驟可概括為：

①求取GPS點高精度的WGS84平面坐標及大地高H84。在實際工程中，可以通過建立GPS網(wǎng)，并與IGS站聯(lián)測來實現(xiàn)。

②利用EGM2008地球重力場模型求解由重力場模型所計算的各GPS點高程異常的長波項ξGM。

③利用GPS網(wǎng)中部分已知正常高的GPS點，通過GPS高程擬合的方法擬合出未知正常高的GPS點正常高Hr。

4　算例分析

某新建鐵路工程位于天山南麓，線路長度約130.5 km。線路起點、終點位于兩條國道附近，能收集到已知正常高的GPS點；其間線路穿越約80 km的無人區(qū)，無控制點資料。測區(qū)最低處海拔約1 600 m，最高處海拔約3 700 m，平均海拔約2 500 m。勘察初期，沿線路方向建立了CPI控制網(wǎng)[5](帶狀的四等GPS網(wǎng))，全網(wǎng)共計40個CPI點；并在線路約0 km、12 km、100 km及130 km處聯(lián)測了已知正常高的GPS點；勘察中期，采用16臺Dini12電子水準儀，按三等水準作業(yè)要求對全部CPI點進行了水準測量；工程控制網(wǎng)及聯(lián)測情況如圖2所示。通過傳統(tǒng)的GPS高程擬合與基于EGM2008地球重力場模型的GPS高程擬合實驗，根據(jù)擬合的高程結(jié)果和幾何水準成果進行統(tǒng)計分析，得出了一些結(jié)論。

圖2　工程控制網(wǎng)及聯(lián)測情況示意

4.1　傳統(tǒng)GPS高程擬合

先按照傳統(tǒng)的GPS高程擬合的方法進行高程擬合。由于4個已知正常高的GPS點位于GPS網(wǎng)兩端，網(wǎng)段中間大于80 km無約束點。采用首末兩個已知正常高的GPS點約束擬合時，離起算點遠的CPI點擬合高程與幾何水準高程之差最大達1 137 mm；用四個已知正常高的GPS點約束擬合時，離起算點遠的CPI點擬合高程與幾何水準高程之差最大達1 654 mm。傳統(tǒng)的GPS高程擬合不能滿足測量要求，且離起算點越遠差值越大，約束點間距越遠差值越大(如圖3所示)。

圖3　傳統(tǒng)GPS高程擬合結(jié)果與幾何水準高程比較

4.2　基于EGM2008地球重力場模型的GPS高程擬合

采用基于EGM2008地球重力場模型的GPS高程擬合時，先根據(jù)GPS網(wǎng)測得的高精度WGS84平面坐標及高精度的大地高H84，解算出高精度的WGS84三維坐標，故能得到高精度的WGS84橢球大地高。利用EGM2008地球重力場模型求解由重力場模型所計算的長波項ξGM(如表1所示)。

表1　各CPI點高程異常長波項　m

圖4　基于EGM2008模型的GPS擬合高程與幾何水準高程比較

5　結(jié)束語

利用基于EGM2008模型的GPS高程擬合方法進行高程擬合可以明顯提高高程擬合的精度，擬合精度可以達到四等精度要求。

利用基于EGM2008模型進行GPS高程擬合，已知水準點的間距可以放寬至80 km左右。不同地區(qū)的已知水準點間距有待進一步研究。

采用基于EGM2008模型進行GPS高程擬合時，需要采集較高精度擬合點的WGS84橢球大地高。

[1]Pavlis N K，S A Holmes，S C Kenyon，et all.An Earth GravitationalModeltoDegree 2160: EGM2008[R]. Presented atthe 2008 General Assembly of the European Geosciences Union， Vienna， Austria， April， 2008:13-14

[2]章傳銀，郭舂喜，陳俊勇，等.EGM2008地球重力場模型在中國大陸適用性分析[J].測繪學報，2009，38(4)：283-289

[3]魏子卿，王剛.用地球位模型和GPS/水準數(shù)據(jù)確定我國大陸似大地水準面[J].測繪學報，2003，32(1)：1-5

[4]陳俊勇.高精度局域大地水準面對布測GPS水準和重力的要求[J].測繪學報，2001，30(3)：189-192

[5]TB10101—2009鐵路工程測量規(guī)范[S]

[6]胡明城.現(xiàn)代大地測量學的理論及其應用[M].北京：測繪出版社，2003

[7]武鵬.GPS擬合高程應用于鐵路定測的研究[J].鐵道勘察，2009，35(6)：11-13

[8]孫風華，吳曉平，張傳定.我國1′×1′高程異常模型的建立及應用[J].地球物理進展，2005，20(4)：986-990

[9]夏哲仁，石磐，李迎春.高分辯率區(qū)域重力場模型DQM2000[J].武漢大學學報：科技信息版，2003，28：124-128

[10]左虎，范東明.GPS高程轉(zhuǎn)換的新方法[J].鐵道勘察，2006，32(6)：6-8