戎貴文,楊緒婷,薄 飛,郭魯楠,賀 甜

(1.安徽理工大學(xué)地球與環(huán)境學(xué)院,安徽 淮南 232001;2.安徽省引江濟(jì)淮集團(tuán)有限公司淮南建管處,安徽 淮南 232001)

水電開(kāi)發(fā)在促進(jìn)社會(huì)發(fā)展的同時(shí)也對(duì)水生態(tài)環(huán)境產(chǎn)生不利影響,河流筑壩導(dǎo)致河流縱向連通性降低,魚(yú)類洄游受阻引發(fā)洄游性魚(yú)類的數(shù)量減少,甚至種群滅絕[1]。為減輕水利水電工程建設(shè)對(duì)魚(yú)類洄游的不利影響,通常在河流上修建過(guò)魚(yú)設(shè)施,輔助魚(yú)類完成生命周期。過(guò)魚(yú)設(shè)施對(duì)于有效恢復(fù)水生生物多樣性以及推進(jìn)生態(tài)文明建設(shè)具有現(xiàn)實(shí)意義。

魚(yú)道是一種供魚(yú)類洄游的過(guò)魚(yú)設(shè)施,豎縫式魚(yú)道因其能夠更好適應(yīng)水位變化、符合多類型魚(yú)的上溯需求,在國(guó)內(nèi)外被廣泛應(yīng)用。魚(yú)道流速是影響魚(yú)類通過(guò)效率的主要水力學(xué)指標(biāo),通過(guò)對(duì)豎縫式魚(yú)道結(jié)構(gòu)體型的研究,可以優(yōu)化魚(yú)道內(nèi)部的流速場(chǎng)分布。許多學(xué)者對(duì)豎縫式魚(yú)道作了大量研究,通過(guò)觀察魚(yú)道水流,發(fā)現(xiàn)過(guò)魚(yú)效率較高的魚(yú)道都具有滿足魚(yú)類輕松上溯的流態(tài)特征[2],而且主流位于中心位置、回流區(qū)面積接近的魚(yú)道水流對(duì)魚(yú)類洄游更為友好,使得魚(yú)類上溯消耗能量較少[3];魚(yú)道池室長(zhǎng)寬比為1.25時(shí),魚(yú)道池室具有較為穩(wěn)定的流態(tài)和可供魚(yú)類休息的回流區(qū)[4];魚(yú)道低坡度條件下,池室寬度與豎縫寬度之比為7時(shí),魚(yú)道池室會(huì)出現(xiàn)主流撞壁的不利流態(tài)現(xiàn)象[5];魚(yú)道內(nèi)出現(xiàn)回流區(qū)過(guò)于靠近下游豎縫的流態(tài),將會(huì)增加魚(yú)類被沖刷至下游的風(fēng)險(xiǎn)[6]。豎縫式魚(yú)道隔板型式對(duì)魚(yú)道水力特性影響十分顯著,在魚(yú)道池室水位差較大的情況下,能夠通過(guò)調(diào)整魚(yú)道隔板改變流速[7];豎縫式魚(yú)道池室采用流線型墩頭能夠減小主流偏轉(zhuǎn)角度,有利于主流流速沿程衰減[8];豎縫式魚(yú)道中的豎縫寬度明顯影響池室場(chǎng),豎縫寬度與池室寬度的比值在0.15～0.20范圍內(nèi),魚(yú)道內(nèi)縱向流速更加均勻且橫向流速值較小、流速條件較好[9]。水流的湍流特性對(duì)魚(yú)類上溯影響較大,過(guò)大的湍流動(dòng)能(Turbulent Kinetic Energy,TKE)會(huì)增加上溯阻力,影響魚(yú)類洄游的成功率[10-11]。豎縫寬度及位置變化對(duì)魚(yú)道流場(chǎng)和湍流動(dòng)能影響顯著,對(duì)魚(yú)道水力特性(流態(tài)、流速、紊動(dòng)能等)的影響機(jī)理仍有待進(jìn)一步研究。本文以江濟(jì)淮某工程魚(yú)道為研究對(duì)象,研究豎縫位置變化對(duì)魚(yú)道湍流動(dòng)能、流場(chǎng)和回流區(qū)的影響,探討池室水流結(jié)構(gòu)最佳時(shí)的魚(yú)道豎縫相對(duì)位置,以期為同類型魚(yú)道工程設(shè)計(jì)提供參考。

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

