復雜海域高精度海底地形測量關鍵問題研究

王雙喜，肖 強，孫雪潔

（92899 部隊，浙江 寧波 315200）

蘇北淺灘海域輻射狀沙洲群，灘闊槽深、地形獨特，屬于強潮區。影響此區域的潮流系統有東海的前進波、黃海的旋轉波和蘇北的沿岸潮流，這些不同流系在前進的過程中達到某一諧振狀態，剛好與淺灘沙脊地形配合，在潮流水道里迅速形成海水堆積，堆積速度過快，超出某一限度時，極易引起局部流速驟增、水位急漲的現象，簡稱“怪潮”[1]。“怪潮”極易引發嚴重的海難傷亡事故，對蘇北沿海地區的生產安全構成了非常嚴重的威脅[2]。“怪潮”海域的預警和研究離不開高精度海底地形數據，但在蘇北淺灘這種潮汐性質比較復雜的區域，應用常規水深測量方法無法進行有效的潮位改正，主檢比對、相鄰圖幅拼接誤差會非常明顯。進行常規水深測量時，潮位處理的精度取決于驗潮站的布設密度和分布情況。為了應對“怪潮”，原國家海洋局建立了多個海洋觀測平臺和樁基潮位站[3]。而在這種潮汐性質非常復雜的區域，這些驗潮站很難滿足高精度水深測量的要求。

在遠距離進行高精度海底地形測量研究方面，國內外相關研究主要集中在基于全球導航衛星系統（Global Navigation Satellite System，GNSS） 動態后處理技術（Post Processed Kinematic，PPK）進行潮位觀測和無驗潮水深測量上。趙建虎等[4]基于PPK技術對潮位觀測和垂直基準面轉換方法進行了研究，在基線距離100 km 范圍內，基于深度基準的全球定位系統（Global Positioning System，GPS）潮位精度優于10 cm。魏榮灝等[5]基于PPK 定位技術對杭州灣河口水下地形測量進行了實際應用，定位精度滿足實際工作需要。汪連賀[6]聯合GPS 和慣性測量單元（Inertial Measurement Unit，IMU） 姿態，得到瞬時水位的精確高程，驗潮距離100 km 內，精度優于10 cm。但在蘇北淺灘這些遠離海岸區域實施測量，面臨更多的挑戰。短時間GNSS PPK 潮位無法實現當地多年平均海面的精確傳遞，故不能直接使用GNSS 驗潮，需要在蘇北淺灘這種無依托的軟質海底設立驗潮站，采用自動驗潮儀和GNSSPPK 同步的方法進行深度基準面大地高的測量。此外，各種測量手段和方法獲取的海底地形大地高需要轉換到深度基準面上，因此還需要在蘇北淺灘“怪潮”海域建立一個無縫的垂直基準轉換模型，通過海洋垂直基準轉換獲得高精度的海底地形。

1 關鍵問題研究

1.1 軟質海底無依托驗潮站的布設

自容式壓力式驗潮儀以其輕便靈活、無須設置驗潮井等優勢在海洋測繪中得到了廣泛應用。以TGR-2050 驗潮儀為例，其觀測精度可達1 cm 以內。自動驗潮儀通過檢測海水的壓力變化來計算海面的波動。實測水下固定位置的壓力，將絕對壓力值再減去同一時刻的大氣壓值就可以得到該驗潮點的水壓值[7]。

蘇北淺灘底質大部分是屬于沙泥性質，相對比較松軟，而且遠離海岸無依托。在無依托海域設置驗潮站，針對不同的海底底質和水深情況往往有著不同的驗潮站布設方法。但需要慎重選擇拋投地點，并采取必要的保護措施，保證驗潮站的自身安全穩定和潮位數據的準確獲取。

