做综合能源系统优化,最麻烦的不是设备建模,而是“不确定”三个字。风电、光伏出力说变就变,绿证价格和碳配额又给运行成本叠加了一层隐形成本。我最近把一个面向绿证-碳交易机制的综合能源系统两阶段鲁棒优化方法用Python完整实现了一遍,求解核心是C&CG(列与约束生成)算法。这套代码能处理风光出力的区间不确定性,在保证系统安全约束的前提下给出最经济的日前调度方案,同时把碳排放履约成本和绿证消纳权重都折算进优化目标。如果你正在做低碳园区、多能互补系统,或者想在电力市场环境下做鲁棒调度,这篇内容应该能帮你省掉不少折腾时间。
这类项目网上能搜到不少源码,但很多版本不是对偶变换符号写反,就是大M取值离谱导致结果完全不可信。所以这篇文章我会把建模思路、CCG迭代逻辑、Python实现细节和踩坑记录全部摊开来讲,尽量让拿到代码的人不仅能跑通,还能看懂为什么这么写。
1. 项目场景与建模对象拆解
1.1 综合能源系统到底在优化什么
我以一个典型的园区级综合能源系统(IES)作为算例。系统外部接入电网和天然气网,内部包含光伏、风电、热电联产机组(CHP)、燃气锅炉、电锅炉、电储能和热储能,同时供应电负荷和热负荷。优化目标很直白:在满足24小时冷热电平衡和设备运行约束的前提下,确定各设备的出力计划、购电购气量、储能充放策略,以及绿证购买量和碳配额履约量,使总运行成本最低。
这套结构是综合能源系统优化里最经典的配置。CHP同时产电和产热,能显著提高天然气利用效率,但它的电出力会受热出力约束,存在“以热定电”的耦合关系。燃气锅炉作为热力备用,响应速度快但效率不如CHP。电锅炉的价值在于消纳风电或光伏高峰时段的多余电力,把电转换成热,相当于给系统增加了一条“电转热”的灵活性通道。储能的作用更直白,低谷时段蓄能、高峰时段放能,平抑负荷波动和新能源的不确定性。
这里要注意,IES优化和普通电力系统优化的关键差异在于“多能互补”四个字。电能、热能、天然气在每个时段都要同时满足平衡方程,而且设备之间的出力还会互相制约。比如CHP多产电,则余热必然增加;燃气锅炉多产热,则购气量必然上升。这些耦合关系在建模时要全部写进约束矩阵,否则结果会出现“同时买高价电和烧气产热”这种不合常理的情况。
1.2 绿证和碳成本是怎么叠进目标函数的
绿证-碳交易机制是最近几年IES调度的热点约束,但如果只看论文里的目标函数,其实一点都不复杂。碳排放成本的核心是计算实际排放量与免费配额之间的差额,再乘上碳价。在日内调度模型中,我按常用做法把年度配额拆到当日,得到逐时段或整日的免费配额量。
外购电力对应的间接排放等于购电量乘电网排放因子,天然气燃烧对应的直接排放等于购气量乘天然气排放因子。两者相加得到总排放量,再减去免费配额,差额乘碳价就是碳履约成本。如果实际排放低于配额,差额为负,相当于
