基于ＶＯＦ?ＤＰＭ 模型的明渠含沙數(shù)值模擬及量水方案比選

摘 要：為探究含沙水流對(duì)明渠流量測(cè)量精度的影響，基于Ｆｌｕｅｎｔ 軟件采用ＶＯＦ－ＤＰＭ 模型對(duì)明渠進(jìn)行數(shù)值模擬，分別模擬不同進(jìn)口流速和含沙量組合條件下渠道內(nèi)流量、液位及泥沙顆粒分布，并將數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比分析，計(jì)算出各量水設(shè)施的流量誤差與液位誤差。結(jié)果表明：采用ＶＯＦ－ＤＰＭ 模型可以較好地模擬含沙水流在明渠的流態(tài)，通過(guò)開(kāi)展不同流量級(jí)的流量測(cè)試可知，在不同含沙量條件下，巴歇爾槽流量計(jì)誤差最?。ǎ常ァ保ィ?，超聲流量計(jì)（－４％～３％）和電磁流量計(jì)（ －９％～２％）量測(cè)效果均較好，雷達(dá)流量計(jì)誤差較大（－２０％～４％），說(shuō)明非接觸式流量計(jì)在高含沙量、大流量的情況下需進(jìn)一步提高測(cè)量精度；超聲液位、電磁液位與雷達(dá)液位誤差均較大，試驗(yàn)中尾部巴歇爾槽壅水以及泥沙在渠道底部淤積是造成較大液位誤差的原因之一。

關(guān)鍵詞：數(shù)值模擬；ＶＯＦ－ＤＰＭ 模型；含沙量；測(cè)流精度

中圖分類號(hào)：ＴＶ９３ 文獻(xiàn)標(biāo)志碼：Ａ ｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２４．１１．０２０

引用格式：唐淵，馬文波，溫紅艷，等．基于ＶＯＦ－ＤＰＭ 模型的明渠含沙數(shù)值模擬及量水方案比選［Ｊ］．人民黃河，２０２４，４６（１１）：１２６－１３２，１４８．

０ 引言

灌區(qū)量水是灌區(qū)水資源優(yōu)化配置和現(xiàn)代化管理的基本手段，是提高灌溉水有效利用系數(shù)的重要途徑［１］ 。目前的量水方法包括建筑物和設(shè)備型量水兩大類。各種量水方法都有其特定的適用范圍和條件，例如：巴歇爾槽量水精度高、適用范圍大、觀測(cè)方便，但結(jié)構(gòu)復(fù)雜、造價(jià)較高，且會(huì)影響渠道過(guò)水能力；超聲流量計(jì)測(cè)量精度高、操作簡(jiǎn)便、安裝方便、智能化程度高，但水位變化較大、水深較深、流態(tài)復(fù)雜的斷面布設(shè)多層測(cè)速聲路成本較高；電磁流量計(jì)精度高、測(cè)量范圍廣、壓力損失小、結(jié)構(gòu)簡(jiǎn)單、易維護(hù)，但測(cè)量精度易受外部磁場(chǎng)干擾，同時(shí)對(duì)液體介質(zhì)的導(dǎo)電性有一定要求；雷達(dá)流量計(jì)適用范圍廣、抗干擾能力強(qiáng)，不受流體性質(zhì)、溫度、壓力等因素影響，但安裝和使用復(fù)雜且成本較高。因此，應(yīng)根據(jù)灌區(qū)基本情況、量水精度要求、渠道水力邊界條件，確定相應(yīng)的量水方法和量水設(shè)施具體型號(hào)［２］ 。在明渠測(cè)流中，水體含沙是一個(gè)關(guān)鍵問(wèn)題，水體中的泥沙會(huì)對(duì)量水設(shè)備的測(cè)量精度產(chǎn)生不同程度的影響，因此開(kāi)展泥沙對(duì)量水設(shè)備測(cè)量精度的影響研究、優(yōu)化量水設(shè)備性能、提升測(cè)量精度有重要意義。

１ 概述

