泵閘聯(lián)合體中雙向引排水泵站進水系統(tǒng)的水力優(yōu)化研究

摘要：為了改善側(cè)向進流式泵站前池內(nèi)的繞流和回流等不良流態(tài)，基于 RNG k-ε模型對泵閘聯(lián)合體泵站運行的流態(tài)進行數(shù)值模擬，提出適用于側(cè)向收縮進流式泵站前池的整流措施。數(shù)值計算結(jié)果表明：泵站前池內(nèi)存在由于側(cè)向收縮引起的繞流，并向下旋流在進水池內(nèi)形成回流，致使泵站邊機組流道的軸向流速均勻度較低，進流條件不佳。在進流收縮區(qū)布設(shè)導(dǎo)流墩可有效改善繞流現(xiàn)象，抑制回流形成，提高進流均勻性；通過對導(dǎo)流墩長度進行優(yōu)化，并采取水工模型試驗驗證了優(yōu)化方案的整流效果。研究成果可為泵閘聯(lián)合體側(cè)向收縮進流式泵站進水系統(tǒng)的流態(tài)分析和水力優(yōu)化提供技術(shù)參考。

關(guān)鍵詞：泵閘聯(lián)合體；泵站進水系統(tǒng)；水力優(yōu)化；計算流體力學(xué)；水工模型試驗

中圖分類號：TV67 文獻標識碼：A文章編號：1001-9235（2024）06-0065-08

Hydraulic Optimization of Inflow System of Bidirectional Drainage Pumping Station in Pump-sluice Complex Engineering

XIAO Zhiqiao1， FEI Zhaodan2， QIAN Shangtuo3

（1. Shanghai Puhe Engineering Design Co.， Ltd.， Shanghai 200333， China;2. Institute of Water Science and Technology， HohaiUniversity， Nanjing 211100， China;3. College of Agricultural Science and Engineering， Hohai University， Nanjing 211100， China）

Abstract： To improve the adverse flow pattern such as circling and recirculation in the lateral forebay of the lateral inflow pumping station， numerous simulations of the flow during the operation of the pumping station in the pump-sluice complex engineering were conducted based on the RNG k-ε model， and the rectification measures for the lateral forebay of the lateral contraction inflow pumping station were proposed. The numerical results indicate that： the circling flow occurred in the forebay due to the effect of lateral contraction， resulting in recirculation flow forming in the suction sump due to the combination of circling flow and downward swirling flow. This leads to a lower axial-flow uniformity in the flow channels of the pumping station， indicating a poor inflow condition. Diversion piers in the inflow contraction area can effectively improve the flow pattern， inhibit the formation of recirculation flow， and increase the uniformity of inflow. The length of the diversion pier was optimized and the rectification of the optimized scheme was verified by the hydraulic model test. The research results can provide technical references for the flow analysis and optimization of the lateral inflow system of the pumping station in the pump-sluice complex engineering.

Keywords： pump-sluice complex engineering; inflow system of the pumping station; hydraulic optimization; computational fluid dynamics; hydraulic model test

平原地區(qū)河網(wǎng)交錯，為開發(fā)和利用水資源，經(jīng)常在靠近主干河道的支流口門上修建水閘和泵站等建筑物，一般采用水閘與泵站組合方式，稱為泵閘聯(lián)合體［1-2］。泵閘聯(lián)合體在低潮位時采用水閘自排，高潮位時采用泵站強排水［3］。泵閘聯(lián)合體中泵站和水閘共用進出水流道，相對于分建工程減少了占地面積和工程量，降低了工程造價，因而得到廣泛應(yīng)用［4-6］。泵閘聯(lián)合體中的流道共用，使得其流道寬度遠比單獨的泵站或水閘來得大。當泵站單獨運行時，泵站進水池主流寬度小于實際的泵閘聯(lián)合體流道寬度，水流難以在較短距離內(nèi)通過自身調(diào)整而充分收縮，將在上游形成較大范圍的回流和偏流等惡劣流態(tài)。尤其當泵閘聯(lián)合體的泵站需要實現(xiàn)雙向引排水時，進出水流動更為復(fù)雜。因此，研究泵閘聯(lián)合體雙向運行泵站進/出水系統(tǒng)的水力特性，探索有效的整流措施改善泵閘聯(lián)合體泵站運行時的流態(tài)，對保證實際工程安全運行具有重要意義［7-9］。

