企业官网建设流程全解析

简介：这份资源是面向电气工程及其自动化专业学生的毕业设计参考资料，聚焦某220kV电网的潮流计算与输电线路继电保护配置，适合正在做电力系统方向毕设或需要复习相关计算的读者。压缩包内仅含1个doc文档，大小约1.43MB，内容按设计流程组织，涵盖网络参数计算、PQ分解法潮流计算源程序及结果、对称与不对称短路电流计算程序、线路保护整定与选型，并附有系统接线图、零序网络图、潮流分布图和线路保护配置图等图纸。读者可据此掌握变压器与线路阻抗标幺值计算、潮流分布求解、短路电流程序实现以及保护方式选择与整定计算的完整思路，对理解220kV电网从参数建模到保护配置的全过程有较高参考价值。目前已有303人学习下载。

1. 从一份 220kV 电网毕设拆起：潮流计算与继保整定到底怎么落地

带过几届电力系统方向的毕设之后，我发现一个规律：真正卡住人的从来不是“潮流计算是什么”，而是拿到一份 220kV 电网的原始参数后，不知道从哪一步开始算、算完怎么验证、短路电流出来之后怎么配保护。这份资源就是一份完整的 220kV 电网潮流计算及输电线路继电保护配置设计文档，覆盖了从网络参数折算、PQ 分解法潮流计算、对称与不对称短路电流计算，到线路保护选型、整定计算和原理图绘制的全链路。它适合电力系统及其自动化专业的本科生做毕设参考，也适合刚入行的变电/线路设计人员拿来对照自己的手算流程。整份文档的核心工具链是 C 语言编写的潮流计算程序和短路计算程序，配合手工整定计算，最终产出潮流分布图、短路电流表和保护配置图。下面我按实际复现的顺序，把每一步拆开讲。

2. 网络参数折算：标幺值、变压器等值电路和线路参数怎么算不出错

2.1 标幺值折算的基准选取逻辑

整个计算的第一步是把所有元件的有名值折算成统一基准下的标幺值。这份设计选的是 SB = 100MVA，VB 取各电压等级的平均额定电压。为什么选 100MVA？因为电力系统计算中 100MVA 是惯例基准容量，方便和设备铭牌参数对照。平均额定电压取 230kV（220kV 侧）和 115kV（110kV 侧），而不是标称的 220kV 和 110kV，原因是平均额定电压更接近实际运行电压，折算误差更小。

标幺值折算的核心公式就一条：

X* = X(Ω) × SB / VB²

其中 X(Ω) 是元件的实际电抗有名值，SB 是基准容量（MVA），VB 是基准电压（kV）。注意这里 VB 要取元件所在电压等级的平均额定电压，不是标称电压。很多新手在这里翻车——220kV 侧用 220 而不是 230 去算，结果标幺值偏大 9% 左右，后面潮流迭代收敛变慢甚至不收敛。

2.2 双绕组变压器参数计算

双绕组变压器的等值电路需要四个参数：电阻 RT、电抗 XT、电导 GT、电纳 BT，加上变比 K。这些参数全部从铭牌上的短路损耗 ΔPS、短路电压 VS%、空载损耗 ΔP0、空载电流 I0% 算出来。

以文档中 B1~B3（型号 SFP7-150000/220）为例，铭牌参数为 ΔP0 = 140kW，ΔPS = 450kW，VS% = 14%，I0% = 0.8%，接线方式 YN,d11。计算过程如下：

# 双绕组变压器参数计算（以 B1~B3 为例） SB = 100.0 # 基准容量 MVA VB_220 = 230.0 # 220kV 侧基准电压 kV SN = 150.0 # 变压器额定容量 MVA VN = 220.0 # 变压器额定电压 kV dPS = 450.0 # 短路损耗 kW VS_pct = 14.0 # 短路电压百分比 I0_pct = 0.8 # 空载电流百分比 # 电阻 RT（有名值，Ω） RT = dPS * VN**2 * 1e3 / SN**2 / 1e6 # RT = 450 * 220^2 * 1000 / 150000^2 = 0.968 Ω # 电抗 XT（有名值，Ω） XT = VS_pct * VN**2 * 1e3 / 100 / SN / 1e3 # XT = 14 * 220^2 * 1000 / 100 / 150000 = 45.173 Ω # 电纳 BT（有名值，S） BT = I0_pct * SN * 1e-3 / 100 / VN**2 # BT = 0.8 * 150000 * 1e-3 / 100 / 220^2 = 0.0227e-3 S # 折算到标幺值 RT_star = RT * SB / VB_220**2 # 0.00186 XT_star = XT * SB / VB_220**2 # 0.0854 BT_star = BT * VB_220**2 / SB # 0.012