國(guó)內(nèi)外學(xué)者過(guò)去主要集中在清水條件下對(duì)量水設(shè)備測(cè)流精度進(jìn)行研究，關(guān)于泥沙或水中雜質(zhì)對(duì)測(cè)流設(shè)備潛在影響的研究相對(duì)較少，例如：范恬等［３］ 利用橫向擺桿提出一種基于流速面積法的明渠測(cè)流方法，結(jié)果表明該測(cè)流方法的測(cè)流精度符合取水計(jì)量的規(guī)定且橫向擺桿裝置對(duì)流場(chǎng)的擾動(dòng)較小；鮑子云等［４］ 將不同流量測(cè)箱與標(biāo)準(zhǔn)流量計(jì)在不同工況條件下的測(cè)量結(jié)果進(jìn)行對(duì)比分析，結(jié)果表明超聲、壓力液位流量測(cè)箱能在不同工況下工作且精度較高；萬(wàn)里等［５］ 設(shè)計(jì)了一種基于彈簧形變量與渠道過(guò)水?dāng)嗝嫠矔r(shí)流量之間關(guān)系的明渠測(cè)流裝置，擬合得到彈簧板式測(cè)流裝置的流量公式，計(jì)算流量與試驗(yàn)時(shí)的實(shí)測(cè)流量吻合；王寶賀等［６］ 提出了一種非接觸式流量計(jì)算方法，結(jié)果表明基于該流量系數(shù)的計(jì)算流量與實(shí)測(cè)流量十分接近；Ｚｈｏｕ 等［７］ 研究發(fā)現(xiàn)提高電磁流量計(jì)傳感器磁場(chǎng)激勵(lì)的均勻性可以提高電磁流量計(jì)的測(cè)量精度；王寶賀等［８］ 研究發(fā)現(xiàn)ＧＡＢＰ神經(jīng)網(wǎng)絡(luò)僅需４ 個(gè)特定位置（測(cè)點(diǎn)）流速即可得到滿足精度的計(jì)算結(jié)果，實(shí)現(xiàn)對(duì)明渠測(cè)點(diǎn)數(shù)量的優(yōu)化，保證測(cè)流精度；Ｓｈｉ 等［９］ 研究提出了一種基于ＡＩＮ 壓電微機(jī)械超聲換能器（ＰＭＵＴｓ）的高精度雙通道超聲流量計(jì)，用于測(cè)量小口徑管道中的流量。近年來(lái)，人們開(kāi)始重點(diǎn)研究泥沙或水中雜質(zhì)對(duì)設(shè)備測(cè)流的影響，例如：段炎沖等［１０］ 研究了含沙量對(duì)明渠超聲時(shí)差法測(cè)流精度的影響，結(jié)果表明超聲時(shí)差法流量計(jì)的測(cè)量值偏大，且偏差隨含沙量增大呈增大趨勢(shì)。

Ｆｌｕｅｎｔ 是一個(gè)用于模擬和分析在復(fù)雜幾何區(qū)域內(nèi)的流體流動(dòng)的專用ＣＦＤ 軟件［１１］ 。在明渠水流動(dòng)模擬中通常使用的多相流模型是ＶＯＦ 模型，模擬水中泥沙顆粒則需要借助離散相模型（Ｄｉｓｃｒｅｔｅ Ｐｈａｓｅ Ｍｏｄｅｌ）在Ｌａｇｒａｎｇｉａｎ 坐標(biāo)下模擬流場(chǎng)中離散的第三相。本文基于Ｆｌｕｅｎｔ 軟件和ＶＯＦ－ＤＰＭ 模型對(duì)明渠進(jìn)行數(shù)值模擬的結(jié)果（不同入口流速和含沙量組合條件下渠道內(nèi)流量、液位及泥沙顆粒分布等），將數(shù)值模擬數(shù)據(jù)與試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比分析，計(jì)算得出各量水設(shè)施的流量誤差與液位誤差，分析其分布規(guī)律，以期為明渠含沙水流條件下量水裝置的選擇提供技術(shù)支撐。

２ ＶＯＦ－ＤＰＭ 模型

ＶＯＦ－ＤＰＭ 模型依賴于兩種或多種流體互不滲透這一前提，對(duì)模型中的每一相引入體積分?jǐn)?shù)，并以此作為求解變量。在每個(gè)控制體積內(nèi)，所有相的體積之和等于１。全部特性參數(shù)和變量的場(chǎng)由各相共享，且這些場(chǎng)的值代表了整個(gè)系統(tǒng)內(nèi)體積平均值［１２］ 。根據(jù)控制體積內(nèi)各相的體積分?jǐn)?shù)，控制體積被賦予適當(dāng)?shù)奶匦院妥兞恐怠?/p>

