小容量超临界二氧化碳燃煤发电系统的一维设计及回热器紧凑性分析

徐进良1,2,田 贺1,2,孙恩慧1,2,常 诚1,2,郭源东1,2,王兆福1,2

(1.华北电力大学 能源动力与机械工程学院,北京 102206;2.华北电力大学 低品位能源多相流与传热北京市重点实验室,北京 102206)

摘 要:超临界二氧化碳(S-CO2)循环因其高效、紧凑、灵活的优势,有望成为新一代发电动力循环。选用三压缩末级部分压缩的循环结构,将其应用到20 MW S-CO2燃煤发电系统中,并对发电系统的关键部件进行一维设计。设计结果表明,透平和压缩机尺寸整体较小,尺寸优势明显,相比水机组具备明显的竞争力;基于不同容量锅炉的体积变化规律,合理推断由于炉内CO2工质换热系数较小,使得S-CO2锅炉相比水机组锅炉体积有所增加;基于印刷电路板换热器的回热系统占地较大,在本文设计条件下回热器总体积约为18.95 m3,远超现有水机组回热器总体积,问题显著。因此,进一步从回热器换热效果角度探究了回热系统特性,讨论了现有技术面临的困难与瓶颈,发现循环热效率与回热系统体积存在着较强的制约关系,进而对回热系统进行了成本分析,在考虑基础建设费及人工费后,回热系统的投资成本约为9 347.47万元,初始投资成本较高。同时发现,印刷电路板换热器在缩小回热系统体积方面的潜力有限,因此仍需研发新的换热器形式,实现换热系数的量级性提升。最后,针对回热系统体积问题,提出了S-CO2燃煤发电系统的一体化布置方案,通过空间布置的优化实现了系统的整体紧凑性。我国正在推进中等容量S-CO2燃煤发电系统示范,此项研究为示范电站建设提供了新的技术思路。

关键词:超临界二氧化碳;燃煤发电;S-CO2循环;系统一维设计;回热系统;成本分析

我国发布的《“十四五”现代能源体系规划》[1]指出,在未来一段时间内,煤炭仍是我国能源结构主体,仍需大力推动煤电灵活性改造,根据发展需求合理建设先进煤电,充分发挥现有煤电机组应急调峰能力。因此,探索更加高效、低碳的先进煤电系统是我国实现“碳达峰、碳中和”目标的有力支撑。S-CO2燃煤发电系统相对水工质燃煤发电具有效率高、部件紧凑、负荷响应迅速等优势[2],适合应用于燃煤发电领域。

在S-CO2燃煤发电领域,国内外诸多研究机构、学者围绕系统性能优势、关键部件特性等进行了深入研究。法国电力公司MOULLEC[3]率先将S-CO2循环应用于燃煤发电领域,并耦合了碳捕集,热力系统采用二次再热再压缩循环,利用烟气冷却器法吸收烟气余热,系统总利用效率可达41.35%,同时通过与水机组的比较,发现将S-CO2循环引入燃煤发电系统是提升电厂效率的有效途径。MOULLEC提出的系统流程也为后续多数研究者提供了参考。MECHERI等[4]进行了换热器夹点、部件压损、锅炉和循环布置对于燃煤S-CO2循环系统的敏感性分析,在主气参数为30 MPa、620 ℃条件下,通过优化,使得容量为1 009 MW的S-CO2燃煤电厂的总净效率可达47.80%。LI等针对1 000 MW[5]和300 MW[6]燃煤S-CO2热力系统进行了深入研究,并对锅炉[6-7]、热交换器[8-9]以及旋转机械[10-11]等进行概念设计,目前已经建成了5 MW级的S-CO2燃气发电示范装置[12]

现有围绕S-CO2燃煤发电系统的研究多基于再压缩循环开展,相比于水机组能得到一些效率优势,但考虑到计算条件较为理想化,实际效率优势仍未明确。针对该问题,SUN等[13]基于拆分法在再压缩循环基础上提出更为共性的多级压缩循环,并指出三压缩循环可用于燃煤发电领域,在30 MPa、620 ℃主气参数条件下,三压缩循环(Tri-Compression,TC)热效率可达52.54%,相比再压缩循环效率提升0.99%,显著提高了系统效率潜力。进而XU等[14]考虑广义S-CO2循环和耦合锅炉热源2个方面,探索燃煤S-CO2效率极限,综合权衡循环复杂性和效率提升贡献度,将三压缩+二次再热+间冷的循环结构应用于1 000 MW S-CO2燃煤发电系统,在35 MPa、630 ℃主气参数条件下,发电效率51.03%,这是目前大容量S-CO2燃煤发电系统效率的最高水平。

围绕S-CO2燃煤发电,在“十四五”科技部重点研发计划煤炭专项与国家能源局“十四五”能源领域科技创新规划中均布局了20 MW等级的工业示范项目。20 MW容量的示范机组在锅炉、透平、回热系统等方面已初步具备大容量机组各部件的关键特征,同时小容量机组也更能发挥循环紧凑、灵活、高效等特点,具备充当综合能源系统中调节性电源的条件,能够更好地提升系统灵活性,因此有必要针对该容量的燃煤机组进行研究。目前针对小容量系统的研究较少,未能充分发挥系统效率优势,孙恩慧[15]对三压缩循环进行改进,提出末级部分压缩循环,同时系统更为简单,适合应用于燃煤发电领域,然而对于该循环的详细设计还未开展,系统整体的特性仍需挖掘。

基于此,笔者提出基于三压缩末级部分压缩(PTC)+再热(RH)的循环结构,并将其应用到小容量20 MW S-CO2燃煤发电系统中,进而对系统各关键部件进行了一维设计与分析,相对全面地揭示了系统整体特性,并对回热系统存在的紧凑性问题进行讨论,并提出了一种合理的解决方案。

1 末级部分压缩S-CO2燃煤发电系统

1.1 系统流程

基于三压缩末级部分压缩循环(PTC)构建了20 MW燃煤发电系统,图1为PTC+RH循环系统示意,图2为20 MW燃煤发电系统示意。系统采用一次再热布置、烟气冷却器法吸收尾部烟道余热,采用烟气再循环缓解炉膛内冷却壁超温爆管问题,并考虑了透平压缩机泄漏、回热器压降等。

