植被莖稈傾伏角變化對(duì)明渠水流結(jié)構(gòu)的影響

劉明，張升堂，張景洲，趙慶星，孫榿汶，王磊

植被莖稈傾伏角變化對(duì)明渠水流結(jié)構(gòu)的影響

劉明，張升堂*，張景洲，趙慶星，孫榿汶，王磊

（山東科技大學(xué) 地球科學(xué)與工程學(xué)院，山東 青島 266590）

【】分析植被莖稈傾伏角變化對(duì)明渠水流結(jié)構(gòu)流態(tài)的影響，充分發(fā)揮河道水生植被功能以減緩河堤的侵蝕和水土流失的速率。建立明渠植被水流結(jié)構(gòu)模型將水流液層在空間高度上分為植被層和自由層，在莖稈不同的變形狀態(tài)下及流量范圍內(nèi)開(kāi)展試驗(yàn)。①在同一莖稈傾伏角下，植被層v隨水深增加呈增加減少再增加最終趨于穩(wěn)定的變化規(guī)律，自由層平均速度s、雷諾數(shù)s、弗勞德數(shù)s隨水深的增大迅速減小后趨于穩(wěn)定。②不同莖稈傾伏角下，植被層水力參數(shù)v、v、v在空間上表現(xiàn)出各向異性，自由層s、s、s隨植被傾伏角增大而減小，植被淹沒(méi)臨界水深存在“滯后”現(xiàn)象。自由層的水流流速s、雷諾數(shù)s、弗勞德數(shù)s均大于植被層；植被層水流流態(tài)為緩紊流，自由層水流流態(tài)從急紊流迅速過(guò)渡到緩紊流。

傾伏角；植被層；自由層；流態(tài)

0 引言

【研究意義】全球洪澇災(zāi)害頻發(fā)、水土流失加劇環(huán)境污染日趨嚴(yán)重威脅到人類(lèi)生存，維護(hù)生態(tài)系統(tǒng)平衡至關(guān)重要[1-2]。水生植被作為河流生態(tài)系統(tǒng)的不可替代的一部分，一方面，水生植被不僅可為水生生物提供營(yíng)養(yǎng)物質(zhì)和生長(zhǎng)環(huán)境，而且植被的根莖具有吸附固定作用抑制泥沙懸浮和改善水質(zhì)；另一方面，植被可增加河道水流阻力、流速減緩以及水位升高，從而間接地控制河道泥沙輸移、堤岸侵蝕、航運(yùn)效率、水質(zhì)凈化及水環(huán)境等方面，所以植被是影響河道水流特性的關(guān)鍵因子[3-4]。因此研究河道植被水流流動(dòng)變化規(guī)律在水土保持、穩(wěn)定河堤及環(huán)境工程等領(lǐng)域具有重要意義?！狙芯窟M(jìn)展】針對(duì)植被莖稈對(duì)明渠水流水力特性的研究，國(guó)內(nèi)外研究者提出了明渠植被水流結(jié)構(gòu)“分層理論”方法，在明渠植被水流中沿水深分為若干層以研究植被對(duì)水流水動(dòng)力學(xué)特性的影響。Huai等[5]基于普朗特混合長(zhǎng)度理論進(jìn)行分析，通過(guò)三層水流結(jié)構(gòu)模型預(yù)測(cè)明渠植被水流的垂向時(shí)均流速分布和紊流特性。王雯等[6]研究了剛性雙層淹沒(méi)植被水流特性，依據(jù)垂向流速的拐點(diǎn)劃分為4個(gè)區(qū)域，在植被層的頂端雷諾應(yīng)力和紊動(dòng)能分布明顯增強(qiáng)。Yang等[7]將雙層水流結(jié)構(gòu)模型運(yùn)用于研究淹沒(méi)植被明渠水流，假設(shè)植被層流速分布均勻以及自由層流速遵循對(duì)數(shù)分布，結(jié)果證實(shí)了實(shí)測(cè)數(shù)據(jù)和理論值高度吻合。Nepf等[8]分析淹沒(méi)植被冠層頂部和內(nèi)部水流紊動(dòng)特性，冠層頂部產(chǎn)生的KH渦流控制垂直方向上動(dòng)量交換，冠層內(nèi)部的湍流動(dòng)能與莖稈尺度相關(guān)。Li等[9]研究淹沒(méi)植被的不同排列方式、間距以及密度對(duì)邊界底層、植被層、非植被層水流結(jié)構(gòu)的影響，結(jié)果表明時(shí)均流速在植被的頂部發(fā)生顯著的變化。Yang等[10]采用塑料圓柱模擬河道剛性植被，并建立二維格點(diǎn)玻爾茲曼模型用來(lái)反應(yīng)植被覆蓋下河道水流流場(chǎng)分布。Liu等[11]通過(guò)多普勒測(cè)速儀測(cè)得高矮剛性植被組合水流流速在斷面的分布特征。Tan等[12]通過(guò)水槽試驗(yàn)研究了剛性植被群對(duì)水流流速分布規(guī)律、紊動(dòng)特性和消波效果的影響?！厩腥朦c(diǎn)】以往研究大多集中在探究剛性植被對(duì)水流斷面點(diǎn)流速、雷諾應(yīng)力等水力特性的影響，忽略了生長(zhǎng)在河道中水生植被形態(tài)、種類(lèi)、屬性以及生長(zhǎng)周期等特征并不相同，水生植被在河道中受到水流的沖擊呈現(xiàn)不同程度的傾伏變形?！緮M解決的關(guān)鍵問(wèn)題】比較分析不同植被莖稈傾伏角下植被層、自由層水流水力特性的變化，研究不同植被莖稈傾伏角對(duì)明渠水流特性的影響。為河道輸配水和水土保持提供理論指導(dǎo)。

