潘凌 武治國(guó) 張春萍 沈海超 陳銀 張家銓 毛紅梅

摘要：選擇實(shí)用的測(cè)量方法對(duì)減少城市河渠剖面流速的測(cè)量誤差十分關(guān)鍵。多普勒流量計(jì)采用多點(diǎn)測(cè)量的方法，對(duì)城市河渠的剖面流速進(jìn)行分析，研究城市河渠剖面垂向流速的實(shí)際分布規(guī)律。經(jīng)研究發(fā)現(xiàn)，針對(duì)不同的應(yīng)用場(chǎng)合，垂向流速的分布適合采用對(duì)數(shù)和拋物線兩種曲線進(jìn)行的分段擬合，且分界點(diǎn)相對(duì)水深宜在0.6-0.85范圍內(nèi)，此時(shí)相關(guān)系數(shù)可達(dá)到0.95。采用兩點(diǎn)法與預(yù)測(cè)法相結(jié)合的方式獲取垂向平均流速，兩種方法都具有較好的精度，可獲得較為準(zhǔn)確的明渠剖面流速，對(duì)工程實(shí)際應(yīng)用具有一定的指導(dǎo)和參考價(jià)值。

關(guān)鍵詞：流速測(cè)量;城市河渠;剖面流速;多普勒流量計(jì)

中圖法分類(lèi)號(hào)：TP391.44

文獻(xiàn)標(biāo)志碼：A

DOI：10.15974/j .cnki.slsdkb.2019.08.005

1 研究背景

目前，針對(duì)河渠剖面流速分布規(guī)律的理論研究較多，如常用的剖面流量測(cè)量方法有堰槽法、水位流量關(guān)系法、流速儀法、電磁感應(yīng)法等[1]。部分學(xué)者依據(jù)動(dòng)能方程推導(dǎo)出流量計(jì)算公式[2]，但主要是針對(duì)自然河道，并且參數(shù)繁多，不易應(yīng)用。對(duì)于剖面流速分布規(guī)律與相關(guān)經(jīng)驗(yàn)公式及計(jì)算方法，韓金旭等[3-4]提出了“水面流速法”，可根據(jù)水面單點(diǎn)的流速測(cè)量推算出剖面流速分布規(guī)律。此方法雖然簡(jiǎn)便，但平均流速的誤差對(duì)測(cè)量點(diǎn)的測(cè)量誤差有較大依賴。王二平等[5-7]對(duì)剖面流速的分布規(guī)律提出了垂線拋物線、橫向乘冪函數(shù)算法公式，對(duì)流速剖面分析有一定參考價(jià)值，但所得出的剖面中心區(qū)域的垂向流速誤差較大。由此可見(jiàn)，上述方法在應(yīng)用到河渠的流速測(cè)量時(shí)均存在不同程度的問(wèn)題。隨著城市水環(huán)境問(wèn)題的凸顯，亟需解決城市河渠的剖面流速測(cè)量問(wèn)題，以實(shí)現(xiàn)水資源的合理化管理，因此，探索便于實(shí)際應(yīng)用且較為精確的剖面流量測(cè)量技術(shù)迫在眉睫。

本文以云南省玉溪市內(nèi)河道為研究對(duì)象，通過(guò)現(xiàn)場(chǎng)實(shí)測(cè)與數(shù)據(jù)擬合等方式，分析出剖面流速經(jīng)驗(yàn)擬合公式，結(jié)合在線多普勒流量檢測(cè)設(shè)備，建立城市規(guī)則河道剖面流速分析方法。將理論與實(shí)際相結(jié)合，總結(jié)出一種可應(yīng)用于城市河渠流速、流量在線測(cè)量的方法。

2 城市河渠剖面流速分布規(guī)律

2.1 河渠剖面流速測(cè)點(diǎn)的選擇及布置

為盡可能得到準(zhǔn)確的剖面平均流速，在玉溪市內(nèi)各種河渠流域中，分別挑選不同流域的7個(gè)監(jiān)測(cè)斷面進(jìn)行人工剖面流速測(cè)量，斷面編號(hào)分別為H1、H3、H5、H10、H11、H13、H15（見(jiàn)圖1）。在監(jiān)測(cè)過(guò)程中，H5（玉溪大河流出示范區(qū)）斷面水位較高，其他監(jiān)測(cè)斷面水位較低。玉溪市河道流域分布見(jiàn)圖1。