图1 PTC+RH循环系统
Fig.1 Diagram of PTC+RH cycle system

图2 20 MW级S-CO2燃煤发电系统
Fig.2 20 MW coal-fired S-CO2 power plant

图1中,C1、C2、C3为压缩机;T1、T2为透平;LTR、MTR、HTR、HTR2分别为低温回热器、中温回热器、高温回热器、高温回热器;1~14为状态点;x为分流系数。循环部分工作流程为:S-CO2工质在中温回热器MTR低压侧出口分流,一部分工质进入低温回热器LTR低压侧,另一部分进入压缩机C3,经C3压缩后依次进入高温回热器HTR2高压侧和加热器3,进入低温回热器LTR的工质在LTR低压侧出口进一步分流,其中一部分依次进入冷却器、压缩机C1后进入LTR高压侧,另一部分进入压缩机C2,经C2压缩后与LTR高压侧出口工质汇流依次进入MTR、高温回热器HTR、加热器1、透平 T1,在加热器2入口处与进入C3的工质汇流并进入T2做功后再次进行分流,一部分进入HTR 低压侧,另一部分进入HTR2低压侧,2股工质在MTR低压侧入口处汇流后进入MTR,至此完成一个循环。

图2中,UCW为上冷却壁;LUW为下冷却壁;SH1为过热器;RH1为再热器;FGC 为尾部烟气余热吸收器;AP1、AP2为空气预热器;1~17、6′、8′、ab为状态点;Tfg,iTfg,oTfg,ex分别为尾部烟气进入FGC及AP2入口、出口温度。烟气工作流程为:煤粉在炉膛内燃烧,火焰与烟气主要通过辐射方式与下冷却壁LCW和上冷却壁UCW管内的S-CO2工质换热后,烟气依次与布置于锅炉烟道中的过热器SH1、再热器RH1、加热器3以及尾部烟气余热吸收器FGC内的S-CO2工质换热;空气预热器AP2与FGC并列布置在锅炉尾部,同时吸收烟气热量;随后烟气进入空气预热器AP1内与空气换热,最后经过多个热交换过程的低温烟气排出锅炉。图1中加热器1由LCW和SH1组成,加热器2由UCW和RH1组成。

1.2 循环热力学特性

本文计算包括循环热力计算、锅炉热力计算以及部件的一维设计计算,计算程序通过Matlab软件编制而成,CO2的物性参数通过调用REFPROP 9.0获得,表1为系统热力计算的设计条件。

表1 系统设计条件
Table 1 Design condition of system

参数数值透平入口温度T6/℃605透平入口压力P6/MPa27压缩机C1入口温度T1/℃32压缩机C1入口压力P1/MPa8.10空预器入口一次风温度/℃31空预器出口一次风温度/℃320一次风率/%19空预器入口二次风温度/℃23二次风率/%81炉膛出口过量空气系数1.20尾部烟道排烟温度/℃125

表2为循环部件热力计算模型。其中,m为质量流量;T为温度;h为状态点的焓;下角数字为各状态点;下角标s为等熵状态,可由REFPROP物性函数调用;ηC,s为压缩机的等熵效率;ηT,s为透平的等熵效率;ΔTLTR、ΔTMTR、ΔTHTR、ΔTHTR2分别为低温回热器、中温回热器、高温回热器、高温回热器2的夹点温差;xC2xC3分别为分流进入C2、C3的分流系数;xFGC为烟气冷却部分的分流系数;WT为透平所做的功;WC为压缩机耗功。循环的热效率ηth计算公式为

表2 循环部件热力计算模型
Table 2 Thermodynamic models of components

部件循环热力学计算公式C1ηc,s=h2s-h1h2-h1,WC1=(1-xC2-xC3)(h2-h1)mC2ηc,s=h3s-h12h3-h12,WC2=xC2(h3-h12)mC3ηc,s=h13s-h11h13-h11,WC2=xC3(h13-h11)mT1ηr,s=h6-h7h6-h7s,WT1=(h6-h7)mT2ηr,s=h8-h9h8-h9s,WT1=(h8-h9)mLTRT12=T2+ΔTLTR(1-xC2-xC3)(h3-h2)=(1-xC3)(h11-h12)MTRT11=T3+ΔTMTR(1-xC3)(h4-h3)=(h10-h11)HTRT10=T4+ΔTHTR(1-xFGC-xC3)(h5-h4)=(1-x)(h9-h10)HTR2T10=T13+ΔTHTR2xC3(h14-h13)=x(h9-h10)

(1)

式中,WT,total为透平做的总功;WC,total为压缩机的总耗功;Wnet为总净功;Qin为工质吸入的总热量。

S-CO2循环因其高效紧凑的优势,因此有望成为未来更具竞争力的发电循环,但S-CO2循环的部件之间具有很强的耦合效应,循环效率高度依赖部件性能[16]。相关文献研究指出当压缩机效率和透平效率变化2%时,循环效率将分别变化0.5%和1.0%[17],回热器的效能降低5%或压降增加2.5%可能导致循环效率降低1.0%[17-18]。因此,若想对循环性能进行准确分析就需要建立更为详细的部件一维计算模型。此外,更具体的S-CO2燃煤发电系统布局和关键部件设计评估也是证明其商业应用的技术可行性的必要条件[16]

2 末级部分压缩S-CO2燃煤发电系统一维设计

2.1 锅炉与旋转机械一维设计

S-CO2燃煤发电系统中,采用四切角圆燃烧方式π型锅炉,设计煤种性质见表3。

表3 设计煤种特性
Table 3 Properties of the designed coal

元素分析/%工业分析/%CarHarOarNarSarAarMarVdafQnet,ar/(kJ·kg-1)58.563.367.280.790.6319.779.6132.3122 441

锅炉热力计算的模型参考《锅炉机组热力计算方法》第73版,在S-CO2发电系统中,锅炉内的CO2工质接近定物性,因此可采用D-B公式[19]计算其对流换热系数,得到的锅炉参数见表4。

表4 锅炉尺寸计算结果
Table 4 Calculation results of boiler size

锅炉结构参数数值炉膛宽度/m3.84炉膛深度/m3.84炉膛周界/m15.36主炉膛高度/m15.71折焰角区高度/m1.55水平烟道高度/m4.00炉膛有效高度/m20.45燃烧器标高/m8.00折焰角深度/m1.00折焰角上下倾角/(°)30冷灰斗倾角/(°)55下炉膛冷却壁面平均热负荷/(kW·m-2)74.86

