基于凸规划形式与启发式方法的N-1输电线路开断下电力系统电压稳定裕度估计

《Modelling》：Estimation of the Voltage Stability Margin in Power Systems Under Transmission Line Contingencies Using a Convex Formulation and a Heuristic Approach

　　

摘要：输电断线N-1故障下的电压稳定性是现代电力系统面临的关键问题，日益增长的电力需求和大规模可再生能源并网使系统运行安全受到挑战。本文研究了N-1输电线路开断条件下的电压稳定裕度(voltage stability margin / loadability factor λ)估计问题，提出三种求解方法：基于内点算法(IPOPT)配合多重初始点策略的非线性规划(NLP)模型、基于逐次牛顿–拉夫森(Newton–Raphson, NR)潮流与渐进负荷缩放的递归启发式方法，以及二阶锥规划(SOCP, Second-Order Cone Programming)凸松弛形式。所提方法在IEEE 9?、14?、30?及57?节点测试系统上验证，涵盖不同拓扑复杂度与冗余度的网络。通过与NLP参考解对比精度，并在全量N-1故障场景下评估计算效率。结果表明，启发式方法精度更高，凸规划形式求解速度显著更快，且两种方法在NLP方法不收敛的情形下均表现出鲁棒性。

论文解读：基于凸规划形式与启发式方法的N-1输电线路开断下电力系统电压稳定裕度估计

一、研究背景与意义

电压稳定性(voltage stability)是指电力系统在遭受扰动后维持所有母线(st bus)可接受稳态电压幅值的能力。随着用电需求增长及非传统可再生能源大规模接入，系统惯量降低且出力波动增大，使电压稳定问题尤为突出，其中未计划输电线路开断(N-1 contingency)是导致电压失稳甚至电压崩溃(voltage collapse)的重要诱因。工程中通常用负荷能力系数(loadability factor, λ)——可同时按比例放大有功和无功负荷的最大倍数——来量化系统在给定约束下距离电压崩溃点的远近。精确求取λ需解决含交流(AC)潮流方程的非线性非凸优化问题，大规模系统下计算困难且易陷入局部最优或不收敛。现有方法如连续潮流(continuation power flow, CPF)、灵敏度分析、启发式算法各有局限。因此，研究人员开展了基于凸松弛(convex relaxation)与递归启发式(recursive heuristic)两种替代方案的N-1故障下电压稳定裕度估计研究，并与内点法NLP参考解对比，以期为在线预想故障筛选与离线规划提供兼顾精度与效率的工具。该文发表于《Modelling》。

二、主要关键技术方法

研究人员采用MATPOWER标准算例(IEEE 9?、14?、30?、57?bus)，分别构建三种模型：①参考非线性规划(NLP)模型——目标函数为max λ，约束为节点有功/无功平衡(含λ均匀缩放负荷)、电压幅值上下限，采用IPOPT求解并施加multi-start策略(5组不同|V_i|初值)以避免局部最优；②递归启发式方法——从基态λ=1开始按固定步长Δλ递增负荷并逐次求解NR交流潮流，记录最后一次收敛的λ为近似最大负荷能力；③二阶锥规划(SOCP)凸松弛模型——引入平方电压幅值W_ii=|V_i|²及两两电压乘积变量W_ij，将原非凸潮流等式松弛为二阶锥约束(‖·‖₂≤W_kk)，目标仍为max λ。排除导致系统孤岛(islanding)的开断线路，NLP与SOCP用Julia+JuMP实现，启发式用MATLAB实现，收敛判据为功率不平衡范数＜10^-8p.u.。

三、研究结果

Mathematical Modeling（数学模型）

研究人员将N-1故障下电压稳定裕度估计建模为最大化λ的NLP问题：目标函数max λ_c（开断线路c时的负荷能力系数），约束包括各节点有功功率平衡P_i(V,θ)=P_G,i?(1+λ)P_L,i^base、无功功率平衡Q_i(V,θ)=Q_G,i?(1+λ)Q_L,i^base、电压幅值V_i,min≤|V_i|≤V_i,max及λ≥0。该AC潮流方程为强非凸，引出后续两种替代解法。

Reference Nonlinear Programming Model（参考非线性规划模型）

研究人员以IPOPT求解上述NLP并配合multi-start(初值|V|=0.95, 1.0, 1.02, 1.04, 1.06 p.u.)取最大可行λ，作为精度参照基准。结果表明，部分复杂N-1场景(如IEEE 57?bus中线路34–35)即便multi-start也无法收敛，说明原始NLP方法存在局限性。