1 試驗(yàn)設(shè)計(jì)與原理

試驗(yàn)在室內(nèi)長(zhǎng)5 m、寬0.4 m、高0.3 m的矩形明渠水槽中進(jìn)行（圖1），水槽從上游到下游依次分為上游平水段（1 m）、試驗(yàn)鋪設(shè)段（3 m）、下游量水段（1 m），在水槽試驗(yàn)鋪設(shè)段內(nèi)設(shè)置2個(gè)斷面，2個(gè)斷面間隔1.5 m，且每個(gè)斷面均設(shè)有測(cè)壓玻璃管以觀測(cè)斷面水位的變化；試驗(yàn)鋪設(shè)段底部放置有機(jī)玻璃板用來(lái)模擬河道水流的下墊面，有機(jī)玻璃板表面呈現(xiàn)相同孔徑的鉆孔，相鄰鉆孔縱橫間距為60 mm×60 mm，在有機(jī)玻璃板的鉆孔中放置鋁棒以替代水生植被，圓柱鋁棒長(zhǎng)為10 mm、直徑為4 mm，其與豎直面形成不同的夾角（圖2）。水泵將自來(lái)水抽送至蜂窩處，水流經(jīng)過(guò)電磁流量計(jì)和閥門(mén)，進(jìn)入試驗(yàn)鋪設(shè)段。流量和水深分別通過(guò)電磁流量計(jì)、測(cè)壓玻璃管測(cè)量，電磁流量計(jì)的量程為0～0.018 m3/s，測(cè)壓管記錄斷面1、斷面2水深數(shù)據(jù)值，分別記作1、2，將斷面1、斷面2水深取平均值得到平均水深，利用水力學(xué)理論公式計(jì)算平均流速弗勞德數(shù)雷諾數(shù)等相關(guān)水力要素。為研究植被莖稈傾伏角變化對(duì)水流結(jié)構(gòu)水力特性的影響，試驗(yàn)中采用控制變量法，即底坡與莖稈直徑等，其他因素保持不變，僅通過(guò)改變莖稈的傾伏角進(jìn)行量化分析。試驗(yàn)設(shè)置了4組莖稈傾伏角，分別為0°、20°、40°、60°，在坡度1%下開(kāi)展試驗(yàn)（表1）。試驗(yàn)過(guò)程中記錄流量水深以及試驗(yàn)現(xiàn)象。

圖1 試驗(yàn)裝置

圖2 植被變形下雙層水流結(jié)構(gòu)

表1 試驗(yàn)設(shè)置

