利用衛(wèi)星測高技術(shù)確定浙江海域大地水準面

梁子亮,岳建平,呂志才,吉淵明

(1 河海大學(xué) 地球科學(xué)與工程學(xué)院，江蘇 南京 210098；2.浙江省測繪科學(xué)技術(shù)研究院，浙江 杭州 310012)

利用衛(wèi)星測高技術(shù)確定浙江海域大地水準面

梁子亮1,岳建平1,呂志才1,吉淵明2

(1 河海大學(xué) 地球科學(xué)與工程學(xué)院，江蘇 南京 210098；2.浙江省測繪科學(xué)技術(shù)研究院，浙江 杭州 310012)

利用多顆衛(wèi)星的測高數(shù)據(jù)，經(jīng)共線平均及交叉點平差，建立浙江深海海域2.5′×2.5′格網(wǎng)分辨率的平均海面高模型，在扣除海面地形影響后得到海域的大地水準面起伏，并與EGM2008所得計算結(jié)果進行對比；利用移去-恢復(fù)技術(shù)及SVR方法，聯(lián)合驗潮站GPS/水準數(shù)據(jù)與EGM2008大地水準面模型，計算浙江近岸海域大地水準面起伏；最終建立浙江海域2.5′×2.5′格網(wǎng)分辨率的大地水準面模型。

衛(wèi)星測高；大地水準面；平均海面；EGM2008；移去恢復(fù)

大地水準面的概念是隨著測量精度的提高而改變的。在相當長的一段時間內(nèi)，平均海面被視為重力等位面，并且與大地水準面重合。近20年來的研究表明，平均海面并不具有這種性質(zhì)，其相對于大地水準面的起伏被稱為海面地形，全球變化幅度為±(1～2) m[1],在平均海面中扣除海面地形的影響即可得到海域大地水準面。

衛(wèi)星測高以其能提供豐富海面高信息的特點，在平均海面的研究中被廣大學(xué)者廣泛采用，如何聯(lián)合處理多代測高衛(wèi)星數(shù)據(jù)，成為計算高精度、高分辨率平均海面高模型的關(guān)鍵技術(shù)[2]。本文聯(lián)合Envisat數(shù)據(jù)、Envisat新軌道數(shù)據(jù)、ERS-1數(shù)據(jù)、ERS-2數(shù)據(jù)、T/P數(shù)據(jù)、T/P新軌道數(shù)據(jù)、Jason-1數(shù)據(jù)、Jason-1新軌道數(shù)據(jù)、Jason-2數(shù)據(jù)和GFO數(shù)據(jù)，建立浙江海域深海2.5′×2.5′格網(wǎng)分辨率的平均海面高模型，在減去MDT_CNES-CLS09 (http://www.aviso.oceanobs.com/)全球海面地形模型值后，得到深海海域的大地水準面起伏。

在近岸區(qū)域，陸地對測高衛(wèi)星的影響較大，其測高數(shù)據(jù)含有較大誤差，由衛(wèi)星測高數(shù)據(jù)確定的平均海面的可靠性不佳，且目前已有的全球海面地形模型在近岸海域存在系統(tǒng)差的幾率較大[1]，因此近岸海域更難從平均海面中分離出海面地形。本文聯(lián)合近海驗潮站的GPS/水準數(shù)據(jù)，利用移去-恢復(fù)及SVR方法計算近岸海域的大地水準面起伏。

1 平均海面的計算

1.1 測高數(shù)據(jù)

衛(wèi)星測高數(shù)據(jù)采用了Aviso發(fā)布的DT CorSSH產(chǎn)品，該數(shù)據(jù)由 Ssalto/Duacs開發(fā)，并由AVISO發(fā)布，產(chǎn)品提供了包括波高和后向散射系數(shù)等多項參數(shù)，且所有衛(wèi)星測高數(shù)據(jù)的基準均已統(tǒng)一到與T/P、Jason-1、Jason-2測高衛(wèi)星相同的坐標框架上。

由圖1可見，所采用的測高數(shù)據(jù)密度較高，且分布也較均勻，可以滿足模型2.5′×2.5′分辨率的要求。

圖1 測高數(shù)據(jù)分布圖

1.2 共線平均

共線平均是一種減小衛(wèi)星軌道誤差并確定平均海面的方法，其基本原理為通過固定的參考軌跡來確定其它周期相對應(yīng)弧段上同緯度點的經(jīng)度及其海面高。經(jīng)過推導(dǎo)可以得到一組上升及下降弧段統(tǒng)一的計算公式[3]為

(1)

(2)

式中:λP為P點經(jīng)度，φP為P點緯度，λQ為Q點經(jīng)度，φQ為Q點緯度，φO為O點緯度，HP為P點海面高，HQ為Q點海面高，D為共線弧的斜率，λ為O′點經(jīng)度，H為O′點海面高。

參考式(1)、式(2)，共線平均的具體步驟如下[4]：①分別對各衛(wèi)星的測高數(shù)據(jù)進行共線平均；②剔除與平均值之差的絕對值大于1 m的海平面觀測值；③重新計算新的平均海面高，并形成各測高衛(wèi)星的平均軌道。對測高數(shù)據(jù)進行共線平均不僅可以減弱海面高的時變影響，還可以減弱在某一特定時期發(fā)生的大范圍海洋學(xué)異?，F(xiàn)象所引起的海平面異常變化。

1.3 交叉點平差

本文利用驗后平差方法進行多種測高數(shù)據(jù)聯(lián)合交叉點平差。平差中，除徑向軌道誤差以外，將其他誤差源也包含在假定的誤差模型中。由于本文研究區(qū)域較小，各個弧段上交叉點較少，因此采用線性函數(shù)和傅里葉函數(shù)之和作為回歸數(shù)學(xué)模型?？紤]到十類測高數(shù)據(jù)基準均統(tǒng)一在Jason-1、Jason-1新軌道以及Jason-2所在的坐標框架，且該框架RMS較小，本文認為：Jason-1、Jason-1新軌道以及Jason-2所在平均框架的測高值含有十類數(shù)據(jù)中最少的中長波誤差，因此，采用Jason-1、Jason-1新軌道以及Jason-2的框架作為平均框架來改正其他測高數(shù)據(jù)平均框架，從空域上削弱中長波誤差。平差前后交叉點不符值精度(RMS)如表1所示(“-”代表該衛(wèi)星未形成交叉點)。

