吳寒,吳燕蘋,鄧明鏡

(1. 重慶交通大學(xué) 土木工程學(xué)院,重慶 400074;2. 重慶工程職業(yè)技術(shù)學(xué)院,重慶 402260)

0 引 言

隨著全球衛(wèi)星導(dǎo)航系統(tǒng)(GNSS)技術(shù)的快速發(fā)展,GNSS高精度定位已廣泛應(yīng)用于工程建設(shè),其中GNSS高程測量由于其不同于傳統(tǒng)水準(zhǔn)作業(yè)方式,極大地提高了高程測量效率．但是GNSS測量的高程為大地高,并不能直接應(yīng)用于工程實(shí)際需要,如何將大地高轉(zhuǎn)換為工程建設(shè)所需的正常高是一個重要研究內(nèi)容,其關(guān)鍵在于高程異常的解算[1]．

目前,高程異常解算的方法主要有物理重力法和幾何擬合法[2-3],幾何擬合法不深入考慮異常值造成的各種因素,統(tǒng)一采用數(shù)學(xué)逼近方法進(jìn)行解決,比較適用于小區(qū)域或地形平坦區(qū)域的工程項(xiàng)目,方法簡單,應(yīng)用廣泛,但在地形起伏較大或大范圍區(qū)域,僅使用幾何擬合法不能滿足實(shí)際需要．隨著衛(wèi)星重力測量的發(fā)展,利用衛(wèi)星重力資料并結(jié)合地面重力數(shù)據(jù)獲得的全球重力場模型的精度和分辨率在逐步提高,特別是EGM2008模型的公布,使得利用重力場模型進(jìn)行GNSS高程轉(zhuǎn)換的精度和可靠性獲得了大幅度提升．但是單一的擬合方法僅在一定的區(qū)域范圍和地形條件下適用,為提高高程擬合的適用性、可靠性和精度,本文綜合利用EGM2008模型和移動曲面模型對GNSS高程進(jìn)行擬合．

1 原理與方法

1.1 整體流程

首先將高程異常分為兩個部分,第一部分是由重力場模型求得的高程異常,第二部分是剩余高程異常部分,具體為

ζ=ζEGM+ζRES，

(1)

式中:ζ為總的高程異常;ζEGM為利用EGM2008模型求解所得的高程異常值;ζRES為剩余高程異常值．

選擇“移去-擬合-恢復(fù)”法求解高程異常,即根據(jù)少量已知水準(zhǔn)點(diǎn),利用EGM2008模型求得對應(yīng)的高程異常ζEGM,然后結(jié)合已知的高程異常ζ求得剩余高程異常值ζRES．考慮到待定點(diǎn)高程異常受到周邊區(qū)域影響,采用以待定點(diǎn)為中心建立剩余高程異常值的移動曲面模型,求出待定點(diǎn)的剩余高程異常值ζRES,再加上對應(yīng)的ζEGM即可求出待定點(diǎn)的高程異常[4-5]．

1.2 EGM2008重力場模型求解高程異常

EGM2008重力場模型是美國國家地理空間情報(bào)局在充分利用地面重力、衛(wèi)星測高、衛(wèi)星重力等多種數(shù)據(jù)的基礎(chǔ)上研制發(fā)布的2 190次球諧系數(shù)的全球超高階地球重力場模型,該模型的基本格網(wǎng)分辨率為5′,地面數(shù)據(jù)覆蓋率達(dá)83.8%,該模型能夠精確表達(dá)高程異常的長波項(xiàng)．利用EGM2008重力場模型計(jì)算高程異常的計(jì)算公式

(2)

章傳銀等[6]對EGM2008地球重力場模型在中國大陸適用性進(jìn)行分析,結(jié)果表明EGM2008模型具有很高的精度,在我國大陸的精度與全球范圍內(nèi)的精度相當(dāng),總體精度為20 cm,華東華中地區(qū)12 cm,華北地區(qū)達(dá)到9 cm,西部地區(qū)為24 cm．束蟬方等[7]進(jìn)一步對EGM2008模型適用性分析認(rèn)為將EGM2008模型和GPS結(jié)合代替低等級的水準(zhǔn)測量是可行的,特別是對于長距離的跨障礙物的高程傳遞更為經(jīng)濟(jì)有效．榮敏等[8]利用中國854個GPS/水準(zhǔn)點(diǎn)數(shù)據(jù)對EGM2008重力場模型的精度進(jìn)行分析,由模型計(jì)算得到的大地水準(zhǔn)面高與實(shí)測值之間存在明顯的系統(tǒng)偏差．綜合分析,僅使用EGM2008模型求解高程異常不能滿足實(shí)際工程需要,還需對剩余的高程異常值進(jìn)行處理．