目前，學(xué)者們針對泵閘聯(lián)合體泵站進/出水流道流態(tài)及水流特性開展了大量研究。周楊等［10］對涵閘樞紐引水工況進行了數(shù)值研究，認為導(dǎo)流墻長度對引航道區(qū)域的流速分布影響很大。馮建剛等［11］分析了某水源泵站前池流態(tài)以及產(chǎn)生不良流態(tài)的原因，提出了消渦板、非全段式底坎和八字形導(dǎo)流墩等改善不良進水流態(tài)的工程措施，并對底坎和八字形導(dǎo)流墩措施進行比較分析。成立等［12］探討了 Y 形導(dǎo)流墩的分流效果，通過在前池加設(shè) Y形導(dǎo)流墩可顯著改善流態(tài)，并對 Y形導(dǎo)流墩的位置、長度等參數(shù)給出了建議值。張睿等［13］、李志祥等［14］通過數(shù)值模擬分析了泵站進水系統(tǒng)的水流流態(tài)和流道配水的均勻性，針對泵站進水系統(tǒng)流量分配不均提出組合式整流措施。資丹等［15］采用數(shù)值模擬和現(xiàn)場測試相結(jié)合的方法，設(shè)計了八字型導(dǎo)流墩、川字型導(dǎo)流墩和十字型消渦板相結(jié)合的組合式整流措施，發(fā)現(xiàn)該組合式整流措施減弱前池大尺度表面旋渦，提高了水流均勻性。羅燦等［16］分析了小型閘站式側(cè)向進水泵站進水流態(tài)，對比不同整流措施效果，以泵前縱剖面軸向速度分布均勻度、速度加權(quán)平均角和軸向平均流速為優(yōu)化目標，確定了合適的整流措施。上述研究大多都是對單獨的泵站進水系統(tǒng)進行分析，本文的泵閘聯(lián)合體的泵站運行為雙向運行，具有獨特性，且研究成果較少。

本文通過數(shù)值模擬結(jié)合試驗驗證的方式，對泵閘聯(lián)合體中雙向引排水泵站運行時的流態(tài)進行分析并開展水力優(yōu)化研究。通過數(shù)值模擬對泵站初步設(shè)計方案在排水和引水工況運行的進水系統(tǒng)的流態(tài)進行分析，針對初步設(shè)計方案中存在的不良流態(tài)進行水力優(yōu)化；采取布設(shè)導(dǎo)流墩的方式來進行整流研究，并通過水工模型試驗進一步驗證整流效果。研究所得成果可為泵閘聯(lián)合體側(cè)向收縮進流式泵站進水系統(tǒng)的流態(tài)分析與優(yōu)化提供技術(shù)參考。

1計算模型

1.1模型構(gòu)建

圖1所示為某泵閘聯(lián)合體工程的初步設(shè)計方案平面，泵閘聯(lián)合體的主要計算域包括了內(nèi)河部分河道、泵站、節(jié)制閘、套閘、外河部分河道。其中，泵站位于河道左側(cè)（順水流方向），順水流方向總長為234.00 m，設(shè)3臺雙向潛水貫流泵，安裝高程為-2.00 m，單機流量為13.33 m3/s，主要包括泵房、前池、進出水池、海漫段等；套閘位于河道右側(cè)（順水流方向），順水流方向總長為234.00 m，主要包括內(nèi)外閘首、內(nèi)外河消力池、內(nèi)外河海漫段、內(nèi)外河防沖槽等；節(jié)制閘位于套閘與泵站之間，順水流方向總長為234.00 m，主要包括閘室、內(nèi)外河消力池、海漫及防沖槽等。泵站運行時，套閘和節(jié)制閘關(guān)閉?；诖诉\行策略，建立泵閘聯(lián)合體泵站運行時的三維幾何模型，主要包括排水工況和引水工況的泵站進水系統(tǒng)模型，見圖2。