1.1 控制方程

RNGk-ε雙方程紊流模型應(yīng)用于魚(yú)道水流模擬已經(jīng)較為成熟,該模型既能提供水流的時(shí)均特性,又能更好地模擬彎曲流線和高應(yīng)變率的流動(dòng)。本文采用RNGk-ε雙方程紊流模型研究魚(yú)道水力特性,采用VOF(Volume of fluid)方法捕捉自由液面,控制方程如下:

連續(xù)方程

(1)

動(dòng)量方程

(2)

k方程

(3)

?方程

(4)

對(duì)自由表面的捕捉采用VOF方法

(5)

式中,ui為各方向速度分量,m·s-1;F為水或空氣所占體積百分比。若F=1,表示水充滿了整個(gè)計(jì)算區(qū)域;若F=0,表示計(jì)算區(qū)域內(nèi)是空氣,水與空氣各自所占百分比的總和為1;函數(shù)F是空間和時(shí)間的函數(shù),即F=F(x,y,z,t)。

1.2 數(shù)值方法及邊界條件

1)數(shù)值方法 采用有限體積法離散控制方程,空間離散采用基于單元體的最小二乘法,動(dòng)量方程、湍流動(dòng)能、湍流耗散率均采用二階迎風(fēng)格式,速度-壓力耦合計(jì)算采用壓力隱式算子分割法(PISO算法)[12]。

2)邊界條件 模型進(jìn)口和出口邊界條件分別設(shè)置為壓力入口和壓力出口邊界,并給定水深條件,以模擬魚(yú)道池室水流。其他壁面設(shè)置為無(wú)滑移邊界條件,模型與空氣交界面均采用壓力入口邊界。

2 計(jì)算區(qū)域及網(wǎng)格劃分

本文以引江濟(jì)淮某工程魚(yú)道為研究對(duì)象,按照1∶1比例進(jìn)行建模,模型共包含15個(gè)魚(yú)道普通池室、1個(gè)休息池。圖1為魚(yú)道池室平面圖,池室長(zhǎng)L為3m、寬度B為2.5m,休息池長(zhǎng)度是魚(yú)道池室長(zhǎng)度的兩倍,豎縫寬度為0.4m,豎縫法線與魚(yú)道中心線的夾角為45°,魚(yú)道內(nèi)部設(shè)同側(cè)隔板,隔板厚0.3m、長(zhǎng)隔板長(zhǎng)1.7m、短隔板長(zhǎng)P=0.65m。

圖1 魚(yú)道池室平面結(jié)構(gòu)圖(單位mm)

本文通過(guò)改變豎縫在池室中的相對(duì)位置,研究短隔板與魚(yú)道寬(P/B)之比為0.18、0.26、0.34、0.42、0.5和0.58的6種工況(見(jiàn)表1)不同豎縫位置對(duì)魚(yú)道池室水力特性的影響。

表1 研究工況

模型采用六面體網(wǎng)格進(jìn)行劃分,網(wǎng)格數(shù)為692 090,節(jié)點(diǎn)為754 344,如圖2所示。模擬過(guò)程迭代計(jì)算的殘差值為10-5,時(shí)間步長(zhǎng)為0.01s。

圖2 網(wǎng)格劃分圖

3 模型驗(yàn)證

利用物理模型試驗(yàn)結(jié)果對(duì)本文的數(shù)學(xué)模型進(jìn)行驗(yàn)證。物理模型試驗(yàn)中,在中間魚(yú)道池室內(nèi)選取7個(gè)橫斷面,每個(gè)橫斷面間隔為0.5m;每個(gè)斷面縱向設(shè)9個(gè)測(cè)點(diǎn),測(cè)點(diǎn)間隔為0.25m。池室共布置測(cè)點(diǎn)63個(gè),如圖3所示。

圖3 測(cè)點(diǎn)分布圖(單位mm)

選取主要斷面Ⅰ、Ⅲ、Ⅵ的流速實(shí)測(cè)值與數(shù)值模擬結(jié)果進(jìn)行驗(yàn)證(見(jiàn)圖4)。從圖4中可以看出,所有測(cè)點(diǎn)中實(shí)測(cè)值的最大流速為0.68m/s,最小流速為0.03m/s,模擬結(jié)果的最大流速為0.74m/s,最小流速為0.04m/s,表明數(shù)值模擬的流速值與試驗(yàn)結(jié)果基本吻合。因此,本文所構(gòu)建的數(shù)學(xué)模型是合理的。