1.3 剩余高程異常值的處理

在實(shí)際的工程應(yīng)用中,高程擬合方法很多,當(dāng)擬合區(qū)域地形平坦,范圍較小時,使用二次曲面函數(shù)擬合方法就能滿足一般工程需求,當(dāng)?shù)匦纹鸱兇蠛蛿M合區(qū)域變大,二次曲面擬合函數(shù)已無法滿足工程需求,而其他常用的模型如多面函數(shù)法、神經(jīng)網(wǎng)絡(luò)法等也存在計(jì)算復(fù)雜或結(jié)果隨機(jī)等問題,綜合考慮精度、可靠性和計(jì)算復(fù)雜性等因素,本文采用距離定權(quán)的移動曲面模型進(jìn)行剩余高程異常值的處理,移動曲面模型既能很好地解決擬合面的整體性問題,抵償部分系統(tǒng)偏差,同時考慮到待定點(diǎn)受到周邊地形環(huán)境的影響情況,利用不同權(quán)重能更好地模擬待定點(diǎn)擬合面,適用面較廣,靈活性較強(qiáng),在GNSS高程轉(zhuǎn)換中有著很好的優(yōu)勢[9-12]．具體處理方法如下．

以待定點(diǎn)為中心,建立剩余高程異常ζRES與對應(yīng)已知點(diǎn)坐標(biāo)(x,y)關(guān)系模型如式(3)．

ζRES=f(x,y)+ε，

(3)

式中:ε為兩者的殘差;f(x,y)為剩余高程異常的趨勢面,常用的數(shù)學(xué)函數(shù)有曲面函數(shù),平面函數(shù),多項(xiàng)式曲線函數(shù),在實(shí)際工程應(yīng)用中,一般選擇二次曲面擬合函數(shù),其函數(shù)表達(dá)式為

f(x,y)=a0+a1x+a2y+a3x2+

a4y2+a5xy.

(4)

式(3)寫成矩陣形式:

ζRES=AX+ε.

(5)

根據(jù)上述模型可知,測區(qū)已知點(diǎn)數(shù)量n≥6,為提高模型精度,參與模型解算的已知點(diǎn)可結(jié)合區(qū)域地形變化、區(qū)域大小、已知點(diǎn)分布和數(shù)據(jù)精度等諸多因素綜合考慮．已知點(diǎn)與待定點(diǎn)的相關(guān)性一般采用距離定權(quán)方法[13-16],本文采用向徑距離定權(quán)方式,即P=1/(r2+Δh2),Δh為已知點(diǎn)與待定點(diǎn)的高程差,距離r可以是選擇兩點(diǎn)之間的向量長度或兩點(diǎn)的平面距離．

按最小二乘原理,即可求出六個參數(shù)a0、a1、a2、a3、a4、a5,由此可計(jì)算出待定點(diǎn)的剩余高程異常值．

1.4 幾點(diǎn)說明

GNSS高程擬合精度除受到轉(zhuǎn)換模型影響外,還受到諸多因素影響,在數(shù)據(jù)處理過程中應(yīng)綜合考慮．

1) 已知點(diǎn)的測量誤差,如GNSS觀測誤差和水準(zhǔn)測量誤差．

2) 已知點(diǎn)的分布,相比曲線和曲面擬合等方法,盡管移動曲面擬合模型受點(diǎn)位分布的影響較小,但在實(shí)際工程中盡量保證足夠的已知點(diǎn)分布于擬合區(qū)域范圍內(nèi)提高成果的可靠性．

3) 轉(zhuǎn)換范圍,從全國大范圍來看,不同區(qū)域水準(zhǔn)網(wǎng)本身也存在系統(tǒng)偏差,且差值各不相同,我國部分學(xué)者曾認(rèn)為國家一等水準(zhǔn)網(wǎng)在西部地區(qū)存在分米量級的系統(tǒng)偏差,有的甚至認(rèn)為系統(tǒng)偏差超過0.5 m．因而數(shù)據(jù)處理前需對所在區(qū)域的水準(zhǔn)網(wǎng)成果數(shù)據(jù)進(jìn)行分析,對于采用不同水準(zhǔn)網(wǎng)平差的區(qū)域或跨越不同省市地區(qū)的較大區(qū)域,可采用分區(qū)方式處理提高成果的精度和可靠性．