在河道中阻礙水流流動(dòng)因素眾多，與其他因素相比，植被是制約水流流動(dòng)的主導(dǎo)因素，并且由于植被自身內(nèi)在的特征，如植被形態(tài)、屬性和組合類(lèi)型等特征在空間維度上對(duì)河道植被水流產(chǎn)生了重要影響，使得植被水力參數(shù)的計(jì)算變得復(fù)雜，為了簡(jiǎn)化相關(guān)水力參數(shù)的計(jì)算，將植被視作是均勻分布于河道下墊面的粗糙單元。大多數(shù)學(xué)者采用傳統(tǒng)的一維達(dá)西-魏斯巴赫阻力公式分析，然而植被群作為一個(gè)大尺度的阻水體，在此情況下用達(dá)西-魏斯巴赫阻力系數(shù)及曼寧糙率系數(shù)量化含水植被阻水效果存在一些缺陷[13-14]?？紤]到植被水流相關(guān)水力參數(shù)是由植被的莖稈與水流接觸引起的，一些學(xué)者定義了新的相關(guān)水力參數(shù)以研究植被對(duì)水流流動(dòng)的影響[15-17]。

式中：為淹沒(méi)度；表征植被垂直投影于下墊面的高度與水深之比；v為植被垂直投影于下墊面的高度；w為水深；v為植被相關(guān)的水力半徑（m）；為單位面積上模擬植被的密度；為植被的直徑（mm）；π取值為3.14。

式中：D為植被阻力（kN）；為下墊面單位面積上模擬植被數(shù)量（/m2）；d為植被拖曳力系數(shù)；為水的密度（kg/m3）；v為植被的高度（m）；h為整體水深（m）；v為明渠水流的平均速度（m/s）；g為（9.8 m/s2）；0為能坡（J/m），表征斷面沿程能量損失的比降。

關(guān)于上述式中植被拖曳力系數(shù)的研究，通過(guò)查閱文獻(xiàn)得到v～d的函數(shù)關(guān)系。Ergun[18]和Cheng[19]通過(guò)試驗(yàn)數(shù)據(jù)得到函數(shù)曲線。Cheng[19]分別進(jìn)行了5種不同植被體積分?jǐn)?shù)水流試驗(yàn)，結(jié)果表明同一植被體積分?jǐn)?shù)水流條件下，d隨v的增大呈先遞減變化后趨于穩(wěn)定。Zhang等[20]試驗(yàn)研究不同植被走向和坡度對(duì)拖曳力系數(shù)的影響，隨著水深的增加拖曳力迅速減小后趨于穩(wěn)定。Tanino等[21]研究了不同植被密度對(duì)拖曳力的影響，結(jié)果試驗(yàn)證明v～d關(guān)系與上述等研究結(jié)果相似。在非淹沒(méi)狀態(tài)下，d隨的增加迅速減小后趨于穩(wěn)定，當(dāng)v大于某一值時(shí)，d變化幅度較小趨于穩(wěn)定。Schlichting等[22]給出了不同雷諾數(shù)值下的植被水流拖曳力d值：

弗勞德數(shù)是用來(lái)判斷水流流型的重要參數(shù)，表征動(dòng)能與勢(shì)能之比，植被層v定義為：

式中：v為植被層的弗勞德數(shù)；v為植被層的平均水流速度（m/s）；g為重力加速度（9.8 m/s2）；v為植被層水深（m/s）。

上述已經(jīng)給出植被層水流相關(guān)水力參數(shù)的推導(dǎo)公式，依據(jù)過(guò)水?dāng)嗝婵偟牧髁康扔谥脖粚印⒆杂蓪恿髁恐?，可以解得到自由層其他相關(guān)水力參數(shù)s、s。

式中：e為植被層有效過(guò)水?dāng)嗝婷娣e（m2）；v為植被層平均水流速度（m/s）；s為自由層過(guò)水?dāng)嗝婷娣e（m2）；s為自由層水流平均速度（m/s）；下標(biāo)“s”表征自由層相關(guān)的水力參數(shù)。

2 結(jié)果與分析

2.1 傾伏角θ及水深h對(duì)植被層平均流速Vv的影響

圖3、圖4、圖5分別為植被變形作用下植被層水流平均速度v、雷諾數(shù)v、弗勞德數(shù)v、與水深的關(guān)系。由圖3可知，當(dāng)植被傾伏角保持不變時(shí)，植被層的水流v與水深呈倒“S”形分布[23]。水深在0～0.02 m范圍內(nèi)，水流v隨著水深的增加呈顯著上升的趨勢(shì)，水流v切線斜率逐漸增大，此時(shí)水深是影響植被層水流v的主要因素；水深0.02～0.04 m時(shí)，隨著明渠植被水深逐漸增加，同時(shí)植被與水流接觸的橫截面積增加，水流阻力增大，植被的“阻塞效應(yīng)”大于水深增加對(duì)水流v貢獻(xiàn)值，故植被層水流v與水深呈現(xiàn)負(fù)相關(guān)關(guān)系。水深大于0.04 m時(shí)，植被層水流v隨水深的增加而增加后趨于穩(wěn)定。可見(jiàn)，植被的存在使明渠水流產(chǎn)生植被繞流現(xiàn)象，植被的迎水面阻擋水流流動(dòng)，水流液面雍高，液體層間的黏滯力使水體在植被變形周?chē)蛛x，植被中間的雍水向二側(cè)擠壓，水流向植被之間的空隙流動(dòng)，促使植被層水流v趨向均勻化[24]。

