做电热综合能源系统调度研究的朋友应该都有过这种体验模型框架搭起来很快真正让人熬夜的往往是“为什么结果不合理”和“为什么求解器不收敛”。这段时间我刚好在Matlab环境里完整复现了一个考虑可再生能源消纳的电热综合能源系统日前经济调度模型从数学建模、求解器选型到算例验证都走了一遍踩了不少坑也沉淀了一些经验写出来供大家参考。这个模型的核心价值很明确把电力系统调度和热力系统调度放进同一个优化框架里让电锅炉、储热罐、热泵这些设备参与削峰填谷从而提升风电、光伏等可再生能源的消纳水平。无论你是刚接触综合能源系统方向的研究生还是已经在做园区级多能互补工程的工程师这篇内容都会对你有用。我会围绕模型为什么这么建、Matlab代码怎么组织、求解器怎么选、算例结果怎么分析这几个层次展开最后把我在调试过程中遇到的高频问题一并梳理出来。1. 为什么必须把热网放进日前调度里1.1 “以热定电”是弃风的第一推手北方冬季供暖期热电联产机组CHP的发电出力并不自由。为了保证供热温度机组必须保持一定的抽汽量而抽汽量直接约束了发电出力区间——这就是行业内常说的“以热定电”。简单说当热负荷高的时候CHP机组的电出力下限会被抬高于是系统中其他电源的可用空间就被压缩了尤其是优先性较差但边际成本极低的风电第一个被牺牲掉。我最早做这个项目时第一版模型只考虑电力系统的平衡结果非常典型夜间热负荷不变但电网负荷很低风电出力却达到峰值。因为没有热力环节的调节手段调度结果直接砍掉大量风电弃风率一度超过30%。这时候我才意识到问题不在风电预测不准而在电热两个系统之间存在“权力的不平衡”——热负荷在早上六点提交需求电系统只能在剩下空间里安排风电这本质上是一个系统耦合结构导致的资源错配。所以把热网纳入优化范围目的不是把两个系统强行拼接而是让热负荷也具备“可调节性”为风电让出空间。常见的手段包括电锅炉、吸收式热泵、储热罐等它们共同的作用是拆掉CHP机组的“电热刚性耦合”。1.2 热网本身的“储能属性”被大多数人忽略了很多初学者有个误区觉得储热罐才是唯一的热储能设备。事实上热力管网本身就是一个自带惯性的大规模储能系统。供热管道的温度变化是缓变的热水在管网里循环一圈需要几十分钟甚至数小时这意味着热负荷的供应不必严格实时匹配——只要保证在几个小时的时间窗口内有足够的蓄热量用户端的室温就不会明显波动。日前调度的时间尺度是小时级而且调度视野是24小时这就给了热网温度平抑策略足够的发挥空间。在建模时可以做两种选择一种是忽略管道动态只把热负荷当成普通节点负荷另一种是把管网模型简化为一阶惯性环节或储热等效节点。后者虽然多了一些约束但模型对“风电高发时段提前蓄热”这种操作的表达能力完全不同。我个人的建议是第一版先别追求管网动态精度先用储热罐热泵的配置验证调度逻辑跑通了再往里面加管网细节。1.3 目标函数不是成本最小这么简单经济调度的目标函数看起来非常直白——“总运行成本最小”但展开之后每一个成本项都有自己的脾气。常规机组有煤耗成本、启停成本CHP有电热联合生产成本电锅炉有耗电成本这些都好办。真正需要认真权衡的是“弃风惩罚”怎么定价。如果惩罚系数设得过高模型会为了让风电多发而不惜启动昂贵的电锅炉从全系统经济运行角度反而不划算如果惩罚系数设得过低模型会倾向于弃风。业界通常的做法是把弃风惩罚设置为常规机组发电成本的1.2到1.5倍让它在“多烧煤”和“弃风电”之间形成合理的经济折中。还有个可选加分项是碳交易成本。在双碳背景下很多省份已经把碳配额纳入火电成本这时候目标函数里要加上碳排放量的线性惩罚项。不过要注意碳价参数本身就是不确定的强耦合进调度模型后会给敏感性分析带来额外工作量除非论文或工程确有需要第一版不建议加入。2. 日前经济调度模型的数学表达方式2.1 决策变量和系统边界一定要先画清楚建模之前我习惯先画一张系统能量流图电源侧接风电场、光伏电站、常规火电机组、CHP机组热源侧接CHP余热、电锅炉、储热罐负荷侧是电负荷和热负荷系统边界处是联络线可以与大电网交换功率。这张图决定了所有变量的定义范围画错了后面会很痛苦。具体的决策变量包括每台CHP机组的电出力和热出力、火电机组出力、风电场实际上网功率、光伏实际出力、电锅炉耗电功率和产热功率、储热罐充放热功率和储热 SOC、联络线交换功率以及各机组的启停状态二进制变量。这些都是典型的连续变量加整数变量混合体所以模型天然是一个混合整数优化问题。值得额外注意的是可再能源出力变量的处理方式。很多文章喜欢用“预测出力等于实际上网功率”的等式约束这是很粗糙的。更合理的做法是让实际上网功率成为一个决策变量只约束它不小于0且不大于预测值至于弃风弃光量则通过目标函数的惩罚项来优化决定。这就是“软约束惩罚”的设计思路比硬性约束科学得多。2.2 约束条件的全家福与分类约束条件大致可以分为四层。第一层是平衡约束包括电功率平衡约束、热功率平衡约束核心思想是每个时刻“源侧出力储能释放外购电”必须等于“负荷储能充入”。第二层是设备出力约束包括机组上下限、爬坡速率限制、电锅炉容量、储热罐储热容量和充放热速率限制。第三层是网络约束包括支路潮流限制和节点电压上下限视算例规模和精度需求可以用直流潮流或交流潮流建模。第四层是启停逻辑约束涉及最小开关机时间、启停状态与出力范围联动等。我在建模过程中用到的约束分类可以看下面这个表格约束类型典型表达式建模要点电功率平衡∑P_G ∑P_W P_buy P_load P_EB P_HP外购电变量方向要统一符号别弄反热功率平衡∑H_CHP H_EB H_dis H_load H_ch储热放热为正储热为负机组出力约束P_min·u ≤ P ≤ P_max·u二进制变量配合大M实现逻辑联动爬坡约束-RD ≤ P(t1) - P(t) ≤ RU要考虑启停时的出力量是否跳过爬坡限制储热约束SOC_min ≤ SOC(t) ≤ SOC_maxSOC初值要和24小时末值联合约束形成日循环风电消纳约束0 ≤ P_W(t) ≤ P_W_pre(t)避免写成等式导致强制弃风或无谓惩罚这里面最容易出错的就是储热罐SOCState of Charge储热状态的日循环约束。日前调度一般是周期性调度如果SOC只约束下限上限而不约束一天结束时的状态那么优化结果会倾向于把热量在最后时刻全部放光以便“低成本”满足负荷结果完全失真。解决方法是加一个等式约束SOC_T SOC_0或者至少约束SOC_T ≥ SOC_min_return保证第二天调度有足够连续工况。2.3 非线性项的线性化处理电热联合调度模型里有两个经典的非线性来源一个是火电机组的煤耗曲线通常是凸二次函数直接进优化模型就成了MIQP混合整数二次规划虽然不是不能求解但求解效率要让位于MILP另一个是CHP机组的可行域工程上常用一个梯形或五边形区域来刻画电热出力耦合关系边界是分段线性的但引入二进制变量后同样会变成混合整数结构。我在实现中选择的做法是把煤耗成本曲线按出力区间做分段线性近似PWL。简单说就是把发电出力分成若干段每段对应一个固定斜率整体近似出一条凸分段线性费用曲线。分段数取3到5段时精度已经能控制在1%以内但求解速度会比完全二次模型快数倍。这个线性化操作在Matlab里用矩阵形式写非常简单把每个区间段的出力分解成多个连续变量并设置每段上限再在目标函数里用各段斜率乘上对应段出力即可。核心要义是确保分解后各段的取值范围形成组合约束防止模型“跳过前面的区间直接用后面的高价段”。这里必须强调一个经验分段线性化引入的额外变量和约束并不多但收益很大。我在测试中同一个38节点算例原MIQP在CPLEX里要跑40多秒分段线性化后的MILP只需要6秒代价仅仅是成本误差0.3%在工程上完全可以接受。3. Matlab代码实现中的关键细节3.1 算例结构和数据组织方式项目选用中国电科院一个经典38节点电热耦合算例做骨架。该算例包含常规火电节点、风电场节点、CHP节点和热负荷节点拓扑结构适中非常适合学术验证。数据组织是整个编码阶段的基石我强烈建议把所有原始参数统一放到一个结构体或表格对象里不要在求解脚本里散落硬编码数字。我的数据组织方式是定义一个名为sys的结构体内部字段包含gen发电机参数表、chp热电联产参数表、wind风电预测时序、load电负荷时序、heat_load热负荷时序、grid网络拓扑邻接矩阵等。每个字段的列顺序在注释里写清楚比如gen的第1列到第6列分别是机组编号、最大出力、最小出力、爬坡上限、爬坡下限、煤耗系数。这样做的好处是后边求解器生成约束时可以直接按列索引构建矩阵不容易出现数据对不上号的情况。另外一点容易被忽视的是时间粒度的选择。欧美文献习惯用1小时为调度间隔国内许多园区项目用15分钟一个断面全天96个点。时间粒度选得越小系统动态表达越精细但变量规模翻倍求解时间非线性增长。我的建议是第一版先从24小时粒度跑通逻辑确认结果合理后再切到96点精细优化两个版本之间只改时间轴数组和负荷曲线的插值方式其他代码结构完全一样。3.2 求解器选型CPLEX、BONMIN还是原对偶内点法Matlab环境下求解优化问题有几种路线各有利弊。如果模型最终是一个MILP商业求解器CPLEX效果最稳定YALMIP或MATLAB自带的intlinprog接口都能调用。这里有个容易混淆的点许多文献写的是“用BONMIN求解”那是针对MINLP混合整数非线性规划的求解器。如果你的模型里保留CHP可行域非线性项或交流潮流约束用BONMIN这类外点法求解器是合理的但收敛速度和初值依赖性都是问题。我个人的实操体会是第一版务必把模型改造成MILP用CPLEX跑先确认调度结果在经济学上是合理的如果需要计及真实交流潮流或热网水力工况的非线性约束才考虑引入原对偶内点法或迭代式潮流求解后者我会在下一节单独说。至于有人用智能优化算法比如遗传算法、粒子群求解我理解某些论文需要对比算法性能但工程上不建议作为首选因为可重复性和全局最优性都无法保证。使用YALMIP的建模过程其实可以写得非常简洁。定义一个包含24小时所有时段优化的Constraints单元数组循环叠加平衡约束和设备约束目标函数用sum求和最后调用optimize(Constraints, Objective, options)即可。如果不用YALMIP直接用intlinprog就要把所有时段变量一次性都堆进一个大向量手动构建矩阵和不等式系数代码会多出几百行。这里不是看不起手写矩阵而是明确告诉大家在Matlab里做这类研究YALMIP是生产力工具值得优先掌握。3.3 雅可比迭代与潮流约束的嵌入逻辑当需要在日前经济调度里加入电网潮流约束时一个常见的做法是用直流潮流近似——精度足够且保持线性。但直流潮流丢掉电压和无功信息如果论文里要算概率潮流指标就得回到交流潮流框架。折中方案是把潮流计算放进一个外层迭代循环内层用原对偶内点法求解调度主问题得到各节点注入功率后外层用牛顿-拉夫逊法计算交流潮流并提取雅可比矩阵和灵敏度信息把线路过载或电压越限的约束反馈回主问题通过增广拉格朗日项或者割平面迭代修正。这个过程很像工程里的“分解协调法”虽然数学上不是全程严格最优但在实际应用中收敛性不错计算时间也可控。这种迭代框架在Matlab里用自带优化工具箱就能搭起来并不一定需要CPLEX参与。唯一需要注意的是潮流方程中节点的PV/PQ转换逻辑一定要写对否则雅可比矩阵奇异时会直接报错排查就非常费劲。如果嫌麻烦也可以直接把潮流约束简化为功率传输分布因子PTDF矩阵用直流潮流的灵敏度公式把传输容量约束写成线性不等式。这个方法在38节点规模下完全够用。3.4 代码结构建议五段式布局一份完整的日前调度Matlab程序建议按五个模块组织数据加载模块load_data.m、模型构建模块build_model.m、求解模块solve_model.m、结果后处理模块postprocess.m、绘图输出模块plot_result.m。这套分层设计的核心价值在于调试友好当结果出现反直觉的现象比如弃风率不降反升我能很快定位是数据问题、约束写错还是求解参数不合适而不是在一整片脚本里从头到尾翻。说一个真实的教训我曾经在postprocess阶段复核电量平衡发现每个断面上总有0.5 MWh左右的电量缺口。排查了一整天最后发现是电锅炉的耗电功率变量在模型里被定义成了“热功率输出”而非“电输入”两者之间少乘了一个能效比系数。这种单位换算错误从语法层面完全抓不到只能靠结果复核校验逻辑来发现。所以强烈建议在每个模块间增加中间变量打印或者写入Excel人眼检查几个关键断面。4. 可再生能源的不确定性到底怎么纳入模型4.1 确定性调度的局限性刚才讲的全过程本质上是一个确定性优化风电预测曲线是给定的模型在“完美预测”的前提下做决策。但实际上风电预测误差是非常显著的尤其在日内尺度上。确定性调度结果虽然看着成本低、消纳率高但一旦实际风力和预测值差出30%备用容量不足系统就得实时切负荷或者紧急调整机组实际运行成本会远高于日前计算结果。我见过不少文献在“结论”里直接写“本模型可有效提高可再生能源消纳率”却没有交代预测误差带来的风险补偿成本。这类结论其实是不严谨的。所以如果想证明你搭建的模型有应用价值建议在确定性模型跑通之后至少做一次蒙特卡洛仿真测试把风电预测误差按正态分布或实测误差分布抽样模拟多个场景看看日前调度方案在滚动修正后的实际运行成本分布情况。4.2 随机优化与鲁棒优化两种最常用的路线考虑不确定性的建模路线业界主流有两套。一套是随机规划Stochastic Programming生成若干风电出力场景每场景赋一个概率目标函数变成所有场景运行成本的期望值最小化同时约束分别在各场景下满足。这个方法的好处是信息利用充分、经济性较好坏处是需要提前做场景削减否则模型规模会爆炸。常见的场景削减手段有K-medoids聚类、同步回代削减法Backward Reduction等可以把500个场景缩到20个而误差保持在可接受范围。另一套是鲁棒优化Robust Optimization假设风电出力在一个不确定区间内波动目标是找到即使在最恶劣出力组合下也可行的调度方案必要时代价为让渡部分经济性。鲁棒优化的关键是不确定集合的设计一般会用盒式集合加预算约束避免所有节点同时达到最恶劣出力这种保守场景。从计算角度说鲁棒对偶后模型仍可保持为MILP很适合在Matlab里用CPLEX求解。我的建议是如果你的研究重点在消纳机制和系统耦合先做随机优化如果你的研究重点在极端天气下系统的安全韧性那鲁棒优化更合适。两者并不冲突也可以做成两阶段模型日前阶段做鲁棒决策日内阶段利用滚动预测做随机调整这种模型写起来复杂但能同时覆盖经济性与保守性。4.3 概率潮流与机会约束的实操接口在计及风电不确定性时还有一个工程实现细节值得单独说就是概率潮流。传统潮流计算给定一组注入功率输出一组确定的潮流分布概率潮流则是给注入功率一个概率分布输出节点电压和线路潮流的概率分布比如越限概率为5%时对应的最低电压。这一指标特别适合用来生成机会约束比如约束“线路传输功率越限的概率不能超过5%”通过概率潮流的灵敏度分析可以得到线性化的近似不等式。我在Matlab里的实现方式是分三明治结构外层抽样生成风电出力场景中间层每个场景跑一次直流潮流或交流潮流统计线路流过功率的分布内层再把分位数约束线性化回到主问题。这个结构的计算量主要在外层场景数量上一般500个场景就足够稳定。需要注意的是当变量维数较高时直接保存所有场景的潮流结果矩阵可能会让内存吃紧建议边抽样边累积统计量不要一次性存完整场。4.4 改进智能算法适合在什么位置我身边不少朋友拿到模型第一反应是“用改进灰狼优化算法IGWO跑”这个思路在论文里做算法性能对比没有问题但如果模型本身是一个大规模MILP用智能算法求解实际上是拿自己的短处和商用求解器的长处硬碰硬。我测试对比过CPLEX和IGWO在同一个38节点确定性模型上的表现CPLEX得到全局最优解耗时6秒IGWO在2000次迭代后得到的解仍比CPLEX成本高约2.3%并且每次运行结果有波动。所以我的建议是当模型保持MILP结构时优先用CPLEX当模型引入高度非线性的热网水力学或交流潮流约束导致无法直接交给商业求解器时再考虑IGWO等启发式算法作为替代求解手段而不是把它当“默认解法”。5. 算例结果分析、敏感性测试与常见故障排查5.1 典型结果怎么看消纳率提升是表象日循环是本质我在上述38节点算例上设定了四种运行场景做对比纯电调度不考虑热网、电热联合且配置电锅炉、电热联合且配置储热罐、电热联合同时配置电锅炉和储热罐。其他参数保持一致风电装机容量占比设为总负荷峰值的25%。结果对比如下运行场景弃风率%总运行成本/万元储热罐日循环次数纯电调度28.7132.6无电热联合电锅炉14.2127.9无电热联合储热罐9.6125.40.9电热联合电锅炉储热罐4.8122.11.2从结果里能读出几个有价值的信息。第一单纯加电锅炉可以降弃风但电锅炉本身耗电在高电价时段运行反而推高成本所以经济性提升不显著。第二储热罐的加入让热电联产机组可以在夜间多储热、白天少产热实际上是把CHP的发电出力从夜间高位移到白天高位更好匹配了电网负荷曲线成本下降更明显。第三两个手段组合起来成本降幅并非简单叠加而是有协同效应因为储热罐给了电锅炉更灵活的运行区间电锅炉不必在风电出力高峰一股脑全开。5.2 敏感性分析到底要测哪些变量仿真模型跑通之后除了展示“分时电价下系统经济性最优”这种结果我建议做两个方向的敏感性分析。第一个是风电渗透率从10%逐步增加到40%观察弃风率拐点在哪里。很多算例在渗透率20%以下时弃风率增长平缓一旦超过某个阈值弃风率会陡增。这时候你需要调整储热容量或增加电锅炉容量观察扩容对拐点的推动作用。这个结果对系统规划非常有参考价值能为储能配比提供量化依据。第二个是热负荷水平变化的影响。北方供暖期的热负荷并非恒定初冬和深冬差得很多。可以在热负荷峰值上下浮动20%做场景扫描看调度成本和弃风率的谱线趋势。我自己实际跑下来发现热负荷越高风电消纳率下降越明显——因为高负荷下CHP电出力下限抬升得更狠这从物理直觉上完全说得通但把这条曲线表绘制出来放在论文里审稿人会明显觉得你的分析有深度。5.3 求解不收敛和不可行问题排查手册下面这部分是我最想写给所有调试者的内容。Matlab环境下做优化模型遇到的报错就那几类每一类的排查路径都很固定。第一类“Infeasible problem”不可行问题。出现这个提示首先查的是约束矛盾而不是求解器参数。我常用的二分法定位思路是这样的先把所有整数变量固定成某组可行值然后依次屏蔽约束集合找到让模型恢复可行的那一组约束再仔细审查这组约束里的边界条件和单位。最常见的坑是热负荷和电负荷的量纲不统一比如一个用了MW一个用了GJ/h两者在热平衡方程里直接相加数值差了近三倍。第二类求解时间爆炸。模型不是不可行就是半天不出结果。这类问题的根源多数是决策变量没有缩放到合理量纲。Chp出力是兆瓦量级储热罐容量是吉焦量级如果变量直接混在一起目标函数里某一项的数值远大于另一项求解器会把大量时间花在克服数值病态上。我的处理办法是统一除以系统基准值比如100 MVA或某个典型负荷值把变量归一化到0.1到10这个舒适区间求解速度往往提升一个数量级。第三类结果震荡或初值敏感。这个问题多出现在用内点法或启发式算法时。内点法对初值很敏感一个不合理的初始点可能导致迭代发散。我的经验是先让所有机组从中间出力比如50%额定起步储热罐SOC从50%起步求解稳定性会好很多。如果用IGWO这类启发式算法结果震荡多是收敛条件设置得太宽松或者种群规模太小需要把最大迭代次数拉长并重复跑多次取最优。第四类雅可比矩阵奇异。这个报错十有八九是节点类型设置错误比如把某个本来应该做PV节点的发电机节点写成了PQ节点导致该节点既没有给定电压也没有给定无功出力潮流方程缺了一个自由度。排查方法是打印节点类型和节点注入功率逐行核对是否有节点未被任何设备连接以及平衡节点是否设置正确。5.4 后处理阶段容易被低估的价值很多人在算例结果出来后就急着截图收工但后处理阶段如果能多做一步“日运行Gantt图”或“电热负荷平衡堆叠图”对理解调度逻辑有很大帮助。我在Matlab里用area函数画96点电功率平衡堆叠图底部是常规机组出力上面叠CHP电出力再叠风电出力最顶部是电锅炉耗电和联络线交换功率。这样一眼就能看出风电机组在半夜是被哪些设备“让路”的也能看出储热罐充放热时段和电价峰谷的关系。我还常用一个交叉验证脚本在求解结束后重新把所有平衡约束计算一遍检查最大不平衡量是否在容差范围内。这个脚本看起来“多此一举”但在我调试过程中至少帮我抓出过两次隐蔽的单位错误和数据错位比事后看曲线的效率高得多。做仿真研究的人和做工程的人有个共同点永远不要完全信任第一次运行的结果数据校验环节绝不能省。6. 我的几点体会与后续扩展建议做完整套模型最大的体会是综合能源系统调度的难点不在于某一类设备的运行约束有多复杂而在于多类设备之间的时间尺度、量纲和可调节性差异。电是秒级传输热能是小时级惯性风电是间歇随机热负荷又有很强的刚性拼劲。把这四个特性拧在一起还能用一套优化模型描述本身就需要反复平衡建模精度和求解可行性。对于刚入手的读者我的建议是先不要碰交流潮流也不要急着加不确定性。老老实实先把确定性MILP跑通电力平衡、热力平衡、储热日循环、风电消纳约束四大块都正确了再逐步扩展。这一步如果走了捷径后面所有分析都会建立在一个不牢固的基础上。另外我个人建议最好给这个模型加一个简单的“数据可视化Dashboard”哪怕只是把所有设备和负荷曲线画在同一个面板上。模型规模越大越需要一眼看清系统状态的工具。我在调试时还发现一个非常实用的小技巧每个求解断点都打印一次目标函数各成本分量的组成比如煤耗成本、购电成本、弃风惩罚各占多少。如果发现弃风惩罚在目标函数里占比异常高那就不是消纳能力问题而是惩罚系数设置不合理需要回到参数标定环节。目前这套模型还有两个方向想继续扩展一是把过去滚动调度机制加进去模拟日前计划和日内调整的双尺度衔接二是引入碳捕集设备的电碳耦合特性把碳流和能流一起优化。两个方向都还在搭模块阶段等代码调通之后我再把实际效果和踩坑经历整理出来继续分享。
