基于EGM2008模型進行蕪湖市似大地水準面精化

李華平，馬新瑩，章緯之，馬祥，謝方

(1.蕪湖市勘察測繪設計研究院有限責任公司，安徽蕪湖 241000； 2.國家測繪局大地測量數據處理中心，陜西西安 710054)

基于EGM2008模型進行蕪湖市似大地水準面精化

李華平1?，馬新瑩2，章緯之1，馬祥1，謝方2

利用高精度的GPS水準數據、重力數據、數字地形模型和EGM 2008地球重力場模型，采用移去-恢復技術進

行蕪湖市1 400 km2似大地水準面精化。分別采用蕪湖市坐標轉換重合點，利用二維高斯平面坐標轉換模型，建立了

蕪湖市2000國家大地坐標系與1980西安坐標系、1954年北京坐標系及蕪湖獨立坐標系的坐標轉換模型，并且新建了

蕪湖2000獨立坐標系。綜合利用多種現代擬合方法完成了分辨率為2.5′×2.5′，高精度、高分辨率的蕪湖市似大地水

準面模型的計算，模型擬合精度達到±1.1 cm，外部檢核精度達到±1.3 cm。

似大地水準面精化；GPS；水準測量；重力

1 引 言

近年來，在信息化浪潮的推動下，我國一些城市已經開始或正在啟動各自的數字城市建設計劃，城市信息化熱潮正在全國范圍內展開。數字城市建設是堅持科學發展觀的重要體現，將有力地推動我國城市信息化在深度和廣度上不斷發展深化，進而從全局上促進我國社會和經濟的可持續發展。

本項目的目的是在蕪湖市及周邊地區的國家連續運行站及蕪湖市已建成的CORS站基礎上，綜合利用地面重力數據、GPS、水準數據、數字地面高程模型數據，基于EGM 2008重力場模型等技術手段，確定蕪湖市區域高精度似大地水準面模型，并應用于蕪湖市的城市規劃、國土測量、大地測量、工程測量、氣象監測、地震監測、地面沉降監測以及城市地理信息系統等領域，同時兼顧社會公共定位服務，以滿足日益增長的城市綜合管理與城市化建設的需求。它是衛星導航定位技術、現代大地測量學與計算機技術的有機結合，可為促進蕪湖市經濟持續、健康、快速發展以及城市規劃、建設、管理與服務水平的提高提供支撐和保障。

EGM2008是近來由NGA釋放的全球超高階地球重力場模型，該模型的階次完全至2 159(另外球諧系數的階擴展至2 190，次為2 159)，相當于模型的空間分辨率約為5′(約9 km)。該模型采用了GRACE衛星跟蹤數據(ITG-GRACE03S位系數信息以及相應的協方差信息)、衛星測高數據和地面重力數據等，該模型無論在精度還是在分辨率方面均取得了巨大進步，采用該模型以及GPS/水準數據有望獲得更高精度的區域似大地水準面。

項目目標是依托國家連續運行參考站與蕪湖市已建立的4個連續運行參考站，利用GPS測量、精密水準測量、及(似)大地水準面精化技術，建立和維持區域大地基準、高程基準及高精度(似)大地水準面，項目包括:

(1)建立和完善蕪湖市GPS水準網；

(2)確定集平面控制和高程控制于一體的蕪湖市三維空間基準網及其應用系統；

(3)確定蕪湖市高精度(似)大地水準面模型。

2 項目利用資料情況

2.1 國家大地控制網資料

在蕪湖市精化區域及周邊有一條國家一等水準路線及多條“華東、華中區域大地水準面精化項目”安徽省二等水準路線通過。蕪湖市周邊有國家二等三角點14個，2000國家GPS大地控制網點1個，安徽省GPS A、B級點2個，安徽省GPSC級網點16個。

(1)蕪湖市精化區域城市控制網資料

①衛星連續運行站:在蕪湖市區域內已建成4個衛星連續運行站；

②GPS點:蕪湖精化區域及周邊已建成一期GPS三等點38個，二期GPS三等點36個，一期GPS四等點40個，二期GPS四等點31個；航測區域GPS四等點50個；

③水準點:蕪湖精化區域及周邊已建成二等水準網點51個，三等水準網點41個。

