宋發(fā)軍

(重慶市奉節(jié)縣水利局，重慶 404600)

我國幅員遼闊，不同地區(qū)地形地貌差異較大，且氣候環(huán)境各異。對于我國云貴川等西南地區(qū)，多為山地和丘陵地貌，地形落差大，年均降雨量充沛，適合建造水電工程。不過，正是由于這些區(qū)域復雜的地形地貌，在修建水電工程時，安全問題也稱為了最受關注問題。

泄水建筑物是用來保證水電工程安全性的設施，其中泄洪道的合理設計能夠有效地保證泄洪水流通暢，降低空化現(xiàn)象。不同消能方式下泄洪道的消能效果各不相同、功能各異[1-5]。例如，底流泄洪的流態(tài)穩(wěn)定，具有較好的消能效果，能夠適應高中低水頭，靈活應對不同的地質(zhì)環(huán)境和水位變化；尾坎式消流也具有較好的消能效果，而且對河道兩岸及下流河道的影響較小；階梯消能可以強化泄洪水流的紊動，有效耗散水流的動能，提高摻氣量，延長溢洪道的使用壽命[6-10]。而且，對于較為復雜的地形地貌來說，常常會結(jié)合多種消能模式，保證泄洪道的消能效果，提高溢洪道的使用壽命[11-15]。

本文對不同消能方式下的彎曲溢洪道水流消能的效果進行了模擬，并對其消能效果進行了分析，對比不同消能方式下的消能能力的差距，在實際工程應用中具有指導意義。

1 消能率算法

溢洪道消能效果通常用消能率來衡量。本文設計彎曲溢洪道在不同消能方式下的效果，消能率越高，則說明該方式下彎曲溢洪道對水流的消能效果越好[16]。對于溢洪道消能率的計算方式各不相同，本文采用公式(1)的方法計算溢洪道的消能率：

(1)

式中，η為絕對消能率；ΔE為溢洪道入口處與出口處水流的能量差；E0為溢洪道入口處的水流能量；E1為溢洪道出口處的水流能量；V0和V1分別為溢洪道入口處和出口處水流的平均流速。

水流能量E可按照式(2)進行計算：

(2)

式中，H為水流斷面高度；α為水流流速系數(shù)；V為水流流速；g為重力加速度。

本文認為在溢洪道入口和出口處，水流斷面高度H和水流流速系數(shù)未發(fā)生明顯改變，因而對于消能率的計算可以簡化為式(1)。

2 數(shù)學模型的建立

2.1 流體體積函數(shù)

本文追蹤水氣界面的方法采用了VOF方法，分別定義水流中水的體積分數(shù)函數(shù)為αw，氣體的體積分數(shù)函數(shù)為αa，二者之和為1，即：

αw+αa=1

(3)

水的體積分數(shù)控制微分方程則可用式(4)進行表示：

(4)

式中，t為時間；ui為水流速度分量；xi為水流坐標分量；其中i可以為1或2。

通過對式(4)進行求解可以對水氣界面進行跟蹤，求解過程需要進行瞬態(tài)求解以獲得非恒定流解，然后通過時間逐步迭代獲得恒定流解。

2.2 紊流模型

溢洪道中尾坎處的水流具有各向異性和高速湍流的特點，為了更貼切真實的水流狀態(tài)，本文采用RNGk-ε模型進行了紊流數(shù)值模擬。該模型的連續(xù)方程、動量方程、k方程分別如式(5)—式(7)所示：

(5)

(6)

(7)

式中，μ為流水的分子黏性系數(shù)；ρ為平均密度；ui為水流速度分量；P為修正壓力；μt為流水的紊動黏性系數(shù)，可按照公式(8)計算；Gk和Cε為由于不同位置平均水流速率的梯度導致產(chǎn)生的紊動能項，分別可按照公式(9)和公式(10)進行計算：

(8)

(9)

(10)

式中，Cμ為0.084 5；C1ε為1.42；C2ε為1.68；σk為0.72；σε為0.72；η0為4.38；β為0.012。

