基于豎縫隔板的大底坡魚道水力優(yōu)化數(shù)值模擬

鄧 斌,黃姣鳳,徐 拓,尹龍斌,張家強(qiáng),蔣昌波

(1.長沙理工大學(xué)水利與環(huán)境工程學(xué)院,湖南 長沙 410114; 2.水沙科學(xué)與水災(zāi)害防治湖南省重點(diǎn)實(shí)驗(yàn)室,湖南 長沙 410114; 3.洞庭湖水環(huán)境治理與生態(tài)修復(fù)湖南省重點(diǎn)實(shí)驗(yàn)室,湖南 長沙 410114)

閘、壩等水利設(shè)施的興建在給防洪、灌溉、緩解缺水、開發(fā)清潔能源等方面帶來巨大效益[1]的同時會截?cái)嗪恿?導(dǎo)致完整的河流碎片化[2],阻斷魚類洄游通道并影響其種群繁殖[3]。魚道作為一種生態(tài)型水利工程,能夠重建河流連通性,幫助魚類克服障礙恢復(fù)自然洄游[4],對保護(hù)魚類資源及水生生物多樣性具有一定的積極作用。魚道主要有仿自然型魚道(仿自然魚梯、魚坡和過魚旁路等)和技術(shù)型魚道(豎縫式魚道、池堰式魚道和丹尼爾魚道等)兩種類型[5],其中豎縫式魚道因結(jié)構(gòu)簡單、能適應(yīng)多變的上下游水位、消能效果顯著等特點(diǎn),成為應(yīng)用最廣泛和最有效的魚道型式之一[6]。

研究發(fā)現(xiàn),我國已建的豎縫式魚道大多底坡較緩,為1/100～1/50[7]。而底坡坡度決定了魚道的總長度,是影響魚類上溯時間、工程建設(shè)成本及規(guī)劃布置的主要因素。一般來說,魚道底坡越陡,總長度越短,占地面積就越小,建設(shè)成本也越低。但是大底坡結(jié)構(gòu)所形成的高速度、高湍動水流流場不利于游泳能力較弱的魚類洄游上溯[8]。因此,考慮河流生態(tài)修復(fù)及工程經(jīng)濟(jì)性,國內(nèi)外學(xué)者針對大底坡下豎縫式魚道的水力優(yōu)化開展了相關(guān)研究,如部分學(xué)者基于物理模型試驗(yàn)[9-10]和數(shù)值模擬[11-14]方法,分別研究了在各種大底坡(1/20～1/7)結(jié)構(gòu)下,豎縫式魚道豎縫下游處設(shè)置不同直徑、個數(shù)及排列方式的圓柱體障礙物對池室水力特性的影響,并通過比較池室流速和池室湍動特性等相關(guān)參數(shù)來評估不同圓柱體配置方式的阻流效果。此外,也有部分學(xué)者研究了在豎縫處增設(shè)隔板對池室水流分布的影響,如Pena等[15]認(rèn)為在底坡為1/10的豎縫式魚道豎縫底部設(shè)置不同長度的隔板可以減小池室最大流速,但同時會產(chǎn)生比較復(fù)雜的流場結(jié)構(gòu)并增加池室湍動能;Mallen-Cooper等[16]基于原型觀測研究了底坡為1/20時在豎縫式魚道豎縫中間位置設(shè)置隔板對魚道過魚效率的影響,發(fā)現(xiàn)與未經(jīng)改造的豎縫式魚道相比,其小魚通過效率提升了6～13倍,這種通過在豎縫式魚道結(jié)構(gòu)水流最集中的部位設(shè)置消能措施以減少池室流速和湍流的方法比較有效、簡單和經(jīng)濟(jì),但是鮮見其相關(guān)的池室水力特性研究報(bào)道,而水力特征分布是探討魚類上溯行為機(jī)理的依據(jù),也亟須對其進(jìn)行相關(guān)研究。

