不同坡度三種渠道下巴歇爾槽泥沙試驗(yàn)研究及數(shù)值模擬

曾鈺峰，何新林，李小龍，楊 廣，喬長錄，金 瑾

（石河子大學(xué)水利建筑工程學(xué)院，新疆 石河子 832000）

0 引 言

巴歇爾量水槽作為灌區(qū)使用最廣泛的量水槽，具有水頭損失小，上游雍水程度低等優(yōu)點(diǎn)。但其在淹沒出流，橫向沉降，縱向沉降［1］等條件下具有較大測流誤差。灌區(qū)常見的3種渠道分別為梯形渠道、矩形渠道、U 形渠道［2，3］。同時(shí)新疆河流含沙量較大，灌區(qū)輸水存在淤積情況［4］。所以探究巴歇爾量水槽在來沙下的適宜渠道及坡度，有利于提高實(shí)踐中的量水精度，促進(jìn)農(nóng)業(yè)水價(jià)改革。

國內(nèi)外學(xué)者對巴歇爾槽也做過一系列研究。KHOSRONEJAD 等［5］通過數(shù)值模擬和現(xiàn)場試驗(yàn)研究，量化了兩個(gè)現(xiàn)場尺度巴歇爾量水槽的誤差范圍和水面波動(dòng)。HEINER 等［6］通過試驗(yàn)發(fā)現(xiàn)在巴歇爾量水槽2∕3A 處進(jìn)行水頭測量時(shí)，誤差可高達(dá)60%，同時(shí)，安裝不正確及缺少翼壁設(shè)計(jì)會(huì)形成駐波，這也會(huì)影響流量測量的準(zhǔn)確性。KIM 等［7］研究發(fā)現(xiàn)巴歇爾量水槽在估算枯水期的流量時(shí)比水位測量站的狀態(tài)流量額定曲線方程更為精確。SAVAGE 等［8］通過物理實(shí)驗(yàn)和Flow-3D 模擬研究發(fā)現(xiàn)，當(dāng)在巴歇爾量水槽規(guī)定以外的位置進(jìn)行水位測量時(shí)，可能會(huì)出現(xiàn)高達(dá)60%的流量誤差。同時(shí)，不同的入口翼墻配置也影響了巴歇爾槽的測流精度。許虎等［9］對不同進(jìn)口連接段形式的巴歇爾槽進(jìn)行了數(shù)值模擬，結(jié)果表明直面過渡段的巴歇爾槽測流精度較高。黃智剛等［10］通過對拋物線形量水槽和巴歇爾槽的水頭損失對比分析，認(rèn)為巴歇爾槽在U 形渠道上會(huì)產(chǎn)生較大的水頭損失。THORNTON 等［11］將巴歇爾量水槽原始的測流方程進(jìn)行了修改，開發(fā)了獨(dú)特的方程用于測量超臨界和亞臨界流態(tài)的流量。前人對巴歇爾量水槽的自身幾何結(jié)構(gòu)和進(jìn)口連接段的研究進(jìn)行的較多，而缺少對巴歇爾槽適用性的研究。所以探究巴歇爾量水槽的適宜渠道及坡度對實(shí)踐具有重要意義。同時(shí)分析巴歇爾槽在各工況下的過沙能力及測流精度，有利于指導(dǎo)新疆的灌溉輸水及農(nóng)業(yè)水價(jià)改革。

1 試驗(yàn)方案

1.1 試驗(yàn)布置

試驗(yàn)在石河子大學(xué)水利大廳進(jìn)行，試驗(yàn)布置如圖1所示，試驗(yàn)系統(tǒng)包括：泵房、電磁流量計(jì)、調(diào)節(jié)閥門、穩(wěn)水池、閘門、U 形渠、矩形渠、梯形渠、巴歇爾槽、回水池、變坡控制系統(tǒng)等。水泵流量調(diào)節(jié)范圍為0～0.1 m3∕s，根據(jù)明渠堰槽流量計(jì)計(jì)量檢定規(guī)程JJG（水利）004-2015［12］及設(shè)計(jì)流量，選用喉道寬為0.152 m 的巴歇爾槽進(jìn)行試驗(yàn)。其測流范圍為0.001 5～0.1 m3∕s，符合設(shè)計(jì)要求。

