張國堙，陶春輝，王奡，3，鄧顯明，何擁華

(1.國家海洋局第二海洋研究所，浙江 杭州 310012；2.國家海洋局海底科學重點實驗室，浙江 杭州 310012；3.中國地質(zhì)大學(武漢) 地球物理與空間信息學院，湖北 武漢 430074)

深水多波束聲吶測深數(shù)據(jù)精度評估

張國堙1，2，陶春輝1，2，王奡1，2，3，鄧顯明1，2，何擁華1，2

(1.國家海洋局第二海洋研究所，浙江 杭州 310012；2.國家海洋局海底科學重點實驗室，浙江 杭州 310012；3.中國地質(zhì)大學(武漢) 地球物理與空間信息學院，湖北 武漢 430074)

多波束測深精度評估是水深測量質(zhì)量控制的重要方面，靜態(tài)精度評估與交叉測線動態(tài)精度評估能夠從不同角度表征測深精度，估計測量樣本的綜合誤差。在實際調(diào)查作業(yè)過程中，由于缺少水深真值，在進行精度估計時缺少可操作性。本文利用Kongsberg EM120型深水多波束系統(tǒng)的測深數(shù)據(jù)，基于某一區(qū)域的重復測量數(shù)據(jù)，應用中央波束的水深數(shù)據(jù)進行靜態(tài)精度分析；通過引入網(wǎng)格化方法，進行動態(tài)水深精度評估分析，并通過偏差分析揭示測量樣本的誤差分布特征。結(jié)果表明，中央波束水深數(shù)據(jù)靜態(tài)精度評估與基于網(wǎng)格化方法的動態(tài)精度評估具有實際可操作性，其結(jié)果能夠有效估計測深的綜合誤差；重復測量數(shù)據(jù)的偏差分析能夠有效展示誤差的分布特征。

多波束聲吶；精度評估；誤差分析

1 引言

海洋水深測量是進行科學研究、資源勘查以及大陸架劃界等的基礎性工作，多波束聲吶是目前開展大面積水深測量的主要手段，也是國際海道測量組織(International Hydrographic Organization， IHO)規(guī)定的在高級別水深測量中必須采用的技術[1]。多波束聲吶獲得的水深值包含了多種因素帶來的測量誤差，其測量值只能無限逼近真值，但不等于水深真值，對測量值進行精度評估是進行水深測量質(zhì)量控制的重要手段。測量精度是表征測量誤差的指標，應用精度表征測量值與真實值之間的偏差，具體描述測深的誤差，表征水深真值的范圍，對于水深測量的質(zhì)量控制具有重要意義。