在較淺（低潮水深在3 m 以內）的軟基海域一般布設樁基式驗潮站。長樁桿埋入海底，在樁桿底部綁定安裝自容式驗潮儀。其他海域可以用拋投式方法布設海上定點驗潮站。拋投式潮位站一般選擇在海底平坦、海流較小的海域，盡量避開航道和漁網區。拋投式驗潮站一般固定安裝在基座上，為了有效保護驗潮站基座安全，基座通過足夠長的纜繩與鐵錨連接，基座沿4 個方向展開分別設置防沉板。此外，拋投式潮位站需設置警示浮筒和網位儀，采用錨拖式方法進行回收。

1.2 GNSS 同步法確定驗潮站深度基準面大地高

高精度的動態大地高測量一般采用實時動態定位技術（Real-Time Kinematic，RTK）、PPK 或精密單點定位技術（Precise Point Positioning，PPP），RTK 需要通過電臺實時傳遞差分信息，作業距離在15 km 以內[8]。遠距離時可以采用PPK 或PPP，PPK測量的基本原理與RTK 相同，只是不需要實時傳遞差分信息，移動站只接收GNSS 信息，事后進行差分處理，作業距離可以達到100 km 甚至更遠[9]。而PPP 測量模式不受作業距離限制，但動態情況下PPP 大地高測量精度低于RTK 和PPK，特別是大地高方向的精度稍差。隨著北斗衛星定位技術的成熟，遠距離PPK 測量時，利用GPS、GLONASS（Global Navigation Satellite System）和北斗三模GNSS接收機，再應用精密衛星軌道和鐘差，就能得到更高精度的三維坐標。

實際海底地形測量應用時以參考橢球面作為測量基準獲取的數據最終應轉換成深度基準面下的水深，這需要構建深度基準面和參考橢球面的垂直基準轉換模型[10-12]。而對于每個驗潮站來說，則需要獲得其深度基準面的大地高。針對蘇北淺灘海域這種遠離岸邊的驗潮站，可以通過測量船錨定方法，在驗潮站附近用GNSS PPK 技術并綜合自動驗潮儀數據進行深度基準面大地高測量。具體過程如下。

（1）自動驗潮儀與長期站同步觀測，推算出當地多年平均海面MSL 和深度基準面L。

（2）錨定情況下通過遠距離GNSS PPK 獲取大地高h。

（3）對測船姿態進行改正，對高程信息進行濾波。

（4）GNSS PPK 測量的瞬時水位HW與自動驗潮儀潮位進行比對，計算得到該驗潮站的深度基準面大地高。

測量船錨定情況下通過姿態傳感器和GNSS 數據進行PPK 解算，獲得該驗潮站的瞬時水面信息。GNSS 同步法如圖1 所示，通過PPK 解算能夠準確測定GNSS 天線的大地高，天線到水面的天線高可以在碼頭靠泊時通過全站儀測量得到。因此瞬時水面的大地高計算公式如下。

圖1 GNSS 同步法確定深度基準面大地高

式中，HW為瞬時海面的大地高；Ght為GNSS 天線大地高；z 為GNSS 天線到海水面的高差。

與自動驗潮儀潮位數據進行對比，可以獲取深度基準面的大地高，計算公式如下。

式中，HCD為驗潮站深度基準面大地高；TCD為自動驗潮儀以深度基準面為基準的潮位。經濾波得到的GNSS 潮位信息能夠正確反映海上定點驗潮站的潮位信息，從而驗潮站就能獲得基于參考橢球面和深度基準面的兩種基面的潮位觀測值。直接比較潮位序列即可確定基面的關系，即深度基準面大地高。

式中，n 表示潮汐序列的總個數；TGNSS表示GNSS 潮位。

為了保證精度，GNSS 同步法需要連續觀測4 h以上，觀測時間至少要有一個低平潮或者高平潮過程，最好進行完整潮汐周期對比。自動驗潮儀布設期間應該進行兩次GNSS 同步觀測，然后將計算得到的深度基準面大地高進行對比，消除各種系統誤差影響。

1.3 復雜海域無縫垂直基準模型的構建