綜上所述,通過在豎縫式魚道結(jié)構(gòu)內(nèi)增設(shè)局部消能措施(如圓柱體、豎縫隔板等)可改變池室流場分布,并且關(guān)于不同配置方式的圓柱體對池室水力特性的影響已進(jìn)行了一定的研究,然而對于豎縫隔板所產(chǎn)生的影響,目前相關(guān)研究尚不夠全面、系統(tǒng),還有待進(jìn)一步探索與分析。因此,本文從水力、工程經(jīng)濟(jì)等角度出發(fā),在Mallen-Cooper等[16]的研究基礎(chǔ)上,通過建立魚道池室三維水動力數(shù)值模型,深入分析了5種大底坡J(1/40、1/30、1/20、1/15、1/10)條件下,豎縫隔板及其不同安裝高度h(0、0.04、0.08m,為豎縫隔板底部到池室底板的距離)對池室水力特性的影響,并闡明了豎縫隔板影響和優(yōu)化大底坡魚道池室流場分布的內(nèi)在原因。通過對池室水力特性的數(shù)值模擬研究,以期獲取關(guān)于這類魚道結(jié)構(gòu)運(yùn)行的基本知識,可為改進(jìn)更有效、經(jīng)濟(jì)的大底坡魚道結(jié)構(gòu)提供相關(guān)技術(shù)依據(jù)。

1 數(shù)學(xué)模型

1.1 控制方程

本研究采用計(jì)算流體力學(xué)商業(yè)軟件FlOW-3D對豎縫式魚道內(nèi)的水動力進(jìn)行模擬,采用有限體積法求解三維雷諾平均納維-斯托克斯(RANS)方程。湍流模型選用重整化群(RNG)k-ε模型,該模型能夠較準(zhǔn)確地模擬復(fù)雜水工結(jié)構(gòu)內(nèi)的流場分布[17]。

采用流體體積(VOF)方法模擬自由液面,即氣相和液相之間的界面。VOF方法通過求解整個計(jì)算域內(nèi)各網(wǎng)格單元的體積分?jǐn)?shù),并進(jìn)行持續(xù)記錄,構(gòu)建出自由液面。其表達(dá)式為

(1)

式中:aw為水的體積分?jǐn)?shù);t為時間;ui為流體速度;xi為坐標(biāo)。當(dāng)aw=0時,計(jì)算域內(nèi)全部充滿空氣;當(dāng)aw=1時,計(jì)算域內(nèi)全部充滿水;當(dāng)0

1.2 計(jì)算域及網(wǎng)格劃分

數(shù)學(xué)模型計(jì)算域包括魚道入流段、出流段及5級池室,如圖1所示。豎縫寬度b為0.05m,單個池室凈長L=10b=0.5m,凈寬B=8b=0.4m,池室長寬比為10/8,底坡J分別取1/40、1/30、1/20、1/15、1/10。豎縫隔板尺寸為0.05m×0.01m×0.05m(長×寬×高),其中高度與豎縫寬度相等,布置方式如圖2所示。為減少進(jìn)出口邊界對池室流場的影響,取第三級池室作為主要研究對象,采用六面體結(jié)構(gòu)化網(wǎng)格對模型進(jìn)行網(wǎng)格劃分,網(wǎng)格單元邊長設(shè)為0.006m。

圖1 魚道數(shù)學(xué)模型平面示意圖

圖2 豎縫隔板布置三維示意圖

1.3 邊界條件及初始條件

考慮到水流流動的物理特性及工程模擬的相似性,魚道上游進(jìn)口及下游出口邊界條件均采用壓力邊界,并設(shè)置相應(yīng)的水深,魚道進(jìn)口水深為0.195m,出口水深為0.245m;模型頂部設(shè)置為壓力進(jìn)口邊界,其相對壓強(qiáng)大小為0,水體體積分?jǐn)?shù)設(shè)置為0,以便空氣可自由進(jìn)入;魚道左右邊墻及底部設(shè)置為固壁無滑移邊界。模型初始水深H設(shè)置為0.2m,以使計(jì)算更快收斂并更加符合實(shí)際情況。模型采用廣義最小殘量(GMRES)隱式求解器計(jì)算,初始時間步長設(shè)置為0.001s,最小時間步長設(shè)置為10-7s。

