做综合能源系统调度优化这几年我最大的感受是模型再漂亮落地时最头疼的往往不是设备本体而是新能源出力那股“说不准”的劲。光伏清晨爬坡、午后飘云风机一阵风过来出力猛涨、一阵风过去近乎归零如果优化模型还是拿一个固定的预测值去排计划实际运行时的功率偏差就会让调度员非常被动。这也是我为什么把“计及新能源出力不确定性的电气设备综合能源系统协同优化”作为核心课题并最终用Matlab把整套算法完整跑通的原因。这篇文章不是教科书式的理论堆砌而是一个可复现的项目经验梳理。我会从模型边界、不确定性建模、协同优化求解、Matlab代码实现、仿真结果分析和踩坑记录六个方面展开。不管你是刚接触综合能源优化的研究生还是已经在做园区级能量管理系统的工程师这篇文章应该都能帮你少走不少弯路。1. 先把模型边界说清楚这个“电气设备综合能源系统”到底优化什么标题里的“电气设备综合能源系统”我倾向于理解为广义的电气化能源设备集合也就是由电力、天然气、热力等多种能源转换与存储设备构成的园区级供能系统。最典型的成员包括光伏、风电、电储能、燃气轮机或热电联产机组、电锅炉、热储能以及电负荷、热负荷、气负荷三类用能需求。调度中心要回答的问题是未来24小时每台设备在每个时刻该发多少电、烧多少气、充多少储、放多少热才能让总运行成本最低同时满足各类负荷的供需平衡和安全约束。1.1 设备构成与决策变量的对应关系先列一个我在建模时一定会先画的表。这张表的作用不是装饰而是把“物理对象”和“数学变量”一一对应起来防止后面写约束时张冠李戴。设备核心决策量典型出力范围主要约束来源光伏/风电实际并网功率0到预测上限出力不确定区间、弃电控制电储能充电/放电功率-额定功率到额定功率SOC上下限、充放功率限、周期回位CHP热电联产机组发电功率、供热功率电出力上下限、热出力上下限热电比可行域、爬坡约束电锅炉产热功率0到额定功率电热转换效率、用电上限热储能充/放热功率-额定功率到额定功率储热罐容量约束外购电力从上级电网购入功率0到购电上限购电合同上限、功率平衡外购天然气从气网购气流量0到购气上限气网供应上限这张表列完之后整个优化问题的骨架基本就出来了。决策变量是一个多维时变序列优化维度通常是“设备数量 × 调度时段数”。如果做24小时日前调度、时间粒度取1小时那么一个中等规模园区的决策变量就有几百个再加上0-1启停变量规模达到数千维并不稀奇。1.2 为什么“协同优化”而不是各设备独立调度这是我自己体会最深的一点。很多人一开始会把问题拆成“电系统一个优化、热系统一个优化”然后单独排计划最后发现根本合不上。原因在于综合能源系统里电和热是强耦合的CHP机组发一度电的同时会带出一部分热电锅炉用电烧热又把电负荷抬高而热储能则可以把“热电联产”和“电转热”两条供热路径在时间上解耦。所以必须把电力平衡、热力平衡、燃气供应放在同一个优化框架里同时求解才能找到真正的全局最优。举个实际例子。在一个冬季典型日里如果只盯着电力成本调度方案会让你在电价低谷时段多买电、少开CHP但如果把热负荷也纳入目标你会发现在电价低谷时段用电锅炉烧热存进储热罐白天高峰时段再由储热罐放热反而比单纯调节CHP更省钱。这类跨能源品种的套利机会只有协同优化才能捕捉到。1.3 数学模型的总体框架用一个概括性的表述来说这个优化问题的标准形式是目标函数最小化总运行成本包括购电费用、购气费用、设备启停成本、运维成本以及弃新能源惩罚和失负荷惩罚。约束条件包括电功率平衡约束、热功率平衡约束、气源购气约束、各设备出力上下限、储能SOC转移方程、机组爬坡约束、外购电量上限等。我通常把这个模型写成混合整数线性规划MILP来处理。虽然原问题里有一些非线性环节例如设备效率随负载率变化但通过分段线性化之后可以在不牺牲太多精度的前提下大幅降低求解难度。这一点我在第3章还会详细展开。2. 新能源出力的不确定性建模场景法、鲁棒法与模糊机会约束的选择逻辑明确了优化问题之后接下来最关键的一步就是回答“不确定的新能源出力到底怎么放进优化模型”这一步选型直接决定了求解难度和结果的保守程度需要在建模阶段就考虑清楚。2.1 不确定性从哪来、影响有多大光伏和风电出力的不确定性主要来自预测误差。我们通常有一个预测曲线比如“明天9点到15点光伏预测功率约80%额定容量”但真实出力可能因为云层遮挡、风速变化、温度漂移而偏离预测。在确定性优化里调度员把这个预测值当作硬约束只要实际出力低于预测系统就会面临功率缺额如果实际出力高于预测又可能出现弃电或潮流越限的问题。这里要特别说一个容易被忽视的细节不确定性不只是“出力大小”的变化还涉及“时间相关性”。比如午后的光伏爬坡事件如果预测模型假设出力是平滑的而实际在10分钟内跌掉一半那么只做单一时刻的静态优化是完全不够的。所以我在建模型时习惯于把所有不确定量的时序特征保留住把场景作为一个完整的时间序列来生成而不是一时刻一时刻孤立抽样。2.2 三条技术路线的对比目前主流的建模方法可以归纳为三类我在实际项目中都试过各自的逻辑和适用场景差别很大。场景法随机规划的思路是用概率场景集合刻画可能的光伏/风电出力曲线每个场景有一个概率权重目标函数取所有场景下的期望成本。这种做法信息量最丰富结果也最贴近真实运行但计算代价随场景数量上升非常快。鲁棒优化法的思路是不关心具体概率分布只给定一个不确定集合比如出力在预测值上下浮动的区间优化目标是在最恶劣情况下的成本最小化。它的优点是对分布误差不敏感、求解相对高效缺点是非常保守有时为了应对极低概率场景会把调度成本推高10%以上。模糊机会约束法则是把约束满足度做成概率形式例如“功率平衡被满足的概率不低于95%”然后用清晰等价类把概率约束转换成确定性的形式。它比鲁棒法灵活比场景法计算量小但对分布参数的建模能力要求较高。三种方法的对比我整理成了表格方便你根据自己的数据条件做选择题。方法所需数据计算量保守程度模型复杂度工程落地难度场景法历史出力数据、预测误差分布中高中低低低鲁棒法出力上下界区间低高中低模糊机会约束预测误差分布参数中低中较高需要概率转化2.3 场景生成与削减的实操做法如果选了场景法接下来就要解决“场景怎么来”以及“场景太多怎么办”这两个问题。我推荐的做法是误差分布采样加场景削减两步走。第一步根据历史数据统计每个时段新能源出力的预测误差假设误差服从正态分布或拉普拉斯分布。用拉丁超立方采样代替普通蒙特卡洛采样可以保证覆盖概率空间的均匀性同样的样本数量下效果明显更好。第二步对生成的大规模场景集做场景削减常用手段包括同步回代削减法和聚类削减法。我自己试下来先用K-means聚类把几百个场景聚成30到50个典型场景再对每个聚类分配概率权重平衡效果和计算速度都很理想。2.4 我在代码里实际采用的方案在这个项目里我最终采用了“场景法为主、鲁棒边界为辅”的混合思路。核心框架是随机规划的场景法同时给最重要的功率平衡约束加上一个鲁棒化的备用容量项。这么做的好处是典型的调度结果兼顾了经济性和健壮性不会因为某几个极端场景而把成本撑得过高。准确说目标函数仍然是所有场景下的期望运行成本但约束里多了一条“系统上调备用和下调备用分别不低于预测误差区间”这等于给不确定性加上了一道安全底垫。3. 协同优化的数学建模与求解器配置目标函数、约束条件与线性化处理模型框架定下来后真正动手写Matlab代码前必须先把目标函数和约束条件在纸面上列清。这一步如果偷懒后面调试代码时一定会付出好几倍的代价。3.1 目标函数运行成本构成的拆解我采用的日前调度目标函数为最小化总运行成本第一项是外购电力费用电价采用分时电价曲线有的园区还有容量费或需量费我通常把容量费转化成对购电功率峰值的惩罚项一并放入目标第二项是天然气购气费用第三项是设备启停成本主要来自于CHP这一类的大型机组频繁启停不仅磨损设备还会在热力学上造成额外燃料消耗第四项是设备运维成本可以简化成出力或功率的线性函数第五项是弃新能源惩罚项用于软化“新能源出力必须全额消纳”的硬约束惩罚系数需要设置得比购电和购气成本高一些才能让优化器在几乎所有情况下都优先消纳新能源。这里有一个很容易犯的错误惩罚系数设置不当。如果弃电惩罚设置得太低优化器会在某些时段宁可通过弃电来避免调节其他设备导致新能源消纳率异常。如果设置太高又会迫使系统在极端场景下花高昂代价去消纳边际效率极低的新能源出力抬升总成本。我通常的做法是把弃电惩罚设为最高电价水平的两到三倍再根据结果做敏感性分析微调。3.2 核心约束条件的写法约束条件里电功率平衡是最基础的一条各发电设备出力加储能放电功率加外购电减去储能充电功率减去用电负荷要等于零。热功率平衡类似CHP供热加电锅炉产热加热储能放热减去热负荷和热储能充热功率等于零。设备出力上下限和爬坡约束也需要结合实际情况考虑。比如CHP机组不能简单把它当作一个可发任意功率的电源它存在一个热电耦合可行域在纯凝工况下发电效率高、供热为零在抽凝工况下发电供热同向变化。我在建模时用多边形可行域近似这个热电耦合关系用几条线性不等式把可行域包络出来效果比单一固定热电比要好得多。储能约束更是重头戏。电储能和热储能的SOC状态转移方程是典型的动态约束SOC(t1) SOC(t) 充电功率 × 充电效率 - 放电功率 / 放电效率这块必须注意效率的放置位置放错位置会直接导致“越充越亏”或“越放越赚”的荒谬结果。另外一个工程上常见的硬性要求是调度周期末SOC要回到初始值附近否则电池无法做循环往复的日常调度会持续漂移到满电或空电状态。3.3 非线性环节的线性化处理原模型中存在两类常见的非线性环节一是设备效率随负载率变化二是包含0-1整数变量的启停逻辑。第一类可以通过分段线性化逼近把每台设备的电热转换关系用多个线性段表示第二类则需要引入大M法或特殊顺序集来线性化。虽然这会增加约束数量但换来的是可以使用商用线性求解器求解速度快且解的质量有保证。我在实际代码里面用Yalmip做建模层求解器选用CPLEX或Gurobi。对于30个场景、96个时段的典型算例模型规模通常在数万约束量级求解时间在几十秒到几分钟之间完全可以满足日前调度的计算需求。3.4 求解器选型和运行环境配置如果你已经有了Matlab环境请务必安装Yalmip工具箱并接好一个可用的混合整数线性规划求解器。环境配置上有一个容易踩的坑Yalmip里指定求解器的方式是sdpsettings(solver, gurobi)如果你安装的是CPLEX要对应改成cplex。另外还要注意求解器许可证的位数必须与Matlab位数一致64位Matlab配32位求解器接口时经常出现无法加载的动态库报错。4. Matlab代码实现的运行流程与核心细节模型数学形式明确之后代码实现的难度已经降了大半。这一章我就直接分享我的代码架构以及几个真正影响运行效果的关键细节。4.1 代码模块怎么划分我的项目代码文件组织如下main.m总入口负责读取参数、调用各模块、汇总结果并绘图。load_data.m读取负荷曲线、新能源预测曲线、设备参数、电价和气价。scenario_gen.m生成新能源不确定性场景集并完成场景削减。build_model.m用Yalmip定义决策变量、目标函数和约束条件。solve_model.m调用求解器并处理求解结果。plot_results.m绘制电功率平衡、热功率平衡、储能SOC、新能源消纳率等曲线。模块划分最核心的原则就是数据生成和模型求解分离。这样做的好处非常实际当你想调整场景数量、修改预测误差参数时只需要改scenario_gen.m不需要动build_model.m里的任何约束当你想增加一个新设备时也只需要在build_model.m里追加变量和约束不用重写数据预处理逻辑。4.2 场景生成模块的代码骨架场景生成模块我用拉丁超立方采样加K-means削减的写法。下面这段代码是核心逻辑% scenario_gen.m 核心片段 Nscenes 1000; % 初始抽样场景数 K 30; % 削减后的典型场景数 T 96; % 调度时段数 mu Ppv_forecast; % 光伏预测曲线 sigma 0.1 * mu; % 标准差按预测值10%设置 % 拉丁超立方抽样 lhs_samples lhsdesign(Nscenes, T); scenes zeros(Nscenes, T); for t 1:T scenes(:, t) norminv(lhs_samples(:, t), mu(t), max(sigma(t), 0.01)); end scenes max(scenes, 0); % 出力非负 % 限制在装机容量内 scenes min(scenes, Ppv_cap); % K-means场景削减 [idx, C] kmeans(scenes, K, MaxIter, 500, Replicates, 5); prob histcounts(idx, K) / Nscenes;这段代码里有个细节值得注意标准差按预测值的比例设置因为光伏出力的预测误差通常随着出力水平变化而不是恒定绝对量。夜间光伏出力趋近于零时误差也相应变小所以我会加一个0.01的保底下限避免极端抽样值频繁出现。4.3 模型构建模块的Yalmip写法build_model.m中用Yalmip定义变量的方式是整个代码的核心。我习惯用三维变量第一维是设备编号第二维是时段第三维是场景索引。如果一个变量不随场景变化那就定义成两维。这个区分非常重要日前调度中启停计划和购电决策可以是“非预期”的也就是在不确定性实现前就要确定下来而机组实际出力、储能充放电可以随场景自适应调整。用专业术语说这是两阶段决策结构第一阶段变量与场景无关第二阶段变量依赖场景。% build_model.m 核心片段 % % 变量定义 % g(i, t) 表示设备i在t时段的出力s表示场景 g sdpvar(repmat(Ndev, 1, 1), T, K); % 依赖场景 u binvar(Ndev, T); % 0-1启停变量不依赖场景 % 目标函数场景期望成本 Objective 0; for k 1:K Objective Objective prob(k) * ... (购气成本(k) 购电成本(k) 启停成本 运维成本 弃电惩罚(k)); end % 约束条件拼接 Constraints []; for k 1:K Constraints [Constraints, 电功率平衡表达式 0]; Constraints [Constraints, 热功率平衡表达式 0]; % 设备上下限、储能SOC、爬坡、备用等约束 end实际代码中不会把所有成本项真的写成“购气成本”这种中文变量这里只是表达结构。我需要强调的一点是第一阶段变量和第二阶段变量的区分直接决定了模型是否符合“非预期性”原则。如果不做这个区分把所有变量都设置成随场景变化的求解器就会获得一种信息上的“先见之明”得到过度乐观的结果。这是随机优化里最容易错的模型设定没有之一。4.4 求解结果的处理与可行性检查求解完成之后不要急着画图先做一步基础检查查看求解器的返回状态确认问题被求解到最优或至少在可接受的gap范围内。我习惯用以下代码sol optimize(Constraints, Objective, options); if sol.problem ~ 0 warning(求解状态异常: %s, yalmiperror(sol.problem)); else g_opt value(g); u_opt value(u); end如果返回值不是0最常见的两个问题是“不可行”和“数值问题”。前者说明约束之间互相矛盾比如某时段热负荷太高而所有供热设备上限加起来都不够后者通常是量纲不一致导致的比如功率用kW能量用MWh储能SOC方程里差了个倍数。我会在下一章详细讲这类坑的排查方法。5. 仿真结果分析不确定性水平对调度策略与经济性的影响代码跑通之后最需要回答的问题是这个模型相对于确定性模型到底带来了什么价值我通过一组对照实验来说明结果差异实验中只改变新能源预测误差的标准差从5%逐步提高到15%观察成本和调度策略的变化。5.1 实验设置与经济性对比基准场景是某园区冬季典型日包含光伏100MW、风电50MW、电储能20MW/40MWh、两台50MW级CHP机组、电锅炉30MW和热储能50MWh。负荷曲线来自该园区实测数据的典型剖面电价采用峰平谷三段式结构。我在不同误差水平下分别跑确定性模型和场景随机优化模型得到的总运行成本如下误差水平确定性模型成本万元场景法模型成本万元成本增幅5%128.6129.80.9%10%128.2132.43.3%15%127.9137.67.6%这个表格里有一个反直觉的结果确定性模型的成本随着误差增大反而下降。这不是因为它变得更好而是因为确定性模型使用的是预测出力均值误差越大它就越敢减少备用、少开机组把成本“虚低”了。一旦实际场景偏离预测均值这个低成本的计划根本没有可执行性实际运行会发生大量切负荷或弃风弃光。换句话说确定性模型的成本是“数据偏差掩盖出来的好消息”而场景法的成本才是真正能落地的期望值。5.2 调度策略的差异在哪进一步看设备出力曲线有两处规律非常明显。第一不确定性越大CHP机组的出力水平越高运行方式更稳定。原因是CHP被视为园区内部的可靠电源不仅发电还能供热等于同时为电系统和热系统提供支撑。优化结果倾向于利用这种“电热双保险”来对冲新能源出力波动。第二电储能在不确定性增大时不再只盯着峰谷价差套利而是更多地保留“备用余量”。典型表现是确定性模型里储能往往在谷电时段充满、峰电时段放空表现出很强的套利特征而场景模型给储能设置了一个最低SOC门槛例如不低于20%或更高以保证在出现低光伏场景时仍有足够的上调备用来平衡系统。储能从“纯套利工具”变成了“套利加备用”的双重角色这是不确定性计入模型后最本质的变化。5.3 从结果中得到的运行启示这对实际运行有三点启示。一是调度计划必须预留足够的备用容量不能把新能源预测值当作实际值来排计划。二是热储能的价值在不确定性环境下被放大因为它既能参与热平衡调节又能间接为CHP“解耦”提供缓冲让机组在电供给富余时段仍然保持较高的可用性。三是备用容量的分布不能只盯着电网侧分布式发电和储能同样应该承担一部分备用责任否则所有风险都压在电网主变上容易造成购电费用和功率越限同时飙升。6. 写码和调参过程中的踩坑记录最后分享几个我在这个项目里真实踩过的坑。这些坑的共性是模型在纸面上完全正确但进了代码就出问题而且报错信息往往不会直接指向病因。6.1 场景数贪多导致求解时间爆炸最开始我把场景数设成了200个想着场景越密结果越精确结果模型求解时间直接飙到几十分钟以上而且求解器内存占用接近失控。后面做了场景削减实验发现把200个场景削减到30个成本计算结果只偏差了大约1%但求解时间从几十分钟降到了两分钟以内。这个性价比非常惊人。我的建议是场景数用“阶梯式增加”的方式试探从20个开始逐步增加到50个观察目标函数值的变化当增加场景数量带来的成本变化小于0.5%时就没有必要继续增加了。6.2 储能SOC约束写错引发的“永动机”这是最典型的一个模型错误。我在储能约束里曾把充电效率和放电效率的作用位置搞反导致约束在数学上允许系统在同一时段同时充电和放电通过来回充放凭空产生能量。这类错误在结果图上非常隐蔽因为整体功率平衡仍然成立但仔细观察单个设备的功率曲线会发现储能充放电经常同时大于零。排查方法是写一个后处理校验脚本对所有时段检查储能SOC变化是否严格等于“充电量乘效率减放电量除效率”一旦发现偏离就立即定位到约束公式。6.3 量纲不一致的灾难有一版代码里负荷数据来自现场采集系统单位是kW电价数据单位是元/MWh燃气热值单位是kJ/m^3我在目标函数里直接相加导致求解结果完全不可解释。后来我强制要求所有输入参数进入模型之前统一转换为标幺值或统一量纲——我习惯用MW、MWh、元/MWh作为基准单位。在load_data.m里做一步单位换算并打印日志每类参数换算完都复核一遍峰值基本能拦住这类问题。6.4 “可弃电”约束的建模选择新能源出力“可弃”这件事在数学上处理起来比想象中麻烦。很多初学者会写一个等式并网功率 弃电量 预测值。这本身没错但弃电量是连续变量且非负有时候还需要走0-1逻辑来避免“预测为0时弃电量为正”的荒谬结果。我的做法是把弃电量定义为连续非负变量并允许其等于预测值和并网功率的差值再让惩罚系数足够大。只要目标函数里惩罚够高优化器不会故意弃电不需要额外引入整数变量这样既保证了模型简洁也没有陷入过量0-1变量的泥潭。6.5 求解器不可行时的标准排查流程当求解器返回不可行时我的标准流程是先把所有场景去掉只留下确定性版本看约束是否本身矛盾。如果确定性版本可解就逐个场景加回来用二分法定位是哪个场景引起的冲突。还有一种高频原因是储能的周期回位约束与末端时段出力限制冲突比如最后两个小时热负荷很低但热储能必须回到初值此时很可能出现热功率不足或过剩的不可行。解决方法是给周期回位约束加一个松弛变量并让目标函数对松弛量做高额惩罚成本比反复手工调参数高得多。提示模型不可行时不要急着改约束的数值范围先检查是不是有两个约束在边界处交叉锁死。最常见的位置就是“储能SOC下限”与“末端SOC回位”的组合。最后再分享两点操作体会第一个体会是项目启动时一定先跑通一个确定性模型把所有设备的物理约束理顺之后再引入不确定性建模。直接一步到位把不确定性和多设备耦合同时上马一旦结果异常根本分不清是场景生成的问题还是约束写错的问题。分阶段调试每一阶段都保存一份可以运行的版本这样回退和定位都很方便。第二个体会是代码里把数据和模型严格分离会让后续的实验效率翻倍。我在scenario_gen.m里把不确定性水平做成一个参数变量然后在main.m里用循环遍历5%、10%、15%三档误差水平所有结果自动保存成结构体并绘图。换一个负荷剖面、调整设备容量全程不需要碰任何模型约束代码。这听起来像是老生常谈但实际项目里能一开始就按这个习惯做的人真的不多。这个项目后续还有很多可以扩展的方向比如引入碳交易成本、加入多时间尺度协调、考虑设备随机故障等。但就目前这套“计及新能源出力不确定性的电气设备综合能源系统协同优化”Matlab实现而言它已经足够支撑园区级能量管理的实际需求也能够作为进一步研究的可靠基线。
