压缩空气储能系统动态运行特性.pdf

第 12 卷 第 6 期2023 年 6 月Vol.12 No.6Jun.2023储能科学与技术Energy Storage Science and Technology压缩空气储能系统动态运行特性孙晓霞1，桂中华1，高梓玉2,3，周冰倩2,3，刘夏2，张新敬2,3,4，郭欢2,3，李文2,3，盛勇2,3，朱阳历2,3，周健1，徐玉杰2,3（1国网新源控股有限公司抽水蓄能技术经济研究院，北京 100761；2中国科学院工程热物理研究所，北京 100190；3中国科学院大学，北京 100049；4中科南京未来能源系统研究院，江苏 南京 211135）摘要：发展基于可再生能源为主体的新型电力系统，支撑“碳达峰、碳中和”战略目标的实现，由于风、光等可再生能源的间歇性、波动性、周期性等特点，需要集成大规模长时储能系统，提升风光等可再生能源发电的品质与可控性，压缩空气储能具有效率高、成本低、环境友好等优点，被认为是最具发展潜力的大规模长时储能技术。压缩空气储能系统通常为定容储气，因此其储能(储气)过程与释能(释气)过程处于动态，本工作围绕储/释能过程的压力变化，开展了压缩空气储能系统不同运行模式特性研究，建立了部件的动态模型，通过仿真获得了系统主要部件的工作特性，以及系统的总体性能。研究结果表明，在释能过程采取定压和滑压结合模式和扩大储气室压力变化范围可以提高TS-CAES系统效率和能量密度。释能时间为6 h，系统效率和能量密度分别为 73.98%、26.49 MJ/m3。关键词：压缩空气储能系统；动态特性；定压运行；滑压运行；系统性能doi：10.19799/ki.2095-4239.2023.0181 中图分类号：TK 02 文献标志码：A 文章编号：2095-4239（2023）06-1840-14Dynamic characteristics of compressed air energy storage systemSUN Xiaoxia1,GUI Zhonghua1,GAO Ziyu2,3,ZHOU Bingqian2,3,LIU Xia2,ZHANG Xinjing2,3,4,GUO Huan2,3,LI Wen2,3,SHENG Yong2,3,ZHU Yangli2,3,ZHOU Jian1,XU Yujie2,3(1Pumped-storage Technological&Economic Research Institute,Beijing 100761,China;2Institute of Engineering Thermophysics,Chinese Academy of Sciences,Beijing 100190,China;3University of Chinese Academy of Sciences,Beijing 100049,China;4Nanjing Institute of Future Energy System,Nanjing 211135,Jiangsu,China)Abstract:Renewable energy is aimed to be the main part of a new electrical system to support the strategic goal of“Carbon Peak,Carbon Neutrality”;however,due to the drawbacks of intermittence,fluctuation,and periodicity of renewable energy,large-scale,long-duration energy storage systems urgently need to improve the quality and flexibility of renewable energy.Compressed air energy storage(CAES)is considered one of the most promising large-scale long-duration energy storage technologies with high efficiency,low cost,and environment-friendly merits.Generally,the CAES system utilizes constant-volume storage 储能系统与工程收稿日期：2023-03-27；修改稿日期：2023-04-03。基金项目：北京市自然科学基金（JQ21010），国家杰出青年科学基金（51925604），国家电网公司科技项目（525730210003）。第一作者：孙晓霞（1990），女，工学硕士，工程师，研究方向为抽水蓄能及新型储能技术应用，E-mail：；通讯作者：张新敬，博士，研究员，研究方向为压缩空气储能，E-mail：。引用本文：孙晓霞,桂中华,高梓玉,等.压缩空气储能系统动态运行特性J.储能科学与技术,2023,12(6):1840-1853.Citation：SUN Xiaoxia,GUI Zhonghua,GAO Ziyu,et al.Dynamic characteristics of compressed air energy storage systemJ.Energy Storage Science and Technology,2023,12(6):1840-1853.第 6 期孙晓霞等：压缩空气储能系统动态运行特性caverns.Thus,the charging and discharging processes are under dynamic conditions,especially the storage pressure.Various CAES operation modes,including dynamic component features,are investigated due to the dynamic pressure conditions and system modeling.Similarly,the operation characteristics and performance of both component-level and system-level are analyzed.The results show that the combination of constant pressure and sliding pressure mode in the discharging process,and enlarging the pressure range of the air chamber,can improve the round-trip efficiency and energy density of the TS-CAES system,which are 73.98%and 26.49 MJ/m3,respectively,at the discharging time of 6 hours.Keywords:compressed air energy storage;dynamic characteristics;constant-pressure operation;sliding-pressure operation;system performance为了推动能源低碳绿色发展，应对全球气候变化，我国提出了“碳达峰、碳中和”的发展战略1。实施“双碳”战略是保障能源安全和实现能源强国的重要举措，2060年实现碳中和，要求我国非化石能源比例由目前的20%上升到80%2，未来能源系统将是以新能源为主体、多种形式能源共同构成的多元化能源系统。然而风力发电、光伏发电本身的波动性和间歇性决定了灵活性将是新能源系统必不可少的组成部分3。现有电力系统无法适应大规模可再生能源的接入，难以满足新形势下的电网运行需求，储能技术则是解决可再生能源大规模消纳和稳定输出、提高电力系统效率与安全性的关键技术，储能技术通过人为地释放和存储能量弥补能量在供需之间的差异性，高效利用能源，是支撑和落实能源革命的关键环节，也是实现我国“碳达峰、碳 中 和”目 标 的 关 键 技 术3-5。图 1 为 20232030年全球储能装机预测值，可以看出，储能行业发展迎来重大机遇，2030年全球储能装机预测将达到150 GWh左右6。国内外大规模新型储能技术研发与示范应用项目陆续启动，储能技术主要包括液流电池、钠离子电池、锂离子电池、压缩空气储能、铅碳电池、储热技术等3，7-8。其中压缩空气储能(CAES)具有效率高、容量大、寿命长、成本相对低廉、对环境友好等优点，被认为是最具有发展前景的大规模储能技术之一9-10。图 2 为美国能源部(DOE)对 20232030年全球CAES系统装机预测值6，可以看出，2030年CAES系统装机预测将达到60 GWh左右，约占全球储能装机的40%；CAES将快速发展，成为大规模储能技术的主力军之一6。截至目前，世界上仅有德国 Huntorf 和美国McIntosh两座投产的商业化传统CAES电站10，由于传统CAES系统存在依赖化石燃料、效率低、能量密度低等缺点，近年来形成了蓄热式压缩空气储能(TS-CAES)、等温压缩空气储能、液态空气储能、超临界压缩空气储能、水下压缩空气储能、湿空气透平压缩空气储能以及CAES耦合系统等多种新型CAES系统10。以上新型CAES系统中，TS-CAES 系统的技术最为成熟，具有不消耗化石燃图2全球CAES装机预测(DOE)Fig.2Global CAES installed capacity forecast(DOE)图1全球储能装机预测(DOE)Fig.1Forecast of global energy storage installations(DOE)18412023 年第 12 卷储能科学与技术料、流程简单、效率高等优点，受到国内外学者的广泛关注。中国科学院工程热物理研究所研发的山东省泰安市肥城10 MW盐穴CAES电站于2021年9月并网发电成功11、张家口市张北县的国际首套100 MW先进CAES国家示范项目已建设安装完成，并于2022年9月顺利实现发电并网12，TS-CAES系统成为最具有产业化前景的大规模CAES储能技术。CAES系统主要包括发电机、压缩机、燃烧室、储气室、膨胀机和电动机等关键部件，分为储能与释能两个过程。储能过程，利用风电、光电等可再生能源电力驱动压缩机压缩空气，将高压空气存储于储气室中；释能过程，储气室中的高压空气驱动膨胀机做功进行发电。CAES系统的运行过程中，由于可再生能源电力输入波动、储气室压力变化、负荷调节需要以及环境温度、压力变化等原因，其一直处于变工况条件下运行，如图3所示13-14。通常CAES系统为定容储气，储能过程中储气室内的压力不断升高；释能过程中储气室内的压力不断降低。目前，多采用节流阀节流的方式稳定压缩机出口压力和膨胀机入口压力，即定压运行方式。该方式使压缩机和膨胀机在接近设计条件下运行，但节流阀会产生节流损失，可能会降低系统效率。因此为减少节流损失，滑压运行方式开始获得广泛关注。郭欢等15对比了定压和滑压运行策略对 TS-CAES系统储/释能变工况特性和质量流量、功率、各级压比/膨胀比等关键参数的影响，结果表明：与定压运行相比，滑压运行可使系统效率提高2.08%。Han等16考虑了压缩和膨胀过程压力调控的不同模式，对比了定压-定压、定压-滑压、定压-滑压、滑压-滑压运行模式下的储释能特性。在4种运行模式中，滑压-滑压运行模式的系统效率和储能密度最大。但以往研究仅针对同一储/释能时间开展不同运行模式对比分析，且未关注在同一过程(如储能过程或释能过程)定压与滑压结合的运行模式对CAES系统性能的影响，以及未考虑动态运行调控对部件性能的影响。鉴于CAES系统变工况特性研究的重要性及对TS-CAES系统变工况特性研究不足的问题，本工作通过对CAES系统储、释能过程变压力工作过程的仿真分析与性能计算，能够获得CAES在变工况条件下的系统效率，同时，也可以获得不同工作参数下，各个部件的工作性能，指导部件的设计与调节。1 蓄热式压缩空气储能系统工作原理TS-CAES系统流程图如图4所示，其工作原理为：储能过程中电动机驱动压缩机将空气压缩，储存在储气室中，同时，来自冷罐的换热流体经间冷器将压缩过程产生的压缩热吸收并储存到热罐，从而将电能转化为压力能和热能；释能过程，储气室的高压空气释放，再被热罐的换热流体经再热器加热后进入膨胀机做功发电，将压力能和热能转化为电能，同时，被高压空气冷却后的换热流体经后冷器散热后，流入冷罐。2 蓄热式压缩空气储能系统热力学模型2.1压缩机离心式压缩机变工况运行时，其压比、效率的通用特性曲线可以近似表示为相对折合质量流量和相对折合转速的函数17：c=c1(n c)G2c+c2(n c)Gc+c3(n c)(1)c=1-c4(1-n c)2(n c/Gc)(2-n c/Gc)(2)c1=n c/a(1-b/n c)+n c(n c-b)2(3)GM输入功率波动压力上升(储能)压力下降(释能)负荷需求波动压缩机膨胀机发电机电动机储气室电网燃烧室风能太阳能燃料图3CAES系统动态运行示意图Fig.3Dynamic operation diagram of CAES system1842第 6 期孙晓霞等：压缩空气储能系统动态运行特性c2=(a-2bn 2c)/a(1-b/n c)+n c(n c-b)2(4)c3=-(abn c-b2n 2c)/a(1-b/n c)+n c(n c-b)2(5)c4=0.3(6)式中，Gc为相对折合空气质量流量，n c为相对折合转速，c为相对折合压比，c为相对折合等熵效率。各级压缩机出口温度、出口压力以及压缩功分别为：Tc，out=Tc，in(1+ck-1k-1c)(7)pc，out=pc，inc(8)wc=hc，out-hc，in=cp(Tc，out-Tc，in)(9)式中，Tc,out、pc,out分别为压缩机的出口温度、出口压力；hc,in、hc,out分别为压缩机的入口比焓、出口比焓。2.2膨胀机透平膨胀机变工况运行时，其膨胀比、效率随转速、流量等参数变化而变化，根据弗留格尔公式得到透平膨胀机的通流特性公式17：Ge=1.4-0.4n e()1/2e-1/()1/2e0-1(10)e=1-t4(1-n e)2(n e/Ge)(2-n e/Ge)(11)t4=0.3(12)式中，e、e0分别为膨胀比、设计膨胀比，Ge为相对折合空气质量流量，n e为相对折合转速，e为相对折合等熵效率。各级膨胀机出口温度、出口压力、膨胀功以及等熵膨胀功分别为：Te，out=Te，in 1-(1-k-1ke)e(13)pe，out=pe，in/e(14)we=he，in-he，out=cp(Te，in-Te，out)(15)we,iso=we/e(16)式中，Te,out、pe,out分别为膨胀机的出口温度、出口压力，he,in、he,out分别为膨胀机的入口比焓、出口比焓，e为膨胀机的效率。因此实际膨胀过程的膨胀效率定义为：=1t14we，i/1t14we，iso，i(17)式中，i、t分别代表膨胀机级数和释能时间。2.3换热器根据能量守恒，空气释放的热量和循环水吸收的热量相等：m hot(hin，hot-hout，hot)=m cold(hout，cold-hin，cold)(18)式中，下标cold和hot分别表示冷流体和热流体，in和out分别表示进口和出口。采用效能-传热单元数(-NTU)法建立换热器变工况数学模型：=()m cp|Tin-Touthot()m cpmin()Tin，hot-Tin，cold=()m cp|Tout-Tincold()m cpmin()Tin，hot-Tin，cold(19)(m cp)min=min(m cp)hot，(m cp)cold(20)=()m cpmin()m cpmax(21)=1-exp-NTU()1-1-exp-NTU()1-(22)NTU=UhAh()m cpmin(23)压缩机1压缩机2压缩机3压缩机4膨胀机1膨胀机2膨胀机3膨胀机4后冷器后冷器再热器1再热器2再热器3再热器4间冷器1间冷器4间冷器3间冷器2冷罐热罐压缩机5压缩机6间冷器5间冷器6储气室空气空气图4TS-CAES系统流程图Fig.4TS-CAES system flowchart18432023 年第 12 卷储能科学与技术式中，为换热器效能，NTU为换热单元数，是反映换热器结构和流体最小热容率的物理量，mcp为流体比热容，为两种流体比热容的比值，Uh为换热器传热系数，Ah为换热器传热面积。当两种流体比热容相等时，换热器效能表示为17：=NTU1+NTU(24)由于换热器的面积和空间有限，空气通过换热器时，存在不可忽视的流动阻力，采用以下压损模型：p/pin()p/pind=()m T/p2in()m T/p2in，d(25)式中，T表示温度；p表示压力；p为空气经间冷器/再热器的压力损失；下标in表示入口参数，d表示设计值。从各级换热器出来的水混合满足能量守恒：i=1Nmihi，in=hin，tanki=1Nmi(26)式中，mi为各级换热器的循环水流量；hi,in为各级换热器出口循环水的比焓；hin,tank为热罐中热水的比焓。2.4储气室储气室与压缩机组/膨胀机组之间的物质流动和能量流动关系，属于非稳定的热力学开口系统问题。根据质量守恒定理和能量守恒定理可得：dmdt=m in-m out(27)d()mudt=m inhin-m outhout-KwAw(T-Tenv)(28)式中，m in、m out为储气室充、放气时的空气质量流量，m为储气室的空气质量，t为时间，u为空气的内能，h为空气的比焓，T为储气室内空气温度，Tenv为环境温度，Kw为储气室内空气与环境的换热系数，Aw为储气室的换热表面积。进一步结合理想气体微分方程关系式：dpp+dVV=dTT+dmm(29)可得，储气室内部空气压力、温度随时间变化的微分方程式：dTdt=1m(m incpcvTin+(1-cpcv)m out-m in T)-KwAwmcv(T-Tenv)(30)dpdt=RVcpcv(m inTin-m outT)-KwAwVRcv(T-Tenv)(31)2.5系统性能评价指标采用储能效率与储能密度来衡量先进CAES系统能量转化利用情况，储能效率为系统总膨胀功与总压缩功的比值，衡量系统的能量转换情况；储能密度为系统总膨胀功与储气室体积的比值，衡量系统在单位储气室体积下的做功能力。ESE=WeWc=0tdischargewedt0tchargewcdt(32)ESD=WeV=0tdischargewedtV(33)3 结果与讨论本工作以功率等级为10 MW的TS-CAES系统作为研究对象，六级压缩、四级膨胀，采用定容储气室储气，以及主动蓄热方式进行蓄热，蓄热介质为水，储能过程换热器为间冷器进行级间冷却，释能过程换热器为再热器进行级间再热。TS-CAES系统设计参数如表12所示。不同释能时间下储气室容积的计算流程如图5所示，以满足储气室最终压力和释能时间为条件来确定储气室容积。3.1储能过程3.1.1压缩机定压运行储能过程采取定压运行策略，即在压缩机组与储气室之间通过安装调节阀抵消由于储气室充气导致的背压升高，此种运行方式使得压缩机组在设计工况下运行，系统稳定性增强。储气室初始压力为7 MPa、初始温度为298 K，终止条件为储气室压力达到设计值，最终获得不同释能时间下储气室空气压力、温度、质量、空气入口流量和总功率等参数随时间的变化规律。图 6 为在不同释能时间下储能过程储气室压表1TS-CAES系统设计参数Table 1Design parameters for the TS-CAES system指标环境温度环境压力储气室体积释能时间储气室最高压力数值298 K0.1 MPa11200 m3/16500 m3/21700 m3/26800 m34 h/6 h/8 h/10 h10 MPa1844第 6 期孙晓霞等：压缩空气储能系统动态运行特性力、温度与质量变化曲线，可以看出，储气室压力和质量随时间线性增加，温度随时间缓慢增加最终趋于平缓。原因是储能过程中空气不断流入储气室，且由于调节阀的存在，空气质量流量恒定，随着储能过程的进行，定容储气室内的空气被压缩，压力和内能均增加，同时储气室内外温差增大，换热量增加，导致储气室温度最终将趋于稳定，而压力和质量受流量影响较为显著，都近似于线性增长。不同释能时间并不改变储气室压力、温度和质量的变化趋势，只是由于工作时间不同，变化斜率不同。图7为释能10 h时定压运行下的质量流量和压缩机总功率变化曲线，可以看出，由于定压运行使用节流阀，从而使得空气质量流量恒定并使压缩机入口压力恒为10 MPa，所以压缩机总功率也保持设计值10 MW不变。3.1.2压缩机滑压运行储能过程采取滑压运行策略，即末级压缩机出口压力随储气室压力变化而变化。储气室初始压力为7 MPa、初始温度为298 K，终止条件为储气室压力达到设计值，最终获得不同释能时间下储气过程中储气室压力、温度、质量、压缩机压比、等熵效率和总功率等参数随时间的变化规律。图8为储能过程储气室压力、温度和质量随时间的变化曲线，可以看出，储气室压力和质量随时间线性增加，温度随时间缓慢增加最终趋于平缓。原因是储能过程中空气不断流入储气室，储气室内空气压力、温度均上升，储气室内的空气被压缩，压力和内能均增加，同时储气室内外温差增大，换热量增加，导致储气室温度最终将趋于稳定。不同释能时间并不改变储气室压力、温度和质量的变化趋势，只是由于工作时间不同，变化斜率不同。图9为储能过程滑压运行时空气质量流量和压缩机总功率随时间的变化曲线，可以看出：空气质量流量随时间逐渐减小，压缩机总功率随时间先增大后减小。原因是随着储能过程的进行，储气室压力逐渐增加至设计值，压缩机背压增加，压比增大，由压缩机流量特性曲线可知，流量随压比的增大而减小；压缩机总功率主要影响因素为压缩机压比和质量流量，储气室空气压力较小时，流量大而压比小，储气室空气压力较大时，压比大而流量小，因此，压缩机总功率最大值为中间某一时刻的值。储能过程结束，储气室压力上升至 10 MPa，此时，空气质量流量和压缩机总功率均达到设计值。不同释能时间并不改变空气质量流量和压缩机总功率的数值和变化趋势，只是由于工作时间不同，曲线斜率在改变。图5储气室容积的计算流程Fig.5The flowchart for determining the volume of the air chamber表2储能过程与释能过程设计参数Table 2Design parameters for the charge and discharge process储能过程压缩机背压压缩机总功率压缩机等熵效率间冷器两端换热温差/压损空气质量流量压缩机压比热罐入口热水混合温度数值10 MPa10 MW0.845 K/0.02 MPa18.0427 kg/s2.2191389.2 K释能过程膨胀机入口压力膨胀机总功率膨胀机等熵效率再热器两端换热温差/压损空气质量流量膨胀机膨胀比再热器出口热水混合温度数值7 MPa10 MW0.885 K/0.02 MPa28.7062 kg/s2.8192302.7 K18452023 年第 12 卷储能科学与技术图10为储能过程滑压运行时各级压缩机压比随时间的变化曲线，可以看出：各级压缩机压比随时间逐渐增加至设计值，且级数越靠后的压缩机压比越小。原因是储气室压力越小于设计压力，总压比就越小；通过各级压缩机的空气质量流量相同，而各级入口压力的变化较大，导致级数越靠后的压缩机入口压力减小值越大，相对折合质量流量增加越多，压比就越小。图11为储能过程滑压运行时各级压缩机等熵效率随时间的变化曲线，可以看出：压缩机等熵效率变化趋势与压比相近。原因是入口压力变化越大，则相对折合流量增加越多，导致等熵效率减小越多，级数越大，等熵效率越低。3.2释能过程3.2.1膨胀机定压运行以储能过程结束时储气室温度为释能过程初始温度条件，释能过程采取定压运行策略，即在储气室与膨胀机组之间通过安装节流阀抵消由于储气室放气导致的压力下降，保证膨胀机入口压力为设计入口压力7 MPa，获得不同释能时间下储气室压力、温度、质量、空气入口质量流量和总功率随时间的变化规律，系统运行终止条件为储气室压力减小至7 MPa。图 12 为在不同释能时间下定压运行时储气室压力、温度与质量变化曲线，可以看出，储气室压力和质量随时间线性减少，温度随时间缓慢减少最终趋于平缓。释能过程中空气不断从储气室流出，节流阀使得空气质量流量恒定，随着释能过程的进行，储气室内压力和内能均降低，同时储气室内外温差增大，换热量增加，所以储气图6不同释能时间下储气室压力、温度和质量变化曲线(a)4 h；(b)6 h；(c)8 h；(d)10 hFig.6Pressure,temperature and mass curve of the air storage chamber in different charging time图7储能过程定压运行质量流量及膨胀机总功率变化曲线(10 h)Fig.7Variation curve of air mass flow and total output power in charging process under constant pressure(10 h)1846第 6 期孙晓霞等：压缩空气储能系统动态运行特性室温度趋于平缓。不同释能时间并不改变储气室压力、温度和质量的变化趋势，只是由于工作时间不同，变化斜率不同。图 13 为释能 10 h 时定压运行下的质量流量和膨胀机总功率变化曲线，图8不同释能时间下储气室压力、温度和质量变化曲线(a)4 h；(b)6 h；(c)8 h；(d)10 hFig.8Pressure,temperature and mass curve of the air storage chamber in different charging time图9不同释能时间下质量流量与压缩机总功率变化曲线(a)4 h；(b)6 h；(c)8 h；(d)10 hFig.9Variation curve of air mass flow and total output power in different charging time18472023 年第 12 卷储能科学与技术可以看出，由于定压运行使用节流阀，从而使得空气质量流量恒定并使膨胀机入口压力恒为7 MPa，所以膨胀机总功率也保持设计值10 MW不变。图10不同释能时间下滑压运行各级压缩机压比变化曲线(a)4 h；(b)6 h；(c)8 h；(d)10 hFig.10Variation curve of each stages compression ratio in different discharging time图11不同释能时间下滑压运行压缩机各级等熵效率变化曲线(a)4 h；(b)6 h；(c)8 h；(d)10 hFig.11Variation curve of each stages insentropic efficiency in different discharging time1848第 6 期孙晓霞等：压缩空气储能系统动态运行特性3.2.2膨胀机滑压运行释能过程采取滑压运行策略，即储气室压力为第一级膨胀机入口压力，膨胀机入口压力随储气室压力的减小而减小，获得不同释能时间下储气室压力、温度、质量、空气入口质量流量、膨胀机膨胀比、等熵效率和总功率等参数随时间的变化规律，终止条件为储气室压力减小至7 MPa。图14为不同释能时间下滑压运行时储气室压力、温度和质量随时间的变化曲线，可以看出：储气室内空气压力和质量随时间线性下降，温度随时间快速降低后再增加。原因是释能过程，空气不断流出，储气室内的压力和温度均下降，同时由于储气室内外温差逐渐增大，换热量逐渐增加，导致储气室温度减小一定时间后逐渐上升。不同释能时间并不改变储气室压力、温度和质量的变化趋势，只是由于工作时间不同，变化斜率不同。图15为释能过程滑压运行时不同释能时间下空气质量流量和膨胀机总功率变化曲线，可以看出，空气质量流量和膨胀机总功率均随时间呈下降趋势。原因是随着释能过程的进行，储气室压力逐渐下降，即第一级膨胀机入口压力降低，由膨胀机特性曲线可得，空气质量流量下降，同时，膨胀机总功率随之降低。不同释能时间并不改变空气质量流量和膨胀机总功率的数值和变化趋势，只是由于工作时间不同，变化斜率在改变。图16为不同释能时间下各级膨胀机膨胀比的变化曲线，可以看出，第一至三级膨胀机膨胀比随时间变化不太明显，第四级膨胀机膨胀比变化显著，且随时间逐渐下降。原因是前三级的膨胀机入口压力折合流量变化不大，而靠近大气压侧的膨胀图13释能过程定压运行质量流量及膨胀机总功率变化曲线(10 h)Fig.13Variation curve of air mass flow and total output power in discharging process under constant pressure(10 h)图12不同释能时间下储气室压力、温度和质量变化曲线(a)4 h；(b)6 h；(c)8 h；(d)10 hFig.12Pressure,temperature and mass curve of the air storage chamber in different discharging time18492023 年第 12 卷储能科学与技术机，其相对折合流量变化较大，所以膨胀比随时间变化明显。不同释能时间并不改变空气质量流量和膨胀机总功率的数值和变化趋势，只是由于工作时间不同，变化斜率在改变。图14不同释能时间下储气室压力、温度和质量变化曲线(a)4 h；(b)6 h；(c)8 h；(d)10 hFig.14Pressure,temperature and mass curve of the air storage chamber in different discharging time图15不同释能时间下质量流量与膨胀机总功率变化曲线Fig.15Variation curve of air mass flow and total output power in different discharging time1850第 6 期孙晓霞等：压缩空气储能系统动态运行特性图17为不同释能时间下滑压运行时各级膨胀机等熵效率随时间的变化曲线，可以看出，四级膨胀机的等熵效率随时间基本不变，但相对而言第四级膨胀机与前三级相比，变化较为明显。原因是前图16不同释能时间下滑压运行各级膨胀机膨胀比变化曲线(a)4 h；(b)6 h；(c)8 h；(d)10 hFig.16Variation curve of each stages expansion ratio in different discharging time图17不同释能时间下滑压运行膨胀机各级等熵效率变化曲线(a)4 h；(b)6 h；(c)8 h；(d)10 hFig.17Variation curve of each stages insentropic efficiency in different discharging time18512023 年第 12 卷储能科学与技术三级膨胀机膨胀比变化较小，由膨胀机特性曲线可得其等熵效率变化较小，第四级膨胀比变化大导致等熵效率变化也较大，同样，释能时间也不影响各级膨胀机效率。3.3不同运行模式对比TS-CAES系统储-释能过程(压缩机、膨胀机)分别采用定压-定压、定压-滑压、滑压-定压和滑压-滑压四种运行模式下TS-CAES系统的储能效率的对比如图 18 所示，储气室压力由 10 MPa 到7 MPa变化，可以看出：释能时间越长，储能效率越大，但相差不大，不同模式下的TS-CAES系统储能效率随时间分别增长了1.42%、1.38%、1.45%、1.40%。对比四种不同运行模式，可以看出，储能过程滑压运行和释能过程定压运行时，TS-CAES系统储能效率最高，在释能时间为10 h时，储能效率为72.11%。由上述可知，储能过程滑压运行可以提高性能，因此为进一步探究储气室压力变化范围对系统性能影响，针对释能时间6 h的TS-CAES系统进行对比分析，运行模式分为以下四种：滑压-滑 压(10 MPa-7 MPa)；滑 压-定 压(10 MPa-7 MPa)；滑压-滑压(10 MPa-4 MPa)；滑压-定压+滑压(10 MPa-7 MPa-4 MPa)，不同运行方式的系统性能对比如表3所示。由表可知，扩大储气室压力变化范围可以提高系统储能效率和储能密度，最多分别提高了 3.22%、13.53 MJ/m3。但是由于释能时间增加，会增加膨胀过程能量损失，从而降低膨胀效率。对于储气室压力在10 MPa至4 MPa范围，模式4的储能效率和储能密度最大，分别提高了 0.88%、0.31 MJ/m3。这说明在释能过程只采用滑压运行虽然可以减少节流损失，但是其也会使膨胀机无法在设计工况下工作而对系统性能产生负影响。然而，若释能过程采取定压与滑压运行相结合的模式，即在高于膨胀机设计入口压力时定压运行，低于膨胀机入口压力时滑压运行，则既可以让膨胀机在设计工况下工作，又可以减少定压节流损失，从而提高TS-CAES系统性能。4 结论蓄热式压缩空气储能(TS-CAES)系统是一种大规模、低成本、高效率的长时储能技术，将其集成于大规模可再生能源系统，能够有效提升其发电质量及并网率，支撑新型电力系统发展，促进“双碳”目标的实现。在高比例可再生能源电力系统中，为了平衡“源-荷”的不稳定性，需要TS-CAES具有较好的动态运行性能。本研究围绕TS-CAES系统的动态运行工况，开展了基于不同运行模式的系统性能研究，建立了详细的关键部件模型，以10 MW TS-CAES系统为对象，揭示了储能过程定压运行和滑压运行不同释能时间的储气室内气体状态变化规律，滑压运行模式不同释能时间下各级压比和各级压缩效率的变化曲线；释能过程定压运行和滑压运行不同释能时间的储气室内气体状态变化规律，滑压运行模式不同释能时间下各级膨胀比和各级膨胀效率的变化曲线。探究了四种运行模式：定-定、定-滑、滑-定和滑-滑对不同释能时间下系统性能的影响，进而图18不同运行模式下TS-CAES系统的储能效率Fig.18Energy storage efficiency of the TS-CAES system in different operating modes表3不同运行方式下系统性能对比Table 3Performance comparison of system in two different operation modes运行模式模式1模式2模式3模式4储能时间/s30769307696172961729释能时间/s17873218054811552208储能效率/%70.7571.5473.1073.98膨胀效率/%88.0087.9984.0587.97储能密度/(MJ/m3)12.9613.2226.1826.491852第 6 期孙晓霞等：压缩空气储能系统动态运行特性针对释能时间 6 h 的 TS-CAES 系统提出释能过程滑压和定压结合的运行模式。研究结果表明：TS-CAES系统在不同释能时间下均是滑-定运行下的储能效