泵站直邊正向前池流態(tài)模擬與泥沙淤積預(yù)防措施

徐存東,劉璐瑤,王國(guó)霞,田子荀,劉麗霞,程 慧

(1.華北水利水電大學(xué)水資源高效利用與保障工程河南省協(xié)同創(chuàng)新中心,河南 鄭州 450045;2.甘肅農(nóng)業(yè)大學(xué)工學(xué)院,甘肅 蘭州 730000)

泵站直邊正向前池流態(tài)模擬與泥沙淤積預(yù)防措施

徐存東1,劉璐瑤1,王國(guó)霞1,田子荀1,劉麗霞2,程 慧1

(1.華北水利水電大學(xué)水資源高效利用與保障工程河南省協(xié)同創(chuàng)新中心,河南 鄭州 450045;2.甘肅農(nóng)業(yè)大學(xué)工學(xué)院,甘肅 蘭州 730000)

針對(duì)從高含沙河道取水的大型泵站直邊正向前池泥沙淤積導(dǎo)致前池調(diào)節(jié)能力下降、泵站運(yùn)行效率低下等問(wèn)題,利用ICEM-CFD軟件構(gòu)建原型前池結(jié)構(gòu)三維模型,采用FLUENT軟件基于Standardk-ε湍流模型對(duì)前池流態(tài)進(jìn)行模擬分析,并針對(duì)不良水流流態(tài)提出增設(shè)導(dǎo)流墩與壓水板2種優(yōu)化措施。結(jié)果表明,數(shù)值模擬結(jié)果與實(shí)際調(diào)查情況相一致,前池水流流態(tài)紊亂,泥沙淤積較嚴(yán)重,八字形導(dǎo)流墩可在一定程度上削減主回流區(qū)范圍,但對(duì)前池水流僅起到分配作用,流場(chǎng)分布無(wú)明顯改善;45°壓水板能有效改善前池主流效應(yīng),兩側(cè)回流區(qū)明顯減弱,在改善前池泥沙淤積問(wèn)題的同時(shí)保證了水泵有利的進(jìn)水條件。

泵站前池;泥沙淤積;流態(tài)模擬;流態(tài)改善;八字形導(dǎo)流墩;45°壓水板

甘肅省景泰川電力提灌工程灌區(qū)(簡(jiǎn)稱(chēng)景電灌區(qū))是位于我國(guó)西北干旱區(qū)的大型引黃灌區(qū)。該灌區(qū)內(nèi)建有泵站43座,裝機(jī)容量25.97萬(wàn)kW,灌溉面積6.51 萬(wàn)hm2。由于泵站引水含沙量高,加之泵站設(shè)計(jì)不合理造成池內(nèi)流態(tài)不良,而流態(tài)又是決定挾沙能力的重要因素,故目前灌區(qū)內(nèi)泵站前池淤積問(wèn)題嚴(yán)重,影響了泵站的正常運(yùn)行,制約了灌區(qū)工程效益的充分發(fā)揮[1]。因此,必須深入探索泵站前池泥沙淤積的流態(tài)成因并采取適當(dāng)工程整流措施對(duì)前池的水力條件進(jìn)行改善。

