黃 煒，唐運憶，趙德友，丁昌言，李國棟

（1. 江蘇省水文水資源勘測局，江蘇 南京 210029；2. 江蘇省水文水資源勘測局揚州分局，江蘇 揚州 225000；3. 江蘇省水文水資源勘測局鹽城分局，江蘇 鹽城 224051）

基于聲學(xué)多普勒流速儀測流的斷面流速分布研究

黃 煒1，唐運憶1，趙德友1，丁昌言2，李國棟3

（1. 江蘇省水文水資源勘測局，江蘇 南京 210029；2. 江蘇省水文水資源勘測局揚州分局，江蘇 揚州 225000；3. 江蘇省水文水資源勘測局鹽城分局，江蘇 鹽城 224051）

選擇涇河、射陽河閘 2 個水文站的走航式聲學(xué)多普勒流速儀測流數(shù)據(jù)，進行測流斷面上的垂向、水平層流速分布及垂向平均流速的橫向分布研究。研究結(jié)果表明：對于順直河段，垂線點流速分布基本符合理論流速分布；對于受感潮影響的測流斷面，相對流速脈動和均化影響較大，難以找到流速分布規(guī)律。通過基于斷面流速分布特點的斷面有效測速單元確定，可以為測速垂線和垂線上測速點的布設(shè)、固定式聲學(xué)多普勒流速儀安裝位置選擇等提供技術(shù)支撐。

聲學(xué)多普勒流速儀；測流斷面；垂向流速分布；水平層流速分布；垂向平均流速橫向分布；有效測速單元

0 引言

斷面上的流速分布是指流速在斷面上沿水深和橫斷面方向的變化規(guī)律[1]，研究測流斷面上的流速分布對分析河床演變和泥沙運動規(guī)律，以及開展水文測驗等具有重要意義。傳統(tǒng)的流速垂線分布測量方法采用單點式流速儀逐一對各條流速垂線上不同深度測點進行點測量。但這種垂線逐點測量方法費時、費力，已不能滿足水文測驗現(xiàn)代化的發(fā)展要求。近年來，基于聲學(xué)多普勒原理的多普勒剖面流速測試技術(shù)得到了快速發(fā)展，并已普遍應(yīng)用于流速、流量測驗[2-5]。相比于傳統(tǒng)的固定式垂線測量，具有速度快，精度高，不擾動流場等優(yōu)點，且聲學(xué)多普勒流速儀（ADCP）空間采樣密度高，測速單元較多，測流時可以覆蓋全部過水?dāng)嗝?，能夠提供更為精?xì)的全斷面二維流速數(shù)據(jù)，為斷面流速分布模型的研究提供數(shù)據(jù)支撐。

選擇涇河、射陽河閘 2 個水文站的走航式 ADCP測流數(shù)據(jù)，進行測流斷面上的流速分布研究，旨在通過摸清流速分布進而指導(dǎo)測速垂線和垂線上測速點的布設(shè)、固定式 ADCP 安裝位置的選擇，以及為流速、流量測驗中的異常數(shù)據(jù)的剔除，斷面特性對流速的影響分析等提供依據(jù)。

1 ADCP 測流原理介紹

