企业官网建设流程全解析

简介：《电力系统分析：CHAPTER 11稳定性》是电力工程专业核心课件，面向电气工程学生、电网工程师及备考人员，系统讲解电力系统稳定性的基本原理与分析方法。内容沿着“稳定性概述→摆动方程→同步电机稳态模型→静态稳定与暂态稳定→等面积准则→三相短路影响→非线性方程数值解法→多机系统暂态分析”逐层展开，既涵盖理论推导，也结合故障案例，使读者能完整理解单机到多机系统的稳定性分析框架。资源为单个PPT文件，大小约11.46MB，已有64人学习。课件详细演示了转子运动方程、功角概念，以及用欧拉法、龙格-库塔法求解摆动方程的具体步骤，并针对三相短路等严重扰动给出等面积准则判定方法与暂态稳定判据，对理解继电保护配合、断路器开断时间选择以及电网安全运行均有直接帮助。此外，多机系统与多机暂态分析的讲解，也适合作为电气工程高年级本科生与研究生的教学参考。

1. 电力系统稳定性分析到底在分析什么

电力系统稳定性问题并不是一个纯理论概念，它直接与电网调度、继电保护定值、机组励磁参数甚至新能源场站的控制策略挂钩。第 11 章这类内容通常出现在《电力系统分析》教材的后半部分，意味着读者已经具备潮流计算、短路计算、元件数学模型的基础，接下来要面对的是系统受到扰动之后能不能继续同步运行的问题。简单说，稳定性分析回答三件事：扰动后系统是否还能回到稳定运行点、如果能需要多长时间、边界条件在哪里。对从事电网规划、保护整定、发电厂电气设计或者新能源并网研究的工程师来说，这一章是把前面学到的稳态分析推向动态分析的关键节点，也是后续学习电力系统仿真软件时的理论基础。

2. 稳定性的分类与数学模型：静态、暂态、动态的区别

2.1 先看数学模型：从单机无穷大系统切入

绝大多数稳定性教材在讲透理论时会从单机无穷大系统（Single Machine Infinite Bus，SMIB）开始，因为这个模型把问题压缩到最核心的物理关系上：一台发电机通过输电线路连接到一个电压和频率恒定的系统。此时发电机的电磁功率可以写成：