精度是對誤差的估計，多波束聲吶是由多個傳感器構(gòu)成的組合測量系統(tǒng)，測量誤差來源于系統(tǒng)的多個方面，其水深測量精度的評估包含了各方面誤差的綜合影響。目前，尚未有比多波束測深系統(tǒng)更高數(shù)量級精度的測深設備誕生，即難以使用另一種專門的更高精度的測量器具來評估多波束系統(tǒng)的水深測量精度；另外，在水下建立高精度的測量控制點比較困難，因此，對多波束聲吶測深的絕對精度評估難以實現(xiàn)[2]。國際海道測量組織為了規(guī)范測量數(shù)據(jù)具有足夠精度，根據(jù)不同的測量等級，制定了相應的精度容限要求。我國在《海洋調(diào)查規(guī)范：海底地形地貌調(diào)查》(GB/T 12763.10-2007)中，要求多波束測深作業(yè)須布置垂直于主測線的檢查測線，利用冗余數(shù)據(jù)進行數(shù)據(jù)質(zhì)量控制，具體開展數(shù)據(jù)精度評估，可利用定點法、重復測線法、交叉測線法和相鄰測區(qū)拼接的重疊區(qū)的重合點水深數(shù)據(jù)進行比對[3]。在行業(yè)標準《海洋多波束測量技術規(guī)程》(DZ/T 0292-2016)中，使用均方根誤差(Root Mean Square Error， RMSE)計算水深測量準確度，分別對定點測量、平行測線以及交叉測線的準確度計算作了相關規(guī)定[4]。目前，海洋水深測量的技術指標存在較大的模糊性，對實際測量作業(yè)沒有全面的參考意義[5]。多波束測深的誤差分析與精度評估一直受到學者們的關注，朱慶和李德仁[6]系統(tǒng)分析了多波束測深數(shù)據(jù)的誤差來源，并開展了相關的誤差改正處理模型，以保證測深系統(tǒng)必要的精度；李家彪等[7]探討了多波束測深的精度影響因素，并提出了控制精度改正的措施；吳自銀等[8]分析了不同誤差對測量精度的影響并進行了綜合校正；陳非凡[9]提出了一種多波束測深動態(tài)測量誤差評估方法，應用了相對誤差的指標開展了分析；李宜龍等[10]針對描述精度的數(shù)理統(tǒng)計指標進行了分析；吳英姿等[11]對交叉測線重疊區(qū)域的數(shù)據(jù)從不同層面進行了統(tǒng)計分析處理；吳超等[12]基于IHO規(guī)范的不確定度指標分析了多波束測深數(shù)據(jù)的質(zhì)量評估方法，討論了各分量精度對測量結(jié)果的影響；高君等[13]通過設計具體的檢定系統(tǒng)，進行了精度評估試驗；馬延霞[14]從具體理論上對精度評估進行分析，開展模擬仿真研究，討論多波束測深的相對精度和絕對精度的系統(tǒng)精度評估方法；楊廣義[15]針對基于多波束測深數(shù)據(jù)建立的海底數(shù)字地形模型，開展精度評估研究，優(yōu)選地形建模算法；唐秋華等[16]應用不同的多波束測深系統(tǒng)進行測深精度比較分析，基于均方根誤差開展內(nèi)符合與外符合精度評估研究；肖波等[17]應用均方根誤差，分別針對定點與重復測量區(qū)域進行了精度評估研究；黃辰虎等[18]針對不同的誤差源，建立了質(zhì)量評價指標；總之，國內(nèi)對多波束測深誤差的來源以及改正方法進行了深入研究并有較為完善的方法，在精度評估方面開展了多層次的探討，較多地使用均方根誤差作為精度指標進行分析研究[16-17]，但并未形成完善統(tǒng)一的方法。國外學者也使用均方根誤差表征精度，并且具體針對不同波束范圍內(nèi)的測深精度進行了評估[19]；對多波束聲吶測量的水平定位與垂直測深精度進行聯(lián)合評估[20]；較早對多波束測深的誤差來源進行了深入研究，并推導了總傳播誤差[21]；應用船舶在動力定位狀態(tài)下，保持艏向不變，采集多波束靜態(tài)水深數(shù)據(jù)，對比分析不同性能指標的運動傳感器與深度誤差的相關性[22]；通過與單波束測深數(shù)據(jù)進行比對分析，探討多波束測深的絕對偏差，并且比較邊緣波束與中央波束的測量精度[23]；基于實際測深數(shù)據(jù)進行了精度評估分析，并與IHO規(guī)范作比較研究[24]；基于IHO規(guī)范中的不確定度評估指標，在實際數(shù)據(jù)質(zhì)量評估中開展相關的應用研究[25]；對除導航定位、運動傳感器以及聲速剖面等輔助傳感器之外，針對聲學測量技術本身的誤差評估開展了研究[26]；利用多次重復測量數(shù)據(jù)，針對水平定位誤差與深度測量誤差開展了研究[27]。國際與國內(nèi)均發(fā)布了較多的多波束水深測量相關規(guī)范以及相關精度評估方法研究，但是對于多波束測量成果精度評估尚未形成一個統(tǒng)一的精度評估模型，并且，大部分測量規(guī)范強調(diào)了均方根誤差的計算，但是對于如何利用均方根誤差沒有給出進一步的詳細說明[28]。目前，通過對儀器本身的多次測量結(jié)果進行統(tǒng)計分析是評估其精度的主要手段，即相對精度評估。相對精度評估包含靜態(tài)精度評估與動態(tài)精度評估，靜態(tài)精度評估即對多波束系統(tǒng)在固定點作業(yè)的測深數(shù)據(jù)進行精度評估，動態(tài)精度評估即對多波束系統(tǒng)在動態(tài)掃測作業(yè)的測深數(shù)據(jù)進行精度評估。

多波束測深系統(tǒng)最大的特點是其一次測量能夠獲得上百個測點，多個測點對應不同位置的地形，針對定點靜態(tài)精度評估，本文應用測量船在開啟動力定位時，船舶最大限度保持靜態(tài)時所采集的多波束水深數(shù)據(jù)，通過統(tǒng)計分析，獲得靜態(tài)作業(yè)狀態(tài)下的定點測量精度；對于動態(tài)測量數(shù)據(jù)，選取同一區(qū)域交叉測線數(shù)據(jù)，應用相同參數(shù)進行網(wǎng)格化處理，獲得同一位置的兩次測量數(shù)據(jù)，進行系統(tǒng)動態(tài)精度評估。同時，基于IHO測量規(guī)范，針對動態(tài)測量數(shù)據(jù)，計算相應測量等級的誤差容限，與動態(tài)精度評估結(jié)果比較，評價測量等級。

2 多波束測深系統(tǒng)的誤差來源與測深精度評估指標

2.1 多波束測深系統(tǒng)原理

多波束測深系統(tǒng)的基本原理，是通過同時獲得的多個波束在海水中的旅行時和到達角，結(jié)合聲速剖面，應用斯涅耳法則進行聲線跟蹤，計算波束腳印在船體坐標系下的水深，再應用實時三維姿態(tài)數(shù)據(jù)和全球?qū)Ш叫l(wèi)星系統(tǒng)的定位數(shù)據(jù)，把船體坐標系下的水深點轉(zhuǎn)換到地理坐標系下，最終獲得波束腳印對應的水深及其地理坐標。多波束聲吶系統(tǒng)一次發(fā)射接收聲脈沖可以形成上百個波束，即可獲得地理坐標系下的上百個水深數(shù)據(jù)，極大地提高了水深測量的作業(yè)效率。在水深的計算過程中，多個傳感器參與其中，誤差也伴隨在整個測量過程中。具體涉及了換能器的安裝、回波信號的檢測、姿態(tài)數(shù)據(jù)的融合、聲速剖面的計算、不同坐標系之間的轉(zhuǎn)換，所有這些過程都將對測深結(jié)果引入不同程度的誤差。

