礦井多水源復(fù)雜并聯(lián)水泵的工況點簡易數(shù)解法研究

汪永彪

（中煤西安設(shè)計工程有限責(zé)任公司,陜西 西安 710054）

1 概述

在泵站工程設(shè)計中,水泵工況點的確定對水泵的效能分析及調(diào)速具有重要的意義[1-2]。水泵工況點是水泵加壓供水系統(tǒng)中在水泵穩(wěn)定運行狀態(tài)下的供需能量平衡點,在泵站設(shè)計后期運行中,需要將水泵工況點控制在水泵高效區(qū),以提高水泵的總效率,降低運行費用[3]。水泵的選型及后期泵站自動化運行管理中都需要確定水泵的工況點,以便于對水泵的適用性及效能分析做出準(zhǔn)確的評價,同時為水泵工況點調(diào)節(jié)提供有效依據(jù)[4]。在工程設(shè)計階段需要根據(jù)不同水泵及管路特性,計算確定水泵工況點,選擇最合適的水泵配置方案。在后期泵站管理和運行中,針對智慧供水及無人泵站的要求,為達(dá)到自動化及節(jié)能的目標(biāo),提出有一種適合于計算機應(yīng)用,并且簡單快速的工況點數(shù)解法,顯的尤為重要。

目前確定泵工況點的方法有圖解法和數(shù)解法,圖解法的優(yōu)點是直觀,圖解法適用于單個水泵或一座泵站內(nèi)多臺水泵并聯(lián)運行、水泵型號一致、泵后管道較短、連接簡單的情況。當(dāng)多個泵站并聯(lián)運行,并聯(lián)水泵型號差異較大、泵后管道較長、連接復(fù)雜時,圖解法存在最終結(jié)果誤差大及繪圖復(fù)雜的問題[5],圖解法需要人工手動繪制圖表,在方案選擇時需將每一個可行的水泵型號組合工況點求解圖紙全部繪制出來,需耗費較大時間,在后續(xù)自動化運行管理中無法程序化。

現(xiàn)有的數(shù)解法是將水泵特性曲線擬合為數(shù)學(xué)公式,聯(lián)合管道特性曲線求解工況點,可通過計算機來實現(xiàn),適用于單泵運行時計算工況點,對于多泵并聯(lián)情況,水泵特性曲線擬合復(fù)雜,其迭代次數(shù)多,運算量大[6],如涉及到泵出口管道連接復(fù)雜,不同型號的水泵并聯(lián),則現(xiàn)有的數(shù)解法無法求解。本文從水泵供水系統(tǒng)能量平衡角度進(jìn)行分析研究,給出一種簡單有效、適用范圍廣泛,針對復(fù)雜并聯(lián)運行水泵工況點的數(shù)解法。同時本文提出的數(shù)解法同樣適用于泵站運行時的變頻調(diào)速、工況點優(yōu)化改造、信息化提升等工作[7]。

在礦井等采用管井抽取地下水作為生產(chǎn)、生活用水時,往往需要建設(shè)多個水源井,利用深井泵將多個水源井內(nèi)的水抽出后,并聯(lián)至一根主干管輸送至用水點。每個水源井的出水量、水位及出水管道的長度、管徑等均不同,需根據(jù)不同情況選擇合適的水泵。在這種情況下找出一種多水源復(fù)雜并聯(lián)水泵工況點的計算方法,對于確定最優(yōu)水泵配置方案及后期泵站實現(xiàn)智慧化運行十分的必要性。

2 水泵工況點求解方法的探討

在泵站設(shè)計中,通過初步的水力計算結(jié)果,選擇某一型號的水泵后,便可得到擬選水泵的流量、揚程參數(shù),利用最小二乘法對水泵的流量、揚程參數(shù)進(jìn)行曲線擬合得到水泵特性曲線[8],擬合后的水泵特性曲線用式（1）表示：

根據(jù)泵后輸水管道參數(shù),可得到輸水管道的特性曲線,用式（2）來表示：

2.1 單泵運行工況點求解

當(dāng)只有單臺水泵運行時,水泵特性曲線與管道特性曲線只有一個交點,該交點便是水泵運行的工況點[9],工況點的流量和揚程需同時滿足式（1）與式（2）,則聯(lián)立式（1）和式（2）,便可求得工況點（H況,Q況）,見式（3）：

對式（3）的求解,在泵設(shè)計工況點Q況處用采用泰勒公式展開,并忽略高次項,將擬定工況點（H況, Q況）作為初始值,求得一個q,然后在初始分配流量下添加校正流量Δq,得到新的泵站輸送流量Δqi（k+1）=qi（0）±Δqi（k）,多次迭代后即可求得水泵工況點（H況, Q況）。

2.2 三臺水泵并聯(lián)運行工況點求解

