沈清華,趙薛強(qiáng)
(中水珠江規(guī)劃勘測(cè)設(shè)計(jì)有限公司,廣東 廣州 510610)
跨海高程傳遞是海島礁測(cè)繪首先需要解決的問(wèn)題,一直以來(lái)深受海洋測(cè)繪學(xué)者的關(guān)注。隨著GNSS技術(shù)、重力測(cè)量技術(shù)和衛(wèi)星重力探測(cè)技術(shù)的迅速發(fā)展,重力場(chǎng)模型以及區(qū)域大地水準(zhǔn)面的精度不斷提升,GNSS跨海高程傳遞技術(shù)成為最為活躍的研究領(lǐng)域之一。2001年,李建成、姜衛(wèi)平等在比較分析靜力水準(zhǔn)法、動(dòng)力水準(zhǔn)法、GNSS水準(zhǔn)法及常規(guī)大地測(cè)量法等目前常用的跨海高程基準(zhǔn)傳遞方法的優(yōu)缺點(diǎn)的基礎(chǔ)上,利用GNSS定位技術(shù)將黃海高程由陸地傳遞到距離上海蘆潮港30km的洋山島上[1]。2013年,王建成、沈清華等,利用GNSS水準(zhǔn)法采用基于EGM2008的“移去-擬合-恢復(fù)”理論與技術(shù)成功將1985國(guó)家高程基準(zhǔn)傳遞到珠江河口以外30km的萬(wàn)山島和50km的擔(dān)桿島上[2]。
由于上述GNSS水準(zhǔn)法是基于大地高差與大地水準(zhǔn)面差距之差進(jìn)行高程傳遞,需要聯(lián)測(cè)足夠數(shù)量的高等級(jí)水準(zhǔn)點(diǎn),同時(shí)對(duì)陸地、海洋重力資料及其他重力場(chǎng)信息和數(shù)據(jù)成果要求較高,在一定程度上制約了此種方法在實(shí)際生產(chǎn)中的應(yīng)用。而基于GNSS連續(xù)運(yùn)行基準(zhǔn)站的同步觀(guān)測(cè)以及區(qū)域似大地水準(zhǔn)面精化模型的擬合推算技術(shù)(以下簡(jiǎn)稱(chēng)CORS同步技術(shù)),由于無(wú)需聯(lián)測(cè)任何已知點(diǎn),將逐漸被廣泛采用,該方法在陸域部分的擬合高程可達(dá)到四等以上水準(zhǔn)精度,但在遠(yuǎn)距離海島區(qū)域,由于基準(zhǔn)站的覆蓋不足,高程擬合精度相對(duì)較低,考慮到伶仃洋海域?qū)俚透叱坍惓^(qū)域,且高程異常變化平緩,而且桂山島架設(shè)了GNSS連續(xù)運(yùn)行基準(zhǔn)站,適合于采用CORS同步技術(shù)進(jìn)行高程擬合從而實(shí)現(xiàn)遠(yuǎn)距離跨海高程傳遞。
我國(guó)高程系統(tǒng)采用的是基于似大地水準(zhǔn)面的正常高系統(tǒng),而GNSS測(cè)量得到的是相對(duì)于WGS84橢球面的大地高,大地高只具有幾何意義,沒(méi)有水準(zhǔn)意義,正常高與大地高H之間的差值稱(chēng)為高程異常,如圖1所示。
圖1 高程異常示意圖
由于大地水準(zhǔn)面是個(gè)不規(guī)則面,因此各點(diǎn)的高程異常是一個(gè)變量,大地高H、高程異常ζ、正常高Hγ三者的關(guān)系如下式。
Hγ=H-ζ
(1)
通過(guò)基于GNSS連續(xù)運(yùn)行基準(zhǔn)站的高等級(jí)GNSS控制網(wǎng)的同步觀(guān)測(cè)和嚴(yán)密解算,經(jīng)過(guò)CGCS2000框架下的三維約束平差,可獲取GNSS控制點(diǎn)的大地坐標(biāo)以及大地高,再通過(guò)區(qū)域似大地水準(zhǔn)面精化模型進(jìn)行擬合與推算,求取對(duì)應(yīng)GNSS點(diǎn)的高程異常從而獲得各點(diǎn)的正常高。
由于伶仃洋海域內(nèi)的桂山島、東澳島、大萬(wàn)山島、外伶仃島四個(gè)海島,距離大陸海岸的距離在20~45km之間,考慮到高程傳遞及成果檢核需要,陸域部分往海岸以?xún)?nèi)等距離范圍進(jìn)行GNSS點(diǎn)布設(shè),按C級(jí)GNSS網(wǎng)要求,共布設(shè)15個(gè)GNSS-C級(jí)點(diǎn),其中陸域10個(gè)點(diǎn),海島5個(gè)點(diǎn),如圖2所示。
圖2 GNSS控制網(wǎng)示意圖
控制網(wǎng)觀(guān)測(cè)采用3臺(tái)Trimble R8和4臺(tái)Leica GS15 GNSS接收機(jī)按國(guó)家C級(jí)網(wǎng)GNSS要求連續(xù)同步觀(guān)測(cè),共5個(gè)時(shí)段,每個(gè)時(shí)段觀(guān)測(cè)8小時(shí)以上。
基線(xiàn)處理時(shí),將沿海及海島區(qū)域的QIAO、ZHGT、GUIS、JWGT、NLGT等5個(gè)GNSS連續(xù)運(yùn)行基準(zhǔn)站納入共同處理,如圖3所示,共組成110個(gè)閉合環(huán),通過(guò)環(huán)檢驗(yàn)的基線(xiàn)128條,平均邊長(zhǎng)30km。
圖3 納入連續(xù)運(yùn)行基準(zhǔn)站的GNSS控制網(wǎng)
控制網(wǎng)平差采用武漢大學(xué)研制的PowerNet軟件,首先在CGCS2000框架下進(jìn)行三維無(wú)約束平差;然后進(jìn)行CGCS2000框架下的三維約束平差,獲取GNSS控制點(diǎn)的大地坐標(biāo)以及大地高;最后通過(guò)區(qū)域似大地水準(zhǔn)面精化模型進(jìn)行擬合與推算,求取對(duì)應(yīng)GNSS點(diǎn)的高程異常從而獲得各點(diǎn)的正常高。
從高程傳遞原理可知,高程傳遞精度取決于解求大地高的GNSS網(wǎng)的精度以及推算高程異常以及正常高的區(qū)域似大地水準(zhǔn)面精化模型的精度。
經(jīng)三維無(wú)約束平差的GNSS網(wǎng)最弱邊相對(duì)中誤差優(yōu)于1/400000,經(jīng)三維約束平差的GNSS網(wǎng)最弱邊相對(duì)中誤差優(yōu)于1/380000,同步環(huán)總長(zhǎng)相對(duì)中誤差優(yōu)于1/4160000,異步環(huán)總長(zhǎng)相對(duì)中誤差優(yōu)于1/190000,達(dá)到國(guó)家C級(jí)GNSS網(wǎng)精度要求;用于高程異常以及正常高擬合與推算的區(qū)域似大地水準(zhǔn)面精化模型內(nèi)符合精度±4.1cm,外符合精度±4.8cm。
為進(jìn)一步分析模型精度,采用四等幾何水準(zhǔn)聯(lián)測(cè)了陸地部分的其中5個(gè)GNSS點(diǎn),聯(lián)測(cè)高程與通過(guò)模型解算獲得的正常高比較,差值最大42mm,最小11mm,可滿(mǎn)足四等水準(zhǔn)要求,詳見(jiàn)表1。
表1 解算高程與水準(zhǔn)高程差值
(1)幾何水準(zhǔn)高差與所求得的高程差比較
對(duì)傳遞到外伶仃島的兩個(gè)GNSS點(diǎn)GC13和GC15的高差按四等水準(zhǔn)要求進(jìn)行了幾何水準(zhǔn)檢測(cè),所得的獨(dú)立高差,與本次傳遞的高差比較,相差41mm,見(jiàn)表2。
表2 傳遞高差檢測(cè)情況
(2)重合點(diǎn)高程比較
對(duì)傳遞到萬(wàn)山島的GNSS點(diǎn)GC14和2013年傳遞的水文基點(diǎn)高程[2]按三等水準(zhǔn)要求進(jìn)行幾何水準(zhǔn)引測(cè),求得的重合點(diǎn)高程,與本次傳遞的高程比較,相差46mm。
(3)用基于EGM2008模型的“移去-擬合-恢復(fù)”法檢核
“移去-擬合-恢復(fù)”法原理是在利用函數(shù)模型(如二次曲面模型)進(jìn)行高程轉(zhuǎn)換前,首先移去用地球重力場(chǎng)模型計(jì)算得到高程異常的長(zhǎng)波部分或者移去地形改正的短波部分,或者移去二者之和,然后對(duì)剩余高程異常進(jìn)行擬合和內(nèi)插,在內(nèi)插點(diǎn)上再利用重力場(chǎng)模型或地形改正公式把移去的部分恢復(fù),最終得到該點(diǎn)的高程異常得,從而可按式(1)求得未知點(diǎn)上的正常高。
表3給出了采用基于EGM2008模型的“移去-擬合-恢復(fù)”法對(duì)四個(gè)海島的5個(gè)GNSS點(diǎn)高程傳遞成果CORS同步技術(shù)傳遞方法所得的成果比較,其中,較差最大42mm、最小僅3mm。
表3 兩種傳遞方法所求的高程差值
(4)同步水位推算高程比較
通過(guò)收集三灶水文站、馬騮洲水文站、九州水文站和大萬(wàn)山島水文站近5年的水位觀(guān)測(cè)資料,利用三灶水文站、馬騮洲水文站和九州水文站三個(gè)長(zhǎng)期驗(yàn)潮站推算出大萬(wàn)山水文站水尺零點(diǎn)的高程,并用三等幾何水準(zhǔn)聯(lián)測(cè)獲得GC14的水準(zhǔn)高程,與CORS同步技術(shù)傳遞方法所得的成果比較相差36mm。
(1)CORS同步技術(shù)是目前最簡(jiǎn)便易用的高程傳遞方法。研究表明,在低高程異常且高程異常變化平緩區(qū)域,采用CORS同步技術(shù)進(jìn)行遠(yuǎn)距離跨海高程傳遞可達(dá)到四等水準(zhǔn)精度要求。
(2)研究有效解決了珠江河口區(qū)域長(zhǎng)距離海島礁高程傳遞的關(guān)鍵技術(shù)難題,獲得了伶仃洋海域桂山島、東澳島、大萬(wàn)山島、外伶仃島四個(gè)遠(yuǎn)距離海島的85基準(zhǔn)高程成果,可為萬(wàn)山島、桂山島、東澳島和外伶仃島等海島海洋水文觀(guān)測(cè)提供基于國(guó)家85高程基準(zhǔn)下的驗(yàn)潮零點(diǎn)數(shù)據(jù)。
(3)隨著GNSS連續(xù)運(yùn)行基準(zhǔn)站覆蓋區(qū)域的不斷擴(kuò)大以及區(qū)域似大地水準(zhǔn)面模型的進(jìn)一步精化,采用此方法有望獲得精度更高的遠(yuǎn)距離跨海高程傳遞成果。