王巍, 劉曉薇, 孫偉文, 王曉放
(大連理工大學 能源與動力學院,遼寧 大連 116024)
在我國以煤炭為主的能源消費模式的背景下,郭烈錦等[1]創新性的提出了超臨界水煤氣化制氫技術,以實現煤的清潔利用。然而,氣化物由CO2+H2O+H2等組成,當其中的氫氣經過氧化形成水蒸氣之后,其混合物主要由H2O+CO2組成,CO2含量較大,且在凝汽器內不發生凝結,將使其凝結傳熱規律不同于純蒸汽,這給混合工質氣輪機最佳排氣壓力的確定以及熱力發電循環系統的設計提出了新的挑戰。
在蒸汽輪機凝汽器真空度的相關研究中,李勇等[2-3]提出了考慮節水因素、水側管壁清潔程度、抽氣設備工作性能的凝汽器最佳真空的確定方法,實現電廠綜合效益最大化。王保力[4]考慮了排汽阻力和鍋爐補充水對排汽壓力的影響,提高了計算最佳排汽壓力的精度。劉克為等[5]整體考慮凝汽器特定管束形式及喉部流場的分布探索了降低凝汽器汽側阻力合理有效手段。趙宏偉等[6]提出基于實時數據的改進型最佳真空計算模型,計算最佳循環水流量和對應的發電凈增功率,降低了電廠煤耗。李建強等[7]基于PSO-SVR模型預測凝汽器真空目標值,對不同負荷下凝汽器真空目標值進行了預測,并通過凝汽器真空實測值與模型預測目標值對比,實現對凝汽器真空運行狀態的評估。Anozie等[8]利用MS Excel宏指令對冷卻水流量和排汽壓力進行優化,使得循環效率提升2%,燃料節約3.8%。Attia[9]研究了冷卻水初溫對某核電站熱力性能的影響,發現冷卻水入口溫度上升1 ℃會使電站輸出功降低0.444%、熱效率降低0.152%,為核電站選址提供參考。Reguckia等[10]對電站冷卻系統的冷卻水流量進行研究,結果表明對冷卻水量分布的掌控可提升凝汽器性能。
可見,蒸汽輪機最佳排汽壓力的研究比較成熟。然而,由于超臨界水煤氣化制氫工藝的獨創性,國內外幾乎沒有關于H2O+CO2混合工質氣輪機最佳排氣壓力和凝汽器內凝結過程方面的研究成果;另一方面,在混合工質循環系統中混合工質氣輪機排氣壓力的選取時,未考慮大量不凝氣體CO2的分壓影響,仍參考傳統水蒸汽汽輪機排氣壓力的常用取值,使得所取排氣壓力偏低,同時受到冷卻水溫所限,導致混合工質在凝汽器內凝結水量不足,為凝汽器設計帶來問題。
因此,本文基于H2O+CO2混合工質組分實驗數據,以50 MW混合工質氣輪機系統凈增功為目標函數,分析了混合工質氣輪機所適用的最佳排氣壓力取值方法。為方便區分,本文中水蒸汽工質一律使用“汽輪機”,混合工質一律使用“氣輪機”。無論是水蒸汽還是混合工質,均只有水蒸汽部分發生凝結,因此一律使用 “凝汽器”。
在超臨界水煤氣化制氫工藝中,水煤漿在氣化爐超臨界水中完成水解、熱解等一系列反應,氣化爐出口為H2、CH4等可燃氣體。在目前情況下,為保證運行的安全,氣化爐出口含氫混合氣體不能直接進入氣輪機做功,需要在氧化爐內進行溫和氧化過程后進入氣輪機,同時,實現氧化爐放熱量、氣化爐吸熱量、超臨界水加熱量等系統熱量的平衡。
表1為煤流量1 kg/s情況下氣化爐出口混合氣體各組分含量的實驗數據,經過熱平衡分析計算,氣輪機入口混合工質中CO2摩爾分數為8.47%,質量分數為18.45%。

表1 氣化爐出口和氣輪機入口各組分含量(以煤流量1 kg/s為準)
基于ASPEN PLUS化工流程軟件,開展實驗工況熱平衡模擬計算,表2為氣化爐吸熱量、超臨界水吸熱量和氧化爐放熱量的模擬結果與實驗測量結果的對比。氣化爐吸熱量和氧化爐放熱量與實驗數據吻合,發電系統實現熱量平衡。建立圖1所示的循環系統圖,氣輪機熱力系統采用四級不調整抽氣,回熱加熱器出口疏水逐級自流進入凝汽器的熱井。若不考慮制氧裝置耗能和CO2捕集耗能,為實現系統熱量自給自足,不外加鍋爐而引入新的污染,需要將氧化爐的放熱量提供給氣化反應爐,并保證給水終溫加熱到318 ℃,以實現發電系統凈輸出功50 MW,由于混合工質汽輪機還未運用到工程實際中,熱力結構的設計尚未完成,因此假定汽輪機內效率0.9,計算得系統效率43.72%,此時混合工質汽機排汽流量為34.78 kg/s,排汽干度為0.882,循環各個物流溫度壓力、干度等值均已在圖中標出。

