杜家臺分洪口閘室消能防沖刷控制方法研究

吳佳倫

（湖北省漢江河道管理局杜家臺分洪閘管理分局，湖北 仙桃 433000）

引 言

消能是指在水利工程或泄水構筑物的建設過程中，為降低或防止由于上游水流對下游建筑物、河渠等造成的沖刷破壞問題，從而采取的消耗水流攜帶沖刷能量、分散水體向下沖擊力而采取的控制措施。由于閘室前后水位高程不等，水閘上、下游落差會將水體的勢能大量轉化為動能，導致過閘的水流流速大大提高，從而對下游河床及兩岸產生嚴重的機械沖刷和側向侵蝕，特別是底部沖刷引起的溯源侵蝕，更會使沖刷范圍從下游不斷后退到達閘室地基，進而造成水閘空懸甚至發生崩塌，危及水利工程或分洪口的安全[1]。因此，在水閘閘室的設計方案中加入控制下的消能防沖刷技術與措施，對于確保水利樞紐的工程安全至關重要。

根據現階段相關工作的實施進展可知，目前的設計成果還不能在水閘工程中發揮更好的作用，所采取的措施對消能防沖刷控制效果較差[2]。同時，在此相關工作的研究中發現，現有的消能防沖刷控制的設計成果較少，盡管相關消能防沖刷控制的物理模型與觀測模型已基本建成，但只適用于某地區特定的水閘工程，并不適合在更大范圍的水利工程領域推廣使用。為進一步深化消能防沖刷控制效果，本文以杜家臺分洪口閘室為例，設計一種全新的消能防沖刷控制方法，保證下游水力設計的合理規范。

1 閘室消能防沖刷控制目標函數設計

目前，閘室消能防沖刷采用的技術措施主要包括底流消能、挑流消能以及面流消能。就具體水利工程而言，適用的技術方案選擇主要取決于過閘的水流流量、上下游水位差、水利工程工況條件以及工程地質狀況等因素。閘室消能防沖刷的計算是水閘工程設計中的一項重要環節，其中消力池、海漫段及防沖刷槽等消能防沖刷設施的結構尺寸都必須經過嚴格計算才能確定。

為實現對杜家臺分洪口閘室消能防沖刷的控制，需要從計算水位組合的選取、河道不沖流速以及允許流量的計算、閘門開度及過閘流量的確定、消能防沖刷計算過程、消能防沖刷設施布置等方面進行控制。通過引進紊流模型，利用紊流模型中的雷諾時均法，建立控制目標函數[3]，對水閘工程消能防沖刷設施進行設計。在此過程中，根據杜家臺分洪口閘室運行工況，建立上游水體連續方程（1）與動量方程（2），方程表達式如下：

式中 ν——紊流黏性系數；

j——張量。

上述計算公式中，i 的取值為正整數，包括1，2，3等。在上述設計內容的基礎上，建立杜家臺分洪口閘室水面模擬模型，模型需要參照VOF 模型構建[4]。在此模型中，假設水流中單相流呈現相同狀態，此時造成水流出現流速差的主要原因是流體密度與分子黏性系數差異。因此，可以根據流體的體積分數，通過對流體體積的加權平均計算，得到流體密度的體積分數函數[5]。函數表達式如下：

式中 α——水體的體積分數；

ρ1——水體密度；

ρ2——空氣密度。

在上述設計內容的基礎上，將流體流速作為杜家臺分洪口閘室消能防沖刷控制目標，根據閘室中間段的支孔數量及其空間幾何尺寸參數，建立閘室三維結構模型[6]。同時，使用六面圖與五面體網格劃分成三維結構模型，分為若干個網格區域，再將杜家臺分洪口閘室消能防沖刷控制工作作為單元控制目標，建立如下計算公式所示的目標函數。函數表達式如下：

式中 f（ν）——杜家臺分洪口閘室消能防沖刷流速目標；

C——微分方程；

μ——有限容積法；

k——單元控制目標；

ε——網格區域。

按照上述方法，完成杜家臺分洪口閘室消能防沖刷控制目標函數的設計。

2 設計消力池控制條件

在上述設計內容的基礎上，將杜家臺分洪口閘室消能防沖刷消力池作為目標函數的約束條件，通過設計控制條件以明確函數的約束方式。根據杜家臺分洪口閘室中閘孔寬度，參照閘孔凈孔寬度與閘室的分布，掌握閘室中總孔數[7]。判定并識別杜家臺分洪口閘室的下游段是否需要在運行中布置消能裝置和并確定消能裝置的布置數量。如果明確分洪口閘室需要布置消力池作為消能裝置，可以根據消能控制情況設計閘道的交合水位[8]。同時，將設計水位作為分洪口閘室閘門全開時泄洪水位，當出現淹沒出流情況時，定義此時為無需消能情況?？紤]到杜家臺分洪口閘室在正常使用條件下，上游水位為河道的最高蓄水位，只有閘門出現局部開啟行為時才會出現閘門水流流速過大的問題，從而造成較為嚴重的河床沖刷現象。將出現此種現象作為前提條件，在杜家臺分洪口閘室對應布置消能設施，設計消能設施的約束條件，約束條件表達式如下：

