重力數據誤差對大地水準面模型建立的影響

李姍姍　曲政豪

1　信息工程大學地理空間信息學院，鄭州市科學大道62號，450001

?

重力數據誤差對大地水準面模型建立的影響

李姍姍1曲政豪1

1信息工程大學地理空間信息學院，鄭州市科學大道62號，450001

基于Stokes理論建立的大地水準面模型，其精度受重力數據誤差，即重力數據分辨率、精度以及積分范圍的影響。針對這一問題，通過重力場譜特征分析，給出不同地形區域重力數據分辨率以及積分半徑造成的大地水準面高頻截斷誤差的量級大小，計算平均重力異常誤差對大地水準面建模精度的影響。研究成果對不同地形區域cm級大地水準面模型的建立具有理論與指導意義。

截斷誤差；大地水準面；平均重力異常；代表誤差；積分半徑

建立高精度高分辨率的大地水準面模型需要有地面上的重力異常。因為測量總是在地面離散點上進行，不能滿足理論上連續分布的數據要求，因此以重力異常計算大地水準面時，都是使用數值化的方法來實現。最為常見的是將積分的球面劃分成若干個面積，以該面積內或其周圍已知的重力異常推求出該面積的平均重力異常，再用求和的方法完成積分的數值計算。由上述或等價的其他數值方法確定大地水準面將產生來自兩個方面的誤差：一是模型誤差，即解式本身的誤差(如地形起伏影響、球近似的影響等)，可通過提高解式的嚴密性來減弱[1-4]；二是源于重力數據的誤差，體現在平均重力異常的分辨率(即面積大小)、積分半徑的選取及其精度兩個方面。本文擬針對1′×1′ cm級大地水準面模型的建立，研究分析不同地形區域對積分半徑選取以及平均重力異常精度的要求。

1　重力數據的分辨率與精度

利用Stokes-Pizzetti公式計算大地水準面可表示為[5-6 ]：

(1)

(2)

在實際應用中，顯然不可能采用全球積分。對于局域大地水準面模型的建立，一般采用求和的方法完成積分的數值計算，即將(1)式表示成如下的形式：

(3)