１）體積分?jǐn)?shù)方程：

式中：ρｑ 和αｑ 分別為第ｑ 相的密度和體積分?jǐn)?shù)（ｑ ＝１，２，…，ｎ）；ｔ 為時(shí)間；ｕ 為混合流體速度；?為散度算子；ｎ 為多相流中流體種類數(shù)目，即相數(shù)；Ｓｑ 為源項(xiàng)，可以通過(guò)用戶自定義方式指定，缺省條件下為０；ｍｐ１ｑ 為從相ｐ１ 到相ｑ 的傳質(zhì)；ｍｑｐ１為從相ｑ 到相ｐ１ 的傳質(zhì)。

２）動(dòng)量方程：

式中：ρ 為混合流體密度，μ 為混合流體動(dòng)力黏度，ｐ２為混合流體壓力，ｇ 為重力加速度，ＦＳＶ 為表面張力的等價(jià)體積力。

３）能量方程（密度ρ、有效導(dǎo)熱系數(shù)ｋｅｆｆ、有效黏度τｅｆｆ各相共用，按各相的體積平均計(jì)算）：

式中：Ｓｈ 為源項(xiàng)（包括熱輻射及其他體積熱源的貢獻(xiàn)），Ｊｊ 為組分ｊ 的擴(kuò)散通量，Ｅ 為各相能量的質(zhì)量加權(quán)平均能量，Ｔ 為溫度，ｈｊ 為組分ｊ 的焓。

離散相模型（Ｄｉｓｃｒｅｔｅ Ｐｈａｓｅ Ｍｏｄｅｌ， ＤＰＭ）屬于歐拉－拉格朗日方法的一種，用于求解帶有顆粒（泥沙顆粒、水滴或氣泡）的流體流動(dòng)問(wèn)題。由于不考慮顆粒與顆粒之間的相互作用，因此一般用于顆粒體積濃度低于１０％、顆粒外觀較清晰，顆粒與流體的相互作用比較明確的場(chǎng)合。顆粒相控制方程為

式中：ｐ 為顆粒的參數(shù)，ｍｐ 為顆粒質(zhì)量，ｕｐ 為顆粒速度，ＦＤ 為曳力，ＦＢ 為重力造成的浮力，Ｆ 為顆粒受到的除曳力和重力造成的浮力外的其他力，ＦＶＭ 為虛擬質(zhì)量力，ＦＰ 為壓力梯度力，ＦＲ 為旋轉(zhuǎn)系統(tǒng)中存在的哥氏力和離心力，ＦＭ 為馬格努斯升力，ＦＳ 為薩夫曼升力，ＦＢＡ為Ｂａｓｓｅｔ 力。

３ 試驗(yàn)研究

３．１ 試驗(yàn)系統(tǒng)

本次測(cè)流試驗(yàn)在室內(nèi)自循環(huán)明渠試驗(yàn)水槽中進(jìn)行，水槽由入口消力池、矩形渠道、尾水池和輸水管道４ 個(gè)部分構(gòu)成（見(jiàn)圖１）。入口消力池可以使下泄急流迅速變?yōu)榫徚?，消力池之后為矩形渠道（渠道長(zhǎng)３３ ｍ、寬０．６ ｍ、深０．７ ｍ，坡降為１‰）。為保證水流流動(dòng)穩(wěn)定，最上游的箱式電磁流量計(jì)布置在距入口１４ ｍ 處，箱式電磁流量計(jì)之后沿渠道每隔２ ｍ 依次布置超聲流量計(jì)與雷達(dá)流量計(jì)。在明渠末端布設(shè)標(biāo)準(zhǔn)小型５ 號(hào)巴歇爾槽，矩形渠道之后為尾水池，尾水由離心泵通過(guò)輸水管道回流至入口消力池。為減少泥沙在池中淤積，尾水池中放置沖沙泵來(lái)清理泥沙?；厮艿罏橹睆剑埃?ｍ 的圓管，管道上布設(shè)內(nèi)徑為０．４ ｍ的管道式電磁流量計(jì)，經(jīng)檢測(cè)在清水條件下示值誤差為０．４２％、測(cè)量精度０．５ 級(jí)。研究發(fā)現(xiàn)管道式電磁流量計(jì)不受介質(zhì)及雜質(zhì)的影響并可保證高精度量測(cè)，因此可將清水標(biāo)定結(jié)果應(yīng)用于渾水流量測(cè)試［１２］（試驗(yàn)通過(guò)變頻器控制離心泵的轉(zhuǎn)速，從而控制和調(diào)節(jié)流量）。

