企业官网建设流程全解析

研究超临界二氧化碳压缩储能系统的人越来越多，但真正把文献里的理论搬到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 PRO，Ebsilon的界面逻辑会稍微有点不同——它更强调可视化的流程图连接关系，而不是设备树式的输入方式。这个特点对做系统方案的人很友好，因为你随时能从画面上看到工质流向和状态点变化。

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质量流量和总储气量。

这里把简单循环的设计工况汇总表列出来，方便对照复现：

4. 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）打开，逐项核对哪个流股的数值偏离了物理常识。很多所谓的不收敛，本质上是某个初值或者边界条件给错了。这套“先查物理、再调数值”的顺序，帮我省下的时间远超想象。