廖智穎，黃春華

（珠江水利委員會(huì)珠江水利科學(xué)研究院,廣東 廣州 510611）

1 引言

長塘水庫位于廣西壯族自治區(qū)桂林市境內(nèi)、洛清江一級(jí)支流西河干流上,擬建壩址在永?？h永福鎮(zhèn)長塘村西河峽谷出口附近,距桂林市約54 km,距永福縣城約7 km。長塘水庫工程由水庫工程、輸水工程和灌區(qū)工程三大部分組成,是一座以供水、灌溉為主,結(jié)合防洪,兼顧發(fā)電等綜合利用的大（2）型水庫工程。

水利水電工程的開發(fā)建設(shè)在帶來了巨大的經(jīng)濟(jì)效益的同時(shí)[1],也破壞了河流連通性[2],必須采取針對(duì)魚類的生態(tài)補(bǔ)償措施[3]。豎縫式魚道作為協(xié)助魚類上溯的生態(tài)補(bǔ)償措施,近年來得到了廣泛的應(yīng)用[4-5]。

2 長塘水利樞紐短魚道工程

2.1 目標(biāo)魚類基本情況

長塘水利樞紐工程確定大刺鰍、圓吻鲴、大眼華鳊、大鰭鳠、波紋鱖、黃顙魚、光倒刺鲃等為主要的過魚對(duì)象,個(gè)體大小約為10 cm～30 cm。根據(jù)過魚種類的洄游特征[6],過魚設(shè)施運(yùn)行時(shí)間暫定過魚季節(jié)主要為4月～7月。

魚類突進(jìn)游速是豎縫式魚道設(shè)計(jì)中的重要參數(shù)[7],是推測(cè)目標(biāo)魚類能否通過魚道豎縫上溯的關(guān)鍵,魚道豎縫設(shè)計(jì)流速應(yīng)小于魚類的突進(jìn)游速[8]。根據(jù)現(xiàn)場(chǎng)魚類游泳能力測(cè)試試驗(yàn)結(jié)果,大刺鰍、圓吻鲴、大眼華鳊、大鰭鳠、波紋鱖的突進(jìn)游速分別位于 0.33 m/s～0.49 m/s、1.09 m/s～1.32 m/s、0.85 m/s～1.11 m/s、1.34 m/s～1.70 m/s、1.10 m/s～1.50 m/s。根據(jù)汪玲瓏、田凱等人[9-10]的研究,光倒刺鲃的臨界游速為（1.13±0.25）m/s,瓦氏黃顙魚臨界游速約為0.98 m/s,一般而言突進(jìn)游速大于臨界游速。因此7種過魚對(duì)象的突進(jìn)游泳速度范圍在0.42 m/s～1.50 m/s間,建議長塘水庫的短魚道豎縫流速設(shè)計(jì)值取1 m/s。

2.2 魚道設(shè)計(jì)方案及相關(guān)參數(shù)

綜合考慮西河魚類需求及各類過魚措施的過魚效果[11-13],本工程采用短魚道+固定門架啟閉機(jī)提升過壩方案來實(shí)現(xiàn)過魚。

短魚道由誘魚系統(tǒng)、進(jìn)口、普通池、休息池、趕魚柵池、集魚池和消力池等組成,見圖1。短魚道總長58.40 m,進(jìn)口閘門控制段長5.4 m（i=0）,直線段長21.6 m（i=2.5%）,后接132°的轉(zhuǎn)彎段,隨后上游集魚槽均為平坡段（i=0）。消力池池寬2 m,池長7.4 m。

圖1 長塘水庫工程位置示意圖

魚道池室結(jié)構(gòu)采用同側(cè)豎縫式,豎縫寬0.365 m。根據(jù)本工程目標(biāo)魚類情況,魚道寬度取2 m,魚室長度為2.7 m,沿程共設(shè)7個(gè)魚室。

整個(gè)系統(tǒng)的運(yùn)行過程為：潛水泵引水管從尾水渠引水,水流順魚道流向下游,最終由魚道誘魚口流出。魚類進(jìn)入誘魚口沿魚道上溯游入趕魚柵池；通過趕魚柵將魚類趕入集魚箱；利用集魚池上部排架的電動(dòng)葫蘆帶動(dòng)自動(dòng)抓取設(shè)備,將集魚箱提升并運(yùn)至運(yùn)魚軌道始端的運(yùn)魚車上,運(yùn)魚車沿著運(yùn)魚軌道將魚運(yùn)至壩腳,再由壩頂?shù)墓潭ㄩT架啟閉機(jī),吊運(yùn)集魚箱至壩前運(yùn)魚船,由運(yùn)魚船將魚運(yùn)送至上游庫尾,將魚類放生,完成魚類洄游過壩的需求。

圖2 短魚道立面圖

2.3 設(shè)計(jì)運(yùn)行水位