圖3可以看出，由0°增加到60°時(shí)，植被層水流v表現(xiàn)為60°>40°>20°>0°，當(dāng)植被傾伏角為60°時(shí)，植被與水流流向之間的夾角較小，順于水流流向，對(duì)水流流速影響較弱；當(dāng)傾伏角為0°時(shí)，植被相當(dāng)于垂直下墊面的“小閘門(mén)”串聯(lián)成一排，水流流經(jīng)“小閘門(mén)”時(shí)，與“小閘門(mén)”發(fā)生強(qiáng)烈的正面碰撞，會(huì)產(chǎn)生較大的機(jī)械能損失，因此植被層水流流速受傾伏角的影響，可提高河道行洪能力和削弱堤岸侵蝕速率。

圖3 植被層不同傾伏角下h～Vv關(guān)系曲線

2.2 傾伏角θ及水深h對(duì)植被層水流流態(tài)的影響

由圖4可知，在～v關(guān)系曲線中，植被層雷諾數(shù)v變化區(qū)間（1 400～65 000），水流流態(tài)均為紊流。當(dāng)試驗(yàn)植被恒定時(shí)，植被層v隨著水深增加呈現(xiàn)線性增加的趨勢(shì)。由于水深的增加，水體與植被粗糙單元接觸的面積加大，植被的存在增加了對(duì)水流的擾動(dòng)性。同一水深試驗(yàn)下，當(dāng)水深大于0.04 m時(shí)，不同傾伏角的形態(tài)植被在水中，植被層水流v隨增加而增加。因此，傾伏角是植被層水流v影響因子之一；當(dāng)水深小于0.04 m時(shí)，不同傾伏角度下植被層水流雷諾數(shù)差異較小，液體層間的黏滯力及水槽邊壁是影響雷諾數(shù)的主導(dǎo)因素。

圖4 植被層不同傾伏角下h～Rev關(guān)系曲線

圖5為植被層弗勞德數(shù)v與水深之間變化曲線。在不同植被傾伏角試驗(yàn)下，植被層水流弗勞德數(shù)v變化區(qū)間為（0.29～0.92），水流流型均為緩流。同一傾伏角下，隨水深的增大呈先減小后趨于穩(wěn)定，隨著植被層水深增加，水流v大體呈增加的趨勢(shì)，但是速度的增長(zhǎng)速率逐漸減小，而水深增加的效果凸顯出優(yōu)勢(shì)，水流的重力勢(shì)能增加，故植被層v隨水深增加呈顯著減小的趨勢(shì)。當(dāng)明渠水流整體水深繼續(xù)增加，淹沒(méi)狀態(tài)下植被層水深保持不變，水流流速趨于穩(wěn)定，植被層水流v隨水深增加呈不變趨勢(shì)；同一水深不同植被傾伏角試驗(yàn)下，植被層v隨增大而增大。

圖5 植被層不同傾伏角下h～Frv關(guān)系曲線

2.3 傾伏角θ及水深h對(duì)自由層平均流速Vs的影響

圖6、圖7分別為自由層流速s與水深關(guān)系曲線圖和流態(tài)流型分區(qū)。從圖6可以看出，相同傾伏角條件下，如在60°植被傾伏角～s曲線中，隨著水深的增加，自由層的平均水流流速先急劇下降后趨于穩(wěn)定。植被初始淹沒(méi)狀態(tài)下，植被層和自由層的液面之間發(fā)生密切的水力，水力要素急劇變化，水流紊動(dòng)增強(qiáng)，自由層水流s迅速減小。同一水深不同植被試驗(yàn)下，植被傾伏角越大，其對(duì)應(yīng)植被淹沒(méi)臨界水深值越小，不同傾伏角下植被淹沒(méi)臨界水深值出現(xiàn)“滯后”現(xiàn)象，自由層平均流速s隨傾伏角的增大逐漸減小。