目前國(guó)內(nèi)外學(xué)者針對(duì)泵站前池泥沙淤積防治措施開(kāi)展了大量研究,取得了豐富的研究成果。Constantinescu等[2]利用Standardk-ε方程模擬了泵站前池的漩渦情況,模擬獲得的漩渦結(jié)構(gòu)與模型試驗(yàn)結(jié)果相似,而湍流數(shù)學(xué)模型的選擇和邊界條件的處理方式會(huì)對(duì)漩渦出現(xiàn)的位置和強(qiáng)度造成影響。劉新陽(yáng)等[3]針對(duì)田山一級(jí)泵站采用雷諾時(shí)均N-S方程,結(jié)合Standardk-ε湍流模型進(jìn)行數(shù)值模擬,分析并選擇了非連續(xù)底坎、非連續(xù)挑流坎與壓水板3種整流措施相結(jié)合的整流工程措施。竇國(guó)仁[4]基于潮流與波浪基本方程及其相似條件,并結(jié)合泥沙運(yùn)動(dòng)方面的多年研究成果,構(gòu)建了針對(duì)懸沙和底沙模擬的物理模型相似理論。Nakato[5]構(gòu)建了Chicot湖取水泵站前池的1∶24的室內(nèi)試驗(yàn)?zāi)Ｐ?通過(guò)對(duì)比現(xiàn)場(chǎng)調(diào)查資料和模型試驗(yàn)結(jié)果,確定泥沙淤積分布和平均水深、平均流速及水流流態(tài)關(guān)系密切。但是,采用直邊正向前池的結(jié)構(gòu)形式開(kāi)展池內(nèi)流態(tài)模擬和泥沙淤積防治專(zhuān)項(xiàng)技術(shù)研究的還不多見(jiàn)。本文選取甘肅省景電灌區(qū)一期西干三泵站典型直邊正向前池為模擬對(duì)象,采取不同整流方案分別進(jìn)行研究分析,以期為同類(lèi)泵站的前池流態(tài)改善的數(shù)值模擬提供參考依據(jù)。

1 泵 站 概 況

1.1泵站參數(shù)

景電灌區(qū)一期西干三泵站共布置5臺(tái)套機(jī)組,1～3號(hào)機(jī)組流量均為0.85 m3/s,吸水管直徑為800 mm,4號(hào)機(jī)組流量為0.28 m3/s,吸水管直徑為500 mm,5號(hào)機(jī)組流量為1.6 m3/s,吸水管直徑為1 000 mm;泵站吸水管安裝高程為1 656.40 m,前池設(shè)計(jì)水位為1 658.67 m,前池?cái)U(kuò)散角為36°。為了保證前池流態(tài)數(shù)值模擬的真實(shí)性,選取的計(jì)算區(qū)域除了前池部分(水平段和斜坡段),還包括泵站前池前的10 m引水渠、漸變段、閘室段和進(jìn)水池以及吸水管道,構(gòu)成較為完整的流場(chǎng)范圍[6],其結(jié)構(gòu)形式如圖1所示,圖中,y1-1～y1-5、z1-1～z1-3為觀測(cè)斷面。

圖1 直邊正向前池結(jié)構(gòu)示意圖(單位：m)Fig.1 Structural diagram of straight-edge front inflow forebay (units: m)

1.2現(xiàn)場(chǎng)調(diào)查

現(xiàn)場(chǎng)淤積程度調(diào)查結(jié)果如圖2所示,圖中虛線范圍內(nèi)為泥沙淤積范圍。

由圖2可知,正向泵站前池泥沙淤積形態(tài)近似軸對(duì)稱(chēng)分布,兩側(cè)泥沙淤積頂部范圍大小和厚度接近,厚度約2.6 m,淤積坡面平順,形成了新的過(guò)水流道。淤積形態(tài)的尾部呈喇叭口形,兩側(cè)底腳泥沙淤積至進(jìn)水池,影響泵站兩側(cè)水泵的取水。

圖2 直邊正向前池淤積程度示意圖Fig.2 Diagram of sediment deposition in straight-edge front inflow forebay

2 水流挾沙理論與數(shù)學(xué)計(jì)算模型

2.1水流挾沙能力結(jié)構(gòu)公式

舒安平、費(fèi)祥俊基于固液兩相挾沙水流紊動(dòng)能量方程得出水流挾沙能力的結(jié)構(gòu)公式,證明了水流流態(tài)是決定挾沙能力的重要因素[7],其公式為

(1)

式中:S*——挾沙能力;f(μr)——待定函數(shù);μr——挾沙水流的相對(duì)黏度;κ——水流流態(tài)卡門(mén)常數(shù);fm——流動(dòng)阻力;ρm——挾沙水流密度;ρs——泥沙密度;u——流速;ω——泥沙沉速;g——重力加速度;q、N——經(jīng)驗(yàn)參數(shù)；R——水利半徑。

