王 鑫,韓偉實,何藝峰,2,幸奠川
(1.哈爾濱工程大學(xué)核科學(xué)與技術(shù)學(xué)院,黑龍江哈爾濱150001;2.西北核技術(shù)研究所,陜西西安710024)
離心泵是給水系統(tǒng)中的重要設(shè)備,在實際應(yīng)用中為了保證連續(xù)供水、提高給水系統(tǒng)穩(wěn)定性和可靠性,經(jīng)常設(shè)置兩臺相同的離心泵,一臺運行,另一臺備用。當(dāng)一臺給水泵停機時另一臺能夠迅速投入使用,保證系統(tǒng)繼續(xù)穩(wěn)定供水[1]。在運行泵發(fā)生故障停機、不停止供水檢修某一臺給水泵以及均衡兩臺給水泵壽命等條件下,需要進行給水泵切換。
給水泵切換過程中會產(chǎn)生給水流量波動,尤其對于某些需要嚴(yán)格穩(wěn)定流量輸出的系統(tǒng),這種流量波動將產(chǎn)生較大影響。目前很多學(xué)者已對離心泵瞬態(tài)理論分析及數(shù)值模擬進行了深入研究:田文喜、蘇光輝等針對雙泵并聯(lián)系統(tǒng)進行了水錘因素分析[2-3];吳大轉(zhuǎn)等數(shù)值模擬了離心泵在瞬態(tài)操作條件下內(nèi)部流動,與試驗結(jié)果吻合[4];楊洪波等分析了核電廠常規(guī)電動給水泵切換瞬態(tài)特性[5];張龍飛等研究了轉(zhuǎn)動慣量對船用核動力主泵瞬態(tài)特性的影響[6]。
影響切換過程中流量特性的因素較多,主要包括切換時的運行工況、給水泵轉(zhuǎn)動慣量、切換方式等,本文著重研究不同的切換方式對雙泵切換過程流量特性的影響。建立雙泵并聯(lián)給水系統(tǒng)模型,通過實驗數(shù)據(jù)對模型進行驗證,利用驗證后的模型對不同切換方式下流量特性進行研究,為實際工程應(yīng)用提供參考。
本試驗中設(shè)計雙泵并聯(lián)給水系統(tǒng),為防止支路倒流造成給水泵倒轉(zhuǎn),在每個支路給水泵下游都安裝有一個止回閥。
試驗裝置結(jié)構(gòu)如圖1所示。以常溫下水為工質(zhì),在接近大氣壓力條件下進行試驗研究。兩臺給水泵的額定轉(zhuǎn)速1 450r/min。電磁流量計最高流速15m/s,精確度為0.3級(流速≥1m/s),環(huán)境溫度-10~+60℃。采用3595系列IMP數(shù)據(jù)采集器作為動態(tài)數(shù)據(jù)采集系統(tǒng),對回路流量參數(shù)進行檢測。
按照圖1將試驗裝置連接并檢查,保證一切設(shè)備正常工作。試驗中在1#泵運行、2#泵備用條件下開始,先關(guān)閉1#離心泵,再啟動2#離心泵。
建立模型需要做適當(dāng)?shù)暮喕图僭O(shè):①給水系統(tǒng)工質(zhì)為常溫常壓下水;②流體處于單相狀態(tài);③流體沿回路看作一維流動??紤]到離心泵和止回閥的特性和計算精度,通過計算將時間步長定為0.1s。
圖1 試驗裝置結(jié)構(gòu)示意圖Fig.1 Scheme of experimental device
為進行數(shù)值模擬,需要對整個給水系統(tǒng)中的關(guān)鍵設(shè)備建立模型。
1.2.1 離心泵模型
離心泵是整個回路中流體流動的動力來源。目前離心泵瞬態(tài)過程的相關(guān)研究較多[7-9],在泵開啟過程中,流體流動行為主要是由轉(zhuǎn)動部件的慣性及流體慣性控制。離心泵瞬態(tài)方程可以寫成[10]:
式中:g為重力加速度,m/s2;A為流通截面,m2;De為等效管道直徑,m;q為質(zhì)量流量,kg/s,t為時間,s;hp為泵揚程,m;f為摩擦因子;De為等價管道直徑,m;Ae為等價流通截面,m2;下標(biāo)i表示第i段控制體。
1.2.2 止回閥模型
將這兩個數(shù)組放在軟件數(shù)據(jù)顯示功能模塊中的組織數(shù)據(jù)步驟里,運行效果見圖5。結(jié)果表明數(shù)據(jù)顯示功能模塊能夠?qū)⒔邮艿降臄?shù)據(jù)正常地顯示在布局中。
止回閥穩(wěn)態(tài)與瞬態(tài)兩種狀態(tài)進行模擬。1.穩(wěn)態(tài)情況,壓降公式為:
2.瞬態(tài)情況,采用以下公式:
式中:Cd為系數(shù),不同時間下,補償系數(shù)Cd采用不同的計算方法;A0為管道截面積,m2;ΔP為閥門入口和出口之間的壓降,Pa。
1.2.3 阻力件模型
阻力件一般采用如下公式計算:
式中:ΔP0為阻力件前后壓降,Pa;K為損耗系數(shù);A0為過水?dāng)嗝婷娣e,m2;Q為體積流量,m3/s。
利用建立的雙泵并聯(lián)系統(tǒng)模型,模擬1#運行泵切換到2#備用泵的瞬態(tài)過程,模擬結(jié)果與試驗結(jié)果比較如圖2所示。從圖2可以看出試驗結(jié)果與模擬結(jié)果的偏差在誤差允許范圍內(nèi),因此,本模型是合理的。
在雙泵切換過程中,系統(tǒng)流量曲線近似呈W形。圖2中,在AB段,1#離心泵由開啟狀態(tài)變?yōu)殛P(guān)閉狀態(tài),支路流量減少,但是由于離心泵惰轉(zhuǎn)作用,支路仍然會產(chǎn)生正向流量,此時2#離心泵處于關(guān)閉狀態(tài),總回路流量下降;在BC段,1#泵惰轉(zhuǎn)產(chǎn)生正向流量,2#泵由關(guān)閉轉(zhuǎn)為開啟,2#泵支路流量增加,導(dǎo)致整個回路流量增加,到達(dá)C點時,1#泵支路正向流量為零,整個回路流量達(dá)到局部峰值;在CD段,1#泵支路出現(xiàn)倒流,2#泵支路正流,主回路系統(tǒng)總流量下降;在DE段,1#止回閥逐漸關(guān)閉,倒流減小,2#泵流量逐步增大,系統(tǒng)總流量增加,到達(dá)E點時1#止回閥全部關(guān)閉,系統(tǒng)重新達(dá)到穩(wěn)定狀態(tài)。
圖2 雙泵切換流量特性曲線Fig.2 Flow curves in the process of switching pumps
雙泵切換方式包括以下三種:
①先啟動備用泵后關(guān)閉運行泵(方式一);②先關(guān)閉運行泵后啟動備用泵(方式二);③同時開啟備用泵和關(guān)閉運行泵(方式三)。
實際工程中不同工況下對應(yīng)不同的泵切換方式:在運行泵發(fā)生故障停機,備用泵啟用時,雙泵切換采用方式一,當(dāng)時間間隔很小時也可視為方式三;在不停止供水條件下檢修某一臺給水泵或者均衡兩個給水泵壽命等計劃內(nèi)切換時,可根據(jù)需要采用三種啟動方式的任意一種。
利用已經(jīng)驗證的系統(tǒng)模型對三種不同泵切換方式下系統(tǒng)流量波動情況進行模擬,結(jié)果如圖3所示。
圖3 不同泵切換方式下系統(tǒng)流量波動曲線Fig.3 Flow curves in different switching mode
(1)先啟后閉,在AB段,備用泵首先啟動,此時運行泵(1#)滿功率運行,備用泵(2#)流量增加,1#、2#支路正流,給水系統(tǒng)流量增加,流量超過正常供給值,到達(dá)B點時系統(tǒng)流量最大;此后1#泵關(guān)閉,1#支路流量減少,但是1#泵在慣性惰轉(zhuǎn)作用下繼續(xù)正向給水,2#泵流量繼續(xù)增加,系統(tǒng)流量減小;之后1#泵停止,所在支路出現(xiàn)逆流,2#泵支路正流,直到C點系統(tǒng)流量達(dá)到最??;在CD段,1#止回閥逐漸關(guān)閉,倒流減小,系統(tǒng)總流量增加,到達(dá)D點時1#止回閥全部關(guān)閉,系統(tǒng)流量達(dá)到穩(wěn)定。由于采用先開啟后關(guān)閉的切換方式,系統(tǒng)流量在波動階段前期大于穩(wěn)定輸出流量值,切換過程流量輸出總和大于穩(wěn)定輸出時的流量總和。
(2)先閉后啟與同時啟閉切換方式,雙泵切換過程中,系統(tǒng)流量特性曲線近似成W型,其形成原理見前文1.3節(jié)中W型曲線形成機理分析。
對于計劃內(nèi)雙泵切換,三種切換方式均可采用。方式一中系統(tǒng)失水量小,最低流量較高,系統(tǒng)穩(wěn)定時間短,但流量波動較大,在切換瞬態(tài)過程中系統(tǒng)總流量大于穩(wěn)定輸出時流量。與方式一相比,方式二、三在雙泵切換過程中失水量較大,系統(tǒng)總流量低于穩(wěn)定輸出時供水量。
對于非計劃內(nèi)事故工況雙泵切換,可采用方式二、三,即運行泵先停,備用泵迅速啟動。比較發(fā)現(xiàn)方式三在雙泵切換過程中系統(tǒng)失水量小,最低流量較高,系統(tǒng)恢復(fù)穩(wěn)定時間短。在實際工況中,采用同時啟閉方式,即盡量減少雙泵切換時間間隔將更易于緩解流量波動問題。
不同的切換方式具有不同的流量特性曲線,對于切換方式一和方式二,切換過程時間間隔也發(fā)揮著明顯的作用。通過模擬不同時間間隔下系統(tǒng)流量的波動,比較時間間隔大小對系統(tǒng)流量特性的影響。
先啟后閉雙泵切換方式在不同時間間隔下系統(tǒng)流量波動如圖4所示,切換時間間隔為分別0.5s、1s和1.5s。圖4表明,對于不同時間間隔,切換過程中流量曲線形狀相近,AB段內(nèi)2#備用泵開啟,流量曲線相重合;時間間隔大小不同,系統(tǒng)流量峰值出現(xiàn)的時間和位置也不同,即間隔時間越長,峰值出現(xiàn)越晚、位置也越高,總體上系統(tǒng)流量波動也將愈加明顯。不同時間間隔直接影響1#運行泵關(guān)閉時刻,時間間隔越長,運行泵關(guān)閉延遲越長,具體表現(xiàn)為流量曲線峰值位置延時并且上移,造成后續(xù)的流量波動階段以及流量恢復(fù)穩(wěn)定階段流量波動更加明顯。因此,對于先啟后閉切換方式,切換間隔時間越長,系統(tǒng)流量波動也更大,給水穩(wěn)定性受到更大的影響。
圖4 先啟后閉方式切換下系統(tǒng)流量波動Fig.4 Flow curves in the switching mode of closing one pump after starting another one
先閉后啟雙泵切換方式在不同時間間隔下系統(tǒng)流量波動如圖5所示,切換時間間隔分別為0.5s、1s和1.5s。由圖5可知,對于三種不同時間間隔切換方式,切換過程中流量曲線形狀近似,AB段內(nèi)1#運行泵關(guān)閉,流量曲線相重合;間隔時間長短影響了系統(tǒng)流量波谷出現(xiàn)的時間和位置,間隔時間越長,系統(tǒng)流量波谷出現(xiàn)越晚、出現(xiàn)位置越低,總體上系統(tǒng)流量波動也將愈加明顯。不同時間間隔直接影響2#備用泵開啟時刻,時間間隔越長,備用泵開啟延遲越長,具體表現(xiàn)為流量曲線波谷位置下移,造成后續(xù)的流量波動階段以及流量恢復(fù)穩(wěn)定階段流量波動更加明顯。因此,對于先閉后啟切換方式,切換間隔時間越長,系統(tǒng)失水量越大,流量波動也更大,給水穩(wěn)定性受到更大的影響。
圖5 先閉后啟泵切換方式下系統(tǒng)流量波動圖Fig.5 Flow curves in the switching mode of starting one pump after closing another one
綜上,無論采取哪種切換方式,縮短泵切換時間間隔,能減小系統(tǒng)流量波動,使系統(tǒng)在較短時間內(nèi)恢復(fù)穩(wěn)定工作。
對雙泵并聯(lián)給水系統(tǒng)進行了建模,通過試驗證明了模型的準(zhǔn)確性。模擬不同雙泵切換方式條件下系統(tǒng)流量波動過程,得到以下結(jié)論:(1)模擬結(jié)果與實驗結(jié)果偏差在允許的范圍內(nèi);(2)先啟后閉切換方式下系統(tǒng)失水量小,最低流量大,同時啟閉切換方式系統(tǒng)流量波動最小,系統(tǒng)恢復(fù)穩(wěn)定時間短;(3)采用同時啟閉方式在泵切換過程中系統(tǒng)失水量小,最低流量較高,系統(tǒng)恢復(fù)穩(wěn)定時間短;(4)雙泵切換時間間隔越短,系統(tǒng)流量波動越小,系統(tǒng)穩(wěn)定時間越短,在實際工況中,要盡量減小切換時間間隔。
[1] 郭強,趙新文,蔡琦.反應(yīng)堆冗余泵組的共因失效分析[J].核動力工程,2009,30(5):89-92.
[2] Wenxi Tian,GH Su,Gaopeng Wang,et al.Numerical simulation and optimization on valve-induced water hammer characteristics for parallel pump feedwater system[J].Annals of Nuclear Energy,2008,35(12):2280-2287.
[3] TIAN Wenxi,SU Guanghui,WANG Gaopeng,et al.Mitigating check valve slamming and subsequentwater hammer events for PPFS using MOC[J].Nuclear Science and Techniques,2009,20(2):118-123.
[4] 吳大轉(zhuǎn),許斌杰,李志峰,等.離心泵瞬態(tài)操作條件下內(nèi)部流動的數(shù)值模擬[J].工程熱物理學(xué)報,2009,30(5):781-783.
[5] 楊洪波,吳煥云.核電廠常規(guī)島電動給水泵切換瞬態(tài)分析[J].核電技術(shù)論壇熱力發(fā)電二〇二,2011,40(1):82-83.
[6] 張龍飛,張大發(fā),王少明.轉(zhuǎn)動慣量對船用核動力主泵瞬態(tài)特性的影響研究[J].船海工程,2005,(2):55-57.
[7] Kolev N,Petrov N,Ivanov B,et al.Simulation of the VVER-1000pump start-up experiment of the OECD V1000CT benchmark with CATHARE and TRAC-PF1[J].Progress in Nuclear Energy,2006(48):922-936.
[8] Rohani M,Afshar M H.Simulation of transient flow caused by pump failure:Point-Implicit Method of Characteristics[J].Annals of Nuclear Energy,2010(37):1742-1750.
[9] 郭玉君,張金玲,秋穗正,等.反應(yīng)堆系統(tǒng)冷卻劑泵流量特性計算模型[J].核科學(xué)與工程,2005,15(3):220-231.
[10] Kazem Farhadi,Anis Bousbia-salah,F(xiàn)ranscesco D’Auria.A model for the analysis of pump start-up transients in Tehran Research Reactor[J].Progress in Nuclear Energy,2007(49):499-510.