[ P_e = \frac{E'V}{X_{\Sigma}} \sin\delta ]

其中 (E') 是发电机暂态电抗后的电动势，(V) 是无穷大母线电压，(X_{\Sigma}) 是暂态电抗与线路电抗之和，(\delta) 是发电机功角。这个公式是整个稳定性分析的出发点，后续所有关于功角稳定、暂态稳定、小信号稳定的讨论都从这里延伸。工程上做等值计算时，经常把多机系统等值为两机或单机无穷大系统来快速评估，虽然是近似，但能抓住主要矛盾。

转子运动方程是另一个核心表达式，通常写成：

[ \frac{2H}{\omega_0} \frac{d^2\delta}{dt^2} = P_m - P_e ]

(H) 是惯性时间常数，(\omega_0) 是同步角速度，(P_m) 是机械功率。这个方程说明了功角变化的根本原因：机械功率与电磁功率不平衡时，转子转速就会变化，功角随之改变。编写仿真程序时，通常把二阶方程改写为一阶微分方程组，用状态空间形式表示，方便数值积分。

2.1.1 为什么用暂态电抗而不是同步电抗

这里有一个容易混淆的点：手工计算短路电流时用 (X_d'') 或 (X_d')，但稳定性分析中转子运动方程里的电磁功率表达式用的是暂态电抗 (X_d') 后的电动势。原因在于稳定性分析的时间尺度是秒级，励磁绕组产生的磁链在这个时间尺度内基本保持不变，用暂态电抗可以避开励磁绕组的电磁暂态细节，把注意力集中在转子机械运动上。如果分析的是次同步振荡或电磁暂态，那就要用更详细的 Park 方程模型，这是不同章节的分工。

2.2 稳定性的分类界限

工程上通常把电力系统稳定性分成三类：静态稳定、暂态稳定、动态稳定（小信号稳定）。静态稳定研究系统在平衡点附近受微小扰动后的恢复能力，暂态稳定研究系统在遭受大扰动后能否保持同步，动态稳定则关注系统在控制器作用下的振荡收敛特性。实际电网中这三者不是独立的，大扰动后系统先经历暂态过程，随后进入动态过程，最终回到新的稳态。

2.2.1 静态稳定的判定条件

静态稳定的判据在功率特性曲线上看得很直观：当 (\delta < 90^\circ) 时，(\frac{dP_e}{d\delta} > 0)，系统静态稳定；当 (\delta > 90^\circ) 时，斜率变负，系统无法维持稳定运行。也就是说，静态稳定极限功率对应的功角是 (90^\circ)，对应的功率极限是：

[ P_{max} = \frac{E'V}{X_{\Sigma}} ]

实际运行中需要留裕度，电力系统安全稳定导则规定正常运行方式下功角一般不允许超过这个极限，这是调度部门安排送端出力时的理论依据。这里要注意，“静态稳定”和“电压稳定”容易混淆，静态稳定特指功角稳定问题，电压稳定讨论的是负荷节点的电压可维持性，两者判据完全不同。

表 2-1 整理了三种稳定性的对比，方便后续讨论时对照。

2.3 暂态稳定分析中常用的简化假设

做暂态稳定分析时，为了在可接受的计算量下获得可靠结果，通常会做几个假设：忽略定子电磁暂态，只用代数方程描述定子磁链；用暂态电抗后的电动势表示发电机；负荷简化成恒定阻抗。这些假设在传统同步机系统中误差不大，但面对双馈风机和光伏逆变器时明显不够用，新能源机组的控制响应时间在毫秒级，必须用详细的电磁暂态模型或动态相量模型。因此现在做电网稳定性分析的工程师往往要同时掌握机电暂态和电磁暂态两种工具，遇到新能源比例高的系统时不能只依赖教材里的经典模型。

3. 暂态稳定计算的落地路径：等面积法则与数值积分

3.1 等面积法则：不写代码也能判断稳定性

等面积法则（Equal Area Criterion）适用于单机无穷大系统，它的核心思想是：当系统受到扰动时，加速面积与减速面积的相对大小决定了机组能否回到稳定运行点。假设输电线路发生单相接地短路，故障期间电磁功率下降，发电机输入机械功率大于电磁功率，转子加速；故障切除后线路阻抗变化，电磁功率上升，转子开始减速。如果故障切除时的功角使得减速面积大于或等于加速面积，系统稳定；否则失步。

手算等面积法则的步骤非常固定：

1. 计算正常运行时的功率特性曲线 (P_{e1}(\delta))
2. 计算故障期间的功率特性曲线 (P_{e2}(\delta))
3. 计算故障切除后的功率特性曲线 (P_{e3}(\delta))
4. 根据故障切除角 (\delta_c) 比较加速面积 (A_1) 和减速面积 (A_2)

[ A_1 = \int_{\delta_0}^{\delta_c} (P_m - P_{e2}) d\delta ]

[ A_2 = \int_{\delta_c}^{\delta_{max}} (P_{e3} - P_m) d\delta ]

如果 (A_2 \geq A_1)，暂态稳定；否则失稳。这个方法虽然简单，但能帮人建立“面积”这个直观概念，后面看仿真曲线上功角来回摆动的现象时更容易理解原因。实际电网分析中不直接用它算结果，而是拿它来估算临界切除时间的大致范围，再结合时域仿真精确定位。

3.2 用 Python 求临界切除时间的数值实现

工程上精确计算临界切除时间（Critical Clearing Time，CCT）需要用数值积分法逐步求解转子运动方程。最常见的做法是用四阶龙格-库塔法（RK4）对功角和角速度进行时域积分，直到系统失稳或恢复到稳定范围。下面给出一个最小可运行的 Python 示例，针对单机无穷大系统在短路故障情况下的摇摆曲线计算。

import numpy as np # 系统参数 H = 5.0 # 惯性时间常数，单位秒 X1 = 0.3 # 正常运行总电抗（标幺值） X2 = 0.5 # 故障期间总电抗（线路部分被短路旁路后变大） X3 = 0.4 # 故障切除后总电抗（单回线运行） E = 1.0 # 发电机暂态电动势 V = 1.0 # 无穷大母线电压 Pm = 0.8 # 机械功率 # 不同运行方式下的电磁功率函数 def power_angle(delta, X): return E * V / X * np.sin(delta) # 故障期间电磁功率 def pe_fault(delta): return power_angle(delta, X2) # 故障切除后电磁功率 def pe_post(delta): return power_angle(delta, X3) # 转子运动方程的一阶状态空间形式 def rotor_eq(state, t, X_mode): delta, omega = state if X_mode == 'fault': pe = pe_fault(delta) else: pe = pe_post(delta) ddelta = omega domega = (Pm - pe) * (np.pi * 50 / H) # 50Hz系统，标幺值转矩转换为角加速度 return np.array([ddelta, domega]) # 四阶龙格-库塔积分 def rk4_step(f, state, t, dt, X_mode): k1 = f(state, t, X_mode) k2 = f(state + 0.5 * dt * k1, t + 0.5 * dt, X_mode) k3 = f(state + 0.5 * dt * k2, t + 0.5 * dt, X_mode) k4 = f(state + dt * k3, t + dt, X_mode) return state + (dt / 6.0) * (k1 + 2*k2 + 2*k3 + k4) # 计算某一故障切除时间下的功角响应 def simulate(tc, t_total=3.0, dt=0.001): t_clear = tc steps_total = int(t_total / dt) steps_clear = int(t_clear / dt) state = np.array([0.5, 0.0]) # 初始功角约28.6度，初始角速度为零 delta_hist = [] for i in range(steps_total): if i < steps_clear: mode = 'fault' else: mode = 'post' state = rk4_step(rotor_eq, state, i * dt, dt, mode) delta_hist.append(state[0]) # 功角超过180度视为失稳，提前终止 if state[0] > np.pi or state[0] < -np.pi: break return delta_hist # 用二分法逼近临界切除时间 def find_cct(low, high, precision=0.001): while high - low > precision: mid = (low + high) / 2 hist = simulate(mid) if max(hist) > np.pi: high = mid # 失稳，降低切除时间 else: low = mid # 稳定，可尝试更大切除时间 return (low + high) / 2 cct = find_cct(0.05, 1.0) print(f"临界切除时间约为: {cct:.4f} 秒")

代码中rk4_step函数实现了标准四阶龙格-库塔单步推进，rotor_eq返回状态向量的一阶导数。simulate函数在故障切除时刻切换电磁功率表达式，模拟故障和切除后的两段过程。找到临界切除时间用的是二分法：给定区间内，如果模拟结果功角超过 180 度（(\pi)）则判断失稳，说明切除时间过长；否则认为在当前切除时间下系统能维持同步。

这里有一个需要注意的数值细节：步长dt的选择直接影响结果精度。步长太大（比如 0.01 秒）时，功角接近 180 度附近的积分误差会放大，可能导致误判；步长太小则计算量增加。实际工程仿真软件（如 BPA、PSASP）内部采用变步长积分器，输出步长与积分步长分离，读者手工编程时固定步长建议不大于 0.005 秒。

3.3 从摇摆曲线判断暂态稳定与失稳

运行上面的代码后，可以分别用tc = 0.1和tc = 0.3调用simulate，绘制功角随时间的变化曲线。稳定情况下，功角曲线在故障切除后先上升到一个最大值，然后回落并围绕某个新的平衡点振荡，振幅逐渐衰减（如果模型中考虑了阻尼）；失稳情况下，功角单调上升穿越 180 度，机组与系统失去同步。

实际电网的暂态稳定分析不止看功角一个指标，还会看母线电压、线路功率、机组转速等物理量，但功角是判断同步性的最直接依据。电力系统安全稳定导则对暂态稳定的判据描述是：在同一时间内，系统中任意两台机组之间的功角差不超过某个限值（通常 180 度或 360 度，按系统规模确定）。这句话翻译到仿真分析里，就是计算所有机组相对参考机组的最大功角差，观察其是否在仿真时段内先增后减。

4. 小信号稳定性分析：状态空间与特征值

4.1 从非线性方程到线性化状态空间

暂态稳定分析关心的是“大扰动后系统能不能撑住”，小信号稳定性分析关心的是“系统在运行点附近受到小扰动后是收敛还是发散”。两者的数学工具不同：暂态稳定用时域积分，小信号稳定用李雅普诺夫线性化方法。把转子运动方程在某运行点附近线性化，得到：

[ \Delta\dot{x} = A \Delta x ]

其中矩阵 (A) 是状态方程的雅可比矩阵。系统的稳定性由 (A) 的特征值决定：所有特征值实部均为负时，系统在该运行点渐近稳定；只要有特征值实部为正，系统不稳定；实部为零时处于稳定边界。对于 50Hz 电力系统，特征值通常被转换成模式的形式来描述振荡频率和阻尼比：

[ \lambda = \sigma + j\omega ]

阻尼比的计算公式是：

[ \zeta = \frac{-\sigma}{\sqrt{\sigma^2 + \omega^2}} ]

工程上对阻尼比有明确要求：电力系统稳定器（PSS）投运后，系统主要振荡模式的阻尼比一般应大于 0.03 到 0.05，低于这个值就要调整 PSS 参数或增加额外的阻尼控制。这个值比特征值本身更容易被工程人员理解，所以小信号稳定分析报告里通常同时给出特征值、振荡频率和阻尼比。

4.1.1 单机系统的小信号模型搭建

单机无穷大系统的小信号分析只需要两个状态变量：功角偏差和角速度偏差。线性化后的状态矩阵为：

[ A = \begin{bmatrix} 0 & 1 \ -\frac{K_S}{M} & -\frac{D}{M} \end{bmatrix} ]

其中 (K_S = \frac{dP_e}{d\delta}) 是同步转矩系数，(M = \frac{2H}{\omega_0}) 是惯性常数，(D) 是阻尼转矩系数。用 Python 可以很方便地计算这个矩阵的特征值：

import numpy as np H = 5.0 D = 1.5 omega0 = 2 * np.pi * 50 M = 2 * H / omega0 # 在功角约70度处的同步转矩系数 delta_s = 70 * np.pi / 180 E = 1.0 V = 1.0 X = 0.4 KS = E * V / X * np.cos(delta_s) A = np.array([[0, 1], [-KS / M, -D / M]]) eigvals = np.linalg.eigvals(A) for lam in eigvals: sigma = lam.real omega = lam.imag zeta = -sigma / np.sqrt(sigma**2 + omega**2) freq = omega / (2 * np.pi) print(f"实部={sigma:.4f}, 虚部={omega:.4f}, 频率={freq:.2f}Hz, 阻尼比={zeta:.4f}")

这段代码直接构造了状态矩阵并计算特征值。注意KS的计算依赖于当前运行点的功角，运行点不同，同步转矩系数不同，特征值也不一样。这也是小信号分析必须针对具体运行方式逐个校核的原因。实际系统中要针对多种典型运行方式分别计算特征值，找出阻尼比最小的模式，再判断是否需要投运或调整 PSS。

4.2 特征值分析结果怎么看

表 4-1 列出了一个典型单机系统在不同运行点下的特征值变化趋势。功角从 30 度增加到 80 度时，同步转矩系数 (K_S) 逐渐下降，特征值实部向正方向移动，阻尼比随之降低。

阻尼比随功角增大的下降趋势说明：系统重负荷时小信号稳定性变差，这与工程经验一致。这就是为什么调度运行中要限制重要断面的传输功率，并且在线监测系统中会实时计算关键模式的阻尼比。当阻尼比低于阈值时，系统会发出告警，提示运行人员调整出力或采取其他措施。

小信号分析相对时域仿真有一个明显优势：特征值直接给出稳定裕度和振荡特性的量化指标，而且计算结果不依赖仿真时长。缺点是特征值分析基于线性化模型，当系统中存在大量电力电子装置时，控制器的开关特性和限幅环节会使线性化失去精度，这时需要结合时域仿真和阻抗分析法做交叉验证。

5. 工程中校验稳定性结果的具体技巧

5.1 判断仿真结果是否可信的三个核心检查

做完稳定性仿真，第一步不是看功角曲线漂不漂亮，而是先校验模型的初始条件是否正确。潮流计算得到的发电机有功出力、机端电压、功角初值必须输入到动态仿真程序中作为初始状态，否则转子运动方程在 (t=0) 时刻就有不平衡转矩，仿真的第一条曲线就会偏离正确轨迹。检查方法是看仿真前 0.1 秒内各状态变量是否基本保持不变，如果有明显漂移，说明初值没有对齐。

第二步是检查故障切除时间与保护动作时间是否匹配。实际系统中线路主保护动作时间约 20 到 40 毫秒，后备保护最长可达数百毫秒。仿真时设置的故障切除时间如果短于实际保护动作时间，计算结果会偏乐观，得到的稳定结论没有工程意义。一般做法是先按主保护时间做一种工况，再按主保护拒动、后备保护动作的极端工况做另一种，两者都稳定才算满足要求。

第三步是看曲线末端是否真正进入稳态。有些工况下功角曲线在前 2 秒内看起来收敛，但之后又缓慢发散，这是因为模型里包含了调速器和励磁系统的慢动态环节。判断稳定性的时间窗口要覆盖这些慢动态的周期，不能只观察故障后 1 秒。GB 38755 中要求暂态稳定仿真的时间长度要达到能够正确反映系统动态特性的程度，工程上通常取 10 到 20 秒。

5.2 提高暂态稳定性的参数调整方向

从工程角度看，提高稳定性的思路比公式更实用。快速切除故障是最有效的手段，缩短保护动作时间和断路器开断时间直接减小加速面积。发电机的快速励磁系统可以提高同步转矩，但过强的励磁调节可能带来负阻尼，所以现代大型机组普遍配置电力系统稳定器（PSS）来弥补励磁系统对低频振荡的不利影响。这是“暂态稳定”和“小信号稳定”在工程措施上交叉的典型场景。

另一个常用措施是电气制动和快关汽门。在送端电网发生故障时，通过投入制动电阻消耗过剩的电磁功率，或者快速关闭汽轮机调节汽门减少机械功率输入，都能有效减少加速面积。新能源场站参与稳定控制则更多依赖于逆变器的暂态控制和频率响应功能，这是新一代电力系统稳定性分析的热点方向。

5.3 用输出结果反向检查模型参数的技巧

最后分享一个实用技巧：把仿真得到的功角振荡频率和现场实测数据做对比。对于一个双机系统，振荡频率约等于 (\sqrt{K_S / M})，因此可以粗略用摇摆曲线的主频率反推系统的等值惯量或等值电抗。如果仿真频率与实测频率偏差超过 20%，大概率是模型参数有问题，比如惯性时间常数设置错误、线路电抗取值与实测不符、励磁系统模型的时间常数不合理。用这个思路可以在没有完整实测数据的情况下快速定位模型中的明显错误。

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