帶狀廣域高程擬合方法的探討

李明軍，楊國東

(1.鐵道第三勘察設計院集團有限公司航遙測繪分院，天津 300251； 2.吉林大學地球探測科學與技術學院，吉林長春 130026)

帶狀廣域高程擬合方法的探討

李明軍1?，楊國東2

(1.鐵道第三勘察設計院集團有限公司航遙測繪分院，天津 300251； 2.吉林大學地球探測科學與技術學院，吉林長春 130026)

結合高速鐵路建設中C級GPS控制網和二等水準數據，對石鄭客運專線沿線143 km帶狀區域，分別構建多面函數模型、二次多項式和EGM2008重力場模型進行高程擬合分析。在特定區域，確定帶狀廣域高程擬合最佳模型，驗證獲取厘米級甚至更高精度的正常高，以滿足鐵路工程、公路工程以及其他帶狀工程施工要求的可行性。

高程擬合模型；正常高；大地高；高程異常

1 引 言

隨著科學技術的不斷發展，GPS定位技術已廣泛用于各種工程平面控制測量，其平面控制精度已可達到10－9數量級上[1]。但由于參考橢球體與似大地水準面之間的關系難于準確確定，致使GPS大地高同水準高轉換模型中高程異常無法用帶有固定參數的數學模型精確求解。造成了GPS高程在實際工程應用中，利用率相對比較低，甚至出現不被利用的現象。為此，本文以石鄭客運專線為例，采用數學模型構建C級GPS控制網與二等水準數據之間的關系，分幾個方案研究帶狀區域高程擬合方法，以獲取高精度高程異常，減少外業水準工作量，從而提高大地高的利用率，進而提出鐵路工程、公路工程以及其他廣域帶狀工程最佳高程擬合方式。

2 高程擬合原理

高程擬合是指根據測區內若干個既進行了GPS測量又聯測了水準高程的控制點，利用大地高、水準高和高程異常之間關系:

(式中，H為大地高，H正常為地面點沿正常重力垂線到似大地水準面的距離－正常高，也就是這里說的水準高。)推算出各點上的高程異常值，利用這些離散點上的高程異常值和GPS點的平面坐標之間存在的關系，構建測區內高程擬合隨坐標不斷變化的數學模型，從而可以根據控制點坐標計算出未知點上的高程異常，實現橢球高向正常高的轉換。高程異常值和GPS點的平面坐標之間存在的關系可由以下方程式表示:

ξ(x，y)為GPS測點的高程異常值，x，y為測點坐標，α為待定參數。f(x，y，α)為所選的函數模型。選擇一定的函數模型，利用已知量和推求出待定參數，可得到該測區似大地水準面的一個逼近，從而可以內插出其余GPS點的正常高值[2]。

3 高程異常解算模型

3.1 多項式模型

似大地水準面的擬合也可采用一次或多次曲面擬合法，即對于公共點上的高程異常與平面坐標之間存在以下數學關系:

一般根據聯測水準的GPS點的個數來選取多項式的模型，比如當聯測水準的GPS點的個數n＝3，4，5時，可選用一次多項式，6以上可以選用二次或更高次多項式。但一般需要擬合的大地水準面區域較小，不會進行大量水準聯測，再加之多項式擬合自身系統精度的限制，一般只應用一、二和三次多項式擬合。

3.2 EGM2008重力場模型

地球重力場模型是通過最新衛星跟蹤數據、地面重力異常數據及衛星測高等重力場信息，建立的重力位的球諧函數級數展開系數，簡稱位系數(Cnm，Snm)。利用位系數，根據數學模型:

可解算地面任意一點P的高程異常。式中:ρ，φ，λ分別為地面點P的向徑、地心緯度和精度，GM為引力常數與地球質量的乘積，γP為地面點P的正常重力值，a為參考橢球的長半軸，在實際計算中采取相對于GRS80的值，即a＝6 378 137。為完全規格化位系數，其中偶次帶諧系數代表實際引力位與正常引力位之差，(sinφ)為完全規格化Legendre函數，N則為地球展開的最高階數。

