泵站流量衰減問題分析與進水系統(tǒng)改造

白綿綿 肖若富 黨照 姚志峰 馮虎虎

：為解決某抽黃灌區(qū)樞紐泵站機組流量大幅衰減問題，采用計算流體動力學(xué)（ＣＦＤ）方法對其衰減原因進行分析，提出進水系統(tǒng)改造方案，并對改造前后的進水系統(tǒng)流態(tài)進行數(shù)值模擬對比。結(jié)果表明：改造前，在箱涵進口產(chǎn)生明顯吸入渦，大量氣體和雜質(zhì)經(jīng)箱涵高速進入吸水池后，造成池內(nèi)水流上下翻滾，導(dǎo)致水泵吸水性能極差，運行一段時間后流量迅速衰減；改造后，進水系統(tǒng)流態(tài)得到有效改善，進水池低速區(qū)變小，吸水池水流上下翻滾現(xiàn)象消失，水泵吸水性能提升較大，機組流量衰減現(xiàn)象得到明顯改善。通過實測改造前后泵站的運行流量，驗證了ＣＦＤ 用于泵站內(nèi)部流場模擬的可靠性；在后續(xù)對原有工程進行改造時，可先采用ＣＦＤ 對原始系統(tǒng)進行流場模擬，有針對性地進行結(jié)構(gòu)改造。該泵站流道改造前后效果驗證，進水系統(tǒng)設(shè)計時，應(yīng)盡量保持各進水建筑物壓力流或明流相同。

關(guān)鍵詞：進水流道；吸水池；流量衰減；表面吸入渦

中圖分類號： ＴＶ２２２．１ 文獻標(biāo)志碼：Ａ ｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２３．１０．０２３

引用格式：白綿綿，肖若富，黨照，等．泵站流量衰減問題分析與進水系統(tǒng)改造［Ｊ］．人民黃河，２０２３，４５（１０）：１２６－１３０．

０ 引言

我國大型灌區(qū)多建于２０ 世紀(jì)六七十年代［１］ ，受當(dāng)時條件限制，進水系統(tǒng)設(shè)計在很大程度上依賴規(guī)范和經(jīng)驗，造成水泵在實際運行時吸水性能差、效率較低［２－３］ 。某大型抽黃灌區(qū)多座泵站就存在嚴重的流量衰減問題，即機組剛開機時各項參數(shù)正常，但運行一段時間后（幾天甚至幾小時），水泵流量逐漸減小，配套電機電流隨之逐漸降低，之后穩(wěn)定在某個遠小于開機流量值運行，無法滿足灌溉高峰期用水要求，對灌區(qū)效益和社會穩(wěn)定造成一定影響［２］ 。

經(jīng)試驗發(fā)現(xiàn)，引起流量衰減的原因為：１） 進氣。當(dāng)雙吸泵吸程為正時，進氣在啟動很短時間內(nèi)發(fā)生，且不太可能通過關(guān)機重啟而消失。２）表面旋渦。表面旋渦與進水池、吸水池、進水管道結(jié)構(gòu)密切相關(guān)，通常伴有明顯的振動和壓力脈動。３）汽蝕。汽蝕一般發(fā)生在機組選型或泵站設(shè)計不合理時，伴有異常高頻聲音，產(chǎn)生大量堵塞流道的氣泡，外觀常見明顯破壞。４）流道堵塞。流道堵塞分為管道堵塞和水泵堵塞兩類，前者流量衰減前后的工作點都在水泵特性曲線上，停機重啟不易恢復(fù)；后者流量衰減前后的工作點不在水泵特性曲線上，通過關(guān)機重啟可恢復(fù)。５）其他不確定的原因。

經(jīng)反復(fù)觀察，該抽黃灌區(qū)多座泵站機組出現(xiàn)流量衰減后，關(guān)機重啟就可恢復(fù)正常運行流量，排除進氣因素；葉輪運行２０ 多ａ 無明顯磨損和汽蝕破壞，排除汽蝕因素；吸水池、流道等無堵塞雜物，排除流道堵塞因素；進水池和吸水池內(nèi)旋渦、回流及翻滾現(xiàn)象嚴重，初步懷疑表面旋渦等進水流場紊亂是引起該灌區(qū)機組流量衰減的主要原因。