1.4 數(shù)值模擬與物理模型試驗(yàn)結(jié)果對比分析

為了評估數(shù)學(xué)模型模擬魚道池室水動力特性的能力,通過調(diào)整數(shù)學(xué)模型中魚道邊墻及底部的表面粗糙度rs(0.010、0.012、0.015、0.018),采用與數(shù)學(xué)模型尺寸相一致的物理模型對各表面粗糙度下的流速進(jìn)行了對比驗(yàn)證。物理模型試驗(yàn)在長沙理工大學(xué)水利實(shí)驗(yàn)中心的變坡水槽中進(jìn)行,水槽尺寸為16m×0.4m×0.5m(長×寬×高),水槽底面及兩側(cè)均為透明的鋼化玻璃。為得到充分發(fā)展的水流,將模型布置在水槽中后段部分,并通過多普勒聲學(xué)流速儀測量了第三級池室z=0.5H水深平面處的速度場分布,池室流速測點(diǎn)布置如圖3所示。對比分析了底坡為1/20且無豎縫隔板情況下,距上游長擋板背水面距離0.22m斷面線處的物理模型與數(shù)學(xué)模型流速分布(圖4(a))及不同水深豎縫斷面中心點(diǎn)連線處的流速分布(圖4(b)),并通過決定系數(shù)R2對數(shù)學(xué)模型計(jì)算精度進(jìn)行評價。驗(yàn)證結(jié)果表明,當(dāng)rs=0.015時,數(shù)學(xué)模型結(jié)果與物理模型結(jié)果吻合較好,R2都在0.9以上,數(shù)學(xué)模型能夠較準(zhǔn)確地模擬魚道池室內(nèi)的水力特性。

圖3 第三級池室測點(diǎn)布置(單位:m)

2 結(jié)果與分析

根據(jù)上述數(shù)學(xué)模型,分別對各工況進(jìn)行了模擬計(jì)算,為減少魚道進(jìn)出口邊界條件對計(jì)算結(jié)果的影響,以中間第三級池室作為主要研究對象進(jìn)行水力分析與討論。

2.1 不同底坡池室水力參數(shù)對比

在豎縫式魚道的水力設(shè)計(jì)中,池室最大流速vmax和池室湍動能k是影響過魚效果最主要的水力參數(shù)[18]。由于豎縫式魚道在相鄰池室間會產(chǎn)生水位落差,vmax一般出現(xiàn)在豎縫位置附近,而豎縫是魚類進(jìn)行洄游的必經(jīng)之地,因此,為使魚類能夠順利通行,vmax必須小于魚的爆發(fā)速度[19]。k表示脈動流速的動能,是衡量水流湍動程度的重要指標(biāo)[20];若魚類上溯路徑中的水流k過大,會導(dǎo)致魚類特別是游泳能力較弱的魚類在上溯過程中難以保持方向,同時增加魚類游泳耗能甚至對魚類造成不同程度的傷害[21],這表明池室中適當(dāng)?shù)膋分布也是魚類成功上溯的關(guān)鍵。因此,優(yōu)化大底坡魚道池室水力分布的主要目標(biāo)是減小池室最大流速和池室湍流。

表1 不同底坡最大流速和平均湍動能下降率

2.2 豎縫隔板對池室水力特性的影響

為了進(jìn)一步闡明豎縫隔板對魚道池室流場分布的影響機(jī)理,選取第三級池室距底板0.03m(豎縫隔板下方水深處,除h=0m時)及0.15m(豎縫隔板上方水深處)的水平面作為主要研究對象,具體分析了底坡為1/20時,不同豎縫隔板配置方式在Q=19.36m3/h和Q=10.23m3/h(Q為魚道入口流量)下的池室水力特性分布。