2.4 傾伏角θ及水深h對(duì)自由層流態(tài)的影響

圖7為不同傾伏角下自由層流態(tài)分區(qū)。自由層水流流態(tài)劃分為6個(gè)區(qū)：緩層流、緩過(guò)渡流、緩紊流、急層流、急過(guò)渡流、急紊流。不同植被下自由層水流流速數(shù)據(jù)點(diǎn)大部分落在緩紊流區(qū)，僅在淹沒(méi)臨界水深附近處，少數(shù)試驗(yàn)點(diǎn)分布于急紊區(qū)，其他4個(gè)區(qū)均未分布。受下層植被頂部擾動(dòng)的影響，水流紊動(dòng)加劇。不同植被傾伏角下自由層值s及s存在差異，同一傾伏角工況下，s及s隨水深增加迅速減小后趨于穩(wěn)定。由于自由層水流慣性力減弱且重力勢(shì)能增加，水流流態(tài)從急紊流迅速過(guò)渡到緩紊流，削弱了對(duì)河岸的侵蝕力度。

圖6 自由層不同傾伏角下的h～Vs關(guān)系曲線

圖7 自由層不同傾伏角下流態(tài)分區(qū)圖

3 討論

張景洲等[25]研究了植被的不同排列方式對(duì)坡面糙率及水流特性的影響，采用塑料棒模擬陸地植被與水流呈不同走向角，且植被垂直于水槽底板，結(jié)果表明相同水深工況下，植被走向角越小，植被莖稈與水流的接觸面積越大，相應(yīng)的水流阻力系數(shù)越大，流速越小。這與本試驗(yàn)不一致，考慮植被莖稈傾伏角的變化對(duì)明渠水流特性的影響，相同水深情況下，莖稈傾伏角越小，其莖稈的固液接觸面積愈小，水流阻力愈小，水流流速越小。究其原因?yàn)椋褐脖凰髁魉偈苤饕蛩赜绊懠此髯枇颓o稈與水流碰撞動(dòng)能的損失，莖稈傾伏角越小，水流阻力越小，傾伏角小的莖稈與水流發(fā)生正面碰撞損失的動(dòng)能較多，故本試驗(yàn)植被層的平均流速隨著莖稈傾伏角增大而增大。自由層的水流受莖稈傾伏角的影響，淹沒(méi)臨界水深值出現(xiàn)“滯后”現(xiàn)象，在淹沒(méi)臨界水深處，自由層的平均流速受植被頂端的擾動(dòng)，水流流速迅速減小。隨著水深的增加，不同傾伏角下自由層的流速趨于某一值。

本試驗(yàn)未從微觀斷面點(diǎn)流速進(jìn)行探討，在莖稈傾伏角的變化下，斷面點(diǎn)流速以及水質(zhì)點(diǎn)劇烈摻混引起的附加應(yīng)力如何變化有待進(jìn)一步開(kāi)展試驗(yàn)研究。

本試驗(yàn)在室內(nèi)開(kāi)展，受某些因素的制約，比如試驗(yàn)水槽的寬度、坡度以及坡長(zhǎng)等因素。即使試驗(yàn)是在等比例尺水力模型下進(jìn)行，但仍和實(shí)際自然情況有差距。實(shí)際自然條件的坡度較陡，植被段的寬度較大。因此試驗(yàn)有待繼續(xù)進(jìn)一步優(yōu)化探索。

4 結(jié)論

采用“分層理論”方法，將水流液層在空間高度上分為植被層和自由層，植被莖稈傾伏角變化對(duì)植被層、自由層水力特性影響如下：

①植被層：不同莖稈傾伏角下，植被層v、v、v與傾伏角均正相關(guān)，傾伏角越大，動(dòng)能損失越小，對(duì)應(yīng)的植被層流速越大，水流流態(tài)均為緩紊流，植被層水流流速受傾伏角的影響較大。

②自由層：不同莖稈傾伏角下，自由層s、s、s隨傾伏角的增大而減少，水流流態(tài)從急紊流迅速過(guò)渡到緩紊流，水深淹沒(méi)臨界值出現(xiàn)“滯后”現(xiàn)象；隨著水深進(jìn)一步增加，莖稈的傾伏角變化對(duì)自由層水流影響減弱。

[1] DEWAN T H. Societal impacts and vulnerability to floods in Bangladesh and Nepal[J]. Weather and Climate Extremes, 2015, 7: 36-42.

