1. 项目概述为什么要做这个综合能源调度模型先聊聊这个题目本身。把光热电站、有机朗肯循环ORC、P2G电转气放在同一个调度框架里做联合优化是我在实际项目里接触过的典型场景——它本质上属于综合能源系统Integrated Energy System, IES的优化调度问题研究的对象不是单一的电或热而是电、热、气多种能量形式在时间尺度上的协同与互济。这两年做能源系统优化很多课题方向都往“多能互补”上靠原因是单靠一种技术太容易被天气和负荷牵着走光伏晚上没有、风电大风的白天可能反而不缺、传统火电调峰潜力又有限。而如果在一个园区或局域电网里把带储热的光热电站、能把低品位热能转化成电的ORC机组、以及能把不确定性电力转成氢/天然气的P2G设备衔接起来那情况就完全不一样了——光热电站可以利用储热罐把中午的太阳能挪到晚高峰释放ORC可以回收工业生产过程中的废热并二次发电P2G则可以在风光出力大的低谷时段把多余电能转化成气态燃料存储或外送。三者结合实际上是在不同时间尺度和不同能级品位之间搭了一座桥。我最初接触这类模型是在搞某区域能源互联网规划时当时Matlab代码就是最常用的落地工具。为什么用Matlab而不是Python原因很现实优势比较大的求解器生态在Matlab里都有成熟接口比如Yalmip作为建模层底层接Gurobi或Cplex对混合整数线性规划MILP问题的求解非常友好。在电力行业里Matlab做数据分析、曲线绘制和优化求解的案例积累也更多身边很多同事、高校研究者的代码库里都有大量现成的Matlab函数可以直接复用。这篇博客我想用通俗点的口吻把整套思路拆开来讲含光热电站、有机朗肯循环、P2G的综合能源优化调度模型怎么搭、目标函数怎么写、约束怎么列、Matlab代码怎么一步步实现、运行中又会遇到哪些坑。内容适合正在做综合能源系统课题的研究生、刚入行做园区能源规划的工程师也适合一切对“多能互补调度”感兴趣但对求解落地还比较陌生的人。目的是让你看完之后不需要再去翻一堆零散教程也能在Matlab里把这套东西跑起来。2. 系统架构与核心设备建模范式2.1 三种设备的角色辨析光热电站Concentrating Solar Power, CSP、ORC与P2G三种设备在系统里的“段位”是完全不一样的只有先理解它们的互补逻辑下一步建模才不会写成一锅粥。光热电站的本质是通过聚光集热装置把太阳能转化成高温热能再由热力循环发电整套系统最大的特点是带储热。实际工程中最常见的是双罐熔盐储热白天光照充足时多余的热量用来加热低温罐里的熔盐并存入高温罐晚上光照不足或负荷高峰时再用高温熔盐去加热蒸汽维持汽轮发电机组稳定出力。从调度视角来看光热电站几乎可以等价成一个“热储能可控发电机组”的复合体——它的电出力上限是由储热罐当前可放热量和汽轮机额定容量共同决定的这个特性让它在调度模型里天然具备灵活性。有机朗肯循环Organic Rankine CycleORC解决的问题则是另一个层面的——把中低温热源通常在80℃到300℃之间转化成电能。常规水蒸气朗肯循环在这种温度区间里经济性太差而有机物工质如R245fa、R134a等沸点低、气化潜热小能更适合低温余热回收。实际场景里ORC的输入端往往接的是工厂废热、地热、生物质锅炉烟气或者光热电站低温旁路发电规模从几千瓦到几十兆瓦都有。调度层面上ORC的出力不完全可控它与上游热源的连续性和温度直接相关在模型里如果把ORC简化成一个固定效率的线性转换单元会丢失很多实际运行特征。P2GPower to Gas电转气则属于长周期储能和跨网耦合的范畴电解水制氢氢气后续可以直接利用、可以掺入天然气管网也可以进一步甲烷化。P2G解决了新能源消纳里一个特别棘手的问题就是“电不好存”电网侧电池最多储几个小时但燃气储存的时间和规模大得多。在综合能源系统里建P2G等于给了调度人员一个“把电转成气送出去或存起来”的出口。从时间尺度上看光热储热解决的是小时级的移峰填谷ORC解决的是持续稳定热源的低品位回收P2G则把目光放到日间甚至周级的能量搬移。把它们放进同一套模型本质上是把这些不同时间常数的灵活性资源统一用数学语言描述出来然后在24小时或更长的调度周期里找一组“总运行成本最低”的出力决策。2.2 为什么把光热简化成“储热电站在调峰”建模仿真里经常要做一个选择模型复杂度与求解可行性之间的权衡。有人喜欢把光热电站内部的每个换热环节都详细建模甚至逐段模拟熔盐储罐的温度分层但这样得到的非线性模型在优化调度里几乎没法全局求解只能交给启发式算法慢慢迭代。我的做法和多数调度文献保持一致把光热电站做“黑箱-适度白箱”的聚合建模。所谓“黑箱”就是不过分关心场内的集热管路导热油流向和换热器面积只看宏观能量流。太阳能输入通过光场效率折算成热功率P_sf_t约束有集热场的热功率上限、熔盐储热罐的容量上下限与充放热速率上限、蒸汽发生系统的热输入与发电功率的转换关系。这样就得到一个包含“储热罐soc”和“储热罐充放热响应”的调度状态模型CSP能够在夜间也保持稳定出力的特性就体现出来了。再看ORC我不建议把它的热力学循环蒸发温度、冷凝温度、工质流量在调度里展开做。因为调度的时间颗粒度通常是1小时优化要的是宏观上的“在某时段从某热源取了多少热、产了多少电、热转电效率是多少”这样的数量关系。实际上ORC的发电效率与热源温度存在耦合如果热源来自光热低温侧或工业余热可以在不同时段给ORC赋不同的效率比如设定三档变工况效率查表这样既不引入非线性也能保留机组的部分变工况特征是一种工程上折中、但足够实用的处理办法。P2G也一样电解槽的功率-产氢率曲线本质上是非线性的但实际项目里通常在一定负荷区间内做线性假设给出额定功率和额定效率部分论文用分段线性或带整数变量的启停模型表示电解槽的最小技术出力限制。在MILP框架里这是完全可以实现的只需引入一个二进制变量表示启停状态。2.3 多能源耦合的瓶颈环节在哪里真正让这套系统复杂起来的不是单个设备建模而是耦合约束。电、热、气不是三条独立的管道它们通过设备互相转换光伏/风电在电力母线上汇集供给电负荷与P2GP2G产气汇入气负荷或储气罐光热的储热罐同时服务汽轮机发电和热负荷ORC又把系统中的某些中低温热流转化为额外的电流。因此调度模型里会同时出现如下交叉变量电母线平衡方程里等式右侧既有常规机组出力、光热电站电出力、ORC电出力、新能源出力左侧则有电负荷、P2G耗电热母线平衡里光热电站可调节热输出、储热罐充放热、ORC热端取热等形成热量分配关系气母线中P2G产气量、外部购气量、气负荷等构成平衡。这些耦合关系最终全部落地为优化模型里的等式约束或不等式约束。写代码的时候必须规划好变量命名空间和维度否则后面出错的概率会非常大。我见过很多同学在Matlab里建模型时过一个月回来看代码就看到变量顺序混乱这个在后面第4部分通过代码案例详细展开。3. 目标函数与约束体系设计从物理到数学3.1 两个层次的目标函数怎么取舍优化调度按决策目标不同通常分成两大类。类一以运行成本最小为目标适用于园区、局域电网等经济性敏感的运营场景目标函数由“购电费用 燃料成本 启停成本 弃风弃光惩罚 设备运维成本 购气费用”构成。类二以新能源消纳最大或碳排放最小为目标在政策驱动场景下用得比较多。实际工程项目里经济目标与环境目标往往通过加权系数合在一起。不过对于初学者我建议先做成本最小模型因为单位统一元结果容易量化检查等代码跑通了再加碳约束并不难。目标函数写出来后需要给每一项定义清楚系数的物理意义比如新能源弃电惩罚是不是应该高于其他发电成本如果惩罚设得太低优化结果里可能出现“为了让成本表好看而主动废掉本来可用的新能源出力”的情况这明显不符合调度初衷不过惩罚设得太高又可能扭曲正常购电/购气决策。这类系数的敏感性分析要多跑几次才能找到合理的区间。我常用的经验是单位弃电惩罚电价取系统最高购电电价的1.2~1.5倍左右可以保证新能源优先消纳又不至于让结果失真。3.2 设备约束要细化到什么程度设备约束是整个模型里最琐碎也最考验细心的地方。系统里每一个可调度设备都要回答四个问题出力上下限是多少爬坡速率限制是多少运行会不会有最小启停时间限制储能相关量的状态转移方程和容量约束怎么表达有一说一很多刚接触优化的同学容易犯一个错误就是只写上下限和平衡约束漏掉储能状态量相邻时段耦合的约束导致储热罐在凌晨1点放的热可能凭空“瞬移”到下午3点再用——这在数学里可能没问题物理上是假的。具体到CSP储热罐我认为最关键的约束是储热罐当前蓄热量的状态转移方程S_t S_{t-1} η_ch * P_ch_t - P_dis_t / η_dis - S_loss其中S_t表示t时段末储热罐内的热能量MWhη_ch和η_dis分别是充放热效率P_ch_t和P_dis_t是t时段的充放热功率S_loss是散热损失。注意这里的充放热效率和功率都是写在“热侧”的光热电站的电出力由“蒸汽发生系统从储热罐取热集热场直接供热”驱动汽轮机热-电转换关系也要相应地建模。ORC机组的约束要分清它的热源来源。在综合能源系统里通常存在多种可能的余热源我们要给ORC配置一个“热源取热上限”同时发电功率P_orc_t η_orc * Q_orc_in_t其中Q_orc_in_t为t时段ORC从热源取的热功率这样约束简单但足以体现ORC运行的可行性边界。还有需要注意的是ORC冷凝器排出的余热是低温热不太可能回收利用在建模时不需要画蛇添足把“冷凝废热”当成产出除非你的系统专门设计了梯级利用链条。P2G的模型相对直接一些输入是电功率输出是氢气流量或计入气网的天然气等效流量。有些模型还要加一个氢气储罐用储氢罐状态来平衡P2G产氢和下游用氢/售氢之间的时移。在这套系统中P2G更常见的定位是消纳晚间风电或午间光伏高峰把本来要被限制的电力变成气体储存起来——这也意味着它的运行特征与新能源出力曲线之间存在强烈的时段耦合特性只能通过优化求解来确定具体的运行时刻。3.3 新能源出力与负荷曲线怎么处理调度问题永远都是一个“给定外界条件求最优决策”的过程风电、光伏的预测出力以及电/热/气负荷曲线属于调度模型的外生输入不是优化变量。但外生输入的质量直接影响模型结果。工程实操中需要把一整年的新能源出力和负荷数据聚类成几个典型日场景比如春夏秋冬四季或者丰水期/枯水期。每个典型日取24时间断面作为输入向量。如果没有实测数据可以用Matlab里的现有工具生成模拟曲线光照强度按正弦曲线加云层波动风速用韦伯分布抽样再做时序平滑负荷数据可以参照典型用户曲线叠加随机扰动。原始数据输入前必须做归一化或者去量纲处理。在优化求解层面我注意到如果不把数量级差异控制在合理范围比如电负荷是几百兆瓦时储能罐容量只有几百千瓦时Cplex/Gurobi这种求解器很容易报数值问题。经验做法是乘以一个基准值将功率统一折算到标幺值域中比如以1MW为基准那么100MW负荷变量量纲变为100储热罐蓄热量可能变成35之类数值范围就不会差别太大。求解器的数值稳定性会好很多这一步对后期调参来说绝对是事半功倍的。4. Matlab代码实现与调度模型落地指南4.1 代码框架搭建与Yalmip建模流程Matlab里建立一个确定性优化模型当前最主流的方式是先装好Yalmip工具箱再加一个求解器。Yalmip并不真正执行求解它像一个伟大的翻译官把人类易于理解的建模语句翻译成求解器能读懂的数学规划标准格式。只要你的模型是线性、混合整数线性、二次约束甚至非线性Yalmip都有对应接口去对接合适的求解器。在这个场景里我推荐用Gurobi或Cplex因为白天遇到MILP问题时它们的内点法和分支切割法表现都很成熟。需要提醒的是Yalmip的安装并不复杂简单来说就是下载yalmip压缩包后把文件夹加到Matlab搜索路径里。真正要花时间的是求解器license的配置尤其很多新手同学在学校电脑上用Cplex经常报出Failed to get license或Error: Unable to solve because solver is not instlled这种时候通常不是模型写错了而是求解器接口没配对。应该逐个检查求解器是否安装、环境变量是否设置、Matlab里调用solvesdp或optimize时是否明确指定了参数。避免因许可证或者路径问题影响进度。Yalmip建模的典型代码结构我通常按下面几个区段来组织参数初始化区设置系统参数、读取负荷和新能源预测曲线决策变量定义区用sdpvar定义连续变量用binvar定义二进制变量约束条件声明区把设备模型中的物理公式一行一行写进去目标函数声明区按成本项累加求解与结果输出区调用optimize并绘图。最忌讳的是从头写到尾只有一个大脚本变量作用域混在一起出了错不好定位。我自己实践下来比较舒服的方式是每个功能模块单独写一个函数但必须注意函数里创建的Yalmip变量返回后仍然保持一致因为它本质上是Matlab对象引用关系比较特殊容易踩坑。4.2 光热储热、ORC和P2G的代码怎么组织下面展示几个核心设备的建模片段来说明关键约束在Yalmip里的表达方式。先看光热储热罐与电出力部分的基础逻辑%% 光热电站相关变量 P_sf sdpvar(1, T); % 集热场接收的太阳热功率, MW P_gross sdpvar(1, T); % 汽轮机发电功率, MW P_ch sdpvar(1, T); % 储热罐充热功率, MW P_dis sdpvar(1, T); % 储热罐放热功率, MW S_soc sdpvar(1, T1); % 储热罐蓄热量MWh取时段初/末状态 P_res sdpvar(1, T); % 实际集热场输入到蒸汽发生系统的热功率 % 集热场输入可用太阳热功率受DNI影响上限随时间变化 Constraints [Constraints, 0 P_sf P_dni]; % P_dni是t时段光场可收集的热功率上限 Constraints [Constraints, 0 P_res P_dis Q_sg_max]; % 蒸汽发生系统吸热上限 % 发电功率与吸热功率成正比例 Constraints [Constraints, P_gross eta_pb * (P_res P_dis)]; % 储热罐状态约束 Constraints [Constraints, S_soc(1) S_init]; for t 1:T Constraints [Constraints, S_soc(t1) S_soc(t) * (1 - loss_rate) eta_ch * P_ch(t) - P_dis(t) / eta_dis]; Constraints [Constraints, 0 P_ch(t) P_ch_max * u_ch(t)]; Constraints [Constraints, 0 P_dis(t) P_dis_max * u_dis(t)]; Constraints [Constraints, u_ch(t) u_dis(t) 1]; % 不能同时充放热 Constraints [Constraints, 0 S_soc(t1) S_max]; end这里引入了一个必要的二元变量u_ch和u_dis表示充/放热状态“不能同时充放热”是物理层面的必需约束。在建模时如果省略那么优化可能会倾向于在同一时段用高温熔盐去充热同时从低温罐取出熔盐放热造成无意义的热量循环却不产生电能导致目标函数值失真。再来一个ORC的简单模型假设热源来自系统中某一固定废热流分时温度变化导致效率变化%% ORC 有机朗肯循环约束 P_orc sdpvar(1, T); % ORC电出力 Q_orc_in sdpvar(1, T); % ORC从热源侧取热功率 eta_orc_t interp1(temp_profile, eff_table, temp_source_t, linear); % 由热源温度查表 Constraints [Constraints, Q_orc_in 0]; Constraints [Constraints, Q_orc_in Q_orc_in_max]; Constraints [Constraints, P_orc eta_orc_t .* Q_orc_in];如果ORC的热源来自光热电站的低温旁路那么取热功率也会出现在光热系统热平衡里需要把这项加到集热场/储热罐的 热量分配关系中。P2G模块的核心代码相对简洁不过当电解槽需要启停状态时还得加上最小出力约束这一点往往被不少研究者忽略%% P2G 电制氢模块 P_p2g sdpvar(1, T); u_p2g binvar(1, T); % P2G启停标志 F_h2 sdpvar(1, T); % 氢气流量kg/h或m3/h for t 1:T Constraints [Constraints, P_p2g_min * u_p2g(t) P_p2g(t) P_p2g_max * u_p2g(t)]; Constraints [Constraints, F_h2(t) eta_p2g * P_p2g(t)]; % 线性效率 end需要注意的是如果有储氢罐那么储氢罐的状态约束与储热罐类似会增加一组S_h2的状态变量和充放氢变量。如果模型只允许P2G产出的氢气直接进入氢负荷或天然气管网那么可以不引入储氢罐平衡方程会简单很多——具体选择取决于你研究的问题是否要含氢气长周期存储。4.3 连接能量母线系统级平衡约束怎么写设备建模完成后必须将所有设备的端口变量汇聚到电、热、气三条能量母线的平衡方程中才能让模型真正形成一个“系统”。在我的代码里这些平衡约束的形式通常非常直观%% 电力母线平衡 for t 1:T Constraints [Constraints, ... P_buy(t) P_wt(t) P_pv(t) P_gross(t) P_orc(t) ... P_load(t) P_p2g(t) P_eb(t)]; % P_eb是电锅炉或其他电负荷 end %% 热力母线平衡此处以区域热网所需热量为热负荷光热供汽轮机后剩余热量供应热负荷 Constraints [Constraints, P_heat_supply P_gb P_eb H_load]; %% 天然气/氢气母线平衡 Constraints [Constraints, F_h2 F_buy_gas G_load F_boiler_gas];这里需要特别强调一点电平衡里各项的正负号非常容易出错。买电、新能源出力、发电设备输出是在方程左侧的“源”项电负荷、P2G的耗电则应放在右侧即“荷”项。写反了很容易造成优化结果产生难以置信的“凭空发电”或“凭空消失的电”遇到这种问题可以先检查各项的正负符号。目标函数部分按成本最小形式配置Objective sum(P_buy .* price_buy * dt) ... % 购电成本 sum(F_buy_gas .* price_gas * dt) ... % 购气成本 sum(P_gross .* c_om_csp * dt) ... % 光热运维成本 sum(P_orc .* c_om_orc * dt) ... sum(P_p2g .* c_om_p2g * dt);dt是调度时段长度如果时间颗粒度为小时那么dt1功率乘以时间等于能量乘以电价才是费用。很多第一次做调度优化的同学会在这一步漏算时间颗粒度导致成本低了一倍或高了一倍最后得到的结果无论解释成什么都别扭。4.4 求解器配置与小规模算例结果解读代码写完直接调用optimize时建议显式指定求解器参数。例如用Gurobi就设置options sdpsettings(solver, gurobi, verbose, 2, showprogress, 1); result optimize(Constraints, Objective, options);若模型是纯线性的输出里的gap在零附近solvertime不超过几秒如果条件表达式里埋伏了不必要的非线性函数Yalmip会把它转成二阶锥或非线性模型此时Gurobi往往无法直接求解Yalmip会提示“No suitable solver”。遇到这种情况最简单的办法是检查使用的函数是否引入了非线性比如用乘法c * x没问题这两个是1*n向量点乘但你如果写成x * c且c也是变量就是二次项了。另外绝对值如果不想引入非光滑项可以用变量对偶法或引入辅助变量线性化表达。小算例做出来后通常可以画出三种图一是典型日全天各电源的电力出力堆叠图二是储热罐24小时蓄热状态与P2G运行功率的对比图三是在不同场景下成本项的柱状图。通过这些图就能够验证模型行为是否合理。如果发现储热罐只在午后充热、晚高峰放热那么说明储热的经济价值被正确捕捉了如果发现储热罐压根没被调度使用SOC恒定不变大概率是目标函数里某些成本项设置过大或储热容量本身过大/过小需要回到参数设置排查。5. 常见问题与排查技巧实录5.1 模型不可行或不收敛先查这几处在调度优化里最崩溃的报错就是Infeasible problem。由于约束太多到底哪个限制条件出了问题求解器往往不会直接指出来。我常用的排查思路是按“从简到繁”的顺序来定位第一步先把所有储能相关的约束注释掉看系统本身平不平衡。如果不行大概率是功率平衡方程或热平衡方程中某项不能同时满足单位或许有问题kW/MW混淆第二步只保留储热约束检查储热罐在初始和终末SOC相同的条件下是否能够在各个时段完成充放。这是验证调度模型正确性的关键一步。很多模型强制要求S_{T1} S_1以方便进行多日滚动调度但在给定初始状态的情况下容易导致无解建议先把循环调度约束放宽成不等式比如末时刻SOC要高于某阈值即可第三步再把所有约束放回来从目标函数中逐项删除来观察删除哪一项后模型可行通常能找到是哪类量卡的太紧。另外把所有设备启停的二元变量都设为常数1看连续松弛问题是否可行可以快速判断是否存在零/一整数的组合死角。比如你要求储热罐不能同时充放但一天中储热罐需要先在一定时段满放再满充而约束又把充放热上限限制在一个极小的值就可能找不到离散方案的可行解。5.2 Yalmip建模常见报错自查表我把平时见到的高频报错整理成一个表格大家可以直接对照排查。虽然不是绝对覆盖但至少覆盖了80%以上的入门问题常见现象可能原因解决建议No suitable solver模型包含非线性项但求解器不支持用yalmiptest检查已安装求解器用export导出模型查看类型Infeasible problem平衡约束无法满足且模型被过约束逐模块放行排查不可行约束来源用relax看对偶变量找出弱约束Index exceeds matrix dimensions读入的负荷曲线长度不足24点观察输入数据维度用size(load)检查求解时间异常长二元变量太多分支定界计算量爆炸削减整数变量的范围把某些离散状态处理为设计参数结果出现“充放电循环但不产生真实吞吐”储能约束松弛或漏约束确认加入了同时充放不可行约束及能量状态转移方程目标函数值是NaN初始值含NaN或除以零清洗所有输入数据检查效率参数是否为零5.3 数值调优与加速的个人心得解决完“能跑”的问题之后接下来就是“跑得快”和“算得稳”的环节了。这里分享几个最直接的加速技巧避免在约束循环里重复创建新的Yalmip变量。比如Constraints [Constraints, ...]在一个大循环里反复拼接其实效率不差但如果是多维大矩阵建议一次性把整个变量矩阵定义出来用repmat或矩阵运算一次性生成约束能大幅缩短建模时间变量的取值范围尽量收紧。给每个决策变量写上下限边界Gurobi/Cplex做预处理时会把冗余约束删除加速比较明显。比如储热罐SOC如果明确知道最高只有12MWh就不要定义到100MWh免得数值比例失调对电力平衡和热平衡这类核心约束尽量整行矩阵化写不做逐一循环。例如电平衡可以用sum(A,2)这种向量化操作直接形成等式约束Yalmip内部会把它识别为一个大稀疏矩阵约束求解时的预处理开销小得多。我实际跑过的模型规模大概是24时段、十余台设备、储能罐两个、带少量启停二元变量总变量数量在四五百个左右Gurobi解决时间大约是几百毫秒到几秒。如果你的模型跑了几分钟还没结果一般情况下不是设备数量多而是有地方产生的二次约束或整数变量较多。这种模棱两可的问题建议先冷静分析模型本身而不是盲目增大求解器时间限制。5.4 调度结果里那些“看起来奇怪”但真实合理的现象有时跑完结果会有一些输出让人觉得哪里不对劲。比如P2G在某个时段明明是电价低谷却没有开起来。原因可能是P2G的技术最小出力限制导致这一个小时内制出的氢气没有地方可去——要么是因为下游气负荷此时段很低要么储氢罐已满即使电价低强开P2G只会增加弃氢风险。这就提醒我们看到结果不能只凭直觉判断要结合存储状态的完整曲线来分析。还有一种情况是光热电站傍晚把储热基本放完但夜间1点的热负荷靠电锅炉顶着看起来好像不如“晚点放热”更划算。为什么优化结果是这样因为模型如果在傍晚低负荷时段不发电就需要从外部电网买电购电价格可能远高于储热发电的边际成本所以它宁可傍晚集中放热把节省的购电费放到最大化。对这类现象需要用“边际成本”的语言去解释以验证结果的合理性。6. 从案例到实践整套系统的可扩展方向与个人经验代码跑通的那一刻是很开心的但综合能源系统的优化调度从来不是从零开始完整复现一遍就结束的研究课题。在实际工程或论文延展里真正的价值在于把基础模型往不同方向扩展。对这个含光热、ORC和P2G的模型来说最容易也最有价值的扩展方向有这样几个多时间尺度滚动调度。上面阐述的是单典型日24小时的确定性优化在实际运行调度时需要做日内滚动修正。例如用第1~4小时实际数据校正新能源预测误差然后每4小时滚动优化一次重新求解未来24小时的决策。这个在Matlab里可以通过循环调用同一套模型实现对代码架构的复用很有要求如果你的约束和变量封装成了独立函数扩展滚动调度并不复杂。考虑不确定性。光伏和风功率预测不可能绝对准确常见的处理是采用场景法把风电场景集生成多个可能序列再把目标函数改成期望值最小或者采用鲁棒优化描述成一个min-max-min三层结构。在Matlab里随机规划通常和scenario配合鲁棒优化则用uncertain对象配合Yalmip。光热辐射本身的波动如天气变化引发DNI突变也值得纳入不确定性集合。我做过的模型里有一个例子把DNI风速负荷构建为盒式不确定集在Yalmip里实现两阶段鲁棒优化很考验建模和数值功底但是做完后对整个系统的弹性会有更深的理解。引入碳交易机制。如果目标是“双碳”相关的综合能源系统研究那么把碳排放配额和碳交易成本加入目标函数是自然的扩展。光热和P2G在这方面有天然优势光热电站比燃煤机组省大量燃料碳排放P2G制氢在部分场景下替代灰氢也有减排价值。碳价价格波动可能反而会改变调度结果中光热和ORC的优先序列这个方向特别适合做成敏感性分析。再有一个小而实用的小建议在代码开发过程中及早形成标准的绘图模板。我之前吃过亏——因为Data结果导来导去找不到哪个变量对应哪条曲线耽误了不少时间。建议把所有结果变量整理成一个结构体变量result.P_csp、result.S_soc、result.P_orc、result.F_h2画图也写成统一函数后期加任何设备只要在结构体里补一个字段就能自动出图作为校核。这样整个项目几个人协作也会比较顺畅。最后分享几次调模型积累下的体会综合能源优化调度的核心难点从来不在某一个公式的复杂程度而在能否把不同时间尺度、不同能量品位的物理过程统一到一个边界清晰、数值稳定的数学框架里。面对这类模型与其贪多求全把每个部件都写到极致不如先把关键设备的物理约束抓准、系统能量平衡关系写对、目标函数经济含义理清。保持“先跑通、再跑稳、后跑快”的推进节奏往往比一次想写出完美代码更有效。如果在这过程中遇到“反直觉”的调度结果不用急着怀疑模型出错了——先在设备约束、能量平衡和价格机制几个维度对照分析很多时候会得到比书本案例更有意思的发现。
1. 项目概述为什么要做这个综合能源调度模型先聊聊这个题目本身。把光热电站、有机朗肯循环ORC、P2G电转气放在同一个调度框架里做联合优化是我在实际项目里接触过的典型场景——它本质上属于综合能源系统Integrated Energy System, IES的优化调度问题研究的对象不是单一的电或热而是电、热、气多种能量形式在时间尺度上的协同与互济。这两年做能源系统优化很多课题方向都往“多能互补”上靠原因是单靠一种技术太容易被天气和负荷牵着走光伏晚上没有、风电大风的白天可能反而不缺、传统火电调峰潜力又有限。而如果在一个园区或局域电网里把带储热的光热电站、能把低品位热能转化成电的ORC机组、以及能把不确定性电力转成氢/天然气的P2G设备衔接起来那情况就完全不一样了——光热电站可以利用储热罐把中午的太阳能挪到晚高峰释放ORC可以回收工业生产过程中的废热并二次发电P2G则可以在风光出力大的低谷时段把多余电能转化成气态燃料存储或外送。三者结合实际上是在不同时间尺度和不同能级品位之间搭了一座桥。我最初接触这类模型是在搞某区域能源互联网规划时当时Matlab代码就是最常用的落地工具。为什么用Matlab而不是Python原因很现实优势比较大的求解器生态在Matlab里都有成熟接口比如Yalmip作为建模层底层接Gurobi或Cplex对混合整数线性规划MILP问题的求解非常友好。在电力行业里Matlab做数据分析、曲线绘制和优化求解的案例积累也更多身边很多同事、高校研究者的代码库里都有大量现成的Matlab函数可以直接复用。这篇博客我想用通俗点的口吻把整套思路拆开来讲含光热电站、有机朗肯循环、P2G的综合能源优化调度模型怎么搭、目标函数怎么写、约束怎么列、Matlab代码怎么一步步实现、运行中又会遇到哪些坑。内容适合正在做综合能源系统课题的研究生、刚入行做园区能源规划的工程师也适合一切对“多能互补调度”感兴趣但对求解落地还比较陌生的人。目的是让你看完之后不需要再去翻一堆零散教程也能在Matlab里把这套东西跑起来。2. 系统架构与核心设备建模范式2.1 三种设备的角色辨析光热电站Concentrating Solar Power, CSP、ORC与P2G三种设备在系统里的“段位”是完全不一样的只有先理解它们的互补逻辑下一步建模才不会写成一锅粥。光热电站的本质是通过聚光集热装置把太阳能转化成高温热能再由热力循环发电整套系统最大的特点是带储热。实际工程中最常见的是双罐熔盐储热白天光照充足时多余的热量用来加热低温罐里的熔盐并存入高温罐晚上光照不足或负荷高峰时再用高温熔盐去加热蒸汽维持汽轮发电机组稳定出力。从调度视角来看光热电站几乎可以等价成一个“热储能可控发电机组”的复合体——它的电出力上限是由储热罐当前可放热量和汽轮机额定容量共同决定的这个特性让它在调度模型里天然具备灵活性。有机朗肯循环Organic Rankine CycleORC解决的问题则是另一个层面的——把中低温热源通常在80℃到300℃之间转化成电能。常规水蒸气朗肯循环在这种温度区间里经济性太差而有机物工质如R245fa、R134a等沸点低、气化潜热小能更适合低温余热回收。实际场景里ORC的输入端往往接的是工厂废热、地热、生物质锅炉烟气或者光热电站低温旁路发电规模从几千瓦到几十兆瓦都有。调度层面上ORC的出力不完全可控它与上游热源的连续性和温度直接相关在模型里如果把ORC简化成一个固定效率的线性转换单元会丢失很多实际运行特征。P2GPower to Gas电转气则属于长周期储能和跨网耦合的范畴电解水制氢氢气后续可以直接利用、可以掺入天然气管网也可以进一步甲烷化。P2G解决了新能源消纳里一个特别棘手的问题就是“电不好存”电网侧电池最多储几个小时但燃气储存的时间和规模大得多。在综合能源系统里建P2G等于给了调度人员一个“把电转成气送出去或存起来”的出口。从时间尺度上看光热储热解决的是小时级的移峰填谷ORC解决的是持续稳定热源的低品位回收P2G则把目光放到日间甚至周级的能量搬移。把它们放进同一套模型本质上是把这些不同时间常数的灵活性资源统一用数学语言描述出来然后在24小时或更长的调度周期里找一组“总运行成本最低”的出力决策。2.2 为什么把光热简化成“储热电站在调峰”建模仿真里经常要做一个选择模型复杂度与求解可行性之间的权衡。有人喜欢把光热电站内部的每个换热环节都详细建模甚至逐段模拟熔盐储罐的温度分层但这样得到的非线性模型在优化调度里几乎没法全局求解只能交给启发式算法慢慢迭代。我的做法和多数调度文献保持一致把光热电站做“黑箱-适度白箱”的聚合建模。所谓“黑箱”就是不过分关心场内的集热管路导热油流向和换热器面积只看宏观能量流。太阳能输入通过光场效率折算成热功率P_sf_t约束有集热场的热功率上限、熔盐储热罐的容量上下限与充放热速率上限、蒸汽发生系统的热输入与发电功率的转换关系。这样就得到一个包含“储热罐soc”和“储热罐充放热响应”的调度状态模型CSP能够在夜间也保持稳定出力的特性就体现出来了。再看ORC我不建议把它的热力学循环蒸发温度、冷凝温度、工质流量在调度里展开做。因为调度的时间颗粒度通常是1小时优化要的是宏观上的“在某时段从某热源取了多少热、产了多少电、热转电效率是多少”这样的数量关系。实际上ORC的发电效率与热源温度存在耦合如果热源来自光热低温侧或工业余热可以在不同时段给ORC赋不同的效率比如设定三档变工况效率查表这样既不引入非线性也能保留机组的部分变工况特征是一种工程上折中、但足够实用的处理办法。P2G也一样电解槽的功率-产氢率曲线本质上是非线性的但实际项目里通常在一定负荷区间内做线性假设给出额定功率和额定效率部分论文用分段线性或带整数变量的启停模型表示电解槽的最小技术出力限制。在MILP框架里这是完全可以实现的只需引入一个二进制变量表示启停状态。2.3 多能源耦合的瓶颈环节在哪里真正让这套系统复杂起来的不是单个设备建模而是耦合约束。电、热、气不是三条独立的管道它们通过设备互相转换光伏/风电在电力母线上汇集供给电负荷与P2GP2G产气汇入气负荷或储气罐光热的储热罐同时服务汽轮机发电和热负荷ORC又把系统中的某些中低温热流转化为额外的电流。因此调度模型里会同时出现如下交叉变量电母线平衡方程里等式右侧既有常规机组出力、光热电站电出力、ORC电出力、新能源出力左侧则有电负荷、P2G耗电热母线平衡里光热电站可调节热输出、储热罐充放热、ORC热端取热等形成热量分配关系气母线中P2G产气量、外部购气量、气负荷等构成平衡。这些耦合关系最终全部落地为优化模型里的等式约束或不等式约束。写代码的时候必须规划好变量命名空间和维度否则后面出错的概率会非常大。我见过很多同学在Matlab里建模型时过一个月回来看代码就看到变量顺序混乱这个在后面第4部分通过代码案例详细展开。3. 目标函数与约束体系设计从物理到数学3.1 两个层次的目标函数怎么取舍优化调度按决策目标不同通常分成两大类。类一以运行成本最小为目标适用于园区、局域电网等经济性敏感的运营场景目标函数由“购电费用 燃料成本 启停成本 弃风弃光惩罚 设备运维成本 购气费用”构成。类二以新能源消纳最大或碳排放最小为目标在政策驱动场景下用得比较多。实际工程项目里经济目标与环境目标往往通过加权系数合在一起。不过对于初学者我建议先做成本最小模型因为单位统一元结果容易量化检查等代码跑通了再加碳约束并不难。目标函数写出来后需要给每一项定义清楚系数的物理意义比如新能源弃电惩罚是不是应该高于其他发电成本如果惩罚设得太低优化结果里可能出现“为了让成本表好看而主动废掉本来可用的新能源出力”的情况这明显不符合调度初衷不过惩罚设得太高又可能扭曲正常购电/购气决策。这类系数的敏感性分析要多跑几次才能找到合理的区间。我常用的经验是单位弃电惩罚电价取系统最高购电电价的1.2~1.5倍左右可以保证新能源优先消纳又不至于让结果失真。3.2 设备约束要细化到什么程度设备约束是整个模型里最琐碎也最考验细心的地方。系统里每一个可调度设备都要回答四个问题出力上下限是多少爬坡速率限制是多少运行会不会有最小启停时间限制储能相关量的状态转移方程和容量约束怎么表达有一说一很多刚接触优化的同学容易犯一个错误就是只写上下限和平衡约束漏掉储能状态量相邻时段耦合的约束导致储热罐在凌晨1点放的热可能凭空“瞬移”到下午3点再用——这在数学里可能没问题物理上是假的。具体到CSP储热罐我认为最关键的约束是储热罐当前蓄热量的状态转移方程S_t S_{t-1} η_ch * P_ch_t - P_dis_t / η_dis - S_loss其中S_t表示t时段末储热罐内的热能量MWhη_ch和η_dis分别是充放热效率P_ch_t和P_dis_t是t时段的充放热功率S_loss是散热损失。注意这里的充放热效率和功率都是写在“热侧”的光热电站的电出力由“蒸汽发生系统从储热罐取热集热场直接供热”驱动汽轮机热-电转换关系也要相应地建模。ORC机组的约束要分清它的热源来源。在综合能源系统里通常存在多种可能的余热源我们要给ORC配置一个“热源取热上限”同时发电功率P_orc_t η_orc * Q_orc_in_t其中Q_orc_in_t为t时段ORC从热源取的热功率这样约束简单但足以体现ORC运行的可行性边界。还有需要注意的是ORC冷凝器排出的余热是低温热不太可能回收利用在建模时不需要画蛇添足把“冷凝废热”当成产出除非你的系统专门设计了梯级利用链条。P2G的模型相对直接一些输入是电功率输出是氢气流量或计入气网的天然气等效流量。有些模型还要加一个氢气储罐用储氢罐状态来平衡P2G产氢和下游用氢/售氢之间的时移。在这套系统中P2G更常见的定位是消纳晚间风电或午间光伏高峰把本来要被限制的电力变成气体储存起来——这也意味着它的运行特征与新能源出力曲线之间存在强烈的时段耦合特性只能通过优化求解来确定具体的运行时刻。3.3 新能源出力与负荷曲线怎么处理调度问题永远都是一个“给定外界条件求最优决策”的过程风电、光伏的预测出力以及电/热/气负荷曲线属于调度模型的外生输入不是优化变量。但外生输入的质量直接影响模型结果。工程实操中需要把一整年的新能源出力和负荷数据聚类成几个典型日场景比如春夏秋冬四季或者丰水期/枯水期。每个典型日取24时间断面作为输入向量。如果没有实测数据可以用Matlab里的现有工具生成模拟曲线光照强度按正弦曲线加云层波动风速用韦伯分布抽样再做时序平滑负荷数据可以参照典型用户曲线叠加随机扰动。原始数据输入前必须做归一化或者去量纲处理。在优化求解层面我注意到如果不把数量级差异控制在合理范围比如电负荷是几百兆瓦时储能罐容量只有几百千瓦时Cplex/Gurobi这种求解器很容易报数值问题。经验做法是乘以一个基准值将功率统一折算到标幺值域中比如以1MW为基准那么100MW负荷变量量纲变为100储热罐蓄热量可能变成35之类数值范围就不会差别太大。求解器的数值稳定性会好很多这一步对后期调参来说绝对是事半功倍的。4. Matlab代码实现与调度模型落地指南4.1 代码框架搭建与Yalmip建模流程Matlab里建立一个确定性优化模型当前最主流的方式是先装好Yalmip工具箱再加一个求解器。Yalmip并不真正执行求解它像一个伟大的翻译官把人类易于理解的建模语句翻译成求解器能读懂的数学规划标准格式。只要你的模型是线性、混合整数线性、二次约束甚至非线性Yalmip都有对应接口去对接合适的求解器。在这个场景里我推荐用Gurobi或Cplex因为白天遇到MILP问题时它们的内点法和分支切割法表现都很成熟。需要提醒的是Yalmip的安装并不复杂简单来说就是下载yalmip压缩包后把文件夹加到Matlab搜索路径里。真正要花时间的是求解器license的配置尤其很多新手同学在学校电脑上用Cplex经常报出Failed to get license或Error: Unable to solve because solver is not instlled这种时候通常不是模型写错了而是求解器接口没配对。应该逐个检查求解器是否安装、环境变量是否设置、Matlab里调用solvesdp或optimize时是否明确指定了参数。避免因许可证或者路径问题影响进度。Yalmip建模的典型代码结构我通常按下面几个区段来组织参数初始化区设置系统参数、读取负荷和新能源预测曲线决策变量定义区用sdpvar定义连续变量用binvar定义二进制变量约束条件声明区把设备模型中的物理公式一行一行写进去目标函数声明区按成本项累加求解与结果输出区调用optimize并绘图。最忌讳的是从头写到尾只有一个大脚本变量作用域混在一起出了错不好定位。我自己实践下来比较舒服的方式是每个功能模块单独写一个函数但必须注意函数里创建的Yalmip变量返回后仍然保持一致因为它本质上是Matlab对象引用关系比较特殊容易踩坑。4.2 光热储热、ORC和P2G的代码怎么组织下面展示几个核心设备的建模片段来说明关键约束在Yalmip里的表达方式。先看光热储热罐与电出力部分的基础逻辑%% 光热电站相关变量 P_sf sdpvar(1, T); % 集热场接收的太阳热功率, MW P_gross sdpvar(1, T); % 汽轮机发电功率, MW P_ch sdpvar(1, T); % 储热罐充热功率, MW P_dis sdpvar(1, T); % 储热罐放热功率, MW S_soc sdpvar(1, T1); % 储热罐蓄热量MWh取时段初/末状态 P_res sdpvar(1, T); % 实际集热场输入到蒸汽发生系统的热功率 % 集热场输入可用太阳热功率受DNI影响上限随时间变化 Constraints [Constraints, 0 P_sf P_dni]; % P_dni是t时段光场可收集的热功率上限 Constraints [Constraints, 0 P_res P_dis Q_sg_max]; % 蒸汽发生系统吸热上限 % 发电功率与吸热功率成正比例 Constraints [Constraints, P_gross eta_pb * (P_res P_dis)]; % 储热罐状态约束 Constraints [Constraints, S_soc(1) S_init]; for t 1:T Constraints [Constraints, S_soc(t1) S_soc(t) * (1 - loss_rate) eta_ch * P_ch(t) - P_dis(t) / eta_dis]; Constraints [Constraints, 0 P_ch(t) P_ch_max * u_ch(t)]; Constraints [Constraints, 0 P_dis(t) P_dis_max * u_dis(t)]; Constraints [Constraints, u_ch(t) u_dis(t) 1]; % 不能同时充放热 Constraints [Constraints, 0 S_soc(t1) S_max]; end这里引入了一个必要的二元变量u_ch和u_dis表示充/放热状态“不能同时充放热”是物理层面的必需约束。在建模时如果省略那么优化可能会倾向于在同一时段用高温熔盐去充热同时从低温罐取出熔盐放热造成无意义的热量循环却不产生电能导致目标函数值失真。再来一个ORC的简单模型假设热源来自系统中某一固定废热流分时温度变化导致效率变化%% ORC 有机朗肯循环约束 P_orc sdpvar(1, T); % ORC电出力 Q_orc_in sdpvar(1, T); % ORC从热源侧取热功率 eta_orc_t interp1(temp_profile, eff_table, temp_source_t, linear); % 由热源温度查表 Constraints [Constraints, Q_orc_in 0]; Constraints [Constraints, Q_orc_in Q_orc_in_max]; Constraints [Constraints, P_orc eta_orc_t .* Q_orc_in];如果ORC的热源来自光热电站的低温旁路那么取热功率也会出现在光热系统热平衡里需要把这项加到集热场/储热罐的 热量分配关系中。P2G模块的核心代码相对简洁不过当电解槽需要启停状态时还得加上最小出力约束这一点往往被不少研究者忽略%% P2G 电制氢模块 P_p2g sdpvar(1, T); u_p2g binvar(1, T); % P2G启停标志 F_h2 sdpvar(1, T); % 氢气流量kg/h或m3/h for t 1:T Constraints [Constraints, P_p2g_min * u_p2g(t) P_p2g(t) P_p2g_max * u_p2g(t)]; Constraints [Constraints, F_h2(t) eta_p2g * P_p2g(t)]; % 线性效率 end需要注意的是如果有储氢罐那么储氢罐的状态约束与储热罐类似会增加一组S_h2的状态变量和充放氢变量。如果模型只允许P2G产出的氢气直接进入氢负荷或天然气管网那么可以不引入储氢罐平衡方程会简单很多——具体选择取决于你研究的问题是否要含氢气长周期存储。4.3 连接能量母线系统级平衡约束怎么写设备建模完成后必须将所有设备的端口变量汇聚到电、热、气三条能量母线的平衡方程中才能让模型真正形成一个“系统”。在我的代码里这些平衡约束的形式通常非常直观%% 电力母线平衡 for t 1:T Constraints [Constraints, ... P_buy(t) P_wt(t) P_pv(t) P_gross(t) P_orc(t) ... P_load(t) P_p2g(t) P_eb(t)]; % P_eb是电锅炉或其他电负荷 end %% 热力母线平衡此处以区域热网所需热量为热负荷光热供汽轮机后剩余热量供应热负荷 Constraints [Constraints, P_heat_supply P_gb P_eb H_load]; %% 天然气/氢气母线平衡 Constraints [Constraints, F_h2 F_buy_gas G_load F_boiler_gas];这里需要特别强调一点电平衡里各项的正负号非常容易出错。买电、新能源出力、发电设备输出是在方程左侧的“源”项电负荷、P2G的耗电则应放在右侧即“荷”项。写反了很容易造成优化结果产生难以置信的“凭空发电”或“凭空消失的电”遇到这种问题可以先检查各项的正负符号。目标函数部分按成本最小形式配置Objective sum(P_buy .* price_buy * dt) ... % 购电成本 sum(F_buy_gas .* price_gas * dt) ... % 购气成本 sum(P_gross .* c_om_csp * dt) ... % 光热运维成本 sum(P_orc .* c_om_orc * dt) ... sum(P_p2g .* c_om_p2g * dt);dt是调度时段长度如果时间颗粒度为小时那么dt1功率乘以时间等于能量乘以电价才是费用。很多第一次做调度优化的同学会在这一步漏算时间颗粒度导致成本低了一倍或高了一倍最后得到的结果无论解释成什么都别扭。4.4 求解器配置与小规模算例结果解读代码写完直接调用optimize时建议显式指定求解器参数。例如用Gurobi就设置options sdpsettings(solver, gurobi, verbose, 2, showprogress, 1); result optimize(Constraints, Objective, options);若模型是纯线性的输出里的gap在零附近solvertime不超过几秒如果条件表达式里埋伏了不必要的非线性函数Yalmip会把它转成二阶锥或非线性模型此时Gurobi往往无法直接求解Yalmip会提示“No suitable solver”。遇到这种情况最简单的办法是检查使用的函数是否引入了非线性比如用乘法c * x没问题这两个是1*n向量点乘但你如果写成x * c且c也是变量就是二次项了。另外绝对值如果不想引入非光滑项可以用变量对偶法或引入辅助变量线性化表达。小算例做出来后通常可以画出三种图一是典型日全天各电源的电力出力堆叠图二是储热罐24小时蓄热状态与P2G运行功率的对比图三是在不同场景下成本项的柱状图。通过这些图就能够验证模型行为是否合理。如果发现储热罐只在午后充热、晚高峰放热那么说明储热的经济价值被正确捕捉了如果发现储热罐压根没被调度使用SOC恒定不变大概率是目标函数里某些成本项设置过大或储热容量本身过大/过小需要回到参数设置排查。5. 常见问题与排查技巧实录5.1 模型不可行或不收敛先查这几处在调度优化里最崩溃的报错就是Infeasible problem。由于约束太多到底哪个限制条件出了问题求解器往往不会直接指出来。我常用的排查思路是按“从简到繁”的顺序来定位第一步先把所有储能相关的约束注释掉看系统本身平不平衡。如果不行大概率是功率平衡方程或热平衡方程中某项不能同时满足单位或许有问题kW/MW混淆第二步只保留储热约束检查储热罐在初始和终末SOC相同的条件下是否能够在各个时段完成充放。这是验证调度模型正确性的关键一步。很多模型强制要求S_{T1} S_1以方便进行多日滚动调度但在给定初始状态的情况下容易导致无解建议先把循环调度约束放宽成不等式比如末时刻SOC要高于某阈值即可第三步再把所有约束放回来从目标函数中逐项删除来观察删除哪一项后模型可行通常能找到是哪类量卡的太紧。另外把所有设备启停的二元变量都设为常数1看连续松弛问题是否可行可以快速判断是否存在零/一整数的组合死角。比如你要求储热罐不能同时充放但一天中储热罐需要先在一定时段满放再满充而约束又把充放热上限限制在一个极小的值就可能找不到离散方案的可行解。5.2 Yalmip建模常见报错自查表我把平时见到的高频报错整理成一个表格大家可以直接对照排查。虽然不是绝对覆盖但至少覆盖了80%以上的入门问题常见现象可能原因解决建议No suitable solver模型包含非线性项但求解器不支持用yalmiptest检查已安装求解器用export导出模型查看类型Infeasible problem平衡约束无法满足且模型被过约束逐模块放行排查不可行约束来源用relax看对偶变量找出弱约束Index exceeds matrix dimensions读入的负荷曲线长度不足24点观察输入数据维度用size(load)检查求解时间异常长二元变量太多分支定界计算量爆炸削减整数变量的范围把某些离散状态处理为设计参数结果出现“充放电循环但不产生真实吞吐”储能约束松弛或漏约束确认加入了同时充放不可行约束及能量状态转移方程目标函数值是NaN初始值含NaN或除以零清洗所有输入数据检查效率参数是否为零5.3 数值调优与加速的个人心得解决完“能跑”的问题之后接下来就是“跑得快”和“算得稳”的环节了。这里分享几个最直接的加速技巧避免在约束循环里重复创建新的Yalmip变量。比如Constraints [Constraints, ...]在一个大循环里反复拼接其实效率不差但如果是多维大矩阵建议一次性把整个变量矩阵定义出来用repmat或矩阵运算一次性生成约束能大幅缩短建模时间变量的取值范围尽量收紧。给每个决策变量写上下限边界Gurobi/Cplex做预处理时会把冗余约束删除加速比较明显。比如储热罐SOC如果明确知道最高只有12MWh就不要定义到100MWh免得数值比例失调对电力平衡和热平衡这类核心约束尽量整行矩阵化写不做逐一循环。例如电平衡可以用sum(A,2)这种向量化操作直接形成等式约束Yalmip内部会把它识别为一个大稀疏矩阵约束求解时的预处理开销小得多。我实际跑过的模型规模大概是24时段、十余台设备、储能罐两个、带少量启停二元变量总变量数量在四五百个左右Gurobi解决时间大约是几百毫秒到几秒。如果你的模型跑了几分钟还没结果一般情况下不是设备数量多而是有地方产生的二次约束或整数变量较多。这种模棱两可的问题建议先冷静分析模型本身而不是盲目增大求解器时间限制。5.4 调度结果里那些“看起来奇怪”但真实合理的现象有时跑完结果会有一些输出让人觉得哪里不对劲。比如P2G在某个时段明明是电价低谷却没有开起来。原因可能是P2G的技术最小出力限制导致这一个小时内制出的氢气没有地方可去——要么是因为下游气负荷此时段很低要么储氢罐已满即使电价低强开P2G只会增加弃氢风险。这就提醒我们看到结果不能只凭直觉判断要结合存储状态的完整曲线来分析。还有一种情况是光热电站傍晚把储热基本放完但夜间1点的热负荷靠电锅炉顶着看起来好像不如“晚点放热”更划算。为什么优化结果是这样因为模型如果在傍晚低负荷时段不发电就需要从外部电网买电购电价格可能远高于储热发电的边际成本所以它宁可傍晚集中放热把节省的购电费放到最大化。对这类现象需要用“边际成本”的语言去解释以验证结果的合理性。6. 从案例到实践整套系统的可扩展方向与个人经验代码跑通的那一刻是很开心的但综合能源系统的优化调度从来不是从零开始完整复现一遍就结束的研究课题。在实际工程或论文延展里真正的价值在于把基础模型往不同方向扩展。对这个含光热、ORC和P2G的模型来说最容易也最有价值的扩展方向有这样几个多时间尺度滚动调度。上面阐述的是单典型日24小时的确定性优化在实际运行调度时需要做日内滚动修正。例如用第1~4小时实际数据校正新能源预测误差然后每4小时滚动优化一次重新求解未来24小时的决策。这个在Matlab里可以通过循环调用同一套模型实现对代码架构的复用很有要求如果你的约束和变量封装成了独立函数扩展滚动调度并不复杂。考虑不确定性。光伏和风功率预测不可能绝对准确常见的处理是采用场景法把风电场景集生成多个可能序列再把目标函数改成期望值最小或者采用鲁棒优化描述成一个min-max-min三层结构。在Matlab里随机规划通常和scenario配合鲁棒优化则用uncertain对象配合Yalmip。光热辐射本身的波动如天气变化引发DNI突变也值得纳入不确定性集合。我做过的模型里有一个例子把DNI风速负荷构建为盒式不确定集在Yalmip里实现两阶段鲁棒优化很考验建模和数值功底但是做完后对整个系统的弹性会有更深的理解。引入碳交易机制。如果目标是“双碳”相关的综合能源系统研究那么把碳排放配额和碳交易成本加入目标函数是自然的扩展。光热和P2G在这方面有天然优势光热电站比燃煤机组省大量燃料碳排放P2G制氢在部分场景下替代灰氢也有减排价值。碳价价格波动可能反而会改变调度结果中光热和ORC的优先序列这个方向特别适合做成敏感性分析。再有一个小而实用的小建议在代码开发过程中及早形成标准的绘图模板。我之前吃过亏——因为Data结果导来导去找不到哪个变量对应哪条曲线耽误了不少时间。建议把所有结果变量整理成一个结构体变量result.P_csp、result.S_soc、result.P_orc、result.F_h2画图也写成统一函数后期加任何设备只要在结构体里补一个字段就能自动出图作为校核。这样整个项目几个人协作也会比较顺畅。最后分享几次调模型积累下的体会综合能源优化调度的核心难点从来不在某一个公式的复杂程度而在能否把不同时间尺度、不同能量品位的物理过程统一到一个边界清晰、数值稳定的数学框架里。面对这类模型与其贪多求全把每个部件都写到极致不如先把关键设备的物理约束抓准、系统能量平衡关系写对、目标函数经济含义理清。保持“先跑通、再跑稳、后跑快”的推进节奏往往比一次想写出完美代码更有效。如果在这过程中遇到“反直觉”的调度结果不用急着怀疑模型出错了——先在设备约束、能量平衡和价格机制几个维度对照分析很多时候会得到比书本案例更有意思的发现。