魚道進(jìn)口底板高程為142.0 m。為保證過魚效果,魚道內(nèi)需保持一定水深,143.0 m為最低運(yùn)行水位。泄洪期間魚道不運(yùn)行,故機(jī)組滿發(fā)時(shí)水位即為設(shè)計(jì)最高運(yùn)行水位,根據(jù)壩址流量-水位關(guān)系及機(jī)組滿發(fā)流量,確定最高運(yùn)行水位為143.64 m。

3 模型構(gòu)建

3.1 控制方程

本文研究的問題屬于不可壓縮流體的湍流流動(dòng),數(shù)值模擬選用RNG k-ε湍流模型,采用VOF模型追蹤自由液面,控制方程如下：

式中：ρ為流體密度；k為湍流動(dòng)能；ui為速度張量；ε為湍動(dòng)能耗散率；eff為擴(kuò)散系數(shù)；Gk為湍動(dòng)能 k的產(chǎn)生項(xiàng)；C1ε、C2ε均為經(jīng)驗(yàn)常數(shù)；ak、aε為常數(shù) 1.39；t為時(shí)間；xi、xj為坐標(biāo)張量。

B：VOF方程：

式中：F為流體體積分?jǐn)?shù)；Ax、Ay、Az分別為對(duì)應(yīng) x、y、z 方向的微元面積；u、v、w對(duì)應(yīng) x、y、z 方向的流體速度；為坐標(biāo)系數(shù),且= 0；FDIF為體積分?jǐn)?shù)擴(kuò)散項(xiàng)；FSOR為體積分?jǐn)?shù)密度源項(xiàng)；R為系數(shù),采用直角坐標(biāo)時(shí)為常數(shù)；UF為擴(kuò)散系數(shù)；CP為常數(shù),是紊流施密特?cái)?shù)的倒數(shù)；是動(dòng)能消散系數(shù)。

3.2 計(jì)算典型區(qū)段及網(wǎng)格劃分

本次僅對(duì)過魚設(shè)施的涉水建筑物部分進(jìn)行模擬計(jì)算,該模擬區(qū)段包括短魚道7個(gè)過魚池、趕魚柵池、集魚池和消力池,見圖3。三維模型X方向0～50.61 m,Y方向0～28 m ,Z方向0～17.79 m,用有限差分法對(duì)模型進(jìn)行離散求解,能大大減少網(wǎng)格數(shù)量,采用六面體結(jié)構(gòu)化網(wǎng)格,設(shè)2個(gè)網(wǎng)格區(qū)塊（圖3）,區(qū)塊 1 單元尺寸 0.05m, x∶y∶z=1∶1∶1,共計(jì) 9004800 個(gè)網(wǎng)格；區(qū)塊 2 單元尺寸 0.1m, x∶y∶z=1∶1∶1,共計(jì) 6048000 個(gè)。

圖3 短魚道+升魚機(jī)涉水部分三維整體模型

3.3 邊界及初始條件

區(qū)塊1的Xmin、Xmax、Ymin、Ymax均采用對(duì)稱S邊界,Zmin為墻體邊界,Zmax為空氣壓力邊界,給定大氣壓標(biāo)準(zhǔn)值；區(qū)塊Xmin為壓力出流邊界,根據(jù)不同工況控制水深分別為1.64 m和1.00 m。Zmin為墻體邊界,Zmax為空氣壓力邊界,Xmax、Ymin、Ymax均采用對(duì)稱S邊界。進(jìn)水端,補(bǔ)水管采用DN500管,在該處設(shè)置源,給定入流流量。為加快模型計(jì)算收斂,在區(qū)塊2范圍設(shè)置初始水體區(qū)域,初始水體水深與魚道進(jìn)口處水深一致。

圖4 模型邊界設(shè)置

3.4 模擬工況

工程為從電站尾水處抽取水體對(duì)魚道內(nèi)進(jìn)行補(bǔ)水。為形成魚道內(nèi)較好的水流條件,補(bǔ)水流量應(yīng)根據(jù)魚道進(jìn)口處水位進(jìn)行調(diào)整。本次主要分析魚道進(jìn)口處水位分別為機(jī)組滿發(fā)時(shí)水位143.64 m與魚道最低運(yùn)行水位時(shí)143.00 m的魚道內(nèi)部水流情況,具體見表1。

表1 三維模擬計(jì)算工況

4 結(jié)果分析

為詳細(xì)了解魚道內(nèi)部流場(chǎng)情況,根據(jù)水深情況本項(xiàng)目對(duì)工況1、2按Z=4.8 m進(jìn)行流場(chǎng)進(jìn)行剖切,對(duì)工況3按Z=4.4 m進(jìn)行剖切。短魚道、趕魚槽、集魚槽及消力池三維整體流速分布見圖5。

圖5 工況1三維整體流速分布