EGM2008是NGA(US National Geospatial Intelligence Agency)釋放的全球超高階地球重力場模型。該模型階次完全至2 159(另外球諧系數的階N擴至2 190，次為2 159)，模型分辨率相當于5′(約9 km)[3]。

3.3 多面函數模型[4]

多面函數的主要思想則是根據“任何數學表面和任何不規則的圓滑表面，總可以利用一系列的有規則的數學表面的總和，以任意精度逼近”原理，建立如下函數模型:

其中xi和yi為中心點，Q(x，y，xi，yi)為核函數，為了簡單一般采用:

或者:

本文則采用式(7)，多面函數圓滑因子取δ＝0.001。

4 實際案例分析

本文針對帶狀線路工程的特征，實施GPS水準點高程擬合均勻性最佳優化設計，對石鄭客運專線沿線143 km，嚴格執行《客運專線無碴軌道鐵路工程測量暫行規定》(鐵建設[2006]189號)，獲取的100個GPS水準點成果，根據上述各種模型的特點，分方案進行了數據處理和分析，其具體案例如圖1所示。

圖1 擬合區域高程異常概算擬合曲面線狀圖

4.1 多面函數擬合案例分析

針對多面函數模型，對現有GPS水準點分別以平均4.5 km、6 km、7.5 km、11.5 km和15 km的主體間距進行數據提取，并依次進行高程擬合，其單點高程誤差統計表及相鄰點高差誤差統計圖如表1和圖2所示。

多面函數擬合單點高程誤差統計表 表1

圖2 多面函數擬合相鄰點高差誤差統計圖

根據四等水準往返高差閉合差限差公式，mh限＝20L，L為路線長度。

根據多面函數擬合成果可以看出，其單點求解高程精度很高，在現有數據最優一組成果中，其最大高程偏離值僅為2.3 cm，而其高差主體滿足四等水準要求。不足的是:在4.5 km處，只有一個略微超限，超限2 mm；11.5 km和15 km處已知GPS水準點包含范圍內，有6個超限，其主要原因是由于離它們的最近兩個GPS水準已知點概算高程異常值出現突變所致；6 km～11.5 km，已知GPS水準點縱向包含范圍外各有幾個超限，并且隨著縱向距離越遠其超限值越大；在間距為15 km時，考慮了縱向范圍包含性問題，在縱向最外延加入已知GPS已知點，上述因外延引起的超限問題得到了解決，并取得了比較理想的成果。

4.2 多項式擬合案例分析

本次多項式擬合實驗中，分別選用了相鄰的24個和63個GPS水準點作為基準數據，然后均按平均4.5 km的間距進行了二次多項式高程擬合。24點擬合相鄰點高差誤差主體滿足四等水準要求，只有一個點超限，超出限差6 mm；63個點擬合相鄰點高差共有5個點超限(其中一個點也是因為外延點而引起的超限)，其他4個中最大超限值為24 mm。

4.3 EGM2008高程擬合案例分析

同樣，本次實驗選擇98個GPS水準點作為基準數據，分別以4.5 km、6 km、7.5 km和15 km的主體間距提取GPS水準數據作為已知數據，然后依次實施EGM2008平面高程擬合，全線單點高程誤差和擬合相鄰點高差誤差統計如表2、圖3所示。

從以下結果可以看出，EGM2008重力高程擬合成果整體比較均一，只要已知GPS水準點布點均勻，其高程誤差以及相鄰點高差誤差基本不受已知點個數限制，其相鄰點高差誤差最大只為0.23 cm，精度均滿足四等水準要求，只是其單點高程擬合精度只是局限于10.5 cm，有待進一步提升。

EGM2008平面高程擬合單點高程誤差統計表 表2

