何金橋，肖波，湯盛萍，王運(yùn)民，曹雄，孫振宇,

(1. 長沙理工大學(xué) 能源與動(dòng)力工程學(xué)院，湖南 長沙，410004；2. 中國廣東核電集團(tuán)有限公司，廣東 深圳，518031)

汽輪機(jī)冷端是指從汽輪機(jī)末級葉片出口→冷凝器→循環(huán)水系統(tǒng)范圍內(nèi)構(gòu)成熱力循環(huán)有效冷源的設(shè)備及其管路系統(tǒng)，其中的主要設(shè)備有凝汽器、循環(huán)水泵、真空泵及冷卻塔等[1]。相關(guān)研究表明，冷端處汽輪機(jī)背壓升高1 kPa，相應(yīng)發(fā)電熱耗率增加0.7%～ 0.9%[2]。電站汽輪機(jī)組凝汽器殼側(cè)的蒸汽流動(dòng)特性及凝結(jié)換熱特性是決定冷端各設(shè)備運(yùn)行工況的重要因素。隨著大功率汽輪機(jī)組的發(fā)展，凝汽器結(jié)構(gòu)越來越復(fù)雜，汽側(cè)乏汽凝結(jié)換熱時(shí)的流動(dòng)特性也更加復(fù)雜[3-5]，使得冷端汽輪機(jī)背壓在變工況負(fù)荷下發(fā)生大范圍波動(dòng)，這將直接影響到整個(gè)火力發(fā)電機(jī)組的熱經(jīng)濟(jì)性和安全性[1,6]。多孔介質(zhì)模型將管子看作多孔骨架，降低了網(wǎng)格生成的復(fù)雜度并大大減少了網(wǎng)格數(shù)量和計(jì)算負(fù)荷，有效反映了多孔介質(zhì)換熱時(shí)的流動(dòng)特性。在目前管殼式換熱器數(shù)值研究中，一般采用Patankar等[7-8]提出的分布阻力方法，對管殼式換熱器建立二維或三維的數(shù)值模型并求解[9-12]。我國學(xué)者在應(yīng)用多孔介質(zhì)模型進(jìn)行數(shù)值模擬方面進(jìn)行了探索性研究，但對管束區(qū)換熱時(shí)流動(dòng)特性研究不多[13-14]，這為研究結(jié)構(gòu)復(fù)雜、管數(shù)眾多的凝汽器管束區(qū)換熱時(shí)的流動(dòng)特性創(chuàng)造了有利條件。

1 火電廠凝汽器數(shù)學(xué)模型分析

1.1 凝汽器結(jié)構(gòu)及工作特性分析

株洲電廠310 MW火電機(jī)組采用的凝汽器為單背壓、單殼體、對分式、雙流程、表面式、冷凝管束橫向布置的凝汽器。

該種形式的凝汽器因管內(nèi)冷卻水流阻小，不易堵塞，易加工和大型化，故得到廣泛應(yīng)用。但是，殼側(cè)內(nèi)乏汽因復(fù)雜的布管導(dǎo)致氣流進(jìn)入管束后會(huì)產(chǎn)生較大波動(dòng)，這將使汽輪機(jī)末級排汽背壓波動(dòng)，導(dǎo)致凝汽式汽輪機(jī)機(jī)組冷端運(yùn)行工況劣化。

圖1 株洲電廠310 MW機(jī)組凝汽器管束布置Fig.1 Condenser tube bundle distribution of 310 MW unit in Zhuzhou Power Plant

圖1所示為該型凝汽器布置管束的管板分區(qū)結(jié)構(gòu)簡圖。由圖1可知：該凝汽器不同管束區(qū)結(jié)構(gòu)存在較大差異，具有分型特性[9-13]，布管數(shù)量并不均勻，相應(yīng)區(qū)間的孔隙率也不一致。因此，應(yīng)根據(jù)分型多孔介質(zhì)模型的特點(diǎn)，將管束區(qū)分為不同孔隙率的3個(gè)區(qū)域，即上管束區(qū)、下管束區(qū)和空冷區(qū)。不同管束區(qū)域的結(jié)構(gòu)特點(diǎn)如表1所示。

表1 管板結(jié)構(gòu)特點(diǎn)Table 1 Structure characteristics of tube plate