4.1 機(jī)組滿發(fā)時(shí)魚道流場(chǎng)分析

（1）工況1

3臺(tái)機(jī)組滿發(fā),當(dāng)補(bǔ)水流量為0.80 m3/s時(shí),模擬段最大流速為4.21 m/s,主要分布在消力池中。該流量下,消力池與集魚槽隔墻上方水深厚約0.10 m,集魚槽部分流速較小約為0.3 m/s～0.4 m/s,紊動(dòng)能為0.02 m2/s2,集魚槽道趕魚槽的突變段紊動(dòng)能為0.06 m2/s2～0.05 m2/s2（圖8）,一般認(rèn)為紊動(dòng)能0.05 m2/s2對(duì)魚類的影響較小。趕魚槽水流較為平順,與魚道連接段水流無明顯突變。而短魚道中豎縫的流速最高達(dá)1.3 m/s,其它多處于1.20 m/s,高于設(shè)計(jì)流速1 m/s,且流線不平順,豎縫處水有對(duì)沖墻體現(xiàn)象,不利于魚類洄游上溯。

圖6 工況1切片流場(chǎng)

圖7 工況1短魚道切片流場(chǎng)

圖8 工況1趕魚槽、集魚槽及消力池切片流場(chǎng)及紊動(dòng)能分布

（2）工況2

3臺(tái)機(jī)組滿發(fā),當(dāng)補(bǔ)水流量為0.65 m3/s時(shí),模擬段最大流速為3.43 m/s,主要分布在消力池中。該流量下,消力池與集魚槽隔墻上方水深厚約0.05 m,有一定的跌水現(xiàn)象。集魚槽部分流速有所增大,約為0.7 m/s～0.8 m/s,紊動(dòng)能為0.02 mm2/s2～0.04 m2/s2,集魚槽道趕魚槽的突變段紊動(dòng)能為0.06 m2/s2～0.05 m2/s2（圖11）。趕魚槽水流較為平順,與魚道連接段水流無明顯突變。短魚道中豎縫的流速最高達(dá)1.10 m/s,其它多處于1.00 m/s,與設(shè)計(jì)流速1m/s相符,流線相對(duì)完整,有利于魚類洄游。機(jī)組滿發(fā)下,補(bǔ)水流量0.65 m3/s較為理想。

圖9 工況2切片流場(chǎng)

圖10 工況2短魚道切片流場(chǎng)

圖11 工況2趕魚槽、集魚槽及消力池切片流場(chǎng)及紊動(dòng)能分布

4.2 最低運(yùn)行水位時(shí)魚道流場(chǎng)分析

下游水位為143 m時(shí)是魚道運(yùn)行的最低水位,通過設(shè)置補(bǔ)水流量0.65 m3/s,分析內(nèi)部流場(chǎng)特征。由圖可知,此時(shí)模擬段最大流速為5.73 m/s,由于流量較小,消力池與集魚槽間出現(xiàn)跌水現(xiàn)象,此處的流速最高。集魚池部分流速有所增大,達(dá)為 2.7 m/s～2.8 m/s,紊動(dòng)能為 0.24 m2/s2～0.28 m2/s2,但主要集中在表面（圖14）。趕魚槽水流較為平順,與魚道連接段水流無明顯突變。短魚道中豎縫的流速最高僅0.8 m/s,其它多處于0.7 m/s,流線相對(duì)平滑,水深也基本滿足魚類洄游需求,該工況下補(bǔ)水流量0.30 m3/s較為理想。

圖12 工況3切片流場(chǎng)

圖13 工況3短魚道切片流場(chǎng)

圖14 工況3趕魚槽、集魚槽及消力池切片流場(chǎng)及紊動(dòng)能分布

5 結(jié)論與建議

（1）通過三維數(shù)值模擬計(jì)算,魚道在電站滿發(fā)水位運(yùn)行時(shí),補(bǔ)水流量宜小于0.65 m3/s,此時(shí)魚道豎縫最大流速為1.10 m/s,多為1.00 m/s,流線較為平順,有利于大部分過魚對(duì)象上溯；魚道在下游最低運(yùn)行水位143 m運(yùn)行時(shí),補(bǔ)水流量宜為0.30 m3/s,此時(shí)豎縫最大流速為0.80 m/s,多為0.70 m/s,流速平順,除大刺鰍外,可兼顧本工程的所有過魚對(duì)象克流上溯需求。綜上,魚道的運(yùn)行補(bǔ)水流量應(yīng)維持在0.30 m3/s～0.65 m3/s間。

（2）同時(shí),因大刺鰍的突進(jìn)游速較小,游泳能力較弱,僅0.42 m/s,與其它過魚對(duì)象的游泳能力差異太大。設(shè)計(jì)工況下的豎縫流速對(duì)其會(huì)形成流速阻礙,建議專門制定大刺鰍過壩技術(shù)研究。