简介：这份《电力系统暂态分析期末复习重点》文档面向电气工程及其自动化专业本科生，针对期末备考场景，帮助读者在有限时间内梳理暂态分析的核心考点与典型题型。内容围绕无限大功率电源三相短路、中性点直接接地系统的单相与三相短路、重合闸对暂态稳定的作用、负序分量与转移阻抗、同步发电机并列运行暂态稳定性、等面积定则、隐极式发电机磁链方程、提高静态稳定性的措施、转移电抗与计算电抗的异同，以及对称分量法计算不对称短路电流的完整步骤等展开，并配有选择题与简答题便于自测。资源包共1个doc文件，约1.25MB，结构紧凑，适合打印或平板批注。目前已有134人学习，可作为考前快速回顾与查漏补缺的参考。
1. 电力系统暂态分析期末复习：为什么背完公式还是算不对
考场上给你一道三相短路题，你记得冲击电流公式，记得派克变换矩阵，甚至记得发电机转子运动方程的每一项系数，但算出来的结果和答案差了一个数量级。这不是你公式背得不够，而是暂态分析这门课考的是“什么时候用哪个模型”的判断力，不是记忆力。电力系统暂态分析的核心就三件事：短路电流怎么算、稳定性怎么判、故障后系统能不能撑住。期末复习的重点不是把课本从头翻一遍，而是把这三条主线上的关键模型、等值电路和计算步骤串起来。适合电气工程本科高年级学生，也适合考研复试要考电力系统分析的考生。下面按短路计算、稳定性分析、避坑排查、进阶技巧四块展开，每一步都落到能动手算的程度。
2. 短路电流计算：从标幺值到冲击电流的完整链路
2.1 标幺值制下的等值电路怎么建
短路计算的第一步不是套公式，是把系统化成标幺值等值电路。很多同学直接跳到冲击电流公式，结果基准值选错，后面全错。常见做法是取基准功率 S_B = 100 MVA 或 1000 MVA，基准电压取各电压等级的平均额定电压 U_av（比如 10.5 kV、37 kV、115 kV、230 kV）。标幺值转换的核心公式是：
# 标幺值转换：已知实际值求标幺值 def to_per_unit(actual_value, base_value): """ actual_value: 实际值（阻抗Ω、电压kV、功率MVA） base_value: 基准值（同量纲） 返回标幺值 """ return actual_value / base_value # 例：一台发电机 X_d'' = 0.2 pu (已标幺)，S_N = 50 MVA，U_N = 10.5 kV # 折算到 S_B = 100 MVA 基准下 S_N = 50 # MVA S_B = 100 # MVA X_d_pp_old = 0.2 # 发电机自身基准下的次暂态电抗 # 电抗折算公式：X_new = X_old * (S_B / S_N) * (U_N / U_B)^2 # 当 U_N = U_B 时，简化为： X_d_pp_new = X_d_pp_old * (S_B / S_N) print(f"折算后次暂态电抗: {X_d_pp_new:.4f} pu") # 输出 0.4000 pu这段代码的逻辑是：发电机铭牌参数是在自身容量基准下的标幺值，必须统一到系统基准才能和其他元件串联。参数说明：S_B 一般选 100 MVA 方便手算，U_B 取平均额定电压而不是变压器额定电压，这是最容易翻车的地方。变压器变比不匹配时，还要做变比修正，否则等值电路里会出现莫名其妙的环流。
建好等值电路后，用网络化简求短路点看进去的等值电抗 X_Σ。手算时用串并联逐步化简，机算时列节点导纳矩阵求逆。期末考一般给 3~5 个元件的系统，手算完全够用。
2.2 三相短路电流周期分量与冲击电流的计算步骤
等值电路建好后，短路电流计算本身不复杂，但步骤不能跳。以三相短路为例，周期分量有效值 I_p 的计算：
import math def three_phase_short_circuit(E_eq, X_sigma, S_B, U_B): """ E_eq: 等值电源电势（标幺值，通常取 1.0） X_sigma: 短路点看进去的等值电抗（标幺值） S_B: 基准功率 MVA U_B: 短路点基准电压 kV 返回：周期分量有效值 kA、冲击电流 kA """ # 周期分量标幺值 I_pu = E_eq / X_sigma # 基准电流 I_B = S_B / (math.sqrt(3) * U_B) # kA # 周期分量实际值 I_p = I_pu * I_B # kA # 冲击电流：K_sh 取 1.8~1.9（发电机母线取 1.9，远离发电机取 1.8） K_sh = 1.8 i_sh = math.sqrt(2) * K_sh * I_p # kA return I_p, i_sh # 例：E_eq=1.0, X_sigma=0.5, S_B=100 MVA, U_B=10.5 kV I_p, i_sh = three_phase_short_circuit(1.0, 0.5, 100, 10.5) print(f"周期分量: {I_p:.2f} kA, 冲击电流: {i_sh:.2f} kA")逻辑说明：先算标幺值下的电流，再乘基准电流还原为实际值。冲击电流系数 K_sh 的选取是考点——发电机机端短路取 1.9，远离发电机取 1.8，差 0.1 在 10 kA 量级上就是几百安的差距。参数说明：E_eq 一般取 1.0，但如果题目给了发电机次暂态电势 E''，要用 E'' 而不是 1.0。X_sigma 是化简后的总电抗，包含发电机、变压器、线路的所有电抗。
2.3 不对称短路：正负零序网络与复合序网
不对称短路是期末必考但最容易丢分的部分。核心思路是对称分量法：把不对称故障分解为正序、负序、零序三个对称网络，在故障点用边界条件连接成复合序网。
正序网络和前面建的等值电路一样，电源电势为正序电势。负序网络没有电源，只有阻抗，发电机负序电抗 X_2 一般取 X_d'' 的 0.8~1.0 倍。零序网络最特殊：零序电流只能通过接地中性点形成回路，所以零序网络的结构和正序完全不同，变压器接线组别（Ynd、YNyn 等）直接决定零序电流能不能流通。
def single_line_to_ground(E_eq, X1_sigma, X2_sigma, X0_sigma): """ 单相接地短路复合序网计算 E_eq: 正序等值电势 X1_sigma: 正序等值电抗 X2_sigma: 负序等值电抗 X0_sigma: 零序等值电抗 返回：正序、负序、零序电流标幺值 """ # 单相接地：三个序网串联 X_total = X1_sigma + X2_sigma + X0_sigma I1 = E_eq / X_total # 正序电流 I2 = I1 # 负序电流等于正序 I0 = I1 # 零序电流等于正序 # 故障相电流 = 3 * I1 I_fault = 3 * I1 return I1, I2, I0, I_fault # 例：E_eq=1.0, X1=0.3, X2=0.3, X0=0.9 I1, I2, I0, I_f = single_line_to_ground(1.0, 0.3, 0.3, 0.9) print(f"正序: {I1:.4f}, 负序: {I2:.4f}, 零序: {I0:.4f}, 故障相: {I_f:.4f} pu")逻辑说明：单相接地时三个序网串联，故障相电流是正序电流的 3 倍。两相短路时正序和负序并联，没有零序。两相接地短路是正序和负序并联后再与零序串联。参数说明：零序电抗 X0 通常比正序大很多，因为零序磁通路径经过空气或油箱壁，磁阻大。如果变压器中性点不接地，零序网络在该变压器处断开，X0 趋于无穷大，单相接地电流为零。
注意：复合序网连接方式取决于故障类型，单相接地是串联，两相短路是并联，两相接地是混联。考前把三种故障的序网连接图默画一遍，比背公式管用。
3. 暂态稳定性分析：等面积定则与摇摆曲线
3.1 等面积定则的物理含义与手算方法
暂态稳定性分析的核心问题是：系统发生故障后，发电机转子会不会失步。等面积定则是最直观的判据——加速面积小于等于减速面积，系统稳定。
功角特性曲线 P_e = (E' * U / X) * sin(δ) 是分析的基础。故障前、故障中、故障后三种状态对应三条不同的功角特性曲线，因为故障期间线路阻抗变化导致传输电抗 X 变化。加速面积是故障期间机械功率 P_m 大于电磁功率 P_e 的部分，减速面积是故障切除后 P_e 大于 P_m 的部分。
import numpy as np from scipy.integrate import quad def equal_area_criterion(E_prime, U, X_pre, X_fault, X_post, P_m, delta_0, delta_clear): """ 等面积定则计算极限切除角 E_prime: 发电机暂态电势 U: 无穷大母线电压 X_pre: 故障前传输电抗 X_fault: 故障中传输电抗 X_post: 故障后传输电抗 P_m: 机械功率 delta_0: 初始功角（弧度） delta_clear: 故障切除角（弧度） 返回：加速面积、减速面积、极限切除角 """ # 三条功角特性曲线的幅值 P_max_pre = E_prime * U / X_pre P_max_fault = E_prime * U / X_fault P_max_post = E_prime * U / X_post # 加速面积：故障期间 (P_m - P_e_fault) 对 delta 积分 def accel_integrand(delta): return P_m - P_max_fault * np.sin(delta) A_accel, _ = quad(accel_integrand, delta_0, delta_clear) # 减速面积：故障切除后 (P_e_post - P_m) 对 delta 积分 # 需要先求最大摇摆角 delta_max，使减速面积等于加速面积 def decel_integrand(delta): return P_max_post * np.sin(delta) - P_m # 用数值方法求 delta_max from scipy.optimize import brentq def area_diff(delta_max): A_decel, _ = quad(decel_integrand, delta_clear, delta_max) return A_decel - A_accel # delta_max 上限取 pi - arcsin(P_m/P_max_post) delta_max_limit = np.pi - np.arcsin(P_m / P_max_post) delta_max = brentq(area_diff, delta_clear, delta_max_limit) A_decel, _ = quad(decel_integrand, delta_clear, delta_max) return A_accel, A_decel, delta_max # 例：E'=1.2, U=1.0, X_pre=0.5, X_fault=1.0, X_post=0.6, P_m=0.8 A_acc, A_dec, d_max = equal_area_criterion(1.2, 1.0, 0.5, 1.0, 0.6, 0.8, 0.5, 1.0) print(f"加速面积: {A_acc:.4f}, 减速面积: {A_dec:.4f}, 最大摇摆角: {np.degrees(d_max):.2f}°")逻辑说明：加速面积和减速面积的积分本质是能量守恒——转子动能增加等于减少。参数说明：X_fault 通常比 X_pre 大，因为故障时线路跳开或阻抗增加，导致 P_max_fault 降低，加速面积增大。极限切除角是使减速面积恰好等于加速面积时的切除角，超过这个角度系统必然失稳。
3.2 摇摆曲线的数值求解与稳定判据
等面积定则只适用于简单系统，多机系统必须数值求解转子运动方程。转子运动方程是二阶微分方程：
dδ/dt = ω - ω_0 M * dω/dt = P_m - P_e - D * (ω - ω_0)
其中 M 是惯性时间常数，D 是阻尼系数。用欧拉法或龙格库塔法求解，得到每台发电机的功角随时间变化曲线，即摇摆曲线。
import numpy as np import matplotlib.pyplot as plt def swing_equation_rk4(M, P_m, E_prime, U, X, D, delta_0, omega_0, t_end, dt): """ 龙格库塔法求解单机无穷大系统摇摆曲线 M: 惯性时间常数 P_m: 机械功率 E_prime: 暂态电势 U: 无穷大母线电压 X: 传输电抗 D: 阻尼系数 delta_0: 初始功角 omega_0: 同步角速度 (rad/s) t_end: 仿真时长 dt: 时间步长 """ t = np.arange(0, t_end, dt) delta = np.zeros(len(t)) omega = np.zeros(len(t)) delta[0] = delta_0 omega[0] = omega_0 def derivatives(delta, omega): P_e = E_prime * U / X * np.sin(delta) ddelta = omega - omega_0 domega = (P_m - P_e - D * (omega - omega_0)) / M return ddelta, domega for i in range(len(t)-1): k1_d, k1_w = derivatives(delta[i], omega[i]) k2_d, k2_w = derivatives(delta[i] + 0.5*dt*k1_d, omega[i] + 0.5*dt*k1_w) k3_d, k3_w = derivatives(delta[i] + 0.5*dt*k2_d, omega[i] + 0.5*dt*k2_w) k4_d, k4_w = derivatives(delta[i] + dt*k3_d, omega[i] + dt*k3_w) delta[i+1] = delta[i] + dt/6 * (k1_d + 2*k2_d + 2*k3_d + k4_d) omega[i+1] = omega[i] + dt/6 * (k1_w + 2*k2_w + 2*k3_w + k4_w) return t, delta, omega # 例：M=6.0, P_m=0.8, E'=1.2, U=1.0, X=0.6, D=0.1 t, delta, omega = swing_equation_rk4(6.0, 0.8, 1.2, 1.0, 0.6, 0.1, 0.5, 314.16, 5.0, 0.001) plt.plot(t, np.degrees(delta)) plt.xlabel('时间 (s)') plt.ylabel('功角 (°)') plt.title('摇摆曲线') plt.grid(True) plt.show()逻辑说明：RK4 比欧拉法精度高，适合刚性不强的转子方程。参数说明：M 越大，功角变化越慢，系统越稳定。D 提供阻尼，D 太小会导致等幅振荡，D 太大则可能过阻尼。稳定判据：功角随时间单调回到稳态值或衰减振荡，则稳定；功角持续增大超过 180°，则失稳。
注意：数值求解时 dt 不能太大，一般取 0.001~0.01 s。dt 太大会导致数值不稳定，误判为系统失稳。
4. 暂态分析期末复习避坑：5 个血泪教训
4.1 标幺值基准混乱导致全盘皆错
现象：短路电流算出来是正确答案的 10 倍或 0.1 倍。原因：不同电压等级没有折算到同一基准，或者变压器变比没有参与标幺值计算。解决：建等值电路前先画一张基准值表，列出每个电压等级的 S_B 和 U_B，所有元件参数统一折算后再化简。变压器用 π 型等值电路时，非标准变比要加理想变压器。
4.2 冲击系数选错导致冲击电流偏差
现象：冲击电流和答案差 5%~10%。原因：K_sh 取值没有根据短路点位置判断。解决：发电机机端短路取 1.9，远离发电机取 1.8，如果题目给了具体参数就按公式 K_sh = 1 + exp(-0.01/T_a) 计算，T_a 是衰减时间常数。
4.3 零序网络结构画错
现象：单相接地电流算出来为零或无穷大。原因：零序网络没有根据变压器接线组别正确断开或连通。解决：画零序网络时逐台变压器检查——YNd 接线零序电流在三角形侧环流，不流出；YNyn 接线零序电流可以流通；中性点不接地则零序网络断开。零序电抗查表时注意区分变压器类型。
4.4 等面积定则中功角特性曲线用错
现象：加速面积和减速面积算反。原因：故障中、故障后、故障前的 X 搞混了。解决：记住故障中 X 最大（线路跳开或阻抗增加），故障后 X 次之，故障前 X 最小。P_max 与 X 成反比，所以故障中 P_max 最小，加速面积最大。
4.5 摇摆曲线数值求解步长过大
现象：仿真结果显示功角振荡发散，但理论上应该稳定。原因：RK4 步长 dt 太大，数值误差累积导致虚假失稳。解决：dt 取 0.001 s，同时检查 M 和 D 的单位是否一致。M 的单位是秒，D 的单位是标幺值/弧度/秒，单位不统一会导致数量级错误。
5. 从期末题到工程判断：用 MATLAB 复现一道综合题
期末复习的终点不是背完公式，而是拿到一道综合题能自己判断用什么模型、算什么量、结果合不合理。下面用 MATLAB 复现一道典型综合题：单机无穷大系统，机端发生三相短路，0.15 s 后切除故障，判断系统暂态稳定性。
% 单机无穷大系统暂态稳定分析 % 参数设置 M = 6.0; % 惯性时间常数 (s) Pm = 0.8; % 机械功率 (pu) E_prime = 1.2; % 暂态电势 (pu) U = 1.0; % 无穷大母线电压 (pu) X_pre = 0.5; % 故障前电抗 (pu) X_fault = 1.0; % 故障中电抗 (pu) X_post = 0.6; % 故障后电抗 (pu) D = 0.1; % 阻尼系数 delta_0 = 0.5; % 初始功角 (rad) omega_0 = 2*pi*50; % 同步角速度 (rad/s) t_clear = 0.15; % 故障切除时间 (s) t_end = 5.0; % 仿真时长 (s) dt = 0.001; % 步长 (s) % 初始化 t = 0:dt:t_end; delta = zeros(size(t)); omega = zeros(size(t)); delta(1) = delta_0; omega(1) = omega_0; % 分段仿真：故障前、故障中、故障后 for i = 1:length(t)-1 if t(i) < t_clear X = X_fault; % 故障中 else X = X_post; % 故障后 end Pe = E_prime * U / X * sin(delta(i)); % RK4 求解 ddelta1 = omega(i) - omega_0; domega1 = (Pm - Pe - D*(omega(i)-omega_0)) / M; delta_mid = delta(i) + 0.5*dt*ddelta1; omega_mid = omega(i) + 0.5*dt*domega1; Pe_mid = E_prime * U / X * sin(delta_mid); ddelta2 = omega_mid - omega_0; domega2 = (Pm - Pe_mid - D*(omega_mid-omega_0)) / M; delta(i+1) = delta(i) + dt*ddelta2; omega(i+1) = omega(i) + dt*domega2; end % 绘图 figure; plot(t, rad2deg(delta), 'LineWidth', 1.5); xlabel('时间 (s)'); ylabel('功角 (°)'); title('摇摆曲线 - 故障切除时间 0.15 s'); grid on; % 判断稳定性 delta_max = max(rad2deg(delta)); if delta_max < 180 fprintf('系统稳定，最大功角 %.2f°\n', delta_max); else fprintf('系统失稳，最大功角 %.2f°\n', delta_max); end这段代码的逻辑是分段仿真：故障期间用 X_fault，故障切除后用 X_post。参数说明：t_clear 是故障切除时间，是稳定性分析的关键变量。如果 t_clear 增大到 0.25 s，加速面积增大，系统可能失稳。实际工程中，继电保护动作时间、断路器分闸时间共同决定 t_clear，一般要求小于极限切除时间。
我自己的习惯是：算完摇摆曲线后，把 t_clear 从 0.1 s 扫到 0.3 s，画出最大功角随切除时间的变化曲线，找到临界切除时间。这比单次仿真更能说明问题，也是期末大题常见的考法。另外，MATLAB 的 ode45 比手写 RK4 更方便，但期末考如果要求手算，等面积定则还是得练熟。
希望帮到你。
本文还有配套的精品资源，点击获取