2.2 多波束測深系統(tǒng)的誤差

多波束測深的誤差可以分為系統(tǒng)誤差和隨機誤差，系統(tǒng)誤差通過各種校正方法進行補償消除；對于隨機誤差，其伴隨整個測量過程，通過精度評估評定其大小，確定測量水深的偏差范圍，進而估計水深的真值。系統(tǒng)誤差與系統(tǒng)的各個單元相關，多波束系統(tǒng)的聲學換能器陣安裝誤差會對水深測量結(jié)果產(chǎn)生影響，由于實際安裝無法與設計的位置精確一致，造成實際船體坐標系與設計坐標系存在偏差，運動傳感器的安裝誤差也會對測量結(jié)果產(chǎn)生影響，針對聲學換能器和運動傳感器安裝產(chǎn)生的誤差，在多波束系統(tǒng)作業(yè)時，一般會先采用參數(shù)校準的方式，測得船體坐標系在3個方向(橫搖、縱搖以及艏搖)的偏差，通過實時補償或者數(shù)據(jù)后處理，消除聲學換能器和運動傳感器安裝帶來的誤差；潮位和換能器吃水誤差也會導致水深誤差，潮位誤差與潮位測量和計算方法有關，而換能器吃水誤差除了靜態(tài)吃水測量誤差外，還有動態(tài)吃水帶來的誤差，原則上，在水深大于200 m的海域進行水深測量時，不考慮潮位變化，不需要進行潮位改正[3]；聲速剖面是進行聲線跟蹤的重要參數(shù)，其會影響波束腳印的水平位置和測點水深，嚴重的聲速剖面誤差會引起真實海底的彎曲，聲速剖面誤差對測深精度影響最大，通過溫鹽深儀或者拋棄式溫深計進行實際聲速剖面測量，用于聲線跟蹤，獲得真實的波束腳印位置。這些誤差是相對固定的，屬于系統(tǒng)性誤差。還有一些在實時測量中由于信號檢測以及作業(yè)環(huán)境引入的隨機誤差，多波束測量系統(tǒng)通過振幅檢測與相位檢測兩種方式對回波信號進行處理，在平坦海底的條件下，相位檢測可以使得入射角計算誤差達到0.05°，在近似垂直入射情況下，振幅檢測可以使得射程計算誤差達到約1～2個脈沖波長，該誤差會隨著入射角的增加而增加，這兩項誤差均影響水深測量中波束腳印的水平位置與水深，但總體來說對多波束測量精度影響不大[29]。

多波束測深系統(tǒng)的不同類型的誤差有對應的處理方法。對于系統(tǒng)誤差，通過參數(shù)校準、吃水改正、聲速剖面校正等方法進行消除[30]，水深數(shù)據(jù)經(jīng)過對系統(tǒng)誤差處理改正后，理論上應該僅存在隨機誤差并達到較好的精度；對于隨機誤差，由于其無規(guī)律性，僅能通過對測深數(shù)據(jù)的精度評估評價其誤差大小。而對于測量中的錯誤，也稱為粗差，經(jīng)過數(shù)據(jù)后處理過程的濾波進行消除。事實上，系統(tǒng)參數(shù)校正和姿態(tài)補償并不完全，各輔助參數(shù)的測量以及上述各項改正模型也存在誤差，這些誤差對水深數(shù)據(jù)的影響呈現(xiàn)系統(tǒng)性，并沒有徹底消除，不同種類的誤差對測深數(shù)據(jù)精度的影響程度不一致，但是系統(tǒng)誤差已經(jīng)較小，而粗差也經(jīng)后處理剔除。隨機誤差受船舶和環(huán)境噪聲等外界因素影響，呈現(xiàn)無規(guī)律性，無法通過系統(tǒng)校正補償或者數(shù)據(jù)后處理改正進行消除，是精度評估的主要對象。測深數(shù)據(jù)在經(jīng)過精細處理后，已經(jīng)消除了系統(tǒng)誤差與粗差，系統(tǒng)的綜合誤差中，主要部分可以認為是隨機誤差，因此，綜合誤差近似服從高斯分布。

2.3 水深精度的評估指標