３．２ 試驗(yàn)步驟

為了較全面地分析泥沙對(duì)明渠量水設(shè)備的影響，開(kāi)展了不同含沙量下不同流量級(jí)的測(cè)量試驗(yàn)。試驗(yàn)采用寧夏引黃灌區(qū)的天然沙樣作為試驗(yàn)用沙，試驗(yàn)用沙粒徑為０．０１９ ｍｍ。在此基礎(chǔ)上進(jìn)行了５ 個(gè)量級(jí)含沙量試驗(yàn)，含沙量Ｓ 為０（清水）、５、１０、１５、２０ ｋｇ／ ｍ３，每級(jí)含沙量下均開(kāi)展５ 級(jí)流量（Ｑ≈３７０、５５０、７１０、８７０、１ ０１０ ｍ３ ／ ｈ）測(cè)試，利用水泵變頻器調(diào)節(jié)頻率來(lái)實(shí)現(xiàn)５個(gè)流量變化（管道式電磁流量計(jì)為滿管流，將各瞬時(shí)實(shí)測(cè)流量等價(jià)為相應(yīng)入口流速進(jìn)行后續(xù)Ｆｌｕｅｎｔ 軟件速度入口設(shè)置）。

為了便于試驗(yàn)，按以下步驟進(jìn)行：１） 清水試驗(yàn)。通過(guò)調(diào)節(jié)水泵變頻器頻率實(shí)現(xiàn)５ 個(gè)不同量級(jí)流量變化，每級(jí)流量設(shè)定試驗(yàn)時(shí)間６０ ｍｉｎ。２） 含沙水流試驗(yàn)。依據(jù)之前設(shè)計(jì)好的４ 個(gè)量級(jí)含沙量（不含清水），將提前稱量好的泥沙分批次倒入水池中，待泥沙攪拌均勻后開(kāi)始測(cè)流，依次調(diào)節(jié)水泵變頻器頻率完成不同流量級(jí)的試驗(yàn)。３）重復(fù)步驟２）完成不同含沙量、流量級(jí)的試驗(yàn)。每個(gè)流量級(jí)設(shè)定試驗(yàn)時(shí)間６０ ｍｉｎ，每級(jí)沙含量對(duì)應(yīng)５ 個(gè)流量級(jí)，共計(jì)２０ 組試驗(yàn)。計(jì)算機(jī)自動(dòng)采集各設(shè)備的實(shí)時(shí)水位、流量、流速等數(shù)據(jù)，采樣頻率為２ ｓ 一次。含沙水流試驗(yàn)每級(jí)流量均進(jìn)行含沙水樣采集，以計(jì)算水體含沙量。尾水池的沖沙泵始終保持運(yùn)行以保證不同流量級(jí)泥沙含量數(shù)據(jù)的一致性。

４ 數(shù)值模擬

４．１ 幾何建模

通過(guò)試驗(yàn)資料可查得試驗(yàn)斷面幾何形狀為矩形渠道，渠道高為０．７ ｍ、寬為０．６ ｍ、長(zhǎng)為１４ ｍ。根據(jù)《明渠實(shí)流法流量比對(duì)現(xiàn)場(chǎng)檢測(cè)規(guī)程》（Ｔ／ ＣＩＤＡ ００１４—２０２２）的規(guī)定，上游順直段長(zhǎng)度應(yīng)大于被測(cè)斷面最大液位的２０ 倍，由實(shí)測(cè)數(shù)據(jù)得被測(cè)斷面最大液位（水深）為０．６２３ ｍ，由此可知渠道長(zhǎng)?。保?ｍ 滿足水流流態(tài)穩(wěn)定的要求；入口管道為直徑０．４ ｍ 的圓管，管道內(nèi)為滿管流，因此將其簡(jiǎn)化為一個(gè)面。