從式(3)可以看出，大地水準面模型建立的精度主要受兩部分影響：一方面，平均重力異常的格網面積越小，即數值積分的單元面積劃分越細，積分范圍越大，則數值積分的結果就越精確；另一方面，水準面的精度還受到網格平均重力異常誤差的影響，這一誤差主要因為觀測值稀疏與數據處理不完善造成?？蓪⒋蟮厮疁拭娼Ｕ`差分為由平均重力異常分辨率以及積分半徑引起的截斷誤差部分m1和平均重力異常誤差所引起的誤差部分m2，即

(4)

2　截斷誤差計算

取計算點周圍(球冠半徑為ψ0)1′×1′分辨率的重力異常數據，結合EGM2008模型，采用移去-恢復技術計算大地水準面，位系數誤差的影響可通過具體的模型誤差來分析，這里暫且將它忽略。則截斷誤差主要包含兩部分：

第一部分，是僅采用1′×1′重力數據計算大地水準面所造成的影響，即在整個球面上忽略了頻段n=10 801～的重力數據而導致的高頻截斷誤差：

(5)

其中，誤差估計為：

(6)

第二部分，由于積分半徑(球冠半徑ψ0)以外區域采用的是2 160階的EGM2008模型，相當于5′×5′的重力數據，也就是說，在該積分區域外還存在n=2 161～10 800的截斷誤差：

(7)

其中，誤差估計為：

(8)

式中，Qn為截斷系數，可表示為[7]：

(9)

研究推導Qn的文獻很多。忽略其推導步驟，其計算式如下[7]：

(10)

式中，

Jn(ψ0)=

J0(ψ0)=I1(ψ0)

Kn+1(ψ0)-2Kn(ψ0)+Kn-1(ψ0)=

利用Rapp-Tscherning模型估計得重力異常階方差為[8]：

(11)

式中，A=425.28，S=0.999 623，B=24，C2=7.5 (mGal)2。

根據式(6)和式(11)，計算得到1′×1′重力數據分辨率造成的高頻截斷誤差為0.63 cm。然后取不同的積分半徑，依據式(8)～(11)計算積分半徑外的截斷誤差，計算結果如表1所示。

表1　不同積分半徑造成的截斷誤差

3　不同地形區域重力場的譜特征

當利用式(6)、式(8)來估計分辨率以及遠區高頻截斷部分造成的誤差并以此確定平均重力異常積分半徑時，是根據重力異常階方差Cn估算的。而由式(11)給出的階方差模型(以及其他一些已知模型)是國外學者通過全球重力異常的統計分析得出的經驗公式，因此代表了全球重力異常的一般特征，但各個不同地區重力場的具體情況是有差異的。數據分辨率誤差的本質是平均重力異常與其代表面積內各點重力異常的差異，即代表誤差。對于局部地區，通常采用下面的代表誤差模型[8]：

E2=4C2D

(12)

其中，D是平均重力異常矩形網格的邊長，以km為單位；C稱為代表誤差系數，它與各地區的地形情況(反映了重力場的復雜度)有關。

由于

(13)

式中，角標K代表不同的地形類別。由代表誤差的經驗公式可見，不同地形類別的代表誤差僅相差一個常數因子，故可?。?/p>

(14)

即

(15)

從而，可由C求得相應的常數AK：

(16)

表2　不同地形區域AK統計

因而，考慮到不同的地形類別，式(6)與式(8)的高頻截斷誤差分別變為：

(17)

(18)

根據式(17)計算的不同地形區域重力數據1′×1′分辨率造成的高頻截斷誤差如表3所示。

表3　不同地形區域重力數據1′×1′分辨率

從表3可以看出，1′×1′數據分辨率造成的高頻截斷誤差，即使在特大山區，也僅為1.88 cm，因此1′×1′重力數據分辨率能滿足cm級大地水準面建模的精度要求。根據式(18)計算的不同地形區域不同積分半徑造成的截斷誤差示于表4。

從表4可以看出，當積分半徑取至20′時，在平原、丘陵和小山區均能滿足cm級大地水準面建模的要求；但對于中山區與大山區而言，則積分半徑至少需要取至30′和70′；而對于特大山區，積分半徑即使取至2°也無法滿足cm級建模的精度要求。繼續擴大特大山區積分區域，結果示于表5。當積分半徑取至3°時，在特大山區積分半徑造成的截斷誤差為9.65 cm。

表4　不同積分半徑在不同地形區域造成的截斷誤差

表5　不同積分半徑在特大山區造成的截斷誤差

4　重力數據誤差影響

僅考慮近區平均重力異常誤差對大地水準面的影響，將式(3)表示為：

(19)

式中，S為線性化Stokes積分算子，δΔg為重力異常值向量。設平均重力異常的方差為CδΔg，根據誤差傳播定理，可得大地水準面的方差CδN為：

(20)

因此，平均重力異常中誤差與大地水準面高中誤差的關系為：

(21)

考慮到中央區S(ψ)的計算存在奇異問題，因此將式(20)表示為中央區σ0與扣除中央區σ-σ0的組合：

(22)

式中，MδN(σ-σ0)計算如式(21)。而在中央區，由于[6]

(23)

代入到式(1)，得到平面近似下中央區大地水準面的計算公式[1]：

(24)

其中，

(25)

通過式(25)，依據誤差傳播定理，將Δg(x,y)的誤差傳播給擬合系數aij，再通過式(24)計算MδN(σ0)。

選擇不同積分半徑計算3mGal、5mGal、10mGal以及20mGal平均重力異常誤差對大地準面計算精度的影響，計算點位置(33°，110°)，計算結果如表6所示(分辨率1′×1′，單位cm)。

表6　不同積分半徑平均重力異常數據誤差

從表6可以看出，隨著積分半徑的增大，平均重力異常誤差對大地水準面建模精度的影響逐漸增大。綜合表3、表4與表6，我們可以根據局域大地水準面建模精度的要求，選擇合適的積分半徑，并對1′×1′平均重力異常精度提出要求。例如假定某計算區域為一大山區，要求大地水準面精度為10 cm，我們選定積分半徑為80′，該區域1′×1′平均重力異常精度為10 mGal，則大地水準面建模精度為11.23 cm；如果該區域1′×1′平均重力異常精度為3 mGal，則大地水準面建模精度為10.12 cm。選定積分半徑為120′，平均重力異常精度為10 mGal，則大地水準面建模精度為9.21 cm；如果該區域平均重力異常精度為3 mGal，則大地水準面建模精度為8.04 cm。顯然，積分半徑越大，平均重力異常精度越高，建模精度就越高?？紤]到重力測量的工程量，在建立局域大地水準面模型之前，需要根據建模精度的需求選擇合適的積分半徑與平均重力異常精度，盡可能地滿足建模精度需求，且測量任務又不至于過于繁重。

5　結　語

本文研究分析了1′×1′cm級大地水準面建模精度對不同地形區域積分半徑選取以及平均重力異常精度的要求。實驗表明，當積分半徑取至20′時，在平原、丘陵和小山區均能滿足cm級大地水準面的建模要求；對于中山區與大山區而言，則積分半徑至少需要取至30′和70′；而對于特大山區，當積分半徑取至3°時，其造成的截斷誤差為9.65 cm。同時，給出了3 mGal、5 mGal、10 mGal以及20 mGal平均重力異常誤差對大地準面建模精度影響的量級大小。綜合計算結果，我們可以選擇合適的積分半徑以及平均重力異常精度，使之既滿足建模精度需求，又避免冗余繁重的工程測量任務。本文的研究成果對于計算確定局域cm級大地水準面模型具有參考價值與指導意義。

[1]李姍姍，吳曉平，張傳定，等．顧及地形與完全球面布格異常梯度項改正的區域似大地水準面精化[J]．測繪學報，2012, 41(4)：510-516(Li Shanshan, Wu Xiaoping, Zhang Chuanding,et al. Regional Quasi-Geoid Refining Considering Corrections of Terrain and Complete Spherical Bouguer Anomaly’s Gradient Term[J]．Acta Geodaetica et Cartographia Sinica, 2012, 41(4):510-516)

[2]陳俊勇，李建成，寧津生，等．我國大陸高精度、高分辨率大地水準面的研究和實施[J]．測繪學報，2001，30(2)：95-99(Chen Junyong, Li Jiancheng, Ning Jinsheng, et al. On a High Resolution and High Accuracy Geoid in China Mainland[J]．Acta Geodaetica et Cartographia Sinica, 2001, 30(2): 95-99)

[3]陳俊勇，李建成，寧津生，等．中國新一代高精度、高分辨率大地水準面的研究和實施[J]．武漢大學學報:信息科學版，2001，26(4)：283-289(Chen Junyong, Li Jiancheng, Ning Jinsheng, et al. A New Chinese Geoid with High Resolution and High Accuracy[J]. Geomatics and Information Science of Wuhan University, 2001, 26(4): 283-289)

[4]寧津生，羅志才，楊沾吉，等．深圳市1 km高分辨率cm級高精度大地水準面的確定[J]．測繪學報，2003，32(2)：102-107(Ning Jinsheng, Luo Zhicai, Yang Zhanji, et al. Determination of Shenzhen Geoid with 1 km Resolution and Centimeter Accuracy[J]．Acta Geodaetica et Cartographia Sinica, 2003, 32(2): 102-107)

[5]Moritz H. Physical Geodesy [M]. New York:Springer-Verlag Wien, 2005

[6]郭俊義．物理大地測量學基礎[M]．武漢：武漢測繪科技大學出版社，1994(Guo Junyi. Foundation of Physical Geodesy[M]. Wuhan: Wuhan Techonology University of Surveying and Mapping Press, 1994)[7]Zhang C D , Lu Z L, Wu X P. Truncation Error Formulae for the Disturbing Gravity Vector[J]．Journal of Geodesy, 1998, 72: 119-123

[8]陸仲連．地球重力場的理論與方法[M]．北京：解放軍出版社，1996(Lu Zhonglian. Theory and Method of the Earth’s Gravity Field[M]. Beijing: PLA Publishing House, 1996)

About the first author:LI Shanshan, PhD, professor, PhD supervisor, majors in physical geodesy, E-mail: zzy_lily@sina.com.

Effect of Error of Gravity Data on Geoid Determination

LIShanshan1QUZhenghao1

1Institute of Surveying and Mapping, Information Engineering University, 62 Kexue Road, Zhengzhou 450001, China

Determining accurate regional geoid with Stokes theory can be affected by gravity data errors, such as resolution of gravity, precision of mean gravity and integral radius. In order to solve this problem, the magnitude of truncation errors of geoids in different terrain areas, which are caused by gravity resolution and integral radius, are given based on the analysis of characteristics of the gravity frequency spectrum. Also, the effect of error of mean gravity anomalies on geoid construction is calculated. Therefore, the results put forward in this paper contain great theoretical value for centimeter geiod calculation in different terrain areas.

truncation error; geoid; mean gravity anomaly; representation error; integral radius

National Natural Science Foundation of China, No.41274029; National High Technology Research and Development Program, No.2013AA122502.

2015-10-30

李姍姍，博士，教授，博士生導師，主要從事物理大地測量學研究，E-mail:zzy_lily@sina.com。

10.14075/j.jgg.2016.10.001

1671-5942(2016)010-0847-04

P223

A

項目來源：國家自然科學基金(41274029)；國家863計劃(2013AA122502)。