计及光伏快速无功响应特性的分布式电源优化配置方法

说句实在话，配电网分布式电源规划这个方向，我见过太多人一开始就把“优化配置”理解成单纯的选址定容——目标函数列个年综合费用，约束里加个潮流方程，再套个智能算法跑几十轮就交差了。这种思路在几年前还行得通，因为那时候DG渗透率低、电网对并网点的要求也宽松。可现在光伏电站大规模接入，调度对快速无功响应能力的要求已经写进了并网规程，如果你还在用“功率因数为固定值”甚至“单位功率因数”的静态模型去做配置优化，算出来的方案投运后很可能会在电压支撑能力上掉链子——要么无功备用不足，要么响应速度跟不上电压跌落的节奏。

这篇就围绕“计及光伏电站快速无功响应特性的分布式电源优化配置方法”展开。文章的主要内容包括快速无功响应到底指的是什么物理过程、怎么把这种动态特性转化为可在优化模型中使用的约束、用Matlab实现时目标函数与变量怎么搭、以及算例结果如何对比和分析。适合用Matlab做分布式电源规划的电气工程学生、做新能源接入评估的工程技术人员，以及对配电网电压控制感兴趣、正在找论文思路的研究者参考。

1. 为什么忽略无功响应特性的配置方案容易“纸上谈兵”

1.1 传统DG优化配置模型的简化假设

先说个常见的现象。翻很多学术论文和开源代码，分布式电源优化配置的模型长得很像：决策变量是待选节点的安装容量（连续变量或离散档位）和安装位置（0-1变量），目标函数要么是年综合费用最小，要么是系统网损最小，约束条件则是节点功率平衡、节点电压上下限、支路潮流限额、DG总安装容量约束这几条。

这类模型看似完整，但内在假设几乎都是静态稳态的：把每个DG当作一个恒功率源（PQ节点），或者最多给定一个固定的功率因数，比如0.95滞后，让它在全天候运行中都按这个功率因数吸收或发出无功。这种简化在DG渗透率不高的时候误差可接受，因为单点注入功率对系统电压分布的影响有限，固定功率因数也有一定的工程可操作性。

可当光伏电站容量动辄几兆瓦、多点接入配电网之后，问题就变了。白天大发时段线路潮流重载，若此时系统发生短路故障或大型负荷突切，光伏电站并网点电压会出现明显的暂态跌落，传统静态配置模型并没有考虑“光伏电站在这个瞬间能不能顶上去”这个能力维度。说白了，配置出来的是一个“平时正好够用、出事没有余量”的方案。

1.2 静态无功平衡与快速无功响应的本质差异

静态无功平衡解决的是稳态运行点上的无功供求关系，对应的时间尺度是分钟级甚至小时级，对应电气量是潮流方程解出的节点电压幅值。而快速无功响应解决的是暂态过程中电压支撑不到位的问题，时间尺度是毫秒到秒级，靠的是逆变器控制环路的快速调节能力，对应的是动态电压支撑特性。

我把两者的差异用一张表来说明，平时讲给学生听也特别直观：

从这个表可以看出，单靠稳态潮流约束只能保证“系统里有多少无功设备”，没法保证“故障瞬间这些无功设备能不能按要求顶上来”。光伏逆变器的快速无功响应如果被配置模型忽略，本质上是把一个动态能力问题和静态规划问题强行解耦了。

1.3 动态响应缺口引发的电压支撑失效场景

我再举一个实践中很典型的现象。某配电网接入了一批分布式光伏，规划阶段是按固定功率因数0.9来做的，正常工作状态下电压确实维持在合格范围。然而一回10kV馈线发生瞬时短路故障，主网电压快速跌落，按理说光伏电站应该进低压穿越并输出无功支撑电压，但由于并网点配套的无功容量是按稳态最大无功需求量的1.1倍设计，动态支撑能力并不可靠，故障期间实际能顶上去的无功只有铭牌容量的60%左右。

结果是：故障切除后电压恢复时间被拉长，局部节点甚至出现低电压脱网，光伏电站因保护动作全部退出，进一步加剧了功率缺额。这种连锁反应在任何静态规划模型里都是看不见的，因为模型里所有变量都处于同一个平衡点，不存在“故障中”“故障后”这些时序状态。

2. 光伏电站快速无功响应特性的物理边界与建模方法

2.1 逆变器有功-无功容量约束：不是简单的PQ上下限

要在优化配置模型里正确体现快速无功响应，第一步是搞清楚光伏电站并网逆变器的无功调节能力边界。很多人随手写一个约束“Q_min ≤ Q ≤ Q_max”，把Q当成独立变量处理，这在数学上很省事，但物理上不够准确。

光伏逆变器的容量一般用视在功率S限制，有功出力P由光照条件决定，无功出力Q必须满足：

[
P^2 + Q^2 \le S^2
]

也就是说，在P比较大的时候，可用的无功调节范围实际上是收窄的。比如一台1MW逆变器，额定视在功率1MVA，正午满发时P接近1MW，Q几乎只能为0，谈不上向电网提供多少无功支撑；反而是早晚P较小、或阴天P受限时，有比较大的无功裕量。所以在优化模型中，无功出力的允许范围应表达为随有功出力变化的运行域，写成约束形式是：

[
-\sqrt{S_i^2 - P_i^2} \le Q_i \le \sqrt{S_i^2 - P_i^2}
]

这类约束是非线性的，直接放进潮流优化会造成凸性变差，后面我介绍用旋转锥约束做二阶锥松弛的办法，可以在保证求解精度的同时顺利把模型交给Cplex/Gurobi这类商业求解器去解。

2.2 电压跌落识别与无功支撑的毫秒级调节过程

快速无功响应到底“快”到什么程度？工程上光伏逆变器通常具备低压穿越（LVRT）能力，在并网点电压跌至0.2倍额定电压甚至更低时仍保持并网，并通过无功电流注入帮助电压恢复。控制上，逆变器检测到电压跌落后的无功支撑动作时间一般要求在几十毫秒以内，这是由瞬时电流内环和电压外环的响应速度决定的。

很多人会问一个实际问题：既然这种动态过程发生在毫秒级，我的优化模型是稳态模型，怎么把这部分考虑进去？答案是把“是否具备快速无功响应能力”当成设备选型或DG类型的属性，而不是把电磁暂态过程原样搬到优化配置里。换句话说，我们在规划阶段的目标是：配置的光伏DG必须有一个足够大的动态无功可行域，这个可行域在系统发生典型电压跌落场景下仍能给出所需的支撑电流/无功量，大于阈值才算具备合格的无功支撑能力。

这种“规划-校验”的思想处理起来很顺手：先用稳态优化模型确定位置和容量，再用各节点的故障电压跌落程度反推需要的最小动态无功备用，如果配置结果不满足就迭代调整。后来我把校验约束直接塞进优化问题，省掉了迭代流程，见下一节。

2.3 把动态特性折算为稳态规划约束的三种思路

要把毫秒级的快速无功特性纳入分钟级的稳态规划模型，业内常见三种折算方式：

*方式一：附加无功备用约束。