[2] 史志華, 楊潔, 李忠武, 等. 南方紅壤低山丘陵區(qū)水土流失綜合治理[J]. 水土保持學(xué)報(bào), 2018, 32(1): 6-9.

SHI Zhihua, YANG Jie, LI Zhongwu, et al. Soil conservation in the hilly red soil region of Southern China[J]. Journal of Soil and Water Conservation, 2018, 32(1): 6-9.

[3] 王子建, 丁雪, 吉慶豐. 含剛性沉水植物明渠水流結(jié)構(gòu)的試驗(yàn)研究[J]. 灌溉排水學(xué)報(bào), 2020, 39(2): 78-83.

WANG Zijian, DING Xue, JI Qingfeng. Experimental study on flow structure of open channels with rigid submerged plants[J]. Journal of Irrigation and Drainage, 2020, 39(2): 78-83.

[4] FRINGS R M, D?RING R, BECKHAUSEN C, et al. Fluvial sediment budget of a modern, restrained river: The lower reach of the Rhine in Germany[J]. Catena, 2014, 122: 91-102.

[5] HUAI W X, CHEN Z B, HAN J, et al. Mathematical model for the flow with submerged and emerged rigid vegetation[J]. Journal of Hydrodynamics, Ser B, 2009, 21(5): 722-729.

[6] 王雯, 槐文信. 剛性淹沒(méi)雙層植被明渠水流紊流特性研究[J]. 四川大學(xué)學(xué)報(bào)(工程科學(xué)版), 2014, 46(1): 61-67.

WANG Wen, HUAI Wenxin. Turbulence structure of open channel flow through double layer rigid submerged vegetation[J]. Journal of Sichuan University (Engineering Science Edition), 2014, 46(1): 61-67.

[7] YANG W, CHOI S U. A two-layer approach for depth-limited open-channel flows with submerged vegetation[J]. Journal of Hydraulic Research, 2010, 48(4): 466-475.

[8] NEPF H, GHISALBERTI M. Flow and transport in channels with submerged vegetation[J]. Acta Geophysica, 2008, 56(3): 753-777.

[9] CHEN S C, KUO Y M, LI Y H. Flow characteristics within different configurations of submerged flexible vegetation[J]. Journal of Hydrology, 2011, 398(1/2): 124-134.

[10] YANG Z H, BAI F P, HUAI W X, et al. Modelling open-channel flow with rigid vegetation based on two-dimensional shallow water equations using the lattice Boltzmann method[J]. Ecological Engineering, 2017, 106: 75-81.

[11] LIU D, DIPLAS P, HODGES C C, et al. Hydrodynamics of flow through double layer rigid vegetation[J]. Geomorphology, 2010, 116(3/4): 286-296.

[12] TAN C, HUANG B S, LIU D, et al. Effect of mimic vegetation with different stiffness on regular wave propagation and turbulence[J]. Water, 2019, 11(1): 109.

[13] FATHI-MOGHADAM M, KASHEFIPOUR M, EBRAHIMI N, et al. Physical and numerical modeling of submerged vegetation roughness in rivers and flood Plains[J]. Journal of Hydrologic Engineering, 2011, 16(11): 858-864.

[14] BANERJEE S, NAIK B, SINGH P, et al. Flow resistance in gravel bed open channel flows case: Intense transport condition[J]. ISH Journal of Hydraulic Engineering, 2019, 25(3): 298-309.

[15] CHENG N S, NGUYEN H T. Hydraulic radius for evaluating resistance induced by simulated emergent vegetation in open-channel flows[J]. Journal of Hydraulic Engineering, 2011, 137(9): 995-1 004.

[16] D’IPPOLITO A, LAURIA A, ALFONSI G, et al. Investigation of flow resistance exerted by rigid emergent vegetation in open channel[J]. Acta Geophysica, 2019, 67(3): 971-986.

[17] KOTHYARI U C, HAYASHI K, HASHIMOTO H. Drag coefficient of unsubmerged rigid vegetation stems in open channel flows[J]. Journal of Hydraulic Research, 2009, 47(6): 691-699.

[18] ERGUN S. Fluid flow through packed columns[J]. Journal of Materials Science and Chemical Engineering, 1952, 48(2): 89-94.

[19] CHENG N S. Calculation of drag coefficient for arrays of emergent circular cylinders with pseudofluid model[J]. Journal of Hydraulic Engineering, 2013, 139(6): 602-611.