精度是對綜合測量誤差的估計，誤差是表征測量值與真實值之間的偏差，真值是客觀存在的，但無法精確測得，測量值僅能無限逼近真值，通過精度描述測量值與真值的偏離程度，從而獲得真值的范圍。通過對多波束測深數(shù)據(jù)的精度評估，可以描述測量水深與真實水深的偏差，評價測量數(shù)據(jù)的可信度。精度的度量指標一般采用偏差、標準差和均方根誤差等評估。偏差是觀測值的數(shù)學期望與真值之間的差別，反映的是觀測量的精確度；標準差是表示觀測值與其數(shù)學期望的差別，反映的是觀測值的精密度；均方根誤差是觀測值與真值之間的差別，反映的是觀測量的精確度。均方根誤差是使用較為廣泛的精度模型，它從整體意義上描述了水深測量值與真值的離散程度，其真正價值在于能提供真值可能存在的范圍[31]。國際海道測量組織的誤差估計要求具有95%的置信度，本文選擇95%置信度的均方根誤差作為精度評估指標。水深真值是一個客觀存在的值，但是其又是不可知的，通常，利用多次測量的數(shù)學期望即算術平均值近似地作為真值，測量的次數(shù)越多，數(shù)學期望就越逼近真值。

多波束測深精度評估分為絕對精度評估與相對精度評估。絕對評估即采用另外一種水深測量設備進行精度評估，測深儀是一種高精度的水深測量設備，可應用測深儀衡量多波束系統(tǒng)的測深精度。具體方法是，在同一水域同時應用多波束聲吶與測深儀進行測量，然后對兩種測量設備的結(jié)果進行偏差評估，以測深儀的數(shù)據(jù)為基準，對多波束的測深結(jié)果進行精度評估，與測深儀比對進行絕對精度評估，這種方法存在一些固有缺陷，測深儀僅僅能夠測得正下方的水深數(shù)據(jù)，僅能對多波束聲吶的正下方數(shù)據(jù)進行比對，無法對多波束的全覆蓋測量數(shù)據(jù)進行全面的評估，另一方面，測深儀與多波束聲吶的測深精度處于同一量級，難以真正進行絕對精度評估。由于尚無更高量級精度的全覆蓋測深設備，對多波束的測深數(shù)據(jù)需通過對其自身測量數(shù)據(jù)進行統(tǒng)計分析，即進行相對精度評估。相對精度評估是應用多波束系統(tǒng)本身的重復測深數(shù)據(jù)進行相對的精度評估，具體包含兩種方式，一種是定點靜態(tài)精度評估，采用船舶在定點狀態(tài)下，例如開啟船舶動力定位系統(tǒng)的情況下，對經(jīng)過后處理的測深數(shù)據(jù)進行精度評估；另一種是往返或交叉測線數(shù)據(jù)的精度評估，對交叉點的重復測深數(shù)據(jù)進行統(tǒng)計分析。

靜態(tài)精度評估是在多波束測深系統(tǒng)靜止的條件下考核其對同一位置測量深度的誤差，其本質(zhì)是一種系統(tǒng)穩(wěn)定性測試，用來評價聲吶系統(tǒng)本身的測量精度，但其無法暴露系統(tǒng)的多種誤差源；動態(tài)精度評估方法利用不同測線的多次測量數(shù)據(jù)進行統(tǒng)計分析，可以表征系統(tǒng)的綜合誤差；根據(jù)誤差理論，兩種精度評估方法分別從系統(tǒng)穩(wěn)定性和自符合性方面評估系統(tǒng)測深精度，是檢驗多波束測深系統(tǒng)精度是否符合海道測量標準的有效方法，可有效評估測量結(jié)果[11]。本文分別采用靜態(tài)精度評估和交叉測線精度評估兩種方法對Kongsberg EM120型多波束聲吶系統(tǒng)采集的測深數(shù)據(jù)進行精度評估，揭示誤差的分布與量級，表征水深真值的可能存在區(qū)間。

3 Kongsberg EM120深水多波束系統(tǒng)測深精度評估

3.1 系統(tǒng)概況與數(shù)據(jù)來源

EM120多波束聲吶是挪威Konsberg公司生產(chǎn)的全海深多波束測深系統(tǒng)。其工作主頻為12 kHz，最大ping率為5 Hz，測深范圍為20～12 000 m，最大發(fā)射扇面開角為150°，一次發(fā)射接收可以形成191個波束，波束角1°×1°，工作模式為等距作業(yè)模式或者等角模式。系統(tǒng)的標稱精度(RMS，95%置信度)為：在0°～±45°覆蓋區(qū)域內(nèi)，精度為0.2%倍水深；在±60°覆蓋區(qū)域內(nèi)，精度為0.3%倍水深；在60°～70°與-60°～-70°覆蓋區(qū)域內(nèi)，精度為0.5%倍水深[32]。

本文的數(shù)據(jù)來源于中國大洋調(diào)查航次，在多波束測量作業(yè)時，工作模式為等距模式，扇面開角為120°，因此本文中數(shù)據(jù)對應的系統(tǒng)標稱精度為0.3%倍水深(95%置信度)。在系統(tǒng)作業(yè)前，進行了參數(shù)校正工作，把獲得的橫搖偏差、縱搖偏差以及艏搖偏差引入系統(tǒng)的測深計算中，消除了系統(tǒng)安裝誤差，由于系統(tǒng)采用了秒脈沖信號同步，無需進行延時校準。在獲得多波束原始數(shù)據(jù)后，采用實測的聲速剖面數(shù)據(jù)，對數(shù)據(jù)進行精細后處理，剔除了粗差，獲得有效的水深數(shù)據(jù)，本文進行精度評估的對象即是已經(jīng)過系統(tǒng)安裝誤差改正和粗差濾波，主要包含隨機誤差的水深數(shù)據(jù)，其誤差近似服從高斯分布。