由表1可知，交叉點平差后，各衛(wèi)星交叉點的精度(RMS)均有所提高。由此可知，通過交叉點平差，不僅提高了各衛(wèi)星數(shù)據(jù)的徑向軌道精度，而且也較系統(tǒng)地統(tǒng)一了各種數(shù)據(jù)的參考基準。

表1 平差前后交叉點不符值RMS m

1.4 格網(wǎng)化方法比較

為了獲得較高的格網(wǎng)化精度，利用Surfer 11軟件選取多種方法進行了測高數(shù)據(jù)的格網(wǎng)化，并進行交叉驗證，交叉驗證結(jié)果見表2。

表2 交叉驗證統(tǒng)計

由表2可知，最小曲率法格網(wǎng)化精度最好，其絕對值偏差平均值、RMS、STD均為最小。故最終選取最小曲率法對測高數(shù)據(jù)進行格網(wǎng)化。為了避免近岸與深海大地水準面模型拼接時再次進行格網(wǎng)化帶來的誤差，在對深海大地水準面進行插值的同時，計算了近岸的格網(wǎng)坐標。

2 大地水準面的計算

2.1 深海大地水準面的計算

EGM2008是美國國家地理空間情報局(NGA)發(fā)布的超高階全球重力場模型，該模型在計算時綜合考慮了最新的衛(wèi)星測高、衛(wèi)星重力和地面重力數(shù)據(jù)，模型的階次完全至2159階[5]，近年來的研究表明，該模型在我國深海海域具有較高的精度[6]。

深海大地水準面由格網(wǎng)點平均海面高減去格網(wǎng)點海面地形模型值計算，其中海面地形模型采用CNES-CLS09模型。為了檢核計算結(jié)果，將所得大地水準面起伏與EGM2008計算所得大地水準面起伏進行比較，結(jié)果見表3。

表3 數(shù)據(jù)比較結(jié)果

從表3可以看出，計算所得的浙江深海海域大地水準面起伏與EGM2008所得數(shù)值差值的標準差是±0.085 39 m，RMS為±0.085 52 m。由此可見，本文所建立的浙江深海大地水準面起伏模型是正確的，精度可靠。

2.2 近岸海域大地水準面的計算

本文采用移去-恢復(fù)方法計算近岸大地水準面起伏，參考模型為EGM2008重力場模型，具體做法為[9]：①以EGM2008重力場模型作為參考場，采用移去-恢復(fù)原理[10-11]從驗潮站數(shù)據(jù)中扣除EGM2008模型大地水準面起伏，得到剩余大地水準面起伏Nres；②采用SVR插值方法，訓(xùn)練剩余大地水準面起伏Nres，并以局部擬合半徑為2倍格網(wǎng)間距在格網(wǎng)點上進行插值，若在擬合的范圍內(nèi)少于1個觀測值，則舍去該點；③將網(wǎng)格剩余大地水準面起伏和EGM2008重力場模型的格網(wǎng)值Nmod相加，恢復(fù)得到浙江近岸大地水準面模型。

利用Surfer11軟件將沿海與深海大地水準面模型進行拼接，并繪制成等高距為1 m的等值線柵格圖(L，B，N)，如圖2所示。