[20] ZHANG S T, LIU Y C, ZHANG J Z, et al. Study of the impact of vegetation direction and slope on drag coefficient[J]. Iranian Journal of Science and Technology, Transactions of Civil Engineering, 2018, 42(4): 381-390.

[21] TANINO Y, NEPF H M. Closure to “l(fā)aboratory investigation of mean drag in a random array of rigid, emergent cylinders” by yukie tanino and heidi M. nepf[J]. Journal of Hydraulic Engineering, 2009, 135(8): 693-694.

[22] SCHLICHTING B H, KESTIN T B J. Boundary-layer theory[J]. Fluid Mechanics, 1990, 20(1): 48-89.

[23] STOESSER T, KIM S J, DIPLAS P. Turbulent flow through idealized emergent vegetation[J]. Journal of Hydraulic Engineering, 2010, 136(12): 1 003-1 017.

[24] MAYAUD J R, WIGGS G F S, BAILEY R M. Characterizing turbulent wind flow around dryland vegetation[J]. Earth Surface Processes and Landforms, 2016, 41(10): 1 421-1 436.

[25] 張景洲, 張升堂, 劉元臣, 等. 植被排列方式對(duì)坡面糙率系數(shù)的影響研究[J]. 人民黃河, 2017, 39(5): 93-96, 111.

ZHANG Jingzhou, ZHANG Shengtang, LIU Yuanchen, et al. Studyon impact of vegetation arrangement on slope surface roughness coefficient[J]. Yellow River, 2017, 39(5): 93-96, 111.

The Effects of Plant Stem Inclination on Water Flow in Open Channel

LIU Ming, ZHANG Shengtang*, ZHANG Jingzhou, ZHAO Qingxing, SUN Qiwen, WANG Lei

( College of Earth Science and Engineering, Shandong University of Science and Technology, Qingdao 266590, China)

【】Most open channels are planted with vegetation to avoid erosion of channel beds and slopes while in the meantime facilitating decontamination of organic and inorganic pollutants in the water. However, the vegetation could detrimentally affect water flow. How the depth and growing direction of the vegetation affect water flow across the channel is an issue not well understood. The objective of this paper is to plug this gap.【】A physical model was established to simulate water flow in open channel with the vegetation taken into account. Water in the channel was divided into two layers: A vegetation layer and an overlaid free layer. The impact of the growing direction of the vegetation, characterized by their stem inclination angle, was studied experimentally.【】①When the inclination angle of the stems was the same, the water velocity in the vegetation layer,v, tended to level off as the water depth in the channel,, increased. In the free layer, the average velocity,s, the Reynolds numbers, Froude numbersincreased asymptotically asincreased. ②We found thatv,Reandvin the vegetation layer were spatially anisotropic varying with the stem inclination angle, while in the free layers,sandsdecreased as the stem inclination angle increased. There was a delay in the critical water depth above which the vegetation was submerged. 【】Water velocitys,sandsin the free layer were higher than that in the vegetation layer. Water flow in the vegetation layer was slow and turbulent, while in the free layer the water flow was rapid, in a transition from fast turbulence to slow turbulence. Our results can help open channel design on how to reduce erosion of its bed and slope by growing vegetation without risking flooding.

pitching angle; vegetation layer; free water flowing layer; water flow regime

TV14

A

10.13522/j.cnki.ggps.2020523

1672 – 3317（2021）08 - 0088 - 06

劉明, 張升堂, 張景洲, 等. 植被莖稈傾伏角變化對(duì)明渠水流結(jié)構(gòu)的影響[J]. 灌溉排水學(xué)報(bào), 2021, 40(8): 88-92, 128.

LIU Ming, ZHANG Shengtang, ZHANG Jingzhou, et al.The Effects of Plant Stem Inclination on Water Flow in Open Channel[J]. Journal of Irrigation and Drainage, 2021, 40(8): 88-92, 128.

2020-09-18

國(guó)家自然科學(xué)基金項(xiàng)目（41471025）；山東科技大學(xué)研究生科技創(chuàng)新項(xiàng)目（SDKDYC180320）

劉明（1994-），男。碩士研究生，主要從事水文水資源方面研究。E-mail: lm0029@126.com

張升堂（1970-），男。教授，主要從事水文水資源方面研究。E-mail: zst0077@163.com

責(zé)任編輯：白芳芳