在利用上述模型進行數(shù)值模擬計算時，將模型進行簡化，利用二維平面模型進行計算。邊界條件包括入口邊界和出口邊界。其中，入口邊界條件包括水流入口邊界和空氣入口邊界；水流入口邊界條件為入口處平均流速，為給定值；空氣入口邊界條件為大氣壓，亦為給定值。出口邊界條件為非給定值，用函數(shù)形式進行表示，認定計算域內(nèi)充滿水，對溢洪道底層水流進行標準壁函數(shù)處理，將邊墻水流定義為無滑移邊界條件。計算過程進行離散計算，采用控制體積法與非結(jié)構(gòu)網(wǎng)格進行結(jié)合。在每個控制體積中進行微分方程的積分，得到積分方程，然后進行線性化處理，得到各個變量。壓力—速率的耦合求解利用適用于瞬態(tài)求解的PISO算法進行求解。

3 尾坎挑流消能研究

本文以某水庫為模擬的工程模型，所模擬的洪峰為1 450 m3/s，進水口的流量為125 m3/s，泄洪沖沙閘流量為490 m3/s，泄洪閘流量為445 m3/s。當河道流量超過1 140 m3/s時，開始溢洪。本模型中溢洪道尾坎的平均流速見表1。尾坎挑流消能的模式種類繁多，不同的尾坎挑流消能效果也各不相同。本文通過數(shù)值模擬的方法研究差動式尾坎挑流消能和連續(xù)式尾坎挑流消能2種方式下的彎曲溢洪道水流消能效果。為了更為準確地對比2種消能方式下的消能效果，本文不僅對其各自的挑流消能率進行了計算，還對比了不同消能方式下水流的挑流水舌挑距、溢洪道各部位的水流速率。

表1 溢流道不同位置流速Tab.1 Flow velocities at different positions of overflow passage

本文采用流體分析軟件Flow-3D建立了溢洪道模型，采用網(wǎng)格劃分方法將模型分為958 127個網(wǎng)格，采用紊流模型，進口邊界條件采用進水量，出口邊界為自由邊界，邊界條件為大氣壓力，本文中的彎曲溢洪道為混凝土結(jié)構(gòu)，因此洪道邊界為無滑移速度邊界條件，數(shù)值模型如圖1所示。

圖1 計算模型Fig.1 Calculation model

2種尾坎挑流消能方式下，彎曲溢洪道水流的水舌挑距見表2。

表2 不同消能方式下水舌挑距Tab.2 Launching tongue spacing in different energy dissipation modes

由于挑流尾坎方式的影響，水流的形式發(fā)生了變化，水舌挑距也出現(xiàn)了明顯的差別。差動式尾坎挑流方式下挑流水舌挑距為86 m，連續(xù)式尾坎挑流方式下挑流水舌挑距為90 m。由于發(fā)生水流挑射的主要原因是溢洪道中水流的動能所導致的，差動式尾坎挑流方式下水流受到差動式尾坎的影響，水流發(fā)生碰撞，消耗掉部分水流動能，因此導致挑距有所降低；相對來說，連續(xù)式尾坎方式下，水流碰撞的幾率降低，未發(fā)生明顯的動能消耗，挑距較大。

為了研究不同尾坎挑流方式下水流平均速率的變化，采用數(shù)值分析的方法得到了差動式尾坎挑流方式下溢洪道中水流平均流速的等值線如圖2所示，連續(xù)式尾坎挑流方式下溢洪道中水流平均流速的等值線如圖3所示。2種尾坎挑流方式下，水流的平均流速變化趨勢沒有明顯差別，均是在溢洪壩堰頂出現(xiàn)最小的流速，在挑流水舌處出現(xiàn)了最大流速。不同的尾坎挑流方式下，溢洪壩堰頂處的流速未出現(xiàn)明顯變化，均為5 m/s，也是2種方式下的水流最小流速；但是，2種方式下挑流水舌入水處的流速個不相同，其中差動式尾坎挑流消能方式下的最大流速為35 m/s，連續(xù)式尾坎挑流消能方式下的最大流速為36 m/s。

圖2 差動式尾坎挑流消能方式下的平均流速等值線Fig.2 Contour map of average velocity under mode of differential tail sill flip flow energy dissipation

