魯 洪

（貴州省水利水電勘測設(shè)計研究院有限公司，貴州 貴陽 550002）

1 概述

沉沙池乃是水利樞紐工程中重要排沙、沖沙的水利設(shè)施[1-2］，其運營可靠性對水利樞紐水沙演變、泥沙淤積等均有影響，而沉砂池內(nèi)由于水沙演變關(guān)系，池內(nèi)滲流場特征受結(jié)構(gòu)設(shè)計影響顯著[3-4］，因而研究沉砂池的設(shè)計方案，對推動沉砂池流場穩(wěn)定性與沉沙運營可靠性兼具價值。呂會嬌等[5］、職承杰等[6］基于水工模型復(fù)制比例理論，設(shè)計了沉砂池、消能池等水工結(jié)構(gòu)模型，開展了相應(yīng)的滲流試驗，分析了池內(nèi)滲流場特征，評價水工結(jié)構(gòu)設(shè)計參數(shù)下體型優(yōu)化。高興輝[7］、劉許超等[8］采用數(shù)值計算方法，對沉沙池、閘門等水工結(jié)構(gòu)的靜、動力特性開展計算分析，探究不同運營工況荷載下水工結(jié)構(gòu)應(yīng)力、位移變化特征，豐富了水工結(jié)構(gòu)設(shè)計優(yōu)化研究成果。對于水工結(jié)構(gòu)三維流場計算問題，范海東等[9］、郭鑫等[10］采用Flow 3D、Fluent 等計算手段開展了消能池、溢洪道及沉沙池滲流特征分析，探討了三維流場內(nèi)流速、水面線等特征參數(shù)影響特性。本文為研究桂北水利樞紐沉沙池結(jié)構(gòu)翼片設(shè)計優(yōu)化問題，采用ANSYS-Fluent 開展沉沙池斜板翼片設(shè)計改進，從水力特征角度遴選最優(yōu)方案。

2 工程建模分析

2.1 工程概況

北江上游桂北地區(qū)水資源豐富，計劃分期建設(shè)一水電站樞紐工程，一期規(guī)劃投入運營工程包括有發(fā)電廠房、防洪主壩及沉沙池與沖沙閘等排沙設(shè)施。由于上游北江、珠江等支流流域內(nèi)人口活動，溢洪道、發(fā)電引水口等區(qū)段泥沙懸浮、沉降較多，水位受泥沙影響，削弱水利工程運營效率，因而工程管理部門計劃在一期工程中建設(shè)沖沙閘與沉沙池。沉沙池設(shè)計剖面見圖1，包括有引水渠、擴散段及池內(nèi)沉沙段及下游泄流通道，沖沙閘建設(shè)在沉沙池尾端。該沉沙池整體高程分布為345.822 m～346.968 m，池內(nèi)坡度為1／100，全長為12 m，擴散段寬度為7 m。在確保沉沙主池運營安全可靠的前提下，計劃對該沉沙池進行設(shè)計改進，提升其沉沙效率與泄流期安全綜合穩(wěn)定性，采用斜板沉砂池作為結(jié)構(gòu)設(shè)計體型，池內(nèi)配置有斜板構(gòu)件，覆蓋有翼片等擋流結(jié)構(gòu)，設(shè)計長度為0.3 m～0.45 m，斜板與翼片為平行相貼合布置，水平夾角為60°，結(jié)構(gòu)設(shè)計概化見圖2。為確保該斜板沉砂池改進后不影響池內(nèi)滲流場穩(wěn)定性，考慮以沉砂池斜板翼片高度H 為研究對象，探究其對沉沙池水力特性影響。

圖1 沉砂池總體設(shè)計剖面圖

圖2 沉砂池結(jié)構(gòu)設(shè)計

2.2 設(shè)計模擬

為分析斜板式沉沙池翼片高度參數(shù)H 對池內(nèi)滲流場及水力特性影響，對池內(nèi)沉沙模型、泥沙懸浮運動狀態(tài)進行概化分析，圖3 為池內(nèi)泥沙顆粒懸浮遷移運動演變概化示意。斜板式沉沙池運營期內(nèi)上游引水明渠中水流進入擴散段后削弱水動力勢能，沉積一部分泥沙，而在擴散段其泥沙顆粒沿著斜板表面運動，顆粒碰撞的過程中沉積至斜板與翼片下方。沉沙池進、出水段分別設(shè)置為弧形觸水面，與水流方向分別呈95°、65°。

圖3 泥沙沉降概化圖

采用ANSYS 完成該斜板式沉沙池結(jié)構(gòu)幾何建模，見圖4（a），由于斜板構(gòu)件為沉沙池重要組成部分，因而獨立提取其幾何模型，見圖4（b）。幾何模型簡化了進、出水段，且省略尾端沖沙閘結(jié)構(gòu)，重點分析流體運動下斜板翼片設(shè)計對池內(nèi)滲流場影響。經(jīng)ANSYS 幾何模型完善后，導(dǎo)入至CFD 計算Fluent 流場軟件中[4，9］，采用不規(guī)則四邊體為沉沙池網(wǎng)格劃分單元體，劃分后三維流場計算模型見圖5。該模型共有88274個微單元，75486 個節(jié)點，網(wǎng)格單元最小體積為6.3528×10-9m3，在圖4（b）所示的斜板構(gòu)造處進行加密劃分，該部分區(qū)域網(wǎng)格覆蓋密度較之其他非加密區(qū)域增大了33.5%。模型中頂、底面分別氣液二相交流邊界與液相獨立邊界。本文模型中X～Z正向分別為出水渠方向、斜板傾斜向及水體豎直向下。

圖4 斜板式沉砂池模型

圖5 沉砂池三維計算模型

改進后沉沙池鋪設(shè)有間距為1 m 的翼片，其高度按照不超過池寬度的3／4 設(shè)定，池深寬度為8 m，因而本文設(shè)定翼片高度H 計算方案為1 m～6 m，各方案間階次為1 m，且設(shè)定有無翼片斜板式沉沙池方案，以翼片高度為0 m 指代。計算工況中初始泥沙含量輸入為1.5 kg／m3，流速為1.2 m／s?；贔luent 計算不同設(shè)計方案下斜板式沉砂池滲流場特征，探討斜板上翼片設(shè)計參數(shù)的影響變化特性。

3 沉沙池水力特征影響

3.1 流速

由于斜板式沉沙池翼片高度設(shè)計參數(shù)差異，池內(nèi)流速會有相應(yīng)影響，圖6 為翼片高度不同方案下池內(nèi)流速等值線與斷面變化特征。

從圖6（a）中可知，不論翼片高度為何值，流速等值線走向與水流方向均一致，且均以水面線流速最高，即翼片高度設(shè)計參數(shù)對流速等值分布影響較小。從池內(nèi)流速等值線分布來看，翼片高度5 m、6 m 方案中池內(nèi)中、下游等區(qū)段聚集有渦旋流等現(xiàn)象，高、低流速在該區(qū)域內(nèi)發(fā)生水力動勢能的“碰撞”，極易引發(fā)池內(nèi)水流的翻滾、擾動及漩渦等現(xiàn)象，不利于水體沖砂、排沙[11］，故從沉沙池運營考慮，翼片高度不應(yīng)過大。

分析圖6（b）可知，斜板增設(shè)翼片方案下，其流速水平均高于無翼片方案，翼片高度1 m、4 m 方案下流速水平較之無翼片方案分別增長了38%～68.9%、1.75～3.3 倍，即增設(shè)翼片高度對沉沙池流速提升具有顯著效果。當(dāng)改變翼片高度時，池內(nèi)沿程斷面流速具有差異性，且翼片高度與池內(nèi)流速水平為正相關(guān)關(guān)系，在池內(nèi)斷面2.4 m 處翼片高度1 m 方案下流速值為0.03 m／s，而高度3 m、4 m、6 m 方案下流速較之分別增大了35.7%、65.4%、85.6%，表明增大翼片高度，沉沙池內(nèi)流速可增大，對水體排沙、沖沙具有正面作用。在翼片高度1 m 方案下沿程斷面平均流速為0.025 m／s，而隨高度方案階次遞增1 m 時，沿程平均流速的增幅達18.4%，但增幅最大集中在翼片高度1 m～4 m 區(qū)間，該區(qū)間內(nèi)平均流速的最大增幅為26.8%，平均增幅為22.3%，而在翼片高度超過4 m 方案平均增幅僅為9.8%，表明翼片高度設(shè)計參數(shù)對沉沙池斷面流速促進作用為減弱態(tài)勢。對比沿程流速變化可知，當(dāng)翼片高度超過4 m 后，受限于局部漩渦流影響，其峰值流速斷面具有差異性，高度5 m、6 m 方案中峰值流速分別為0.066 m／s、0.074 m／s，位于斷面8.4 m 處，且在該斷面周圍具有流速波動性。高度1 m～4 m 方案內(nèi)峰值流速分布為0.033 m／s～0.054 m／s，均為斷面3.6 m，沿程流速穩(wěn)定性均較佳，從“上-中-下游”沿程斷面，流速為先增后減變化。綜合分析可知，翼片高度超過4 m 后，流速穩(wěn)定性欠佳，且漩渦等非穩(wěn)定滲流易引起流速波動段。

圖6 翼片高度影響下流速變化（圖（a）從下至上分別為高度0cm、2m、4m、5m、6m）

3.2 壓強

圖7 為斜板翼片高度設(shè)計參數(shù)影響下池內(nèi)時均壓強變化特征。分析壓強變化可知，壓強水平最高為翼片高度4 m 方案，隨翼片高度增大，池內(nèi)壓強為先增后減變化；對比沿程平均壓強也可知，無翼片方案中沿程平均壓強為6.6 kPa，而增設(shè)翼片后，高度1 m 方案下平均壓強為9.92 kPa，隨高度方案階次每增長1 m 下，在高度1 m～4 m 區(qū)間內(nèi)沿程平均壓強的增幅為28.1%，而超過高度4 m 區(qū)間內(nèi)平均壓強具有降幅17.3%。從壓強量值對比來看，應(yīng)控制翼片高度在低于4 m，更有利于池內(nèi)氣液固三相分布，限制泥沙顆粒運動，降低水體泥沙懸浮量[10］。

圖7 沉沙池沿程斷面壓強特征

池內(nèi)沿程斷面上壓強均為遞增，以出渠口處壓強值最高，在高度1 m～4 m 方案內(nèi)沿程壓強增幅均較穩(wěn)定，各方案斷面間平均增幅分布為5.7%～8.2%，以高度4 m 方案下壓強受影響敏感度最高。而在高度超過4m 方案內(nèi)， 壓強在斷面0～7.2 m 內(nèi)具有遞增特性，而在臨近出渠口斷面處壓強呈穩(wěn)態(tài)分布，壓強遞增效應(yīng)受限制，甚至高度6 m 方案在斷面9.6 m 后高度6 m 方案的壓強水平低于高度3 m 方案。綜合壓強特征可知，翼片高度4 m 下，沉沙池內(nèi)水力特征技術(shù)優(yōu)勢最大，對沉沙、排沙運營具備更高的可靠性。

4 沉沙池水沙特征分析

圖8 為本文模擬計算運營期500 d 時池內(nèi)含沙量變化特征。

圖8 沉沙池運營期500d 含沙量變化特征

分析圖中含沙量變化可知，翼片高度愈大，則含沙量水平愈低，翼片高度1 m 方案下斷面平均含沙量為3.47 kg／m3，而隨翼片高度每米遞增，其平均含沙量的降幅為13.1%，由此表明，翼片高度愈大，對池內(nèi)排沙、沉沙效率具有促進作用。另一方面，池內(nèi)斷面含沙量均為遞增特征，以出渠口含沙量水平最大，翼片高度1 m 方案中各斷面間含沙量的增幅為4.8%，在翼片高度3 m、4 m 方案中相應(yīng)的增幅分別為4.2%、3.2%，即翼片高度愈大，斷面含沙量水平控制效果愈佳，但是在翼片高度5 m、6 m 方案中含沙量的增幅分別“慢-快”兩階段，兩個方案均在斷面4.8 m 后出現(xiàn)含沙量的快增變化，該斷面后含沙量的平均增幅分別可達7.2%、11.6%，即翼片高度愈大，由于對水體控制性欠弱，因而含沙量水平在一定斷面區(qū)間上會“脫離”斜板沉沙、排沙控制[12］，導(dǎo)致含沙量出現(xiàn)高增幅階段。綜合池內(nèi)水力特性及水沙演變，認為翼片高度4 m 時沉沙池結(jié)構(gòu)設(shè)計水平最優(yōu)。

5 結(jié)論

（1）翼片高度對池內(nèi)流速等值線分布影響較小，但高度5 m、6 m 方案內(nèi)流速等值線分布存在渦旋，且在峰值流速斷面周圍區(qū)域存在流速波動段；翼片高度愈大，流速愈高，但增幅集中在高度1 m～4 m 方案，超過高度4 m 后平均增幅僅為9.8%。

（2）各翼片方案中以高度4 m 下壓強水平最高，在翼片高度1 m～4 m 與4 m～6 m 梯次方案內(nèi)，分別具有平均增幅28.1%與降幅17.3%；蓋度4 m 方案下斷面間壓強增幅最高，而高度5 m、6 m 方案臨近出渠口壓強增長停滯。

（3）翼片高度與含沙量水平為負相關(guān)，翼片高度方案每梯次下，池內(nèi)斷面含沙量的降幅為13.1%；池內(nèi)斷面間含沙量為遞增變化，但高度5 m、6 m 方案在斷面4.8 m 前、后的遞增呈“慢-快”兩階段特征。

（4）綜合流場計算，沉沙池翼片高度4 m 方案下排沙、沖沙設(shè)計最優(yōu)。