劉桐渤，閆 濱，王鐵良

(1.大連市莊河市水務(wù)局,遼寧大連 116400；2.沈陽農(nóng)業(yè)大學水利學院,沈陽 110866)

環(huán)形波紋鋼管涵洞式魚道水力特性數(shù)值模擬

劉桐渤1，閆 濱2，王鐵良2

(1.大連市莊河市水務(wù)局,遼寧大連 116400；2.沈陽農(nóng)業(yè)大學水利學院,沈陽 110866)

環(huán)形波紋鋼管涵洞式魚道內(nèi)波紋的粗糙度可以使邊界附近產(chǎn)生足夠低的流速,魚類可以由此游向上游。對3種流量0.05,0.07,0.09 m3/s,3種埋深0D,0.1D,0.2D,坡度0.4%工況下,環(huán)形波紋鋼管涵洞式魚道內(nèi)水面線、流速場和紊流場分布分別進行了數(shù)值模擬研究。結(jié)果表明:環(huán)形波紋鋼管涵洞中,高流速區(qū)(無量綱流速大于0.9)位于涵洞過水斷面中心區(qū)域,高紊流區(qū)(紊流強度大于0.2)位于水面中心附近。在涵洞底部邊壁附近較大區(qū)域內(nèi),由于流速較小,紊流強度也較低(紊流強度小于0.1),魚類可以由此順利完成上溯。相比無埋深式涵洞,嵌入式涵洞內(nèi)平均流速較小,紊流強度變化較為平緩,更有利于魚類的洄游。

數(shù)值模擬；水力特性；魚道；波紋鋼管涵洞；粗糙度；高流速區(qū)；高紊流區(qū)

2015,32(11):25-29,39

1 研究背景

大多數(shù)魚類都有在一定季節(jié)溯流而上遷移洄游的習性,這對于維持魚類種群和所處生態(tài)系統(tǒng)的穩(wěn)定至關(guān)重要。因而,在魚類棲息的河流與道路、堤防或其他溪流交匯處,常修建涵洞等交叉建筑物。

國外涵洞的設(shè)計一方面要求滿足輸水需求,另一方面必須保證魚類順利通過涵洞[1]。在加拿大,涵洞式魚道設(shè)計準則是獲得設(shè)計條件下涵洞內(nèi)特定的平均流速,且該流速不得超過當?shù)佤~類物種的游泳能力[2]。波紋鋼管涵洞過水斷面的大部分區(qū)域水流流速低于平均流速,適合魚類上溯洄游[3]。特別是波紋鋼管中波紋的粗糙度可以在邊界附近產(chǎn)生足夠低的流速,即使涵洞內(nèi)的平均流速超出了魚種的游泳能力,魚類仍然可以由此游向上游[4]。因而,波紋鋼管涵洞式魚道在國外應(yīng)用十分廣泛,過魚效果也比較理想。

適宜的水深、流速和紊流條件是輔助魚類順利通過魚道完成洄游的必要條件。將涵洞進口內(nèi)底置于正常河床高程以下(嵌入式涵洞)或使涵洞在回水條件下運行(淹沒式涵洞)可大大降低涵洞內(nèi)的水流流速,顯著增加涵洞內(nèi)低于平均流速的過水斷面面積,改善過魚條件[1]。若魚道內(nèi)水流紊動過大,魚類將會迷失方向,體能消耗過快,從而無法順利上溯[5]。常用紊動能、紊流強度等來描述紊流。Hunt等[6](2010)分析了環(huán)形波紋鋼管涵洞內(nèi)的紊動能分布,發(fā)現(xiàn)橫斷面內(nèi)最大紊動能位置與順水流向低流速區(qū)域的位置密切相關(guān)；涵洞進口附近紊動能顯著大于流場充分發(fā)展區(qū)域的紊動能。Abbs等[1](2007)研究發(fā)現(xiàn),環(huán)形波紋鋼管涵洞在邊界附近紊流強度較高,而水流中心區(qū)域紊流強度較低。最大紊流強度出現(xiàn)在涵洞進口過水斷面的上角偏下游處[7]；盡管這些區(qū)域的局部紊流強度值很高,當水流加速流入涵洞時,仍有一個很大的低紊流強度的中央核心區(qū)域；涵洞內(nèi)底嵌入河床不會降低涵洞進口處的紊流強度峰值,但是在距離進口約一倍洞徑以外紊流強度有所減小。Enders等[8]發(fā)現(xiàn)紊流強度會影響魚類游泳的消耗。試驗表明,紊流強度的影響約占幼鮭游泳總消耗的14%；隨著水溫、魚體質(zhì)量、平均流速和流速標準差的增加,幼鮭游泳總消耗量也在增加[9]。Smith等[10]利用幼年虹鳉魚進行試驗研究,發(fā)現(xiàn)在一定平均流速范圍內(nèi),魚類樂于聚集在低紊流區(qū)域；魚類常避開紊流強度較高的低流速區(qū)域；只有出現(xiàn)過大的平均流速時,魚類才會占據(jù)高紊流區(qū)。

