移动储能预布局与动态调度：两阶段鲁棒优化提升配电网韧性

做配电网韧性的同学应该都有这种体会：台风、冰灾这类极端事件一来，静态的应急电源车和固定储能往往只能“头疼医头”，要么提前蹲错了位置，要么抢修时利用率不高。前一阶段我完整复现了一套移动储能预布局与动态调度策略，走的是两阶段鲁棒优化的路子，测试算例用的IEEE33节点系统。整套代码跑通之后最大的感受是：移动储能能不能真正撬动配电网恢复能力，关键在于“预测—布局—响应”这条链路能不能闭环。这篇内容就是把我复现的模型、求解思路、代码结构和踩坑记录全部梳理一遍，给正在做配电韧性、储能调度相关方向的朋友一个完整参考。

1. 项目背景与核心问题拆解

1.1 为什么移动储能是韧性提升的关键角色

配电网韧性问题，说白了就是电网在极端灾害下“扛得住、恢复快”。传统配电网里，应急手段无非是固定储能、柴油发电机、抢修队伍这些。但固定储能有一个天然问题：位置锁死了，灾害发生路径和负荷损失位置不匹配时，储能就算容量再大也送不出去；柴油发电机响应慢、污染大，调度灵活性也差。移动储能（MES）就不一样了，它本质上是“储能电池+运输载体”，既可以像固定储能一样提供电压支撑和功率补偿，也可以在灾前根据灾害预报移动到关键负荷附近，在灾后根据实际故障动态转移位置。

这套策略的核心价值在“预布局”这三个字上。灾害之前，移动储能的位置和容量分配还有足够的时间窗口，可以利用气象预报信息和历史负荷数据去做优化决策。等到灾害真正发生、故障线路暴露出来之后，再去调度移动储能，其实已经晚了半拍。所以预布局解决的是“把正确的资源放在正确的地方”，动态调度解决的是“资源到位之后怎么最大化发挥恢复效果”。两者耦合在一起，就是一个典型的两阶段决策问题。

1.2 两阶段决策的结构逻辑

第一阶段（预布局）发生在灾害尚未到来时，决策变量是移动储能从哪些初始位置出发、移动到哪些节点、每个节点配置多少容量。这个阶段的决策要考虑最坏情况下的故障场景，因为灾害实际造成的线路故障是不可预知的，但可以通过鲁棒优化把最恶劣的场景纳入考量。

第二阶段（动态调度）发生在灾害发生、故障场景清楚之后，决策变量是具体的放电功率、充电功率、移动路径、节点连接状态。这时候目标很直接：以最快的速度恢复关键负荷，最小化停电损失，同时保证系统运行约束不越限。

这种“先预布局、后动态调度”的结构，和我们平时处理不确定性问题时的思路高度一致。电网不只配电领域，包括输电线路抢修、燃气系統应急、供水管网调度，凡是涉及“灾前准备+灾后响应”的资源分配问题，几乎都可以套用这个两阶段框架。它最大的价值是把不确定性的影响通过场景和最坏情况的方式显式建模，而不是拍脑袋决定储能放哪里。

2. 模型构建与数学表达

2.1 不确定性建模方式

这套策略里，不确定性主要有两个来源。一个是负荷波动，尤其是关键负荷在灾害前后的变化；另一个是故障场景，也就是哪些线路会断开、哪些区域会失电。

处理故障不确定性的常见方式有两种：场景法和鲁棒优化。场景法是把历史台风、冰灾数据整理成几个代表场景，每个场景配一个概率，缺点是要覆盖足够多的极端情况，场景数往往要很大；鲁棒优化则是用一个不确定集合来刻画故障可能出现的范围，不需要概率信息，更适合数据不足但需要保证安全性的场景。我复现的这套策略采用的是鲁棒优化的思路，用盒式不确定集合来描述故障线路的变化范围，保证预布局方案在最恶劣情况下都能实现负荷恢复。

这种建模方式的物理意义很清晰：我们不需要精确预知“哪条线会断”，只需要知道“最多断几条”“哪些区段受损风险高”，系统就能给出一个备选方案。实际做韧性分析时，这正是工程上最容易获得的信息。

2.2 目标函数设计

目标函数衡量的是韧性提升效果，通常用灾害后系统总失负荷量最小化来表示。形式上是一个最小化问题，分为两个阶段累加：

第一阶段的目标是让预布局方案在后续最坏场景下具备最优的恢复潜力，所以目标函数内嵌了第二阶段的最小化失负荷。第二阶段的目标才是真正执行层面的最小化。

