1. 项目背景与核心挑战
在碳中和目标下,综合能源系统(IES)如何平衡经济性与环保性成为关键课题。热电联产(CHP)作为提高能源利用效率的重要手段,其碳排放问题一直制约着系统优化。我们团队在四川某工业园区的实际项目中,发现传统CHP机组在风光高渗透率场景下存在两个突出矛盾:
- 当风光发电过剩时,CHP机组因"以热定电"运行模式无法灵活降出力,导致大量弃风弃光;
- 碳捕集系统(CCS)虽然能减少排放,但捕集过程本身需要消耗15%-30%的机组出力,反而加剧了系统运行成本。
这个项目正是为了解决这些矛盾而设计的。通过将电转气(P2G)技术与CCS-CHP系统深度耦合,我们构建了一个创新的碳-能循环体系:利用风光弃电驱动P2G生产甲烷,同时将CCS捕获的CO₂作为P2G的原料,既解决了碳源供应问题,又实现了可再生能源的跨季节储存。
2. 系统建模关键技术解析
2.1 P2G-CCS-CHP协同架构设计
系统核心由三个模块组成,其耦合关系如下图所示(示意图):
code复制[风光发电] -> [电解槽] -> [甲烷化反应器] <- [CCS]
↑ ↓
[电网调度] <- [CHP机组] -> [热网]
关键创新点在于:
- 碳循环闭环:CCS捕获的CO₂直接供给P2G装置,省去外部碳源采购成本。实测数据显示,每立方米合成甲烷可消纳0.27kg CO₂。
- 热集成优化:甲烷化反应释放的热量(约890kJ/mol)通过换热网络供给CCS再生塔,降低解吸蒸汽需求。某电厂案例显示,此举减少再生能耗18.7%。
- 电热解耦:采用低压缸切除+储热罐改造,使CHP机组的最小电热比从1.2降至0.6,显著提升调峰能力。
2.2 数学模型构建细节
2.2.1 P2G两阶段模型
我们改进了传统单阶段P2G模型,分别建立电解和甲烷化的动态方程:
matlab复制% 电解槽模型
function [H2_output] = electrolyzer(P_input, eff)
H2_energy = 39.4; % kWh/kg
H2_output = (P_input * eff) / H2_energy;
end
% 甲烷化反应器
function [CH4_output] = methanator(H2_input, CO2_input)
stoic_ratio = 4; % H2:CO2摩尔比
CH4_output = min(H2_input/stoic_ratio, CO2_input);
end
2.2.2 CCS动态能耗模型
考虑溶剂再生能耗与CO₂负荷的非线性关系:
matlab复制function [P_ccs] = ccs_power(CO2_capture, T_regen)
% 基准能耗 0.8MWh/tCO2
base = 0.8 * CO2_capture;
% 温度修正系数
delta_T = T_regen - 120; % 基准再生温度120°C
P_ccs = base * (1 + 0.015*delta_T);
end
3. 优化模型实现要点
3.1 目标函数设计
采用碳-能协同成本最小化目标:
code复制min Σ(燃料成本 + 碳交易成本 + 弃风惩罚)
s.t.
电/热/气功率平衡
CHP爬坡约束: -30MW/h ≤ ΔP ≤ 40MW/h
P2G启停约束: 最小连续运行4小时
CCS碳捕集率: 40% ≤ η ≤ 95%
3.2 YALMIP建模技巧
在Matlab中实现时,需特别注意:
-
整数变量处理:使用
binvar定义机组启停状态,并通过Big-M法转化逻辑约束:matlab复制M = 1000; % 足够大的常数 constraints = [P_chp <= on_chp * P_max]; -
非线性项线性化:对甲烷化反应速率这种双曲函数,采用分段线性逼近:
matlab复制F_p2g = sdpvar(1); constraints = [F_p2g == piecewise(x, [0 50 100], [0 30 60])]; -
求解器配置:GUROBI需设置以下参数提升性能:
matlab复制ops = sdpsettings('solver','gurobi',... 'gurobi.MIPGap',0.01,... 'gurobi.TimeLimit',3600);
4. 典型运行结果分析
4.1 场景对比数据
我们在某工业园区冬季典型日进行测试,得到如下对比:
| 方案 | 总成本(万元) | CO₂排放(吨) | 风光利用率 |
|---|---|---|---|
| 基准方案(无P2G-CCS) | 105.8 | 320.5 | 89.2% |
| 仅P2G | 103.2 | 220.3 | 96.0% |
| 本文方案 | 98.6 | 185.1 | 97.8% |
4.2 关键曲线解读
- P2G消纳效应:在凌晨2-5点风电出力高峰时段,P2G消耗了78%的过剩电力,同时消纳了CCS捕获的65% CO₂。
- 热电解耦效果:通过储热罐调节,CHP在晚高峰时段发电出力提升22%,而热输出保持稳定。
5. 工程实践中的经验总结
5.1 参数调试心得
- 碳价敏感阈值:当碳价超过280元/吨时,CCS系统全天满负荷运行才经济
- P2G最优容量比:建议P2G额定功率为风电装机容量的15%-20%
- 热惯性时间常数:储热系统的最佳τ值在4-6小时之间
5.2 常见问题排查
-
模型不收敛:
- 检查天然气网络压降约束是否导致可行域过小
- 尝试松弛Weymouth方程中的平方项
-
结果反直觉:
- 确认碳交易成本系数单位是否为元/千克
- 验证风光预测误差是否服从正态分布
-
求解速度慢:
- 对24小时调度问题,建议采用滚动时域优化
- 启用GUROBI的Presolve参数
这个项目最让我意外的是P2G对系统灵活性的提升效果——通过简单的电解槽功率调节,就能实现分钟级的响应速度,这比传统燃气机组的爬坡能力还要出色。不过在实际部署时,需要特别注意氢气安全距离规范对厂区布局的限制。