ADCP 是利用聲學(xué)多普勒效應(yīng)進行測流的[6-7]。換能器以固定的頻率向水中發(fā)射聲波，當(dāng)聲波碰到水中的聲散射體（如浮游生物、懸浮粒子等）后，少量脈沖返回并被 ADCP 接收，反射回來的小數(shù)量的聲脈沖頻率發(fā)生改變，與發(fā)射聲波的頻率之間會產(chǎn)生一個頻差，即多普勒頻移，通過測量得到多普勒頻移即可得到散射粒子的速度。由于水體中散射粒子懸浮于水中，它們的速度通常與水流速度保持一致，故顆粒物的速度即為水流的流速。ADCP 可以接受經(jīng)懸浮顆粒反射后多個時段的回波信號，從而能夠?qū)Σ煌钐幍乃魉俣冗M行測量，可同時獲得垂向不同水深處的流速值[8]。ADCP 在移動過程中所測的垂線可以是非常多的條數(shù)，在每條“垂線”上設(shè)置不同厚度的深度單元，相應(yīng)測到幾十甚至上百個不同深度點的流速值。通過軟件提取各個深度單元流速，能夠提供更為精細(xì)的全斷面二維流速數(shù)據(jù)，將全斷面二維流速數(shù)據(jù)組合成流速剖面，進而得到測驗斷面的流速分布。通過對流速分布情況的分析，可為斷面流速分布模型的研究提供數(shù)據(jù)支撐，并指導(dǎo)固定式 ADCP 安裝和垂線等位置的選擇。

2 ADCP 數(shù)據(jù)采集與處理

2.1 測站選擇

江蘇省內(nèi)水系發(fā)達，航運頻繁，且水利工程眾多，大多測驗河段受航運和水利工程調(diào)度影響較大。

大運河涇河站為南水北調(diào)市際控制站，測流斷面河寬 180 m，常水位時水面寬約 140 m。測流斷面的上、下游數(shù)千米范圍內(nèi)，河道基本順直。測流斷面處的河道，為人工開挖的規(guī)則梯形斷面，河床土質(zhì)以細(xì)顆粒粘土為主，兩岸為塊石護坡，斷面沖淤年度變化率不明顯，河勢比較穩(wěn)定。該河道為北煤南運的重要通道，航道等級為二級，年通航能力約5 億 t，給該站水文測驗工作帶來了一定影響。射陽河閘為江蘇省沿海最大的擋潮閘壩，該閘有 35 孔，每孔凈寬 10 m，閘底板高程 -3.50 m，閘身全長410.1 m，是集擋潮、御鹵、蓄淡、抗旱、排洪于一體的綜合型閘壩。測站測驗河段順直，水面寬約270 m，高水時無岔流、串溝，河床及河岸均為砂土，閘上游河段沖淤變化甚微，主流一般穩(wěn)定，但受閘門開啟位置不同而略有擺動。

為此，選擇涇河、射陽河閘 2 個水文站的走航式 ADCP 測流數(shù)據(jù)，做為江蘇省內(nèi)典型測驗斷面進行測流斷面上的流速分布研究。

根據(jù)不同測驗環(huán)境共選擇 4 份測驗樣本，各樣本基本情況如表 1 所示。

表 1 測驗數(shù)據(jù)情況一覽表

2.2 數(shù)據(jù)處理

在 ADCP 測流中，每完成 1 次采樣，即 ADCP每完成 1 次信號發(fā)射、接受和處理，則視為 1 個呯，多個呯采樣的平均值即為呯集合。根據(jù)各測次呯集合對應(yīng)的起點距數(shù)據(jù)，涇河站各測次呯集合均化為 1 00 根垂線，射陽河閘站各測次呯集合均化為1 000 根垂線，便于進行樣本內(nèi)各測次數(shù)據(jù)的疊加分析。應(yīng)用 WinriverⅡ 軟件回放每半測回測量數(shù)據(jù)，輸出 ASCⅡ 數(shù)據(jù)文件，提取測流斷面上的二維流速、流向數(shù)據(jù)，以及每個呯集合對應(yīng)的起點距數(shù)據(jù)。

走航式 ADCP 單次獲取的垂線流速為瞬時的紊動流速，為了求取平均流速分布，需要將走航式ADCP 各次重復(fù)測量獲得的多組瞬時流速作為平均流速，然而，垂線組數(shù)的選取對平均流速的值會產(chǎn)生影響。以涇河為例，以某一測點、某一層數(shù)據(jù)提取垂線流速，對采用不同垂線組數(shù)建立平均流速與垂線組數(shù) N 的關(guān)系曲線，發(fā)現(xiàn)用于計算平均流速的垂線組數(shù) ＞ 70 組（約 70%）時，平均后得到的流速值已趨于穩(wěn)定。因此，在對各測點的流速數(shù)據(jù)進行分析時，選取的瞬時垂線組數(shù)均不小于全部砰集合的 70%。