挑選的河渠寬深比均在15以上，選擇區(qū)域剖面規(guī)則、水流相對(duì)平穩(wěn)。每個(gè)測(cè)量剖面布置9條垂線，每條垂線設(shè)置7個(gè)測(cè)點(diǎn)（部分水位較淺的河道只設(shè)置了3-5個(gè)測(cè)點(diǎn)）。9條垂線的相對(duì)起點(diǎn)間距分別為0.1，0.2，0.3，0.4，0.5，0.6，0.7，0.8，0.9（個(gè)別剖面根據(jù)實(shí)際情況有所調(diào)整），7個(gè)測(cè)點(diǎn)的相對(duì)水深分別為0.95，0.75，0.50，0.375. 0.25，0.125，0.05（某些垂線根據(jù)實(shí)際情況有所調(diào)整），每個(gè)點(diǎn)位進(jìn)行分時(shí)測(cè)量。

現(xiàn)場(chǎng)測(cè)量使用的流速檢測(cè)設(shè)備為便攜式流速儀，型號(hào)為NF-ZBX-1，測(cè)量范圍不大于5 m/s，測(cè)量精度1%。

2.2 河渠剖面流速整體分布規(guī)律

一般河渠的截面形狀會(huì)對(duì)水流的狀態(tài)產(chǎn)生較大影響，進(jìn)而影響斷面流速分布。對(duì)于寬淺的河渠，流速在垂線方向上的分布梯度較大，垂線平均流速沿河寬方向（橫向）的分布較均勻;對(duì)于較窄深的河渠，側(cè)壁對(duì)流速分布的影響較大，流速沿垂線分布較均勻，但流速在河寬方向的分布梯度較大[8]。因此，對(duì)于窄深的河道必須考慮側(cè)壁對(duì)流速分布的影響，而對(duì)于寬淺的河道側(cè)壁影響則較小。在一定的寬深比范圍內(nèi)，河渠的剖面流速分布在垂向或橫向上存在一定規(guī)律。一般當(dāng)剖面規(guī)則且寬深比在15以上時(shí)，可以認(rèn)為是寬淺類(lèi)型的河渠。

以玉溪城市河道現(xiàn)場(chǎng)的兩處勘測(cè)斷面為例，對(duì)兩處斷面的流速數(shù)據(jù)進(jìn)行分析，繪制的仿真圖如圖2～3所示。橫坐標(biāo)軸為測(cè)點(diǎn)距河岸的距離，縱坐標(biāo)軸為測(cè)點(diǎn)距水面的距離，垂向坐標(biāo)軸表示流速。

由圖2-3可以看出，雖然是不同的河道，但剖面流速分布情況近似。最大流速出現(xiàn)在接近水面的位置，由于靠近側(cè)邊的水流受到側(cè)壁的阻力，導(dǎo)致水流速度小于河道中間的水流速度，底部水流速度小于水面的水流速度。表1對(duì)不同河渠剖面的橫向與垂向流速的平均變化率（流速變化率=兩測(cè)點(diǎn)的流速差/兩點(diǎn)間的距離）進(jìn)行了對(duì)比，所有被測(cè)量河道的垂向流速平均變化率遠(yuǎn)遠(yuǎn)大于橫向流速平均變化率，這與盧金友等[8]研究得出規(guī)律相符。

由此可見(jiàn)，對(duì)于規(guī)則的城市河道或明渠，有必要明確垂向的流速分布規(guī)律，其直接影響河渠剖面的流速測(cè)量精度。

2.3 垂向流速分布規(guī)律

2.3.1 現(xiàn)有分布規(guī)律算法討論及存在的問(wèn)題

在城市河道或明渠環(huán)境流速數(shù)據(jù)研究中（寬深比15以上），在河道中某一相對(duì)確定位置，測(cè)線離渠壁距離D在某一范圍內(nèi)時(shí)，發(fā)現(xiàn)各垂線上測(cè)點(diǎn)流速與該點(diǎn)所在位置有著緊密、一致的對(duì)應(yīng)關(guān)系，無(wú)量綱化的相對(duì)流速u(mài)/V（u為點(diǎn)流速，V為垂線平均流速）與y/H（y為測(cè)點(diǎn)離渠底距離，H為水深）的關(guān)系在明渠底部?jī)?nèi)區(qū)符合對(duì)數(shù)分布規(guī)律[9]。