通过计算可得,炉膛容积热负荷为155.0 kW/m3,炉膛截面热负荷为3 400 kW/m2,其中锅炉效率为92.99%,燃煤量Bcal为2.17 kg/s,尾部烟道排烟温度设置为125 ℃,根据计算结果可知,在设计条件下锅炉炉膛尺寸为3.84 m×3.84 m×20.45 m。

在S-CO2循环系统中,涡轮机械是系统的重要组成部件,对系统性能的影响显著。对于压缩机以及透平,本文选择三级轴流透平,采用平均线模型[20]设计,压缩机采用单级离心式压缩机模型,采用子午线模型设计,其中关于透平以及压缩机的一维设计计算,具体计算流程如图3、4所示。参考表1给出的设计计算条件,可以得到关于透平以及压缩机的一维设计结构参数,其中表5为单级离心式压缩机的一维设计结构参数,表6为轴流透平的各级结构参数。

图3 透平一维设计计算流程
Fig.3 Flow chart of turbine 1D design calculation

图4 压缩机一维设计计算流程
Fig.4 Flow chart of compressor 1D design calculation

表5 压缩机一维计算结果
Table 5 One-dimensional calculation results of compressors

叶轮一维设计参数C1C2C3压比3.443.392.13转速n/(r·min-1)16 50025 00025 000叶轮进口轮毂直径D1h/mm151520叶轮进口叶尖直径D1s/mm83.5087.3086.20叶轮出口直径D2/mm237.60259.80260.20叶轮出口叶片高度b2/mm8.105.507.60叶顶间隙cl/mm0.400.400.40叶片数Z8+88+88+8叶轮进口叶片安装角/(°)59.9559.5854.01叶轮出口叶片安装角/(°)505050无叶扩压器进口直径D3/mm249.50272.70273.20无叶扩压器出口直径D4/mm285.10311.70312.30

由表5可知,压缩机C1、C2、C3叶轮直径最大处即扩压器出口直径分别约为0.29、0.31、0.31 m,压缩机的叶片直径相对较小,紧凑性较高。根据设计参数可得到压缩机的进口条件以及等熵条件下压缩机的出口焓,参照文献[16,21-23],通过对叶轮以及扩压器的几何设计以及损失模型计算,最终计算得到压缩机的出口的实际焓,进而根据等熵效率定义[18],计算得到循环中压缩机C1、C2、与C3的等熵效率分别为80.0%、78.4%、78.4%,C2与C3效率较低的原因在于2个压缩机的流量偏小,因此需通过变速箱调整其转速以适应运行要求,变速箱效率取98%。

表6为三级轴流透平T1与T2的结构计算结果,对于透平的等熵效率计算,通过透平的入口参数条件可得到透平出口在等熵条件下的焓,参照文献[20,24-26],通过设置各级透平焓降相等的原则,结合损失模型逐级校核计算各级的出口压力以及焓,三级计算结束,最终可得到透平出口的实际焓,由等熵效率计算公式,得到循环中2个透平的等熵效率均为87%[18],其中透平T1与T2的体积分别为0.15、0.68 m3,与现役的水机组汽轮机相比,S-CO2密度大且比容小,相同质量流量情况下S-CO2透平的体积以及叶片参数尺寸紧凑性更高,竞争优势更明显。

表6 透平一维计算结果
Table 6 One-dimensional calculation results of turbines

级数轴向弦长/mm弦长/mm叶片数量轮毂/mm叶高/mm叶尖半径/mmT1第1级静叶7833464894第1级动叶1718373962101第2级静叶7849435396第2级动叶171837425598第3级静叶7849405999第3级动叶171837455095T2第1级静叶787911535150第1级动叶32353611339152第2级静叶789211143154第2级动叶32353610947156第3级静叶789210752158第3级动叶32353610457161

2.2 回热器一维设计

印刷电路板式换热器(Printed Circuit Heat Exchanger,PCHE)是一种微通道热交换器,由大量刻蚀微通道板在高温高压下堆叠组成,工质在刻蚀通道中流通,其基本结构如图5所示,其中,t1为通道距离边缘位置;t2为板厚;t3为通道节距;Dh为刻蚀通道直径;Lw为Z型半波长;α为通道波纹角度;LxLyLz分别为PCHE芯体的宽度、高度、长度。

图5 PCHE几何示意、通道尺寸及形状示意
Fig.5 PCHE geometry diagram,channel size and shape diagram

PCHE中有多种通道类型,如波形、Z字型、直线型、S型等,其中Z型通道具有更加优越的流动传热性能及加工性能[27]。因此,本文系统中所用回热器和冷却器均采用Zigzag型PCHE,通道波纹角度为40o,PCHE流通通道横截面为半圆形[28]。印刷电路板换热器计算预设的几何参数见表7。在PCHE的热力计算中,由于CO2在临界点和近临界点附近的比热变化较大,因此采用分段设计。分段计算条件:设置每个子换热段内的压力恒定;假设与外界无热量交换;计算长度时,并不考虑半圆型通道的曲率,在横向方向上,忽略温度梯度与交叉流的影响的方法[29]。如图6所示(Th,inPh,in分别为热侧入口温度、压力;Tc,outPc,out分别为冷侧出口温度、压力;分别为第i个子单元上热侧平均温度、压力;Ti-1hPi-1h分别为第i-1子单元上热侧温度、压力;Ti-1cPi-1c分别为第i-1子单元上冷侧温度、压力;分别为第i个子单元上热侧温度、压力;分别为第i个子单元上冷侧温度、压力;分别为第i个子单元上冷侧平均温度、压力;Th,outPh,out分别为热侧出口温度、压力;Tc,inPc,in 分别为冷侧入口温度、压力;N为总子单元数;i为第i个子单元;下角h、c分别为热侧和冷侧)。将单个换热器分成N段,本文选取子单元数为20。假设已知热流体入口温度、冷流体入口温度、指定冷流体的出口温度、子单元换热系数Ki(i=1,2,…,n),子单元上热流体释放热量全部被该段冷流体吸收[29]。应用分段设计方法,对换热器进行分段,取分段数为N,子换热单元的换热量为q