上、下管束區(qū)通過凝汽器殼體外部的彎頭連接形成雙流程的冷卻效果，管束區(qū)中部的空冷區(qū)在管束邊緣的上、下管束區(qū)之間存在明顯空隙，這會(huì)使得一部分乏汽直接進(jìn)入到空冷區(qū)，因此，空冷區(qū)中也布置一定數(shù)量的冷卻水管，以進(jìn)一步提高凝汽器的凝結(jié)效率。

1.2 管束區(qū)乏汽換熱特性分析

該凝汽器采用水平管外凝結(jié)的工作形式，即在凝汽器內(nèi)部，冷卻水在循環(huán)水泵的驅(qū)動(dòng)下連續(xù)不斷地流過冷卻管，乏汽不斷地在冷卻管外壁凝結(jié)。由于乏汽溫度不高，一般不考慮凝汽器與外界的散熱，故凝汽器的傳熱系數(shù)主要受凝汽器水側(cè)污垢、汽側(cè)凝結(jié)液膜厚度、管子的排列、蒸汽中空氣含量等影響。顯然，凝結(jié)液膜的厚度主要取決于材料的表面特性，水側(cè)污垢通過定期清洗也可以消除其對換熱的影響，所以，影響凝汽器管束區(qū)不同位置換熱性能的因素在于管子的排列和乏汽中空氣的含量。凝汽器內(nèi)少量從外界漏入的不凝結(jié)氣體隨著乏汽的凝結(jié)，其濃度逐漸增大，這將導(dǎo)致乏汽的分壓力降低，難以凝結(jié)在冷卻管表面，對凝結(jié)換熱過程產(chǎn)生重要影響[2,10]，但不凝結(jié)氣體含量小，因此，主要影響空冷區(qū)的換熱凝結(jié)。

1.3 多孔介質(zhì)模型分析

凝汽器殼側(cè)蒸汽流動(dòng)是發(fā)生在冷卻管束內(nèi)外的多組分同時(shí)伴有相變的復(fù)雜過程。合理研究凝汽器內(nèi)發(fā)起流動(dòng)特性需對該凝汽器內(nèi)部的流動(dòng)與傳熱過程進(jìn)行以下幾方面簡化和假設(shè)：忽略溫度及液相體積的影響；空氣泵具有足夠抽除空氣和未凝結(jié)蒸汽的能力；凝汽器外殼絕熱，蒸汽凝結(jié)釋放的熱量全部傳給冷卻水；忽略管板、支撐件及抽氣導(dǎo)出管等的影響，各汽室的流動(dòng)可近似看成是乏汽橫掠管束的二維定常流動(dòng)。

此外，對凝汽器的工作條件也進(jìn)行如下假設(shè)：蒸汽在飽和狀態(tài)下凝結(jié)；冷卻水流量在同一流程的各冷卻管中始終相等；忽略不凝結(jié)氣體對換熱的影響，凝汽器內(nèi)的空氣和蒸汽均勻混合；殼側(cè)管束外各處壓力相等。

根據(jù)上述假設(shè)，影響汽輪機(jī)末級排汽壓力的因素主要集中在凝汽器管束區(qū)，因此，只需將管束區(qū)當(dāng)作多孔介質(zhì)區(qū)域，將蒸汽在其中的流動(dòng)看作是在多孔介質(zhì)內(nèi)的流動(dòng)，就可以建立具有分布阻力和分布質(zhì)量匯的多孔介質(zhì)模型。

質(zhì)量連續(xù)性方程：

X方向動(dòng)量守恒方程為

Y方向動(dòng)量守恒方程為

式中：ε為孔隙率，主要取決于管束區(qū)布管特性；μ′為慣性阻力系數(shù)；eμ為當(dāng)量黏性系數(shù)；fx和fy為冷卻管束間的摩擦阻力系數(shù)，均取決于凝汽器負(fù)荷特性；Q，Qu和Qv為源項(xiàng)。

2 凝汽器管束區(qū)乏汽流場特性分析

2.1 流動(dòng)特性分析

圖2所示為300 MW負(fù)荷下管束區(qū)乏汽的速度分布矢量圖。由圖2可知：下管束區(qū)氣流速度分布較均勻，且較大；上管束區(qū)頂部氣流速度較大，下部氣流速度較低，且低于空冷區(qū)及下管束區(qū)中大部分區(qū)域的氣流速度；空冷區(qū)中僅有部分氣流流向不凝結(jié)氣體排氣口，大量的氣流又回流到上、下冷凝管束區(qū)。其中，回流到下管束區(qū)中乏汽氣流速較大且較均勻；回流到上管束中的乏汽流速較小且主要回流到上管束區(qū)側(cè)向區(qū)域，頂部區(qū)域回流量較小。