2.2.1流態(tài)

圖5為不同流量條件下豎縫隔板不同配置方式的三維流線分布:當(dāng)豎縫處無隔板時(圖5(a)(e)),水流通過上游豎縫后,以射流的形式進(jìn)入池室,并沿池室隔板導(dǎo)角方向向前流動,射流到達(dá)池室中心后受下游擋板墩頭影響,其斷面逐漸收縮并同時折向下游豎縫進(jìn)入下一級池室,呈現(xiàn)出反S形的主流,并在主流兩側(cè)產(chǎn)生了兩個尺度相當(dāng)?shù)幕亓鲄^(qū);而對于大底坡和流量較大的情況,這種流型的主流流速(0.5～0.6m/s)過大,會導(dǎo)致魚類上溯耗能的增加并不利于游泳能力較弱的魚類通行。豎縫處設(shè)置隔板后(圖6(b)(c)(d)(f)(g)(h)),水流也以射流形式流入池室,但整體流動模式發(fā)生了顯著變化:豎縫隔板的存在改變了射流流向,形成了利用池室整體容積進(jìn)行低湍流能量轉(zhuǎn)換、自身回轉(zhuǎn)的主流,使得主流路徑更加彎曲并呈Ω形,該流型通過射流混合消能的方式,可有效延長大底坡、大流量條件下的水流能量耗散距離,在池室內(nèi)形成面積較大的低流速區(qū)(0.25～0.35m/s),使魚類能以持續(xù)游泳速度和耐久游泳速度相交替的游泳方式[22]游過各級魚道池室,有利于減小魚類上溯耗能;同時,在池室的上游角和主流的凸出部分產(chǎn)生了兩個尺度不一的回流區(qū),給魚類洄游提供休息區(qū)域。而改變豎縫隔板安裝高度,流態(tài)呈現(xiàn)相似的分布規(guī)律。

圖6 不同水深處池室流速分布

圖6為Q=19.36m3/h時豎縫隔板不同安裝情況下不同水深處的池室流速分布。與無隔板相比,3種隔板安裝高度下的池室流速都顯著下降,且具有相似的分布情況;當(dāng)h為0時,在0.03m水深處(圖6(b)),由于豎縫隔板的阻擋,其豎縫下游處的流速分布較緩;而豎縫隔板上下水深處的流速分布(圖6(c)(d)(f)(g)(h))大體保持一致,這表明改變豎縫隔板安裝高度,對其上下水深處的流速分布無顯著影響。

2.2.2主流區(qū)水力特性

圖7 最大流速軌跡線

圖8 最大流速沿程分布曲線

2.2.3豎縫斷面流速

分析發(fā)現(xiàn),不同流量條件下豎縫斷面流速分布規(guī)律具有較好的相似性。因此,圖9僅給出了Q=19.36m3/h時H=0.03m和H=0.15m處豎縫斷面流速沿豎縫寬度b方向的分布規(guī)律。在豎縫處設(shè)置隔板使豎縫中間位置(b′/b=0.3～0.9,其中b′為豎縫上的點(diǎn)到隔板的距離)的流速顯著下降,而對靠近兩側(cè)擋板處的流速影響較小;同時,當(dāng)h=0m時,魚道底部由于豎縫隔板的阻擋,使得魚道底部(H=0.03m)豎縫流速為0,因此圖9(a)未給出其流速分布,當(dāng)h=0.04m時,對豎縫隔板上方水深處(H=0.15m)的流速影響較小,當(dāng)h=0.08m時,豎縫隔板上下水深處(H=0.03m、0.15m)的流速都有明顯降低,這表明豎縫隔板能顯著減小其安裝位置附近水深處的豎縫斷面流速,且h越大,這種影響范圍也越大?；诖?初步認(rèn)為可將隔板安裝在豎縫的中上方水深處,以保障魚類能夠更加順利地通過豎縫完成洄游。

圖9 不同水深處豎縫斷面流速分布

2.2.4湍動能