有点绕，但捋清楚之后逻辑就顺了：这是一个 min-max-min 结构。外层 min 是第一阶段的预布局决策，中间的 max 是寻找最恶劣故障场景，内层 min 是在最恶劣场景下求最小失负荷。目标值越小，说明预布局方案越能有效应对极端故障。

除了失负荷量，部分方案也会把移动储能投资成本、运输成本、网络损耗纳入目标函数做多目标加权。实际复现时，如果追求的是复现论文效果，保持纯失负荷量最小化就够了；如果想迁移到实际项目，再考虑加成本项不迟。

2.3 关键约束条件

约束条件是这套模型最核心也最容易出错的部分，主要包括五类：

• 配电网潮流约束：采用DistFlow线性化模型，处理有功、无功、电压幅值之间的关系。这是可以在YALMIP里快速实现的线性约束集合，也是很多配电网优化问题的标准基础。
• 移动储能运行约束：包括充放电功率上下限、SOC更新方程、容量限制、以及“同一时刻只能在一个位置”的逻辑约束。这里特别要处理移动时间约束——储能从节点i移动到节点j需要时间，在灾害窗口内不是瞬间完成的。
• 网络拓扑约束：预布局完成后，移动储能接入节点和配电网拓扑结构要保持一致，不能出现孤岛和断开节点在电气上储能却连不上的矛盾。
• 不确定性集合约束：定义故障场景的变化范围，例如“最多断开的线路数不超过某阈值”或“故障线路集合属于某候选集合”。
• 线性化处理引入的辅助变量约束：涉及二进制变量（0/1变量），用于表示储能是否接入某节点、线路是否故障、负荷是否恢复等。

这些约束需要特别注意SOC更新方程中时间步长的处理。我做复现的时候，一开始把移动过程的时间尺度和调度周期（通常1小时或15分钟一个时段）用同一个时间步长来建模，但实际移动时间可能不是整倍数，就会导致SOC跳变出现错误。后来改成在移动过程中加入修正项，才把物理过程对齐了。

3. 求解算法与Matlab代码架构

3.1 两阶段鲁棒优化与C&CG求解框架

两阶段鲁棒优化问题直接求解很困难，工程上通常用列约束生成算法（C&CG）来迭代求解。核心思想是：把问题拆成一个主问题和一个子问题，主问题求解预布局方案，子问题固定预布局方案之后去寻找最恶劣故障场景，然后把找到的恶劣场景以“列”的形式加回主问题，再重新求解。这样反复迭代，直到主问题和子问题的最优值收敛到同一个目标值。

用大白话解释这个过程：先随便给一个储能位置方案，然后问“这种情况下最差会断哪几条线”，找到了最差断线组合之后，重新优化“储能应该挪到哪些位置才能对付这个最差组合”，挪完之后再看“新的最差组合是什么”，以此循环。C&CG的优势是收敛快、割平面有效性强，在配电网这类中等规模问题上表现非常稳定。

不过这里要提醒一句：如果故障不确定性集合是非离散的，比如考虑线路故障概率和持续时间并用连续区间描述，C&CG的子问题可能就不是单纯的线性规划，而需要加入整数变量来处理线路投切状态。这套策略里线路故障通常建模成0/1变量，所以子问题实际是一个混合整数线性规划，需要调用Gurobi或Cplex这类商业求解器来解。

3.2 代码模块划分

整个Matlab工程的代码结构，我复现时按功能拆成了六个模块，维护起来很舒服：

• 数据初始化模块：定义IEEE33节点的线路参数、节点负荷、网络拓扑、移动储能数量/容量/速度/充放电功率参数。
• 主程序模块：搭建C&CG主循环，设置收敛精度和最大迭代次数。
• 主问题建模模块：用YALMIP建立预布局决策变量和目标函数，加入第一阶段及与第一阶段相关的约束。
• 子问题建模模块：固定第一阶段结果，建立最恶劣故障场景搜索模型，求出故障集合对应的最小失负荷量。
• 结果输出模块：输出两阶段最优目标值、储能位置方案、各时段充放电功率、SOC变化曲线、节点电压与负荷恢复状态。
• 绘图模块：绘制配电网结构图、韧性曲线、储能转移路径图、失负荷柱状图。

3.3 求解器配置与调试要点

YALMIP是Matlab下最常用的优化建模工具箱，配合Gurobi求解MILP表现很好。我的建议是优先用Gurobi，因为C&CG迭代过程中需要反复求解子问题，Gurobi的求解速度比Cplex快不少，而且在Windows系统下部署也更省心。安装好了之后在Matlab里跑一下yalmiptest确认求解器能被识别。

