蔣建平，朱漢華

（長(zhǎng)江水利委員會(huì)水文局長(zhǎng)江下游水文水資源勘測(cè)局，江蘇 南京 210011）

0 引言

聲學(xué)多普勒流速剖面儀（ADCP）具有不擾動(dòng)流場(chǎng)、連續(xù)施測(cè)、測(cè)流歷時(shí)短的特點(diǎn)[1]，目前在水文測(cè)驗(yàn)中已普遍應(yīng)用，但設(shè)備主要依賴進(jìn)口，存在價(jià)格高、采購(gòu)周期長(zhǎng)、維修困難等特點(diǎn)。國(guó)內(nèi)近些年在 ADCP 的研發(fā)和產(chǎn)品化方面取得了一些進(jìn)展，推出了一系列國(guó)產(chǎn) ADCP 產(chǎn)品，2016 年 9 月RIV-600 型 ADCP 通過(guò)了水利部新產(chǎn)品鑒定。為了推進(jìn)國(guó)產(chǎn) ADCP 應(yīng)用進(jìn)程，針對(duì)國(guó)外同頻率的ADCP 底跟蹤施測(cè)流量出現(xiàn)偏小現(xiàn)象，如流量在 70 000 m3/s 時(shí)底跟蹤施測(cè)會(huì)偏小近 20%，以及鐵質(zhì)測(cè)船導(dǎo)致 ADCP 內(nèi)部羅經(jīng)出現(xiàn)偏差，影響到流量計(jì)算的準(zhǔn)確性等技術(shù)問(wèn)題，長(zhǎng)江水利委員會(huì)水文局開(kāi)展了國(guó)產(chǎn) ADCP 多傳感器應(yīng)用的適應(yīng)性測(cè)試。

適應(yīng)性測(cè)試采用一種新的方法，分別將國(guó)產(chǎn)RIV ADCP 和同頻率的進(jìn)口 RDI ADCP（作為標(biāo)準(zhǔn)）安裝在測(cè)船同一側(cè)，均外接 GNSS（全球?qū)Ш叫l(wèi)星系統(tǒng)）和羅經(jīng)，并同步進(jìn)行數(shù)據(jù)的采集，分析國(guó)產(chǎn)RIV ADCP 與進(jìn)口 RDI ADCP 之間的誤差及相關(guān)性，檢驗(yàn)國(guó)產(chǎn) ADCP 在長(zhǎng)江應(yīng)用的適應(yīng)性及流量等測(cè)量數(shù)據(jù)的精度。

1 測(cè)試

1.1 流量測(cè)試設(shè)備

流量測(cè)試標(biāo)準(zhǔn)采用 RDI WH600 型 ADCP，因該設(shè)備在應(yīng)用于流量基本資料收集前作了大量的比測(cè)，精度滿足規(guī)范要求。國(guó)產(chǎn) ADCP 設(shè)備為RIV-600 型 ADCP，主要參數(shù)如表1 所示。考慮到測(cè)試用船均為鐵質(zhì)，斷面流量施測(cè)時(shí)有“動(dòng)底”的存在，RDI ADCP 配備了 Trimble SPS461 GNSS 一體化羅經(jīng)，并開(kāi)通了 OmniSTAR HP 星站差分服務(wù)，水平位置偏差小于 0.15 m。2 套衛(wèi)星天線安裝間距為 1.0 m，羅經(jīng)定向精度優(yōu)于 0.25°。RIV-600 型 ADCP 采用的外接定位和定向數(shù)據(jù)與 RDI ADCP 相同，分頻接入。

1.2 流量施測(cè)及采用模式選擇

流量適應(yīng)性采用走航式測(cè)試，即測(cè)船沿測(cè)流斷面往返橫渡，通過(guò)施測(cè)各點(diǎn)的流速、水深和測(cè)船的移動(dòng)方式測(cè)定流量，每測(cè)次包含 2 個(gè)來(lái)回。長(zhǎng)江干流斷面一般在中高水都存在“動(dòng)底”現(xiàn)象，且施測(cè)均采用鐵質(zhì)船作為平臺(tái)，為確保流量施測(cè)的精度，外接的 GNSS 和外部羅經(jīng)施測(cè)前進(jìn)行了校正[2-3]，并分別施測(cè)了底跟蹤和 GNSS GGA 跟蹤模式下不同的流量，通過(guò)計(jì)算統(tǒng)計(jì)用于國(guó)產(chǎn) ADCP 不同模式的適應(yīng)性，以及流量、面積、最大流速和水深等特征值的精度評(píng)估。當(dāng) GGA 跟蹤模式流量值大于底跟蹤流量值 1% 時(shí)，認(rèn)為有“動(dòng)底”現(xiàn)象，采用 GGA 模式流量值，反之則采用底跟蹤模式值。

