南寧站定點(diǎn)式ADCP指標(biāo)流速與斷面平均流速關(guān)系模型的建立和應(yīng)用

韋廣龍

（廣西水文水資源局南寧水文站，南寧 530007）

南寧站定點(diǎn)式ADCP指標(biāo)流速與斷面平均流速關(guān)系模型的建立和應(yīng)用

韋廣龍

（廣西水文水資源局南寧水文站，南寧 530007）

南寧站通過(guò)定點(diǎn)式與走航式聲學(xué)多普勒流速儀（ADCP）比測(cè)，經(jīng)率定后，建立定點(diǎn)式H-ADCP指標(biāo)流速與斷面平均流速關(guān)系模型，以此關(guān)系模型應(yīng)用到流量自動(dòng)監(jiān)測(cè)，實(shí)現(xiàn)南寧水文站定點(diǎn)式H-ADCP流量自動(dòng)監(jiān)測(cè)。

南寧站；定點(diǎn)式ADCP；流速關(guān)系

1 概述

應(yīng)用定點(diǎn)式H—ADCP自動(dòng)監(jiān)測(cè)流量，是解決南寧（三）站受水利樞紐蓄水回水影響所造成水位流量關(guān)系不穩(wěn)定、測(cè)流多、定線難等問(wèn)題的新技術(shù)，在定點(diǎn)式H—ADCP應(yīng)用前，要選好定點(diǎn)式H—AD?CP、優(yōu)選測(cè)點(diǎn)、安裝設(shè)備，通過(guò)定點(diǎn)式H—ADCP與走航式ADCP流量比測(cè)，進(jìn)行指標(biāo)流速Vi～斷面平均流速Vm關(guān)系率定，建立指標(biāo)流速Vi～斷面平均流速Vm關(guān)系模型，達(dá)到定點(diǎn)式H—ADCP正式投入生產(chǎn)應(yīng)用的目的，實(shí)現(xiàn)南寧（三）站流量自動(dòng)監(jiān)測(cè)。

（1）河流特性。南寧（三）站為郁江上游控制站，集水面積72 656 km2，河段順直，水流含沙量小，為懸移質(zhì)泥沙，水位年最大變幅17.5m，流量年最大變幅13 400 m3/s。因受下游西津電站調(diào)蓄水影響，無(wú)穩(wěn)定水位流量關(guān)系。

（2）斷面特征。南寧（三）站斷面為梯形（見(jiàn)圖1），河床穩(wěn)定，河底為沙卵石，無(wú)水草，流向與斷面垂直，無(wú)回流，適合使用ADCP測(cè)流。

（3）水文特征。南寧（三）站多年最低、最高水位60.65 m～77.95 m，多年最小、最大流速0.045 m/s～2.45 m/s（2001年），多年最小、最大流量85.0 m3/s～13 400 m3/s，多年最小、最大含沙量0.0 kg/m3～5.71 kg/m3。

圖1 南寧（三）站流速斷面及水平ADCP位置示意圖

2 聲學(xué)多普勒測(cè)流原理

ADCP一般配有4個(gè)（或3個(gè)）換能器，換能器與ADCP軸線成一夾角，每個(gè)換能器即是發(fā)射器也是接收器。換能器發(fā)射的聲波具有指向性。換能器某一固定頻率的聲波，然后接收被水體中顆粒物散射回來(lái)的聲波。一般認(rèn)定顆粒物的運(yùn)動(dòng)速度與水體流速相同，當(dāng)顆粒物的運(yùn)動(dòng)方向接近換能器時(shí)，換能器聽(tīng)到的回波頻率高；當(dāng)顆粒物的運(yùn)動(dòng)方向是背離換能器時(shí)，換能器聽(tīng)到的回波頻率低。聲學(xué)多普勒頻移，即回波頻率與發(fā)射波頻率之差[1]，由下式確定：

式中：FD——聲學(xué)多普勒頻移；

F——發(fā)射波頻率；

V——顆粒物沿聲速方向的移動(dòng)速度；

C——聲波在水中的傳播速度。