表面旋渦分為自由表面旋渦和次表面旋渦兩種。表面吸入渦是自由表面旋渦６ 個階段中的最強階段［４－５］ ，旋渦從吸水池自由水面一直延伸到進水喇叭口，形成連續(xù)的空氣吸入渦，夾帶大量氣體進入水泵，會在一定程度上堵塞流道，同時，吸入渦的強旋轉(zhuǎn)特性導(dǎo)致水泵進口流速變得極不均勻，造成水泵出水流量減小、運行效率下降、壓力脈動和振動加劇。

資丹［６］ 采用耦合Ｌｅｖｅｌ－ｓｅｔ 和ＶＯＦ 方法對表面旋渦的產(chǎn)生機制進行了研究，認為主導(dǎo)渦量分量為垂直渦量，而流向速度沿展向的波動是垂直渦產(chǎn)生的原因；流向進水管的主流將垂直渦輸送到進水管附近，在進水管上方區(qū)域的垂直渦被拉伸，其強度隨著拉伸而增大，當(dāng)增大到一定值時，空氣會被垂直渦吸入到水下或進水管中。多位學(xué)者通過數(shù)值模擬，研究了進水系統(tǒng)幾何形狀、參數(shù)與進水流場的關(guān)系，并提出多種改善進水流態(tài)的措施。賈君德等［７－１１］ 研究了箱形、簸箕形、開敞式等不同形式進水部分的水力特性，通過對長度、寬度、高度、擴散角、底坡等系列參數(shù)進行分析，認為各參數(shù)之間存在最優(yōu)比。朱紅耕等［１２－１６］ 對進水部分進行了整流優(yōu)化研究，提出設(shè)置底坎、導(dǎo)流墩、消渦板等多種控制渦流的措施和“八”字形導(dǎo)流墩＋“川”字形導(dǎo)流墩＋“十”字形消渦板的組合控渦方案。陸林廣等［１７］通過物理模型試驗，驗證了數(shù)值模擬的準(zhǔn)確性。針對情況更為惡劣的低水位工況，Ｍｏｌｌｅｒ 等［１８－２１］ 提出在吸水池設(shè)置導(dǎo)流臺＋水下消渦板＋Ｗ 形后墻導(dǎo)流墩的組合整流方案，并采用三維模擬和模型試驗對整流效果進行了驗證，結(jié)果表明，該組合整流方案可以完全消除吸水池表面旋渦，吸水池內(nèi)水流流態(tài)和水泵進口的軸向速度分布都得到了明顯改善。

然而，目前沒有對泵站流量衰減與進水流場之間的具體相關(guān)性研究。筆者針對上述表面旋渦等進水流場紊亂導(dǎo)致的流量衰減問題，結(jié)合某抽黃灌區(qū)樞紐泵站（以下簡稱某泵站），對進水部分進行計算流體動力學(xué)（ＣＦＤ）數(shù)值模擬，分析流場特點，研究造成流量衰減的原因，通過改造進水流道來降低流量衰減率；并通過現(xiàn)場試驗，對ＣＦＤ 模擬結(jié)果和改造效果進行了驗證。

１ 研究對象及存在問題

某泵站于１９９７ 建成運行，安裝８ 臺１２００ ＬＷ－６０型立式離心泵，單泵設(shè)計流量４．０５ ｍ３ ／ ｓ、設(shè)計揚程６０．３ ｍ。泵站進水采用單機單池單流道形式，自上級渠道引水，水流先后通過進水渠道、進水池、壓力箱涵進水流道、吸水池，進入水泵吸入口（見圖１），其中壓力箱涵進水流道相關(guān)參數(shù)為底寬Ｂ×高度Ｈ ＝ ２．５ ｍ×２．０ ｍ、長１８．７ ｍ、縱坡比１ ∶ ５．８６（壓力箱涵出口與橢圓形吸水池底部連接）。

測點設(shè)置在水泵出口與電動檢修閥之間，流量計型號ＡＤＭ６７２５，測量原理為時差相關(guān)原理，測量精度±０．５％，分辨率０．０２５ ｃｍ／ ｓ。

２．２ ＣＦＤ 模型及邊界條件

以７＃、８＃機組為例，考慮設(shè)計水位為３８２．５ ｍ 時雙機運行工況，對其進水系統(tǒng)進行真實尺寸的數(shù)值模擬。原始流道計算模型包括進水渠道、進水池、壓力箱涵進水流道、吸水池及進水管路（見圖２）。