經(jīng)測試，上、下管束區(qū)之間間隙的垂直距離為600 mm，空冷區(qū)經(jīng)該間隙直接與凝汽器殼體中的乏汽聯(lián)通，必然導(dǎo)致大量乏汽直接流入空冷區(qū)，加之空冷區(qū)中冷卻管數(shù)量少且主要集中在排氣口周圍及其頂部軸線上，因此，乏汽一旦流入空冷區(qū)，其流動(dòng)方向?qū)⒅饕Q于排氣口初始壓力及上、下冷卻管束的冷凝作用所導(dǎo)致的真空效應(yīng)。圖2所示的速度場表明：空冷區(qū)中的乏汽大量流入到下管束區(qū)，而只有一小部分進(jìn)入空冷區(qū)，使得下管束區(qū)形成多面進(jìn)氣，因此，其中氣流速度分布較均勻；而上管束雖然在與空冷區(qū)的界面處也有乏汽流入，但由于排氣口置于下管束區(qū)中，且上管束與空冷區(qū)之間的內(nèi)邊界較狹長，因此，上管束區(qū)中間除頂部外的管束中基本上沒有氣流，也難以達(dá)到較好的冷凝作用。

圖2 300 MW負(fù)荷下管束區(qū)乏汽速度分布Fig.2 Velocity of flow distribution of tube bundle under load of 300 MW

2.2 流量特性分析

圖3所示為300 MW負(fù)荷下管束區(qū)乏汽的流量分布等值線圖。乏汽在3個(gè)不同管束區(qū)中的流量存在顯著差異，上管束區(qū)中流量極大值出現(xiàn)在頂部轉(zhuǎn)彎處，即圖1所示的D和E處；下管束區(qū)的則分布在其與空冷區(qū)的界面處，即圖1中所示的N和H處?？绽鋮^(qū)中靠近排氣口處的氣流存在顯著的回流現(xiàn)象，這使得排氣口處的乏汽量大大減少，有利于減少乏汽中可用工質(zhì)的損失，但是，由于汽輪機(jī)末級工質(zhì)已經(jīng)大大低于環(huán)境大氣的壓力，必然會(huì)有極少量的不凝性空氣泄露進(jìn)入到乏汽中，若不能及時(shí)排除，則會(huì)影響冷凝管束外側(cè)的換熱。

上管束區(qū)中部的乏汽流量等值線表明其中乏汽流量較低，但變化梯度不大，這主要由冷凝管束的冷凝作用及空冷區(qū)中的乏汽熱浮升作用共同導(dǎo)致。下管束底部中間位置處的流量等值線表明該處乏汽流量較少而且變化梯度較大，這說明乏汽自汽輪機(jī)末級排入到凝汽器后，經(jīng)過凝汽器管束區(qū)上部不同位置的冷凝和流動(dòng)，到達(dá)底部中間位置處的乏汽量很少，其中冷凝管束難以起到冷凝作用。

圖3 300 MW負(fù)荷下的管束區(qū)乏汽流量分布Fig.3 Rate of flow distribution of tube bundle under load of 300 MW

圖4所示為不同負(fù)荷下不同管束區(qū)中的最大流量分布情況。由圖4可知：下管束區(qū)中的最大流量最大，空冷區(qū)中的最大流量雖然最小，三者隨負(fù)荷的增大基本上近似平行增大，只是在由 300 MW 增加到 310 MW的負(fù)荷時(shí)，個(gè)管束區(qū)中的最大流量變化較大。由此可見：導(dǎo)致該種結(jié)構(gòu)形式的凝汽器冷凝管束的換熱不均的首要原因應(yīng)在于管板管束區(qū)內(nèi)邊和外邊緣的進(jìn)汽不均勻，內(nèi)邊緣需通過空冷區(qū)到達(dá)主凝結(jié)區(qū)凝結(jié)換熱，這將增加流動(dòng)阻力，導(dǎo)致流量減少；此外，管束底仍有部分區(qū)域氣流難以到達(dá)，也降低了其冷凝作用。