４．２ 網(wǎng)格劃分

本次模擬使用ＡＮＳＹＳ 的ｍｅｓｈ 軟件進(jìn)行網(wǎng)格劃分，并對(duì)幾何的各部分進(jìn)行定義。本次數(shù)值模擬對(duì)象為矩形渠道，結(jié)構(gòu)簡(jiǎn)單，使用結(jié)構(gòu)化網(wǎng)格劃分，提高精度，減少網(wǎng)格數(shù)量，提高計(jì)算效率。采用Ｍｕｌｔｉｚｏｎｅ 方法進(jìn)行均勻化網(wǎng)格劃分，網(wǎng)格數(shù)量一共為１３３ ７５６ 個(gè)，網(wǎng)格質(zhì)量較好。

４．３ Ｆｌｕｅｎｔ 計(jì)算

為了使結(jié)果更貼合實(shí)際，采用Ｆｌｕｅｎｔ 軟件模擬明渠內(nèi)部流場(chǎng)及外部特性。將模型導(dǎo)入Ｆｌｕｅｎｔ 軟件中，設(shè)置多相流模型為ＶＯＦ 模型（加水和空氣），并添加離散相模型（研究氣液固三相流），含沙量分別為０、５、１０、１５、２０ ｋｇ／ ｍ３，顆粒直徑為０．０１９ ｍｍ，其分別對(duì)應(yīng)的體積分?jǐn)?shù)見(jiàn)表１。

湍流模型設(shè)為標(biāo)準(zhǔn)ｋ－ｅｐｓｉｌｏｎ 模型，此模型在滿足計(jì)算效率的前提下滿足相應(yīng)的計(jì)算精度，廣泛應(yīng)用于黏性模擬，壁面函數(shù)選取標(biāo)準(zhǔn)壁面函數(shù)（ＳＷＦ）。定義邊界條件，流體在豎直ｙ 方向的重力加速度為９．８１ ｍ／ ｓ２，進(jìn)口設(shè)置為速度入口，入口速度根據(jù)實(shí)測(cè)數(shù)據(jù)分別設(shè)定，離散相模型的壁面條件（Ｂｏｕｎｄａｒｙ Ｃｏｎｄｉｔｉｏｎ Ｔｙｐｅ）選擇ｅｓｃａｐｅ 顆粒逃逸；上表面設(shè)置為壓力進(jìn)口，壓力為一個(gè)大氣壓；出口設(shè)置為壓力出口，渠道邊壁設(shè)置為標(biāo)準(zhǔn)固定無(wú)滑移壁面，粗糙度常數(shù)設(shè)置為０．０１３。求解算法設(shè)置為基于壓力基的ＰＩＳＯ 算法，其包括預(yù)測(cè)—修正—再修正３ 個(gè)步驟，在本次模擬明渠含沙水流瞬態(tài)三維流動(dòng)過(guò)程中，無(wú)論計(jì)算精度還是計(jì)算時(shí)間都具有一定優(yōu)勢(shì)。空間離散中的壓力差值采用Ｂｏｄｙ Ｆｏｒｃｅ Ｗｅｉｇｈｔｅｄ（體積力加權(quán)壓力）插值，其余差值均使用精度更高的二階迎風(fēng)格式，從全局對(duì)其進(jìn)行初始化設(shè)置：步數(shù)為６ ０００ 步，時(shí)間步長(zhǎng)為０．０５ ｓ。

５ 結(jié)果分析

５．１ 同一含沙量不同流量

以含沙量１０ ｋｇ／ ｍ３為例，分別對(duì)實(shí)測(cè)入口速度０．８１６、１．１９９、１．５６０、１．９１６、２．２３０ ｍ／ ｓ 進(jìn)行數(shù)值模擬，并以速度１．１９９、１．９１６ ｍ／ ｓ 為例研究渠道水體積分?jǐn)?shù)、流速、壓強(qiáng)及泥沙濃度的變化（見(jiàn)圖２）。