1.2數(shù)值計算方法

泵閘聯(lián)合體的水流流態(tài)復(fù)雜，通常處于湍流狀態(tài)，容易產(chǎn)生回流、旋渦、二次流等不良流態(tài)。在泵站工程水力學(xué)中，k-ε模型被廣泛應(yīng)用于相關(guān)研究中［17-19］?；?k-ε模型的 RNG k-ε模型則更多地考慮了湍流漩渦的影響［20-21］，可提高計算精度，因此選用 RNG k-ε模型進行數(shù)值模擬計算。采用有限體積法對控制方程進行離散，對流項采用高分辨率格式，流場求解采用全隱式多重網(wǎng)格耦合方法。

1.3邊界條件設(shè)定

泵閘聯(lián)合體泵站排水工況中，進口為內(nèi)河河道側(cè)，設(shè)置為壓力進口；出口為泵站進水管，設(shè)置為均勻流速出口，流速為設(shè)計流量與進水管橫截面積的比值；壁面設(shè)置為無滑移壁面。泵閘聯(lián)合體泵站引水工況中，進口為外河河道側(cè)，設(shè)置為壓力進口；出口為泵站出水管，設(shè)置為均勻流速出口，流速為設(shè)計流量與出水管橫截面積的比值；壁面設(shè)置為無滑移壁面，近壁面流動采用標準壁面函數(shù)進行處理。

由于本泵閘聯(lián)合體寬度較大、河床平穩(wěn)，其河道內(nèi)和進出水池的水面平穩(wěn)，波動較小，因此基于剛蓋假定將自由水面設(shè)置為對稱面。本研究中，選取泵閘聯(lián)合體泵站運行時的最不利水位進行流態(tài)分析，其中排水工況時水位為2.50 m，引水工況時水位為1.58 m。

1.4網(wǎng)格劃分

泵閘聯(lián)合體的幾何結(jié)構(gòu)復(fù)雜，不利于六面體網(wǎng)格的劃分，因此選用適應(yīng)性更好的四面體網(wǎng)格對三維幾何模型進行劃分，并對泵站進出水池等重點關(guān)注部位進行加密處理。網(wǎng)格數(shù)量對數(shù)值模擬的計算精度至關(guān)重要，有必要對其進行網(wǎng)格無關(guān)系分析。水力損失和流速是泵閘工程的重點關(guān)注因素，因此本研究選用前池和進水池段的水力損失及進水池進口斷面的平均流速進行網(wǎng)格無關(guān)性分析。

針對排水引水2種工況，各選取4套方案對網(wǎng)格劃分策略進行分析。表1所示為水力損失和平均流速隨網(wǎng)格數(shù)量的變化，結(jié)果表明當網(wǎng)格數(shù)量超過250萬（排水工況）、245萬（引水工況）時，誤差不超過1%，可認為超過此數(shù)量的網(wǎng)格方案滿足計算分析要求。綜合考慮計算精度和計算效率，本研究選取的計算網(wǎng)格數(shù)量為255.11萬（排水工況）、247.16萬（引水工況），并且基于此網(wǎng)格劃分策略進行水力優(yōu)化方案，不再進行相應(yīng)的網(wǎng)格無關(guān)性分析。此數(shù)量下的網(wǎng)格劃分策略見圖3。

1.5數(shù)值計算驗證

數(shù)值計算通過模型試驗結(jié)果進行計算準確度驗證。泵閘聯(lián)合體水工模型按照重力相似準則設(shè)計，模型線性比尺λL =25，建筑物采用透明有機玻璃制作。圖4a、4b分別為排水和引水工況的流態(tài)對比。試驗觀測（左）和數(shù)值模擬（右）的流態(tài)相似，可認定本數(shù)值模擬可準確模擬出泵閘聯(lián)合體泵站運行時的流態(tài)，具有可靠性。