(a)斷面Ⅰ (b)斷面Ⅲ (c)斷面Ⅵ

4 討論

采用所建立的數(shù)學(xué)模擬對(duì)6種工況分別進(jìn)行模擬計(jì)算,得到魚(yú)道池室的流場(chǎng)和湍流動(dòng)能,對(duì)不同工況下池室內(nèi)的水力特性進(jìn)行分析。

4.1 池室流場(chǎng)特性

豎縫位置的改變對(duì)魚(yú)道池室流場(chǎng)有直接影響[13]。圖5為6種工況水深h=0.8 m剖面流場(chǎng)。

(a)工況1 (b)工況2

從圖5可知,主流從豎縫以射流形式流入池室,射流沿?cái)嗝娉尸F(xiàn)先擴(kuò)散后收縮的趨勢(shì),流速也隨之先降低后增高,兩側(cè)都有流速較低的回流區(qū)。主流區(qū)的流速分布和衰減率與魚(yú)類在上游過(guò)程中的體能消耗、上溯時(shí)間和洄游成功率密切相關(guān)[14]。對(duì)比發(fā)現(xiàn),6種工況的主流距池室左壁面漸遠(yuǎn),橫向擴(kuò)散趨勢(shì)更明顯,彎曲程度加大。其中,工況1主流較為順直,偏轉(zhuǎn)角度小,主流靠近池室左側(cè);工況2、3主流彎曲角度比工況1大,主流呈“S”形;工況4、5、6的主流經(jīng)過(guò)上一個(gè)豎縫后向右側(cè)偏轉(zhuǎn)的角度較大,橫向?qū)挾葦U(kuò)大,到達(dá)下一豎縫處收縮,部分水流貼近池壁沖擊下一池室擋板,工況6水流貼壁現(xiàn)象更加明顯。由此可見(jiàn),改變擋板的長(zhǎng)度會(huì)對(duì)主流軌跡產(chǎn)生影響。這與文獻(xiàn)[15]洄游魚(yú)類通常沿著主流向上游游動(dòng)相一致。工況3中主流位于池室中央,這種流態(tài)對(duì)于魚(yú)類上溯有利,而工況5、6出現(xiàn)主流貼壁現(xiàn)象,在魚(yú)類上溯過(guò)程中可能會(huì)對(duì)魚(yú)體造成損傷,降低魚(yú)類存活能力。

根據(jù)以往的研究,洄游魚(yú)類需要以爆發(fā)泳速連續(xù)通過(guò)魚(yú)道豎縫,因此會(huì)消耗大量體能,需要適當(dāng)休息后繼續(xù)上溯[16]?；亓鲄^(qū)能夠?yàn)轸~(yú)類提供低流速休息環(huán)境,因此魚(yú)道回流區(qū)面積分布對(duì)魚(yú)類持續(xù)上溯十分重要。工況1中池室右側(cè)回流區(qū)面積較大,約為左側(cè)回流區(qū)的3倍;工況2、3左側(cè)回流區(qū)面積增加,工況3主流兩側(cè)回流區(qū)面積接近,分布較合理;工況4、5、6屬于同類型流場(chǎng),兩側(cè)的回流區(qū)分布與工況1相反,主流左側(cè)有大面積回流區(qū),且隨著豎縫位置右移,左側(cè)回流區(qū)面積依次增加,主流右側(cè)回流區(qū)面積被主流壓縮。由此可見(jiàn),豎縫位置距池室邊壁的距離應(yīng)該控制在一定范圍內(nèi),否則會(huì)導(dǎo)致兩邊回流區(qū)面積分布差異很大,魚(yú)類容易在大面積回流區(qū)中迷失方向,降低魚(yú)類的通過(guò)效率。因此,P/B的取值范圍控制在0.26～0.34之間時(shí),池室流態(tài)分布更適宜。

4.2 池室流速云圖