圖3 連續(xù)式尾坎挑流消能方式下的平均流速等值線Fig.3 Contour map of average velocity under condition of continuous wake bucket flip bucket energy dissipation

本文還利用公式(1)對2種尾坎挑流消能模式下的消能率進行了計算，結(jié)果見表3。2種方式下，溢洪道入口處水流的能量相同，但溢洪道出口處的水流能量出現(xiàn)了明顯的差別。差動式尾坎挑流消能方式下的消能率為21.7%，連續(xù)式尾坎挑流消能方式下的消能率為17.9%。相比之下，差動式尾坎挑流消能方式下的消能率較高，比連續(xù)式尾坎挑流消能方式高3.8%。這是由于，在水流進入到溢洪道中后，流水與溢洪道壩面之間發(fā)生碰撞，從而造成了動能的消耗。相比連續(xù)式尾坎挑流消能方式來說，差動式尾坎挑流消能方式下流水與壩面碰撞更為強烈，消能效果較好，且對下游和兩岸的影響較小。

表3 不同消能方式下消能率Tab.3 Energy dissipation rate under different energy dissipation modes

由表3可知，差動式尾坎挑流方式消能率更高，因此對該消能方式下溢洪道中心沿程壓力分布從模擬結(jié)果中進行了提取，壓力分布結(jié)果如圖4所示。從圖4可知，通過計算得到的壓力曲線與數(shù)值模擬曲線較為接近，二者的壓力最小值在出口附近，壓力最大值在進口附近。

4 工程實踐應用效果

重慶李家口水電站位于重慶市巫山縣的清溪口河，水電站上距孔巖廟6.5 km，下離張皮村11.7 km。工程河段控制流域面積約253.9 km2，年均徑流量約1.32億m3，設計洪水位(根據(jù)水電樞紐等級取1%)為547.84 m，正常蓄水位為546.00 m。該水電站溢洪道下泄水流流速過大，對左岸邊壁有較明顯的沖擊。采用差動式尾坎挑流消能方式后，有效改變彎道底部部分水流的方向，由于動量交換的機械作用，底部水流動量變化傳給上部水流，使上部水流動量在整個方向上改變一個角度，部分水流從凹岸流向凸岸，對凹岸和凸岸水流水深平衡有顯著促進作用。挑流消能效果分析見表4。本工程挑流消能段的消能效率高達25.7%，有效消耗了下泄水流能量，減小了下泄水流流速，消能效率良好，對保證下游河床穩(wěn)定有積極作用。

圖4 壓力分布Fig.4 Pressure distribution

表4 挑流消能效果Tab.4 Energy dissipation effect of flip flop

5 結(jié)論

本文利用VOF方法對溢洪道中水流的水氣界面進行了追蹤，采用RNGk-ε模型進行了紊流數(shù)值模擬，并基于Flow-3D軟件對水流平均流速進行了模擬，研究并對比了差動式尾坎挑流和連續(xù)式尾坎挑流兩種方式下，彎曲溢洪道水流消能的效果。結(jié)果表明，相比連續(xù)式尾坎挑流方式，差動式尾坎挑流方式下，水流之間以及水流與壩面之間的碰撞更為激烈，導致這種方式下挑流水舌挑距更低，水舌入水處的最大流速更小，消能率更大。研究結(jié)果表明，差動式尾坎挑流方式下，挑流水舌挑距為86 m，最大流速為35 m/s，消能率為21.7%，是更為有效的消能方式。

參考文獻(References)：

[1] 楊金孟.溢洪道泄槽彎道透水丁壩布置試驗研究[D].泰安：山東農(nóng)業(yè)大學，2019.

[2] 高學平，井書光，賈來飛.溢洪道彎道水流影響因素研究[J].水力發(fā)電學報，2014，33(4)：132-138.

Gao Xueping，Jing Shuguang，Jia Laifei.Study on influencing factors of flow in spillway curve[J].Journal of Hydroelectric Engineering，2014，33(4)：132-138.

[3] 張峰，劉韓生.臺階式溢洪道消能特性的研究[J].水利與建筑工程學報，2012，10(4)：11-13.