波紋鋼管涵洞在我國高速公路(如青藏公路)、鐵路、堤壩輸水和舊橋涵加固中已有應(yīng)用,但尚未考慮兼做魚道使用。我國常用混凝土管或鋼筋混凝土管作為河道交叉建筑物,涵洞設(shè)計僅以滿足過流能力為準則,尚未考慮過魚需求。此種涵管的長期使用會對水體產(chǎn)生腐蝕作用；且其內(nèi)壁粗糙度不夠,流速過大,必將嚴重阻礙魚類洄游。而波紋鋼管涵洞造價低、材質(zhì)輕、強度高、壽命長、耐腐蝕、低碳環(huán)保、施工安裝便捷,且適于通車,因此在歐美等國家應(yīng)用非常廣泛。我國魚道的主要過魚對象為珍貴魚類、鯉科魚類和蝦蟹等。波紋鋼管涵洞內(nèi)波紋附近流速較低,適合我國魚類通行。隨著國家對生態(tài)環(huán)境保護措施扶持力度的不斷增強,波紋鋼管涵洞式魚道在我國必將有廣闊的應(yīng)用前景。

目前,國內(nèi)對波紋鋼管涵洞式魚道水力特性的研究比較少見。本文利用數(shù)值模擬方法,分析直徑500 mm的環(huán)形波紋鋼管涵洞在3種流量、3種埋深情況下,隨著水流向下游推進,魚道內(nèi)的水面線變化,以及無量綱流速和紊流強度的沿程分布,藉此深入剖析波紋鋼管涵洞內(nèi)水流的水力特性和魚類的生存環(huán)境,以便采取適宜的措施改善魚類生活環(huán)境,為維持水生態(tài)系統(tǒng)平衡,實現(xiàn)人與自然的和諧共處創(chuàng)造條件。

2 數(shù)值模擬控制方程及模型

本文采用標準k-ε紊流模型建立波紋鋼管涵洞式魚道的數(shù)學模型,對波紋鋼管涵洞式魚道內(nèi)的水流進行三維數(shù)值模擬。

2.1 標準k-ε紊流模型控制方程組

采用標準k-ε模型求解流動及換熱問題時,控制方程包括連續(xù)性方程、動量方程、能量方程、k方程和ε方程等。若不考慮熱交換的單純流場計算問題,則不需要包含能量方程。若考慮傳質(zhì)或有化學變化的情況,則應(yīng)再加入組分方程。這些方程都可表示成如下通用形式:

使用散度符號,上式記為

式中:?為通用變量,可以代表u,v,w,T等求解變量；Γ為廣義擴散系數(shù)；S為廣義源項。

三維直角坐標系下,與通用形式(1)對應(yīng)的k-ε模型的控制方程見表1。

2.2 VOF模型

水流進入涵洞后,由于過水面積變小,水流流態(tài)發(fā)生變化,水流在橫向、縱向和垂向都具有水流流速,水流與空氣的接觸面積特別大,流動屬于典型的水氣兩相流。VOF模型通過求解單獨的動量方程和處理穿過區(qū)域的每一流體的體積分數(shù)來模擬2種或3種不能混合的流體。因此在本文計算中采用VOF模型來處理水氣交界面。VOF模型的控制方程如下所示。

通用控制方程:

連續(xù)性方程:

動量方程:

能量方程:

表1 通用形式k-ε模型的控制方程Table 1 Control equations of k-ε model corresponding to the general form

2.3 模型網(wǎng)格與邊界條件

該數(shù)值模擬模型采用直徑為500 mm,長度為8 m,波紋尺寸為68 mm×17 mm的波紋鋼管,在3種埋深,1種坡度和3種流量下分別進行數(shù)值模擬。即:①坡度為0.4%；②埋深為無埋深(0D)、0.1D埋深、0.2D埋深,其中D為涵洞直徑；③流量為0.05, 0.07,0.09 m3/s。

魚道的結(jié)構(gòu)布置見圖1。魚道頂端有開孔,因此可以滿足通行魚類對光線的要求。

圖1 魚道結(jié)構(gòu)布置Fig.1 Structural layout of the fishway