以蘇北淺灘為例，淺灘面積大，局部地形水深相對較淺，低潮時大面積干出。如果只采用水深測量方法則效率較低，只能搶在高潮時進行測量，還容易造成船只擱淺產生事故，人員安全無法保證。在大面積干出灘進行測量時，可以采用氣墊船和機載LIDAR（Light Detection And Ranging）測量手段，并輔以人工測量進行驗證。在重點航道，地形復雜的深水區域可以采用多波束無驗潮水深測量，其他區域可以高潮時段進行單波束無驗潮水深測量。這些測量手段和方法的共同點是實際獲得的測量數據為海底地形的大地高，因此需要構建連續無縫的垂直基準模型[13]，以保證每個測量點的大地高能夠轉換到深度基準面上，而且還可以進行無縫連接。

如圖2 所示，測量數據的參考基面為參考橢球面，根據平均海面與地球重力場模型EGM2008 大地水準面模型相近似的特征，可以采用移去恢復法，根據EGM2008 模型，通過驗潮站深度基準面大地高構建出復雜海域無縫垂直基準模型，具體步驟如下。

圖2 垂直基準關系

（1）先確定復雜海域包括驗潮站在內的所有測量范圍。

（2）在測量區域建立深度基準面規則格網，格網大小根據實際的潮汐性質來確定，通過各驗潮站的深度基準面L 值在網格點進行插值。

（3）測區范圍內建立一個大地水準面差距N 的格網文件。

（4） 建立海面地形的插值模型，通過GNSS PPK 測量和當地多年平均海面計算，求出當地多年平均海面的大地高HMSL，并與大地水準面差距N 進行比較，得到海面地形HTSS如下。

（5）建立整個測區的深度基準面大地高模型，深度基準面大地高HCD的計算公式如下。

2 實例分析

2.1 技術方案

實驗區域選擇蘇北淺灘海域，本區域包含各種海底地形，包括“怪潮”容易產生海域、大面積干出灘、港口和航道等。為了得到高精度的海底地形，本文通過以上關鍵問題技術研究，制定了一套針對“怪潮”海域的海底地形測量作業的技術方案。為了構建大范圍測量海域的垂直基準轉換模型，共布設6 個拋投式無依托海上驗潮站，測區已有6 個長期驗潮站。驗潮站都采用TGR-2050 自容式驗潮儀，并與長期驗潮站同步觀測1 個月。驗潮儀布設24 h 后采用錨定法進行第一次GNSS 同步觀測，驗潮儀回收前再進行第二次GNSS 同步觀測，觀測時長4 h 以上，其中包含低平潮或高平潮，長期站也進行兩次GNSS 同步觀測。自動驗潮儀回收后，所有數據都穩定有效，沒有出現明顯的沉降和位移。

2.2 精度分析

GNSS 接收機采用雙頻多模接收機，GNSS PPK測量采用雙基站模式，最遠基線長度為94 km。為了驗證遠距離GNSS PPK 解算的精度，設計了兩種驗證方案：一是用兩個基準站（相距62 km）中的一臺當流動站進行處理；二是利用三臺接收機，流動站與基站相距92 km 左右，與另一個已知控制點相距5 km 以內。控制點的坐標先通過靜態相對測量進行解算獲得毫米級的三維定位精度，以此控制點作為基準站進行短距離動態解算，解算結果認為是真值，再評估遠距離GNSS PPK 解算結果的精度。觀測數據采用Waypoint 軟件進行處理，遠距離時下載IGS（International GNSS Service） 精密星歷參與解算。兩種方案計算出的大地高對比情況如圖3 和圖4 所示。

圖3 遠距離靜態測量時動態處理結果對比

圖4 遠距離動態測量結果對比

以初始坐標為參考，評估靜態情況下的動態處理的精度，統計結果顯示大地高差值Δh 的誤差區間為[-0.08 m，0.12 m]，集中分布在5 cm 以內，中誤差為±0.032 m；以近距離解算結果為真值，評估動態模式的精度，統計結果顯示大地高差值Δh 的誤差區間為[-0.13 m，0.14 m]，集中分布在10 cm 以內，中誤差為±0.045 m。兩種情況的大地高精度都在10 cm 以內，80%以上在5 cm 以內。靜態測量時精度要明顯高于動態測量，這是由于靜態測量時的PPK 解算中的濾波算法精度更高，此外遠距離動態測量時只能與近距離測量結果進行對比，近距離PPK 測量本身也會存在一定的誤差。