2.2控制方程

泵站前池水流屬于復(fù)雜的三維湍流運(yùn)動(dòng),水流視為不可壓縮流體(密度為常數(shù)),流體的熱交換量較小,故模擬不考慮能量方程[8]。因此,泵站前池的水流連續(xù)性方程、動(dòng)量守恒方程分別為[9]

(2)

(3)

式中:ui、uj——流速矢量;ρ——流體密度;t——時(shí)間;p——包含湍流動(dòng)能的靜壓力;μeff——流體運(yùn)動(dòng)的有效黏滯性系數(shù);gi——i方向的重力分量。

2.3Standardk-ε模型

Standardk-ε模型是基于湍流動(dòng)能耗散率的半經(jīng)驗(yàn)公式,是針對(duì)湍流充分發(fā)展的高雷諾數(shù)流動(dòng)建立的[10],模型的穩(wěn)定性和精確性在科學(xué)界得到了廣泛驗(yàn)證[11]。泵站前池水流流態(tài)屬于高雷諾數(shù)湍流,因此采用Standardk-ε模型進(jìn)行計(jì)算。其湍動(dòng)能方程、湍動(dòng)能耗散率ε方程分別為[9]

(4)

(5)

式中:μ、μt——層流黏滯性系數(shù)和湍流黏滯性系數(shù);Pk——湍流能生成項(xiàng);C1ε、C2ε、Cμ、σk、σε——經(jīng)驗(yàn)常數(shù),分別取值1.44、1.92、0.09、1.0和1.3,即FLUENT中的默認(rèn)設(shè)置。

圖3 直邊正向前池計(jì)算區(qū)域網(wǎng)格劃分Fig.3 Meshing of computing area of straight-edge front inflow forebay

2.4邊界條件及網(wǎng)格劃分

根據(jù)泵站前池進(jìn)口斷面處的平均流速,選用速度進(jìn)口邊界條件。出口邊界水流為完全發(fā)展流動(dòng),將出口邊界條件設(shè)置為outflow邊界條件[12]。固體壁面邊界處理選擇FLUENT軟件常用的壁面函數(shù)法。泵站前池自由表面受外界環(huán)境影響小,浮動(dòng)變化不大,采用剛蓋假定法模擬計(jì)算自由表面,即在FLUENT中自由表面條件設(shè)置為symmetry[13]。

泵站前池邊界條件較為復(fù)雜,池內(nèi)有吸水管及防淤設(shè)施等,因此采用非結(jié)構(gòu)化網(wǎng)格進(jìn)行模型網(wǎng)格劃分[14]。根據(jù)泵站前池的實(shí)際尺寸,利用ICEM-CFD軟件對(duì)其進(jìn)行三維幾何模型構(gòu)建和網(wǎng)格劃分,局部位置進(jìn)行了適當(dāng)加密,調(diào)整網(wǎng)格數(shù)量進(jìn)行試算,結(jié)果誤差在5%以?xún)?nèi),即滿足網(wǎng)格無(wú)關(guān)性檢驗(yàn),采用的計(jì)算網(wǎng)格單元總數(shù)為1 074 457。直邊正向前池三維模型和網(wǎng)格劃分結(jié)果如圖3所示。

表1 觀測(cè)斷面參數(shù)

注:y是觀測(cè)斷面到渠道中心線的距離,順流向左為正,右為負(fù);z是觀測(cè)斷面距水面線以下距離,負(fù)號(hào)表示水面線以下。

3 模 擬 分 析

將ICEM-CFD軟件生成的網(wǎng)格文件導(dǎo)入FLUENT軟件進(jìn)行模擬分析。模擬計(jì)算精度均取10-4,計(jì)算過(guò)程中,連續(xù)性收斂速度較慢,但并非每項(xiàng)殘差都收斂,模擬過(guò)程才收斂,對(duì)于定常流,通過(guò)比較進(jìn)、出口流量以及觀察其他參數(shù)的殘差波動(dòng)性即可判斷連續(xù)性是否收斂。當(dāng)模型進(jìn)、出口流量誤差幾乎為零,其余參數(shù)殘差達(dá)到精度要求且波動(dòng)不大,故認(rèn)為模擬過(guò)程收斂[15]。