由圖２（ａ）可知，紅藍(lán)交界處為水氣交界面，水面上方藍(lán)色區(qū)域?yàn)榭諝猓嫦路郊t色區(qū)域?yàn)樗?，隨著入口流速的增大，穩(wěn)定后的液位上升，模擬水流流態(tài)符合實(shí)際水流流態(tài)，水面線平穩(wěn)，水位與實(shí)際水位一致；由圖２（ｂ）可知，入口后方水低流速區(qū)面積增大，渦流現(xiàn)象加??；由圖２（ｃ）（ｅ）可知，隨著入口流速的增大，水中泥沙濃度下降，底部泥沙濃度比上方泥沙濃度大，入口處泥沙沉積較多；由圖２（ｄ）（ｆ）可知，隨著入口流速的增大，底部壓強(qiáng)增大，渠底壓強(qiáng)分布較為均勻。以標(biāo)準(zhǔn)表流量與巴歇爾槽液位為基準(zhǔn)，將模擬流量、液位等數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比（見(jiàn)表２）；同時(shí)，對(duì)含沙量０、１０、２０ ｋｇ／ ｍ３三種工況計(jì)算不同量水設(shè)備實(shí)測(cè)流量誤差與液位誤差（見(jiàn)圖３）。

由圖３（ａ）（ｃ）（ｅ）可知，在含沙量一定的情況下，隨著入口流速的增大，超聲流量計(jì)流量測(cè)量誤差在±３％以內(nèi)，量測(cè)精度較高；電磁流量計(jì)測(cè)量誤差在－９％～２％范圍內(nèi)波動(dòng)，量測(cè)效果較好；雷達(dá)流量計(jì)誤差則較大，誤差為－２０％～４％，且隨著入口流速的增大呈逐漸遞增的趨勢(shì)，這說(shuō)明非接觸式流量計(jì)在高含沙量、大流量的情況下需要進(jìn)一步提高測(cè)量精度；巴歇爾槽流量計(jì)總體誤差最小，誤差為－３％ ～１％，量測(cè)效果較好，驗(yàn)證了傳統(tǒng)量水堰槽的可靠性。

由圖３（ｂ）（ｄ）（ｆ）可知，在含沙量一定的情況下，隨著入口流速的增大，超聲液位、電磁液位與雷達(dá)液位總體誤差均較大，電磁液位誤差為－１２％ ～ －４％，超聲液位與雷達(dá)液位誤差為－２３％ ～ －８％。試驗(yàn)中尾部巴歇爾槽壅水以及泥沙在渠道底部的淤積是造成較大液位誤差的原因之一。

５．２ 同一流量不同含沙量

由于實(shí)測(cè)數(shù)據(jù)是基于不同含沙量分批進(jìn)行試驗(yàn)所得，無(wú)法在同一時(shí)刻對(duì)多種含沙水流工況進(jìn)行試驗(yàn)，因此此處以近似流量７１０ ｍ３ ／ ｈ 為例，分別對(duì)含沙量０、５、１０、１５、２０ ｋｇ／ ｍ３ 進(jìn)行數(shù)值模擬， 并以含沙量５、１５ｋｇ／ ｍ３為例研究渠道水體積分?jǐn)?shù)、流速、壓強(qiáng)及泥沙濃度的變化，如圖４ 所示。