3.2 靜態(tài)精度評估方法

靜態(tài)精度評估是對系統(tǒng)在定點作業(yè)狀況下，采集的數(shù)據(jù)進行統(tǒng)計分析。由于海洋環(huán)境的影響，船舶無法保持絕對靜止狀態(tài)，本文采用2014年調(diào)查船在南海作業(yè)時，船舶在開啟動力定位情況下，對船舶在定點作業(yè)時采集的數(shù)據(jù)進行分析。船舶作業(yè)時，動力定位的定位點選取為船尾部，因此，船艏并非靜態(tài)，而是持續(xù)往復偏轉(zhuǎn)，多波束系統(tǒng)的波束穩(wěn)定單元為保持波束腳印與航跡垂直，會依據(jù)艏向參數(shù)進行波束偏轉(zhuǎn)，在靜態(tài)作業(yè)時，這樣就造成中央波束腳印位置保持近似不變，而邊緣波束腳印旋轉(zhuǎn)(圖1a)，邊緣波束的測量對象在不斷變化，不具備統(tǒng)計分析的條件，中央波束近似保持指向海底同一位置，故本文提取中央波束數(shù)據(jù)進行統(tǒng)計分析。

船舶保持在固定位置連續(xù)測量，對中央波束的水深測量結(jié)果進行分析統(tǒng)計。在誤差理論中，符合同一測量條件下，對同一對象多次重復測量的誤差統(tǒng)計條件，可以采用常規(guī)誤差統(tǒng)計分析理論。本文選用1 h內(nèi)采集的多波束數(shù)據(jù)，共發(fā)射與接收1 517ping，即采集獲得1 517個中央波束數(shù)據(jù)(圖1b)，采用均方根誤差對該中央波束數(shù)據(jù)進行精度評估。

(1)

式中，σ為均方根差；di為每一ping測得的水深值；d為所有水深值的平均值即數(shù)學期望，用來指代水深真值；n為測深數(shù)據(jù)的個數(shù)。經(jīng)計算，2σ=5.528 0(95%置信度)。

如果用相對誤差δ來表示，則有

(2)

得到，δ=0.46%。

系統(tǒng)在中央鏡像區(qū)域的標稱精度為0.2%倍水深，該靜態(tài)數(shù)據(jù)的樣本統(tǒng)計結(jié)果顯示，其實際探測精度低于系統(tǒng)標稱精度。

3.3 交叉測線精度評估方法

交叉測線對海底同一區(qū)域進行不同方向的測量，對兩次測量的水深數(shù)據(jù)進行統(tǒng)計分析，可以描述系統(tǒng)的自符合性，獲得系統(tǒng)全覆蓋測量的測深精度。本文選取的數(shù)據(jù)為南北向與東西向的兩條交叉測線，對重疊區(qū)域的測深數(shù)據(jù)進行精度評估(圖2)。

圖1 1小時各個波束靜態(tài)數(shù)據(jù)(a)和中央波束數(shù)據(jù)(共1 517個點)(b)Fig.1 Beam bathymetry data in one hour (a)，and central beam bathymetry data in one hour (b)

圖2 交叉測線區(qū)域地形圖(紅框內(nèi)為評估數(shù)據(jù))Fig.2 Topographic map of the survey line (evaluation data in the red box)

圖3 選定重疊區(qū)域的水深點(a)和經(jīng)過網(wǎng)格化處理的水深數(shù)據(jù)點(b)Fig.3 Bathymetry data of the study area(a)， and grid bathymetry data of the study area(b)

圖4 往返測量重疊區(qū)域的水深數(shù)據(jù)偏差(白色線為兩條交叉測線)Fig.4 Difference between the two survey line (the two cross survey lines in white colour)

測區(qū)海底較為平坦，對該區(qū)域的交叉測線數(shù)據(jù)進行精細后處理，消除了粗差。為了便于與系統(tǒng)指標進行對比，選取發(fā)射扇面開角-60°～60°區(qū)域數(shù)據(jù)進行統(tǒng)計分析，重疊區(qū)域的水深點見圖3a，其中，紅色點為南北向測線數(shù)據(jù)，黑色點為東西向測線數(shù)據(jù)。由于兩條測線的波束腳印地理位置并不一致，為了使得數(shù)據(jù)統(tǒng)計分析結(jié)果可靠，對原始的水深數(shù)據(jù)進行網(wǎng)格化處理，獲得同一地理位置的兩次測量水深值。原始數(shù)據(jù)中，波束腳印之間的距離約為100 m，采用100 m網(wǎng)格間距對數(shù)據(jù)網(wǎng)格化，網(wǎng)格化后的數(shù)據(jù)見圖3b。

利用均方根誤差進行計算：

(3)

式中，σ為均方根差；di為兩次測深值的偏差；n為一條測線的測深數(shù)據(jù)個數(shù)。經(jīng)計算，得2σ=35.974(95%置信度)，相對誤差δ=0.84%。

