趙 磊,孫洪瑞

(中交第三航務(wù)工程勘察設(shè)計院有限公司，上海 200038)

隨著中國“一帶一路”倡議的推進(jìn)，海外水運(yùn)工程項目日益增加。由于一些國家的經(jīng)濟(jì)、技術(shù)水平較低，沒有完善的大地控制網(wǎng)，平面控制點和高程控制點資料稀少，甚至個別國家連國家控制網(wǎng)都沒有。此外，一些擬建地區(qū)地質(zhì)穩(wěn)定性較差，比如西非海岸以沙灘為主，南亞地區(qū)有大面積的沖積平原，南太平洋島國經(jīng)常發(fā)生火山爆發(fā)、地震等情況，造成建成的控制點經(jīng)常發(fā)生沉降位移，給工程測量控制網(wǎng)的建立和后期維護(hù)帶來很大的困難，而且項目建設(shè)設(shè)計與建設(shè)周期較長，前后期高程基準(zhǔn)的統(tǒng)一問題亟待解決。

1 工程控制網(wǎng)布設(shè)現(xiàn)狀

1.1 有官方控制點的情況

以收集到的官方控制點為起算點，與項目新布設(shè)的工程控制點進(jìn)行聯(lián)測，獲取工程平面和高程數(shù)據(jù)。

1.2 無官方控制點的情況

平面控制點采用獨立坐標(biāo)系統(tǒng)，以GPS單點定位坐標(biāo)起算，利用GPS同步觀測，建立工程控制網(wǎng)，然后以WGS-84橢球為參考面進(jìn)行高斯投影，中央子午線根據(jù)工程位置進(jìn)行選擇，建立平面坐標(biāo)系統(tǒng)。高程是通過工程所在海域設(shè)立臨時潮位站，通過30 d潮位數(shù)據(jù)計算出平均海平面，按三等水準(zhǔn)精度把臨時潮位站高程與各GPS控制點聯(lián)測。

1.3 存在的問題

工程區(qū)域的平面和高程坐標(biāo)完全依賴于測量控制點而存在，一旦布設(shè)的控制點被毀損或發(fā)生沉降位移，這些重要的測繪資料將不復(fù)存在,后期使用很難重現(xiàn)或恢復(fù)。

平面坐標(biāo)數(shù)據(jù)還有可能通過重復(fù)GPS測量重新建立起來，高程數(shù)據(jù)將不可能復(fù)現(xiàn)，因為短期臨時潮位站測得的月平均海平面并不穩(wěn)定，會隨時間變化而變化。圖1a)、b)分別顯示2017、2018年我國東海和南海月平均海平面與常年同期月平均海平面變化情況[1]。

圖1 東海、南海沿海月平均海平面變化

2 高程起算基面的確定

隨著GPS單點精密定位技術(shù)的出現(xiàn)，通過GPS測量能提供地面點精確的三維坐標(biāo)值,其高程數(shù)據(jù)依據(jù)地球橢球面，稱之為大地高程；而海外工程偏遠(yuǎn)地區(qū)使用的高程數(shù)據(jù)，來源于前期水文測量人員建立的臨時潮位站一個月的平均海平面測量數(shù)據(jù)。為充分利用GPS所提供的大地高程數(shù)據(jù)，將工程上使用的高程數(shù)據(jù)與GPS所提供的大地高數(shù)據(jù)進(jìn)行關(guān)聯(lián)是使工程上使用的平均海平面數(shù)據(jù)得以長久保存的關(guān)鍵方法。

2.1 精密單點定位的定義

精密單點定位是利用國際GPS地球動力服務(wù)局(IGS)提供的精密星歷和精密鐘差對單臺雙頻GPS接收機(jī)所采集的相位及偽距觀測值進(jìn)行定位解算，獲得高精度控制點的ITRF坐標(biāo)的定位方法,簡稱PPP。

2.2 精密單點定位的基本原理

GPS精密單點定位一般是采用單臺雙頻GPS接收機(jī)采集的載波相位觀測值,利用IGS提供的精密星歷(事后精密星歷或快速精密星歷)和鐘差進(jìn)行的高精度定位。