由圖４（ａ）可知，在初始時(shí)刻入口處水流受到重力的作用使得渠道內(nèi)水的體積分?jǐn)?shù)銳減，水面線逐漸下降并產(chǎn)生略微波動(dòng)；隨著水流繼續(xù)在渠道中流動(dòng)水面線逐漸趨于平穩(wěn)，含沙量為５、１５ ｋｇ／ ｍ３的水面線在穩(wěn)定后高度一致。由圖４（ｂ）可知，含沙量為５、１５ ｋｇ／ ｍ３ 的渠道內(nèi)速度分布較為接近。由圖４（ｃ）（ｅ）可知，隨著含沙量增大，泥沙濃度增加，邊壁與渠底的泥沙濃度大于水中的泥沙濃度。由圖４（ｄ）（ｆ）可知，隨著含沙量增大，渠底部壓強(qiáng)降低，入口處渠底壓強(qiáng)低于渠道中后部的壓強(qiáng)，渠道中后部壓強(qiáng)分布較為均勻。以標(biāo)準(zhǔn)表流量與巴歇爾槽液位為基準(zhǔn)，將近似流量為７１０ ｍ３ ／ ｈ 的模擬流量、液位等數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比，具體見(jiàn)表３。同時(shí)，對(duì)近似流量為３７０、７１０、１ ０１０ ｍ３ ／ ｈ 三種工況計(jì)算不同量水設(shè)備實(shí)測(cè)流量誤差與液位誤差，具體見(jiàn)圖５。由圖５（ａ）（ｃ）（ｅ）可知，在入口近似流量一定情況下，隨著含沙量增大，超聲流量計(jì)流量測(cè)量誤差為－４％～３％，量測(cè)精度較高；電磁流量計(jì)流量測(cè)量誤差在－６％ ～２％之間波動(dòng)，量測(cè)效果較好；雷達(dá)流量計(jì)誤差則較大，誤差為－２０％～４％，且隨著入口流量的增大呈逐漸遞增趨勢(shì)；巴歇爾槽流量計(jì)總體誤差最小，誤差為－２％～０％，量測(cè)效果較好。

由圖５（ｂ）（ｄ）（ｆ）可知，在入口流量一定的情況下，隨著含沙量增大，超聲液位、電磁液位與雷達(dá)液位誤差均較大，電磁液位誤差為－８％ ～０％、超聲液位誤差為－１２％～３％、雷達(dá)液位誤差為－１４％～－３％。

６ 結(jié)論

１）數(shù)值模擬結(jié)果表明，采用ＶＯＦ－ＤＰＭ 模型可以較好地模擬含沙水流在明渠的流態(tài)。

２）通過(guò)開(kāi)展不同流量級(jí)的流量測(cè)量試驗(yàn)可知，在不同含沙量條件下，巴歇爾槽流量計(jì)量總體誤差最小、量測(cè)效果較好、體現(xiàn)了傳統(tǒng)量水堰槽的可靠性，超聲流量計(jì)和電磁流量計(jì)量測(cè)效果均較好，雷達(dá)流量計(jì)誤差較大，說(shuō)明非接觸式流量計(jì)在高含沙量、大流量的情況下需要進(jìn)一步提高測(cè)量精度；超聲液位、電磁液位與雷達(dá)液位誤差均較大，試驗(yàn)中尾部巴歇爾槽壅水以及泥沙在渠道底部淤積是造成較大液位誤差的原因之一，高含沙條件下先進(jìn)行的試驗(yàn)在最大流量情況下將泥沙充分混合，切換到小流量時(shí)導(dǎo)致部分泥沙沉積。

３）本研究可為數(shù)字孿生灌區(qū)建設(shè)、灌區(qū)明渠量水及分離式測(cè)控一體化閘門(mén)量水設(shè)備的合理選取，提供一定的理論依據(jù)和技術(shù)支持。

參考文獻(xiàn)：

［１］ 邵帥超，章少輝，張凱，等．灌區(qū)量測(cè)水方法與技術(shù)研究進(jìn)展［Ｊ］．節(jié)水灌溉，２０２３（９）：１－７．

［２］ 張義強(qiáng)，劉惠忠，付國(guó)義．灌區(qū)量水實(shí)用手冊(cè)［Ｍ］．北京：中國(guó)水利水電出版社，２０１６：４－５．

［３］ 范恬，周義仁．橫向擺桿式明渠測(cè)流方法在矩形渠道中的應(yīng)用研究［Ｊ］．節(jié)水灌溉，２０２２（１）：６２－６６，７３．

［４］ 鮑子云，王潔，李王成，等．明渠非滿管電磁流量測(cè)箱測(cè)試評(píng)價(jià)［Ｊ］．排灌機(jī)械工程學(xué)報(bào)，２０２２，４０（５）：４８２－４８７．

［５］ 萬(wàn)里，李永業(yè)，張港，等．彈簧板式明渠測(cè)流裝置的水力特性［Ｊ］．排灌機(jī)械工程學(xué)報(bào)，２０２２，４０（１２）：１２４１－１２４７．

［６］ 王寶賀，王竹青，蘇沛蘭，等．基于水面一點(diǎn)法的明渠流量計(jì)算方法研究［Ｊ］．節(jié)水灌溉，２０２３（５）：１０９－１１３，１２１．