u/V=a+bIn（y/H）

（1）式中，u、y分別為測(cè)線上任意一點(diǎn)的流速與平均流速;a、b為垂向流速分布系數(shù)。玉溪市玉帶河匯人玉溪大河下游相對(duì)起點(diǎn)距為0.5的垂線上的人工測(cè)量數(shù)據(jù)見(jiàn)表2。

用對(duì)數(shù)分布規(guī)律描述無(wú)量綱化的相對(duì)流速u(mài)/V與y/H的關(guān)系，發(fā)現(xiàn)擬合相關(guān)性系數(shù)僅為0.87。為更好地體現(xiàn)垂向流速分布情況，孫東坡等[10]提出實(shí)際河道明渠流速的垂線分布更接近于二次函數(shù)曲線的特征，垂線上相對(duì)流速與相對(duì)水深的無(wú)量綱函數(shù)關(guān)系可表示為

u/V= a（y/H）2+ b（y/H）+c

（2）式中，c為垂向流速分布系數(shù)。對(duì)試驗(yàn)數(shù)據(jù)分析可知，若使用二次函數(shù)來(lái)描述垂向流速的分布，其相關(guān)性有所提高，但仍然不能精確完整描述垂向流速分布情況。因此，在一條垂線上單獨(dú)使用對(duì)數(shù)函數(shù)或二次函數(shù)分析所得出的結(jié)果誤差較大。

2.3.2 針對(duì)現(xiàn)有算法問(wèn)題的解決辦法

為解決現(xiàn)有分析方法在實(shí)際應(yīng)用中存在的問(wèn)題，且要既能保證實(shí)用性，又能有效提高測(cè)量精度，可參照胡云進(jìn)[9]的方法將垂線分為內(nèi)區(qū)與外區(qū)，分別運(yùn)用不同的公式進(jìn)行擬合計(jì)算，但文獻(xiàn)[9]未提出內(nèi)外區(qū)分界點(diǎn)的位置。本文對(duì)玉溪河道實(shí)測(cè)的數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，以H3斷面（玉帶河匯人玉溪大河下游斷面）數(shù)據(jù)為例，在進(jìn)行數(shù)據(jù)擬合的過(guò)程中，通過(guò)不斷調(diào)整兩種公式擬合的分界點(diǎn)位，找出最終可使擬合結(jié)果與實(shí)測(cè)數(shù)據(jù)相關(guān)性最高的點(diǎn)。通過(guò)對(duì)大量數(shù)據(jù)的對(duì)比分析發(fā)現(xiàn)，此分界點(diǎn)一般位于相對(duì)水深y/H值為0.6-0.85的位置。因此，針對(duì)內(nèi)區(qū)與外區(qū)分別用對(duì)數(shù)函數(shù)與二次函數(shù)進(jìn)行擬合，相關(guān)系數(shù)可達(dá)到0.95以上，如圖4所示。

對(duì)于寬淺的過(guò)水剖面，其中心區(qū)為二維流動(dòng)，各垂線流速分布相同，兩側(cè)為受邊界粗糙率影響的邊壁區(qū)，兩個(gè)區(qū)域的分界與水流強(qiáng)度有關(guān)。沿剖面橫向水流運(yùn)動(dòng)可認(rèn)為是對(duì)稱的。參數(shù)a，b，c可借鑒同類(lèi)河道數(shù)據(jù)。

當(dāng)河道中心垂向上可同時(shí)安裝兩個(gè)測(cè)點(diǎn)時(shí)，還可采用兩點(diǎn)法來(lái)計(jì)算垂線的平均流速。河道中心垂向平均流速可用中心任意兩個(gè)流速值和其對(duì)應(yīng)的水位值來(lái)求得[11]：式中，uδ1和uδ2分別為兩個(gè)測(cè)點(diǎn)的流速;δ1和δ2為兩個(gè)測(cè)點(diǎn)的流速到河底的距離。將式（3）中求得的河道中心平均流速代人式（1），即可求出a，b的值，再將式（1）應(yīng)用于其他垂線上的單個(gè)測(cè)量點(diǎn)。這樣只用一個(gè)點(diǎn)，就可近似計(jì)算出該剖面的其他測(cè)量點(diǎn)所在垂線的平均流速。

3 多普勒河渠流速剖面在線監(jiān)測(cè)方法及驗(yàn)證

