張桂花,胡士輝,張曉雷,李 彬
(1.河南省水利勘測設計研究有限公司,河南 鄭州450016;2.黃河水利委員會 水文局,河南 鄭州450004;3.華北水利水電大學 水利學院,河南 鄭州450046;4.中水北方勘測設計研究有限責任公司,天津300222)
近年來我國碾壓混凝土重力壩建設較多,工程所處區域地形、地貌、地質條件越來越復雜,在壩高、泄流量、河段地形、地貌、地質條件、經濟性等因素影響下無法采用底流消能方式時,一般采用挑流方式消能。挑流消能工挑坎體形設計甚為關鍵,一些學者對其進行了大量研究,如:麥麥提江·艾麥提[1]采用FLUENT軟件對挑流水舌形狀進行了數值模擬,并對連續式挑流鼻坎進行了體形優化;衛勇[2]采用整體水工模式對復式挑流鼻坎的應用進行了研究,得到復式挑流鼻坎對減小岸邊沖刷深度和降低下游淤積高程具有良好效果的結論;李澤等[3]采用整體水工模型研究提出了分條復式扭曲挑流鼻坎;李元杰等[4]采用整體水工模型研究提出了一種新型非對稱轉向收縮差動式挑坎;劉斌[5]通過水力學模型試驗研究提出兩底孔出口可采用矩形不對稱寬尾墩形成的窄縫體形,三表孔可采用左、右表孔差動式齒坎與中表孔舌形鼻坎組合的優化方案來改善水墊塘內的水流流態,提高消能效率,使水墊塘底板所受動水沖擊壓力及脈動壓力達到設計要求。
大石澗水庫工程由于大壩上下游水位落差大,溢洪道單寬流量大,消能段河谷狹窄,加之地形、地貌、地質等因素影響,因此其泄洪消能工設計受到較大的限制,而上述學者所采用的挑流消能工無法解決該工程的泄水消能沖刷問題。筆者結合連續式和差動式挑流鼻坎體形的優點,研究提出了一種新型舌瓣鼻坎消能工,并通過物理模型試驗和數值模擬等對其進行設計優化,最終使該消能工在實際工程中得到成功應用。
大石澗水庫工程位于黃河二級支流渡洋河中上游,是一座以灌溉供水為主并可起到攔洪滯蓄作用的中型水利樞紐工程。其控制流域面積173 km2,總庫容3 309萬m3,大壩為碾壓混凝土重力壩,建筑物級別為3級,設計壩高68.9 m、壩頂長度295.0 m,大壩由擋水壩段和溢流壩段組成,共分14個壩段,其中溢流壩段64.0 m、左岸擋水壩段122.5 m、右岸擋水壩段108.5 m。表1給出了水庫設計及校核洪峰流量。

表1 水庫設計及校核洪峰流量
大石澗水庫壩址區河槽兩岸均為雄厚的中低山山體,山頂高程600~710 m,兩岸岸坡坡度35°~50°,局部近直立,泄洪消能區左岸有一凸出的山包,泄洪消能設施的布置受到較大限制。水庫溢流壩下游河床覆蓋層厚8~10 m,局部達12 m。原設計泄洪消能方案采用連續坎挑流消能方式,其模型試驗結果表明:庫區來流平順,泄流建筑物泄流能力滿足要求,表孔堰面曲線設計合理,溢流表孔出流正對河道,水流沒有沖刷兩岸;在設計工況下,下游河道最大沖刷深度12.8 m,略深于覆蓋層厚度;由于水舌空中擴散程度低,水流落點過于集中,因此兩岸均有不同程度的淘刷,影響岸坡穩定。為達到最佳經濟效果,不宜處理地形地質條件,因此只能對末端挑流鼻坎進行修改和優化[7],考慮對原挑流消能設計進行優化,解決其消能防沖問題。
常用挑流鼻坎主要有連續式和差動式,兩種挑流鼻坎體形各有優缺點[2]:連續式挑流鼻坎構造簡單,不易發生空蝕破壞,但挑流時水股密實集中,高水頭大流量時對下游河床沖刷較深;高低差動式挑流鼻坎的沖淤情況有較大改善,上下兩層水舌挑距差太小,水舌落點仍較為集中,消能效果較連續式好,但易發生空蝕破壞。
基于有關學者已有研究成果[1-6],結合連續式和差動式挑流鼻坎體形的優點,研究提出了一種新型舌瓣鼻坎體形,即在連續式基礎上向前不連續伸出一段舌形底板形成長短差動的舌形挑流鼻坎。
大石澗水庫工程溢流壩總寬度為60.0 m,其新型舌瓣鼻坎消能工采用三瓣式舌形挑流鼻坎,即整個消能工由3個大舌形鼻坎加2個缺口組成,見圖1。經過多種方法多次對比調整,最終確定左右兩側舌形鼻坎對稱布置,均由半徑13.7 m和10.01 m的兩段圓弧組成,鼻坎反弧半徑均為15 m,最大挑射角35.0°;中間舌形鼻坎平面投影為半徑11.79 m的一段圓弧,最大挑射角30.0°;缺口寬度為3 m,挑射角15.0°。