流速是魚(yú)道水力特性的重要參數(shù)之一,若池室流速大于魚(yú)類的爆發(fā)游速,魚(yú)類難以成功上溯。圖6展示的是6種工況下h=0.8m剖面的流速等值線云圖,6種工況均有明顯的流速分區(qū)現(xiàn)象。水流呈現(xiàn)多樣性,但流速分布差異不大,中層流速范圍0.10～1.0m/s,受主流偏轉(zhuǎn)影響,最大流速都出現(xiàn)在豎縫處及隔板前后;最大流速為0.9～1.0m/s左右,其中工況2～4主流流速峰值區(qū)面積較其他工況小。一般情況下,魚(yú)類上溯過(guò)程中會(huì)優(yōu)先選擇低流速區(qū)域通過(guò)[17]。故工況2～4的流速分布符合魚(yú)類洄游偏好。

(a)工況1 (b)工況2

隨著豎縫位置變化,流速較大的面積呈現(xiàn)先下降后增加的趨勢(shì),表明豎縫位置改變對(duì)主流最大流速分布有一定影響,但影響較小。文獻(xiàn)[18]對(duì)豎縫魚(yú)道體型作了相關(guān)研究,發(fā)現(xiàn)豎縫位置對(duì)主流流速的衰減及豎縫處流速分布影響不大,與本文的研究結(jié)果一致。豎縫位置靠近池壁時(shí),最大流速面積增大。同時(shí)6種工況的主流高流速區(qū)面積較小,僅占單個(gè)池室面積的1/4左右,低速區(qū)分布在主流區(qū)兩側(cè),與回流區(qū)重合,為魚(yú)類提供了一個(gè)休息空間。低流速區(qū)分布有所不同,工況1～3池室左側(cè)的低流速區(qū)最大值均為0.2m/s,主流右側(cè)回流區(qū)流速全部低于0.1m/s;工況4左側(cè)回流區(qū)最大流速為0.3m/s,且最大流速區(qū)域面積增加,右側(cè)回流區(qū)面積較小,最大流速為0.2m/s;工況5左側(cè)回流區(qū)最大流速為0.2m/s,右側(cè)回流區(qū)最大流速為0.1m/s,主流左側(cè)低流速區(qū)面積較大;工況6左側(cè)回流區(qū)最大流速為0.3m/s,右側(cè)回流區(qū)最大流速為0.2m/s。相比前3種工況,后3種工況回流區(qū)流速大于0.2m/s的區(qū)域增加。當(dāng)P/B的取值范圍控制在0.26～0.34之間時(shí),池室流速分布更有利于魚(yú)類上溯。對(duì)比不同工況下魚(yú)道池室內(nèi)的流速分布,發(fā)現(xiàn)存在友好流態(tài)時(shí),池室內(nèi)最大流速出現(xiàn)一定程度的降低,流態(tài)與流速之間的變化關(guān)系對(duì)于魚(yú)道設(shè)計(jì)具有參考價(jià)值。

4.3 池室湍流動(dòng)能(TKE)

研究表明,魚(yú)類在上溯中都偏好利用低湍流動(dòng)能區(qū)作為休息區(qū),然后再嘗試穿越更高流速和湍流區(qū)域,研究魚(yú)道池室TKE分布對(duì)提高魚(yú)道過(guò)魚(yú)效率具有重要意義[19-20]。6種工況h=0.8m剖面的TKE等值線云圖如圖7所示。

(a)工況1 (b)工況2

由7圖可知,整個(gè)池室的TKE的范圍在0～0.05m2/s2,隔板導(dǎo)角周圍的TKE值最大,并且沿射流方向衰減。主流區(qū)兩側(cè)的回流區(qū)TKE較小,回流區(qū)TKE范圍在0～0.02m2/s2的面積區(qū)域較大,流速與湍流動(dòng)能成正相關(guān),意味著在低速區(qū)域也應(yīng)提供低TKE值,主流流態(tài)的分布也存在很強(qiáng)的關(guān)聯(lián)。