湍流模型采用Ｒｅａｌｉｚａｂｌｅ ｋ－ε 湍流模型，有關(guān)研究表明該模型能較好地模擬表面旋渦和附壁渦，且預(yù)測的旋渦位置和強度與模型試驗較為吻合，在進水系統(tǒng)流動模擬中具有一定優(yōu)勢［５，１１，１７］ 。

對流項采用二階迎風(fēng)格式進行離散，源項和擴散項采用中心插分格式進行離散，采用ＳＩＭＰＬＥＣ 算法進行求解。邊界條件設(shè)置：１）進口邊界為進水池上游進水渠道６０ ｍ 處斷面，給定流量；２）出口邊界為進水管路出水?dāng)嗝妫o定壓力；３）進水渠道和進水池頂部采用鋼蓋假定；４）吸水池水面給定大氣壓力；５）其他壁面為無滑移壁面。

采用四面體的非結(jié)構(gòu)化網(wǎng)格，在各部位銜接處（進水渠道與進水池銜接處、進水池與進水流道銜接處、進水流道與吸水池銜接處以及吸水池與進水管路銜接處） 等流動變化劇烈的位置對網(wǎng)格進行局部加密，并進行網(wǎng)格無關(guān)性檢查，最終選用網(wǎng)格單元數(shù)為５００ 萬個，其中進水渠道１２４ 萬個、進水池１７３ 萬個、進水流道２４ 萬個、吸水池１６５ 萬個、進水管路１４ 萬個。

３ 原進水系統(tǒng)內(nèi)部流動分析

圖３ 為原始流道進水系統(tǒng)流線圖。從圖３ 可以看出，以進水池靠近箱涵進口處為分界點，前后流場截然不同。分界點前的進水渠道和進水池前半段，流線整體較為平順，流速在進入進水池后平面擴散逐漸減小；分界點后流線紊亂，尤其在進水池靠近箱涵進口處和吸水池內(nèi)流線更加混亂，壓力箱涵內(nèi)流速較吸水池高，水泵進口處流速最高。４ 為原始流道進水系統(tǒng)流場細節(jié)圖，進水池內(nèi)流速均值為５．５ ｍ／ ｓ，流速在垂直方向從池底到水面自下而上逐漸提高，且底部流速梯度較大，最底部流速為０．１ ｍ／ ｓ；自上游至下游流速均勻部分占比越來越小；在末端靠近箱涵進口處，形成了一個流速小于０．２ ｍ／ ｓ 的低速旋渦區(qū)，容易造成泥沙淤積；進入箱涵后，流速迅速提高至９．３ ｍ／ ｓ，挾帶表面吸入渦的水流從箱涵進水流道中高速流出，快速進入吸水池。吸水池流場縱向明顯分為３ 部分，最底部小區(qū)域流速低，往上比箱涵高度較低處流速高，再往上流速又變低；由于下部為高速流動區(qū)、上部為低速流動區(qū)，因此大的流速差造成下部水流自內(nèi)圈流向上部，再在重力作用下自外圈向下流動，形成翻滾式流道，造成水泵吸水性能下降。

原始流道進水池靠近箱涵進口處出現(xiàn)典型的表面吸入渦［６］ ，水流在水面背離軸線（順?biāo)鞣较颍┳酝庀騼?nèi)旋轉(zhuǎn)，水平向隨旋轉(zhuǎn)半徑變小速度逐漸提高，豎直向整個旋渦自水面向下尺寸越來越小、速度越來越高，直至高速進入壓力箱涵，隨水流一直延伸到進水喇叭口。吸水池內(nèi)多處水流上下翻滾，水流在翻滾中心自下而上流動，又沿著翻滾外緣自上而下流動。實拍照片（見圖５）與ＣＦＤ 模擬結(jié)果一致，證實數(shù)值模擬可靠，可用來指導(dǎo)下一步進水系統(tǒng)改造。

４ 流道改造及流動分析

４．１ 進水流道改造