表2 系統主要部分吸/放熱量
由于本文的H2O+CO2混合工質氣輪機研究中所涉及的排氣壓力為5~53 kPa,排氣溫度為23~80 ℃,在此范圍內CO2的露點溫度非常低,不會出現CO2氣體的凝結,因此本文將其作為不凝氣體處理,圖2展示了SRK方程和PR方程計算出的不同壓力所對應的CO2露點溫度數據。
凝汽器內部壓力是蒸汽工質溫度降低,造成蒸汽凝結后形成的真空。在圖1所示的H2O+CO2混合工質氣輪機電站系統中,最佳氣輪機排氣壓力受許多因素的影響。在本研究中,將以系統凈增功為目標函數,通過增加循環水量實現氣輪機排氣壓力的降低,最終提出氣輪機最佳排氣壓力值。

圖1 基于實驗工況設計的熱力發電系統

圖2 5~50 kPa壓力下CO2露點溫度
將氣輪機發電系統相比于設計工況下功率的凈增功ΔW的值作為衡量標準,確定混合工質氣輪機最佳排氣壓力,其計算公式為:
ΔW=ΔWT-ΔWpu
(1)
ΔWT=Q·(hs1-hsc)·ηoi
(2)
ΔWpu=Δqw·(pwi-p0)·vm·ηpu
(3)
式中:ΔWT為氣輪機輸出功增量,kW;Q為氣輪機進氣質量流量,kg/s;hs1為對應某一排氣壓力下氣輪機等熵焓降,kJ/kg;hsc為參考排氣壓力下氣輪機等熵焓降,kJ/kg;ηoi為氣輪機相對內效率;ΔWpu為凝結水泵耗功增量,kW;Δqw為冷卻水流量增量,kg/s;pwi為冷卻水泵出口壓力,kPa;p0為冷卻水泵入口壓力(大氣壓),kPa;vm為冷卻水平均比容,m3/kg;ηpu為冷卻水泵效率。
為保證在排氣壓力變化過程中,氣輪機相對內效率的變化盡可能在小范圍內保持,選取參考排氣壓力值pcc作為比較基準,忽略氣輪機相對內效率的變化。本文在蒸汽輪機最佳排汽壓力確定時,選取參考排氣壓力值pcc=10 kPa。
以系統凈增功為目標函數,降低汽輪機排汽壓力值,最終確定蒸汽輪機最佳排汽壓力值。圖3為相同入口流量的水蒸汽汽輪機發電系統凈增功隨冷卻水流量變化曲線。可以看出,隨著汽輪機排汽壓力逐漸降低,系統凈增功先增大后減小,當冷卻水流量增量為9 170 t/h時,系統凈增功達到最大,此時汽輪機排汽壓力為7 kPa。我國固定式電站凝汽式汽輪機設計背壓一般為4~7 kPa,計算結果符合工程實際[11]。

圖3 純蒸汽汽輪機系統凈增功隨冷卻水流量變化
凝汽器屬于兩相流橫向沖刷管束的管殼式換熱器,并且流向為垂直下降,本文采用Chisholm擬合Baroczy曲線的公式來計算汽阻[12]:
(4)
(5)
(6)
(7)