ADCP每一個(gè)聲速軸線即為一個(gè)聲束坐標(biāo)。每個(gè)換能器測(cè)驗(yàn)的流速是水流沿其聲束坐標(biāo)方向的流速分量。南寧（三）站定點(diǎn)式H-ADCP為二個(gè)換能器。

3 選用設(shè)備

《聲學(xué)多普勒流量測(cè)驗(yàn)規(guī)范》（SL337—2006）2.2.1“定點(diǎn)式流量測(cè)驗(yàn)，應(yīng)根據(jù)所測(cè)斷面的河寬、水深、流速和含沙量等情況選用合適的聲學(xué)多普勒流速儀”。

南寧站使用定點(diǎn)式H-ADCP是由TRDI公司生產(chǎn)的固定式多普勒流速剖面儀600 kHz（90 m）型，為ChannelMaster系列H-ADCP。主要技術(shù)參數(shù)如下：工作頻率600 kHz；換能器發(fā)射（聲速）角20°；測(cè)速范圍±0.01～±20.0 m/s（水平方向）；速度分辨率0.1 cm/s；流速精度0.25%，測(cè)量范圍1.0～20 m；盲區(qū)0.5 m；單元數(shù)128；單元尺寸100 cm；聲速1 500.0 m/s；水體溫度1.0℃～38.0℃；鹽度0ppt；內(nèi)存2 M，內(nèi)部聲學(xué)水位計(jì)和內(nèi)部壓力水位計(jì)。符合該站水文特征。

4 設(shè)備安裝

（1）安裝位置。南寧（三）站定點(diǎn)式聲學(xué)多普勒流速剖面儀（H-ADCP）安裝在流速儀測(cè)流斷面上游左岸10.0 m處，安裝高程為59.20 m，起點(diǎn)距60.0 m，距河底1.5 m。固定安裝在特制的鐵架上。

（2）安裝要求。南寧（三）站定點(diǎn)式H-ADCP，根據(jù)《聲學(xué)多普勒流量規(guī)范》第2.2.4、2.2.5、3.2.2、3.2.3的要求及規(guī)定進(jìn)行安裝，電源采用直流電供電，外加太陽(yáng)能充電，確保電壓保持11 V以上。HADCP的橫搖（度）為0.200，縱搖（度）0.000，波特率9 600，內(nèi)驅(qū)版本：28.36，采樣間隔時(shí)間：5 min。

5 比測(cè)及成果分析