Zhang Feng，Liu Hansheng.Study on energy dissipation characteristics of stepped spillway[J].Journal of Water Resources and Architectural Engineering，2012，10(4)：11-13.

[4] 向蕾.溢洪道挑流消能段水氣兩相流的數(shù)值模擬[J].人民長江，2016(S1)：161-163.

Xiang Lei.Numerical simulation of water-gas two-phase flow in the energy dissipation section of the spillway[J].Yangtze River，2016(S1)：161-163.

[5] 薛宏程，刁明軍，岳書波，等.溢洪道出口斜切型挑坎挑射水舌三維數(shù)值模擬[J].水利學報，2013，44(6)：703-709.

Xue Hongcheng，Diao Mingjun，Yue Shubo，et al.Three-dimensional numerical simulation of slanting water tongue at spillway outlet[J].Journal of Hydraulic Engineering，2013，44(6)：703-709.

[6] 張志昌，徐嘯.臺階式溢洪道摻氣水流水深的試驗研究[J].應用力學學報，2011，28(6)：644-648.

Zhang Zhichang，Xu Xiao.Experimental study on the depth of aerated flow in stepped spillway[J].Chinese Journal of Applied Mechanics，2011，28(6)：644-648.

[7] 補舒棋.階梯溢洪道紊流數(shù)值模擬及消能影響因素分析[D].烏魯木齊：新疆農(nóng)業(yè)大學，2010.

[8] 包中進，王月華.基于FLOW-3D軟件的溢洪道三維水流數(shù)值模擬[J].浙江水利科技，2012，(2)：5-9.

Bao Zhongjin，Wang Yuehua.Three dimensional numerical simulation of spillway flow based on FLOW-3D software[J].Zhejiang Hydrotechnics，2012，(2)：5-9.

[9] 董玉蓮.彎曲溢洪道多級跌水與底流復合消能研究[D].大連：大連理工大學，2013.

[10] 張峰，劉韓生.臺階式溢洪道消能特性的研究[J].水利與建筑工程學報，2012，10(4)：11-13.

Zhang Feng，Liu Hansheng.Research on energy dissipation characteristics of stepped spillway[J].Journal of Water Resources and Architectural Engineering，2012，10(4)：11-13.

[11] 莊佳，張炯，高勝杰.廣東省郁南縣向陽水庫除險加固工程溢洪道消能試驗研究[J].中國水運，2013，13(9)：182-184.

Zhuang Jia，Zhang Jiong，Gao Shengjie.Experimental study on spillway energy dissipation in Xiangyang reservoir reinforcement project of Yunan County，Guangdong Province[J].China Water Transport，2013，13(9)：182-184.

[12] 周亮.壩體溢洪道設計研究[J].水能經(jīng)濟，2017(7)：65.

Zhou Liang.Research on design of spillway in dam[J].Journal of Hydropower Economics，2017(7)：65.

[13] 邱紅.臺階式溢洪道相對消能率研究[J].東北水利水電，2019，37(7)：18-22.

Qiu Hong.Study on the relative energy dissipation rate of step spillway[J].Water Resources & Hydropower of Northeast China，2019，37(7)：18-22.

[14] 李松泉.溢洪道多級消能優(yōu)化研究[J].城市建設理論研究(電子版)，2015，5(36)：4444-4445.

Li Songquan.Multi-level energy dissipation optimization of spillway[J].Urban Construction Theory Research(Electronic Edition)，2015，5(36)：4444-4445.

[15] 龔明正，付旭輝，楊勝發(fā).南川金佛山水電站溢洪道消能結(jié)構(gòu)優(yōu)化試驗研究[J].水電能源科學，2017，35(7)：121-124.

Gong Mingzheng，F(xiàn)u Xuhui，Yang Shengfa.Experimental study on energy dissipation structure optimization of spillway of Jinfo Mountain Hydropower Station in Nanchuan[J].Water Resources and Power，2017，35(7)：121-124.

[16] 田嘉寧.臺階式泄水建筑物水力特性試驗研究[D].西安：西安理工大學，2005.