根據以上方法計算可知,凝汽器汽阻為0.24 kPa,汽阻的存在使得汽輪機排汽壓力提高3.4%。若忽略汽阻對于排汽壓力的影響,可以近似認為蒸汽在凝汽器內為定壓凝結。
同時,若假設H2O+CO2混合工質氣輪機排氣壓力、排氣流量、干度與純蒸汽相同,凝汽器內冷卻水管管徑相同,由于混合工質中含有18.45%的不凝結氣體,一方面,與純蒸汽凝汽器相比,混合工質凝汽器內傳熱系數的大幅減小[11,13],傳熱面積增大,進而使得管束區面積質量流量G降低,摩擦阻力系數fL0減小,另一方面,由于CO2的存在,與純蒸汽工質相比,混合工質氣相密度和氣相粘度增大,液相密度和粘度不變,導致混合工質Chisholm數Y2大于純蒸汽,因此,應用式(4)~(7)計算得出的混合工質在凝汽器內的氣阻均小于純蒸汽汽阻。并且在凝汽器實際運行中,當存在不凝結氣體CO2時,只要保證抽氣泵的能力足夠大, CO2可以快速排出凝汽器,混合工質在凝汽過程中真空壓力可以瞬間建立,此時可以忽略氣阻,基本可以認為混合工質在凝汽器中為定壓凝結,即近似看做凝汽器內工作壓力與氣輪機排氣壓力值相等。
不同于純蒸汽在凝汽器中定溫凝結過程,CO2+H2O混合工質在凝汽器中,隨著水蒸汽凝結量逐漸增加,水蒸汽分壓力逐漸降低,而混合工質的飽和溫度等于蒸汽分壓下水的飽和溫度,因此,混合工質在凝汽器中為降溫凝結過程。
由于存在著CO2分壓力,水蒸汽不可能百分百地凝結,存在小部分蒸汽的不凝結。因此,研究混合工質在凝汽器內部壓力與溫度的關系時,須事先給定蒸汽凝結率。蒸汽凝結率表示混合工質中水蒸汽在凝汽器內的凝結程度,即混合工質凝結完成后,在凝汽器出口,水蒸汽組分中凝結為水的流量占總水蒸汽組分流量的百分數。其計算公式為:
(8)
式中:ξ為蒸汽凝結率;mc為混合工質中水蒸汽部分中凝結為水的質量流量,kg/s;ms為混合工質中水蒸汽部分質量流量,kg/s。凝結率代表了工質水重新進入系統的百分比,凝結率越大,工質回收的水量越多,系統補水量越少。
基于立方形狀態方程RKS,實現了混合工質氣液兩相凝結計算。圖4給出了在凝汽器內部,混合工質壓力、溫度和蒸汽凝結率三者之間的變化規律。針對不同凝汽器內部壓力,較高的蒸汽凝結率對應著較大混合工質溫度的降低,當蒸汽凝結率超過80%以后,混合工質溫度會急劇下降。

圖4 不同壓力下混合工質溫度隨凝結率變化(twi=25 ℃)
然而,凝汽器內混合工質的降溫冷卻過程是通過循環冷卻水換熱實現的,因此,由于環境溫度的限制,混合工質的溫度也不可能無限降低。如圖4所示,當凝汽器內部壓力45 kPa,蒸汽凝結率達到99.4%時,混合工質溫度迅速降低到24.7 ℃,低于環境溫度。因此,當冷卻水入口溫度較高時,很難同時達到較低凝汽器工作壓力和較高蒸汽凝結率。
為方便分析,定義Δt2為凝汽器出口溫差,即凝汽器出口工質與冷卻水入口溫差,如圖5所示。計算公式為:
Δt2=ts2-twi
(9)
式中:ts1為氣輪機排氣溫度,℃;ts2為凝汽器出口工質溫度,℃;two為冷卻水出口溫度,℃;twi為冷卻水入口溫度,℃;Δt1為凝汽器入口端差,℃。當凝汽器出口工質溫度降低,Δt2減小,由于環境溫度限制,Δt2存在最小值,即凝汽器出口端差Δt2 min,混合工質在凝汽器中的換熱過程如圖6。

圖5 混合工質凝汽器進出口溫差示意

圖6 混合工質和冷卻水溫度沿冷卻面積變化
圖7進一步揭示了凝汽器內部混合工質最低壓力的影響因素。分析表明:當冷卻水入口溫度一定時,隨著蒸汽凝結率的提高,凝汽器出口溫差Δt2勢必降低。當蒸汽凝結率繼續增大時,混合工質的出口溫度接近或達到冷卻水入口溫度,造成傳熱面積設計值無限大,這在實際設計中是無法實現的。因此,存在凝汽器出口的最小溫差Δt2 min。若給定Δt2 min=8 ℃,入口水溫twi=25 ℃,則凝汽器出口工質最低溫度ts2=33 ℃。在上述條件限制下,不同的凝汽器內部壓力對應著不同的最高蒸汽凝結率,即等壓線與等溫線交點的橫坐標值。隨著最大蒸汽凝結率的下降,凝汽器壓力降低。而在一定的凝結率要求下,總是可以獲得最小的排氣壓力值Pc min,即等溫線與等壓線交點的縱坐標值,并且隨著蒸汽凝結率的減小,凝汽器內部混合工質的最低壓力降低。
當凝汽器出口端差增大,即等溫線上移,其與等壓線的交點向左上方移動,一定壓力下最大蒸汽凝結率隨之減小,而特定蒸汽凝結率下最低排氣壓力增大;反之,當凝汽器出口端差減小,等溫線下移,其與等壓線的交點向右下方移動,一定壓力下最大蒸汽凝結率隨之增大,而特定蒸汽凝結率下最低排氣壓力降低。