原始進水流場特征表明，某泵站流場紊亂的源頭在進水池靠近箱涵進口處，即明流變壓力流處。本次改造擬將整個進水系統(tǒng)流態(tài)統(tǒng)一變?yōu)槊髁骰驂毫α鳌Ｓ捎诔M水流道外其余進水建筑物皆為明流，考慮投資最省，因此將原來的壓力箱涵進水流道改造為明渠進水流道，整個進水系統(tǒng)變?yōu)槊髁鳎ㄒ妶D６）。改造后明渠橫斷面為矩形，底寬Ｂ 為１．２ ｍ，明渠之間采用塊石混凝土填筑，底坡維持１ ∶ ５．８６ 不變。

４．２ 流道改造后進水系統(tǒng)內(nèi)部流動分析

為研究進水流道改造前后流態(tài)差異，對７＃、８＃機組在設(shè)計水位雙機運行工況下進行模擬。由進水系統(tǒng)流線圖（見圖７）看出，隨著進水流道由壓力箱涵改為明渠，進水池與流道銜接處流線也由紊亂變得極為平滑；從圖８ 可以看出，進水流道、進水池中水流速度較為均勻，與原始流道相比，改造后進水池靠近流道進口處底部的低速流動區(qū)域明顯縮小，不到改造前的１／７，泥沙淤積現(xiàn)象減輕。在進水流道底部出現(xiàn)了一部分低速區(qū)，但此處淹沒深度較大，不會對進水系統(tǒng)的整體流動產(chǎn)生較大影響。

由圖９、圖１０ 可知，流道改造后，吸水池沿縱向流速基本一致，上下翻滾式流動消失；然而，池內(nèi)仍存在較為明顯的吸入渦，原因可能是帶入空氣或者雜質(zhì)到進水管路中，導(dǎo)致水泵吸水性能變差，流量衰減現(xiàn)象雖然有所改善，但沒有完全消除。

５ 流道改造前后試驗對比

為進一步研究流道改造效果，對流道改造前后的實測數(shù)據(jù)進行對比，驗證了數(shù)值模擬的準(zhǔn)確性（見表２）。流道改造后，水泵開機流量增大，流量衰減速度變慢，運行７ ｄ 后流量衰減率由改造前的３３％減小到２０％，衰減幅度變小，衰減率降低了３９％，衰減問題大幅改善。但是，由于吸水池內(nèi)還可能存在其他引起流場紊亂的因素，因此流量衰減并沒有完全消除。

６ 結(jié)論

針對已建泵站進水系統(tǒng)流態(tài)紊亂引起的機組流量衰減問題，采用ＣＦＤ 模擬和試驗測試相結(jié)合的方法，對某抽黃灌區(qū)樞紐泵站進水系統(tǒng)流態(tài)進行研究，提出將進水流道由壓力箱涵改造成明渠的方案，從而將整個進水系統(tǒng)流態(tài)統(tǒng)一變?yōu)槊髁鳎档土髁克p率。

１）改造前，在進水池靠近箱涵進口處出現(xiàn)明顯低速區(qū)，在明流變壓力流的箱涵進口產(chǎn)生明顯的吸入渦，帶入大量氣體和雜質(zhì)隨著箱涵內(nèi)壓力流高速進入吸水池，造成吸水池內(nèi)水流上下翻滾，降低水泵吸水性能，隨開機時間延長，流量衰減嚴重。

２）將進水流道由壓力箱涵改為明渠后，進水系統(tǒng)流場較改造前平穩(wěn)。進水池低速區(qū)明顯變小，進水流道底部出現(xiàn)一部分對整體流動影響較小的低速區(qū)。吸水池內(nèi)翻滾式流動消失，水泵吸水性能有明顯提高。

３）通過對流道改造前后開機流量和運行７ ｄ 后流量實測數(shù)據(jù)進行對比，驗證了數(shù)值模擬的正確性。流道改造后，開機流量變大，流量波動幅度變小，流量衰減現(xiàn)象明顯改善，運行７ ｄ 后流量衰減率由改造前的３３％降低到２０％，衰減率降低３９％。

參考文獻：

［１］ 李琪，許建中，李端明．全國大型泵站更新改造前期工作有關(guān)問題探討［Ｊ］．中國農(nóng)村水利水電，２００７（１２）：４１－４４，４８．

［２］ 丁天壽．１２００ＬＷ－６０ 大型泵低效原因淺析及改善措施探討［Ｊ］．西北水力發(fā)電，２００６，２２（５）：８９－９１．