模型網(wǎng)格采用線、面、體的順序進行劃分, Interval size的值分別設(shè)置為40,1 600,640 000 mm。

在邊界條件設(shè)置時,將波紋鋼管涵洞進口有水流部分設(shè)置成流速進口,涵洞進口空氣部分和涵洞內(nèi)水流上部切面設(shè)置為壓力進口,涵洞出口部分設(shè)置為壓力出口。

3 數(shù)值模擬結(jié)果與分析

依據(jù)上述數(shù)值模擬模型,對不同工況下波紋鋼管涵洞式魚道內(nèi)的水面線、流速場和紊流場進行了數(shù)值模擬,并對模擬結(jié)果進行了分析。

3.1 水面線

流量0.05 m3/s,坡度0.4%,埋置深度分別為0D,0.1D和0.2D情況下,涵洞內(nèi)縱向中心剖面水面線見圖2。

由圖2(a)知,在流量和坡度相同,不同埋置深度的情況下,涵洞內(nèi)水面線的走勢大致是相同的。3種工況下,水流在進入涵洞后,水面線均下降,在距離涵洞進口0.5 m處左右的位置水面線開始回升,并且,在距離涵洞進口4 m處左右,水面線趨于穩(wěn)定,涵洞內(nèi)水流流態(tài)穩(wěn)定,呈均勻流流態(tài),并在涵洞出口處自由出流。

圖2 不同埋深情況下涵洞內(nèi)水面線比較(流量0.05 m3/s，坡度0.4%)Fig.2 Water levels in the culvert of different embedded depths(flow rate 0.05m3/s，slope gradient 0.4%)

由圖2(b)知,在流量和坡度相同的情況下,嵌入式涵洞進口處水位雍高了。由于涵洞內(nèi)水流流量不變,水位雍高后,過水斷面面積增加,相應(yīng)地,涵洞內(nèi)平均流速降低?？梢?流量和坡度一定的情況下,嵌入式涵洞比無埋深式涵洞,更利于魚類上溯洄游。

涵洞無埋深(0D),坡度0.4%,流量分別為0.05, 0.07,0.09 m3/s情況下,涵洞內(nèi)水面線見圖3。其中粗實線,細實線和點劃線分別對應(yīng)流量0.05, 0.07,0.09 m3/s。

圖3 不同流量情況下涵洞內(nèi)水面線比較(埋深0D，坡度0.4%)Fig.3 Water levels in the culvert at different flow rates (non-embedded，slope gradient 0.4%)

由圖3可知,在涵洞坡度和進口埋深一定的情況下,流量越大,涵洞進口處水深也越大,涵洞內(nèi)水面線也越高。

3.2 沿涵洞縱向各斷面無量綱流速等值線云圖

為突出顯示涵洞內(nèi)流速較大(流速大于平均流速)區(qū)域,將斷面內(nèi)的流速值與正常水深下涵洞內(nèi)的流速值相比,就得到涵洞縱向各斷面無量綱流速等值線云圖。

流量0.05 m3/s,坡度0.4%,無埋深和0.1D埋深的涵洞內(nèi),與進口不同距離x處各斷面無量綱流速等值線云圖見圖4。

圖4 不同埋深情況下涵洞沿程橫斷面無量綱流速等值線云圖(0.05 m3/s，0.4%)Fig.4 Contours of dimensionless flow velocity at cross sections along the culvert of different embedded depths (flow rate 0.05m3/s，slope gradient 0.4%)

由圖4可見,相同流量和坡度時,無埋深(0D)和0.1D埋深涵洞對比,后者水位明顯壅高,邊壁附近無量綱流速明顯降低,較高流速(無量綱流速v＞0.9)區(qū)域面積顯著減小。可見,將涵洞進口埋置于天然河床高程以下形成嵌入式涵洞,可以有效地降低涵洞內(nèi)的水流流速,尤其是邊壁附近的流速,改善涵洞內(nèi)的水流流態(tài),有利于魚類的上溯洄游。

對嵌入式涵洞(圖4(a2至d2)),相同流量和坡度情況下,順水流方向,涵洞內(nèi)無量綱流速較大區(qū)域(v＞0.9)的面積逐漸減小,涵洞邊壁附近低無量綱流速區(qū)域面積逐漸增大,魚類可以由此游向上游。

對無埋深涵洞(圖4(a1至d1)),高無量綱流速區(qū)占據(jù)了過水斷面的大部分區(qū)域,僅在邊壁附近還存在較小范圍的低流速區(qū)；從涵洞進口到斷面x=1.5 m處,在水面附近存在流速較低區(qū)域。自斷面x=4 m處之后,無量綱流速等值線云圖幾乎不發(fā)生變化,表明涵洞內(nèi)流態(tài)已趨于穩(wěn)定。