圖1 新型舌瓣鼻坎消能工布置(單位:m)
根據溢流堰堰上水頭及泄流能力,物理模型選取模型比尺為1∶40。物理模型按照重力相似準則設計,幾何比尺λL=40、流速比尺λv=6.32、壓力比尺λP=40、流量比尺
動床模擬:工程沖刷區附近河中覆蓋層厚度約8 m;左岸邊坡較緩,覆蓋層較厚,為10~12 m;右岸邊坡較陡,覆蓋層厚僅2 m左右。為了提高物理模型沖刷與原型的相似性,物理模型將消能區左岸動定床邊界適當外移,擴大動床范圍,將左岸動床模型沙按基巖和覆蓋層分層模擬,其中下層基巖用中值粒徑3 cm的鵝卵石模擬,其抗沖流速折合原型流速為5 m/s左右。
圖2 給出了整體模型試驗平面布置情況。
3.2.1 挑流水流流態
圖3 給出了不同泄洪工況下挑流水舌形態及下游河道水流流態。不同泄洪工況下挑流水舌水流的主入水區均基本位于河床中心位置,避免了挑流水舌對兩岸坡的直接沖刷,同時減輕了對河床的沖刷;校核洪水時,中間水舌最遠落點較左右兩側水舌最遠落點近2~3 m,缺口和左右邊墻處挑流水舌落點距壩腳50 m左右,消能區左右岸最大回流流速為5.53、3.95 m/s,下游河道最大流速為5.87 m/s,平均流速為5.55 m/s;設計洪水時,水舌擴散形態和校核洪水相似,范圍略小,但水舌更薄、更光滑,缺口和左右邊墻處挑流水舌落點較校核洪水時短10 m左右,消能區左右岸最大回流流速分別為4.31、3.96 m/s,下游河道最大流速6.54 m/s,平均流速5.04 m/s;30 a一遇洪水時,水舌在空中呈燈泡形,水舌光滑擴散形態較好,挑流水舌最遠落點處三股水舌落點位置相差不大,消能區左右岸最大回流流速為3.85、1.47 m/s,下游河道最大流速5.44 m/s,平均流速4.89 m/s。

圖2 大石澗水庫工程整體模型試驗平面布置示意

圖3 新型舌瓣鼻坎消能工物理模型試驗挑流水舌形態及下游河道水流流態
3.2.2 下游河道沖淤形態
圖4 給出了不同泄洪工況下下游河道沖淤形態。校核工況消能區兩岸沖刷嚴重,最大沖深6.5 m,仍比壩基礎高出6.5 m,不會影響壩體穩定;設計工況河中沖深比校核工況減少1.5 m,左岸沖刷范圍減小,基巖略有沖刷,右岸基巖沖刷6 m左右;30 a一遇洪水時,河中沖深與設計工況接近,但范圍較小,消能區左岸覆蓋層淘刷后基巖裸露,但基巖沒有產生沖刷,右岸基巖沖刷3~4 m,不會影響岸坡穩定。
綜上分析,新型舌瓣鼻坎消能工起到了明顯的消能效果,減小了下游沖刷范圍及沖坑深度。

圖4 新型舌瓣鼻坎消能工物理模型試驗下游河床沖淤形態
3.2.3 挑坎段流速及壓力分布
為了分析新型挑流鼻坎體形設計的合理性,沿模型鼻坎布置了9個測壓點,各測點流速值見表2;在左側鼻坎上布置了4個測壓點來測試鼻坎段壓力,各測點位置、高程及校核、設計和30 a一遇洪水3種工況下壓力值見表3。
不同泄洪工況溢流表孔泄流時堰面流速呈逐漸增大趨勢,進口流速最小,鼻坎反弧最低點流速達到最大,出口流速略有減小,挑坎段流速分布較為均勻,水流流態較好;校核洪水時,墩頭進流表底流速分別為5.6、8.3 m/s,鼻坎上最大流速達27.29 m/s;設計、30 a一遇洪水時,鼻坎最大流速分別為26.04 m/s和25.80 m/s。

