简介这份资源是面向电气工程及其自动化专业学生的毕业设计参考资料聚焦某220kV电网的潮流计算与输电线路继电保护配置适合正在做电力系统方向毕设或需要复习相关计算的读者。压缩包内仅含1个doc文档大小约1.43MB内容按设计流程组织涵盖网络参数计算、PQ分解法潮流计算源程序及结果、对称与不对称短路电流计算程序、线路保护整定与选型并附有系统接线图、零序网络图、潮流分布图和线路保护配置图等图纸。读者可据此掌握变压器与线路阻抗标幺值计算、潮流分布求解、短路电流程序实现以及保护方式选择与整定计算的完整思路对理解220kV电网从参数建模到保护配置的全过程有较高参考价值。目前已有303人学习下载。1. 从一份 220kV 电网毕设拆起潮流计算与继保整定到底怎么落地带过几届电力系统方向的毕设之后我发现一个规律真正卡住人的从来不是“潮流计算是什么”而是拿到一份 220kV 电网的原始参数后不知道从哪一步开始算、算完怎么验证、短路电流出来之后怎么配保护。这份资源就是一份完整的 220kV 电网潮流计算及输电线路继电保护配置设计文档覆盖了从网络参数折算、PQ 分解法潮流计算、对称与不对称短路电流计算到线路保护选型、整定计算和原理图绘制的全链路。它适合电力系统及其自动化专业的本科生做毕设参考也适合刚入行的变电/线路设计人员拿来对照自己的手算流程。整份文档的核心工具链是 C 语言编写的潮流计算程序和短路计算程序配合手工整定计算最终产出潮流分布图、短路电流表和保护配置图。下面我按实际复现的顺序把每一步拆开讲。2. 网络参数折算标幺值、变压器等值电路和线路参数怎么算不出错2.1 标幺值折算的基准选取逻辑整个计算的第一步是把所有元件的有名值折算成统一基准下的标幺值。这份设计选的是 SB 100MVAVB 取各电压等级的平均额定电压。为什么选 100MVA因为电力系统计算中 100MVA 是惯例基准容量方便和设备铭牌参数对照。平均额定电压取 230kV220kV 侧和 115kV110kV 侧而不是标称的 220kV 和 110kV原因是平均额定电压更接近实际运行电压折算误差更小。标幺值折算的核心公式就一条X* X(Ω) × SB / VB²其中 X(Ω) 是元件的实际电抗有名值SB 是基准容量MVAVB 是基准电压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 450kWVS% 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 是 kWVN 是 kVSN 是 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并联单回路直接乘长度。以 L12×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] 存电导 Gyb[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; } }这段代码的关键在于区分三种支路类型接地支路i0、纯阻抗支路kb0、变压器支路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。迭代步骤给定各节点电压初值PQ 节点 V1.0∠0°PV 节点 V 给定、δ0°平衡节点 V 和 δ 给定计算各 PQ 节点的有功不平衡量 ΔPi 和无功不平衡量 ΔQi判断收敛性如果所有 |ΔP| 和 |ΔQ| 都小于收敛精度 ε通常取 10⁻⁵停止迭代解修正方程 ΔP/V B × Δδ得到相角修正量 Δδ解修正方程 ΔQ/V B × ΔV得到电压幅值修正量 ΔV更新 δ 和 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 i1; 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、9PV 节点若干平衡节点 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 100MVAVB 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.0s5.2 距离保护整定计算距离保护 I 段的整定阻抗Z_set_I K_rel × Z_line其中 K_rel 是可靠系数一般取 0.8~0.85Z_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.33I_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 × IfK_sh 一般取 1.8。继电保护整定中电流速断保护需要躲过冲击电流所以要用 i_sh 而不是 If。6.5 距离保护整定未考虑互感器变比现象保护定值与预期偏差很大或者保护范围明显不对。原因整定阻抗是一次侧值但保护装置接收的是二次侧电压和电流需要按互感器变比折算。解决二次侧整定阻抗 Z_set_sec Z_set_prim × n_TA / n_TV。其中 n_TA 是电流互感器变比如 1200/5 240n_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 程序输出逐行检查快得多。从那以后我每次做完潮流计算都强制走一遍这个验证流程——先看电压范围再看平衡节点功率最后抽查几条线路的功率方向。希望帮到你。本文还有配套的精品资源点击获取