由于 ADCP 可同時獲取不同水深處的瞬時流速值 vi,j，因此，時均流速垂向分布可通過對各水深處瞬時流速的時間序列進行算術(shù)平均得到。不同水深處的時均流速值與標(biāo)準(zhǔn)差 σi可由下式得到：

式中：vi,j為第 i 層水流瞬時流速，i 表示第 i 層，j 表示第 j 根垂線；n 表示層流速單元個數(shù)；m 表示測點分層數(shù)量。

通過整理各樣本的流速測量成果，發(fā)現(xiàn)各垂線上測點流速與其位置均有密切、一致的對應(yīng)關(guān)系[9]。為考察各測點平均流速測量結(jié)果的可靠性及合理性，利用對測點不同水深處流速標(biāo)準(zhǔn)差 σi進行無量綱化，將測點所有水深處的相對標(biāo)準(zhǔn)差進行平均得到平均相對標(biāo)準(zhǔn)差，計算公式如下：

3 斷面流速分布研究

3.1 垂向流速分布

在 ADCP 實測數(shù)據(jù)樣本中，選擇受外界影響較小的某個測次（即樣本 1 中第 8 測次，樣本 2 中第 48 測次，樣本 3 中第 16 測次，樣本 4 中第 29 測次），任意提取河道中泓 3 根垂線的測速單元流速數(shù)據(jù)，給出中泓垂線點流速分布圖，如圖 1 所示。

圖 1 單次測驗數(shù)據(jù)中泓垂線點流速分布圖

理論上垂線上的流速分布可以擬合為拋物線、對數(shù)、橢圓、指數(shù)型等，由圖 1 可以看出，單個測次中中泓垂線點流速分布基本符合理論流速分布，但是分布較為紊亂，流速脈動較大。

因此，為了避免單次測驗成果中任意垂線的數(shù)據(jù)偶然性，將各樣本中多次測驗數(shù)據(jù)進行綜合，取不同測次同一位置點流速的平均值，得到綜合分析結(jié)果，如圖 2 所示。

圖 2 全體樣本測驗數(shù)據(jù)中泓垂線點流速分布圖

通過圖 1 和 2 可知：

1）單次測驗成果同一垂線上不同水深的單元流速脈動較大，多次平均后可減弱脈動影響；

2）全體樣本測驗垂線上的流速分布，樣本 1，2，3 總體上符合拋物線形，流速隨水深增大而減?。?/p>

3）樣本 1 單元尺寸為 0.5 m，規(guī)律較為明顯，樣本 2 和 3 單元尺寸較?。?.2 m），受流速脈動影響較大，規(guī)律相較于樣本 1 不明顯；

4）各樣本流速分布在接近河底處易產(chǎn)生奇異值，因此在選擇有效測速單元時應(yīng)考慮避開河底部分；

5）樣本 4 測流時流速較小，并受感潮影響，相對流速脈動和均化影響較大，因此難以找到流速分布規(guī)律。

3.2 水平層流速分布

取同一深度的單元流速計算平均值，可得到每一水平層的平均流速，層的深度取決于單元尺寸和號數(shù)。

涇河站測流斷面選取部分垂線（20～50 號垂線）的每層平均流速與全部斷面水平層平均流速進行對比分析，射陽河閘站測流斷面取各水平層的最大流速與全部斷面水平層平均流速作對比分析。水平層和部分垂線的平均流速或水平層最大流速的分布規(guī)律如圖 3 所示。由圖 3 可知：

1）斷面水平層的平均和最大流速的垂向分布與圖 2 中垂線上測點流速的垂向分布具有相似性；