確定觀測(cè)斷面,觀察各斷面的流場(chǎng)情況,觀測(cè)斷面位置如圖1所示,各斷面具體參數(shù)見(jiàn)表1。模擬5臺(tái)機(jī)組全部以設(shè)計(jì)流量同時(shí)運(yùn)行,總設(shè)計(jì)流量4.43 m3/s時(shí)泵站前池各斷面的流場(chǎng)分布,結(jié)果如圖4、圖5所示。

圖4 z1-1～z1-3斷面前池流線Fig.4 Pathlines in forebay from section z1-1 to z1-3

圖5 y1-1～y1-5斷面流速分布云圖及前池流線Fig.5 Velocity distribution contour from section y1-1 to y1-5 and pathlines in forebay

圖6 直邊正向前池導(dǎo)流墩布置示意圖Fig.6 Layout diagram of diversion pier in straight-edge front inflow forebay

由圖4和圖5可以看出,水流進(jìn)入直邊正向前池后,平面擴(kuò)散和立面擴(kuò)散效果很不理想,由于平面擴(kuò)散差,使得主流限制在前池中間,兩側(cè)形成的回流區(qū)進(jìn)一步擠壓主流;由于立面擴(kuò)散差,下層水流流速減小,回流現(xiàn)象更加突出?；亓鲄^(qū)水流流速較低,尤其是前池下層水域,當(dāng)水流流速低于其挾沙最低流速時(shí),水流含帶的泥沙便開(kāi)始在回流區(qū)域內(nèi)沉降、淤積,該結(jié)論與現(xiàn)場(chǎng)調(diào)查結(jié)果吻合。

4 防淤措施模擬

根據(jù)GB 50265—2010《泵站設(shè)計(jì)規(guī)范》[16]及直邊正向前池實(shí)際運(yùn)行情況,計(jì)算時(shí)選取以下2種方案分析結(jié)構(gòu)優(yōu)化措施對(duì)其防淤效果的改善:

方案1:采用八字形導(dǎo)流墩進(jìn)行優(yōu)化,八字形導(dǎo)流墩設(shè)在前池斜坡段,采用對(duì)稱(chēng)布置,其主要控制參數(shù)有位置參數(shù)(L1、L2、β)和尺寸參數(shù)(L3)[17],如圖6所示。設(shè)置導(dǎo)流墩長(zhǎng)度為1 000 cm,即L3=1 000 cm,其他各參數(shù)通過(guò)模擬結(jié)果與單一變量最優(yōu)取值選取,見(jiàn)表2。

表2 導(dǎo)流墩位置參數(shù)

通過(guò)數(shù)值模擬發(fā)現(xiàn),八字形導(dǎo)流墩對(duì)泵站前池水流僅起到分流作用,引導(dǎo)部分主流流向前池兩側(cè),對(duì)前池兩側(cè)的回流區(qū)有一定的削減作用,然而在導(dǎo)流墩的作用下,前池水流的平面擴(kuò)散和立面擴(kuò)散效果改善不佳,經(jīng)導(dǎo)流墩分配的原主流變?yōu)?股主流,流場(chǎng)分布不均勻, 水流挾沙能力不足,泥沙在前池淤積。

選取7種參數(shù)組合(表2),分別建立幾何模型、劃分網(wǎng)格,不同組合的網(wǎng)格單元數(shù)在146.4萬(wàn)～151.3萬(wàn)之間,前池下層水流的流速分布矢量圖如圖7所示。