［７］ ＺＨＯＵ Ｆ，ＹＡＮＧ Ｑ Ｆ，ＬＩＮ Ｋ． Ｅｖａｌｕａｔｉｏｎ Ｉｎｄｅｘ ａｎｄ Ｐｅｒｆｏｒｍ?ａｎｃｅ Ｓｔｒｕｃｔｕｒｅ Ｏｐｔｉｍｉｚａｔｉｏｎ ｏｆ Ｍａｇｎｅｔｉｃ Ｆｉｅｌｄ Ｕｎｉｆｏｒｍｉｔｙ ｏｆ Ｃｏｍｐｌｅｘ Ｍｕｌｔｉｐｈａｓｅ Ｆｌｏｗ Ｅｌｅｃｔｒｏｍａｇｎｅｔｉｃ Ｃｕｒｒｅｎｔ Ｍｅｔｅｒ ［Ｊ］．Ｍｅａｓｕｒｅｍｅｎｔ ａｎｄ Ｃｏｎｔｒｏｌ，２０２２，５５（１－２）：６２－７１．

［８］ 王寶賀，蘇沛蘭，吳建華，等．基于ＧＡ－ＢＰ 神經(jīng)網(wǎng)絡(luò)的明渠流速測(cè)點(diǎn)優(yōu)化研究［Ｊ］．人民黃河，２０２３，４５（１２）：１１７－１２３．

［９］ ＳＨＩ Ｌ Ｊ，ＱＵ Ｍ Ｊ，Ｌｙｕ Ｄ Ｚ，ｅｔ ａｌ．Ａ Ｔｗｏ?Ｃｈａｎｎｅｌ Ｕｌｔｒａｓｏｎｉｃ Ｆｌｏｗｍｅｔｅｒ Ｂａｓｅｄ ｏｎ ＡｌＮ Ｐｉｅｚｏｅｌｅｃｔｒｉｃ Ｍｉｃｒｏｍａｃｈｉｎｅｄ Ｕｌ?ｔｒａｓｏｎｉｃ Ｔｒａｎｓｄｕｃｅｒｓ Ａｒｒａｙｓ ｗｉｔｈ Ｉｍｐｒｏｖｅｄ Ｃｒｏｓｓ?Ｃｏｒｒｅｌａｔｉｏｎ Ｍｅｔｈｏｄ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｍｉｃｒｏｍｅｃｈａｎｉｃｓ ａｎｄ Ｍｉｃｒｏｅｎｇｉｎｅｅｒｉｎｇ，２０２３，３３（１２）：７－８．

［１０］ 段炎沖，楊郁挺，王忠靜，等．含沙量對(duì)明渠超聲時(shí)差法測(cè)流精度的影響［Ｊ／ ＯＬ］．水力發(fā)電學(xué)報(bào)，１－１２［２０２４－０５－１１］． ｈｔｔｐ：／／ ｋｎｓ． ｃｎｋｉ． ｎｅｔ／ ｋｃｍｓ／ ｄｅｔａｉｌ／１１． ２２４１． ＴＶ．２０２２１２０８．１８１１．００２．ｈｔｍｌ．

［１１］ 王福軍．計(jì)算流體動(dòng)力學(xué)分析［Ｍ］．北京：清華大學(xué)出版社，２０１４：２１－２２．

［１２］ 陸力，徐洪泉，孟曉超，等．渾水條件下水頭及流量測(cè)量方法探討［Ｊ］．水電站機(jī)電技術(shù)，２０１７，４０（５）：１－６，１００．

【責(zé)任編輯 簡(jiǎn) 群】

基金項(xiàng)目：山西省基礎(chǔ)研究計(jì)劃（自由探索類）項(xiàng)目（２０２１０３０２１２３１６９，２０２２０３０２１２１２２７１）；山西省高等學(xué)?？萍紕?chuàng)新項(xiàng)目（２０２１Ｌ０２０，２０２２Ｌ０３３）；太原理工大學(xué)校級(jí)科學(xué)研究基金項(xiàng)目（２０２２ＱＮ０５５）；２０２３ 年度山西省研究生教育教學(xué)管理改革項(xiàng)目（２０２３ＪＧ０３１）