2）樣本 1，2，3 的全斷面水平層平均流速和 20～50 號部分垂線水平層平均流速形狀基本契合，總體分布上符合拋物線形，且流速隨水深增大而減??；

3）樣本 4 測時流速較小，樣本 2 和 3 單元尺寸較小，受流速脈動影響較大，全部垂線水平層平均流速規(guī)律不明顯；

圖 3 斷面水平層流速垂向分布圖

4）樣本 4 各水平層最大流速與平均流速差異較大，說明測流時受外界影響較大，異常數(shù)據(jù)較多。

3.3 垂線平均流速的橫向分布

根據(jù)每條垂線上的單元流速數(shù)據(jù)計算垂線平均和最大流速并進行分析，得到垂線平均和最大流速的橫向分布規(guī)律，如圖 4 所示。

由圖 4 可知：

1）橫斷面上的流速分布與斷面形狀有關(guān)，樣本數(shù)據(jù)基本符合理論的流速分布規(guī)律，近兩岸邊的流速小于中泓流速；

2）樣本 1，2，3 為同一測流斷面，但不同水位和流速大小造成橫斷面流速分布不同，但仍能從圖中看出，起點距 20～80 號之間垂線平均流速較為穩(wěn)定；

圖 4 垂線流速橫向分布圖

3）樣本 2 和 3 各垂線最大流速與平均流速差異較大，說明測流時受外界影響較大，異常數(shù)據(jù)較多，與水平層平均流速表現(xiàn)一致。

3.4 有效測速單元選擇

固定式 ADCP 比測率定精度與安裝位置有很大關(guān)系。將斷面流速分布中最能代表斷面平均流速或與斷面平均流速標(biāo)準(zhǔn)差最小的測速單元定義為有效測速單元。選擇在有效測速單元附近安裝設(shè)備（橫向式 ADCP）和設(shè)置垂線（垂向式 ADCP），對更有效地率定關(guān)系有很大幫助。由于射陽河閘站受感潮影響，相對流速脈動和均化影響較大，水平層最大流速與平均流速差異較大，流速分布規(guī)律不明顯，有效測速單元難以確定，為此取涇河站多個走航式ADCP 測流數(shù)據(jù)，做均化處理后，分析二維斷面流速分布圖，確定有效測速單元。

由分析可知，涇河站垂線點流速分布基本符合理論流速分布，水平層平均流速和部分垂線的水平層平均流速基本形狀契合，總體分布上符合拋物線形，但在接近河底處流速易產(chǎn)生奇異值，因此在選擇有效測速單元時需避開河底部分。因此，可以認(rèn)為涇河站測流斷面中垂向上 0.20～0.60 相對水深處，橫向上 20～40 號垂線處為測流的有效測速單元。

采用有效測速單元流速均值作代表流速，并對測得的比測率定測驗成果進行分析，得到?jīng)芎诱局笜?biāo)流速 vi和斷面平均流速 vp關(guān)系，即 vp= 0.916 vi，相關(guān)系數(shù) R2= 0.998 9。對測點率定關(guān)系曲線進行三項檢驗，檢驗結(jié)果均小于檢驗臨界值，符合資料整編規(guī)范要求。

4 結(jié)語

選擇涇河站、射陽河閘 2 個水文站的走航式ADCP 測流數(shù)據(jù)，建立了江蘇省內(nèi)典型測流斷面上的垂線、水平層流速分布及垂線流速的橫向分布模型，分析了測流斷面上的有效測速單元分布情況。研究結(jié)果表明，當(dāng)河段較為順直穩(wěn)定時，垂線點流速分布基本符合理論流速分布；對于受感潮影響的測流斷面，相對流速脈動和均化影響較大，難以找到流速分布規(guī)律。通過建立有效測速單元，并在儀器設(shè)備安裝時選擇避開流速分布規(guī)律不強或易產(chǎn)生奇異值的區(qū)域，可以更有效地建立斷面平均流速與指標(biāo)流速之間的率定關(guān)系。由于在固定式 ADCP 安裝和運行過程中，受船行或閘門啟閉影響的情況較多，因此在下一步工作中，應(yīng)對受船行波或水利工程控制影響的測流斷面的流速分布情況進行分析，建立不同測驗環(huán)境下的流速分布模型，以期為固定式 ADCP 的安裝位置選擇、成功應(yīng)用及比測率定提供更全面的幫助。