求解器配置上，有两点实测经验：第一，C&CG的主问题和子问题都需要设置求解器为MILP，并让Gurobi启动多线程，代码里加上ops = sdpsettings('solver', 'gurobi', 'gurobi.Threads', 4, 'verbose', 0)，迭代会明显提速；第二，收敛精度不宜设得太严，通常主问题目标值和子问题目标值之差小于1e-3即可，设成1e-6会让迭代次数增加不少，但收益很小。

4. 仿真案例分析：IEEE33节点复现过程

4.1 算例参数与场景设定

IEEE33节点系统是我在配电网优化仿真中用得最频繁的标准算例。33个节点、32条支路、5个联络开关，基准容量10 MVA，基准电压12.66 kV，总负荷大约3.7 MW。算例规模不大，但拓扑特征丰富，既有辐射结构又有联络通道，非常适合测试移动储能预布局策略。

复现时我按论文常见配置设计了这样一组场景参数：

实际运行之前我先把广度搜索验证了一遍：把4条故障线路的所有组合枚举出来，对比C&CG求出的最恶劣场景是否落在枚举集合内。这一步很关键，因为C&CG本质上是在找最恶劣场景，如果找出的场景不在合理候选集中，说明子问题建模或不确定集合定义有问题。

4.2 运行结果与图表解读

代码跑完之后，C&CG在4次迭代内收敛，目标值在迭代过程中的变化是：主问题目标从初始的失负荷较高值逐步上升（因为每加入一个割平面，第一阶段方案就要多考虑一个恶劣场景），子问题目标从较低值逐步下降，最终两条曲线逼近到同一个数值。从实测结果看，移动储能预布局方案把最恶劣场景下的失负荷量降低了约35%，对比不采用移动储能、仅依赖固定储能的方案提升显著。

这里有一个观察值得分享：在初始迭代中，主问题给出的预布局位置倾向于集中在负荷集中区域附近。但随着迭代进行，子问题不断加入新的故障场景，布局方案会逐渐向网络的关键转供节点移动。这背后的逻辑很直观——台风路径不确定时，储能放在负荷中心虽然第一时间能覆盖最大负荷，但如果故障线路刚好切断负荷中心的上游电源，储能孤立无援，恢复效果反而打折扣；放在转供节点，则能灵活支援多个方向。这就是“鲁棒”两个字在结果上的直接体现。

4.3 方案对比与灵敏度分析

做复现时我额外做了一个灵敏度测试，观察移动储能数量对韧性提升的影响。分别用1台、2台、3台储能重复C&CG求解，失负荷量分别下降了约15%、35%、48%。可以看到边际效益在递减，说明工程上并不是储能越多越好，需要结合投资成本做权衡。

故障线路数量上限从2条逐步调到6条时，预布局方案也从“集中在关键负荷周围”逐步变为“分散在多个转供节点”，失负荷量下降幅度逐渐收窄。这个现象说明模型对故障不确定性的敏感度很高，实际应用中要合理评估故障范围——定得太宽，预布局方案会过度保守；定得太窄，方案可能在更严重的灾害下失效。

5. 实操经验与代码使用指南

5.1 运行环境与依赖安装

我复现时用的环境是Matlab R2022a + YALMIP最新版 + Gurobi 10.0，操作系统Windows 11。如果你手头没有商业求解器，也可以先用intlinprog跑一遍小规模算例验证正确性，但迭代速度会明显慢。正式做仿真建议还是上Gurobi或Cplex，差别在分钟级别和秒级别之间。

安装YALMIP很简单，去官方仓库下载，把整个文件夹加入Matlab路径就好。Gurobi则需要在官网申请个人学术许可，注册一下就能拿到。装好后在Matlab命令行输入yalmiptest，看到Gurobi列在可用求解器里，就可以开始跑了。

5.2 参数调整指南

复现这套代码时，需要重点调整的参数主要在三块：

• 网络参数：IEEE33节点的线路阻抗、负荷数据在数据初始化模块里以表格形式给出，做其他配电网算例时，把这里的表格替换成对应数据即可。
• 移动储能参数：储能容量、充放电功率、SOC初值、移动速度、运输时间。需要注意储能台数改变后，预布局变量的维度也要对应的逻辑同步修改。
• 不确定集合参数：故障线路数量上限、故障候选线路集合。如果你的场景没有线路故障，只有负荷不确定性，可以把这一块改成负荷不确定集合。

调参的时候建议一次只改一个参数，改完跑一次迭代曲线对比。别同时改储能台数和故障数量，否则结果变化时很难定位是哪个参数引起的。

5.3 代码验证与结果一致性校验

很多同学复现论文代码时会被“代码跑通但结果对不上”折磨。我的经验是分三步验证：