图6 PCHE分段设计计算示意
Fig.6 PCHE discrete design calculation diagram

表7 PCHE初始参数设置
Table 7 Initial parameter settings for PCHE

参数数值刻蚀通道直径(dc/dh)/mm[29]2.0Z型半波长Lw/mm4.5芯体宽度Lx/m[29]0.6高度上限值Ly,ml/m[29]8通道距边缘距离t1/m[30]0.01板厚t2/mm0.75通道节距t3/mm0.6选用钢材316L

Q=mhcp,h(Th,in-Th,out)

(2)

(3)

式中,Q为热交换器总换热量;mh为热侧质量流量;cp,h为热侧比容;in为入口;out为出口。

由于高温侧工质密度小,且高温侧压降对循环效率影响更大,因此取高温侧与低温侧通道数量之比为2[29-30],记为Ny,h=2Ny,c,由此可以得到PCHE内高温侧与低温侧的换热面积为

A=(1+0.5π)dNxNyLc

(4)

式中,A为传热面积;d为通道直径,热侧、冷侧设置的通道直径相等,因此均为d;NxNy为PCHE内xy方向通道数量;Lc为通道的实际长度。

子换热单元的总传热系数为K,热导为UA,采用对数平均温度方法(LMTD)计算:

(5)

(6)

(7)

ΔTmax=max(Th,in-Tc,out,Th,out-Tc,in)

(8)

ΔTmin=min(Th,in-Tc,out,Th,out-Tc,in)

(9)

式中,λ为材料导热系数,W/(m·K);te为等效管壁厚;对于子换热单元,hhhc分别为高温侧与低温侧管内对流传热系数:

(10)

式中,Nu为努塞尔数;λhλc分别为冷热侧工质导热系数,W/(m·K)。

对于回热器来说,两侧工质均为二氧化碳,由文献[29]可知,对于回热器,其中努塞尔数Nu及范宁摩擦阻力系数f关联式为

(11)

Nu=0.084 5Re0.721Pr0.333

(12)

f=1.336Re-0.126 8

(13)

其中,G为工质质量流速,kg/(m2·s);d为通道直径;μ为黏度。对于冷却器,热侧为CO2,冷侧为水,其中冷侧努塞尔数Nu及范宁摩擦阻力系数f关联式[31]

Nu=0.122Re0.56Pr0.14

(14)

f=(1.12ln Re-0.126 8+0.85)-2

(15)

微通道内流动摩擦阻力(即为压降)为

(16)

(17)

其中,ρ为密度;d′为水力直径。流道截面为半圆形。

结合循环的热力计算条件,可得到不同热交换器的一维设计计算参数见表8,计算可得LTR、MTR、HTR、HTR2、冷却器体积分别为4.21、4.40、8.29、2.05、3.11 m3。图7为基于上述计算方法得到的LTR内部子单元模块沿程换热系数与温度分布。

表8 换热器一维设计计算结果
Table 8 One-dimensional calculation results of exchangers

参数LTRMTRHTRHTR2冷却器Rp22221Lz/m0.820.651.011.071.03Lx/m0.60.60.60.60.6Ly/m8.6011.3813.803.185.07冷侧压降/kPa10010010010060热侧压降/kPa150150150150100换热面积/m22 139.81 420.12 389.1558.21 548.8热负荷/MW32.3529.2842.709.4724.80冷侧板数1 6601 8812 189494966热侧板数3 3193 7624 3789881 932

图7 低温回热器LTR内部子单元模块沿程换热系数和温度的分布
Fig.7 Distribution of heat transfer coefficient and temperature along the course of the subunit modules inside the LTR

2.3 一维设计结果分析

对系统部件进行一维设计计算后,可得到关于涡轮机械的等熵效率(表9)。由于发电效率为循环热效率、锅炉效率和管道效率的乘积。因此,选取管道效率为99%[32-33],考虑管结构参数、材料耐高温极限以及系统循环换热效率等约束条件,参考现役超超临界水机组主汽参数[34],结合课题组已有的研究成果[2,13,15],本文主气参数选取为27 MPa、605 ℃,综合考虑循环热力计算以及部件一维设计计算,可得到系统循环热效率为43.49%,发电效率为40.20%。

表9 部件性能一维设计结果
Table 9 One-dimensional calculation results of component performance

部件透平T1透平T2压缩机C1压缩机C2压缩机C3等熵效率/%87.087.080.078.478.4

根据部件的一维计算结果,轴流透平T1与T2的体积分别为0.15、0.68 m3,压缩机C1、C2、C3叶轮直径最大处即扩压器出口直径分别约为0.29、0.31、0.31 m。可以看出,S-CO2发电系统透平以及压缩机的尺寸均较小,紧凑性好,与水机组相比,竞争优势明显。

对于S-CO2锅炉,XU等[35]比较了容量同为1 000 MW、蒸汽参数相同的2个发电系统,发现由于S-CO2锅炉受CO2传热能力的影响,受热面用钢等级和厚度均增加,锅炉结构复杂,使得其成本较水蒸气锅炉成本增加了36.3%。因此体积也有所增大。

在西安热工院300 MW燃煤S-CO2锅炉概念设计中[6],经计算得到300 MW级S-CO2锅炉炉膛尺寸约为12.5 m×12.5 m×52.2 m,炉膛体积为8 156.25 m3。依据文献[36],某主汽参数为16.76 MPa、555 ℃的亚临界锅炉,炉膛体积约为7 900.40 m3。北京巴威公司制造的亚临界压力300 MW自然循环锅炉[36],锅炉炉膛尺寸约为13.35 m×12.3 m×46.40 m,炉膛体积仅为7 619.11 m3

陈玮等[37]对50 MWS-CO2燃煤发电系统进行设计与参数分析,S-CO2循环系统锅炉采用倒π形布置,在主气参数为25 MPa、620 ℃的条件下,计算得到锅炉炉膛体积为1 641.99 m3;而某传统50 MW水机组[38],在主蒸汽参数为8.62 MPa、535 ℃条件下,锅炉炉膛体积为688.82 m3