圖1 試驗(yàn)布置圖（單位：m）Fig.1 Test layout

1.2 試驗(yàn)設(shè)計(jì)

查閱規(guī)范，0.152 m 喉道寬的巴歇爾槽臨界淹沒度為0.6。為保證試驗(yàn)為自由出流，先在未安裝巴歇爾槽的渠道上進(jìn)行水位～流量關(guān)系測定，試驗(yàn)最大流量為80 L∕s，梯形渠道、U 形渠道、矩形渠道0.1%坡度時(shí)的上游水深h上分別為：23.18，23.50，20.58 cm。即巴歇爾槽底檻最低高度P為［13］：

式中：P為底檻高度，cm；h上為未安裝巴歇爾槽時(shí)的渠道正常水深；ha為巴歇爾槽上游測點(diǎn)處水深。

所以為保證3種渠道均為自由出流，設(shè)置底檻高度為5 cm。L1=4P，L5=6P，渠道的具體設(shè)計(jì)尺寸及巴歇爾槽尺寸見表1及圖2～圖4。

表1 巴歇爾槽結(jié)構(gòu)尺寸圖mTab.1 Structural dimension drawing of Parshall flume

圖2 渠道及巴歇爾槽縱剖面圖Fig.2 Longitudinal section of channel and Parshall flume

圖3 渠道及巴歇爾槽俯視圖Fig.3 Top view of channel and Parshall flume

圖4 3種渠道橫截面圖（單位：m）Fig.4 Cross section of three channels

渠道為鋼結(jié)構(gòu)，糙率為0.012。流量控制采用流量控制系統(tǒng)（DCMS），如圖5所示。水位檢測采用超聲波探頭及雷達(dá)水位計(jì)，如圖6、圖7 所示。流速測量采用多普勒流速儀，如圖8 所示，變坡控制系統(tǒng)如圖9所示。

圖5 流量控制系統(tǒng)Fig.5 Flow control system

圖6 超聲波水位計(jì)Fig.6 Ultrasonic water level gauge

圖7 雷達(dá)水位計(jì)Fig.7 Radar water level gauge

圖8 多普勒流速儀Fig.8 Doppler current meter

圖9 變坡控制系統(tǒng)Fig.9 Slope control system

圖10 網(wǎng)格劃分示意圖Fig.10 Schematic diagram of grid division

泥沙試驗(yàn)參考瑪納斯河多年平均含沙量2.2 kg∕m3［14］，試驗(yàn)用沙取自瑪納斯河天然河沙，根據(jù)流量的不同，在穩(wěn)水池傾倒不同重量的河沙，攪拌均勻，具體試驗(yàn)?zāi)嗌持亓?，泥沙粒徑?jí)配分別如表2、表3所示。

表2 不同流量下試驗(yàn)?zāi)嗌硟A倒量Tab.2 Test sediment dumping under different flow rates

表3 天然河沙粒徑級(jí)配Tab.3 Grain size grading of natural river sand

2 數(shù)值模擬方法

2.1 控制方程

多相流模型采用歐拉模型，歐拉模型允許對多個(gè)獨(dú)立但相互作用的模型進(jìn)行建模。這些相可以是任何組合的液體、氣體或固體。每個(gè)相位的處理均使用歐拉法。歐拉模型適用于水沙兩相流的固液滑移流動(dòng)模擬。它將主相和次相均看作為一個(gè)連續(xù)相，通過壓力和相間交換系數(shù)實(shí)現(xiàn)耦合。流體—固體流動(dòng)間的交換模型采用Wen and Yu 模型，該模型適用于次相體積濃度明顯低于主相的情況［15］。

