基于Flow-3d 的異側豎縫底孔組合式魚道水力特性研究

范聰喆,劉 科,劉彬俠,陳 卓

（陜西省江河水庫工作中心,陜西 西安 710018）

1 前言

隨著人類社會發展,水利工程在防洪、發電、灌溉等方面發揮著巨大作用,但同時對生態環境也造成了一定影響[1]。特別是水庫大壩的修建阻礙了魚類洄游上溯進行基因交流,破壞了河流生態系統的穩定性[2]。為了緩解水壩對魚類洄游通道的阻隔作用,在水利工程內部設置魚道等過魚設施,形成連接通道,可保障魚類洄游、促進上下游種群和基因交流[3]。近年來,隨著國家對生態環境保護的重視,河流筑壩開發時修建過魚設施成為必不可少的環節[4]。

魚道可分為單一式魚道和組合式魚道,單一式魚道按體型又可分為豎縫式魚道、孔口式魚道及溢流堰式魚道[5]?，F有研究多為單一式魚道,但洄游性魚類眾多,其溯游習性各異,傳統的單一式魚道不能同時滿足多種魚類的溯游習性,研究組合式魚道就越發重要[6]。同時,大多數魚道的水力設計只考慮流速因素,忽略了魚道內水流的紊流特性。Herskin等認為[7],魚類成群溯游時,可從領頭魚產生的旋渦中受益。若魚道中存在這種生態友好型的水流形態,有利于魚類溯游,因此在研究魚道可行性時,應適當考慮紊動能變化。

本文利用Flow-3d 軟件模擬魚道的流態分布規律,分析魚道流速、紊動能與過魚對象的適應性,結合其水力特性,為孔口和豎縫組合式魚道的設計和優化提供參考。

2 研究對象和方法

2.1 研究魚道的基本體型與結構

豎縫與孔口組合式魚道在我國應用比較成功,某魚道全長1443.0 m,上下游水頭差15.0 m,魚道底坡為1∶80,魚道隔板203 塊,單個池室長度為6.0 m,有效水深3.0 m。魚道槽身標準斷面底寬5.0 m,兩側邊墻高2.5m,上接1∶2 斜坡高1 m,共高3.5 m,再接馬道,見圖1。隔板高度3.5 m,厚0.2 m,相鄰兩隔板的過魚孔交叉布置[8]。魚道每隔30 塊隔板設休息池,供魚類上溯時短暫休憩,休息池為平底段,長10 m,斷面與魚道槽身相同。隔板下底邊3.5 m,靠一側邊坡布置,邊坡底部設置1.5 m×1.5 m 過魚底孔,見圖2。

圖1 魚道槽身斷面圖（單位：m）

圖2 魚道槽身隔板斷面圖（單位：m）

2.2 流動控制方程和紊流模型

利用Flow-3d 軟件對魚道進行水動力數值模擬。流體動力學計算（CFD）的實質是在計算域內對水流控制方程進行離散迭代求解,獲得計算域內的水動力因子,數值模擬的理論基礎為流動控制方程。

（1）質量守恒方程

式中：xi、ui和ρ分別表示在i 方向上的空間坐標、流體流速和流體密度。

（2）動量平衡方程

式中：fi、ui分別表示平均壓強、i 方向上的重力分量和流體粘性,-是由Boussinesq 假設推得的紊流雷諾應力,采用k-ε方程,則具有如下表達式：

魚道模擬流體為不可壓縮流體,在選用紊流模型時,要考慮模擬局部回流與漩渦,根據各個紊流模型的特點,考慮采用計算精度和穩定性較好的k-ε模型[9-11]。

3 研究對象和方法

3.1 建立三維模型

由于研究魚道較長,結構相對復雜,為便于后期計算分析,現對魚道模型進行簡化,選取魚道10 個池室構建數值模型,魚道中部設置10 m 長休息室。為保證模擬水流進入魚道的穩定性,分別在魚道進口和出口處各設置10 m 的過渡段,隔板厚0.2 m,簡化魚道基本尺寸見圖3,利用Rhino 軟件建立三維模型見圖4。

圖3 模擬魚道內部結構俯視圖（單位：m）

圖4 魚道三維模型

3.2 網格劃分及邊界條件設置

根據魚道原型概況并結合Flow-3d 軟件模擬特性,為提高計算的穩定性、高效性和準確性,利用模塊中的FAVOR功能進行網格處理,建立由0.1 m×0.1 m×0.1 m 的均勻立方體網格單元構成的三維網格系統,笛卡爾坐標系下沿xyz 坐標軸正方向長度分別92 m、3.5 m 和8 m,其網格系統恰好涵蓋魚道模型,共有2576000 個網格單元。

通過相關原型資料分析確定模型邊界條件為：魚道進口段（上游邊界x-Min）設為流量邊界,流量為6.64 m3/s,水深為3 m；下游出口（下游邊界x-Max）設為壓力邊界,設定水深為3 m 的靜水定壓力；模型上部（上邊界z-Max）設置為壓力邊界,設置數值為1 個大氣壓空氣壓強；模型的底部（下邊界z-Min）設為固壁邊界,魚道壁面邊界條件符合流速不分離和無滑移條件,見圖5。

圖5 魚道模型邊界條件