1.3 適應(yīng)性測(cè)試

RIV-600 型 ADCP 的適應(yīng)性測(cè)試主要在長(zhǎng)江下游干流水文站進(jìn)行，主要測(cè)試 RIV-600 型 ADCP 在有“動(dòng)底”情況下多傳感器接入后的流速、流向及流量等的精度。適應(yīng)性測(cè)試進(jìn)行了 38 個(gè)流量測(cè)次測(cè)試，測(cè)試流量范圍為 12 486～70 443 m3/s，測(cè)試期間最大斷面流速范圍為 0.98～2.92 m/s，斷面平均含沙量范圍為 0.026～0.339 kg/m3，泥沙中數(shù)粒徑范圍為 0.009～0.063 mm。為了驗(yàn)證適應(yīng)性，以大通站 2017 年 7 月 6 日高水流量（斷面平均含沙量為 0.339 kg/m3）為例進(jìn)行分析。

大通水文站的一次 RIV 與 RDI ADCP 流量比測(cè)，斷面施測(cè) 2 個(gè)測(cè)回，比測(cè)統(tǒng)計(jì)如表2 所示。該測(cè)次 RDI ADCP 斷面平均施測(cè)了 471 呯數(shù)據(jù)（每呯 1.54 s）；RIV ADCP 則施測(cè)了 1 195 呯數(shù)據(jù)（每呯 0.60 s），數(shù)據(jù)量是 RDI ADCP 的 2.5 倍。RDI ADCP 施測(cè)底跟蹤模式流量值為 56 733 m3/s，GGA 跟蹤模式流量值為 70 342 m3/s，底跟蹤模式流量偏小 13 609 m3/s，相對(duì)偏小 19.3%，說(shuō)明該斷面有明顯的“動(dòng)底”現(xiàn)象。RIV 與 RDI ADCP 底跟蹤和 GGA（VTG）跟蹤模式航跡線比較圖如圖1 所示，從圖1 可以看出，底跟蹤和 GGA（VTG）跟蹤模式的航跡線明顯不重合，而且底跟蹤航跡線明顯偏離在斷面上游，說(shuō)明斷面有“動(dòng)底”現(xiàn)象，底跟蹤模式下流量會(huì)偏小。而圖1 a 和 b 坐標(biāo)和圖形的差異，主要是 2 個(gè)軟件顯示的設(shè)置參數(shù)導(dǎo)致。在底跟蹤模式下，RIV 與 RDI ADCP 流量最大相對(duì)誤差為 1.72%，平均相對(duì)誤差為 0.43%； 在 GGA 跟蹤模式下，兩者最大相對(duì)誤差為 -0.61%，平均相對(duì)誤差為 -0.28%，2 種模式下平均流量相差均小于 ±1%。

表2 RIV 與 RDI ADCP 流量比測(cè)統(tǒng)計(jì)表

圖1 RIV 與 RDI ADCP 底跟蹤和 GGA (VTG) 模式航跡線比較圖

對(duì)測(cè)次施測(cè)的斷面流向、面積及平均流速統(tǒng)計(jì)分析其相關(guān)性，比測(cè)統(tǒng)計(jì)如表3 所示。在 GGA 模式下 RIV 和 RDI ADCP 施測(cè)的斷面流向半測(cè)回最大差值的絕對(duì)值為 0.45°，測(cè)次平均流向差值的絕對(duì)值為 0.21°；施測(cè)的斷面面積半測(cè)回最大相對(duì)誤差絕對(duì)值為 0.46%，測(cè)次平均相對(duì)誤差絕對(duì)值為 0.25%；施測(cè)的斷面平均流速半測(cè)回最大相對(duì)誤差為 0.22%，平均相對(duì)誤差絕對(duì)值為 0.03%，說(shuō)明 2 種設(shè)備施測(cè)的結(jié)果基本一致。

表3 RIV 與 RDI ADCP 斷面流向、面積和平均流速比測(cè)統(tǒng)計(jì)表