2.2 體積分?jǐn)?shù)方程

體積分?jǐn)?shù)表示每個(gè)相所占的空間，每個(gè)相分別滿足質(zhì)量和動(dòng)量守恒定律。

相位的體積Vq由以下公式定義［16］：

其中：

相q的有效密度表達(dá)式為：

式中：ρq為相q的物理密度。

2.3 質(zhì)量守恒連續(xù)性方程

2.4 動(dòng)量守恒方程

式中：μq和λq分別為相的剪切和體積黏度；為外部體積力；為浮托力；為虛擬質(zhì)量力；為相之間的相互作用力；→Fwl，q為壁面滑移力；為湍流擴(kuò)散力；為相間速度；Kpq為相間交換系數(shù)；p為兩相間共享的壓力。

2.5 網(wǎng)格劃分

結(jié)構(gòu)化網(wǎng)格具有計(jì)算效率高，但幾何適應(yīng)性差的特點(diǎn)。而非結(jié)構(gòu)化網(wǎng)格，幾何適應(yīng)性較好，但計(jì)算較慢［17］。經(jīng)過兩者權(quán)衡，梯形和矩形渠道采用結(jié)構(gòu)化網(wǎng)格，U形渠道采用結(jié)構(gòu)化和非結(jié)構(gòu)化網(wǎng)格的混合。網(wǎng)格質(zhì)量較好，均在0.8 以上，滿足計(jì)算精度要求。同時(shí)加密氣液交界處網(wǎng)格數(shù)量，提高對液面捕捉的精度。梯形、矩形、U形網(wǎng)格數(shù)量分別為127萬、105萬、138萬。

2.6 數(shù)值求解方法與邊界條件

湍流方程采用RNGk-ε雙方程模型［18，19］，算法采用基于壓力耦合的coupled 算法。使用瞬態(tài)算法，時(shí)間步長為0.001 s。邊界條件具體設(shè)置如圖11 所示，自由表面邊界采用壓力進(jìn)口，水相和沙相采用速度進(jìn)口，沙相占0.014 67，壁面采用標(biāo)準(zhǔn)無滑移壁面，材質(zhì)為鋼，出口為充分發(fā)展的自由出流。

圖11 邊界條件示意圖Fig.11 Schematic diagram of boundary conditions

3 結(jié)果與分析

3.1 模型驗(yàn)證

為驗(yàn)證模型可靠性，將59.881 456 L∕s流量0.1%坡度下3種渠道巴歇爾槽水位及流速的試驗(yàn)值與模擬值進(jìn)行對比。如圖12～15 所示，梯形渠道、U 形渠道、矩形渠道巴歇爾槽上游測點(diǎn)水位的試驗(yàn)與模擬相對誤差分別為3.90%、1.31%、4.23%，最大誤差不超過5%。由水面線可知，巴歇爾槽在收縮段為緩流，流經(jīng)喉道段，水流由緩流經(jīng)臨界流轉(zhuǎn)變成急流，數(shù)值模擬結(jié)果符合試驗(yàn)量水槽流態(tài)。梯形渠道、U形渠道、矩形渠道巴歇爾槽上游測點(diǎn)流速的試驗(yàn)與模擬相對誤差分別為1.94%、0.82%、2.08%，最大誤差不超過3%。即認(rèn)為該模型滿足模擬精度要求。

圖12 梯形渠道巴歇爾槽水面線數(shù)值模擬與試驗(yàn)對比圖Fig.12 Comparison between numerical simulation and test of water surface profile of trapezoidal channel and Parshall flume

圖13 U形渠道巴歇爾槽水面線數(shù)值模擬與試驗(yàn)對比圖Fig.13 Comparison between numerical simulation and test of water surface profile of U-shaped channel and Parshall flume

圖14 矩形渠道巴歇爾槽水面線數(shù)值模擬與試驗(yàn)對比圖Fig.14 Comparison between numerical simulation and test of water surface profile of rectangular channel and Parshall flume

