⚡ 电力系统分析 · 课程教学辅助平台Power System Analysis Teaching Site
首页/课前・知源/稳定分析/重点难点解析

第4节 · 重点难点解析

Deep Dive

重点难点深度解析

逐层展开:概念 → 原理 → 推导 → 应用,建议按顺序阅读。

01转子运动方程与惯性时间常数▼

(1)从牛顿第二定律到转矩方程。同步发电机的转子是一个巨大的旋转刚体,其机械角加速度由原动机转矩与电磁转矩之差决定。设转子机械角速度为 ωm(rad/s),转动惯量为 J(kg·m²),则

转子转矩方程(物理形式)J dωm/dt = MT − ME − MD

式中 MT 为原动机机械转矩,ME 为电磁转矩(制动性质),MD 为机械阻尼转矩。稳态时三者平衡,转速恒定,发电机与系统严格同步。

(2)为什么要用功角 δ 而不用转速?因为电网频率被系统整体钳制在 50 Hz 附近,绝对转速的变化量很小(±1 Hz 以内),而转子相对同步旋转参考轴的角位移 δ 才是能量不平衡的直接度量。令同步角速度 ωN = 2πfN,则

功角定义dδ/dt = ω − ωN    d²δ/dt² = dω/dt

(3)化为标幺值并统一到电角度。取转矩基准 MB = SN/ωN(SN 为机组容量基准),把 ωm 换成电角速度 ω = pωm(p 为极对数),并定义惯性时间常数

惯性时间常数TJ = Jω²N/SN  (s)    或  H = TJ/2  (s)

其中 ωN 用机械角速度时 TJ = Jω²mN/SN,用电气角速度时须注意 p² 因子,多数教材统一写成 TJ = 2WK/SN(WK 为额定转速下的动能),这也是 TJ = 2H 的来历。TJ 的工程含义:在额定转矩作用下,把转子从静止加速到额定转速所需的时间,它表征机组抵抗转速变化的能力。国产 300 MW 汽轮发电机组 TJ 约 8~10 s,水轮机组约 5~7 s。

转子运动方程(标幺形式,教材主推式)M d²δ/dt² = PT − PE,  M = TJ/ωN  (s²/rad)
等价写法:d²δ/dt² = (ωN/TJ)(PT − PE)

式中 PT、PE 均为标幺值(以机组容量为基准),δ 以弧度为单位,t 以秒为单位;ωN/TJ 的单位是 rad/s²(当 P 用标幺值时),因此用弧度计算最省事。

易错点:① TJ 与 H 差 2 倍,不同教材、不同仿真软件(PSS/E 用 H,BPA 用 TJ=2H)定义不同,套用公式前必须核对;② 惯量基准必须与功率基准一致,多机系统中 TJ 要折算到系统基准 SB:TJ(系统) = TJ(机组)SN/SB;③ 方程中的 δ 必须用弧度,若用度需乘 π/180;④ 转子运动方程只描述机电暂态,不能用来分析电磁暂态(此时要用派克方程)。
拓展:把 M d²δ/dt² = ΔP 两边乘 dδ 并积分,得 ½Mω² = ∫ΔP dδ,这正是等面积定则的能量解释——左边是动能增量,右边是「不平衡功率—功角」平面上的面积。
02功角特性与静态稳定极限▼

(1)单机—无穷大系统的功角特性。发电机经变压器、线路接至无穷大母线(电压 U、频率恒定),设 Eq 为发电机交轴同步电势,XΣ 为发电机电抗与外部电抗之和(含变压器、线路,不计电阻),则输送到系统的电磁功率为

隐极机功角特性PE = (EqU/XΣ) sinδ = Pmax sinδ

对凸极机,还要加上因交直轴磁阻不等而产生的附加分量:

凸极机功角特性PE = (EqU/XdΣ)sinδ + (U²/2)(1/XqΣ − 1/XdΣ)sin2δ

由于 sin2δ 项的存在,凸极机的极限功率对应角略小于 90°(对水轮发电机约 70°~80°),极限功率比隐极机略大;同时 δ = 0 时仍有少量功率输出(磁阻转矩)。