2初步設(shè)計方案的流態(tài)分析

2.1排水工況

圖5所示為排水工況時進水系統(tǒng)的流態(tài)。受到河道進流偏移的影響，泵站前池內(nèi)存在明顯的回流區(qū)，主要集中于靠近節(jié)制閘側(cè)的3號流道上游（圖5a），隨著水流前進，當從前池下坡進入進水池后，回流區(qū)域擴大（圖5b）。這是由于水流自開闊的河道偏移進入泵站前池時，受到側(cè)向收縮的影響，靠近套閘側(cè)的水流產(chǎn)生水流集中現(xiàn)象，并最終形成回流。此外，由于前池底的高程要高于進水池底，水流向下旋滾進入進水池，導(dǎo)致前池內(nèi)的不良流態(tài)惡化，并蔓延至進水池中，形成了較大尺度的回流區(qū)域。

圖6所示為3條進水流道進口斷面的軸向流速分布云圖，3條流道進口斷面的流速均為正值，即與主流方向相同，并且均呈現(xiàn)頂部流速偏大的分布，這是由于流道高程較低導(dǎo)致的，并且可見3號流道斷面頂部的流速較1號和2號流道偏大。為了更好地對比泵站進流條件的優(yōu)劣，采用軸向流速均勻度 Vu 對進流條件進行定量評價［22］，軸向均勻度見式（1）：

Vu =（1-）×100%（1）

式中： a 為平均軸向流速；uai為第i個單元的軸向流速；m 為斷面總單元數(shù)。軸向流速分布越均勻，Vu 越接近于100%。

經(jīng)計算分析，1—3號流道進口斷面的軸向流速均勻度分別為83.10%、87.26%、77.27%?？傮w來看，排水工況1、2號流道的進流條件較好，均勻度可達80%以上，而3號流道受到回流區(qū)的影響，導(dǎo)致流道進口的流速均勻性不良，影響了正常進流。

2.2引水工況

圖7所示為引水工況進水系統(tǒng)水平斷面的流態(tài)。相較排水工況的流態(tài)，引水工況由于節(jié)制閘與套閘之間的分隔墩較長，導(dǎo)致河道主流的偏流程度降低，水流擠壓程度減弱，在套閘側(cè)無明顯的回流區(qū)，但是從圖7b 中仍可見套閘上游部分水流繞流進入2號流道，這是主流受擠壓和向下流動共同導(dǎo)致的結(jié)果。

圖8所示為3條出水流道出口斷面的流速分布云圖。在引水工況下，出水流道具有進流作用，此時的出口斷面為水流進入流道的進口斷面。受到工作閘門和事故檢修閘門之間存在滯水區(qū)的影響，引水工況的進口斷面的頂部存在明顯的負流速區(qū)，導(dǎo)致軸向流速均勻度整體偏小。經(jīng)計算分析，1—3號流道出口斷面的軸向流速均勻度分別為64.10%、62.68%、64.95%，其中，2號流道受到偏移進流的影響，導(dǎo)致軸向流速均勻度偏小。

3整流措施研究

3.1整流方案設(shè)計

為了改善3號流道在排水和引水2種工況下的偏移進流條件，提出在3號流道上游側(cè)增設(shè)導(dǎo)流墩的整流措施，通過對擠壓水流分流的方式改善進流條件。對于每種工況，通過改變墩頭位置的方式共設(shè)計了3套方案進行對比研究，導(dǎo)流墩結(jié)構(gòu)見圖9。

3.2整流措施效果分析