不同參數(shù)組合的導(dǎo)流墩實(shí)施效果差異較大。當(dāng)L1>600 cm時(shí),導(dǎo)流墩距離前池入口過(guò)遠(yuǎn),失去了對(duì)主流的控制作用,主流流道發(fā)生扭曲,如圖7(d)所示。當(dāng)L1<600 cm時(shí),導(dǎo)流墩遠(yuǎn)離前池末端,不能有效改善前池末端回流區(qū)范圍,如圖7(e)所示。當(dāng)L2>300 cm時(shí),導(dǎo)流墩間距過(guò)大,主流區(qū)的回流面積增加。當(dāng)β>12°時(shí),在導(dǎo)流墩背水面極易出現(xiàn)漩渦,且范圍較大,嚴(yán)重壓縮導(dǎo)流墩內(nèi)部主流,如圖7(c)所示。當(dāng)β<12°時(shí),導(dǎo)流墩不能有效削減兩側(cè)回流區(qū),如圖7(a)所示。綜上,參數(shù)組合L1=600 cm、L2=300 cm、L3=1 000 cm、β=12°的八字形導(dǎo)流墩對(duì)直邊正向前池流態(tài)具有一定的改善作用,分配前池水流,導(dǎo)向前池兩側(cè),擠壓和削弱末端回流區(qū),同時(shí)導(dǎo)流墩背水面的漩渦回流現(xiàn)象較弱。

圖7 方案1各組參數(shù)導(dǎo)流墩優(yōu)化前池流線Fig.7 Optimized streamline of the forebay by using different parameters for the diversion pier in solution 1

圖8 直邊正向前池45°壓水板布置示意圖Fig.8 Layout diagram of 45° pressure plate in straight-edge front inflow forebay

方案2:采用45°壓水板對(duì)直邊正向前池進(jìn)行優(yōu)化改造,壓水板主要控制參數(shù)有壓水板下緣入水深度H和到前池入口的距離L,如圖8所示。該方案壓水板參數(shù)取值見(jiàn)表3。

選取9種參數(shù)組合,分別建立幾何模型、劃分網(wǎng)格,不同組合的網(wǎng)格單元數(shù)在117.0萬(wàn)～117.3萬(wàn)之間,各組合下的前池下層水流的流速分布矢量圖如圖9所示。

由圖9可知,泵站前池內(nèi)設(shè)置45°壓水板能夠顯著改善直邊正向前池水流的主流效應(yīng),有效增強(qiáng)前池水流的平面擴(kuò)散和立面擴(kuò)散效果,前池兩側(cè)回流區(qū)明顯減弱,使前池內(nèi)上、下層水流流速分布均勻,上層水流形成穩(wěn)態(tài)擴(kuò)散,改善前池泥沙淤積問(wèn)題的同時(shí)保證了水泵有利的進(jìn)水條件。下層水流流速明顯增加,對(duì)前池底部泥沙具有一定的沖刷作用,同時(shí)水流的挾沙能力也得到大幅度提高,有利于減少前池泥沙淤積。

表3 45°壓水板位置參數(shù)

圖9 方案2各組參數(shù)壓水板優(yōu)化前池流線Fig.9 Forebay optimized streamline by using different parameters for the water pressure plate in solution 2

不同參數(shù)組合的壓水板對(duì)前池水流的作用能力各不相同,通過(guò)對(duì)比模擬結(jié)果可以看出,H和L分別取60 cm和800 cm時(shí),對(duì)前池水流流態(tài)的改善效果最佳。取H<60 cm時(shí),壓水板對(duì)前池水流的作用不夠,導(dǎo)致前池主流的平面擴(kuò)散效果較差,在前池兩側(cè)出現(xiàn)帶狀低速回流區(qū),容易造成泥沙淤積,同時(shí)端機(jī)組附近水流紊亂,水泵進(jìn)水條件較差,如圖9(a)(b)(c)所示。取H>60 cm時(shí),壓水板對(duì)前池水流的作用過(guò)強(qiáng),影響深度過(guò)大,導(dǎo)致下層水流均勻流速區(qū)厚度過(guò)小,甚至低于水泵吸水管口,不利于水泵的取水,造成水泵工作性能下降,工作效率降低,如圖9(g)(h)(i)所示。取L<800 cm時(shí),壓水板的作用未能延伸至前池末端,在端機(jī)組前出現(xiàn)了2個(gè)小范圍低速區(qū),可能會(huì)出現(xiàn)一定的泥沙淤積,如圖9(a)(d)(g)所示。取L>800 cm時(shí),壓水板前空間過(guò)大,水流流速低,在板前出現(xiàn)了漩渦回流,不利于前池的防淤,如圖9(c)(f)(i)所示。