圖10為Q=19.36m3/h時豎縫隔板不同安裝情況在不同水深處的池室湍動能分布:不同工況池室湍動能的分布與其所對應(yīng)的池室流型分布具有較好的一致性;在不同水深處(除豎縫隔板阻擋水深處),湍動能最大值均出現(xiàn)在豎縫附近,同時水流通過豎縫后由于主流區(qū)流速較大,兩側(cè)回流區(qū)流速較小,造成主流邊界流速梯度較大,因此在主流及其左右邊界的剪切區(qū)域亦表現(xiàn)出較大的湍動能,并隨著水深的增加,池室中湍動能的分布規(guī)律大致相同,但相同位置的湍動能有所增大;與無隔板相比,豎縫處設(shè)置隔板可降低池室湍動能的大小,并使其在豎縫下游位置區(qū)域內(nèi)迅速衰減,產(chǎn)生連續(xù)的低湍動區(qū)域(0.007～0.01J/kg),這給魚類特別是游泳能力較弱的魚類提供了一個較好的上溯條件;同流速分布一樣,改變豎縫隔板安裝高度,對豎縫隔板上下水深處池室湍動能的分布規(guī)律無顯著影響。

圖10 不同水深處池室湍動能分布

2.2.5應(yīng)變率

應(yīng)變率(sr)為水流各方向速度梯度的絕對值之和。相關(guān)研究[24-25]認(rèn)為水流速度梯度的大小是影響魚類側(cè)線系統(tǒng)感知水流刺激并產(chǎn)生運(yùn)動的重要因素,因此,應(yīng)變率也是研究魚道水力特性的重要參數(shù)之一。圖11為Q=19.36m3/h時不同隔板安裝情況下不同水深處的池室應(yīng)變率分布。與池室湍動能分布一樣,應(yīng)變率分布也取決于池室流型的分布,并沿主流兩側(cè)邊界有逐漸減小的趨勢;同時在豎縫位置內(nèi)側(cè)以及靠近豎縫處下游位置的兩側(cè)呈現(xiàn)相對較大的數(shù)值,即在這些區(qū)域水流速度變化顯著、流速梯度大(豎縫隔板阻擋水深處除外);設(shè)置豎縫隔板會使應(yīng)變率的分布在豎縫下游處發(fā)生偏轉(zhuǎn),但對應(yīng)變率的大小沒有顯著影響,這表明魚類能在小流速和低湍動的水流條件下感知到相同的水流刺激向前運(yùn)動,可促進(jìn)魚類的洄游上溯;而改變豎縫隔板的安裝高度,對豎縫隔板上下水深處應(yīng)變率的分布規(guī)律無明顯影響。

圖11 不同水深處池室水力應(yīng)變率分布

3 結(jié) 論

a.對于不同的大底坡坡度,豎縫處設(shè)置隔板可有效減小池室最大流速及池室平均湍動能:池室最大流速下降率分布范圍為2.01%～21.49%,池室平均湍動能下降率分布范圍為3.77%～42.19%。

b.底坡為1/20時,池室主流從反S形轉(zhuǎn)變?yōu)槔贸厥艺w容積進(jìn)行低湍流能量轉(zhuǎn)換并自身回轉(zhuǎn)的Ω形,并且其主流最大流速衰減率達(dá)到61.24%～61.82%。

c.改變豎縫隔板安裝高度對隔板上下水深處的豎縫流速分布存在一定的影響,而對流速、湍動能、水力應(yīng)變率的分布規(guī)律無顯著影響,因此為了使魚類能夠更加順利地通過豎縫,可將豎縫隔板安裝在豎縫的中上方水深位置。

d.在豎縫處設(shè)置隔板以優(yōu)化大底坡魚道形成的高流速、高湍動的水流流場具有一定的可行性,但是還需進(jìn)一步結(jié)合過魚試驗(yàn)研究來評估魚類對這些水動力條件的上溯行為反應(yīng),以論證研究結(jié)果的工程應(yīng)用價值。