3.1方法簡(jiǎn)介

在對(duì)試驗(yàn)數(shù)據(jù)的分析與研究的基礎(chǔ)上，結(jié)合在線多普勒超聲波流量計(jì)（型號(hào)為NF-LSX-1，測(cè)量范圍為一5 m/s +5 m/s，精度為±1 cm/s），將上述算法嵌入在線流量測(cè)量設(shè)備中。依據(jù)現(xiàn)場(chǎng)實(shí)際環(huán)境及上述討論的流量計(jì)布點(diǎn)規(guī)律，在測(cè)量區(qū)合理安裝布點(diǎn)。應(yīng)用于城市河道明渠流量監(jiān)測(cè)時(shí)，記錄監(jiān)測(cè)點(diǎn)位河道剖面長(zhǎng)、寬、流量計(jì)安裝位置（流量計(jì)安裝高度、測(cè)點(diǎn)與河邊距離）等信息，并對(duì)在線多普勒超聲波流量計(jì)表頭進(jìn)行配置。在流量計(jì)工況良好情況下，便可實(shí)現(xiàn)河道流速的在線監(jiān)測(cè)。

對(duì)于河渠剖面橫向流速的分布規(guī)律，由于邊壁效應(yīng)，在河道邊壁附近的垂線平均流速會(huì)小于河道中心的垂線平均流速，但在距邊壁一定距離后，垂線平均流速基本保持不變。為此，對(duì)于橫向流速的計(jì)算，不采用已有的橫向流速變化規(guī)律公式，而是用一種更有實(shí)際測(cè)量意義的計(jì)算方法。實(shí)際觀測(cè)發(fā)現(xiàn)，測(cè)試用的室外渠道因受紊流和渠底淤積物的影響橫向流速分布規(guī)律不太明顯，只能看出靠近側(cè)壁的流速較中間偏低。采取多點(diǎn)測(cè)量的方式來(lái)確定橫向流速。測(cè)點(diǎn)位置選在橫向流速變化相對(duì)劇烈的點(diǎn)：渠底規(guī)則、寬度較窄的河道剖面，在橫向上一般布置1-2條垂向測(cè)點(diǎn);渠底規(guī)則、河道有一定寬度、紊流現(xiàn)象不明顯的剖面，一般布置3-4條垂線測(cè)點(diǎn);河面較寬、紊流明顯的截面，布點(diǎn)4-6條垂線測(cè)點(diǎn)，每條垂線測(cè)點(diǎn)上采用1-2個(gè)多普勒流量計(jì)進(jìn)行流速測(cè)量，每個(gè)流量計(jì)的相對(duì)水深分別為0.65-0.80與0.30-0.50（具體位置根據(jù)現(xiàn)場(chǎng)安裝情況進(jìn)行調(diào)整），一般單條垂線的橫向覆蓋區(qū)域范圍是3-6 m。

實(shí)測(cè)發(fā)現(xiàn)，在垂向流速測(cè)量時(shí)，將流量計(jì)置于剖面內(nèi)區(qū)（對(duì)數(shù)區(qū)）可獲得更準(zhǔn)確的測(cè)量數(shù)據(jù)。實(shí)測(cè)時(shí)，將河渠剖面橫向等分為5-10個(gè)區(qū)域（區(qū)域數(shù)量通常大于垂向測(cè)點(diǎn)數(shù)量），將整理后的實(shí)測(cè)數(shù)據(jù)代人式（1）和式（2）可得垂向平均流速。對(duì)于其他未安裝流量計(jì)區(qū)域的垂向平均流速，可通過(guò)對(duì)相鄰區(qū)域垂向平均流速進(jìn)行加權(quán)平均得出。各分區(qū)流速的平均值即為剖面平均流速，再乘以剖面面積即可得到剖面瞬時(shí)流量。

3.2 分析驗(yàn)證

在玉溪市內(nèi)河道中，選取多個(gè)河道斷面規(guī)則度較好的點(diǎn)位，分別用預(yù)測(cè)法和兩點(diǎn)法對(duì)部分河道進(jìn)行剖面流速分析，并將計(jì)算出的平均流速與人工實(shí)測(cè)流速相比較，結(jié)果如表3-4所示。