这里有几个容易搞错的地方。第一，RT 公式里的单位换算：ΔPS 是 kW，VN 是 kV，SN 是 kVA 还是 MVA 要统一。文档里用的是 SN 取 kVA 值（150000），所以分母是 SN² 再乘 10³ 做单位对齐。第二，BT 的标幺值折算方向和阻抗相反——阻抗是乘 SB 除 VB²，导纳是乘 VB² 除 SB。第三，变比 K 对于 YN,d11 接线的变压器，如果两侧都是星形或都是三角形，变比等于匝数比；但星-三角接线时变比要乘 √3。文档中 B1~B3 是 YN,d11，在标幺值等值电路中通常取理想变压器变比来处理，非标准变比放在支路参数里。

2.3 三绕组变压器与线路参数

三绕组变压器（B5~B8，型号 SFPSZ7-120000/220）的参数计算比双绕组多一步：需要从三个短路试验数据中分离出各绕组的短路损耗和短路电压。

# 三绕组变压器各绕组短路损耗分离 # 已知 ΔPS(1-2), ΔPS(2-3), ΔPS(3-1) dPS_12 = 425.0 # kW，示例值 dPS_23 = 310.0 dPS_31 = 380.0 dPS1 = 0.5 * (dPS_12 + dPS_31 - dPS_23) # 247.5 kW dPS2 = 0.5 * (dPS_12 + dPS_23 - dPS_31) # 177.5 kW dPS3 = 0.5 * (dPS_23 + dPS_31 - dPS_12) # 132.5 kW # 各绕组电阻（以高压侧额定电压为基准） VN = 230.0 # kV SN = 120.0 # MVA R1 = dPS1 * VN**2 * 1e3 / SN**2 / 1e6 # Ω R2 = dPS2 * VN**2 * 1e3 / SN**2 / 1e6 R3 = dPS3 * VN**2 * 1e3 / SN**2 / 1e6

短路电压的分离公式和短路损耗完全对称，把 ΔPS 换成 VS% 即可。算出各绕组电抗后，再按双绕组的方式折算标幺值。

线路参数相对简单。文档中线路型号有 LGJ-400、LGJ-185、LGJ-150 等，每公里电阻 r、电抗 x、电纳 b 查表可得。双回路要除以 2（并联），单回路直接乘长度。以 L1（2×LGJ-400/75km）为例：

# L1: 2×LGJ-400, 长度 75km, 双回路 L = 75.0 r_per_km = 0.08 # Ω/km x_per_km = 0.397 # Ω/km b_per_km = 2.92e-6 # S/km R1 = 0.5 * L * r_per_km # 3.0 Ω X1 = 0.5 * L * x_per_km # 14.888 Ω B1 = 2.0 * L * b_per_km # 4.38e-4 S # 折算标幺值（230kV 侧） SB = 100.0 VB = 230.0 R1_star = R1 * SB / VB**2 # 0.00567 X1_star = X1 * SB / VB**2 # 0.02817 B1_star = B1 * VB**2 / SB # 0.2317

注意线路电纳在潮流计算中通常按 π 型等值电路处理，一半挂在首端节点，一半挂在末端节点。所以填入导纳矩阵时，B1_star/2 = 0.1158 分别加到两端节点的自导纳虚部。

提示：线路参数折算时，如果线路跨越两个电压等级（比如 220kV 侧和 110kV 侧通过变压器连接），线路本身的标幺值折算用所在侧的基准电压，不要用对侧的。

3. PQ 分解法潮流计算：从导纳矩阵到 C 程序实现

3.1 导纳矩阵的形成规则

潮流计算的第一步是形成节点导纳矩阵 Ybus。对于 n 个节点的系统，Ybus 是 n×n 的复数矩阵。对角元素 Yii 是节点 i 所有支路导纳之和（包括接地支路），非对角元素 Yij 是节点 i 和 j 之间支路导纳的负值。