5 改造效果分析

景電灌區(qū)一期西干三泵站于2015年1—5月運(yùn)用八字形導(dǎo)流墩對(duì)前池進(jìn)行改造,并進(jìn)行現(xiàn)場(chǎng)測(cè)試與試驗(yàn),后發(fā)現(xiàn)側(cè)邊壁與導(dǎo)流墩兩側(cè)依然淤積較多泥沙。故次年2—6月改用壓水板對(duì)前池進(jìn)行再次改造,現(xiàn)場(chǎng)測(cè)試與試驗(yàn)表明前池兩側(cè)淤積基本清除,并且水泵效率有了明顯提高。表4給出了前池改造前后淤積程度與泵站的相關(guān)性能參數(shù)。

表4 改造前后數(shù)據(jù)對(duì)比

注：裝置效率=電動(dòng)機(jī)效率×水泵效率×管路效率。

通過(guò)現(xiàn)場(chǎng)測(cè)試發(fā)現(xiàn)壓水板整流后前池內(nèi)流態(tài)得到了較好改善,流量明顯提高,泥沙基本消除。并且因水泵效率為(流量×揚(yáng)程)/軸功率,故由表可見(jiàn)改導(dǎo)流墩與壓水板造后水泵效率分別提高了5.71%與8.68%。

6 結(jié) 論

a. 利用FLUENT軟件對(duì)直邊泵站前池水流流態(tài)進(jìn)行模擬可知,前池兩側(cè)中后部回流區(qū)流速較低,平面和立面擴(kuò)散效果較差,易產(chǎn)生泥沙淤積,與現(xiàn)場(chǎng)調(diào)查情況相符,可見(jiàn)該數(shù)值模擬結(jié)果較為精確,可普遍適用于同類(lèi)型泵站前池優(yōu)化改造措施的模擬研究。

b. 前池內(nèi)設(shè)置八字形導(dǎo)流墩對(duì)前池水流有一定的分配和引導(dǎo)作用,但主流效應(yīng)未得到明顯改善,回流區(qū)范圍和流速變化不大。通過(guò)設(shè)置45°壓水板能夠顯著改善直邊正向前池的主流效應(yīng),提高前池底部水流流速,增強(qiáng)水流平面擴(kuò)散,在改善前池泥沙淤積問(wèn)題的同時(shí)保證了水泵有利的進(jìn)水條件,且取H=60 cm、L=800 cm時(shí),對(duì)前池水流流態(tài)的改善效果最佳,水流挾沙能力最強(qiáng),淤積最輕。

c. 泵站前池內(nèi)設(shè)置參數(shù)合理的八字形導(dǎo)流墩可有效改善前池流場(chǎng)分布,但該措施實(shí)施后達(dá)不到預(yù)期效果,僅對(duì)前池水流起到分配作用,分配后原主流變?yōu)槿芍髁?未能得到相對(duì)較好的流態(tài),且不同參數(shù)組合的導(dǎo)流墩實(shí)施效果差異較大。而45°壓水板對(duì)前池流態(tài)改善效果顯著,泥沙淤積明顯減少,但本次研究對(duì)壓水板參數(shù)選擇較少,且只模擬了導(dǎo)流墩和壓水板的單獨(dú)工作狀態(tài),對(duì)于其他形式防淤措施與組合防淤措施的研究還有待進(jìn)一步深入研究。

[ 1 ] 徐存東.高揚(yáng)程梯級(jí)提水灌溉工程建設(shè)監(jiān)理制的應(yīng)用研究[D].咸陽(yáng):西北農(nóng)林科技大學(xué),2003.