2 實(shí)驗(yàn)與分析

為驗(yàn)證本文方法的精度情況和適用范圍,選取4個不同區(qū)域和地形環(huán)境的算例進(jìn)行實(shí)驗(yàn)分析．

算例1:東西長約5 km,南北長約8 km,共有13個聯(lián)測已知點(diǎn)(四等水準(zhǔn)點(diǎn)),選取其中7個點(diǎn)作為參考點(diǎn),6個點(diǎn)作為檢核點(diǎn),最高高程為326 m,最低高程為278 m,平均高程為314 m,地勢平坦．

算例2:東西長約25 km,南北長約35 km,共有27個聯(lián)測已知點(diǎn)(四等水準(zhǔn)點(diǎn)),選取其中7個點(diǎn)作為參考點(diǎn),20個點(diǎn)作為檢核點(diǎn),最高高程為508 m,最低高程為192 m,平均高程為312 m．

算例3:東西南北長約150,共有18個聯(lián)測已知點(diǎn)(四等水準(zhǔn)點(diǎn)),其平均間距為20,選取其中11個點(diǎn)作為參考點(diǎn),7個點(diǎn)作為檢核點(diǎn),最高高程為1 501 m,最低高程為187 m,平均高程為704 m,地勢起伏很大．

算例4:東西南北長約400 km,共有71個聯(lián)測已知點(diǎn)(四等水準(zhǔn)點(diǎn)),其平均間距為35 km,選取其中34個點(diǎn)作為參考點(diǎn),37個點(diǎn)作為檢核點(diǎn),最高高程為1 548 m,最低高程為71 m,平均高程為527 m,地勢起伏很大．

為了比較解算精度,以常用的二次曲面模型作比較,分別采用二次曲面模型(方法1)、移動曲面模型(方法2)、EGM2008模型和二次曲面模型組合(方法3)、EGM2008模型和移動曲面模型組合(方法4)4種方法進(jìn)行計(jì)算,統(tǒng)計(jì)情況如表1所示．

表1中,RMSE為均方根誤差,|ν|max為擬合殘差絕對值最大值,從表1可知,在平坦和較小范圍區(qū)域,采用4種方法的效果相當(dāng),隨著擬合范圍和地形起伏變大(如算例3和算例4),采用方法1和方法2已無法滿足現(xiàn)實(shí)需求,由于EGM2008重力場模型能表達(dá)高程異常長波項(xiàng),即使在大范圍區(qū)域轉(zhuǎn)換時也能有效地處理好高程異常系統(tǒng)偏差．因而,綜合利用EGM2008模型的方法3和方法4表現(xiàn)出明顯的優(yōu)勢．

表1 擬合精度比較mm

方法算例1RMSE|ν|max算例2RMSE|ν|max算例3RMSE|ν|max算例4RMSE|ν|max 方法12.43.314.826.3225.2364.5621.0847.7 方法22.53.415.935.9194.3208.3188.5539.8 方法33.44.66.818.535.156.066.2156.8 方法43.04.27.89.129.247.448.187.4

圖1 檢核點(diǎn)殘差絕對值與四等水準(zhǔn)限差/m

1) 方法4擬合中誤差較方法3明顯降低,精度改善38%．

2) 方法4和方法3在各點(diǎn)的擬合殘差分布不一致,62%的檢核點(diǎn)中方法4擬合殘差小于方法3．

3) 方法4最大擬合殘差為87.4 mm,方法3最大擬合殘差為156.8 mm,按照四等幾何水準(zhǔn)的限差規(guī)范,方法4所有檢核點(diǎn)無殘差超限情況,方法3有4個點(diǎn)超限,方法4可靠性較好．

3 結(jié)束語

本文基于“移去-擬合-恢復(fù)”的思想,提出綜合利用EGM2008模型的移動曲面模型高程擬合方法,該方法對已知點(diǎn)分布要求較低,計(jì)算簡單,適用范圍較廣,可靠性高,該方法在實(shí)際工程項(xiàng)目尤其是大范圍山區(qū)地形GNSS高程擬合中起到一定改進(jìn)作用．