圖3 EGM2008平面高程擬合相鄰點高差誤差統計圖

5 結 論

綜合考慮帶狀廣域線性工程特點，對石鄭線現有GPS水準數據，基于三種擬合方法，分方案實施143 km帶狀區域的高程擬合，通過擬合成果及分析，可以得到如下結論:

(1)三種擬合模型均能獲取高精度高程擬合成果，其擬合精度均能達到四等水準要求，甚至更高。對于廣域帶狀區域而言，采用多面函數模型能獲取高準確度點位高程擬合成果，且其已知GPS點高程擬合殘差均為0，但要獲取高精度擬合成果，它對已知點點位密度和均勻性要求比較高，尤其是在當高程異常變化比較大的情況下，更應想辦法加密控制點的密度，且當未知點超出帶狀控制點縱向范圍，其擬合高程成果不可靠；多項式擬合同樣能獲取高準確度的點位高程成果，但相對于多面函數模型，其對廣域帶狀區域的高程擬合不能建立統一的數學模型來實施對全線其他GPS點正常高的求解，計算過程比較復雜；相對于上述兩種高程擬合模型，EGM2008重力場模型對點位密度和均勻性要求相對較低，能獲取更高準確度的相鄰點高差。但其單點擬合高程準確度略顯偏低，高程誤差有待更進一步提升。因此為獲取高質量廣域帶狀區域高程擬合成果，更適合采用多面函數和EGM2008模型。而在使用多面函數模型時，綜合考慮工程進度和點位高程精度，在高程異常變化比較大的區域要獲取高精度擬合成果，高程已知點點位密度應以5 km/個～8 km/個為宜；

(2)在實際帶狀工程生產應用中，高程擬合已知點個數總是有限，多面函數高程擬合模型的準確度也有待更進一步的提高。通過本文案例分析，是否可以結合多面函數高程擬合成果與EGM2008模型獲取的相鄰點高差數據，在已知水準點的基礎上進行高程平差后的成果，來獲取高準確度的高程成果，可作為高程擬合模型研究的下一課題。

[1] 張冰，馬開鋒，陳南祥.GPS高程擬合模型研究[J].華北水利水電學院學報；2008(2):65～67

[2] 徐紹銓，張華海，楊志強等.GPS原理及應用[M].武漢:武漢大學出版社，2003

[3] 金時華.多面函數擬合法轉換高程擬合模型[J].測繪與空間地理信息，2005年，(6):44～47

[4] 劉成，張幸福.EGM2008重力場模型在GPS高程擬合中的應用分析[J].鐵道勘察；2009(1):1～3

[5] 匡翠林.高精度GPS水準算法研究及其應用[D].長沙:中南大學，2004

[6] 王凱紅，張俊影，趙娟等.GPS高程擬合方法精度分析[J].吉林地質，2007(3):84～87

Discussion on Elevation Fitting Method of Band Wide－area

Li MingJun1，Yang GuoDong2
(1.Department of Remote Sensing and Surveying of theThird Surveying and Design Institute Co.，Ltd.Tianjin 300251，China； 2.College of Geoexploration Science and Technology of JiLin university，ChangChun 130026，China)

In this paper，based on C－class GPS control network and second level data in the high－speed railway construction，researched on 143km band wide－area along passenger dedicated high－speed railway of which from ShiJiazhuang to WuHan，By using of multi－faceted function model and quadratic polynomial model and EGM2008 gravity field model，carried out and analysis elevation fitting。In a specific region determine the best band wide－area fitting elevation model，and verification the feasibility to obtain the normal high posses centimeter－level or even higher accuracy，to meet the railway，road and other ribbon construction requirements。results of study has guiding significance on heoretical research and practical work.

elevation fitting model；normal height；geodetic height；abnormal elevation

1672－8262(2011)02－96－04

P228

B

2010—07—22

李明軍(1984—)，男，助理工程師，主要從事GPS與3S應用研究。