圖10所示為排水工況進水池區(qū)域的流態(tài)對比，對比初步設(shè)計方案，導(dǎo)流墩的整流效果明顯。整流方案 P1的導(dǎo)流墩墩頭與分隔墩墩頭持平，此時導(dǎo)流墩已呈現(xiàn)引流作用，部分水流自導(dǎo)流墩和分隔墩之間的流道流至進水池內(nèi)，3號流道進流流態(tài)得到一定的改善；與此同時，導(dǎo)流墩左側(cè)不可避免地產(chǎn)生低流速區(qū)域，并影響至貫流泵安裝高程。整流方案P2的導(dǎo)流墩墩頭超出分隔墩墩頭0.5 m，此時有更多的水流自導(dǎo)流墩和分隔墩之間的流道流進進水池內(nèi)，導(dǎo)流墩左側(cè)的低流速區(qū)域擴大。整流方案 P3的導(dǎo)流墩墩頭超出分隔墩墩頭1.5 m，此時導(dǎo)流墩左側(cè)低流速區(qū)更大，導(dǎo)致2號流道進流受到明顯影響。

三流道進口斷面的軸向流速均勻度情況見表2。3種整流方案在3號流道的軸向流速均勻度相較初步設(shè)計方案具有一定程度的提升，均大于80%，整流方案 P2提升最大，為87.19%；2號流道相較初步設(shè)計方案有所下降，這是由于2號流道上游存在回流導(dǎo)致，但是影響作用有限，軸向流速均勻度仍然大于80%；1號流道的軸向流速均勻度與初步設(shè)計方案相似，即不受導(dǎo)流墩的整流效果影響。綜上，整流方案 P2的進流均勻性最佳，故選用整流方案 P2作為排水工況的最終整流方案。

對于引水工況的流態(tài)，與排水工況類似，這里不再贅述。引水工況出口斷面的軸向流速均勻度見表3。相較初步設(shè)計方案和其他整流方案，整流方案 Y2在2、3號流道的軸向均勻度最高，并且在1號流道的軸向流速均勻度最高未明顯下降，即此方案的整流效果最佳，故選用整流方案 Y3作為引水工況的最終整流方案。

3.3整流方案驗證

對最終整流方案實施水工模型試驗，進一步驗證整流方案的有效性。圖11為排水和引水工況的模型試驗流態(tài)，圖12為排水和引水工況的前池垂線平均流速曲線。結(jié)合圖11、12，最終整流方案的回流和繞流相較初步設(shè)計方案明顯改善，垂向平均流速分布亦更加均勻，即整流措施效果明顯，也再次驗證了數(shù)值計算的準確性。

通過數(shù)值模擬研究和水工模型試驗驗證，本文所得最終整流方案已滿足泵閘聯(lián)合體安全穩(wěn)定運行的要求，可用作實際工程的最終設(shè)計方案。進一步，基于不同的運行水位，對導(dǎo)流墩的位置和結(jié)構(gòu)進行多參數(shù)的水力優(yōu)化，可獲取最優(yōu)整流設(shè)計方案。本項目將在下一階段進行相應(yīng)研究。

4結(jié)論

a）泵閘聯(lián)合體泵站獨立運行時，無論排水還是引水工況，受到來流側(cè)向收縮的影響，在分隔墩附近產(chǎn)生繞流和回流區(qū)，導(dǎo)致分隔墩處的進水流道流態(tài)較差，流速分布不均。

b）對于排水工況，在回流區(qū)上游增設(shè)導(dǎo)流墩后，回流區(qū)域明顯收縮，進流均勻性顯著提升，整流效果較好。將此整流方案布設(shè)于引水工況的進水池中，可得類似的整流效果。

c）通過水工模型試驗對導(dǎo)流墩整流效果進行了驗證，導(dǎo)流墩整流措施可有效改善泵閘聯(lián)合體側(cè)向泵站進流的擾流和回流現(xiàn)象，提高泵閘聯(lián)合體水泵運行的穩(wěn)定性，可保障泵站的安全高效運行。

參考文獻：

［1］馮建剛，李志祥，錢尚拓.泵閘工程閘下雙底坎整流的水力特性［J］.水利水電科技進展，2020，40（1）：25-31，36.

［2］馮建剛，孟湘云，錢尚拓.分側(cè)式閘站樞紐下游底坎整流特性［J］.水利水電科技進展，2019，39（6）：62-67.