再檢查分析底跟蹤和 GGA 模式下斷面流速等值圖， RIV 與 RDI ADCP 底跟蹤和 GGA 模式下斷面流速等值比較圖如圖2，3 所示。從第 1 半測(cè)回?cái)?shù)據(jù)回放可以看出，RDI ADCP 在施測(cè)時(shí)中間有幾處底跟蹤失效，而 RIV ADCP 相對(duì)數(shù)據(jù)齊全。但 ADCP 底跟蹤模式下，RDI ADCP 施測(cè)的斷面水深存在 4 個(gè)波束施測(cè)相差不一致，導(dǎo)致水深數(shù)據(jù)有鋸齒現(xiàn)象；而 RIV ADCP 施測(cè)的水深變化過(guò)程平緩光滑。GGA 模式下 2 款 ADCP 斷面流速都漸變均勻，說(shuō)明軟件中具有外接設(shè)備校正和解決“動(dòng)底”導(dǎo)致的流量偏小現(xiàn)象功能。RIV 與 RDI ADCP 底跟蹤和 GGA 模式下斷面流速矢量線比較分別如圖4，5 所示，從圖4 和 5 中可以看出，RIV 與 RDI ADCP 底跟蹤速度下流速有一點(diǎn)差異，RDI ADCP 數(shù)據(jù)的缺失主要是底跟蹤有幾處失效導(dǎo)致，但 GGA 模式下施測(cè)的斷面平均流速過(guò)程基本一致，且 RIV ADCP 數(shù)據(jù)明顯密集，這說(shuō)明施測(cè)垂線的數(shù)據(jù)明顯多于 RDI ADCP。

圖2 RIV 與 RDI ADCP 底跟蹤模式下斷面流速等值比較圖

圖3 RIV 與 RDI ADCP GGA 模式下斷面流速等值比較圖

圖4 RIV 與 RDI ADCP 底跟蹤模式下斷面流速矢量線比較圖

圖5 RIV 與 RDI ADCP GGA 跟蹤模式下斷面流速矢量線比較圖

2 斷面流量精度分析

RIV 與 RDI ADCP 共進(jìn)行 38 測(cè)次流量比測(cè)（含底跟蹤和 GGA 2 種模式），流量比測(cè)比較圖如圖6 所示，流量相對(duì)誤差統(tǒng)計(jì)圖如圖7 所示。比測(cè)流量范圍在 12 486～70 443 m3/s 之間，單次流量比測(cè)相對(duì)誤差范圍為 -2.28%～2.86%，系統(tǒng)誤差為 -0.83%，總隨機(jī)不確定度為 2.28%。通過(guò)分析，流量測(cè)驗(yàn)單次流量允許誤差滿足一類站中水位級(jí)總隨機(jī)不確定度小于 6%，系統(tǒng)誤差在 -2%～1% 之間的要求，符合規(guī)范[4]的要求。

3 結(jié)語(yǔ)

通過(guò)斷面流量和流向、面積等特征值比測(cè)，得出 RIV-600 型 ADCP 主要結(jié)論如下：RIV-600 型ADCP 在長(zhǎng)江下游干流斷面流量 70 443 m3/s 以下（斷面平均含沙量為 0.339 kg/m3）出現(xiàn)“動(dòng)底”后，具有與 RDI WH600 型 ADCP 解決流量偏小相同的功能，外接的 GNSS 和羅經(jīng)通過(guò)校正后施測(cè)的流量、流向、斷面和平均流速等特征值數(shù)據(jù)精度與進(jìn)口 RDI WH600 型 ADCP 相當(dāng)，滿足規(guī)范精度要求。RIV-600 型 ADCP 數(shù)據(jù)更新率明顯高于 RDI ADCP 設(shè)備，但沒(méi)有淺水高精度測(cè)量模式，軟件圖表輸出及顯示功能等存在不足和不規(guī)范現(xiàn)象，還需進(jìn)一步完善?？傊?，RIV-600 型 ADCP 主要性能和功能滿足長(zhǎng)江下游流量在 70 443 m3/s（低含沙量）以下測(cè)驗(yàn)，可替代進(jìn)口 ADCP 用于流量測(cè)驗(yàn)基本資料收集。

圖6 RIV 與 RDI ADCP 流量比測(cè)比較圖

圖7 RIV 與 RDI ADCP 流量相對(duì)誤差統(tǒng)計(jì)圖