研究超临界二氧化碳压缩储能系统的人越来越多但真正把文献里的理论搬到Ebsilon 13.02里跑通的人并不多。这篇综述式的项目笔记我围绕简单循环方案做了两件事一是把sCO2储能的核心技术逻辑梳理清楚二是在Ebsilon 13.02里把充电/放电回路完整建模跑通了。不管你是刚接触这个方向的研究生还是想评估新型储能方案的工程师这篇内容应该都能帮你少走不少弯路。1. 超临界二氧化碳压缩储能为什么值得关注1.1 从物理原理看sCO2储能的底层逻辑先说为什么选二氧化碳。CO2的临界点只有7.38 MPa和31.1°C相比水的22.1 MPa和374°C或者空气的临界点约-140°C、3.77 MPa这个参数区间在工程上非常好实现。这意味着我们用常规的压缩机、换热器、储罐材料就能让CO2进入超临界态不需要极端的高温高压条件。超临界态下CO2的密度接近液体粘度又接近气体扩散系数高换热特性好。在压缩储能这个场景里高密度直接带来一个核心收益储罐体积可以大幅缩小。同样是储能用压缩气体sCO2方案的储能密度通常能达到常规压缩空气储能CAES的3到5倍。这条特性让sCO2储能成为长时储能赛道上非常有竞争力的候选技术。从热力学角度看压缩储能本质上做的是“电能→压力能→电能”的转换。充电阶段用电驱动压缩机把低压CO2压缩到高压并储存起来放电阶段高压CO2经过透平膨胀做功带动发电机输出电能。整个过程的核心评价指标是往返效率RTE也就是放电输出电能与充电输入电能的比值。你可能会问这跟常规CAES有什么区别核心区别在于CO2的工作温度区间和储存状态。CAES在储气前通常需要级间冷却储气过程大量热能散失放电时还需要补燃或外部热源加热而sCO2系统可以利用CO2在近临界区物性剧烈变化的特性配合回热和冷热能利用实现更高的系统集成度。即便是不带回热的简单循环因为CO2在超临界态的高密度特性整个系统的占地面积和储罐成本也比CAES更有竞争力。1.2 简单循环的技术定位对比带回热/再热的复杂方案在做文献调研的时候我明显感觉到简单循环在实际项目讨论中经常被低估。很多人一上来就谈“超临界CO2储能带再热、间冷、回热”的复杂构型忽略了简单循环在方案论证阶段的不可替代作用。简单循环的定义很直白充电回路是压缩机冷却器高压储罐放电回路是高压储罐透平发电机中间没有回热器、再热器等附加换热设备。它的优点是系统结构清晰、设备种类少、投资估算容易、控制逻辑简单。在很多预可研阶段简单循环的结果就是项目决策的第一版依据。复杂构型当然效率更高。带回热器regenerator的系统透平出口的高温CO2余热会被用来预热压缩机入口附近的工质这样能显著降低换热器外部热源需求把往返效率从简单循环的50%左右推到65%以上。如果再叠加再热和间冷效率甚至能逼近75%。但代价是换热网络复杂、压损控制难度大、设备投资和占地面积显著增加。我在这个项目里选择先做简单循环就是遵循“先跑通、再优化”的原则。文献综述阶段看到很多研究也是这个思路比如部分文献先给出简单循环的基准参数再以此为参照评估不同改进构型的增益。简单循环的价值不在于效率最优而在于它是一个稳定、可靠的热力学基准能把系统底层的能量转换逻辑讲清楚。2. Ebsilon 13.02选型理由与建模能力盘点2.1 Ebsilon为什么适合做sCO2系统仿真Ebsilon Professional是电站和动力系统领域非常成熟的热力仿真平台我用的13.02版本在电厂热力系统建模中依然有大量存量用户。它的核心优势在两点一是模块化组件的建模方式贴近工程习惯二是内置的流体物性库对超临界流体支持得比较完整。在Ebsilon里压缩机、透平、换热器、储罐、泵、阀门这些设备都是现成的组件直接拖到画布上连线就能搭出系统拓扑。这个体验比从零写MATLAB热力学方程快得多也比用Aspen Plus做动态模拟更贴近“电厂热平衡计算”的语境。Ebsilon的计算引擎擅长稳态工况点求解恰好匹配储能系统充放电两个阶段的稳态设计分析。13.02这个版本在物性计算上的一个特点是支持自定义流体和第三方物性接口。标准库里的CO2物性数据覆盖了从亚临界到超临界的完整区间在临界点附近的计算精度也可以接受。做设计工况计算时Ebsilon的收敛速度比较快配合灵敏度分析parameter variation功能可以在几分钟内扫描完一组压比或温度参数。如果你之前用过Thermoflow或者GT PROEbsilon的界面逻辑会稍微有点不同——它更强调可视化的流程图连接关系而不是设备树式的输入方式。这个特点对做系统方案的人很友好因为你随时能从画面上看到工质流向和状态点变化。2.2 工质物性处理CO2近临界区的建模关键sCO2系统建模最大的难点不在设备而在物性。CO2在临界点附近的热物性变化极其剧烈定压比热容cp会出现尖峰密度随压力和温度的微小变化呈现非线性剧变。这在Ebsilon里会导致两种结果物性计算慢、迭代容易发散。我在建模前先做了一步标定工作在Ebsilon里单独建了一个“流动物性查看”的小模型只放一个管道组件定义流体为CO2然后手动输入不同的压力温度组合对比Ebsilon计算的密度、焓、熵与REFPROP数据的偏差。实测下来13.02版本在7.5~8.5 MPa、30~40°C区间内的密度误差在1%以内焓值误差在0.5%以内这个精度对系统级的效率评估完全够用。这里有个实操建议尽量把CO2流体定义成独立工质流stream而不是借用“空气自定义组分”的方式来拼凑。Ebsilon对纯工质CO2的处理会调用专门的物性子程序收敛稳定性和速度都优于混合气体模型。尤其是做临界点附近的参数扫描时混合气体模型很容易出现反转温度附近迭代失败。另外如果你研究的是跨临界循环即运行范围横跨临界点一定要在流体定义中勾选“包含两相区”选项。这个选项默认可能是关闭的如果不打开Ebsilon在计算透平出口发生部分冷凝的情况时会直接报错或者给出非物理解。这也是我早期踩过的一个坑后面在问题章节会细说。3. 简单循环系统设计与热力学参数整定3.1 系统流程拓扑与充放电回路说明这个项目的简单循环拓扑我最终定为充电侧由低压储罐2.5 MPa初始压力经过压缩机升压到20 MPa高温CO2经冷却器降温至接近环境温度后进入高压储罐储存放电侧高压CO2经调节阀节流到透平入口压力这里保留了微调空间进入透平膨胀做功后低压排气返回低压储罐完成循环。这里补充说明为什么充放电回路要共用一个低压储罐而不是直接向环境排放。CO2不是空气直接排放不仅浪费工质而且涉及捕集和成本问题。闭式循环是sCO2储能的主流方案两个储罐之间来回倒换工质系统总工质装载量固定运行过程中不需要从外界补充CO2。从文献综述的结果看国内外团队采用的储压范围差异不小。有的方案储压上限定在12 MPa左右属于“中压储能”也有研究把压力推到25 MPa以上追求高储能密度。我选择20 MPa是综合考虑了压缩机能耗、储罐壁厚成本和透平膨胀比三方面因素这个数值在大多数已发表的简单循环文献中也处于合理区间。放电回路有一个细节值得注意透平入口需要保持超临界态。如果储罐压力下降过多经过调节阀后压力可能跌到临界点以下透平入口出现两相流这会严重威胁叶片安全。所以在设计中我设置了最低放电压力限值比如14 MPa低于这个值就停止放电宁可牺牲一部分可用能量也要保证设备安全。3.2 关键参数的计算与选取逻辑先看压缩机。Ebsilon里压缩机的输入参数包括入口压力、入口温度、排气压力、绝热效率isentropic efficiency或等温效率。储能系统压缩机一般选用绝热效率作为输入取值0.82~0.87之间比较合理我这里取0.85。根据绝热压缩公式算出来的排气温度大约在155°C左右这个温度对常规碳钢管道和阀门没有问题。关键是冷却器的端差设置。简单循环中没有回热器所以冷却器的任务就是把压缩机排气温度降到接近环境温度减少高温工质在储罐内储存时的散热损失。我在Ebsilon里把冷却器出口温度设定为环境温度10°C即端差10°C环境温度取25°C所以冷却器出口为35°C。这个参数直接影响充电过程的等效电耗冷却温度越低储罐内CO2密度越大同样的储罐容积能存更多工质但同时需要更大的冷却面积和冷却水量。透平侧我取绝热效率0.88。透平进口参数是决定发电量的关键入口温度越高、压比越大理论焓降越大。但因为简单循环没有回热器透平进口温度就取决于储罐内的温度。由于高压储罐在充电末期会积累一定的压缩热虽然经过冷却器但难以完全冷却到环境温度所以放电初期透平入口温度我按40°C代入计算这算一个比较保守的取值。储能容量按10 MW/40 MWh设计即满功率放电4小时。这个容量级别对应的是典型的中小型储能电站换算下来高压储罐体积约在2000~4000 m³取决于储压和温度。Ebsilon里直接给定高压储罐的终态压力和温度就能反算需要的CO2质量流量和总储气量。这里把简单循环的设计工况汇总表列出来方便对照复现参数充电工况放电工况低压储罐压力2.5 MPa2.5 MPa终点高压储罐压力20 MPa14 MPa终点压缩机入口温度25°C-压缩机排气压力20 MPa-压缩机绝热效率0.85-冷却器出口温度35°C-透平入口温度-40°C透平绝热效率-0.88放电时长-4 h系统额定功率-10 MW4. Ebsilon建模实操过程全记录4.1 新建项目与组件拖放从零搭建循环打开Ebsilon 13.02后在“File”菜单新建一个案例case命名好项目路径。我的建议是直接把充放电两个回路建在同一个画布上中间用储罐组件作为连接点。这样做的好处是后续做参数扫描时能同时观察充放电两侧的响应。组件清单如下按Ebsilon默认组件库名称管道Pipe、压缩机Compressor、冷却器Cooler也可以选Heat Exchanger、储罐Tank、透平Turbine、发电机Generator、调节阀Valve、分流/合流器Mixer/Splitter。搭建顺序我习惯按工质流向走先搭充电回路从低压储罐出口拉管道到压缩机入口压缩机出口接冷却器热侧入口冷却器出口接高压储罐入口。放电回路再单独拉高压储罐出口经调节阀到透平入口透平出口接低压储罐入口。连线时Ebsilon会自动识别管口的流动方向如果不小心接反了计算时会提示反向流动错误。我的经验是先给每个组件编好名称前缀C-、T-、HE-等再连线路这样排查错误时能快速定位是哪个组件出了问题。这里特别提醒一个储罐的建模细节Ebsilon的Tank组件有“固定容积”和“固定压力”两种模式。储能系统应该选固定容积模式并且设定初始充装量。如果你误选了固定压力模式储罐就变成了一个无限大的“压力源塔”无法反映充放电过程中储罐压力随工质质量变化的真实物理过程。这个设定我在第一版模型里就选错了导致整个系统能量不平衡。4.2 设计工况计算的收敛调试模型搭完之后按F9或者点击计算按钮进入求解。Ebsilon的求解顺序是迭代式的先假设一组未知量初值根据热平衡方程计算各组件进出口参数再对比残差修正。如果初值给得太离谱就会出现发散。我给初值的经验是所有流体流股先给定一个合理的压力温度猜测值比如均取12 MPa、50°C然后开启“initial values”辅助功能让Ebsilon自己迭代。千万不要让某个流股的初值处在临界点附近比如8 MPa、32°C因为临界点附近物性变化剧烈迭代初期一个微小的初值偏差就能让密度计算翻倍然后整个残差雪崩式放大。计算收敛后第一件事不是看效率而是检查每个流股的压力、温度、焓值是否在物理合理范围。我用一个笨办法把Ebsilon的计算结果导出到Excel同时打开REFPROP交叉验证几个关键状态点的密度和焓。这个步骤看起来很费时间但能帮你尽早发现物性库标定的问题。如果碰到不收敛优先检查三处一是是否有未连接的死支路二是换热器的端差和换热系数是否设置得太激进三是储罐的质量累积是否与能量平衡方程冲突。实际调试中我把冷却器端差从5°C改到10°C后收敛速度明显提升原因也简单——端差越小冷却器需要的UA值越大方程组的数值刚性越强。4.3 参数扫描与离设计工况分析设计工况跑通后我在Ebsilon里启用了“Parameter Variation”工具对压缩机排气压力和透平入口温度做了双参数扫描。这一步的目的是画出简单循环的“效率地图”观察不同运行参数组合下系统往返效率的变化趋势。扫描结果很有代表性在排气压力从12 MPa升到20 MPa的区间内系统往返效率先升后降峰值出现在16~17 MPa附近。原因并不难解释——低压端固定是2.5 MPa压比升高带来的单位工质发电量增加在初期会超过压缩机耗功的增加但当压比超过某个值后压缩过程偏离等熵压缩的损耗增长更快边际收益变为负值。这个峰值点对实际项目选压非常有意义。离设计工况off-design计算也是Ebsilon相对方便的功能。在“Case”里复制设计工况为离设计工况然后把放电功率设定为50%负荷5 MW透平采用等熵效率修正曲线或固定效率计算部分负荷下的效率衰减。我实测下来简单循环在50%负荷时往返效率会比设计工况低约3~4个百分点主要原因是透平和压缩机偏离设计压比后效率下降以及固定端差换热器在部分负荷下传热性能变差。这里有一个关于“往返效率定义”的陷阱Ebsilon默认计算的能效指标是各组件的独立效率它不会自动帮你算出一个“充电耗电vs放电发电”的整站效率。你需要自己在Excel里搭一个折算表把充电过程压缩机的总电耗含电机效率、冷却水泵电耗、控制系统的厂用电都算进去再除以放电过程的发电量扣除变压器损耗。很多文献里的往返效率之所以差距大就是因为口径不统一。5. 常见问题与排查技巧实录5.1 物性突变导致的收敛失败这是我在临界点附近做参数扫描时遇到最频繁的问题。当扫描点落进7.8~8.5 MPa、33~38°C这个区间时Ebsilon要么弹出“iteration failed”的红色警告要么计算结果出现密度跳变到异常值的现象。排查思路是这样先缩小扫描步长确认是否只是在临界点附近那一两个点发散还是整个区间都发散。如果只是个别点可以手动修正这些点的初值把上一步收敛结果的参数作为下一步的初始猜测。如果整个区间都发散那多半是流体物性模型设置的问题回到流体定义里检查是否启用了两相区选项。还有一个土办法很有效把参数扫描的步长从“压力均匀分布”改成“log分布”让临界点附近被更多扫描点覆盖其余区域少放几个点。因为sCO2物性在临界点附近对压力极其敏感线性均匀布点很容易漏掉关键转折。5.2 储罐建模带来的压力漂移问题储罐模型选成“固定容积”后出现了一个新问题充电过程中储罐压力飙升过快跟理论计算不符。排查后发现原因是储罐组件内部默认采用了等温假设也就是储罐内气体温度始终等于环境温度但实际高压CO2充入时焦耳-汤姆逊效应和压缩热的累积会让罐内温度显著升高。解决方式是在储罐组件里选择“绝热”或“多节点温度”模型要么直接给定罐体散热U值。我的做法是取一个折中给定罐体面积和综合散热系数让储罐处于近似绝热到轻微散热的区间。这样算出来的储罐终态温度比纯等温假设高出约15°C对透平入口状态影响显著。这个细节直接关系系统安全性。如果储罐实际温度高于设计值CO2密度会降低同样的质量流量下体积流量增大可能导致透平超速反之如果温度低于预期又可能提前进入两相区。所以储罐热模型的选取要谨慎不能为了省事直接用默认等温假设。5.3 效率评估口径的混乱与修正整理文献时我统计了十几篇sCO2储能的研究论文发现大家报告的往返效率从45%到75%都有初看差异大得离谱。逐一拆解计算口径后发现差异主要来自三个变量是否计入辅助设备电耗、压缩机效率取值高低、以及储能时长内储罐散热的处理方式。Ebsilon本身提供的是设备级热力学效率不是系统级往返效率。我在做这个项目综述时建立了一个统一的评价模板充电电耗 压缩机轴功 / 电机效率 冷却水泵电耗放电输出 透平轴功 × 发电机效率 - 变压器损耗系统往返效率 放电输出 / 充电电耗。采用这个统一口径后简单循环的往返效率大约在51%~54%之间跟文献中同样构型的结果对得上。如果你在写论文或做项目汇报我强烈建议在给出效率数据的同时把计算口径的明细表格附上。这个习惯能避免大量后续评审和返工——至少我是因为这个踩过坑被审稿人追着问“充电电耗里到底含不含泵功”。5.4 透平出口状态点异常的两相区问题最后一个高频问题出现在放电回路的透平计算上。当我尝试把储罐压力从20 MPa降低到12 MPa、透平入口温度又只有35°C时透平出口状态点会越过饱和线进入两相区。Ebsilon算出的透平出口焓值看起来正常但干度小于1意味着有液滴出现。解决思路有两个方向。一是提高透平入口温度通过外部加热或回热方式让膨胀过程远离饱和线二是在简单循环内控制透平膨胀比设置更高的最低储罐压力。我的模型里最终把最低放电压限定在14 MPa同时把透平入口温度最低控制在37°C以上保证透平出口干度不低于0.96。这个问题的工程意义很重要。大量文献在研究简单循环时默认“CO2永远不会液化”但实际只要膨胀比一大、入口温度不够高液滴就会产生。Ebsilon能算出来说明物性处理是到位的。反过来我也建议在审阅别人文献时留意透平出口的干度指标——如果全文没有提及干度大概率是用了理想气体假设效率数据要打个问号。写在最后做这个sCO2压缩储能简单循环的文献综述加Ebsilon仿真前后大概花了三周时间。第一周几乎全部用在啃文献和搭建Ebsilon模型上第二周调试收敛和储罐模型第三周做参数扫描和数据整理。最深的体会是这种新型储能系统的瓶颈往往不在循环思路而在物性处理的细节上。CO2在临界点附近那套“既像气又像液”的行为让所有理想化假设都必须谨慎使用。Ebsilon 13.02作为老牌电厂仿真工具做sCO2系统的稳态设计计算是够用的至少在简单循环这个阶段它的精度和收敛表现都让我满意。后续如果你想做动态仿真或者控制系统设计可能需要转向Modelica或动态流程模拟工具但稳态基准模型的搭建和维护Ebsilon还是能扛大梁。最后分享一个排查问题的小习惯每次计算不收敛先别急着调松弛因子把组件的输入输出表Input/Output List打开逐项核对哪个流股的数值偏离了物理常识。很多所谓的不收敛本质上是某个初值或者边界条件给错了。这套“先查物理、再调数值”的顺序帮我省下的时间远超想象。