因此,在相同容量下,对比燃煤水机组锅炉,S-CO2锅炉炉膛体积均偏大。由于公开发表的20 MW燃煤发电水机组锅炉的文献较少,只能通过其他容量锅炉进行比较,经合理推测后,相比于20 MW燃煤发电水机组,20 MWS-CO2锅炉体积稍有增大。

目前最突出的问题集中在回热系统上,本文所采用的末级部分压缩、三压缩循环共包含4个回热器,其回热系统的总体积约为18.95 m3,回热器系统的占地面积过大,远超现有水机组回热系统,问题明显。

3 基于印刷电路板换热器面临的紧凑性困境

S-CO2循环具有典型气态循环的特点,可通过较高的回热程度显著提高循环效率,然而较高的回热程度也意味着较大的回热面积,进而导致成本增大。因此,开发高效、紧凑、经济的回热器对于S-CO2循环的应用至关重要。

PIERRES[39]、MUSGROVE[40]等研究认为,印刷电路板换热器(PCHE)具有较高的基材强度和耐高压能力,同时具备紧凑性高、换热面积大及换热性能好等特点,适合作为S-CO2循环的换热设备。目前PCHE在S-CO2循环中应用广泛,如韩国先进科学技术研究院[41]、Sandia国家实验室[42]、东京工业大学[43]以及西安热工研究院[9]的实验中均采用此类换热器。

3.1 回热系统体积分析

由2.3节可知,在一维计算中基于PCHE的回热器体积并没有体现出明显的紧凑性优势。在循环层面,回热器对循环效率的影响主要通过压降与夹点温差,这2个参数增大均会使得循环效率显著下降。

对于回热系统的总体积,由图8可知,夹点温差不变,随着压降增大,回热系统的总体积减小;压降不变时,随着夹点温差增大,回热系统总体积减小。夹点温差越小,压降对换热体积的影响越显著,夹点温差为2 ℃时,压降每增加50 kPa,回热系统总体积会显著减小,最大可减少27.02 m3;压降不变时,随夹点温差的逐渐增大,回热系统总体积逐渐减小,在8 ℃以上时变化趋势放缓。

图8 回热系统总体积随夹点温差和压降的变化
Fig.8 Total volume of the regeneration system changed with the pinch temperature difference and pressure drop

经反复计算,本文设计的回热系统夹点温差为10 ℃、压降为150 kPa,回热系统总体积为18.95 m3。由于回热系统性能与循环效率成正比,但与经济性成反比,因此仍需对回热系统进行成本分析。

3.2 回热系统成本分析

WEILAND等[44]认为,回热系统的建设成本约占整个机组建设成本的24%。XU等[35]研究发现,在三压缩S-CO2燃煤发电系统中,由于回热程度的提升,回热系统的投资成本占整个设备投资成本的44%左右。因此回热器的投资费用占总投资费用的比重较大。目前针对换热器成本的计算模型主要分为两大类:一类是基于换热器体积的正比例函数计算模型;另一类是基于换热器热导(UA)的幂函数计算模型,UA=Q/ΔTlm(其中ΔTlm为对数平均温差),计算时,通常将换热器离散为多个单元[45]

ZADA等[46]研究发现,换热器成本模型需要将热导与换热器的效能以及平均传热温差相结合,提出近似于正比例函数的换热器成本计算模型,并给出了单位热导的成本计算参考值,对于高温回热器,每单位热导成本为1.49~2.21$;对于低温回热器,每单位热导成本为1.19~1.77$。

WEILAND等[44]基于回热器热导的幂函数计算模型,提出了同时适用于高温回热器以及低温回热器的计算模型为

Crecup=49.45UA0.754 4fT,recup

(18)

fT,recup=1(Tmax<550 ℃)

(19)

fT,recup=1+0.021 41(Tmax-550),(Tmax≥550 ℃)

(20)

其中,Crecup为回热器成本;fT,recup为温度修正因子;Tmax为回热器沿程或两侧最高温度。该模型的限制条件为换热器热导UA计算范围为1.60×105~2.15×108 W/K,高温回热器的可承受最高温度范围为407~584 ℃,低温回热器为165~252 ℃;回热器允许压降范围为0.07~0.40 MPa,低温侧的最大使用压力为21~32 MPa[44-45]

DOSTAL等[47]基于换热器的质量以及体积等提出了基于换热器体积的正比例函数计算模型,其成本公式为

(21)

(22)

式中,LWH分别为PCHE的长、宽、高;ρ′为换热器材料密度(通常设定为不锈钢7 800 kg/m3);Cm为PCHE单位质量成本(不锈钢成本价为30$/kg);p为通道间距;t为板厚;fm为换热器耐热材料体积分数,当d=2 mm、p=2.60 mm、t=1.75 mm时,PCHE的fm=0.66。

本文选定上述3种换热器成本计算模型,其中Zada模型为case1,Weiland模型为case2,Dostal模型为case3,如图9所示,以LTR的case2成本模型为例,发现压降对于回热系统的投资成本影响并不显著,因此在分析回热系统总投资成本时,认为夹点温差是主导的影响因素。

图9 LTR成本随压降以及夹点温差的变化
Fig.9 LTR cost changed with the pressure drop and the pinch temperature difference

基于3个成本模型在不同夹点温差下计算得到的回热系统总成本,如图10、表10所示。由图10可知,不同的模型成本计算结果差别较大,相比于case3以回热器体积为基础,case1与case2所采用的基于热导的计算模型会更合理。而相比于case1,case2通过线性变化对换热器运行温度的处理更为符合实际情况,因此采用case2的成本模型来计算回热系统的总制造成本。

图10 3个成本模型下回热系统成本随夹点温差的变化
Fig.10 Cost of regeneration system changed with the pinch temperature difference under three cost models

表10 3个成本模型下夹点温差对应的回热系统成本
Table 10 Cost of the regeneration system corresponding to the pinch temperature difference under three cost models

夹点温差/℃回热系统成本/106$case1case2case3255.2832.5314.51624.4917.414.931016.2512.772.881511.359.771.81208.637.981.28

根据3.1节对回热系统的体积分析可知,当夹点温差10 ℃,压降取150 kPa,此时回热系统的总体积维持在18.95 m3左右。按照case2的成本模型计算,由表10回热系统的总制造成本为12.77×106$。此外,公式计算得到的是PCHE的制造成本,在后期的建设过程中还要考虑基础建设成本(约为总成本的2%)以及人工费(约为总成本的3%)[44],因此回热系统的制造成本约为总制造成本的1.05倍,即为13.41×106$,折合人民币约为9 347.47万。

