如何使MATPOWER的runpf函数强制执行支路潮流限值？

在Matlab环境中，我已准备好MATPOWER案例文件格式的电网模型，需验证其AC可行性。为此我使用runpf函数执行计算：

define_constants;
mpc = loadcase(case_file_path);
results = runpf(mpc);

多数案例显示可行，但我发现被标记为可行的案例中仍存在支路潮流越限（从物理角度看不合理）。我通过以下代码计算越限情况：

pfs = results.branch(:, PF);
qfs = results.branch(:, QF);
pts = results.branch(:, PT);
qts = results.branch(:, QT);
Sfs = sqrt(pfs.^2 + qfs.^2);
Sts = sqrt(pts.^2 + qts.^2);
active = mpc_target.branch(:, 11) > 0.5;
rates = mpc_target.branch(:, 6); % the 6th column gives the rateA of the branch 
violations_n = sum((Sfs .* active > rates) | (Sts .* active > rates));

我猜测runpf未强制执行潮流限值，但这正是我需要的功能。
我已尝试的方法：

• 在runpf说明中发现pf.enforce_q_lims选项，将其设为true后结果无变化；
• 尝试runcpf函数的cpf.enforce_flow_lims选项，即使将基准案例与目标案例设为同一案例，越限情况仍未改善。
  我不确定后续如何推进，想了解是否因MATPOWER的未知限制导致无法实现该功能？

1. 明确runpf的核心定位 #

runpf是交流潮流计算工具，仅求解满足基尔霍夫定律和节点功率平衡的潮流解，不会主动强制执行支路潮流限值——这类限值属于**最优潮流(OPF)**的约束范畴，而非常规潮流的求解目标。pf.enforce_q_lims仅针对发电机无功出力限值，和支路潮流无关，所以设置后无变化是正常的。

2. 用runopf替代runpf实现约束校验 #

要验证"满足支路潮流限值的AC可行性"，应使用MATPOWER的runopf函数，它会在求解潮流的同时将支路潮流作为约束加入优化问题：

define_constants;
mpc = loadcase(case_file_path);
% 配置OPF选项，强制执行支路潮流限值
mpc.opf.options.verbose = 0; % 可选，关闭冗余输出
mpc.opf.enforce_flow_lims = true; 
results_opf = runopf(mpc);

• 若results_opf.success == 1，说明电网存在满足所有约束（包括支路潮流限值）的AC可行解；
• 若success == 0，则当前工况下无法满足支路潮流限值，电网不可行。

3. 纠正runcpf的使用误区 #

runcpf是连续潮流计算工具，用于追踪系统从基准工况到电压崩溃的过程。cpf.enforce_flow_lims的作用是在潮流追踪时，遇到支路潮流越限就停止计算，而非主动调整潮流避免越限，因此不适合用于静态可行性验证。

4. 自定义支路约束的细节 #

如果需要针对特定支路设置限值，可手动修改mpc.branch的对应列：

• 第6列(rateA)：支路长期允许视在功率限值
• 第7列(rateB)：短期允许限值
• 第8列(rateC)：紧急状态允许限值
  runopf默认采用rateA作为约束上限，可根据需求调整对应数值。

内容的提问来源于stack exchange，提问作者wittn