圖4 不同負(fù)荷下的管束區(qū)中最大乏汽流量分布Fig.4 The maximum rate of flow distribution of tube bundle under different loads

2.3 壓力分布特性

圖5所示為300 MW負(fù)荷下管束區(qū)的乏汽動(dòng)壓分布等值線圖。由圖5可知：空冷區(qū)中存在一對稱分布的高動(dòng)壓區(qū)。該動(dòng)壓區(qū)緊靠下管束的垂直內(nèi)邊界。結(jié)合圖2所示乏汽速度分布可知：正是由于該動(dòng)壓區(qū)才使得流入空冷區(qū)中的乏汽轉(zhuǎn)向流入到凝結(jié)管束區(qū)中，避免了工質(zhì)直接通過排汽口排放的損失。但是，由于負(fù)荷工況的變化，該動(dòng)壓區(qū)所處的位置以及動(dòng)壓均將發(fā)生變化，若在確定的排汽壓力下，排汽口排出氣體中的可凝結(jié)性工質(zhì)的含量也隨之發(fā)生改變。圖5中，管束區(qū)中的壓力為對應(yīng)負(fù)荷下的最大計(jì)算值；而非管束區(qū)中的乏汽壓力則為實(shí)際測量值。

圖5 300 MW負(fù)荷下的管束區(qū)乏汽壓力分布Fig.5 Steam exhaust dynamic pressure distribution of tube bundle under load of 300 MW

圖5還表明：在整個(gè)管束區(qū)中，位于管板4個(gè)轉(zhuǎn)角及外邊界處的動(dòng)壓較大，而且等值線之間的間距較小，說明動(dòng)壓減少的梯度?。簧?、下管束區(qū)中其他區(qū)域動(dòng)壓等值線之間的間距隨之增大，相應(yīng)地，動(dòng)壓減小得更快，即流動(dòng)阻力增加，這將導(dǎo)致乏汽主要通過上、下管束區(qū)外邊界附近的冷卻管換熱凝結(jié)。

圖6所示為不同負(fù)荷下，不同管束區(qū)外邊界上最大動(dòng)壓。由圖6可知：上、下2管束區(qū)外邊界上的動(dòng)壓變化基本一致，而且均在大于額定負(fù)荷后出現(xiàn)顯著增大；下管束區(qū)中的最大動(dòng)壓大，說明下管束的凝結(jié)換熱效果更好?？绽鋮^(qū)外邊界為圖1所示的CBN及FGH邊所在位置，其上的最大動(dòng)壓在240 MW附近也出現(xiàn)了較大波動(dòng)。由此可見：在240 MW附近的負(fù)荷下，若不及時(shí)調(diào)整排汽口壓力，將明顯導(dǎo)致不良運(yùn)行工況出現(xiàn)。

圖6 不同負(fù)荷下管束區(qū)入口乏汽動(dòng)壓分布Fig.6 Steam exhaust dynamic pressure distribution of inlet of tube bundle under different loads

3 結(jié)論

(1) 在均勻管內(nèi)換熱條件下，該結(jié)構(gòu)類型凝汽器中的凝結(jié)效果主要取決于換熱管束的外邊界附近的冷卻管束的換熱，而且上管束區(qū)域的冷凝效果比下管束區(qū)的冷凝效果差。

(2) 導(dǎo)致該種結(jié)構(gòu)形式的凝汽器冷凝管束換熱不均的首要原因在于管板管束區(qū)內(nèi)邊和外邊緣的進(jìn)汽不均勻，內(nèi)邊緣需通過空冷區(qū)到達(dá)主凝結(jié)區(qū)凝結(jié)換熱，這將增加流動(dòng)阻力，導(dǎo)致流量減少；此外，管束底仍有部分區(qū)域氣流難以到達(dá)，也降低了其冷凝作用。

(3) 空冷區(qū)中存在一對稱分布的高動(dòng)壓區(qū)。該動(dòng)壓區(qū)緊靠下管束的垂直內(nèi)邊界。該動(dòng)壓區(qū)使得流入空冷區(qū)中的乏汽轉(zhuǎn)向流入到凝結(jié)管束區(qū)中，一方面，避免了工質(zhì)直接通過排汽口排放的損失；另一方面，該動(dòng)壓區(qū)所處的位置以及動(dòng)壓隨負(fù)荷均將發(fā)生變化，尤其在240 MW負(fù)荷附近工況下變化顯著，這將導(dǎo)致排氣口帶走更多的工質(zhì)。