(2)静态稳定判据。设原动机机械功率 PT 恒定,运行点在 PE(δ0) = PT 处。若受小扰动使 δ 增大 Δδ,则电磁功率增加量 ΔPE = (dPE/dδ)Δδ,多余部分对转子起制动作用,使 δ 回落。因此稳定条件是同步功率系数大于零:

静态稳定判据dPE/dδ > 0   ⇒   SE = Pmax cosδ0 > 0   ⇒   δ0 < 90°

δ < 90° 为稳定域(P-δ 曲线上升段),δ = 90° 是静态稳定极限(临界点,PE = Pmax),δ > 90° 为不稳定域(下降段,扰动后 δ 单调增大直至失步)。

(3)静态稳定储备系数。工程上不允许运行在临界点附近,必须留有余量:

静态稳定储备系数Kp = (Pmax − P0)/P0 × 100%
与功角的关系:若 P0 = Pmax sinδ0,则 Kp = 1/sinδ0 − 1
运行功角 δ0理论 Kp工程评价
30°100%储备充裕,但线路输送能力未充分利用
45°41.4%较理想的经济—安全平衡点
60°15.5%接近正常方式下限 15%~20%
65°10.3%接近事故后方式下限 10%,需重点关注
70°6.4%低于事故后要求,必须增加无功/串补或限制输送
重点提醒:Kp 的分母是当前输送功率而非额定容量,所以同一条线路在轻载时储备系数很高、重载时急剧下降;做断面校核时必须按「正常运行方式」与「事故后运行方式」分别计算并对照导则限值(15%~20% / 10%)。
03等面积定则的几何与物理含义▼

(1)物理本质是能量守恒。把转子运动方程 M d²δ/dt² = PT − PE(δ) 两边同乘 dδ/dt 并对时间积分:

能量形式(M/2)[(dδ/dt)² − (dδ/dt)²|t=0] = ∫δ0δ[PT − PE(δ)] dδ

扰动前稳态时 dδ/dt = 0;若希望扰动后系统能重新达到 dδ/dt = 0 的新(或原)平衡点,则整个过程中右端积分必须为零。把积分按时间段拆开,就得到「加速面积 = 减速面积」。

(2)三段功率特性。以最常见的「一回线首端三相短路、故障后切除该线路」为例,系统有三个不同的功角特性:

三阶段功率特性故障前  PEⅠ = PmaxⅠ sinδ  (正常双回线)
故障中  PEⅡ = PmaxⅡ sinδ  (单回线+短路点,PmaxⅡ 最小)
故障后  PEⅢ = PmaxⅢ sinδ  (单回线运行)

一般有 PmaxⅠ > PmaxⅢ > PmaxⅡ;三相短路且故障线路被切除时 PmaxⅡ 常常小于 PT,此时转子在故障期间持续加速。

(3)两个面积的几何定义。设初始功角 δ0(满足 PmaxⅠ sinδ0 = PT)、故障切除角 δc、切除后不稳定平衡点 δmax(满足 PmaxⅢ sinδmax = PT):

加速面积与减速面积A加 = ∫δ0δc[PT − PmaxⅡsinδ] dδ = PT(δc − δ0) − PmaxⅡ(cosδ0 − cosδc)
A减 = ∫δcδmax[PmaxⅢsinδ − PT] dδ = PmaxⅢ(cosδc − cosδmax) − PT(δmax − δc)
稳定判据:A加 ≤ A减

(4)几何读图要点。在 P-δ 图上,把水平线 P = PT 画出来:曲线在水平线以下的部分(PE < PT)对应的面积是加速面积(面积上边界为水平线、下边界为曲线);曲线在水平线以上的部分(PE > PT)对应的面积是减速面积。增大 PmaxⅢ(如重合闸成功、投入强励、串补)会抬高减速区曲线,使减速面积变大,从而提高稳定极限。

易错点:① 加速面积与减速面积是面积相等而不是功角相等;② PT > PmaxⅡ 时故障中永远加速,此时不能只比较局部面积,必须用 A加(δc) = A减,max 求极限切除角;③ 曲线交点 δmax 必须用故障后的 PmaxⅢ 计算(不是故障中的 PmaxⅡ);④ 阻抗继电器动作、断路器失灵等使切除时间延长,会直接吃掉稳定裕度。
04极限切除角的推导▼