結合本次模擬的實際情況確定其它相關參數：模擬流體選擇20℃的水,物理模型選擇重力模型（Gravity）、粘度與紊流模型（Viscosityd and Turbulence）,根據魚道坡降,重力加速度在z 軸負方向取-9.809 kg/s3,在x 軸正方向取0.123 kg/s3,粘度與紊流模型選擇無粘度的紊流RNG 模塊（Renomalized group model）,為提高計算效率,初始計算水位設為3 m,運行時間設置為水流狀態趨于穩定時停止模擬,輸出數據添加流速、紊動能等類型,數據輸出間隔為2 s。

3.3 模型驗證

為保證數學模型的有效性及合理性,進行相應水力模型的比對驗證,建立1∶10 幾何比尺模型,見圖6。采用紅外線流速儀對魚道內a-b 監測線（第五級池室的B/2、h/2 位置處,其中B 為池室寬,h 為有效水深）流速進行測定,圖7 為試驗實測流速與模擬流速對比,實測值與模擬值基本吻合、趨勢基本一致,說明所建數學模型合理,模擬成果具有一定可信性。

圖6 魚道模型試驗結構示意

圖7 模擬流速對比

4 結果分析

4.1 流速分析

因魚道內不同水深的流態差異較大,為更好地研究組合式魚道的水力特性,現將魚道沿xy 平面取三個水深剖面進行分析,分別為底層水體H1=0.1H=0.3 m；次底層水體H2=0.25H=0.75 m（底孔中部）；表層水體H3=0.6H=2.4 m；其中H 為自由表面水深,H1、H2、H3為距魚道底部高度,分析平面布設見圖8。

圖8 分析平面布設示意

為提高分析準確性,選取典型段（x=22.1～28.8）進行分析,即魚道第3 池室,三四隔板之間。圖9 為典型段不同水深的x-y 剖面流速分布圖,典型段各層水體流速統計見表1。底層水體（H1=0.3 m）平均流速為0.845 m/s,受過水斷面收縮影響,底孔流速增大,但底孔流速略小于豎縫流速,且底近孔邊壁側流速較近隔板側流速低,符合魚類洄游喜貼壁的生理特性。由于隔板阻擋,隔板后方形成小型低流速回流區（平均流速為0.15 m/s）,為魚類上溯途中提供小型短暫休憩區。次底層水體（H2=0.75 m）平均流速為0.692 m/s,底孔流速整體偏大,流場分布較底層均勻,流體大部分經下一級隔板潛孔流出,在兩級隔板之間生成回流區,在底孔下游也有小型回流區產生。表層水體（H3=2.4 m）平均流速為0.175 m/s,小于底層和中層分析平面平均流速,水體流態相對較為簡單,在隔板后形成較大漩渦,水流呈“s 型”軌跡繞過隔板。圖10 為全魚道x-z 剖面流速分布圖。由圖可知,各池室x-z 剖面流速分布與典型池室基本一致。根據上述分析,組合式魚道的設計能夠很好地滿足大部分魚類上溯。

表1 魚道水力參數統計

圖9 典型段不同水深下x-y 剖面流速分布圖

圖10 全魚道x-z 剖面流速分布圖

4.2 紊動能分析

紊動能可用于度量由于紊動流速波動而導致動能增加的程度,紊動能對魚類洄游影響較大,若紊動能過大影響魚類對于方向的判斷,削弱魚類游泳能力和平衡能力,增大魚類上溯的能量消耗,從而影響魚類上溯率,甚至會對魚類的生理能力造成不可逆損傷。因此合適的紊動能是魚道設計的重要參考指標,可為魚類洄游提供有利流態,保障魚類上溯條件。典型段不同水深下x-y 剖面紊動能分布見圖11。

由圖11 可知,對于底層水體（H1=0.3 m）的紊動能最高值可達到0.162 m2/s2,為三層水體最小,紊動能最大值分布在各級擋板的上游區域,擋板近豎縫側的下游也存在小范圍的紊動能值偏高區域,順水流方向帶狀分布。次底層水體（H2=0.75 m）的紊動能最高值可達到0.225 m2/s2,紊動能的分布規律與底層水體基本相同,較底層水體有所增加,近豎縫側的高紊動能區域較底層水體縮小。表層水體（H3=2.4 m）的紊動能最高值為0.212 m2/s2,高紊動能區域與主流軌跡部分重疊,回流區紊動能較小。魚道池室內的三層水體紊動能基本在0.131 m2/s2以下,符合魚類洄游條件。全魚道x-z 剖面紊動能分布圖見圖12,由圖可知,各池室剖面紊動能流態分布與典型池室基本一致。魚道各層水體紊動能統計見表1。

圖11 典型段不同水深下x-y 剖面紊動能分布圖

圖12 全魚道x-z 剖面紊動能分布圖

5 結論

本文利用Flow-3d 對異側豎縫底孔式組合式魚道水動力特性進行數值計算,分析研究得出以下結論：

（1）池室流速隨水深增加逐漸增大,隔板后存在低流速回流區,為魚類洄游過程提供短暫休憩場所,表層流體流速較低,在隔板后形成較大漩渦,水流“s 型”依次繞過各級隔板流動。

（2）開孔高度內,紊動能隨水深增加而減小,在未開孔高度隨水深增加逐漸降低,魚道池室內大部分區域的紊動能均在0.131 m2/s2以下,符合魚類洄游條件,高紊動能區域與主流軌跡部分重疊,回流區紊動能較小。

綜上所述,組合式魚道能為魚類上溯提供適宜的水力條件,充分發揮魚道功能,對提高魚類的上溯率具有重要作用,本文的研究思路與方法可為魚道的設計及運行提供參考。