圖7 凝汽器內部混合工質最低壓力的影響因素
蒸汽凝結率的大小反映了重新進入到循環系統中水量的多少,為了減少系統補水量,凝結率不宜過小。然而,由于環境溫度和傳熱面積的限制,使得混合工質中蒸汽的凝結率不可能達到100%。當冷卻水入口水溫25 ℃,凝汽器出口端差Δt2 min=8 ℃時,若蒸汽凝結率在96%~99.5%變化,則凝汽器內混合工質的最低壓力在17~48 kPa變化。這與2.1節中分析的純蒸汽凝汽器內部的最佳壓力相差較大,說明盡管在相同的環境溫度、凝汽器端差等條件下,混合工質H2O+CO2在凝汽器內部不可能建立與蒸汽工質相同的真空壓力,只能獲得大于純蒸汽工質的凝汽器壓力。
圖8給出了當冷卻水入口溫度在15 ℃~25 ℃,出口端差為8 ℃時,混合工質在凝汽器內最低壓力隨冷卻水溫和蒸汽凝結率變化的分析結果。

圖8 最低排氣壓力隨蒸汽凝結率變化(Δt2 min=8 ℃)
在給定蒸汽凝結率時,冷卻水入口溫度越低,凝汽器內混合工質壓力越低;在給定冷卻水入口水溫時,凝汽器內混合工質壓力隨凝結率減小而降低。當冷卻水入口水溫15 ℃,凝結率為96%時,混合工質在凝汽器內部壓力可達到9.3 kPa,這為混合工質氣輪機變工況的研究提供了參考。
混合工質在凝汽器內部的最低壓力是否為氣輪機的最佳排氣壓力,還要取決于系統凈增功的大小,圖9給出了最佳排氣壓力的計算過程。

圖9 最佳排氣壓力計算流程
假定冷卻水入口溫度twi=25 ℃,凝汽器出口端差Δt2 min=8 ℃,選取參考排氣壓力pcc=pcmin+5 kPa,系統凈增功隨冷卻水流量增量Δqw的變化趨勢如圖10所示。結果顯示,在一定蒸汽凝結率下,系統凈增功隨冷卻水流量的增大而增大,是單調遞增的,這與圖2蒸汽輪機排氣壓力存在最佳值不同。分析原因是混合工質在凝汽器內部所建立的真空除了與環境溫度、凝汽器出口端差有關外,還與混合工質在凝汽器內的蒸汽凝結率有關,使得混合工質在凝汽器內的傳熱過程更加復雜。由2.2節分析可知,蒸汽凝結率和環境溫度限制了凝汽器內部壓力的進一步降低,存在最低壓力。因此,冷卻水量的增加也受到凝汽器內部壓力限制,圖中給出了不同蒸汽凝結率對應的冷卻水流量增加的極限值,可以看出在有限的冷卻水流量增加的過程中,系統凈增功都是單調遞增的,且無法進一步降低凝汽器內部的壓力。綜上分析,氣輪機最佳排氣壓力應與凝汽器內部混合工質最低壓力值保持相等。

圖10 蒸汽凝結率為96%~99%時ΔW隨Δqw變化(twi=25 ℃, Δt2 min=8 ℃)
在相同的環境溫度和凝汽器出口端差下,提高蒸汽凝結率,氣輪機最佳排氣壓力隨之提高,系統凈增功減小,如表3所示。

表3 不同凝結率對應的氣輪機最佳排氣壓力
1)混合工質在凝汽器內為降溫凝結過程,并受到蒸汽凝結率、環境溫度和凝汽器出口端差Δt2 min的限制。凝結率越大,工質溫度降低越劇烈。在換熱過程中特定的蒸汽凝結率下,存在最低壓力pc min,并且凝結率越大,pc min越高。
2)混合工質氣輪機最佳排氣壓力與凝汽器內最低工作壓力相等,并且在本文所研究工況內,混合工質氣輪機最佳排氣壓力均不同程度地高于目前凝汽式電站水蒸汽汽輪機的設計排氣壓力。
3)混合工質氣輪機最佳排氣壓力隨著環境溫度的升高和凝結率的增大而升高。在所研究的工況范圍內,當環境溫度25 ℃,其他參數不變時,蒸汽凝結率提高3%,混合工質氣輪機最佳排氣壓力增加180%,造成系統凈功率減小了56%,系統效率降低了3.3個百分點;當蒸汽凝結率96%,其他參數不變時,環境溫度提高10 ℃,混合工質氣輪機最佳排氣壓力增大79.6%,造成系統凈增功率減小23%,系統效率降低了2.8個百分點。
因此,根據水蒸汽汽輪機排氣壓力選取混合工質氣輪機排氣壓力是不可取的,混合工質氣輪機排氣壓力的確定需考慮蒸汽凝結率、環境溫度和出口端差等綜合影響。