圖15 3種渠道巴歇爾槽流速數(shù)值模擬值與試驗(yàn)對比圖Fig.15 Comparison of numerical simulation values and test results of flow velocity in three channels

3.2 泥沙分布

圖16 為59.881 456 L∕s 流量下3 種渠道巴歇爾槽分別在坡度0.1%、0.4%、0.8%的泥沙分布模擬圖。圖16中T、U、J分別表示梯形渠道、U 形渠道、矩形渠道。1、4、8 分別表示坡度0.1%、0.4%、0.8%。由圖16 中可知梯形渠道泥沙相對集中在渠道上游與巴歇爾槽連接段底部兩側(cè)及下游渠道底部兩側(cè)，U 形渠道泥沙相對集中在渠道上游中間底部、巴歇爾槽上游連接段及渠道下游中間底部，矩形渠道泥沙分布相對比較分散。從整體來看3種渠道泥沙分布主要位于巴歇爾槽上游連接段位置。而在3 種渠道巴歇爾槽槽身泥沙分布中，梯形渠道巴歇爾槽泥沙積累最少，U形渠道巴歇爾槽泥沙積累最多。隨著坡度的增加，渠道及巴歇爾槽泥沙含量開始逐漸減少，淤積程度逐漸降低。

圖16 59.881 456 L/s流量下3種渠道各坡度巴歇爾槽泥沙分布模擬圖Fig.16 Simulation diagram of sediment distribution in Parshall flume of three channels at the flow rate of 59.881 456 L/s

3.3 測流精度

由于本試驗(yàn)均是在自由出流情況下進(jìn)行，根據(jù)《JJG（水利）004-2015 明渠堰槽流量計(jì)計(jì)量檢定規(guī)程》，巴歇爾槽為單一的水位流量關(guān)系，其流量計(jì)算公式如下［12］：

式中：Q為流量，m3∕s；h為上游實(shí)測水頭，m。

3 種渠道巴歇爾槽在0.1%～0.8%坡度下各流量的測流相對誤差如圖17 所示?？傮w來講，3 種渠道巴歇爾槽的測流誤差為：U形渠道＞矩形渠道＞梯形渠道。隨著坡度的增加，測流誤差逐漸增大；隨著流量的增大，測流誤差逐漸減小。在小流量5.668 975 L∕s時(shí)，3種渠道在0.1%坡度～0.3%坡度測流誤差呈下降趨勢，在0.3%坡度～0.8%坡度測流誤差呈上升趨勢。在0.02～0.04 m3∕s時(shí)，巴歇爾槽測流誤差下降較為明顯。0.02～0.04 m3∕s 對應(yīng)的巴歇爾槽水位為15.5～24 cm，而巴歇爾槽的測量水位范圍為3～45 cm，即在巴歇爾槽最大測量水位的34.44%～53.33%時(shí)，巴歇爾槽的測流誤差下降較快，而超過53.33%后，測流誤差下降較慢，最后曲線趨于平穩(wěn)，測流相對穩(wěn)定，誤差較小。5.668 975 L∕s 對應(yīng)的巴歇爾槽水深為6.97 cm，即當(dāng)巴歇爾槽的水位小于最大測量水位的15.49%時(shí)，測流誤差相對較大。3 種渠道在0.1%～0.4%坡度下的測流誤差相對較小，大部分在規(guī)范規(guī)定的5%以內(nèi)。而在0.4%～0.8%坡度下的測流誤差相對較大，大部分工況逐漸超出5%。

圖17 0.1%～0.8%坡度下各流量3種渠道巴歇爾槽測流相對誤差Fig.17 The relative error of flow measurement in three channels with 0.1% ～ 0.8% gradient

3.4 流速分布

