共计 1829 个字符，预计需要花费 5 分钟才能阅读完成。
火电 - 储能耦合系统的价值与挑战
火电厂耦合储能系统是实现碳中和目标的关键技术路径。对于 660MW 一次再热机组而言，引入超临界 CO₂储能系统可带来三方面优势：
- 调峰能力提升 ：储能系统可存储机组低负荷时的富余电能，在用电高峰时释放，使机组始终运行在最佳效率区间。实测数据表明，耦合系统可使机组调峰范围从 50%-100% 扩展至 30%-110%
- 煤耗降低 ：通过热电解耦，机组避免频繁变负荷导致的效率损失。某 660MW 机组改造后供电煤耗下降 8.2g/kWh
- 碳排放减少 ：系统利用 CO₂作为工质，可实现碳捕集系统的无缝集成，全生命周期碳排放降低 12%-15%
但系统建模面临三大技术难点：
- 超临界 CO₂物性剧烈变化区（近临界点）的热力学计算
- 蒸汽 -CO₂换热器的跨介质传热建模
- 动态过程中机械 - 热力 - 电气的多时间尺度耦合
建模工具对比与选型
Aspen Plus
- 优势 ：
- 内置完善的物性数据库，尤其擅长复杂相态计算
- 提供成熟的电站模块库（如锅炉、汽轮机等）
- 局限 ：
- 动态仿真能力较弱
- 自定义控制逻辑开发困难
EBSILON
- 优势 ：
- 专为热力系统优化设计
- 提供标准火电机组模板
- 局限 ：
- 对新型工质支持有限
- 学习曲线陡峭
MATLAB/Simulink
- 推荐方案 ：
- 灵活构建自定义模型（推荐使用 Simscape Thermodynamics 模块）
- 强大的控制系统开发环境
- 便于与电力系统模型耦合
超临界 CO₂循环建模核心要点
1. 关键设备特性曲线拟合
压缩机等熵效率曲线采用二次多项式拟合：
\eta_{comp} = a_0 + a_1\left(\frac{N}{N_{des}}\right) + a_2\left(\frac{N}{N_{des}}\right)^2其中典型参数范围：
- 设计转速 N_des：6000-8000 rpm
- 系数 a_0：0.82-0.88
- 系数 a_1：0.15-0.25
2. 再热蒸汽 -CO₂换热器建模
采用 NTU- ε 法计算换热效率：
\epsilon = 1 - \exp\left(\frac{-NTU^{0.22}}{C_r}[\exp(-C_r NTU^{0.78})-1]\right)工程注意事项：
- 蒸汽侧压力损失需控制在 0.3MPa 以内
- CO₂出口温度与蒸汽饱和温度最小保持 15℃端差
3. 压力 - 温度协同控制
构建前馈 - 反馈复合控制系统：
% 压力控制器（示例）function dpdt = pressure_control(t,p,setpoint)
Kp = 0.8; % 比例系数（范围 0.5-1.2）Ti = 15; % 积分时间 (s)
error = setpoint - p;
dpdt = Kp*(error + 1/Ti*integral(error));
endMATLAB 建模实战示例
CO₂物性计算
% 使用 CoolProp 获取超临界 CO₂密度
import py.CoolProp.CoolProp.PropsSI;
rho = PropsSI('D','P',20e6,'T',773,'CO2'); % 典型参数：P=20MPa, T=500℃动态系统求解
% 定义微分方程组
function dxdt = system_dynamics(t,x)
% x(1): 储能罐压力 (MPa)
% x(2): 换热器 CO₂温度 (K)
dxdt = zeros(2,1);
dxdt(1) = (m_dot_in - m_dot_out)/V_tank; % 质量守恒
dxdt(2) = (Q_heat - m_dot*cp*(x(2)-T_in))/(m*cp); % 能量守恒
end
% ODE45 求解
[t,y] = ode45(@system_dynamics, [0 3600], [7e6 573]); % 初始条件：7MPa,300℃生产环境避坑指南
材料强度校核
- 在 P =25MPa、T=550℃工况下：
- 管道需采用 P92 级钢材（许用应力≥100MPa）
- 法兰密封面硬度需达 HRC40 以上
收敛性加速技巧
- 采用变步长求解器：
options = odeset('RelTol',1e-6,'InitialStep',0.1); - 对物性计算进行插值预处理
- 合理设置代数环初值
模型验证方法
- 稳态工况误差控制：
- 主要参数相对误差 <3%
- 关键节点温度绝对误差 <5℃
- 动态验证采用阶跃扰动测试
延伸思考方向
- 如何利用机器学习优化压缩机特性曲线？
- 蒸汽 -CO₂换热器能否采用微通道设计提升传热系数？
- 系统惯性时间常数与电网调频需求的匹配关系？
通过本文介绍的建模方法，工程师可在 2 - 3 周内完成基础模型搭建。实际项目中建议先构建简化模型验证控制策略，再逐步增加物理细节。期待与同行探讨更多系统优化可能性。
正文完