計算的樣本精度低于系統(tǒng)的標稱精度。為了直觀表現(xiàn)數(shù)據(jù)誤差的分布，對兩條測深數(shù)據(jù)的偏差進行分析(圖4)，結(jié)果表明，在中央波束兩側(cè)，數(shù)據(jù)偏差較大，邊緣探測誤差應該是精度降低的主要原因。該誤差一直出現(xiàn)在探測結(jié)果中，呈現(xiàn)系統(tǒng)性，在數(shù)據(jù)后處理中，基于地形相關性認為該水深數(shù)據(jù)具有合理性，不認為其是粗差，數(shù)據(jù)后處理時無法消除。通過對交叉數(shù)據(jù)進行處理，重疊區(qū)域的水深偏差可以展現(xiàn)出誤差分布，有助于判讀誤差來源。

4 結(jié)果

系統(tǒng)標稱的精度，是系統(tǒng)在理想的作業(yè)工況條件下才可達到的，是系統(tǒng)精度指標的上限。在實際作業(yè)環(huán)境中，由于各方面的限制，測深精度一般難以達到系統(tǒng)的標稱精度，本文的樣本數(shù)據(jù)精度略低于系統(tǒng)標稱精度。國際海道測量組織為了適應當今測深技術的發(fā)展，更加有效的描述不同測量目的的適用測深精度，對精度要求制定了國際化標準，并給出了不同測量等級對應的精度要求，即規(guī)定了不同測量等級對應要求的精度下限。

IHO S-44標準對各個測量等級的精度要求提出了一套統(tǒng)計指標，即各測量等級的精度容限為：

(4)

式中，a為深度常差，即所有常量誤差的總和；b×d為深度相關誤差，即所有深度相關誤差的總和；b為深度相關誤差因子；d為深度。針對本文中的測量區(qū)

域，各個等級的精度要求(95%置信度)見表1。本文中的d為交叉測線重疊區(qū)域深度數(shù)據(jù)的平均值。

表1 IHO海道測量規(guī)范測量最低標準概要

結(jié)果表明，本文的樣本數(shù)據(jù)精度優(yōu)于IHO-S44要求的一等測量標準的精度容限，測量結(jié)果滿足一等測量標準。

5 結(jié)語

通過系統(tǒng)校準，可以有效消除系統(tǒng)安裝誤差；而通過精度評估與偏差分析，能夠估計隨機誤差的量值，并揭示隨機誤差分布，系統(tǒng)校準與精度評估都是為獲得水深真值服務，互為補充。在多波束水深測量作業(yè)時，從系統(tǒng)參數(shù)校準和精度評估兩個方面對多波束測深質(zhì)量進行控制，可以保證測深的可靠性，提高數(shù)據(jù)質(zhì)量。靜態(tài)精度評估和交叉測線動態(tài)精度評估為多波束測深數(shù)據(jù)的誤差估計提供了有效的表征方法，靜態(tài)精度評估可以表征測量樣本的穩(wěn)定性，而針對交叉測線的測深重疊區(qū)域，基于網(wǎng)格化方法進行動態(tài)精度評估具有實際可操作性，其結(jié)果可以表征測量樣本的自符合性。交叉測線的重疊區(qū)域的數(shù)據(jù)偏差統(tǒng)計分析可以具體的展現(xiàn)誤差量值與分布特征，有利于進行針對性的測深質(zhì)量控制。

[1] International Hydrographic Organization. IHO Standards for Hydrographic Surveys Special Publication N44[S]. 2008．

[2] Dix M, Abd-Elrahman A, Dewitt B, et al. Accuracy evaluation of terrestrial LiDAR and multi-beam sonar systems mounted on a survey vessel[J]. Journal of Surveying Engineering, 2011, 138(4): 203-213．

[3] 國家質(zhì)量技術監(jiān)督局. 中華人民共和國國家標準GB 17378.2-2007 海洋調(diào)查規(guī)范第10部分:海底地形地貌調(diào)查[S]. 2007．

The State Bureau of Quality and Technical Supervision. GB 17378.2-2007 Specification for oceanographic survey—Part 10 Submarine topography and geomorphology[S]. 2007．

[4] 國土資源部. DZ/T 0292-2016 海洋多波束水深測量技術規(guī)程[S]. 2016．

Ministry of Land and Resources of the Peoples’s Republic of China. DZ/T 0292-2016 Technical regulations for application of marine multi-beam bathymetric survey[S]. 2016.

[5] 梁志誠, 黃文騫, 魏帥, 等. 多波束測量的精度控制與規(guī)范指標[J]. 測繪工程, 2016, 25(2): 57-60．

Liang Zhicheng, Huang Wenqian, Wei Shuai, et al. Precision control and standard indexes of multi-beam surveying[J]. Engineering of Surveying and Mapping, 2016, 25(2): 57-60．

[6] 朱慶, 李德仁. 多波束測深數(shù)據(jù)的誤差分析與處理[J]. 武漢測繪科技大學學報, 1998, 21(1): 1-4．