第一步，先把鲁棒优化退化为确定性优化，也就是把故障数量上限设为0或者固定故障集合，这时候两阶段问题退化成一个普通MILP，可以单独验证目标函数和约束是否与其他工具结果一致。

第二步，把移动储能数量设为0，此时模型退化为纯配电网故障恢复问题，可以和IEEE33的标准潮流结果对比电压和功率分布，确认DistFlow方程实现正确。

第三步，再开启完整的C&CG迭代，对比论文报告的目标值。如果前两步都是对的，第三步通常不会有大偏差。如果第三步对不上，大概率是收敛精度设置或不确定集合定义和原论文有出入。

6. 常见问题与排查技巧实录

问题1：C&CG迭代发散或不收敛。

排查方向主要是子问题中加入了未知的二进制变量或非线性约束。如果子问题在固定主问题结果后仍然包含主问题变量，就会导致割平面失效。解决方法是检查子问题约束中是否遗漏了固定主问题变量的赋值语句，尤其是使用YALMIP的value()或assign()函数时，要确保赋值对象和子问题表达式完全匹配。

问题2：SOC越界或充放电功率不平衡。

这种问题多半出在SOC更新方程上。移动储能在节点间转移需要时间，这段时间内SOC可能保持不变，但代码中常常因为没有单独建模移动时段，导致SOC出现不应有的变化。我采用的解决方法是把移动过程单独设为一个不充不放的时段，或者用一个“移动时间修正量”把跨节点转移的时段补偿掉。Debug时把SOC曲线画出来，出现台阶式跳变就说明对时间步长的处理有bug。

问题3：Gurobi求解速度极慢。

先看变量数量是不是被不必要地放大了——比如把上万个辅助变量建立成二元变量，就会拖慢求解器。可以先用size检查YALMIP变量的维度，看看是不是扩展成了预期规模的多倍。如果问题规模没问题，检查收敛精度是否太严密，把gurobi.TimeLimit设为600秒辅助判断问题是否能在可接受时间内收敛。

问题4：结果对故障场景不敏感。

如果不管故障数量上限怎么调，预布局方案都几乎不变，很可能是不确定集合的候选故障线路被限定在了一个很小的范围内，常见原因是在建立不确定集合时只选了与某种负荷模式强相关的少数线路。这时候把候选故障线路的选取逻辑改为基于拓扑脆弱性排序（比如关键联络线、长线路优先），结果就会明显变化。

问题5：画出的配电网结构图节点错位。

IEEE33节点的坐标数据网上有很多版本，部分版本不一致会导致图乱。建议使用标准文献中的坐标值，这类数据在Graphviz绘图时效果最稳定。如果是自己写绘图脚本，注意节点的连接关系要和branch数据完全一致，否则会出现线路交叉和节点悬空。

另外，整个调试过程中我感受最深的一个坑，就是YALMIP中implies函数在整数变量使用时会引入大量的额外二元变量，导致模型规模爆炸。我自己在建模配电网辐射条件时，一开始用了大量的implies来约束线路投切和节点电压，结果模型变量数量翻了好几倍，求解时间从几十秒变成几十分钟。后来把这个逻辑换成用明确的0/1变量配合线性大M约束来表达，问题规模瞬间降下来了。这也是我在复现各类配电网优化论文代码后总结出来的一个通用经验：能用大M法解决的逻辑约束，尽量不要用implies系列函数。

还有一个细节值得单独提醒，就是IEEE33节点系统本身带联络开关，预布局阶段如果没考虑联络开关的作用，储能的位置优化结果会在灾后调度阶段发生较大偏移。因为联络开关在故障后可以提供备用电源通道，移动储能的接入位置应当和这些通道配合起来，而不是孤立地只看故障前拓扑。所以代码里如果有网络拓扑参数表，一定要把联络开关的支路数据检查清楚，别把本来是常开状态的联络线当成了常闭运行线路。

至于这套代码后续的扩展方向，我目前在做的是把单一灾害场景扩展成多时段灾害演进过程，让移动储能既能在灾前移动，也能在灾害持续过程中逐步调整位置。如果你对分布式鲁棒优化更感兴趣，也可以把当前的两阶段鲁棒模型改造成分布鲁棒模型，在不确定集合中加入概率矩约束。这些都是基于同一套代码框架可以快速扩展的方向。

复现一套优化论文的代码，本质上是在还原作者的每一层决策逻辑。这套移动储能预布局与动态调度策略，层次很清楚：先想不确定性是什么，再想资源怎么放，最后想新资源怎么用。把这三步想明白，代码实现只是时间和耐心的问题。希望我的这些实测细节能让你少走几步弯路。