圖18 為沿渠道方向各斷面的位置示意圖，斷面1～14 的相對位置分別為0，1.5，3，3.2，3.403，3.606，3.81，3.912，4.014，4.115，4.318，4.521，5.025，6.325m。其中4～6為巴歇爾槽收縮段斷面，7～9 為喉道段斷面，10～12 為擴(kuò)散段斷面。分析各個(gè)斷面的平均流速，得出如圖17 所示的各斷面沿Z方向的流速變化圖。圖19所示分別為流量5.668 975、59.881 456 L∕s下3種渠道0.1%～0.8%坡度的流速變化。在小流量5.668 975 L∕s下，3種渠道流速相近，巴歇爾槽流速在喉道段前1∕3處（即斷面8處）達(dá)到最大值，由于流量較小，水流在擴(kuò)散段流速較小，導(dǎo)致水流在喉道段末雍水，流速驟降，且泥沙易淤積在喉道末與擴(kuò)散段的連接處，對上游收縮段流速有一定影響，這也是低流量下測流不準(zhǔn)的原因之一。在大流量59.881 456 L∕s 下，3 種渠道流速相差較大，巴歇爾槽的流速最大值發(fā)生在擴(kuò)散段，由于流量較大，水流沖刷較大，在喉道末處雍水現(xiàn)象減弱。從總體來看，隨著坡度的增加，流速逐漸增大。在0.1%坡度～0.4%坡度下，3種渠道巴歇爾槽的流速差距相對不大，而坡度超過0.4%后，流速差距增加，U 形渠道＞矩形渠道＞梯形渠道，流態(tài)越來越急，臨界流發(fā)生位置提前，導(dǎo)致槽內(nèi)流速波動(dòng)加劇，泥沙顆粒碰撞加劇，測流精度降低。

圖18 斷面位置圖Fig.18 Section location map

圖19 5.668 975及59.881 456 L/s流量下沿渠道方向各斷面流速分布Fig.19 Velocity distribution of each section along the channel under the flow of 5.668 975 and 59.881 456 L/s

3.5 挾沙能力

1958 年，張瑞瑾［20］等根據(jù)懸移質(zhì)具有制紊作用的觀點(diǎn)，從挾沙水流的能量平衡原理出發(fā)，建立能量平衡方程，導(dǎo)出了水流挾沙能力公式的結(jié)構(gòu)形式：

式中：S*是懸移質(zhì)水流挾沙能力，kg∕m3；U為斷面平均流速，m∕s；R為水力半徑，m，對于寬淺河道R可近似用斷面平均水深h代替；ω為泥沙沉降速度，m∕s；g為重力加速度，m∕s2；K和m分別為水流挾沙能力系數(shù)和指數(shù)，這些系數(shù)可以通過半理論表達(dá)式和實(shí)測資料確定［21］，K取0.02，m取0.92。

此公式適用于含沙量0.1～100 kg∕m3［18］，其中懸移質(zhì)泥沙粒徑不大于0.062 mm時(shí)，沉速采用斯托克斯公式［20］：

式中：γs為泥沙的干容重，N∕m3；γ為水的容重，N∕m3；d取泥沙的中值粒徑，m；v為運(yùn)動(dòng)黏滯系數(shù)，m2∕s。

當(dāng)粒徑為0.062～2.0 mm 時(shí)，采用沙玉清的過渡區(qū)沉速公式間接推算泥沙顆粒的沉降粒徑：

式中：Sa為沉速判數(shù)；φ為粒徑判數(shù)；ρs為泥沙密度，g∕cm3，ρω為清水密度，g∕cm3；ω沉降速度，cm∕s；d為中值粒徑，mm；g為重力加速度，cm∕s2；ν為黏滯系數(shù)，cm2∕s。