(2)DEM資料

全國已有1∶5萬DEM數據，精度較好，用于蕪湖市似大地水準面確立中修正格網地形改正和均衡改正。

(3)地形改正資料

目前項目區域內已有30″×30″格網地形改正和均衡改正成果，在項目實施過程中利用1∶5萬DEM完成3″×3″、30″×30″、2.5′×2.5′格網地形改正和均衡改正成果。

(4)加密重力資料

本項目可收集蕪湖市用于似大地水準面精化的該區域68 308個重力點。

3 GPS網及水準點的布設

3.1 GPS網

蕪湖市已有CORS站4個，在精化區域范圍內均勻布設GPS水準網點55個，其中利用國家及蕪湖市已有GPS控制點42個、新埋GPS點13個。為了建立蕪湖市已有坐標系統之間的相互轉換關系，同時在蕪湖市保存完好的6個國家二等三角點上實施了GPS觀測。本項目GPS網點共計61個。

GPS網采用與國家連續運行站及蕪湖市CORS基準站同步觀測模式，蕪湖市CORS站及GPS點取得的觀測時間需大于23.5 h，儀器6臺GPS觀測全部采用徠卡雙頻GPS接收機。

①觀測時段數:≥1個時段

②衛星截止高度角:≥10°

③同時觀測有效衛星數:≥4

④有效觀測衛星總數:≥9

⑤數據采樣間隔:30 s

⑥PDOP值:≤6

3.2 GPS水準點

為了滿足精化蕪湖市大地水準面的需要，根據蕪湖市已有國家一等水準網、“華東、華中區域大地水準面精化項目”安徽省二等水準網及蕪湖市已有的GPS及水準網分布情況，及點位普查情況，均勻選取GPS水準點作為精化蕪湖市似大地水準面的控制點，共計55個。其中約40個GPS水準網點用于蕪湖市似大地水準面精化數據處理，15個點用于最終蕪湖市似大地水準面外部檢核點。精化點和檢核點需均勻布設于蕪湖精化范圍區域內，并且需要同步進行組織施測，以保證蕪湖市似大地水準面的一致性。

具體水準線路如圖1所示。

圖1 聯測GPS水準點圖

水準觀測所使用測量儀器設備經儀器檢定部門檢定合格，并在有效期內。采用Ni002(Ni002A)自動安平光學水準儀或DiNi12、LeicaDNA自動安平數字水準儀等DSZ05、DS05類別儀器進行觀測。標尺采用線條式因瓦標尺或條碼式因瓦標尺，尺墊不小于5.0 kg鑄鐵尺墊。

4 數據處理

4.1 數據處理

(1)軟件

采用目前國際上公認的用于GPS數據后處理最優秀、最成熟的軟件美國麻省理工學院(MIT)的GAMIT/GLOBK軟件。

(2)數據處理工作內容

數據處理工作內容包括:蕪湖市及周邊地形與重力資料收集與整理、蕪湖市周邊地區我國連續運行基準站觀測數據收集，2000國家大地坐標系下4個蕪湖市CORS站及67個GPS點的數據處理；2000國家大地坐標系向1980西安坐標系、1954年北京坐標系及蕪湖獨立坐標系進行坐標轉換，建立蕪湖2000獨立坐標系；蕪湖市二等水準網數據處理；蕪湖市似大地水準面精化。

(3)已有資料的收集

項目收集了蕪湖市周圍泰安(TAIN)、上海(SHAO)、武漢(WUHN)、廈門(XIAM)四個國家GPS連續運行站的數據；蕪湖市二等水準網連測的4個國家一等水準點1985國家高程成果(1999年二期復測成果)；蕪湖市似大地水準面精化收集了國家測繪局館藏蕪湖市及周邊地區加密重力點68 308點；蕪湖市及周邊地區的數字地形模型資料收集了國家1∶5萬數字高程模型DEM數據；收集美國最新研制的高精度、高階次的EGM2008模型(2 160階次)。

(4)重力與地形資料的整理與處理

