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

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

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

關鍵詞：進水流道；吸水池；流量衰減；表面吸入渦

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

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

０ 引言

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

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

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

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

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

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

１ 研究對象及存在問題

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

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

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

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

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

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

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

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

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

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

４ 流道改造及流動分析

４．１ 進水流道改造

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

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

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

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

５ 流道改造前后試驗對比

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

６ 結論

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

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

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

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

參考文獻：

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

【責任編輯 簡 群】