(1)极限切除角的定义。当切除角 δc 增大到某一临界值 δcr 时,减速面积恰好等于加速面积,切除后转子角恰好摇摆到 δmax 而速度为 0(临界稳定)。δc < δcr 稳定,δc > δcr 失稳。

(2)推导。令 A加 = A减 且 δc = δcr:

推导第一步:面积相等PT(δcr − δ0) − PmaxⅡ(cosδ0 − cosδcr)
  = PmaxⅢ(cosδcr − cosδmax) − PT(δmax − δcr)

把含 cosδcr 的项移到一边、含 δcr 的项移到另一边:

推导第二步:合并同类项PT(δcr − δ0) + PT(δmax − δcr) = PmaxⅢcosδcr − PmaxⅢcosδmax + PmaxⅡcosδcr − PmaxⅡcosδ0
⇒ PT(δmax − δ0) = (PmaxⅢ + PmaxⅡ)cosδcr − PmaxⅢcosδmax − PmaxⅡcosδ0
极限切除角公式(最终式)cosδcr = [PT(δmax − δ0) + PmaxⅢcosδmax − PmaxⅡcosδ0] / (PmaxⅢ − PmaxⅡ)

由 cosδcr 反解 δcr(取 90°~180° 之间的解,因为切除角总是大于 δ0,多数情况下大于 90°)。

(3)从切除角到切除时间。故障期间 PE = PmaxⅡsinδ,当 PT 恒定时角加速度近似为常数 α = (ωN/TJ)(PT − PmaxⅡsinδ)。在 δ0 到 δcr 区间取平均加速度,或直接用能量法:

极限切除时间tcr = √(2M · A加,max / (PT − PmaxⅡ))  (当 PT > PmaxⅡ 且近似匀加速时)
其中 A加,max = ∫δ0δcr[PT − PmaxⅡsinδ]dδ(以弧度·标幺计)
使用条件:该式只在 PT > PmaxⅡ(故障期间持续加速、加速度符号不变)时成立,此时 A加,max 与匀加速假设自洽。若 PT < PmaxⅡ,则 δ 会在故障中先加速后减速、摇摆回到交点,应用数值积分求 tcr。
易错点:① cosδcr 反解时误取锐角解;② δmax 用的是 arcsin(PT/PmaxⅢ) 的钝角解(δmax = 180° − arcsin(PT/PmaxⅢ));③ 角度与弧度混用导致 tcr 相差 57.3 倍;④ 忘记把 PT 折算到与 Pmax 相同的容量基准。
05小干扰稳定的阻尼力矩分析▼

(1)线性化模型。设 δ = δ0 + Δδ,ω = ωN + Δω,把转矩方程在平衡点线性化,并考虑励磁绕组、阻尼绕组与励磁调节器引入的附加转矩(分解为与 Δδ 同相的同步分量、与 Δω 同相的阻尼分量):

线性化转子方程M d²Δδ/dt² + D dΔδ/dt + SEΔδ = 0
其中 SE = dPE/dδ|δ0(同步功率系数),D 为等效阻尼系数

对应的特征方程为 Mλ² + Dλ + SE = 0,根为

特征根与振荡频率/阻尼比λ1,2 = −D/(2M) ± j√(SE/M − D²/(4M²))
ωn = √(SE/M) (自然振荡角频率,rad/s)  ζ = D/(2√(M SE))

(2)三类稳定性结论。

  • 静态稳定条件:SE > 0,即 δ0 < 90°,与等面积/功角曲线的结论一致。
  • 振荡稳定条件:D > 0。当 D < 0(负阻尼)时特征根实部为正,系统出现自发增幅的低频振荡(0.1~2.5 Hz),即经典的低频振荡(LFO)/弱阻尼振荡。振荡角频率 f = ωn/(2π)。
  • 振荡频率与惯量关系:M 越大(TJ 越大)ωn 越低;δ0 越大(重载)SE 越小、ωn 也越低,因此重载长距离输电断面更容易出现低频振荡。