鉴于目前循环热效率与回热系统经济性的制约关系,在没有更好回热器选择或技术突破的情况下,回热系统的体积难以进一步降低,其投资成本也会处于相对较高的水平。

3.3 回热系统的紧凑性潜力

以低温回热器为例,根据Heatric公司[48]发布的数据,PCHE的刻蚀通道直径最小为0.5 mm,最大为5 mm,基于此计算得到在不同刻蚀通道直径下对流换热系数与回热器压降的关系曲线如图11所示,随着对流换热系数的不断增大,回热器的体积逐渐缩小,呈负相关关系。根据2.3节可知,在现有技术条件下,回热系统的体积优势并不明显,因此需要针对待开发区域进行展望。在待开发领域,根据PCHE换热器的对流换热系数与回热器体积的曲线变化趋势可以看出,如果继续采用PCHE型换热器,即便在未来能够实现更小的通道直径,对流换热系数也不会实现量级的变化,回热器的体积不会得到明显的减小,还是维持在相对较高的水平。因此在表面式换热的形式下,S-CO2燃煤发电的回热系统难以体现出紧凑性的优势,仍需基于新原理开发新型换热技术。

图11 对流换热系数与回热器体积的关系及展望
Fig.11 Relationship between convective heat transfer coefficient and regenerative heat exchanger volume and outlook

3.4 机炉一体化布置方案

受目前技术水平限制,在S-CO2燃煤发电领域,具备应用前景的仍是PCHE,但由上述的分析可知,S-CO2机组回热量大,应用PCHE仍难以体现紧凑性优势,如何解决回热系统的占地是节约投资成本的关键问题。

为此,笔者提出了一种S-CO2燃煤发电系统的一体化布置方案,如图12所示,各回热器均采用串并联模块组合方式围绕布置在炉膛周围,同时透平采用高位布置。回热器采用串并联模块布置,可解决回热系统占地面积过大的问题,更大程度地减小系统压降;透平高位布置可极大缩短输送工质的管道长度,节约投资成本。

图12 S-CO2循环燃煤发电系统一体化布置
Fig.12 Integrated arrangement of S-CO2 cycle coal-fired power generation system

4 结 论

(1)换热器体积与循环热效率存在制约关系,在现有技术水平下,回热系统难以在保证循环高效的同时体现紧凑性与经济性优势。当夹点温差为10 ℃、压降为150 kPa时,回热系统总体积为18.95 m3,系统循环热效率为44.04%,考虑基础建设费及人工费,回热系统的投资成本约为9 347.47万元。

(2)若继续采用印刷电路板换热器,即使能开发出更小的刻蚀通道,回热器体积也不会显著减小,目前仍需研发新的换热器,实现换热系数量级性提升。

(3)针对回热系统体积过大问题,提出了一种高效、紧凑的S-CO2燃煤发电系统一体化布置方案,即采用回热器串并联模块围绕布置在炉膛周围的形式以及透平高位布置。该方案能够很好地解决回热系统占地面积过大的问题,为S-CO2燃煤发电系统的设计提供新的思路。

参考文献(References):

[1] 国家能源局.关于印发《“十四五”现代能源体系规划》的通知[EB/OL](2022-01-29)[2022-12-01].http://zfxxgk.nea.gov.cn/2022-01/29/c_1310524241.htm.

[2] XU Jinliang,SUN Enhui,LI Mingjia,et al.Key issues and solution strategies for supercritical carbon dioxide coal fired power plant[J].Energy,2018,157:227-246.

[3] LE Moullec Y.Conceptual study of a high efficiency coal-fired power plant with CO2 capture using a supercritical CO2 Brayton cycle[J].Energy,2013,49:32-46.

[4] MECHERI M,LE MOULLEC Y.Supercritical CO2 Brayton cycles for coal-fired power plants [J].Energy,2016,103:758-771.

[5] LI Hongzhi,ZHANG Yifan,YANG Yu.Preliminary design assess-

ment of supercritical CO2 cycle for commercial scale coal-fired power plants[J].Applied Thermal Engineering,2019,158:113785.

[6] BAI Wengang,ZHANG Yifan,YANG Yu,et al.300 MW boiler design study for coal-fired supercritical CO2 Brayton cycle[J].Applied Thermal Engineering,2018,135:66-73.

[7] 张一帆,王安,白文刚,等.300 MW超临界二氧化碳锅炉气动力特性及壁温分布[J].中国电机工程学报,2019,39(6):1700-1706,1866.

ZHANG Yifan,WANG An,BAI Wengang,et al.Thermal fluid dynamic characteristics and wall temperature distribution of a 300 MW supercritical CO2 boiler[J].Proceedings of the CSEE,2019,39(6):1700-1706,1866.

[8] 吴家荣,李红智,杨玉,等.超临界二氧化碳动力循环中印刷电路板换热器芯体机械应力和热应力耦合分析[J].中国电机工程学报,2022,42(2):640-650.

WU Jiarong,LI Hongzhi,YANG Yu,et al.Coupling analysis of mechanical stress and thermal stress of printed circuit heat exchanger core in supercritical carbon dioxide power cycle[J].Proceedings of the CSEE,2022,42(2):640-650.

[9] 张一帆,李红智,姚明宇,等.600 MW煤基超临界二氧化碳发电系统回热器和预冷器的概念设计[J].中国电机工程学报,2017,37(24):7223-7229,7436.

ZHANG Yifan,LI Hongzhi,YAO Mingyu,et al.Conceptual design of the recuperator and precooler for a 600 MW fossil-based supercritical CO2 power generation system[J].Proceedings of the CSEE,2017,37(24):7223-7229,7436.

[10] 韩万龙,丰镇平,王月明,等.超临界二氧化碳高压涡轮气动设计及性能[J].哈尔滨工业大学学报,2018,50(7):192-198.

HAN Wanlong,FENG Zhenping,WANG Mingyue,et al.Aerodynamic design and performance of S-CO2 high pressure turbines[J].Journal of Harbin Institute of Technology,2018,50(7):192-198.

[11] 韩万龙,王月明,丰镇平,等.超临界二氧化碳涡轮叶栅端壁附面层分层流动现象研究[J].热力发电,2019,48(2):16-22.