隨著豎縫位置向池室中央靠近,湍流動(dòng)能大于0.02m2/s2的范圍縮小。其中,工況1、2的湍流動(dòng)能分布相似,靠近左側(cè)隔板的水流出現(xiàn)TKE的最大值0.05m2/s2,工況1中TKE≥0.02m2/s2區(qū)域的面積更大,沿豎縫前后延伸,兩種工況中TKE≥0.02m2/s2的面積占整個(gè)池室的1/4左右;工況3中湍流動(dòng)能最大值為0.03m2/s2,對(duì)比前兩種工況,TKE最大值下降約40%,除了擋板兩側(cè)分布有大于0.02m2/s2的區(qū)域,池室中其他區(qū)域TKE值均在0.01m2/s2左右;工況4的TKE整體低于0.02m2/s2,TKE是所有工況中最小的,其TKE極值區(qū)靠池室右側(cè)分布,可能是由于主流與右側(cè)池壁和擋板產(chǎn)生摩擦;工況5、6的TKE最大值為0.05m2/s2,相比前兩種工況有所增加,湍流動(dòng)能最大值位置出現(xiàn)在短擋板導(dǎo)角處,與工況1、2最大值出現(xiàn)位置基本一致,工況6的TKE變化更加劇烈。出現(xiàn)這種狀況是由于豎縫位置改變導(dǎo)致主流偏轉(zhuǎn)角增加,部分水流碰撞短擋板后導(dǎo)致紊流增大。綜上所述,豎縫位置改變對(duì)魚(yú)道池室TKE影響顯著,不僅改變了TKE分布的面積和位置,還影響其大小,豎縫位置位于池室偏中央處,TKE呈下降趨勢(shì);當(dāng)P/B取值范圍為0.34～0.42時(shí),擋板對(duì)池室TKE減弱效果更好。我國(guó)魚(yú)類大部分游泳能力較弱,會(huì)避開(kāi)魚(yú)道高流速和紊動(dòng)能的區(qū)域,對(duì)鳙魚(yú)和草魚(yú)的上溯行為研究發(fā)現(xiàn),它們偏好的TKE范圍為0.02～0.043m2/s2[21],而幼年草魚(yú)在上溯時(shí)會(huì)避開(kāi)TKE大于0.012m2/s2的區(qū)域[22]。本文研究認(rèn)為,P/B取值范圍為0.34～0.42時(shí)能夠減弱魚(yú)道內(nèi)湍流動(dòng)能,將魚(yú)道內(nèi)TKE最大值減小至0.02m2/s2,該結(jié)論對(duì)于魚(yú)道擋板設(shè)計(jì)有一定借鑒意義。

5 結(jié)論與展望

(1)豎縫相對(duì)位置改變對(duì)主流有明顯影響,P/B越大,主流彎曲程度越大。當(dāng)P/B≥0.42時(shí),主流出現(xiàn)貼壁現(xiàn)象,左側(cè)的回流區(qū)面積過(guò)大,不適宜魚(yú)類上溯;當(dāng)P/B在0.26～0.34范圍時(shí),池室內(nèi)流態(tài)分布較好。

(2)豎縫位置變化導(dǎo)致魚(yú)道池室的湍流動(dòng)能、回流區(qū)流速分布產(chǎn)生顯著變化,改變了魚(yú)道池室湍流的分布位置和區(qū)域;當(dāng)0.34≤P/B≤0.42時(shí),湍流動(dòng)能相對(duì)較小,能夠減少魚(yú)類上溯阻力。從湍流動(dòng)能、流速場(chǎng)分布綜合來(lái)看,當(dāng)P/B為0.34時(shí),池室內(nèi)的水流流態(tài)更有利于魚(yú)類洄游,能夠?yàn)轸~(yú)道優(yōu)化設(shè)計(jì)和修復(fù)魚(yú)類生境提供參考。

(3)在魚(yú)道的寬度和長(zhǎng)度確定的情況下,改變P/B比值是調(diào)節(jié)魚(yú)道主要水力學(xué)參數(shù)的一種最經(jīng)濟(jì)有效的方法,營(yíng)造適宜的水流環(huán)境對(duì)提高豎縫式魚(yú)道的運(yùn)行效果和生態(tài)效益非常重要。

對(duì)魚(yú)道水力特性造成影響的因素有很多,本研究立足豎縫相對(duì)位置變化對(duì)水力特性的影響,而實(shí)際天然水流環(huán)境更加復(fù)雜,流量的變化對(duì)魚(yú)道流場(chǎng)和湍流強(qiáng)度影響更大,未來(lái)可通過(guò)原型長(zhǎng)期觀察試驗(yàn)進(jìn)一步揭示魚(yú)道水力特性的影響機(jī)理。