表2 溢流表孔沿程流速分布 m/s

表3 溢流表孔鼻坎段壓力及水流空化數
由表3可知:校核、設計和30 a一遇洪水等工況新型舌瓣鼻坎壓力均較大,不同泄洪工況挑坎段各測點均未出現負壓情況;最大壓力出現在校核工況時鼻坎反弧低點的16號測點,數值為7.93×9.8 kPa;鼻坎段最小空化數0.36,最大水流流速均未超過30 m/s。由于不同泄洪工況挑坎段沿程均未出現負壓,因此挑坎結構不易發生空蝕破壞。
3.2.4 消能效果分析
通過水工試驗對溢流表孔進行了連續挑坎、高低差動挑坎、長短差動挑坎、舌瓣坎等多種方案對比,表4給出了各方案沖刷結果。

表4 各種挑坎沖刷結果 m
由表4可知,6種挑流鼻坎方案河中沖刷深度相差不大,但從左右岸和護坦末段沖刷情況來看,方案6優于其他方案。另外,方案6水舌形態好,左岸回流范圍小、回流強度不大,水舌空中擴散程度高,在入水落點處展開,有效減輕了對河道的沖刷,岸邊沖刷范圍較小。根據下游沖淤結果及水流流態情況,最終采用新型舌瓣挑坎方案。
采用Flow-3D軟件中的泥沙沖刷模型結合RNG模型和FAVOR方法,進行新型舌瓣鼻坎消能工數值模擬試驗,Flow-3D軟件的數學控制方程及其解法詳見文獻[8]。按原型1∶1的比例搭建數學模型,數學模型采用結構化正交網格,其中:一般區域采用網格尺寸為2 m×2 m×2 m、堰體及消能區局部加密網格尺寸為1 m×1 m×1 m,數學模型共有網格600萬個;模型定床糙率設定為0.03,動床均值粒徑取0.02 m;上邊界采用壓力邊界并設定相對應水位,下邊界采用自由出流(Outflow)邊界,其他邊界均采用固壁(Wall)邊界條件。圖5給出了新型舌瓣鼻坎模型效果圖。

圖5 新型舌瓣鼻坎模型效果
圖6 給出了30 a一遇洪水的流速分布、沖淤地形及最大沖坑深度。

圖6 新型舌瓣鼻坎模型挑流數值模擬結果
由圖6可知:溢流堰頂的流速為3.44 m/s,沿著溢流堰下泄過程中水流受到重力、堰面摩擦力和空氣阻力等多種因素影響,流速沿程逐漸增大,在挑坎處達到最大值27.05 m/s,與物理模型流速最大值25.80 m/s相差4.8%;數學模型最遠挑距距壩腳85 m,下游最大沖深12.0 m,與物理模型最大沖深12.5 m相差-4.0%。
為了進一步驗證數學模型的可靠性,將表4中各方案30 a一遇洪水的沖坑試驗結果與數值模擬值進行對比,具體見表5。
由表5可知:數學模型的沖刷模擬結果與試驗結果基本一致,表明Flow-3D軟件可應用于挑流消能沖刷的數值模擬,新型舌瓣鼻坎消能工可以有效減小下游沖刷范圍及沖坑深度,消能防沖效果明顯。
相比于方案1連續坎,方案2、3護坦末端侵蝕嚴重,方案4河中沖刷深點距護坦過近(30 a一遇工況鼻坎缺口挑流落點試驗值距護坦22 m),方案5優于方案1,方案6最佳。

表5 各鼻坎體形30 a一遇沖刷特征值對比
(1)新型舌瓣鼻坎消能工整個形狀同3個獨立擴散的燈泡形水舌結構,缺口處挑流水舌短且較為集中,水流落點遠離護坦且能與舌瓣鼻坎水舌落點錯開,水舌平面擴散范圍大,水舌左右兩側落點前后均勻拉開,可使水舌在空中充分摻氣碰撞消殺能量,從而減小下游河道沖刷范圍及沖坑深度。
(2)不同泄洪工況下挑坎段流速分布較為均勻,水流流態較好,挑坎段沿程未出現負壓,挑坎結構不易發生空蝕破壞。
(3)大壩上下游水位落差大,單寬流量大,消能段河谷狹窄,地形、地貌、地質條件差,泄洪消能工設計受到較大限制時,可采用新型舌瓣鼻坎消能工,其具體設計參數可采用物理模型和Flow-3D軟件搭建數學模型研究確定。