精密單點定位計算的主要過程包括預(yù)處理觀測數(shù)據(jù)、精密星歷和精密衛(wèi)星鐘差擬合成軌道多項式、各項誤差的模型改正及參數(shù)估計等。

2.3 精密單點定位中的坐標(biāo)框架

精密單點定位采用IGS精密星歷，所以精密單點定位解算出的坐標(biāo)與所使用的IGS精密星歷的坐標(biāo)框架(ITRF系列)一致，與常用的WCS-84坐標(biāo)系統(tǒng)有差異，這是因為IGS精密星歷與GPS廣播星歷所對應(yīng)的參考框架不同。另外，不同時期IGS精密星歷所使用的ITRF也不同，所以在進(jìn)行精密單點定位數(shù)據(jù)處理時，需要明確所用精密星歷對應(yīng)的參考框架和歷元，并通過框架和歷元的轉(zhuǎn)換公式進(jìn)行統(tǒng)一。坐標(biāo)計算框架見圖2。

圖2 坐標(biāo)計算框架

2.4 解算軟件及解算精度分析

2.4.1解算軟件

可為用戶提供免費(fèi)在線服務(wù)的解算軟件主要包括加拿大的CSRS-PPP、GAPS-PPP和美國的 APPS-PPP等軟件。

2.4.2解算精度

夏朋飛等[2]通過處理多個IGS跟蹤站的數(shù)據(jù)對這3種精密單點定位軟件進(jìn)行了精度評估，其精度結(jié)果見表1。

表1 3種軟件計算的中誤差 cm

注：E為東向，N為北向，U為垂直方向。

徐永斌[3]于2011年用BJFS IGS站的觀測數(shù)據(jù)，對APPS和CSRS兩種在線軟件解算結(jié)果進(jìn)行分析，結(jié)果與BJFS IGS的真值進(jìn)行比較，其精度結(jié)果見表2。

表2 2種軟件計算的中誤差 cm

可見，在線精密單點定位解算軟件的解算精度可以滿足工程建設(shè)需要。

2.5 高程基準(zhǔn)

2.5.1高程系統(tǒng)

在測量中常用的高程系統(tǒng)有大地高系統(tǒng)、正高系統(tǒng)和正常高系統(tǒng)[4-5]。

1)大地高系統(tǒng)是以地球橢球面為基準(zhǔn)面的高程系統(tǒng)。大地高的定義是：由地面點沿通過該點的橢球面法線到橢球面的距離。大地高也稱為橢球高,一般用符號H表示。利用GPS單點精密定位技術(shù)可以直接測定觀測站在ITRF中的大地高。

2)正高系統(tǒng)是以大地水準(zhǔn)面為基準(zhǔn)面的高程系統(tǒng)。正高的定義是：由地面點沿通過該點的鉛垂線至大地水準(zhǔn)面的距離。正高用符號Hg表示。

3)正常高系統(tǒng)是以似大地水準(zhǔn)面為基準(zhǔn)的高程系統(tǒng)。正常高的定義是：由地面點沿通過該點的鉛垂線至似大地水準(zhǔn)面的距離,正常高用Hγ表示。

2.5.2平均海平面

水運(yùn)工程上經(jīng)常使用的平均海平面是以似大地水準(zhǔn)面為基準(zhǔn)的高程系統(tǒng)，平均海平面與大地高的差值可以視為大地異常，記為ζ。大地高與平均海平面之間的關(guān)系式：

H=Hγ+ζ

(1)

式中：H為平均海平面，Hγ為正常高，ζ為大地異常。

2.6 海外水運(yùn)工程高程基準(zhǔn)的確定

在沒有長期資料的海外國家主要通過潮位觀測確立高程基準(zhǔn)面，潮位站的基本水準(zhǔn)點埋設(shè)應(yīng)利于點位長期保存和水尺零點聯(lián)測。根據(jù)公式(1)可以得到當(dāng)?shù)氐拇蟮禺惓Ｖ郸?，從而將平均海平面歸算至ITRF下的大地高。