文档中的 C 程序用两个实数矩阵分别存储导纳的实部和虚部：yg[i][j] 存电导 G，yb[i][j] 存电纳 B。支路数据输入格式为 I, J, R, X, KB，其中 I 和 J 是支路两端节点号，R 和 X 是支路阻抗（标幺值），KB 是变比或电纳信息。

/* 导纳矩阵形成核心逻辑（从文档 C 程序中提取） */ for (k = 1; k <= ls; k++) { i = z[k][1]; /* 首端节点 */ j = z[k][2]; /* 末端节点 */ r = z[k][3]; /* 电阻标幺值 */ x = z[k][4]; /* 电抗标幺值 */ kb = z[k][5]; /* 变比或电纳 */ if (i == 0) { /* 接地支路：只影响 j 节点的自导纳 */ yg[j][j] += r / (r*r + x*x); yb[j][j] -= x / (r*r + x*x); } else if (kb == 0) { /* 纯阻抗支路（线路） */ g = r / (r*r + x*x); b = -x / (r*r + x*x); yg[i][i] += g; yb[i][i] += b; yg[j][j] += g; yb[j][j] += b; yg[i][j] -= g; yb[i][j] -= b; yg[j][i] -= g; yb[j][i] -= b; } else { /* 变压器支路，kb 为非标准变比 */ /* 等值电路：阻抗在 i 侧，理想变压器在 j 侧 */ g = r / (r*r + x*x); b = -x / (r*r + x*x); yg[i][i] += g; yb[i][i] += b; yg[j][j] += g / (kb*kb); yb[j][j] += b / (kb*kb); yg[i][j] -= g / kb; yb[i][j] -= b / kb; yg[j][i] -= g / kb; yb[j][i] -= b / kb; } }

这段代码的关键在于区分三种支路类型：接地支路（i=0）、纯阻抗支路（kb=0）、变压器支路（kb≠0）。变压器支路的处理最容易出错——非标准变比 kb 放在 j 侧时，j 侧自导纳要除以 kb²，互导纳要除以 kb。如果变比放在 i 侧，公式要相应调整。文档中采用的是“阻抗在 i 侧、理想变压器在 j 侧”的等值电路，这是电力系统计算中的常见约定。

3.2 PQ 分解法的迭代逻辑

PQ 分解法的核心思想是利用高压电网中 R << X 的特点，将有功-相角和无功-电压幅值解耦。修正方程从完整的 Newton-Raphson 方程简化为两个低阶方程：

ΔP/V = B' × Δδ ΔQ/V = B'' × ΔV

其中 B' 和 B'' 是常数矩阵，由导纳矩阵的虚部构成。B' 对应 ΔP-Δδ 关系，B'' 对应 ΔQ-ΔV 关系。文档中的 C 程序用 b1 和 b2 两个矩阵分别存储 B' 和 B''。

迭代步骤：

1. 给定各节点电压初值（PQ 节点 V=1.0∠0°，PV 节点 V 给定、δ=0°，平衡节点 V 和 δ 给定）
2. 计算各 PQ 节点的有功不平衡量 ΔPi 和无功不平衡量 ΔQi
3. 判断收敛性：如果所有 |ΔP| 和 |ΔQ| 都小于收敛精度 ε（通常取 10⁻⁵），停止迭代
4. 解修正方程 ΔP/V = B' × Δδ，得到相角修正量 Δδ
5. 解修正方程 ΔQ/V = B'' × ΔV，得到电压幅值修正量 ΔV
6. 更新 δ 和 V，回到第 2 步

文档中的 C 程序用因子表法（三角分解）求解修正方程。因子表法的好处是 B' 和 B'' 在迭代过程中不变，只需做一次三角分解，后续迭代只需回代，计算量大幅减少。

/* 因子表形成（三角分解）核心逻辑 */ /* 对矩阵 A 进行 LR 分解，A = L × R */ for (i = 1; i <= n; i++) { for (j = i; j <= n; j++) { /* 计算上三角元素 */ sum = a[i][j]; for (k = 1; k < i; k++) { sum -= a[i][k] * a[k][j]; } a[i][j] = sum; } for (j = i+1; j <= n; j++) { /* 计算下三角元素 */ sum = a[j][i]; for (k = 1; k < i; k++) { sum -= a[j][k] * a[k][i]; } a[j][i] = sum / a[i][i]; } }

3.3 节点类型与数据输入