采用同步法進行驗潮站深度基準面大地高推算時，GNSS PPK 采用的兩個基準站分別位于測區南北，兩站相距62 km。解算時為提高精度采用雙基站模式，以其中1 個基站3 次24 h 精密單點定位解算的平均值作為初始坐標，另外一個基準站用靜態相對測量進行解算。12 個驗潮站分別進行了兩次GNSS 同步觀測，經綜合計算后的GNSS PPK 潮位與基于深度基準面潮位進行對比，結果如表1 所示。

表1 GNSS 同步法兩次結果對比

從表1 可以看出，兩次GNSS 同步測量的結果相差在5 cm 以內，具有很好的一致性。同步時間較長的驗潮站，深度基準面大地高的精度也相對較高。此外，驗潮站站間距離對垂直基準關系的計算也有一定的影響，驗潮站與基準站相距較遠時，深度基準面大地高的精度相對較低。但站間距離在100 km 范圍內，內符合精度仍優于5 cm。

為了驗證外符合精度，在小洋口長期驗潮站架設了GNSS 天線進行靜態測量，通過長期站已知的各基面關系，直接推算出長期站的深度基準面大地高，經計算兩者僅相差1.8 cm。結果表明，直接計算的精度與GNSS 同步法基本一致，滿足相應規范精度要求。

為了驗證構建的垂直基準模型精度，小洋口等10 個驗潮站參與垂直基準模型計算，模型格網設置為2 km。利用黃沙洋和火星沙兩個驗潮站作為檢查站。經模型插值計算得到的兩個檢查站的深度基準面大地高與實際測量值進行對比，兩站分別相差3.9 cm 和2.1 cm。結果表明，無縫垂直基準模型的精度在5 cm 以內，考慮到遠距離GNSS 的測量精度，其綜合精度優于10 cm。

為了驗證最終海底地形測量成果的精度，選取測區1∶25 000 的水深成果圖，分別采用常規水位改正方法和本文方法進行水位改正處理，水深成果主檢測線交叉點不符值共檢測1 758 個點，對比如表2 所示。

從表2 可以看出，采用常規水位方法時，水深測量精度不符合規范要求，這是由于潮汐性質的差異較大，而驗潮站布設密度不足所引起。而以參考橢球為測量和數據處理的基準，應用本文方法進行水深改正時得到的最終水深成果具有較高的精度，能夠滿足在復雜海域實施海底地形測量的精度規范要求。

表2 水深成果主檢測線交叉點不符值對比分析

3 結 論

本文通過蘇北淺灘測量實踐進行研究，重點解決了復雜海域潮位改正和垂直基準轉換中存在的問題，結合實際應用得到的主要結論如下。

（1）在軟基海底布設無依托海上驗潮站時，淺水海域采用樁基式，較深海域采用拋投式能夠獲得較為穩定的潮位數據。

（2）GNSS PPK 基線距離在92 km 時，大地高的解算精度優于10 cm，能夠滿足測量規范精度要求。

（3） 通過對12 個驗潮站GNSS 同步法計算出的深度基準面大地高進行對比分析，基線距離小于100 km 時，兩次GNSS 同步測量的結果互差在5 cm以內，兩者精度具有一致性。

（4）通過構建垂直基準模型可以對整個測區水深數據進行無縫轉換，無縫垂直基準模型的精度優于5 cm，考慮到遠距離的GNSS 測量精度，綜合精度優于10 cm。

本文針對蘇北“怪潮”海域這種特殊海底地形測量進行研究分析，技術方案與精度評估結論對其他潮汐復雜海域的海底地形測量也有一定的參考價值。