加密重力成果經過分析、轉換、粗差剔除后用于似大地水準面精化計算。坐標系統的統一采用坐標轉換的方法，將所有加密重力的坐標轉換為CGCS2000坐標系。下面主要介紹重力基準的統一。所有重力點的重力基準全部統一至2000國家重力基本網，重力點坐標全部統一至CGCS2000坐標系，重力點高程全部統一至正常高系統(1985國家高程基準)。

數字高程模型DEM數據生成3″×3″、30″×30″、2.5′×2.5′格網平均高數據。利用3″×3″、30″×30″、2.5′× 2.5′格網平均高數據，完成3″×3″、30″×30″、2.5′×2.5′格網地形改正和均衡改正成果。采用插值計算方法完成離散點重力觀測值的布格、均衡異常歸算；利用離散點的均衡重力異常值，采用線性移動擬合法計算30″× 30″格網平均均衡異常；由30″×30″格網平均均衡異常，采用平均法計算2.5′×2.5′格網平均均衡異常計算；由2.5′×2.5′格網平均均衡異常，計算2.5′×2.5′格網平均空間異常。

(5)GPS數據處理

充分收集蕪湖市周邊的國家GPS連續運行參考站、蕪湖市CORS站的多天觀測數據，采用國際上高精度的數據處理軟件GAMIT/GLOBK最新版本、嚴密的數據處理技術方法，逐級控制，獲得了蕪湖市高精度的2000國家大地坐標系下GPS網數據處理成果。首先，在2000國家大地坐標系下，約束泰安(TAIN)、上海(SHAO)、武漢(WUHN)、廈門(XIAM)4個GPS連續運行站，做三維約束平差，求出蕪湖市CORS站點坐標。其次，在2000國家大地坐標系下，約束繁昌縣站(FCXZ)、南陵縣站(NLXZ)、蕪湖市站(WHSZ)、蕪湖縣站(WHXZ)4個蕪湖CORS站及經分析較可靠的重合安徽省GPSC級點R532、R562、XQDZ做三維約束平差，求出蕪湖市其他GPS網點坐標。

4.2 水準網數據處理

水準網平差以觀測高差為元素，按測站數，采用水準網間接平差方法完成平差計算。高程起算數據采用1999年國家第二期一等水準網復測平差后高程，重力基準為2000國家重力基準；正常重力采用CGCS2000橢球相應公式，二等水準觀測高差中加入標尺尺長改正、正常水準面不平行改正、重力異常改正、固體潮改正四項改正。將成果與前期水準成果進行了分析比較，獲得了現勢性強、精度高、統一的蕪湖市現代高程基準數據。

4.3 大地水準面精化

在完成區域重力似大地水準面計算后，需要利用測區的GPS水準點成果將區域重力似大地水準面擬合適配于該區域的實測似大地水準面。擬合糾正時，利用GPS水準點上的實測似大地水準面ζGPS與由規則格網內插的重力似大地水準面ζgra的差值△ζ，采用多種擬合方法對區域重力似大地水準進行擬合糾正。本項目分別采用最小二乘配置、BP神經網絡、薄板樣條等方法進行擬合糾正，選擇最佳方案確定蕪湖市似大地水準面模型。

4.4 坐標轉換

根據蕪湖市似大地水準面精化GPS點位的分布情況及蕪湖市勘察測繪設計研究院提供的資料，經過綜合分析，反復試算，最終采用重合點11點計算坐標轉換參數。采用四參數坐標轉換模型完成蕪湖市2000國家大地坐標系與1980西安坐標系、1954年北京坐標系、蕪湖獨立坐標系之間的坐標轉換關系，完成各GPS點1980西安坐標系、1954年北京坐標系、蕪湖獨立坐標系成果計算。

5 項目成果精度

(1)項目通過高精度的GPS觀測及數據處理，獲得了蕪湖市4個CORS站點、67個GPS點的2000國家大地坐標系成果。

CORS站點南北方向的精度平均值為±0.4 mm，東西方向的精度平均值為±0.5 mm，高程方向的精度平均值為±1.4 mm。CORS站相鄰點基線南北方向分量測量的精度平均值為±0.4 mm，東西方向分量測量的精度平均值為±0.5 mm；垂直分量測量的精度平均值為±1.8 mm。