Zhu Qing, Li Deren. Error analysis and processing of multi-beam soundings[J]. Journal of Wuhan Technical University of Surveying and Mapping, 1998, 21(1): 1-4．

[7] 李家彪, 鄭玉龍, 王小波, 等. 多波束測深及影響精度的主要因素[J]. 海洋測繪, 2001(1): 26-32．

Li Jiabiao, Zheng Yulong, Wang Xiaobo, et al. Multi-beam sounding and the main factors affecting the depth accuracy[J]. Hydrographic Surveying and Charting, 2001(1):26-32．

[8] 吳自銀, 金翔龍, 鄭玉龍, 等. 多波束測深邊緣波束誤差的綜合校正[J]. 海洋學報, 2005, 27(4): 88-94．

Wu Ziyin, Jin Xianglong, Zheng Yulong, et al. Integrated error correction of multi-beam marginal sounding beam[J]. Haiyang Xuebao, 2005, 27(4): 88-94．

[9] 陳非凡. 多波束條帶測深儀的動態(tài)測量誤差評估[J]. 海洋技術, 1999, 18(1): 42-45．

Chen Feifan. Estimation of dynamic measuring error for multi-beam swath bathymeter[J]. Ocean Technology, 1999, 18(1): 42-45．

[10] 李宜龍, 殷曉冬, 張立華, 等. 交叉檢查在多波束測深精度評估和誤差分析中的應用[J]. 海洋技術, 2006, 25(1): 119-123．

Li Yilong, Yin Xiaodong, Zhang Lihua, et al. The application of crosscheck in multi-beam accuracy evaluation and error analysis[J]. Ocean Technology, 2006, 25(1):119-123．

[11] 吳英姿, 徐新盛, 喬力爭. 多波束測深系統(tǒng)的精度評估方法研究[J]. 海洋技術, 2003, 22(3):65-69．

Wu Yingzi, Xu Xinsheng, Qiao Lizheng. The evaluation method of the precision of the multi-beam bathymetric system[J]. Ocean Technology, 2003, 22(3):65-69．

[12] 吳超, 殷曉冬, 張立華, 等. 基于不確定度的多波束測深數(shù)據(jù)質(zhì)量評估方法[J]. 海洋測繪, 2009, 29(5):11-14．

Wu Chao, Yin Xiaodong, Zhang Lihua, et al. The method of quality estimation in multi-beam sounding data based on uncertainy[J]. Hydrographic Surveying and Charting, 2009, 29(5):11-14．

[13] 高君, 肖付民, 裴文斌, 等. 多波束測深精度檢定方法[J]. 測繪科學技術學報, 2012, 29(3): 179-182．

Gao Jun, Xiao Fumin, Pei Wenbin, et al. Calibrating the depth accuracy of multi-beam sounder[J]. Journal of Geomatics Science and Technology, 2012, 29(3):179-182．

[14] 馬延霞. 多波束條帶測深系統(tǒng)測深精度評估方法研究[D]. 哈爾濱: 哈爾濱工程大學, 2007．

Ma Yanxia. A study on surveying accuracy of estimation in multi-beam sounding system[D]. Harbin: Harbin Engineering University, 2007．

[15] 楊廣義. 基于多波束測深數(shù)據(jù)的海底地形建模技術及精度評估研究[D]. 鄭州：解放軍信息工程大學, 2010．

Yang Guangyi. Research on benthal terrain modeling technology and precision evaluaion base on multi-beam data[D]. Zhengzhou: PLA Information Engineering University, 2010．

[16] 唐秋華, 陳義蘭, 路波, 等. EM1002S與GeoSwath多波束聲吶系統(tǒng)測深精度比較分析[J]. 海岸工程, 2013, 32(4): 56-64．

Tang Qiuhua, Chen Yilan, Lu Bo, et al. Comparison of sounding accuracy between multi-beam sonar systems EM1002S and GeoSwath[J]. Coastal Engineering, 2013, 23(4):56-64．

[17] 肖波, 盛堰, 劉勝旋. 多波束測深質(zhì)量評價方法分析[J]. 海洋測繪, 2013, 33(3): 56-58．

Xiao Bo, Sheng Yan, Liu Shengxuan, et al. Method for multi-beam echosounding quality evaluation[J]. Hydrographic Surveying and Charting, 2013, 33(3):56-58．

[18] 黃辰虎, 陸秀平, 歐陽永忠, 等. 多波束水深測量誤差源分析與成果質(zhì)量評定[J]. 海洋測繪, 2014, 34(2): 1-6．

Huang Chenhu, Lu Xiuping, Ouyang Yongzhong, et al. Analysis of error source and quality assessment about multi-beam sounding product[J]. Hydrographic Surveying and Charting, 2014, 34(2):1-6．

[19] Maleika W. Development of a method for the estimation of multibeam echosounder measurement accuracy[J]. Przeglad Elektrotechniczny, 2012, 88(10):205-208.

[20] Ernstsen V B, Noormets R, Hebbeln D, et al. Precision of high-resolution multi-beam echo sounding coupled with high-accuracy positioning in a shallow water coastal environment[J]. Geo-Marine Letters, 2006, 26(3): 141-149．