［３］ 樊亞龍．三元流葉輪改造技術(shù)在灌溉泵站中的應(yīng)用研究［Ｊ］．陜西水利，２０１８（１）：１７１－１７２，１７６．

［４］ 王福軍．水泵與水泵站［Ｍ］．北京：中國農(nóng)業(yè)出版社，２０２１：２９５－２９７．

［５］ 王福軍，唐學(xué)林，陳鑫，等．泵站內(nèi)部流動分析方法研究進展［Ｊ］．水利學(xué)報，２０１８，４９（１）：４７－６１，７１．

［６］ 資丹．泵站進水池空氣吸入渦水力特性研究［Ｄ］．北京：中國農(nóng)業(yè)大學(xué)，２０１９：３９－４５．

［７］ 賈君德．箱型進水流道水力特性數(shù)值分析研究［Ｄ］．揚州：揚州大學(xué)，２０１９：２１－５２．

［８］ 高傳昌，孫龍月，董旭敏，等．簸箕形進水流道寬度對流道水流流態(tài)的影響［Ｊ］．人民黃河，２０２０，４２（１）：９３－９６，１０４．

［９］ 黃佳衛(wèi)，劉超，嚴天序，等．進水流道對泵裝置性能影響的數(shù)值模擬分析［Ｊ］．水利水電技術(shù)，２０１８，４９（１０）：１１０－１１９．

［１０］ 付建國，李寧．泵站前池和開敞式進水流道的數(shù)值優(yōu)化［Ｊ］．江西科學(xué)，２０２２，４０（１）：１６５－１７０．

［１１］ 董旭敏．引黃低揚程泵站進水流道數(shù)值模擬研究［Ｄ］．鄭州：華北水利水電大學(xué)，２０２０：２２－４９．

［１２］ 朱紅耕，袁壽其．大型泵站進水流道技術(shù)改造優(yōu)選設(shè)計［Ｊ］．水力發(fā)電學(xué)報，２００６，２５（２）：５１－５５．

［１３］ 羅燦，成立，劉超．泵站正向進水前池底坎整流機理數(shù)值模擬［Ｊ］．排灌機械工程學(xué)報，２０１４，３２（５）：３９３－３９８．

［１４］ 王紫陽，楊凱，張興毅，等．整流罩在泵站肘形進水流道中的應(yīng)用效果分析［Ｊ］．水利水電技術(shù)，２０１８，４９（２）：１１２－１１６．

［１５］ 羅海軍，張睿，徐輝．改善城市排水泵站進水流態(tài)的試驗研究［Ｊ］．中國農(nóng)村水利水電，２０１９（１）：１７６－１７９．

［１６］ 資丹，王福軍，姚志峰，等．大型泵站進水流場組合式導(dǎo)流墩整流效果分析［Ｊ］．農(nóng)業(yè)工程學(xué)報，２０１５，３１（１６）：７１－７７．

［１７］ 陸林廣，徐磊，梁金棟，等．泵站進水流道三維流動及水力損失數(shù)值模擬［Ｊ］．排灌機械，２００８，２６（５）：５５－５８．

［１８］ ＭＯＬＬＥＲ Ｇ．Ｖｏｒｔｅｘ?Ｉｎｄｙｃｅｄ Ａｉｒ Ｅｎｔｒａｉｎｍｅｎｔ Ｒａｔｅ ａｔ Ｉｎｔａｋｅｓ［Ｄ］．Ｚｕｒｉｃｈ：Ｓｗｉｓｓ Ｆｅｄｅｒａｌ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｔｅｃｈｎｏｌｏｇｙ Ｚｕｒｉｃｈ，２０１３：２３－２７．

［１９］ 高傳昌，解克宇，黃丹，等．不同水位對泵站進水池流態(tài)影響的數(shù)值模擬［Ｊ］．華北水利水電大學(xué)學(xué)報（自然科學(xué)版），２０１５，３６（６）：１０－１４．

［２０］ 胡良明，柴端伍．不同水位泵房進水流道堵塞對流速分布的影響［Ｊ］．人民黃河，２０１９，４１（８）：１５３－１５６，１６２．

［２１］ 高傳昌，曾新樂，解克宇，等．泵站進水池超低水位下組合整流方案與驗證［Ｊ］．農(nóng)業(yè)工程學(xué)報，２０１７，３３（２３）：１０１－１０８．

【責(zé)任編輯 簡 群】