简介这份《电力系统暂态分析期末复习重点》文档面向电气工程及其自动化专业本科生针对期末备考场景帮助读者在有限时间内梳理暂态分析的核心考点与典型题型。内容围绕无限大功率电源三相短路、中性点直接接地系统的单相与三相短路、重合闸对暂态稳定的作用、负序分量与转移阻抗、同步发电机并列运行暂态稳定性、等面积定则、隐极式发电机磁链方程、提高静态稳定性的措施、转移电抗与计算电抗的异同以及对称分量法计算不对称短路电流的完整步骤等展开并配有选择题与简答题便于自测。资源包共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 MVAU_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_eq1.0, X_sigma0.5, S_B100 MVA, U_B10.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_eq1.0, X10.3, X20.3, X00.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 # 例E1.2, U1.0, X_pre0.5, X_fault1.0, X_post0.6, P_m0.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[i1] delta[i] dt/6 * (k1_d 2*k2_d 2*k3_d k4_d) omega[i1] omega[i] dt/6 * (k1_w 2*k2_w 2*k3_w k4_w) return t, delta, omega # 例M6.0, P_m0.8, E1.2, U1.0, X0.6, D0.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(i1) delta(i) dt*ddelta2; omega(i1) 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 更方便但期末考如果要求手算等面积定则还是得练熟。希望帮到你。本文还有配套的精品资源点击获取