[21] Hare R. Depth and position error budgets for mulitbeam echosounding[J]. The International Hydrographic Review, 1995, 72(2):37-69．

[22] de Moustier C. Field evaluation of sounding accuracy in deep water multi-beam swath bathymetry[C]// Oceans 94 Oceans Engineering for Today’s Technology and Tomorrow’s Preservation’ Proceedings, 2001, 3：1761-1765．

[23] Tonchia H. Evaluation of a deep water wide swath echo sounder for hydrographic surveys[C]//Oceans 94 Oceans Engineering for Today’s Technology and Tomorrow’s Preservation’ Proceedings, 1994, 3：225-230．

[24] Eeg J. On the estimation of standard deviations in multi-beam soundings[J]. The International Hydrographic Review, 1996, 73(2):39-51．

[25] Hare R, Eakins B, Amante C. Modelling bathymetric uncertainty[J]. The International Hydrographic Review, 2011,12(6):31-42．

[26] Lurton X, Augustin J. A measurement quality factor for swath bathymetry sounders[J]. IEEE Journal of Oceanic Engineering, 2010, 35(4): 852-862．

[27] Ernstsen V B, Noormets R, Hebbeln D, et al. Precision of high-resolution multi-beam echo sounding coupled with high-accuracy positioning in a shallow water coastal environment[J]. Geo-Marine Letters, 2006, 26(3): 141-149．

[28] 劉勝旋, 張瑤, 馬金鳳, 等. 多波束測深成果精度評估方法探討[J]. 海洋測繪, 2016, 36(5): 36-39．

Liu Shengxuan, Zhang Yao, Ma Jinfeng, et al. Accuracy assessment of multi-beam echosounding data[J]. Hydrographic Surveying and Charting, 2016, 36(5):36-39.

[29] 李家彪. 多波束勘測原理技術與方法[M]. 北京: 海洋出版社, 1999．

Li Jiabiao. Multi-beam Sounding Principles Survey Technologies and Data Processing Methods[M]. Beijing: China Ocean Press, 1999．

[30] 劉勝旋, 關永賢. 多波束系統(tǒng)的參數(shù)誤差判斷及校正[J]. 海洋測繪, 2002, 22(1): 33-37．

Liu Shengxuan, Guan Yongxian, The determination and calibration for multi-beam system parameters[J]. Hydrographic Surveying and Charting, 2002, 22(1):33-37．

[31] 張云端, 祿豐年. 數(shù)字高程模型DEM精度研究[J]. 測繪與空間地理信息, 2007, 30(3): 120-123．

Zhang Yunduan, Lu Fengnian. Research on precision of digital elevation models[J]. Geomatics & Spatial Information Technology, 2007, 30(3):120-123．

[32] Kongsberg. EM120 multi-beam echo sounder product description[R]. Kongsberg Maritime AS, Norway, 2005.

Accuracy evaluation of multibeam echo sounder bathymetry data

Zhang Guoyin1，2， Tao Chunhui1，2， Wang Ao1，2，3， Deng Xianming1，2， He Yonghua1，2

(1.SecondInstituteofOceanography，StateOceanicAdministration，Hangzhou310012，China; 2.KeyLaboratorySubmarineGeosciences，StateOceanicAdministration，Hangzhou310012，China; 3.InstituteofGeophysicsandGeomatics，ChinaUniverstiyofGeosciences，Wuhan430074，China)

Bathymetry accuracy evaluation of multibeam echo sounder is an important aspect of water depth measurement quality control. It is an effective method using the static data accuracy evaluation and cross line coverage data accuracy evaluation to estimate the depth samples error. In the actual investigation of the water depth， there is no data of the true water depth， so it is difficult to make statistics of the depth error. The depth data collected by Kongsberg EM120 is used to analysis， the static measurement statistics is calculated using the central beam data， then a grid method is applied to make statistics using the cross line coverage depth data. In addition， the distribution of the sample data total error is showed by the cross line coverage depth deviation. In conclusion， it is an effective way to using the statistics central beam data and the cross line coverage depth data for the water depth accuracy evaluation of multibeam echo sounder．

multibeam echo sounder; accuracy evaluation; error analysis

P733.23

A

0253-4193(2017)11-0106-09

張國堙， 陶春輝， 王奡， 等. 深水多波束聲吶測深數(shù)據(jù)精度評估[J]. 海洋學報， 2017， 39(11):106-114，

10.3969/j.issn.0253-4193.2017.11.010

Zhang Guoyin， Tao Chunhui， Wang Ao， et al. Accuracy evaluation of multibeam echo sounder bathymetry data[J]. Haiyang Xuebao， 2017， 39(11):106-114， doi：10.3969/j.issn.0253-4193.2017.11.010

2016-10-27；

2017-06-26。

浙江省自然科學基金項目(LQ16D060008)；國家科技支撐計劃項目(2014BAB14B01)。

張國堙(1984—)，男，山西省朔州市人，工程師，從事海底聲學探測技術及應用研究。E-mail：zgysir@126.com