對嵌入式涵洞,從進口起涵洞邊壁及水面附近水流流態(tài)充分發(fā)展,至斷面x=1.5 m處,水流流態(tài)漸趨平穩(wěn)；自斷面x=4 m處之后,涵洞內(nèi)水流穩(wěn)定,進入均勻流流態(tài)。

3.3 沿涵洞縱向各斷面紊流強度分布

流量0.05 m3/s,坡度0.4%,無埋深和0.1D埋深的涵洞縱向各斷面紊流強度等值線云圖見圖5。

圖5 不同埋深情況下涵洞沿程橫斷面紊流強度等值線云圖(0.05 m3/s，0.4%)Fig.5 Contours of turbulence intensity at cross sections along the culvert of different embedded depths (flow rate 0.05m3/s，slope gradient 0.4%)

由圖5可見,對于嵌入式涵洞,在進口處高紊流區(qū)位于涵洞過水斷面的左上方,隨著水流向下游推移,水流流態(tài)不斷發(fā)展,高紊流區(qū)逐漸向右、向上偏移,面積不斷擴大,到x=1 m斷面以后,高紊流區(qū)穩(wěn)定分布在水面中心位置。

對于無埋深的涵洞,僅在進口附近(x=0.25 m)高紊流區(qū)位于水面中心偏右方,之后,隨著水流向下游推移,高紊流區(qū)穩(wěn)定分布在水面中心位置。

不論進口有無埋深,沿涵洞縱向向下游,涵洞底部低紊流區(qū)(i＜0.1)的面積均在不斷擴大。直至斷面x=4 m處,低紊流區(qū)面積逐步趨于穩(wěn)定。沿涵洞底部邊壁附近較大區(qū)域內(nèi),由于流速較小,紊流強度也較低,魚類可以由此順利完成上溯。

從涵洞進口到斷面x=2 m處,有、無埋深涵洞在過水斷面上部紊流強度變化梯度差別不大；但從斷面x=3 m處起,嵌入式涵洞過水斷面上部紊流強度變化梯度明顯緩于無埋深式涵洞,即紊流變化較為平緩。這表明,涵洞進口嵌入河床高程以下可能會對4倍洞徑以后的涵洞內(nèi)紊流流態(tài)有所改善。

0.1D埋深情況,順水流方向,紊流強度在x=1 m處達到最大,imax=0.23；高紊流區(qū)面積(i≥0.2)在該處也達到最大。之后,較高紊流區(qū)面積逐漸減小,在x=3 m斷面附近高紊流區(qū)面積最小。然后,高紊流區(qū)面積又有所回升,直至斷面x=6 m以后,較高紊流區(qū)面積才逐步趨于穩(wěn)定。

對無埋深(0D)涵洞,從進口到x=0.75 m斷面附近,高紊流區(qū)(i≥0.2)面積逐漸增大。從x= 0.75 m斷面以后,高紊流區(qū)面積逐漸減小,在x=2 m處達到最小值。之后,高紊流區(qū)面積又逐漸擴大,直至x=6 m斷面以后,紊流才逐步趨于穩(wěn)定。

從涵洞進口到斷面x=2 m處,無埋深涵洞(0D)低紊流區(qū)(i≤0.1)面積大于0.1D埋深的涵洞。而在斷面x=3 m以后,正好相反,0.1D埋深涵洞內(nèi)低紊流區(qū)(i≤0.1)面積大于無埋深(0D)涵洞情況,有利于魚類的洄游。這再次表明,嵌入式涵洞對4倍洞徑以后的涵洞內(nèi)紊流流態(tài)會有所改善。

4 結(jié) 語

本文通過數(shù)值模擬分析了不同流量、埋深工況下,環(huán)形波紋鋼管涵洞式魚道內(nèi)的水面線變化、流速沿程分布以及紊流強度在進出口及洞身的變化趨勢。研究結(jié)果表明,涵洞式魚道過水斷面中心區(qū)域流速較高,紊流強度較低；水面中心附近紊流強度較高,水流流速較低。而在涵洞底部邊壁附近,流速較小,紊流強度也較低,最適合作為魚類上溯洄游的通道。與無埋深式涵洞相比,嵌入式涵洞內(nèi)平均流速較小,紊流強度的變化也比較平穩(wěn),有利于魚類的洄游。但嵌入式涵洞僅對4倍洞徑以后的涵洞內(nèi)紊流流態(tài)有所改善。