(3)负阻尼的来源。高放大倍数快速励磁系统(尤其快速励磁+高增益 AVR)在重载时会产生相位滞后的附加转矩,其与 Δω 同相的分量为负,这是低频振荡最主要的机理性原因。此外,电力电子装置(HVDC、FACTS、风机变流器)的快速控制与外送通道的弱阻尼也会引入负阻尼。

(4)PSS 的作用原理。电力系统稳定器(PSS)取发电机转速偏差 Δω(或功率偏差 ΔP)作为输入,经隔直、相位补偿(超前—滞后网络)、放大与限幅后叠加到励磁调节器的电压给定上,把负阻尼补偿为正阻尼,同时尽量不削弱同步力矩。设计要点是相位补偿必须使 PSS 产生的电磁转矩在 0.1~2.5 Hz 频段内与 Δω 接近同相(滞后角补偿到 0°~20°)。

实例感受:设 TJ = 8 s,fN = 50 Hz,δ0 = 45°,Pmax = 2.0,则 M = 8/314.159 = 0.02546 s²/rad,SE = 2.0 cos45° = 1.4142,ωn = √(1.4142/0.02546) = 7.45 rad/s,f = 1.19 Hz,落在典型低频振荡频段内;若 D = 0.02,则 ζ = 0.02/(2√(0.02546×1.4142)) = 0.053(5.3%),属于弱阻尼,若 D 变为 −0.02 则 ζ = −5.3%,系统将振荡发散。
06暂态稳定的计算与仿真流程▼

(1)数值积分法(时域仿真)。多机系统用经典模型时,每台机组的状态方程为

多机经典模型dδi/dt = ωi − ωN
(TJi/ωN) dωi/dt = PTi − PEi(δ1,…,δn)
PEi = E′i²Gii + ∑j≠iE′iE′j(Gijcosδij + Bijsinδij)

用改进欧拉法(预测—校正)或四阶龙格—库塔法逐步积分,可以得到各机组的摇摆曲线 δi(t);两条曲线之间的最大角差 Δδij 超过 180° 即认为失步。数值积分法是目前工程标准方法,可考虑详细发电机模型、励磁/调速系统、负荷特性与保护动作时序,但对计算量与模型精度依赖大。

(2)直接法(暂态能量函数 TEF)。对经典模型构造能量函数

暂态能量函数V(δ,ω) = ½∑iMiω&subi² + ∫δCOIδ∑i(PEi − PTi)dδi = VK + VP

若故障清除时刻的 V 小于临界能量 Vcr(由相关不稳定平衡点 UEP 确定),则系统稳定。实用化方法有 PEBS(势能边界曲面法)、BCU 法,以及在国内广泛应用的 EEAC(扩展等面积定则):先识别临界机群与剩余机群,把多机系统互补群聚合成一个等值单机—无穷大系统,再用等面积定则严格求解,兼具直接法的速度与积分法的精度。

(3)仿真工具与实践流程。国内常用 PSASP、PSD-BPA,国外常用 PSS/E、DIgSILENT PowerFactory 与 MATLAB/Simulink。一次完整的暂态稳定校核流程为:

  1. 建模与数据准备:潮流数据(含发电机出力、负荷水平)、动态数据(TJ、X′d、励磁与调速模型)、保护与安全自动装置定值。
  2. 确定故障集:按导则要求扫描 N-1、N-2 及严重故障(母线三相短路、线路三相短路并切除、单相重合闸失败等),按「重载断面+关键厂站」优先排序。
  3. 计算条件:明确运行方式(正常/检修/事故后)、负荷模型、故障发生时刻、切除时间(近端 0.08~0.10 s、远端 0.12~0.15 s)、重合闸时间。
  4. 仿真与判据:判定标准通常为「任意两台机组最大功角差 < 180° 且能衰减到稳态」「关键断面功率不越限」「母线电压能恢复到 0.75 p.u. 以上」。
  5. 结论与措施:给出极限切除时间、稳定裕度、薄弱环节与提高稳定的具体措施(切机、快关、串补、加装 PSS 等)。
报告要点:一份合格的稳定分析报告必须写清「故障集 + 计算条件 + 判据 + 结论 + 措施」五要素,只给出一张摇摆曲线截图不足以支撑结论——这是工程评审中最常被退回的地方。