［3］謝先坤.物理模型試驗在不對稱泵閘布置驗證及優(yōu)化中的應(yīng)用［J］.中國農(nóng)村水利水電，2022（5）：77-81，86.

［4］戴揚濤.雙向一體化泵閘在東莞市中堂鎮(zhèn)內(nèi)河涌治理中的應(yīng)用［J］.中國水運（下半月），2024，24（4）：93-94，111.

［5］方正杰.上海水利泵閘智慧運維管理模式的探索與思考［J］.珠江水運，2023（5）：23-25.

［6］陳靜，漢京超，周娟娟，等.一體化泵閘的發(fā)展及應(yīng)用研究［J］.中國給水排水，2023，39（18）：49-54.

［7］蘇正洋，陸嘉偉，張志韜，等.泵站前池整流技術(shù)研究綜述［J］.人民珠江，2020，41（1）：69-75.

［8］鄧鵬，張詩媛，王旭瀅.水利泵閘工程安全管理中數(shù)字孿生技術(shù)的應(yīng)用研究［J］.水利信息化，2024（2）：16-20，35.

［9］張新嬌.基于數(shù)值模擬的某彎道河段泵閘站進水流態(tài)改善措施［J］.水利技術(shù)監(jiān)督，2023（11）：119-123.

［10］周楊，楊曉紅，吳佩鋒.涵閘樞紐引水工況水力優(yōu)化數(shù)值模擬研究［J］.人民珠江，2017，38（4）：43-47.

［11］馮建剛，李杰.大型城市水源泵站前池流態(tài)及改善措施試驗［J］.水利水電科技進展，2010，30（2）：70-74.

［12］成立，祁衛(wèi)軍，羅燦，等. Y形導(dǎo)流墩幾何參數(shù)對泵站前池整流效果的影響［J］.水利水電科技進展，2014，34（1）：68-72.

［13］張睿，徐輝，陳毓陵，等.斜向進流城市排水泵站進水箱涵配水均勻性的改善研究［J］.灌溉排水學(xué)報，2017，36（12）：75-80.

［14］李志祥，馮建剛，錢尚拓，等.斜向進流城市雨水泵站引水建筑物流態(tài)改善數(shù)值模擬［J］.中國農(nóng)村水利水電，2020（6）：139-143，147.

［15］資丹，王福軍，姚志峰，等.大型泵站進水流場組合式導(dǎo)流墩整流效果分析［J］.農(nóng)業(yè)工程學(xué)報，2015，31（16）：71-77.

［16］羅燦，劉浩，施偉，等.雙側(cè)向進水泵站前池流態(tài)數(shù)值模擬研究［J］.中國農(nóng)村水利水電，2020（11）：121-128.

［17］CHEN Y L，LIU G Y，QIAN S T，et al. Effects of convergingsidewalls on skimming flow over converging stepped spillway［J］ Applied Science-Basel，2022，12. DOI：10.3390/app12157868.

［18］FEI Z D，XU H ，ZHANG R ，et al. Numerical simulation on hydraulic performance and tip leakage vortex of a slanted axial- flow pump with different blade angles［J］. Proceedings of the InstitutionofMechanicalEngineering ， PartC： Journalof Mechanical Engineering Science，2022，236（6）：2775-2790.

［19］代彬，陳章淼，周維.基于 Realizable k-epsilon 模型的水閘下游水流數(shù)值模擬［J］.水利與建筑工程學(xué)報.2018，16（4）：176-180.

［20］王亞軍.基于 RNG k-epsilon模型的高海拔地區(qū)生態(tài)基流下游階梯消能水力特性數(shù)值研究［J］.水利技術(shù)監(jiān)督，2019（6）：180-183.

［21］馮建剛，曹亞萍，王曉升，等.深層隧道進流式泵站前池流態(tài)及整流措施研究［J］.水利水電科技進展，2024，44（2）：1-6，21.

［22］陸林廣，周濟人.泵站進水流道三維紊流數(shù)值模擬機水力優(yōu)化設(shè)計［J］.水利學(xué)報，1995（12）：67-75.

（責(zé)任編輯：向飛）