HAN Wanlong,WANG Yueming,FENG Zhenping,et al.Stratified flow phenomena of the end wall boundary layer in supercritical CO2 turbine cascades[J].Thermal Power Generation,2019,48(2):16-22.

[12] LI Hongzhi,ZHANG Yifan,YAO Mingyu,et al.Design assessment of a 5 MW fossil-fired supercritical CO2 power cycle pilot loop[J].Energy,2019,174:792-804.

[13] SUN Enhui,XU Jinliang,LI Mingjia,et al.Synergetics:The cooperative phenomenon in multi-compressions S-CO2 power cycles[J].Energy Conversion and Management,2020,X7:100042.

[14] XU Cheng,ZHANG Qiang,YANG Zhiping,et al.An improved supercritical coal-fired power generation system incorporating a supplementary supercritical CO2 cycle[J].Applied Energy,2018,231:1319-1329.

[15] 孙恩慧.超高参数二氧化碳燃煤发电系统热力学研究[D].北京:华北电力大学(北京),2020.

SUN Enhui.Thermodynamic study of ultra-high parameter CO2-fired power generation system[D].Beijing:North China Electric Power Universit(Beijing),2020.

[16] YAO Lichao,ZOU Zhengping.A one-dimensional design methodology for supercritical carbon dioxide Brayton cycles:Integration of cycle conceptual design and components preliminary design[J].Applied Energy,2020,276:115354.

[17] CHENG Wenlong,HUANG Wenxu,NIAN Yongle.Global param-eter optimization and criterion formula of supercritical carbon dioxide Brayton cycle with recompression[J].Energy Conversion and Management,2017,150:669-677.

[18] BRUN K,FRIEDMAN P,DENNIS R,et al.Fundamentals and applications of supercritical carbon dioxide(sCO2)based power cycles[M].Cambridge:Woodhead Publishing an imprint of Elsevier,2017.

[19] 杨世铭,陶文铨.传热学[M].北京:高等教育出版社,2006.

[20] SALAH S I,KHADER M A,WHITE M T,et al.Mean-line design of a supercritical CO2 micro axial turbine[J].Applied Sciences,2020,10(15):5069.

[21] 王彦方.跨临界CO2热泵离心式压缩机研究[D].北京:中国科学院大学,2021.

WANG Yanfang.Research on centrifugal compressor for transcritical CO2 heat pump[D].Beijing:University of Chinese Academy of Sciences,2021.

[22] 黄怡婷.超临界二氧化碳压缩机开发及压缩过程性能分析[D].上海:上海交通大学,2020.

HUANG Yiting.Development on supercritical CO2 compressor and performance analysis on compression process[D].Shanghai:Shanghai Jiao Tong University,2020.

[23] ELKIN I,GUTIERREZ VELASQUEZ.Determination of a suitable set of loss models for centrifugal compressor performance prediction[J].Chinese Journal of Aeronautics,2017,30(5):1644-1650.

[24] DA LIO L,MANENTE G,LAZZARETTO A.Predicting the optimum design of single stage axial expanders in ORC systems:Is there a single efficiency map for different working fluids?[J].Applied Energy,2016,167:44-58.

[25] AUNGIER R H.Axial-flow and radial inflow turbine design and an-alysis[J].Turbine Aerodynamics,2006,215:751-769.

[26] DA LIO L,MANENTE G,LAZZARETTO A.New efficiency charts for the optimum design of axial flow turbines for organic Rankine cycles[J].Energy,2014,77:447-459.

[27] LIU Gangxu,HUANG Yanping,WANG Junfeng,et al.A review on

the thermal-hydraulic performance and optimization of printed circuit heat exchangers for supercritical CO2 in advanced nuclear power systems[J].Renewable and Sustainable Energy Reviews,2020,133:110290.

[28] SOUTHALL D,PIERRES R L,DEWSON S J.Design consider-ati-ons for compact heat exchangers[C]//Proceedings of ICAPP’08,Anaheim.CA USA,2008:8-12.

[29] JIANG Y,LIESE E,ZITNEY S E,et al.Design and dynamic modeling of printed circuit heat exchangers for supercritical carbon dioxide Brayton power cycles[J].Applied Energy,2018,231:1019-1032.

[30] NESTELL J,SHAM T L.ASME code considerations for the compact heat exchanger[R].Oak Ridge National Laboratory,2015.

[31] KIM S G,LEE Y,AHN Y,et al.CFD aided approach to design printed circuit heat exchangers for supercritical CO2 Brayton cycle application[J].Annals of Nuclear Energy,2016,92:175-185.

[32] 石奇光,仝宁,高乃平.火电厂管道热效率及行业标准的探讨[J].华东电力,2004,32(1):44-47.

SHI Qiguang,TONG Ning,GAO Naiping.Exploration of thermal efficiency of thermal power plant pipelines and industry standards[J].East China Electric Power,2004,32(1):44-47.

[33] 中华人民共和国国家发展和改革委员会.火力电厂热平衡导则第三部分:热平衡:DL/T606.3-2006[S].北京:中国电力出版社,2007.

[34] 高子瑜,徐雪元,丘加友,等.上海锅炉厂有限公司百万等级超超临界П型锅炉方案简介[J].锅炉技术,2006(S1):1-7,11.

GAO Ziyu,XU Xueyuan,QIU Jiayou,et al.Introduction to the million-class ultra-supercritical П-type boiler scheme of Shanghai boiler plant company Limited[J].Boiler Technology,2006(S1):1-7,11.

[35] XU Jinliang,WANG Xue,SUN Enhui,et al.Economic comparison between sCO2 power cycle and water-steam Rankine cycle for coal-fired power generation system[J].Energy Conversion and Management,2021,238:114150.

[36] 林荫平.电站锅炉手册[M].北京:中国电力出版社,2013.

[37] 陈玮,罗向龙,梁颖宗,等.50 MW超临界二氧化碳燃煤发电系统设计与参数分析[J].工程热物理学报,2022,43(7):1824-1835.

CHEN Wei,LUO Xianglong,LIANG Yingzong,et al.System design and parametric analysis of 50 MW supercritical carbon dioxide coal-fired power plant[J].Journal of Engineering Thermophysics,2022,43(7):1824-1835.