文档中的系统有 10 个节点，其中 PQ 节点 4 个（节点 4、5、8、9），PV 节点若干，平衡节点 1 个。PQ 节点的有功 P 和无功 Q 给定，节点电压 V 和相角 δ 待求；PV 节点的 P 和 V 给定，Q 和 δ 待求；平衡节点的 V 和 δ 给定，P 和 Q 待求。

数据输入时，PQ 节点的负荷功率要折算成标幺值。以节点 4 为例，S4 = 180 + j100 MVA，折算后 S4* = 1.8 + j1.08。注意负荷功率的符号约定：注入节点的功率为正，从节点流出的负荷功率为负。在潮流计算中，PQ 节点的注入功率通常取负值（表示负荷），但文档中的程序采用负荷功率直接输入正值、在计算不平衡量时做减法的处理方式。

注意：PQ 分解法在 R/X 比值较大的配电网中收敛性会变差，但对于 220kV 输电网，R/X 通常小于 0.3，收敛性很好。如果迭代 10 次以上还不收敛，先检查导纳矩阵是否对称、变比处理是否正确、平衡节点是否设置合理。

4. 短路电流计算：对称与不对称故障的程序实现

4.1 对称短路电流计算

对称短路（三相短路）是最严重的短路类型，也是继电保护整定的基础。计算对称短路电流需要用到正序等效网络。文档中的对称短路程序基于节点阻抗矩阵法：先形成正序导纳矩阵，求逆得到阻抗矩阵，然后根据故障节点 f 的自阻抗 Zff 计算短路电流。

/* 对称短路电流计算核心逻辑 */ /* 故障点 f 的三相短路电流（标幺值） */ /* If = 1 / Zff，其中 Zff 是节点 f 的自阻抗 */ /* 首先形成正序导纳矩阵 Y1 */ /* 然后求逆得到 Z1 = Y1^(-1) */ /* 故障点 f 的短路电流标幺值 */ If_star = 1.0 / Z1[f][f]; /* 折算成有名值 */ If_actual = If_star * SB / (sqrt(3) * VB);

文档中对 5、6、7、8 四个节点分别计算了短路电流。以节点 5 为例，如果 Zff = 0.15（标幺值），SB = 100MVA，VB = 230kV，则：

import math Zff = 0.15 SB = 100.0 VB = 230.0 If_star = 1.0 / Zff # 6.667 If_actual = If_star * SB / (math.sqrt(3) * VB) # 1.673 kA

冲击电流 i_sh = √2 × K_sh × If，其中 K_sh 是冲击系数，一般取 1.8（高压电网）。所以 i_sh = 1.414 × 1.8 × 1.673 = 4.26 kA。

4.2 不对称短路电流计算

不对称短路包括单相接地短路、两相短路、两相接地短路。计算不对称短路需要用到对称分量法，将不对称故障分解为正序、负序、零序三个对称网络。

文档中的不对称短路程序需要输入正序、负序、零序三个网络的节点阻抗矩阵。对于单相接地短路，故障点 f 的短路电流为：

If = 3 × Vf(0) / (Z1ff + Z2ff + Z0ff + 3Zf)

其中 Vf(0) 是故障前故障点电压，Z1ff、Z2ff、Z0ff 分别是正序、负序、零序网络中故障点的自阻抗，Zf 是故障点过渡电阻。

# 单相接地短路电流计算（标幺值） Vf0 = 1.0 # 故障前电压，通常取 1.0 Z1ff = 0.15 # 正序自阻抗 Z2ff = 0.15 # 负序自阻抗（通常等于正序） Z0ff = 0.45 # 零序自阻抗 Zf = 0.0 # 金属性短路，过渡电阻为 0 If_star = 3.0 * Vf0 / (Z1ff + Z2ff + Z0ff + 3*Zf) # If_star = 3.0 / 0.75 = 4.0 # 折算成有名值 SB = 100.0 VB = 230.0 If_actual = If_star * SB / (math.sqrt(3) * VB) # 1.004 kA

两相短路和两相接地短路的公式类似，只是分母中的阻抗组合不同。两相短路：If = √3 × Vf0 / (Z1ff + Z2ff)；两相接地短路：If = 3 × Vf0 × (Z2ff + Z0ff) / (Z1ff×Z2ff + Z1ff×Z0ff + Z2ff×Z0ff)。