對于三臺水泵并聯(lián)運行的情況見圖1。首先對泵A進(jìn)行分析,A號水泵將水從1點吸水后輸送到4點會流,然后再輸送到5點。泵A提供的水流經(jīng)過[1]、[4]號管段,將[1]、[4]號管段稱為水泵A的關(guān)聯(lián)管段,同樣 [4]號管段的水來自A、B、C三個水泵,稱A、B、C三個水泵為[4]管段的關(guān)聯(lián)水泵。從能量平衡的角度分析,當(dāng)水泵在工況點運行時,泵提供的能量一部分轉(zhuǎn)化為水的勢能,剩余的一部分用于抵消管道損失,全部消耗在水的輸送過程中,對于A水泵根據(jù)能量平衡得到式（4）。

圖1 三臺水泵并聯(lián)運行結(jié)構(gòu)圖

對于B、C號泵做同樣的能量平衡分析,可得出式（5）、式（6）。

將能量平衡式（4）、式（5）、式（6）聯(lián)立求解,可得出式（7）。

1）模型坐標(biāo)系：有限元模型坐標(biāo)系：模型建立在直角坐標(biāo)系（X，Y，Z）下，X 軸，沿河流方向并指向下游；Y 軸，垂直河流方向，指向山外；Z軸，豎直向上。

對式（7）分析進(jìn)行,在加壓輸水系統(tǒng)中,須遵守能量與質(zhì)量守恒定律,節(jié)點的流量必須滿足節(jié)點流量平衡,則可得到QA= q1,QB= q2,QC= q3, q4= q1+q2+q3=QA+QB+QC,則式（7）中的變量只有QA、QB、QC。

在實際工程應(yīng)用中,為便于求解,簡化計算,對于結(jié)算結(jié)果影響不大的高次項可忽略,將水泵特性曲線采用二次項擬合,將式（7）做簡化后成為式（8）。通過式（8）可看出,如只有一臺水泵運行,其余兩臺水泵停止的話,式（8）便可推導(dǎo)出式（3）,佐證了上述單泵運行工況點計算的分析方法是正確的。

將式（8）在設(shè)計水泵流量QA（0）、QC（0）、QD（0）處用采用泰勒公式展開,并且忽略展開式中的高次項,得到式（9）。

改寫為矩陣形式方程：

對式（8）中的F1（QA,QB,QC）求一階偏微分并帶入泵站擬定流量QA（0）、QB（0）、QC（0）,得到公式（11）,同樣可得到F2（QA、QB、QC）、F3（QA,QB,QC）的一階偏微分。

在求解線性方程（11）時,可得到ΔqA、ΔqB、ΔqC,將其施加到初始流量QA（0）、QB（0）、QC（0）上便可到了水泵工況點的流量值,帶入式（1）可得到相對應(yīng)的揚程。由于式（9）在忽略了泰勒公式展開的高次項,存在誤差,需要多次迭代,不斷校正誤差。泵站流量的校正公式見式（12）。

2.3 多臺水泵并聯(lián)運行工況點求解

對于L臺水泵通過復(fù)雜管道連接的聯(lián)合運行,在整個加壓供水系統(tǒng)中有k條水泵出水管道,m個局部損失點,可將式（8）擴(kuò)展為式（13）。

式（13）中當(dāng)?shù)趇個管段是第l個泵站關(guān)聯(lián)管段的時候,則管道流量qi不為0,否則為0。當(dāng)將式（13）在初分配流量下按照泰勒公式展開,并忽略高此項轉(zhuǎn)化為矩陣形式,見式（14）：

求出Δql（0）后通過,則可得出每條管段及水泵的新的流量,其中泵站的新的流量計算見式（16）。

式中：Ql（1）為泵站初始流量(水泵設(shè)計流量)。

式（16）中,L個水泵并聯(lián)供水系統(tǒng)中管道流量qi（0）根據(jù)節(jié)點流量平衡可得出,新的管道流量Qi（1）計算見式（17）。管段計流量的初始值根據(jù)節(jié)點平衡計算出后,為了保證節(jié)點流量平衡,在做管道校正流量時,只將與這根管道的關(guān)聯(lián)泵站的校正流量相加,通過多次迭代,便可得到工況點的流量參數(shù),再帶入水泵特性公式,便可得出工況點的流量參數(shù)。

2.4 復(fù)雜并聯(lián)水泵調(diào)速運行工況點數(shù)解

現(xiàn)有水泵的調(diào)速運行情況下,通過水泵性能試驗與相似理論（比例律公式） 得到不同轉(zhuǎn)速下的水泵性能特性曲線[10],其中水泵性能試驗的方法是在實驗室內(nèi)測試水泵在不同轉(zhuǎn)速下的性能曲線,該方法需要將水泵在特定的實驗臺進(jìn)行,在實際工程設(shè)計及水泵運行調(diào)整中可實施性不高。 在實際工程設(shè)計及水泵運行調(diào)整中使用較多的是利用相似理論（比例律公式）進(jìn)行調(diào)速水泵工況點求解,利用比例律公式對水泵調(diào)速運行工況點的求解有圖解法和數(shù)解法。