我國水運(yùn)工程統(tǒng)一采用理論最低潮面為深度基準(zhǔn)面。由于初勘階段時間有限，潮位觀測一般以1個月為主，推算出當(dāng)?shù)卦缕骄Ｆ矫?。再利?個月的潮位數(shù)據(jù)，通過最小二乘法進(jìn)行準(zhǔn)調(diào)和分析，得到11個調(diào)和常數(shù)，利用這11個參數(shù)可以初步計算理論最低潮面，供設(shè)計人員進(jìn)行參考使用。缺少的2 個長周期分潮Sa和SSa待日后進(jìn)行1 a以上潮位觀測獲取。

3 工程實例

3.1 巴布亞新幾內(nèi)亞科考瑞深水港可行性研究項目GPS測量

在巴布亞新幾內(nèi)亞科考瑞深水港可行性研究項目中，利用CSRS-PPP在線服務(wù)，對G002、G004點進(jìn)行了精密單點定位計算。

G002控制點在2019年5月5日使用雙頻GPS接收機(jī)靜態(tài)觀測了4 h21 min,其解算結(jié)果和精度指標(biāo)見表3。

表3 G002控制點三維坐標(biāo) ITRF14(2019)

G004控制點在2019年5月3日使用雙頻GPS接收機(jī)靜態(tài)觀測了6 h25 min,其解算結(jié)果和精度指標(biāo)見表4。

表4 G004控制點三維坐標(biāo) ITRF14(2019)

從表3和表4可知，平面和高程的解算精度均在5 cm以內(nèi)，成果精度均滿足測量規(guī)范要求。

3.2 潮位觀測

巴布亞新幾內(nèi)亞科考瑞深水港可行性研究項目中在河口處建立了T1潮位站，潮位觀測為期1個月，臨時水尺高程按三等水準(zhǔn)精度與水準(zhǔn)點G004接測。潮位觀測采用自容式小型驗潮儀，采用當(dāng)?shù)貢r間(東十區(qū))，高程系統(tǒng)為ITRF14(2019)大地高為基準(zhǔn)，每10 min記錄1次，提供整點潮位資料，其整點潮位觀測曲線見圖3。

圖3 巴布亞新幾內(nèi)亞科考瑞深水港一個月連續(xù)潮位觀測曲線

3.3 高程起算基面確定

在水準(zhǔn)點G004上架設(shè)單臺雙頻GPS 接收機(jī)，連續(xù)觀測6 h,利用PPP技術(shù)得到ITRF下三維絕對坐標(biāo)，其平面精度為2～3 cm,高程精度為2～4 cm。建立在ITRF下的平均海平面、理論最低潮面以及水尺零點高程相互關(guān)系見圖4。

圖4 單點精密定位與潮位觀測聯(lián)測

通過分析臨時潮位站的1個月連續(xù)觀測潮位，可以計算出11個調(diào)和常數(shù)[6](表5)。

進(jìn)一步推算出T1臨時潮位站在ITRF下的各基面換算關(guān)系(圖5)，從而為本項目建立了完整的高程起算面。

表5 主要分潮調(diào)和常數(shù)

圖5 巴新T1臨時潮位站基面關(guān)系(單位： m)

4 結(jié)論

1)通過利用在線PPP服務(wù)，可以獲得厘米級的坐標(biāo)和高程數(shù)據(jù)，為缺乏起算資料的地區(qū)建立平高控制網(wǎng)提供便利。

2)提出一種解決海外水運(yùn)工程高程起算面的方法，將高程控制網(wǎng)建立在ITRF下，解決了海外水運(yùn)工程前后期高程基準(zhǔn)的統(tǒng)一和復(fù)現(xiàn)問題。

3)在小范圍內(nèi)，利用精密單點定位技術(shù)測定的三維坐標(biāo)，可以代替?zhèn)鹘y(tǒng)的水準(zhǔn)測量，節(jié)約大量的人力物力。