由表3-4可知，這兩種方法在實(shí)際規(guī)則剖面城市河道流速測(cè)量中的誤差均在5%以內(nèi)，對(duì)于河道流速測(cè)量要求來(lái)說(shuō)精度相對(duì)較高，且操作方法可行，數(shù)據(jù)可實(shí)時(shí)獲取。

4 結(jié)語(yǔ)

在總結(jié)和分析現(xiàn)有城市河渠剖面流速分布規(guī)律的基礎(chǔ)上，為提高剖面流速的測(cè)量精度，提出了一種基于多普勒流量計(jì)的城市河渠剖面流速分析方法。

城市河道或明渠（寬深比15以上，剖面規(guī)則）剖面流速分布具有較明顯的規(guī)律，相對(duì)流速與相對(duì)水深在剖面內(nèi)區(qū)呈對(duì)數(shù)函數(shù)關(guān)系，在剖面外區(qū)為二次函數(shù)關(guān)系;垂線平均流速沿橫向分布理想情況接近于冪函數(shù)或拋物線函數(shù)分布，但考慮到現(xiàn)場(chǎng)環(huán)境大多數(shù)為非理想，且橫向流速變化關(guān)系大致都符合中間區(qū)域流速比河道邊壁流速大的規(guī)律。因此，可使用多點(diǎn)測(cè)量的方法來(lái)提高橫向平均流速的測(cè)量精度。通過(guò)將實(shí)測(cè)數(shù)據(jù)與現(xiàn)有的經(jīng)驗(yàn)公式進(jìn)行擬合驗(yàn)證，發(fā)現(xiàn)在一條垂線上單獨(dú)使用對(duì)數(shù)函數(shù)或二次函數(shù)分析所得出的結(jié)果誤差較大。通過(guò)對(duì)選取的玉溪市河道不同斷面進(jìn)行多次實(shí)際測(cè)量和數(shù)據(jù)分析，可確定內(nèi)外區(qū)的分界點(diǎn)相對(duì)水深一般位于y/H值為0.6-0.85的位置，此時(shí)相關(guān)系數(shù)可達(dá)到0.95以上。該法可用于城市河道明渠剖面流速實(shí)際監(jiān)測(cè)中，配合多普勒流量計(jì)可實(shí)現(xiàn)河渠的流量在線監(jiān)測(cè)。

參考文獻(xiàn)：

[1]蔡甫款.明渠流量測(cè)量的關(guān)鍵技術(shù)研究[D].杭州：浙江大學(xué)，2003.

[2] 劉星，賀曉春.天然河道斷面流速分布計(jì)算[J].四川聯(lián)合大學(xué)學(xué)報(bào)（工程科學(xué)版），1998，2（4）：40-44.

[3]韓金旭，蔡大應(yīng)，李敬茹，等.明渠水流流層平均流速分布規(guī)律及其應(yīng)用[J].中國(guó)農(nóng)村水利水電，2015（6）：116-128.

[4]韓金旭，黃福貴，卞艷麗，等.明渠水流流速橫向分布拋物線律及其應(yīng)用[J].人民黃河，2013，35（5）： 83-85.

[5]王二平，金輝，王文君，等.矩形明渠流速分布律及流量自動(dòng)化測(cè)量方法.人民黃河，2008，30（3）：31-32.

[6]王二平，金輝，關(guān)靖.矩形明渠流速分布規(guī)律及流量自動(dòng)化量測(cè)方法[J].華北水利水電學(xué)院學(xué)報(bào)，2007，28（6）：13-16.

[7]王二平，金輝，張艷艷，等.矩形明渠流速分布特征及其在流量量測(cè)中的應(yīng)用[J-.灌溉排水學(xué)報(bào)，2008，27（4）：28.

[8]盧金友，詹義正，趙根生，等.河道中受側(cè)壁影響的斷面流速分布規(guī)律研究[J].水利學(xué)報(bào)，2012，43（6）：645-658.

[9] 胡云進(jìn)，萬(wàn)五一，蔡甫款，等.窄深矩形斷面明渠流速分布的研究[J].浙江大學(xué)學(xué)報(bào)（工學(xué)版），2008，42（1）：183-187.

10]孫東坡，王二平，董志慧，等.矩形斷面明渠流速分布的研究及應(yīng)用[J].水動(dòng)力學(xué)研究與進(jìn)展，2004，19（2）：144-150.

[11]J.F.沃克.用通用兩點(diǎn)法測(cè)量明渠流速[J].浙江水利科技，1989（1）：57-59.