圖2 平均海面等值線圖

3 結(jié) 論

本文研究了聯(lián)合處理多代測高衛(wèi)星數(shù)據(jù)的方法。通過共線平均，在時域上削弱了測高數(shù)據(jù)短波誤差的影響，通過以Jason-1、Jason-2平均框架作為參考框架來強制改正其余衛(wèi)星數(shù)據(jù)，從空域上削弱了中長波誤差的影響。研究了數(shù)據(jù)格網(wǎng)化的方法，建立了浙江深海海域2.5′×2.5′格網(wǎng)分辨率的平均海面高模型，并計算了近岸的格網(wǎng)坐標。研究了利用海面高模型與海面地形模型計算深海大地水準面的方法，并將所得大地水準面模型與EGM2008計算所得大地水準面模型進行了比較，所得差值的均方根標準差(STD)為±0.085 39 m，均方根(RMS)為±0.085 52 m。利用移去-恢復(fù)技術(shù)及SVR方法，聯(lián)合驗潮站GPS/水準數(shù)據(jù)與EGM2008大地水準面模型，計算了浙江近岸海域大地水準面，最終聯(lián)合深海數(shù)據(jù)建立了浙江海域2.5′×2.5′格網(wǎng)分辨率的大地水準面模型。

[1]陳俊勇，李建成，晁定波，等. 我國海域大地水準面的計算及其與大陸大地水準面拼接的研究和實施[J]. 地球物理學(xué)報， 2003，46(1): 31-35.

[2]鄧凱亮，暴景陽，許軍，等. 用強制改正法建立中國近海平均海平面高模型[J]. 武漢大學(xué)學(xué)報：信息科學(xué)版，2008，33(12): 1283-1287.

[3]姜衛(wèi)平. 衛(wèi)星測高技術(shù)在大地測量學(xué)中的應(yīng)用[D]. 武漢：武漢大學(xué), 2001.

[4]姜衛(wèi)平，李建成，王正濤. 聯(lián)合多種測高數(shù)據(jù)確定全球平均海面WHU2000[J]. 科學(xué)通報， 2002，47(15):1187-1191.

[5]束蟬方，李斐，郝衛(wèi)峰. EGM2008模型在中國某地區(qū)的檢核及適用性分析[J]. 武漢大學(xué)學(xué)報：信息科學(xué)版，2011，36(8): 919-922.

[6]胡淑梅. 衛(wèi)星測高反演海域重力異常和大地水準面的研究[D]. 阜新：遼寧工程技術(shù)大學(xué), 2012.

[7]陳德忠，趙亮，何書靜，等. BP神經(jīng)網(wǎng)絡(luò)和SVR用于GPS高程擬合研究分析[J]. 測繪地理信息，2012，37(5): 47-49.

[8]翟振和，魏子卿，吳富梅，等. 利用EGM2008位模型計算中國高程基準與大地水準面間的垂直偏差[J]. 大地測量與地球動力學(xué)，2011，31(4): 116-118.

[9]鄧凱亮，暴景陽，章傳銀，等. 聯(lián)合多代衛(wèi)星測高數(shù)據(jù)建立中國近海垂線偏差模型[J]. 海洋測繪，2008，28(3): 15-17.

[10]劉曉剛, 劉雁雨, 曹紀東, 等. GPS水準采用移去恢復(fù)技術(shù)擬合大地水準面方法的研究[J]. 測繪工程, 2008, 17(3): 70-73.

[11]谷延超, 范東明. 顧及EGM2008和殘差地形模型的GPS高程轉(zhuǎn)換方法研究[J].測繪工程, 2013, 22(2): 26-29.

[責任編輯：劉文霞]

Determination of geoid model in zhejiang costal waters by using satellite altimetry

LIANG Zi-liang1,YUE Jian-ping1,Lü Zhi-cai1,JI Yuan-ming2

(1.School of Earth Sciences and Engineering, Hohai University, Nanjing 210098, China；2.Zhejiang Academy of Surveying and Mapping, Hangzhou 310012, China)

Multi-satellite altimeter data is used to establish the Mean Sea Surface (MSS) model over Zhejiang deep waters with a spatial resolution of 2.5′, based on the collinear average method and crossover adjustment method. The geoid undulation is calculated after deducting the impact of sea surface topography compared with the results of EGM2008. The geoid undulation of Zhejiang coastal waters is calculated by using the remove-restore technology and SVR methods with GPS/Leveling data of tide gauges and EGM2008 geoid model. A geoid model over Zhejiang offshore is ultimately established with a spatial resolution of 2.5′.

altimetry; geoid; mean sea surface; EGM2008; remove-restore

2013-08-15

梁子亮(1990-)，男，碩士研究生.

P228

：A

：1006-7949(2014)07-0022-04