零序网络的形成是不对称短路计算中最麻烦的部分。零序电流只能通过变压器中性点接地的路径流通，所以零序网络的结构和正序网络差别很大。文档中专门画了零序网络图（附录 b），需要根据变压器的接线方式和中性点接地情况来确定零序通路。

提示：零序阻抗的计算要注意变压器接线方式。YN,d11 接线的变压器，三角形侧零序电流可以流通但不出现在线路上，星形侧中性点接地时零序电流才能流入大地。如果变压器中性点不接地，零序网络中该变压器支路断开。

5. 线路保护整定与选型：从短路电流到保护定值

5.1 保护方式选择

文档中为线路 6-7 选择保护方式。220kV 线路的常见保护配置是：主保护采用纵联差动保护或距离保护，后备保护采用阶段式距离保护或零序电流保护。

线路 6-7 的具体参数需要从潮流分布图和短路电流表中获取。假设线路 6-7 是双回路 LGJ-400 线路，长度约 73km（参考 L6 参数），则：

• 主保护：光纤纵联差动保护，动作时间 0s，保护线路全长
• 后备保护 I 段：距离保护 I 段，整定范围为线路全长的 80%~85%，动作时间 0s
• 后备保护 II 段：距离保护 II 段，整定范围与相邻线路 I 段配合，动作时间 0.5s
• 后备保护 III 段：距离保护 III 段，整定范围与相邻线路 II 段配合，动作时间 1.0s

5.2 距离保护整定计算

距离保护 I 段的整定阻抗：

Z_set_I = K_rel × Z_line

其中 K_rel 是可靠系数，一般取 0.8~0.85；Z_line 是被保护线路的正序阻抗。

# 距离保护 I 段整定 Z_line = 14.491 # 线路 6-7 的正序阻抗（Ω），参考 L6 参数 K_rel = 0.85 Z_set_I = K_rel * Z_line # 12.317 Ω # 折算成二次侧定值 n_TA = 1200 / 5 # 电流互感器变比 n_TV = 220000 / 100 # 电压互感器变比 Z_set_I_sec = Z_set_I * n_TA / n_TV # 二次侧整定值

距离保护 II 段需要与相邻线路的 I 段配合：