[38] 王志峰.小容量高参数燃煤机组热电联供热力系统设计与研究[D].北京:华北电力大学,2019.

WANG Zhifeng.Cogeneration system design and research of small capacity coal-fired units with advanced steam condition[D].Beijing :North China Electric Power University,2019.

[39] LEPIERRES R,SOUTHALL D,OSBRONE S.Impact of mechanical design issues on printed circuit heat exchangers[C]//Proceedings of sCO2 power cycle symposium 2011.University of Colorado at Boulder,2011:24-25.

[40] MUSGROVE G,SULLIVAN S,SHIFERAW D,et,al.Heat exchan-

gers[M]//BRUN K,FRIEDMAN P,DENNIS R.Fundamentals and applications of supercritical carbon dioxide.(sCO2)based power cycles.Elsevier,2017.

[41] KIM I H,NO H C.Thermal hydraulic performance analysis of a pr-inted circuit heat exchanger using a helium-water test loop and numerical simulations[J].Applied Thermal Engineering,2011,31:4064-4073.

[42] KRUIZENGA A,FLEMING D,CARLSON M,et al.Supercritical

CO2 heat exchanger fouling[C]//Proceeding of the 4th international symposium-supercritical CO2 power cycles.Pittsburgh,PA,2014:9-10.

[43] NGO T L,KATO Y,NIKITIN K,et al.Heat transfer and pressure

drop correlations of microchannel heat exchangers with S-shaped and zigzag fins for carbon dioxide cycles[J].Experimental Thermal and Fluid Science,2008,32(2):560-570.

[44] WEILAN N T,LANCE B W,PIDAPARTI S R.sCO2 power cycle component cost correlations from DOE data spanning multiple scales and applications[C]//Proceedings of the ASME Turbo Expo:Turbomachinery Technical Conference and Exposition.Phoenix,Arizona,USA,2019.

[45] 王雪,孙恩慧,徐进良,等.超临界二氧化碳循环关键部件成本模型研究进展[J].中国电机工程学报,2022,42(2):650-663.

WANG Xue,SUN Enhui,XU Jinliang,et al.Advances in cost modeling of key components of supercritical CO2 cycle[J].Proceedings of the CSEE,2022,42(2):650-663.

[46] ZADA K R,KIM R,WILDBERGER A.Analysis of supercritical CO2 brayton cycle recuperative heat exchanger size and capital cost with variation of layout design[C]//6th International supercritical CO2 power cycles symposium.Pittsburgh,2018.

[47] DOSTAL V,HEJZLAR P,DRISCOLL M J.The supercritical carbon dioxide power cycle:Comparison to other advanced power cycles[J].Nuclear Technology,2006,154(3):283-301.

[48] MEGGITT.Heatric[EB/OL].(2017-08-15)[2022-12-01].www.heatric.com.

One dimensional design of small capacity supercritical CO2 coal fired power generation system and analysis of regenerator compactness

XU Jinliang1,2,TIAN He1,2,SUN Enhui1,2,CHANG Cheng1,2,GUO Yuandong1,2,WANG Zhaofu1,2

(1.School of Energy,Power and Mechanical EngineeringNorth China Electric Power UniversityBeijing 102206,China;2.Beijing Key Laboratory of Multiphase Flow and Heat Transfer for Low Grade Energy UtilizationNorth China Electric Power UniversityBeijing 102206,China)

Abstract:Supercritical CO2(S-CO2)cycle is expected to be the next generation power cycle for electricity generation because of its high efficiency,compactness and flexibility.In this study,a tri-compression with last-stage partial compression S-CO2 cycle was applied to a 20 MW S-CO2 coal-fired power generation system,and the one-dimensional designs of components were conducted.It is shown that the size of turbine is significant smaller than that of the traditional steam turbine,thus it has obvious competitiveness compared with the water steam unit.Based on the volume change law of boilers with different capacities,it is reasonably inferred that the volume of S-CO2 boiler is larger than the steam boiler due to the smaller heat transfer coefficient of the CO2 in the boiler.The regeneration system based on the printed circuit heat exchanger covers a large area with a volume around 18.95 m3,which is far more than the existing regeneration system of the water steam unit and the problem is significant.Therefore,this study further explored the characteristics of the regeneration system from the perspective of the heat transfer of the regeneration,discussed the difficulties faced by the existing technology,and found that there is a strong constraint relationship between the cycle thermal efficiency and the volume of the regeneration system,and then conducted an economic analysis of the regeneration system.After considering the infrastructure and labor costs,the investment cost of the regeneration system is about 93 474 700 yuan,which is a high initial investment cost.It was also found that the printed circuit heat exchanger has limited potential to reduce the size of the regeneration system,thus new heat exchanger types still need to be developed to increase the magnitude in heat transfer coefficient.Finally,to address the size problem of the regeneration system,this study proposed an integrated arrangement scheme for the S-CO2 coal-fired power generation system,which achieves the overall compactness of the system by optimizing the space arrangement.Currently,the demonstration of S-CO2 coal-fired power plants is promoted in China,and this research provides some novel technical ideas for the demonstration system.

Key words:S-CO2 ;coal-fired power plant;S-CO2 cycle;one-dimensional design;regeneration system;cost analysis

中图分类号:TK16

文献标志码:A

文章编号:0253-9993(2023)01-0438-14

移动阅读

收稿日期:2022-12-22

修回日期:2023-01-13

责任编辑:张晓宁

DOI:10.13225/j.cnki.jccs.EX22.1861

基金项目:国家自然科学基金资助项目(52130608,52206010)

作者简介:徐进良(1966—),男,安徽东至人,教授,博士。Tel:010-61771091,E-mail:xjl@ncepu.edu.cn

通讯作者:孙恩慧(1991—),男,黑龙江七台河人,讲师,博士。E-mail:ehsun@ncepu.edu.cn

引用格式:徐进良,田贺,孙恩慧,等.小容量超临界二氧化碳燃煤发电系统的一维设计及回热器紧凑性分析[J].煤炭学报,2023,48(1):438-451.

XU Jinliang,TIAN He,SUN Enhui,et al.One dimensional design of small capacity supercritical CO2 coal fired power generation system and analysis of regenerator compactness[J].Journal of China Coal Society,2023,48(1):438-451.