[1] 徐仿優(yōu), 胡光明, 劉光耀, 等．沙角 A電廠 300 MW機(jī)組冷端系統(tǒng)經(jīng)濟(jì)性診斷及優(yōu)化運(yùn)行[J]. 熱力發(fā)電, 2010, 31(7):7-10.

XU Fangyou, HU Guangming, LIU Guangyao, et al. Operation optimization and economic efficiency diagnosis on cold end system of 300 MW unit Shajiao A power plant[J]. Thermal Power Generation, 2010, 31(7): 7-10.

[2] 王學(xué)棟, 袁建華, 曲建麗, 等. 330 MW機(jī)組凝汽器改造的數(shù)值模擬及性能分析[J]. 汽輪機(jī)技術(shù), 2010, 52(4): 150-154.

WANG Xuedong, YUAN Jianhua, QU Jianli, et al. Numerical simulation and performance analysis of condenser retrofit of 330 MW unit[J]. Turbine Technology, 2010, 52(4): 150-154.

[3] 周蘭欣, 李富云, 李衛(wèi)華. 凝汽器殼側(cè)準(zhǔn)三維數(shù)值研究[J]. 中國電機(jī)工程學(xué)報(bào), 2008, 28(23): 25-30.

ZHOU Lanxin, LI Fuyun, LI Weihua. Quasi-three-dimensional numerical study of shell side of condenser[J]. Proceedings of the CSEE, 2008, 28(23): 25-30.

[4] Zhu Q Y, Li Y. Numerical simulation of the transient heat and liquid moisture transfer through porous textiles with consideration of electric double layer[J]. Int J Heat Mass Transfer, 2010, 53(7): 1417-1425.

[5] 俞茂錚, 姚秀平, 汪國山. 大功率汽輪機(jī)凝汽器汽相流動(dòng)與傳熱特性的數(shù)值模擬[J]. 動(dòng)力工程, 1995, 15(6): 42-49.

YU Maozheng, YAO Xiuping, WANG Guoshan. Numerical analysis of the steam flow behavior and heat transfer performance of condensers of large steam turbines[J]. Power Engineering, 1995, 15(6): 42-49.

[6] 張卓澄. 大型電站凝汽器[M]. 北京: 機(jī)械工業(yè)出版社, 1993:4-5.

ZHANG zhuocheng. Large power station condenser[M]. Beijing:China Machine Press, 1993: 4-5.

[7] Spalding D B. The calculation of free-convertion phenomenon in gas-liquid mixtures[C]//ICHMT Semonar. Dubrovnil, 1967:1-15.

[8] Patankar S V, Spalding D B. Computer analysis of threedimensional flow and heat transfer in a steam generator[C]//Forsch. Ingenieurwes, 1978: 4.

[9] Zhu Q Y, Xie M H, Yang J, et al. Investigation of the 3D model of coupled heat and mass moisture transfer in hygroscopic porous fi brous media[J]. Int J Heat Mass Transfer, 2010, 53(8):3914-3927.

[10] KOU Jianlong, LIU Yang, WU Fengmin, et al. Fractal analysis of effective thermal conductivity for three-phase (unsaturated)porous media[J]. Journal of Applied Physics, 2009(5): 106.

[11] Ono Y, Mayama H, Furo I, et al. Characterization and structural investigation of fractal porous-silica over an extremely wide scale range of pore size[J]. J Colloid Interf Sci, 2009, 336:215-225.

[12] Wu J S, Yu B M. A fractal resistance model for fl ow through porous media[J]. Int J Heat Mass Transfer, 2007, 50(8):3925-3932.

[13] Zhu Q Y, Xie M H, Yang J, et al. A fractal model for the coupled heat and mass transfer in porous fibrous media[J]. Int J Heat Mass Transfer, 2011, 54(7): 1400-1409.

[14] 曾輝, 孟繼安, 李志信. N-17000型凝汽器的熱力性能的數(shù)值分析[J]. 工程熱物理學(xué)報(bào), 2010, 31(7): 1213-1214．

ZENG Hui, MENG Jian, LI Zhixin. Numerical analsis on the thermal performan[J]. Journal of Engineering Thermal Physics,2010, 31(7): 1212-1214.