[1] HUANG H N. Index-velocity rating development for rapidly changing flows in an irrigation canal using broadband StreamPro ADCP and ChannelMaster H-ADCP[C] // Proceedings of Rivers’04, First International Conference on Managing Rivers in the 21st Century: Issues and Challenges. Malaysia:Benang, 2004: 146-154.

[2] 王發(fā)君，黃河寧. H-ADCP 流量在線監(jiān)測指標(biāo)流速法定線軟件“定線通”介紹與應(yīng)用[J]. 水文，2007，27(4): 63-65.

[3] 張孝軍，諶紅信，陳冬，等. 水平式 ADCP 安裝方式探討[J]. 水資源研究，2008 (2): 30-31.

[4] 杜耀東. 現(xiàn)代測流測沙技術(shù)研究與應(yīng)用[D]. 武漢：武漢大學(xué)，2012: 99-100.

[5] 鄭慶濤，唐健奇，張志敏. H-ADCP 在流量自動監(jiān)測系統(tǒng)中的應(yīng)用[J]. 水利科技與經(jīng)濟，2007，13 (3): 173-174.

[6] 李志敏. H-ADCP 在線測流系統(tǒng)在思賢滘的應(yīng)用研究[J].廣東水利水電，2006 (5): 10-11.

[7] 詹智慧，沈鴻金. 哨兵型 H-ADCP 在潮流量實時監(jiān)測中的應(yīng)用研究[J]. 人民珠江，2009，30 (2): 49-51.

[8] 胡雪嬌，李健. H-ADCP 在流域引江濟太工作中的實踐[J]. 水利信息化，2010 (6): 28-31.

[9] 蔣建平. 走航式聲學(xué)多普勒流速儀投產(chǎn)方法探討[J]. 水利信息化，2013 (2): 40-43.

Research on velocity distribution of fl ow measurement based on acoustics doppler current pro fi le measurement method

HUANG Wei1, TANG Yunyi1, ZHAO Deyou1, DING Changyan2, LI Guodong3
(1. Jiangsu Province Hydrology and Water Resources Investigation Bureau, Nanjing 210029, China; 2. Jiangsu Province Hydrology and Water Resources Investigation Bureau Yangzhou Branch, Yangzhou 225000, China; 3. Jiangsu Province Hydrology and Water Resources Investigation Bureau Yangzhou Branch, Yancheng 224051, China)

The research of vertical velocity distribution, horizontal velocity distribution, and transverse distribution of vertical mean velocity, is based on the shipboard-type ADCP discharge measurement data which are measured on Jinghe and Shenyang two hydrological stations. The results show that the vertical velocity distribution satisfies the theoretical velocity distribution on the straight section. It is dif fi cult to fi nd the velocity distribution law for the in fl uence of tidal fl ow on the cross section, which has a great in fl uence on relative velocity fl uctuation and homogenization. The effective velocity measurement unit of the cross section is determined by the velocity distribution characteristics of the cross section, which can provide technical support for the layout of velocity measurement points on the vertical and the installation location of fi xed acoustic Doppler velocity meter.

acoustics doppler current pro fi le; measuring cross-section; vertical velocity distribution; horizontal velocity distribution; transverse distribution of vertical mean velocity; effective velocity measuring unit

P335

A

1674-9405(2017)04-0061-07

10.19364/j.1674-9405.2017.04.012

2017-04-10

黃 煒（1981-），男，湖南永州人，工程師，主要從事水文測驗、水信息技術(shù)及水環(huán)境等研究。