全球重力場模型EGM2008精度分析與改進研究

■陳秋陽，陳尚云

■1.中國水電顧問集團中南勘測設計研究院有限公司，湖南 長沙 410014;2.四川中水成勘院測繪工程有限責任公司，四川 成都 610072

1 前言

隨著美國全球定位系統(GPS)、俄羅斯全球衛星導航系統、中國北斗衛星導航系統(BDS)、和歐盟衛星導航定位系統(GALILEO)的不斷完善和投入民用，GNSS測量的可靠性和精度逐步提高，其應用范圍也越來越廣。

由于GNSS測量得到的成果是WGS84坐標系的地心坐標(X，Y，Z)或大地坐標(B，L，H)，而用戶需要的成果則是參心坐標(如1954北京坐標或1980西安坐標)和正常高，因此，需要將大地坐標(B、L、H)轉換為參心坐標(x、y)和正常高h。平面坐標轉換相對容易實現，而將大地高H轉換為正常高h則需要確定精確的高程異常ζ(即h=H－ζ)，由于似大地水準面的不規則性，大地高H轉換為正常高h的精度較低。為提高GPS大地高H的利用效率，近來由NGA(US National Geospatial－IntelligenceAgency)推出的全球超高階地球重力場模型EGM2008以確定全球似大地水準面，該模型的階次為2159(另外球諧系數的階擴展至2190次)，相當于模型的空間分辨率約為5＇(約9 km)。該模型采用了GRACE衛星跟蹤數據(ITG－GRACE03S位系數信息以及相應的協方差信息)、衛星測高數據和地面重力數據等，該模型無論在精度還是在分辨率方面均取得了巨大進步，采用該模型以及GPS水準數據有望獲得更高精度的區域似大地水準面。

下面我們以工程實例對全球重力模型EGM2008的精度進行分析，并提出精度改進方法。

2 EGM2008模型精度分析

2.1 工程1分析

工程1施測了16個控制點，控制面積約80km2，控制點分布見圖2－1，各點大地坐標與正常高值見表2－1。

圖2 －1

表2 －1 GPS控制點大地坐標與正常高值

點號 大地經度L(°＇″) 大地緯度B(°＇″) 大地高H(m) 正常高h(m)10 104.271805332 E 31.481736606 N 755.431 794.077 11 104.272082132 E 31.503576736 N 624.410 662.895 12 104.283992279 E 31.504462179 N 691.031 729.648 13 104.300153986 E 31.513604053 N 599.711 638.412 14 104.255647407 E 31.505334368 N 776.048 814.352 15 104.272371667 E 31.485384672 N 1081.706 1120.310 16 104.261130931 E 31.504508640 N 705.801 744.148

利用EGM2008模型計算表2－1中各點高程異常值，然后求得其正常高值，并與實測值進行比較，比較結果見表2－2。

表2 －2 GPS點EGM2008高程異常和計算高程成果表

表2－2中，第5列為EGM2008模型計算的正常高HEGM2008=大地高h－高程異常ζEGM2008;第6列為HEGM2008與實測正常高的較差△H=HEGM2008－h，最大值為0.437m，最小值為0.374m，中誤差為±0.411m;然后采用常數擬合的方式，取較差均值對HEGM2008進行改正，得到改正后的模型正常高(見第7列);第8列為改正后的正常高與實測正常高的差值，最大值為0.032m，最小值為－0.031m，中誤差為±0.018m。

2.2 工程2分析

工程2施測了21個控制點，控制面積約247290km2，控制點分布見圖2－2，各點成果見表2－3。

表2 －3 GPS控制點大地坐標與正常高值

圖2 －2

利用EGM2008模型計算表2－3中各點高程異常值，然后求得其正常高值，并與實測值進行比較，比較結果見表2－4。

表2 －4 GPS點EGM2008高程異常和計算高程成果表

根據工程1的計算方法，直接利用EGM2008模型計算的正常高與實測正常高間的較差最大值為2.189嗎，最小值為0.041m，中誤差為±1.067m;若利用較差均值進行修正，則較差最大值為1.074m，最小值為－1.074m，中誤差為 ±0.671m，與工程 1相比，該工程區利用EGM2008獲得正常高精度較低。

3 EGM2008模型正常高精度改進

從上述精度分析中可以看出，若直接利用EGM2008計算高程異常獲取正常高值，其精度較低，如工程1為±0.411m，工程2為±1.067m;若利用較差均值對高程異常進行修正，則得到的正常高精度將得到很大的改善，經改正后工程1的精度為±0.018m，工程2為±0.671m;但工程2經均值修正后精度依然偏低。

為獲取更優的正常高值，針對工程2，利用EGM2008模型正常高和實測正常高，對其差值按10參數多項式模型進行差值擬合，然后對EGM2008模型正常高值進行修正，表3－1是針對工程2進行的擬合計算。

表3 －1 EGM2008正常高優化計算

表3－1是利用21個控制點按10參數多項式擬合計算得到的模型和精度統計，該模型的內符合精度為±0.125m，最弱點正常高較差為0.318m;外符合精度為±0.168m。

若直接利用實測的大地高和正常高按10參數多項式模型進行高程異常擬合，模型成果和精度統計見表3－2。

點名 X(m) Y(m) h(m) H(m) 差值(m)1 3562523.880 208689.936 3742.081 3774.449 －32.368至21 3435605.135 496955.486 865.224 900.822 －35.598

利用該法，由21個點按10參數多項式模型擬合獲取的正常高內符合精度為±1.212m，最弱點正常高較差為2.229m;外符合精度為±1.352，最弱點正常高較差為2.410m。

從上分析可知，利用EGM2008模型獲取的正常高精度比直接利用GPS水準點進行模型擬合獲取的正常高精度高;若利用殘差對EGM2008模型正常高進一步改正，則獲取的正常高精度更高。

4 結束語

本文利用工程實例分析了EGMM2008全球重力場模型的精度，并提出了改進EGM2008模型正常高精度的方法，所形成的工作經驗表現在:

(1)當測區范圍不大于20km×20km，且有一個已知高程點時，則根據已知點的EGM2008正常高與實測正常高的差值，對EGM模型正常高進行修正，可大幅提高EGM2008模型正常高精度，如工程1從±0.411m提高至±0.018m;

(2)當范圍大到500km×500km時或EGM2008模型正常高精度較低時，可根據均勻分布于測區的控制點的EGM2008正常高和實測正常高的差值，采用多項式模型進行擬合，然后對EGM2008正常高進行修正，可大幅提高EGM2008模型正常高的精度，如工程2從±0.671m提高至±0.168m;

(3)根據工程實踐經驗，當測區范圍不大于500km×500km、且有均勻分布的控制點時，可不再作圖根水準測量或圖根三角高程測量，可大大減少外業工作，從而提高經濟效益。

［1］孔祥元，郭際明.劉宗全編著《大地測量學基礎》武漢大學出版社，2005(12).

［2］於宗儔，魯林成.《測量平差基礎》(增訂本)，測繪出版社，1983(6).