式中 h1——分洪口閘室水流躍后水深（m）；

h2——分洪口閘室水流收縮水深（m）；

E0——杜家臺分洪口閘室上游水流總體勢能（J）；

q——水流過閘時，單閘最大流量[m3/（s·m）]；

φ——水流流速控制系數，通常情況下，φ 的取值為常數，此次計算中，取值φ 為0.85；

g——重力加速度（m/s2）。

按照上述方式，完成對杜家臺分洪口閘室消能防沖刷控制中消力池控制條件的設計。

3 分洪口閘室消能防沖刷控制量計算

完成上述設計后，在操作終端創建控制驅動程序，通過對分洪口閘室消能防沖刷控制量的計算，實現對消能防沖刷的控制。在操作終端界面中創建一個窗口，在此窗口中創建并設置4 個驅動程序指令，其中程序1 為參數輸入程序；程序2 為分洪口閘室孔口設計程序；程序3 為消能控制實時情況反饋程序；程序4 為閘口收縮斷面水面試算程序。在此過程中，在程序1 中輸入現有與分洪口閘室相關的參數條件，輸入參數條件應該包括：設計過閘流量、底板高程、單孔閘寬、流速系數、水流動能校正系數、消力池底寬數、邊坡系數以及重力加速度系數。參照《水閘設計規范》的標準要求，明確水流動能校正系數取值為正數，設定此系數的取值在1.05～1.10 之間，根據程序3 反饋的數據計算出泄洪流的控制量，計算公式如下：

式中 Q——泄洪流控制量；

χ——開啟度；

γ′——閘門垂直收縮系數；

b——杜家臺分洪口閘室中閘門寬度；

H0——閘室可開啟孔數。

根據計算的控制量，驅動程序4 控制分洪口閘室消能防沖刷行為，一旦上游出現泄洪行為，可以根據程序3 反饋的實時水位數據，立即驅動控制行為。按照上述方式，完成控制方法的設計，杜家臺分洪口閘室消能防沖刷計算成果見表1。

表1 杜家臺分洪口閘室消能防沖刷計算結果

4 實例應用分析

上文完成了分洪口閘室消能防沖刷控制方法的設計，為實現對此方法在實際應用中效果的校驗，下面以杜家臺分洪口閘室為例，通過實例應用的方式，檢驗本文設計控制方法的有效性。為確?？刂乒ぷ鞯囊幏痘瘜嵤?，在開展相關研究前，應先采用現場勘察的方式，搜集并匯總杜家臺分洪口閘室基本運行情況。本次研究的閘室為開敞式泄洪閘，根據其建設工程信息，統計泄洪閘概況。相關內容見表2。

表2 杜家臺分洪口閘室運行基本概況

掌握杜家臺分洪口閘室的基本情況后，使用本文設計方法啟動控制閘室消能防沖刷行為?？刂七^程中，根據閘室設計運行條件設定分洪口閘室消能防沖刷控制的目標函數，在此基礎上設計消力池控制條件，通過對分洪口閘室消能防沖刷控制量的計算，實現對分洪口閘室消能防沖刷的控制。完成控制后，計算杜家臺分洪口閘室消能防沖刷各個斷面的流速，將其與下游閘道水流流速進行對比，檢驗杜家臺分洪口閘室消能防沖刷控制效果。見圖1。

圖1 杜家臺分洪口閘室消能防沖刷水體流速控制效果

圖1 中：（1）表示下游水體流速；（2）～（4）表示杜家臺分洪口閘室消能防沖刷各個斷面水體的流速。從圖1 所示的結果可以看出，下游水體流速穩定在（1～3）m/s 范圍內，上游水體經過消能防沖刷控制后，流速明顯快速下降，并在12 s 后穩定在3 m/s 范圍內，即設計的方法實現了對上游水流流速的有效控制。在此基礎上，按照下述公式，計算消能控制后的水流動能：

式中 e——消能控制后水流動能；

V——水流流速；

g——重力加速度，計算中g 的取值為常數。

最后統計3 個斷面的水流動能變化情況，其結果見圖2。

圖2 斷面動能變化趨勢圖

5 結 語

文章通過設計杜家臺分洪口閘室消能防沖刷控制目標函數、設計消力池控制條件、分洪口閘室消能防沖刷控制量計算，完成分洪口閘室消能防沖刷控制方法的設計。該方法在通過實踐檢驗后證明具有實際應用價值，利用此方法可以實現對上游水流流速與攜帶動能的有效控制。