Recursive Heuristic Approach（递归启发式方法）

研究人员采用固定步长Δλ=0.01逐级等比放大P_L,i^base、Q_L,i^base，每步做NR潮流，首个不收敛步的前一级λ记为估计值。步长敏感性分析显示Δλ=0.01在误差(<0.2%)与时间间取得最佳折衷。该方法在全部测试系统中给出最接近NLP参考解的λ值。

Second-Order Cone Approximation（二阶锥规划松弛）

研究人员将AC潮流方程经变量代换W_ii=|V_i|², W_ij=V_iV_j^*并添加SOC约束‖[2Re(W_ij),2Im(W_ij),W_ii?W_jj]‖₂≤W_ii+W_jj形成凸优化问题。径向网络下松弛严格等价；网状输电网络中SOCP给出λ的真值下界(偏乐观估计)。研究人员指出可用恢复相角θ检验松弛紧致性(gap＜5%时可接受)。

Test Systems（测试系统）

选用MATPOWER中IEEE 9?(3机9线)、14?(5台发电机)、30?(6台发电机)、57?bus(7台发电机)系统，排除造成电气孤岛的开断线路(如30?bus中25–26，57?bus中32–33、35–36)。

Numerical Validations and Results Analysis（数值验证与结果分析）

• •
  IEEE 9?bus：最严重开断为线路4–9，λ≈2.38；启发式平均误差0.11%，SOCP在此系统中零误差，均能正确识别最严重故障。
• •
  IEEE 14?bus：最严重开断为线路1–2(两发电机节点间)，λ≈2.75；启发式最大误差0.34%、平均0.08%；SOCP最大误差1.37%、平均0.63%，仍能识别最严重故障。
• •
  IEEE 30?bus：最严重开断为线路6–8，λ≈3.05；更网状拓扑使SOCP误差降至最大2.87%、平均＜1%，启发式＜0.1%。
• •
  IEEE 57?bus：最严重开断为线路25–30(近径向区域)，λ≈1.66；参考NLP对线路34–35不收敛而启发式与SOCP分别给出λ≈1.90与1.89 p.u.；启发式整体最准但耗时约7 min，SOCP仅约0.3 min（加速20余倍）。SOCP在大多数N-1场景下松弛间隙＜5%。

Comparative Analysis of the Results（结果对比分析）

综合四系统，启发式方法平均误差最低(最大不超0.34%)，但计算耗时随规模显著上升；SOCP方法误差通常＜3%，计算效率高出一个数量级以上，适合在线N-1筛选；NLP参考方法存在不收敛风险。两法均能正确识别最小λ对应的关键开断线路。网络冗余度高(更网状)时SOCP松弛更接近真值。

四、讨论与结论翻译

结论：本文对比了三种N-1开断下电压稳定裕度估计方法——带multi-start的内点法NLP、基于逐次NR潮流的递归启发式方法及SOCP凸松弛。SOCP在9?bus系统所有故障下零误差，随系统规模增大出现非零误差(IEEE 14?bus最大1.37%，IEEE 30?bus最大2.87%)，网状拓扑改善松弛精度。启发式方法在所有测试中精度最高(9?bus最大误差0.16%，14?bus 0.34%，30?bus 0.09%，57?bus 0.12%)，但计算代价大——57?bus时耗时约7 min，SOCP仅约0.3 min，提速超20倍。SOCP误差一般＜3%，是在线预想故障评估的实用替代方案；启发式更适合离线规划校核。57?bus系统中参考NLP对线路34–35(34–35)无法收敛，而启发式和SOCP给出一致估计(分别≈1.90 p.u.和≈1.89 p.u.)，证明二者鲁棒性。研究还确认节点数多并不保证更高电压稳定——57?bus系统因局部低冗余与近径向结构具最低基态裕度(λ≈2.66)。最关键开断分别为线路9–4(9?bus)、1–2(14?bus)、6–8(30?bus)、25–30(57?bus)，后者仅允许约66%额外负荷增长即达崩溃点。未来研究方向包括扩展至动态电压稳定分析、考虑间歇性可再生电源随机性、采用紧致凸松弛(如半定规划SDP hierarchy)、启发式步长自适应与并行化、孤岛开断专项分析及凸–非凸混合(warm?start)方法以实现大规模系统实时监测。