圖20、圖21 分別為0.1%坡度，流量5.668 975、59.881 456 L∕s下斷面1～14的沿程挾沙力圖。由圖中可知，5.668 975 L∕s流量下巴歇爾槽的最大挾沙力位于喉道段前1∕3處（斷面8），最大值為0.28 kg∕m3；巴歇爾槽最小挾沙力位于擴(kuò)散段進(jìn)口處（斷面10），最小值為0.006 kg∕m3。59.881 456 L∕s 流量下巴歇爾槽最大挾沙力位于擴(kuò)散段后1∕3 處（斷面13），最大值為6.04 kg∕m3；巴歇爾槽最小挾沙力位于巴歇爾槽進(jìn)口，最小值為0.69 kg∕m3。流量越大，挾沙力越大。在大流量的情況下，巴歇爾槽大部分部位的挾沙力均大于2.2 kg∕m3，說明巴歇爾槽在大流量情況下，不存在淤積現(xiàn)象。

圖20 5.668 975 L/s流量0.1%坡度挾沙力沿程分布圖Fig.20 5.668 975 L/s flow 0.1% slope sediment carrying capacity along the road

圖21 59.881 456 L/s流量0.1%坡度挾沙力沿程分布圖Fig.21 59.881 456 L/s flow 0.1% slope sediment carrying capacity along the road

4 討 論

通過計(jì)算流體動(dòng)力學(xué)——CFD 對3 種渠道上巴歇爾槽0.1%～0.8%坡度的水沙試驗(yàn)進(jìn)行了數(shù)值模擬。通過研究發(fā)現(xiàn)，歐拉模型對水沙模擬具有較高的擬合性，更加貼近實(shí)際情況，能夠滿足精度要求。相較于DPM 模型、混合模型，歐拉模型將所有相都看作連續(xù)相，且相與相之間互相耦合作用。而DPM 模型，將顆?？醋麟x散相，相與相之間的交互作用較弱，不適用于河流泥沙模擬。混合模型相當(dāng)于簡化了多相流，雖然其具有計(jì)算較快，收斂較快的特點(diǎn)，但是對于復(fù)雜流動(dòng)，精度較低。

對于巴歇爾槽的研究，本文彌補(bǔ)了前人對巴歇爾槽泥沙及不同坡度不同渠道下研究的缺失。研究發(fā)現(xiàn)，梯形渠道在泥沙情況下具有較高精度。分析是由于梯形渠道將大部分泥沙攔截在巴歇爾槽上游連接段，且梯形渠道泥沙集中在渠底兩側(cè)，對水流影響較小，水與沙的摻互較少，對水流的擾動(dòng)較小。而U形渠道泥沙集中在中間底部，泥沙易排放于巴歇爾槽，對測流精度有一定影響，為含沙量比較高的灌區(qū)提供了理論依據(jù)。同時(shí)，巴歇爾槽的修建應(yīng)該在0.1%～0.4%坡度之間，坡度過大，對巴歇爾槽內(nèi)水流的流態(tài)，流速均會(huì)產(chǎn)生影響，水位處于不穩(wěn)定狀態(tài)。在設(shè)計(jì)流量時(shí)，易將流量設(shè)計(jì)在巴歇爾槽最大測量水位的30%以上，能夠保證測流精度。

但對于不同粒徑的泥沙還有待進(jìn)一步試驗(yàn)與分析，缺少巴歇爾槽對不同粒徑泥沙的過流能力分析，作者將繼續(xù)后續(xù)的研究。

5 結(jié) 論

（1） 在泥沙情況下，梯形渠道的巴歇爾槽具有更高的測流精度和適應(yīng)性，槽身泥沙淤積較少。

（2） 0.1%～0.4%坡度下的巴歇爾槽能夠較好滿足測流精度，超過0.4%坡度后測流誤差逐漸增大。

（3） 在巴歇爾槽最大測量水位的34.44%～53.33%時(shí)，巴歇爾槽的測流誤差下降較快，而超過53.33%后，測流誤差下降較慢，最后曲線趨于平穩(wěn)，測流相對穩(wěn)定，誤差較小。建議在實(shí)踐設(shè)計(jì)中，巴歇爾槽的水位保持在最大水位的30%以上。

（4） 為減少巴歇爾槽上游泥沙淤積，可適當(dāng)提高渠道坡度，但建議不超過0.4%。