Z_set_II = K_rel × (Z_line + K_rel' × Z_next_I)

其中 Z_next_I 是相邻线路 I 段的整定阻抗，K_rel' 是配合系数，一般取 0.8。

5.3 零序电流保护整定

零序电流保护是 220kV 线路接地故障的后备保护。零序 I 段的整定原则是躲过线路末端三相短路时流过保护的最大零序电流：

I_set_0_I = K_rel × 3I_0_max

其中 K_rel 取 1.2~1.3，3I_0_max 是线路末端三相短路时保护安装处的最大零序不平衡电流。

零序 II 段与相邻线路零序 I 段配合，零序 III 段与相邻线路零序 II 段配合。整定计算需要用到短路电流计算的结果，特别是单相接地短路和两相接地短路的零序电流分量。

5.4 保护选型与原理图

保护装置选型需要考虑：额定电压（220kV）、额定电流（根据线路最大负荷电流选择）、保护功能（距离保护、零序保护、纵联差动）、通信接口（光纤纵联需要光接口）。

文档中附录 d 是线路保护配置图，标注了每条线路的保护类型和安装位置。保护原理图（附录中未展开但任务书要求）需要画出电流互感器、电压互感器、保护装置、断路器跳闸线圈之间的电气连接关系。

注意：整定计算中所有短路电流都要折算到保护安装处的互感器二次侧。电流互感器变比和电压互感器变比选错，整定值会差一个数量级，这是血泪经验。

6. 复现这份设计时最容易翻车的五个地方

6.1 标幺值折算基准不统一

现象：潮流计算迭代 20 次以上不收敛，或者收敛后电压幅值明显偏离 1.0。

原因：不同电压等级的元件用了各自的基准电压折算，但导纳矩阵形成时没有按变压器变比归算到同一基准。比如 220kV 侧的线路用 230kV 折算，110kV 侧的线路用 115kV 折算，但变压器支路的变比没有正确处理，导致导纳矩阵不对称。

解决：所有元件的标幺值必须归算到统一基准。变压器的非标准变比要显式体现在导纳矩阵中，不能忽略。检查方法是验证导纳矩阵是否对称：Yij 应该等于 Yji（对于无变压器的纯线路支路）。

6.2 零序网络漏掉中性点接地支路

现象：不对称短路计算中零序电流为零，或者单相接地短路电流等于三相短路电流。

原因：零序网络中漏掉了变压器中性点接地支路，或者把不接地的变压器支路错误地保留在零序网络中。

解决：逐台变压器检查接线方式和中性点接地情况。YN,d11 接线的变压器，星形侧中性点接地时零序通路存在，三角形侧零序电流可以流通但不出现在线路上。中性点不接地的变压器，零序网络中该支路直接断开。

6.3 PQ 分解法迭代发散

现象：前几次迭代不平衡量在减小，但突然增大然后发散。

原因：PV 节点转换为 PQ 节点时无功越限处理不当，或者平衡节点选择不合理。另外，如果系统中存在 R/X 比值较大的线路（比如电缆线路），PQ 分解法的解耦假设不成立，也会导致发散。

解决：检查 PV 节点的无功出力是否在合理范围内。如果无功越限，将该节点转换为 PQ 节点并重新迭代。对于 R/X 较大的线路，改用 Newton-Raphson 法。

6.4 短路电流计算忘记乘冲击系数

现象：短路电流计算结果偏小，继电保护整定值偏低，保护可能误动。

原因：只计算了短路电流的稳态有效值，没有考虑短路瞬间的冲击电流和非周期分量。

解决：冲击电流 i_sh = √2 × K_sh × If，K_sh 一般取 1.8。继电保护整定中，电流速断保护需要躲过冲击电流，所以要用 i_sh 而不是 If。

6.5 距离保护整定未考虑互感器变比

现象：保护定值与预期偏差很大，或者保护范围明显不对。

原因：整定阻抗是一次侧值，但保护装置接收的是二次侧电压和电流，需要按互感器变比折算。

解决：二次侧整定阻抗 Z_set_sec = Z_set_prim × n_TA / n_TV。其中 n_TA 是电流互感器变比（如 1200/5 = 240），n_TV 是电压互感器变比（如 220000/100 = 2200）。折算后 Z_set_sec = Z_set_prim × 240 / 2200 = Z_set_prim × 0.109。

7. 用 Python 做交叉验证：三分钟检查潮流结果是否合理

C 程序跑完之后，怎么快速判断结果对不对？我一般会用 Python 写一个简单的验证脚本，做三件事：检查功率平衡、检查电压范围、检查线路损耗。

import numpy as np # 从 C 程序输出中读取各节点电压和相角 # 假设 10 节点系统，节点 1 为平衡节点 V = np.array([1.05, 1.02, 1.01, 0.98, 0.97, 1.00, 0.99, 0.96, 0.98, 0.97]) delta = np.array([0.0, -2.1, -3.5, -5.2, -6.8, -4.1, -5.5, -7.2, -6.0, -6.5]) # 度 # 检查 1：电压范围是否在合理区间 V_min, V_max = V.min(), V.max() print(f"电压范围: {V_min:.4f} ~ {V_max:.4f} p.u.") # 220kV 电网电压允许范围一般为 0.95 ~ 1.05 p.u. if V_min < 0.95 or V_max > 1.05: print("警告：电压越限，检查无功补偿或变压器变比") else: print("电压范围正常") # 检查 2：平衡节点功率是否合理 # 平衡节点注入功率 = 所有负荷 + 所有线路损耗 # 如果平衡节点功率远大于总负荷，说明线路损耗异常 total_load = 1.8 + 0.72 + 1.2 + 0.36 + 0.4 # 各 PQ 节点有功负荷之和（标幺值） print(f"总负荷: {total_load:.4f} p.u.") # 平衡节点功率应该略大于总负荷（差值即网损） # 如果差值超过总负荷的 10%，检查线路电阻是否偏大 # 检查 3：线路功率方向是否合理 # 功率应该从高电压节点流向低电压节点 # 如果出现功率反向流动，检查相角差是否合理

这个验证脚本的核心逻辑是：电压范围检查能发现无功功率不平衡问题；平衡节点功率检查能发现线路损耗异常；功率方向检查能发现相角计算错误。三分钟跑一遍，比对着 C 程序输出逐行检查快得多。

从那以后我每次做完潮流计算，都强制走一遍这个验证流程——先看电压范围，再看平衡节点功率，最后抽查几条线路的功率方向。希望帮到你。

本文还有配套的精品资源，点击获取