[ 2 ] CONSTANTINESCU G S,PATE V C. Role of turbulence model in prediction of pump-bay vortices[J].Journal of Hydraulic Engineering, 2000,126(5):387-391.

[ 3 ] 劉新陽(yáng),高傳昌,石禮文,等.泵站前池與進(jìn)水池整流數(shù)值模擬[J].排灌機(jī)械工程學(xué)報(bào),2010,28(3):242-246.(LIU Xinyang,GAO Chuanchang,SHI Liwen,et al. Numerical simulation for fluid meliorating in both forebay and suction bay of pumping stations[J].Journal of Drainage and Irrigation Machinery Engineering,2010,28(3):242-246.(in Chinese))

[ 4 ] 竇國(guó)仁.河口海岸全沙模型相似理論[J].水利水運(yùn)工程學(xué)報(bào).2001,1(1):1-11.(DOU Guoren. Similarity theory of total sediment transport modeling for estuarine and coastal regions[J].Hydro-Science Engineering,2001,1(1):1-11.(in Chinese))

[ 5 ] NAKATO T. Application of suction scoops to improve pump-approach flow distributions in three-pump intake bays[C]//Anon.Critical Transitions in Water and Environmental Resources Management.Reston:ASCE,2014:1-10.

[ 6 ] 周大慶,米紫昊,茅媛婷.基于歐拉固液兩相流模型的泵站進(jìn)水側(cè)流場(chǎng)三維模擬[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(1):48-52.(ZHOU Daqing,MI Zihao,MAO Yuanting.3-D numerical simulation of inlet structure flow in pumping station based on Eulerian solid-liquid two-phase flow model[J].Transactions of the Chinese Society of Agricultural Machinery,2013,44(1):48-52.(in Chinese))

[ 7 ] 秦毅,李子文,曹如軒,等.高含沙水流的整體停滯淤積及其影響[J].水利學(xué)報(bào),2011,42(11):1302-1307.(QIN Yi,LI Ziwen,CAO Ruxuan,et al.Whole stagnating boundary siltation of hyperconcentrated flow and its influences[J].Journal of Hydraulic Engineering,2011,42(11):1302-1307. (in Chinese))

[ 8 ] 羅燦,錢(qián)均,劉超,等.非對(duì)稱(chēng)式閘站結(jié)合式泵站前池導(dǎo)流墩整流模擬及試驗(yàn)驗(yàn)證[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(7):100-108.(LUO Can,QIAN Jun,LIU Chao,et al.Numerical simulation and test verification on diversion pier rectifying flow in forebay of pumping station for asymmetric combined sluice-pump station project[J].Transactions of the Chinese Society of Agricultural Engineering,2015,31(7):100-108.(in Chinese))

[ 9 ] 紀(jì)兵兵,陳金瓶.ANSYS ICEM CFD網(wǎng)格劃分技術(shù)實(shí)例詳解[M].北京:中國(guó)水利水電出版社,2012.

[10] 成立,劉超,周濟(jì)人,等.泵站前池底壩整流數(shù)值模擬研究[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2001,29(3):42-45.(CHENG Li,LIU Chao,ZHOU Jiren,et al.Numerical simulation of sill flows in the forebay of pumping station[J].Journal of Hohai University(Natural Sciences),2001,29(3):42-45.(in Chinese))

[11] 徐存東,楊柯,周琳博,等.引黃泵站進(jìn)水池幾何形式的數(shù)值模擬[J].甘肅科學(xué)學(xué)報(bào),2011,23(2):150-154.(XU Cundong,YANG Ke,ZHOU Linbo,et al.Numerical simulation on geometric form of pump sump in the Yellow River pump station[J].Journal of Gansu Sciences,2011,23(2):150-154.(in Chinese))

[12] MAROVIC B,KOLAK D.Combination of CFD dimensions: Industry wide first universal-1D-/3D-computational fluid dynamics simulation software[J].Mechatronik,2012,120(11):34-35.

[13] 曾杭.基于CFD數(shù)值模擬的泵站前池設(shè)計(jì)優(yōu)化[D]. 長(zhǎng)沙:湖南大學(xué),2014.

[14] 唐學(xué)林,王武昌,王福軍,等.泵站前池水沙流的數(shù)值模擬[J].排灌機(jī)械工程學(xué)報(bào),2011,29(5):411-417.(TANG Xuelin,WANG Wuchang,WANG Fujun,et al.Numerical simulation of slit-laden flows in pumping station forebay[J].Journal of Drainage and Irrigation Machinery Engineering,2011,29(5):411-417.(in Chinese))

[15] RAFIEE A,CUMMINS S,RUDMAN M.Comparative study on the accuracy and stability of SPH schemes in simulating energetic free-surface flows[J].European Journal of Mechanics B-Fluids, 2012,36(10):1-16.

[16] 中華人民共和國(guó)水利部.泵站設(shè)計(jì)規(guī)范:GB 50265—2010[S].北京:中國(guó)計(jì)劃出版社,2011.

[17] 資丹,王福軍,姚志峰,等.大型泵站進(jìn)水流場(chǎng)組合式導(dǎo)流墩整流效果分析[J].農(nóng)業(yè)工程學(xué)報(bào), 2015,31(16):71-77.(ZI Dan,WANG Fujun,YAO Zhifeng,et al.Effects analysis on rectifying intake flow field for large scale pumping station with combined diversion piers[J].Transactions of the Chinese Society of Agricultural Engineering,2015,31(16):71-77.(in Chinese))

Flowsimulationoffrontinflowsinthestraight-edgeforebayofpumpingstationandpreventivemeasuresofsedimentdeposition

XUCundong1，LIULuyao1，WANGGuoxia1，TIANZixun1，LIULixia2，CHENGHui1

(1.CollaborativeInnovationCenterofWaterResourcesEfficientUtilizationandProtectionEngineering，NorthChinaUniversityofWaterResourcesandElectricPower，Zhengzhou450045，China；2.InstituteofTechnology，GansuAgriculturalUniversity，Lanzhou730000，China)

Due to the sediment deposition occurring in front inflows of the straight-edge forebay of large pumping station, where its water is taken from the river with high content of sand, the adjustment ability of the forebay is usually reduced and the operation efficiency of the pumping station becomes low. To solve this problem, a prototype 3D model of the forebay structure is established using the ICEM CFD software, with its flow regime being modelled and analyzed by using Standardk-εturbulence model in the FLUENT software. Moreover, two optimization measures are proposed for the bad flow regime by incorporating diversion pier or water pressure plate. The findings show that the computed results are consistent with those from site survey, the flow regime in the forebay is in disorder and the sand sedimentation is quite serious. The eight-shaped diversion pier can reduce to some extent the scope of the main re-circulation region, but only serves as allocating of the forebay flow without obvious improvement in its flow filed distribution. On the contrast, the 45° water pressure plate can effectively improve the mainstream effect in the forebay, and the back-flow region on both sides is obviously reduced. The sand sedimentation problem in the forbay is thus solved with the favorable inflow conditions being ensured for the water pump.

the forebay of pumping station; sediment deposition; flow simulation; flow improvement; eight-shaped diversion pier; 45° water pressure plate

10.3876/j.issn.1000-1980.2017.05.004

2016-10-20

國(guó)家自然科學(xué)基金(51279102;51279064);河南省教育廳科技創(chuàng)新人才支持計(jì)劃(14HASTIT047);河南省高?？萍紕?chuàng)新團(tuán)隊(duì)支持計(jì)劃(14IRTSTHN028)

徐存東 (1972—) 男,甘肅景泰人,教授,博士,主要從事水工結(jié)構(gòu)優(yōu)化設(shè)計(jì)和耐久性研究。E-mail:xcundong@126.com

劉璐瑤,碩士研究生。E-mail:LiuLuyao8888@126.com

TV675

A

1000-1980(2017)05-0398-08