其中利用比例律公式對水泵調(diào)速運行工況點的圖解法,是在繪制水泵標(biāo)準(zhǔn)轉(zhuǎn)速下的水泵特性曲線后,根據(jù)需要調(diào)速后得到的流量、揚程參數(shù),繪制水泵相似工況剖物線,得到標(biāo)準(zhǔn)轉(zhuǎn)速的相似工況點,再根據(jù)相似工況點的流量與調(diào)速后的流量參數(shù)相關(guān)關(guān)系得到需要調(diào)整到的水泵轉(zhuǎn)速。該方法在多臺不同參數(shù)水泵并聯(lián)聯(lián)合調(diào)速運行時,無法使用。

本文以簡單水泵調(diào)速運行的數(shù)解法為基礎(chǔ),根據(jù)調(diào)速后的水泵特性曲線的,結(jié)合能量平衡,則可將式（13）擴(kuò)展為復(fù)雜并聯(lián)水泵調(diào)速工況點求解公式,見式（18）。式中的Ki是指不同調(diào)速水泵在標(biāo)準(zhǔn)轉(zhuǎn)速及調(diào)速后的水泵轉(zhuǎn)速的比指,調(diào)速工況點的求解方法與求解式（13）相同。

3 工程實例應(yīng)用

陜北某礦井生產(chǎn)、生活水源采用地下水,根據(jù)礦井周邊水文地質(zhì)情況及礦井用水情況需建設(shè)4口水源井,3用1備。因礦井工業(yè)場地位于山區(qū)溝道內(nèi),周邊地勢起伏較大,3口水源井建設(shè)位置較遠(yuǎn),沿溝底兩側(cè)布置水源井,每處水源處設(shè)一臺水泵深井泵,設(shè)一條主供水管道,將三個水源井的水分別匯集后輸送到礦井工業(yè)場地。每臺水源取水泵出口管道長度不同,三口水源井按并聯(lián)形式運行,管道布置及參數(shù)見圖2。三臺水泵選用相同型號的離心泵,根據(jù)樣本給出的三組水泵流量、揚程參數(shù),采用最小二乘法擬合的水泵特性曲線為：H=-3121.7Q2+452Q+68.1。對該工程中的多水源并聯(lián)水泵的運行工況點采用本文中的數(shù)解法計算,根據(jù)3次迭代后便可得出結(jié)算結(jié)果,數(shù)解法結(jié)果見表1。

表1 多水源復(fù)雜并聯(lián)水泵工況數(shù)解法結(jié)算結(jié)果

圖2 礦井水多水源取水泵站及輸水管道結(jié)構(gòu)圖

同時采用圖解法求解,根據(jù)水泵出水管道的特性,圖解法計算步驟如下：

（1）繪制出每臺水泵的特性曲線（Q-H）1；

（2）根據(jù)每臺水泵的出水口管道特性,分別繪出三臺水泵的折引曲線（Q-H）'i；

（3）根據(jù)折引曲線繪制出2號及3號水泵的并聯(lián)曲線（Q-H）'2+3；

（4）根據(jù)2號及3號水泵并聯(lián)管道的特性,繪制出2號及3號水泵并聯(lián)曲線的折引曲線（Q-H）''2+3；

（5）繪制出2號及3號水泵并聯(lián)曲線的折引曲線與1號水泵的折引曲線的并聯(lián)曲線（Q-H）''1+2+3；

（6）繪制并聯(lián)干管[5]的特性曲線,讀出并聯(lián)水泵特性點,圖解法結(jié)果見圖3。

圖3 對水源并聯(lián)水泵工況圖解法

從以上的計算過程及結(jié)果中可看到,對于僅有三個水源的并聯(lián)水泵工況點圖解法的計算,如采用圖解法計算,計算過程需多次繪制折引曲線,計算過程復(fù)雜,結(jié)果需人工讀取數(shù)值,誤差較大,當(dāng)并聯(lián)水源泵增加后,計算過程及工程量成倍增加。本文提出的數(shù)解法對于復(fù)雜并聯(lián)泵站工況點求解十分有效,結(jié)果準(zhǔn)確,誤差小。相比于圖解法,過程簡便,計算結(jié)果準(zhǔn)確度高,將計算過程編為程序后計算更加簡便,有利于泵站后期運行的自動控制,選擇最優(yōu)的運行工況,節(jié)約能源。

4 結(jié)論

本文提出的并聯(lián)水泵工況點數(shù)解法,不僅適用于簡單的單座泵站多臺水泵并聯(lián)工況點計算,尤其對于解決多水源供水、不同型號的水泵并聯(lián)運行、每處水泵的出水管道相差較大、并聯(lián)形式復(fù)雜的供水系統(tǒng),具有計算速度快、計算結(jié)果精度高的優(yōu)勢。同時,將原有的靠手工求解的過程通過數(shù)學(xué)方法解決,可將求解過程通過編程自動計算出水泵的工況點,在后期多水源泵站供水系統(tǒng)的控制、調(diào)速運行中,通過自動計算不同情況下的水泵運行工況點,可選擇出最合理、節(jié)能的水泵轉(zhuǎn)速及組合方式,對于智慧泵站建設(shè)及能效提高方面有著重要的意義。