[1]ABBS T J.A Model Study of the Hydraulics Related to Fish Passage Through Backwatered Culverts[C]//Proceedings of the 18th Canadian Hydrotechnical Conference: Changes for Water Resources Engineering in a Changing World.Winnipeg,Manitoba,August 22-24,2007:1-12.

[2]BATES K.Design of Road Culverts for Fish Passage[R]. USA:Washington Department of Fisheries and Wildlife, 2003.

[3]BARBER M E,DOWNS R C.Investigation of Culvert Hydraulics Related to Juvenile Fish Passage[R].USA: Washington State Department of Transportation,1996.

[4]EAD S A,RAJARATNAM N,KATOPODIS C,et al.Turbulent Open-channel Flow in Circular Corrugated Culverts[J].Journal of Hydraulic Engineering,2000, 126(10):750-757.

[5]曹慶磊,楊文俊,陳 輝.2010.同側(cè)豎縫式魚道水力特性的數(shù)值模擬[J].長江科學院報,2010,27(7):26-30. (CAO Qing-lei,YANG Wen-jun,CHEN Hui.Numerical Simulation of Characteristics of Vertical Slot Fishway on Same One Side[J].Journal of Yangtze River Scientific Research Institute,2010,27(7):26-30.(in Chinese))

[6]HUNT M,CLARK S,KEHLER N.A Model Study of the Hydraulics Related to Fish Passage in a CSP Culvert with a Vertical Headwall[C]//Proceedings of the 8th International Symposium on Ecohydraulics.Seoul,Korea,September 12-16,2010:1-8.

[7]GARNER M E,KELLS J A,KATOPODIS C.A Model Study of the Hydraulics Related to Fish Passage through Embedded Culverts[C]//Proceedings of the Joint Conference of IAHR,ASCE and CSCE.Vancouver,August 9-14,2009:1-8.

[8]ENDERS E C,BOISCLAIR D,ROY A G.The Effect of Turbulence on the Cost of Swimming for Juvenile Atlantic Salmon(Salmo salar)[J].Canadian Journal of Fisheries and Aquatic Sciences,2003,60(9):1149-1160.

[9]ENDERS E C,BOISCLAIR D,ROY A G.A Model of Total Swimming Costs in Turbulent Flow for Juvenile Atlantic Salmon(Salmo salar)[J].Canadian Journal of Fisheries and Aquatic Sciences,2005,62(5):1079-1089.

[10]SMITH D L,BRANNON E L,ODEH M.Response of Juvenile Rainbow Trout to Turbulence Produced by Prismatoidal Shapes[J].Transactions of American Fisheries Society,2005,134(3):741-753.

(編輯:劉運飛)

Numerical Simulation on Hydraulic Characteristics of Culvert Fishway with Annular Corrugated Steel Pipe

LIU Tong-bo1,YAN Bing2,WANG Tie-liang2
(1.Water Affairs Bureau of Zhuanghe of Dalian City,Dalian 116400,China；2.College of Water Conservancy,Shenyang Agricultural University,Shenyang 110866,China)

The roughness of annular corrugated steel pipe fishway could generate sufficiently low flow velocities near the boundary,which facilitates fishes to immigrate towards the upstream.In this paper,the water surface line,flow velocity field,and turbulent flow field in the annular corrugated steel pipe culvert fishway are simulated under different conditions(flow rate 0.05 m3/s,0.07 m3/s and 0.09 m3/s,embedded depth 0D,0.1D and 0.2D,and slope gradient 0.4%).Results show that high velocity field(dimensionless velocity greater than 0.9)is located in the central area of the cross section of the culvert,and high turbulent flow field(turbulence intensity greater than 0.2) is located near the central area of the water surface.In large areas near the side walls of culvert bottom,due to relatively low velocity,the turbulence intensity is relatively low(less than 0.1)and fishes could complete migration using this passage.Compared with non-embedded culverts,the embedded culvert has lower average velocity and gentle variation of turbulence intensity,therefore it is more conducive to fish migration.

numerical simulation；hydraulic characteristics；fishway；corrugated steel pipe culvert；roughness；high velocity field；high turbulent flow field

TV135；X171.4

A

1001-5485(2015)11-0025-05

10.11988/ckyyb.20140304

2014-04-21；

2014-07-02

遼寧省農(nóng)業(yè)攻關(guān)計劃項目(2012212001)

劉桐渤(1987-),男,遼寧鞍山人,助理工程師,碩士,從事水工結(jié)構(gòu)工程相關(guān)方面的研究,(電話)15804269289(電子信箱) 15040330981@163.com。