GPS點南北方向的精度平均值為±1.2 mm，東西方向的精度平均值為±1.3 mm，高程方向的精度平均值為±5.4 mm。GPS網相鄰點基線南北方向分量測量的精度平均值為±1.7 mm，東西方向分量測量的精度平均值為±1.8 mm；垂直分量測量的精度平均值為±7.7 mm。

(2)分別采用蕪湖市坐標轉換重合點，利用二維高斯平面坐標轉換模型，建立了蕪湖市2000國家大地坐標系與1980西安坐標系、1954年北京坐標系及蕪湖獨立坐標系的坐標轉換模型，并且新建了蕪湖2000獨立坐標系。其中2000國家大地坐標系向蕪湖市獨立坐標系轉換的精度為:X方向精度為±0.009 m，Y方向精度為±0.016 m，平面精度為±0.018 m。

(3)通過高等級的二等水準連續，建立了的蕪湖市二等水準網，獲得了現勢性強、統一的蕪湖市高程基準，其平差后權(每千米)中誤差為±0.34 mm。

(4)充分利用蕪湖市及周邊地區的重力點成果、數字高程模型、采用EGM2008全球重力場模型、及分布較均勻的GPS水準成果，采用重力法(Molodensky原理)及移去～恢復(remove～restore)技術，綜合利用多種現代擬合方法完成了分辨率為2.5′×2.5′，高精度、高分辨率的蕪湖市似大地水準面模型的計算，模型擬合精度達到±1.1 cm，外部檢核精度達到±1.3 cm。

6 小 結

蕪湖市似大地水準面及其相關成果的確定，是常規大地測量、衛星大地測量、現代地球重力場確定理論與方法的綜合運用，不僅建立了與2000國家大地基準一致的精確的區域大地測量框架，同時也標志著集平面、高程及似大地水準面模型為一體的蕪湖市現代測繪基準體系的建立。

似大地水準面模型成果通過與蕪湖市衛星定位綜合服務系統有機結合使用，可以快速地獲取地面點的高程數據，取代低等級水準測量，極大地改善傳統高程測量作業模式，從而使費用高、難度大、周期長的傳統水準測量工作量減少到最低限度，對滿足目前測繪工程的迫切需要，滿足經濟建設日益發展的需要和加快數字蕪湖建設發展需要，都具有特別重要的現實意義和較高的實用價值。

[1] 黃勁松，魏二虎.GPS測量操作與數據處理[M].武漢:武漢大學出版社，2004.

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

[3] 肖建華，王厚之，李江衛等.武漢市現代測繪基準體系建設[J].城市勘測，2007，6，11～13.

[4] 劉道明.合肥市厘米級大地水準面的確定[J].城市勘測，2011(3).

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

Quasi-Geoid Refinement Based on the EMG2008 M odel of W uhu City

Li Huaping1，Ma Xinying2，Zhang Weizhi1，Ma Xiang1，Xie Fang2
(1.Wuhu Geotechnical and Survey Design Institution Co.，Ltd.Wuhu 41000，China；2.Geodetic Data Processing Certer of Nationl Surveying Bureau，Xi’an 710054，China)

This paper discuss aboutusing high-precision GPS leveling data，gravity data，digital terrainmodels，EGM 2008 Earth gravity field model and removed-recovery techniques in 1400 square kilometers quasi-geoid refinement of Wuhu City.We use coordinates coincident points ofWuhu city and two-dimensional Gaussian plane coordinate transformation model to establish coordinate transformation model of geodetic coordinate system ofWuhu City 2000，Xi'an Coordinate System 1980，Beijing coordinate system 1954 and Wuhu independent coordinate system respectively.We comprehensive utilize a variety ofmodern fittingmethod to complete the high-precision，high-resolution quasi-geoid refinement ofWuhu city calculation with the complete resolution of 2.5′×2.5′，model fitting accuracy of±1.1cm external checking accuracy±1.3cm.

similar geoid refinement；GPS；leveling；gravity

1672-8262(2013)05-101-04

P22

A

(1.蕪湖市勘察測繪設計研究院有限責任公司，安徽蕪湖 241000； 2.國家測繪局大地測量數據處理中心，陜西西安 710054)

2012—11—16

李華平(1970—)，男，高級工程師、國家注冊測繪師、國家注冊監理師，現從事工程測繪及數據處理工作。