南寧（三）站定點(diǎn)式H-ADCP與RS-1000型走航式ADCP進(jìn)行比測(cè)，比測(cè)方法、程序如下：在同一測(cè)流斷面同時(shí)進(jìn)行比測(cè)，走航式ACDP施測(cè)二個(gè)測(cè)回?cái)嗝媪髁浚瑑纱瘟髁空`差應(yīng)小于±5.0%，固定式H-ADCP按單元數(shù)5～65，每5 min采集指標(biāo)流速一次，在比測(cè)時(shí)間內(nèi)的幾次采集指標(biāo)流速的平均值。以走航式ADCP實(shí)測(cè)流量Qm、相應(yīng)水位HC與定點(diǎn)式H-ADCP指標(biāo)流速Vi采用定線通軟件進(jìn)行率定。通過(guò)相關(guān)分析法來(lái)確定斷面平均流速Vm的率定系數(shù)。比測(cè)依據(jù)：根據(jù)《聲學(xué)多普勒流量測(cè)驗(yàn)規(guī)范》定點(diǎn)式流量測(cè)驗(yàn)的要求第2.2.1進(jìn)行比測(cè)。

5.1 比測(cè)方法

（1）比測(cè)前期。比測(cè)前，將走航式ADCP換能器入水深度為0.2 m，入水深度雖較小，但測(cè)流過(guò)程中不發(fā)生換能器露出水面，產(chǎn)生“空蝕”現(xiàn)象，從而保證了測(cè)驗(yàn)精度。走航式ADCP比測(cè)位置在流速儀測(cè)流斷面處，測(cè)船盡可能沿著測(cè)流斷面線勻速橫渡施測(cè)，測(cè)驗(yàn)時(shí)測(cè)船平均速度力求小于或等于水流流速；固定式H-ADCP在比測(cè)前，電腦開(kāi)機(jī)后，運(yùn)行H-ADCP程序，在第一設(shè)置框（見(jiàn)圖2）點(diǎn)擊指定數(shù)據(jù)，點(diǎn)擊OK后，進(jìn)入第二設(shè)置框（見(jiàn)圖3）點(diǎn)擊指定數(shù)據(jù)，導(dǎo)入當(dāng)年大斷面數(shù)據(jù)，再點(diǎn)擊OK，即可進(jìn)行固定式H-ADCP比測(cè)。

圖2 運(yùn)行H—ADCP程序第一設(shè)置框

圖3 運(yùn)行H—ADCP程序第二設(shè)置框

（2）比測(cè)過(guò)程。比測(cè)開(kāi)始：在同一時(shí)間，同一流速儀測(cè)流斷面，進(jìn)行走航式ADCP和定點(diǎn)式H-ADCP比測(cè)。取走航式ADCP施測(cè)開(kāi)始時(shí)間為比測(cè)開(kāi)始時(shí)間，讀取對(duì)應(yīng)施測(cè)開(kāi)始的固態(tài)水位，走航式ADCP采用兩個(gè)測(cè)回，要求每個(gè)測(cè)回流量誤差應(yīng)小于±5.0%，且兩個(gè)測(cè)回流量誤差應(yīng)小于±5.0%，即取兩個(gè)測(cè)回流量的平均值為比測(cè)流量Qm；比測(cè)結(jié)束：取走航式ADCP施測(cè)結(jié)束時(shí)間為比測(cè)結(jié)束時(shí)間，并讀取對(duì)應(yīng)施測(cè)開(kāi)始的固態(tài)水位，取走航式ADCP施測(cè)開(kāi)始、結(jié)束時(shí)間的平均值為比測(cè)時(shí)間Tm，所對(duì)應(yīng)施測(cè)開(kāi)始、結(jié)束的固態(tài)水位平均值為比測(cè)水位Hm。負(fù)責(zé)定點(diǎn)式H-ADCP操作，按走航式ADCP施測(cè)開(kāi)始、結(jié)束時(shí)間內(nèi)，摘取每5 min采集的各個(gè)指標(biāo)流速值，取其平均值為比測(cè)指標(biāo)流速Vi，比測(cè)過(guò)程結(jié)束。

（3）比測(cè)要求。按照《聲學(xué)多普勒流量測(cè)驗(yàn)規(guī)范》第5章第5.2.2的規(guī)定：①相關(guān)關(guān)系的建立應(yīng)以儀器實(shí)測(cè)流速和實(shí)測(cè)斷面平均流速為依據(jù)；②斷面平均流速可用流速儀精測(cè)法或走航式多普勒流速儀測(cè)量；③應(yīng)收集不同水位級(jí)的關(guān)系率定資料，分析各種情況下相關(guān)關(guān)系；④各水流條件關(guān)系率定的樣本數(shù)應(yīng)大于30；⑤定線精度指標(biāo)應(yīng)符合SL247-1999要求；⑥應(yīng)選取相關(guān)關(guān)系最優(yōu)的關(guān)系線推算流量；⑦不同水情可采用不同的相關(guān)關(guān)系。因此，南寧（三）站定點(diǎn)式流量測(cè)驗(yàn)比測(cè)應(yīng)大于30測(cè)次，按不同水位級(jí)、不同的流速進(jìn)行比測(cè)，與走航式ADCP比測(cè)流量誤差應(yīng)小于±5.0%。

5.2 比測(cè)資料

南寧（三）站600 kHz（90 m）型定點(diǎn)式H-ADCP與RS-1000型走航式ADCP于2010年5月19日～9月17日進(jìn)行了39次比測(cè)，比測(cè)水位范圍61.08～69.58 m，流量變幅308～5 490 m3/s;比測(cè)流速范圍0.04～1.15，與走航式ADCP比測(cè)流量誤差應(yīng)小于±5.0%。

5.3 比測(cè)分析

將南寧（三）站600 KHz（90 m）型定點(diǎn)式H-ADCP比測(cè)39次的資料在H-ADCP運(yùn)行程序進(jìn)行分析，單元數(shù)級(jí)別分成：5～50；5～55；5～60；5～65；5～70；5～75等5個(gè)進(jìn)行重新計(jì)算單元數(shù)級(jí)別的流量誤差。計(jì)算方法：將走航式ADCP比測(cè)39次的流量、水位和定點(diǎn)式H-ADCP的指標(biāo)流速Vi，錄入H-AD?CP定線通程序中進(jìn)行計(jì)算。經(jīng)計(jì)算，優(yōu)選單元數(shù)：5～65，所得南寧（三）站2010年斷面平均流速Vm與指標(biāo)流速Vi相關(guān)式及率定指標(biāo)參數(shù)，分別為：

一元線性：

（2）（3）（4）式中：

Vm——斷面平均流速，m/s；

Vi——指標(biāo)流速，m/s；

Hi——相應(yīng)水位，m。

通過(guò)對(duì)39份比測(cè)資料進(jìn)行指標(biāo)流速法率定，采用一元線性、二元線性、一元二次多項(xiàng)式計(jì)算，合格率分別為89.74%、100.0%、97.44%（見(jiàn)表1），各項(xiàng)技術(shù)指標(biāo)符合定點(diǎn)式H-ADCP流量測(cè)驗(yàn)的技術(shù)要求。

經(jīng)分析，采用二元線性法計(jì)算得各種技術(shù)指標(biāo)優(yōu)于前后兩各分析法。因此，優(yōu)選二元線性法的關(guān)系式：Vm=0.1 494+（0.1 979+0.0 108 Hi）Vi為南寧（三）站定點(diǎn)式H-ADCP率定公式（見(jiàn)圖4），用于推算流量。

表1 南寧（三）站2010年定點(diǎn)式H-ADCP指標(biāo)流速率定參數(shù)計(jì)算表

圖4 南寧（三）站2010年斷面平均流速Vm與指標(biāo)流速Vi相關(guān)線

6 應(yīng)用

6.1 應(yīng)用模型

南寧（三）站2010年定點(diǎn)式H-ADCP通過(guò)比測(cè)分析，得到南寧（三）站指標(biāo)流速Vi～斷面平均流速Vm關(guān)系模型，即率定公式：Vm=0.1 494+（0.1 979+0.0 108 Hi）Vi。其中，率定參數(shù)b1=0.1 494，b2=0.1 979，b3=0.0108。對(duì)率定公式：Vm=0.1494+（0.1979+0.0108 Hi）Vi進(jìn)行檢驗(yàn)，相關(guān)線符號(hào)、適線、偏離檢驗(yàn)均合理，曲線標(biāo)準(zhǔn)差2.37%，系統(tǒng)誤差-0.01%，符合《水文資料整編規(guī)范》流量檢驗(yàn)要求。

6.2 數(shù)據(jù)采集

南寧（三）站定點(diǎn)式H-ADCP按照比測(cè)率定得指標(biāo)流速Vi～斷面平均流速Vm關(guān)系模型參數(shù)輸入運(yùn)行，進(jìn)行流量自動(dòng)監(jiān)測(cè)，按每5 min1個(gè)時(shí)段監(jiān)測(cè)1次，每1次指標(biāo)流速監(jiān)測(cè)歷時(shí)為100 s，每2天采集定點(diǎn)式H-ADCP流量數(shù)據(jù)一次，采集的指標(biāo)流速原始數(shù)據(jù)通過(guò)平滑處理后，根據(jù)率定公式：Vm=0.1 494+（0.1 979+0.0 108 Hi）Vi，計(jì)算時(shí)段斷面平均流速Vm和流量Qm，以每5 min讀取1個(gè)流量數(shù)據(jù)，采用實(shí)測(cè)流量過(guò)程線進(jìn)行流量整理，并找出月最大、月最小流量。

采用率定公式：Vm=0.1 494+（0.1 979+0.0 108 Hi）Vi進(jìn)行南寧（三）站2010年7月份流量推算。經(jīng)計(jì)算得，逐日平均流量的平均誤差為0.57%，最大誤差為-2.65%，月平均流量誤差為-0.53%，月最大流量誤差為0，月最小流量誤差為1.4%，符合《水文資料整編規(guī)范》流量整編要求。

6.3 應(yīng)用范圍

南寧（三）站2010年定點(diǎn)式H-ADCP水位范圍：61.08～69.58 m，流量：308～5 490 m3/s;流速范圍：0.04～1.15 m/s。

根據(jù)水利部《聲學(xué)多普勒流量測(cè)驗(yàn)規(guī)范》第2.2.6條:根據(jù)試驗(yàn)資料確定的相關(guān)關(guān)系或模型應(yīng)進(jìn)行校測(cè)，校測(cè)的樣本數(shù)應(yīng)不少于10次。當(dāng)校測(cè)點(diǎn)明顯偏離原率定相關(guān)關(guān)系時(shí)，應(yīng)及時(shí)加密校測(cè)測(cè)次，以滿足重新率定相關(guān)關(guān)系。校測(cè)結(jié)果宜用統(tǒng)計(jì)檢驗(yàn)法進(jìn)行檢驗(yàn)。當(dāng)原相關(guān)關(guān)系與校測(cè)樣本有顯著性差異時(shí)，應(yīng)重新進(jìn)行試驗(yàn)并率定相關(guān)關(guān)系。

7 重新率定相關(guān)線

7.1 校測(cè)檢驗(yàn)

南寧（三）站2010年定點(diǎn)式H-ADCP經(jīng)比測(cè)分析，確定率定公式：Vm=0.1 494+（0.1 979+0.0 108Hi）Vi，上報(bào)主管機(jī)關(guān)審批后，同意從2011年1月1日起，正式應(yīng)用定點(diǎn)式H-ADCP作為流量測(cè)驗(yàn)的常規(guī)方法，其成果用于資料整編。應(yīng)用范圍：水位范圍：61.00～69.50 m；流速范圍：0.05～1.15 m/s。

在采用定點(diǎn)式H-ADCP測(cè)流期間，應(yīng)嚴(yán)格按《聲學(xué)多普勒流量測(cè)驗(yàn)規(guī)范》（SL337-2006）的有關(guān)規(guī)定做好流速關(guān)系校測(cè)工作，每年校測(cè)次數(shù)不少于10次，要求校測(cè)點(diǎn)據(jù)基本均勻分布于各流速級(jí)。并繼續(xù)做好水位69.50 m以上及指標(biāo)流速大于1.15 m/s的比測(cè)率定工作。

南寧（三）站2011年定點(diǎn)式H-ADCP在校測(cè)中，出現(xiàn)了水位71.65 m大于69.50 m，指標(biāo)流速1.22 m/s大于 1.15 m/s，指標(biāo)流速 0.03 m/s小于 0.05 m/s；水位、指標(biāo)流速均超出應(yīng)用范圍，校測(cè)點(diǎn)明顯偏離原率定相關(guān)線。因此，加密校測(cè)測(cè)次，共校測(cè)39次。對(duì)校測(cè)39次指標(biāo)流速按原定率定公式：Vm=0.1 494+（0.1 979+0.0 108Hi）Vi進(jìn)行流量檢驗(yàn)。檢驗(yàn)結(jié)果是：相關(guān)線符號(hào)、適線、偏離檢驗(yàn)均不合理，曲線標(biāo)準(zhǔn)差16.51%大于4.0%，不合理，系統(tǒng)誤差-10.1%大于-2.0%，不符合《水文資料整編規(guī)范》流量檢驗(yàn)要求。

7.2 重新率定

將南寧（三）站2011年定點(diǎn)式H-ADCP校測(cè)39次重新率定。按6.2比測(cè)分析的方法進(jìn)行率定，以率定分析，確定重新率定相關(guān)式：Vm=0.0 473+（1.3 821-0.0 051Hi）Vi（見(jiàn)圖 4），其中，率定參數(shù) b1=0.0 473，b2=1.3 821，b3=-0.0 051。水位范圍：61.00～71.65 m；流速范圍：0.03～1.22 m/s。

圖4 南寧（三）站2011年斷面平均流速Vm與指標(biāo)流速Vi相關(guān)線

7.3 相關(guān)分析

南寧（三）站2011年定點(diǎn)式H-ADCP重新率定相關(guān)線，經(jīng)相關(guān)法檢驗(yàn)得知，相關(guān)線符號(hào)、適線、偏離檢驗(yàn)均合理，曲線標(biāo)準(zhǔn)差1.90%，合理，系統(tǒng)誤差0.001%，符合《水文資料整編規(guī)范》流量檢驗(yàn)要求。

南寧（三）站2011年定點(diǎn)式H-ADCP重新率定相關(guān)線，相關(guān)系數(shù)：0.99 965，合格率100%。

9 結(jié)語(yǔ)

南寧（三）站定點(diǎn)式H-ADCP按照《聲學(xué)多普勒流量測(cè)驗(yàn)規(guī)范》（SL337-2006）的有關(guān)規(guī)定，于2010年與走航式ADCP進(jìn)行比測(cè)，2011年校測(cè)，與原相關(guān)關(guān)系與校測(cè)樣本有顯著性差異，重新率定相關(guān)線得：Vm=0.0 473+（1.3 821-0.0 051Hi）Vi，應(yīng)用范圍，水位范圍：61.00～71.65 m；流速范圍：0.03～1.22 m/s。

南寧（三）站定點(diǎn)式多普勒流速剖面儀600 kHz（90 m）型（定點(diǎn)式H-ADCP）正式應(yīng)用，實(shí)現(xiàn)了流量自動(dòng)化監(jiān)測(cè)，解決了因西津電站調(diào)蓄水影響，造成測(cè)流多，水位流量關(guān)系線多又亂，過(guò)渡線多，定線難諸多問(wèn)題。定點(diǎn)式ADCP應(yīng)用的優(yōu)點(diǎn)多，少測(cè)流，不定線；減輕工作強(qiáng)度，減少測(cè)驗(yàn)成本；提高工作效率，提高監(jiān)測(cè)水平；安全快捷，準(zhǔn)確可靠。

[1] 張?jiān)：＃晫W(xué)多普勒流速儀在浙江水文中的應(yīng)用[J]，水文，2010（1）：51.

[2]SL337-2006，聲學(xué)多普勒流量測(cè)驗(yàn)規(guī)范[S].

（責(zé)任編輯：劉征湛）

Development and application of a model describing the relationship between index velocity from H-ADCP and mean cross-section velocity at Nanning Hydrological Station

WEI Guang-long
（Nanning Hydrological Station of Guangxi Hydrology and Water Resources Bureau，Nanning 530007,China）

At Nanning Hydrologicd Station,comparative measurements between H-ADCP and vessel-mounted AD?CP were conducted.Based on the calibration results of comparative measurement data,a model describing the rela?tionship between index velocity from H-ADCP and mean cross-section velocity was developed.Application of this model permits automatic flow measurement by H-ADCP at Nanning Hydrological Station.

Nanning Hydrological Station；horizontal acoustic Doppler current profiler（H-ADCP）;velocity rela?tionship;

P332

B

1003-1510（2016）02-0017-05

聲 明

《廣西水利水電》編輯部

2016-02-23

韋廣龍（1958-）男（壯族）廣西上思人，廣西水文水資源局南寧水文站工程師，從事陸地水文監(jiān)測(cè)和計(jì)算工作。

本刊已許可中國(guó)學(xué)術(shù)期刊（光盤(pán)版）電子雜志社在中國(guó)知網(wǎng)及其系列數(shù)據(jù)庫(kù)產(chǎn)品中以數(shù)字化方式復(fù)制、匯編、發(fā)行、信息網(wǎng)絡(luò)傳播本刊全文。該社著作權(quán)使用費(fèi)與本刊稿酬一并支付。作者向本刊提交文章發(fā)表的行為即視為同意上述聲明。