简介暂态能量函数法是电力系统暂态稳定分析的重要直接法手段该PPT学习教案系统梳理了方法原理、数学描述与工程应用适合电气工程专业高年级本科生、研究生及从事电网稳定分析的工程技术人员使用。课件以古典力学能量概念切入逐步展开李雅普诺夫函数构造、临界能量判定、李氏定理与引力域等内容并延伸至单机无穷大系统直接法、多机系统特殊问题、两种坐标系、相关不稳定平衡点(RUEP)等专题脉络清晰、层层深入。资源包共1个pptx文件大小约1.23MB、全文共105页公式推导与示意图结合便于课堂授课或自学研读。目前已有161人浏览学习。学习后可掌握暂态能量函数的建模思路与稳定裕度评估方法理解故障轨迹及切除时刻对暂态稳定性的影响机制为电力系统暂态安全分析与控制实践打下扎实基础。1. 暂态能量函数法一套能算出稳定裕度的直接法而不是玄学做暂态稳定分析的人十有八九都有过这样的体会时域仿真曲线拉出来一大把功角摇摆看着是稳了但你要问他“稳定裕度是多少”“再切晚 0.1 秒还稳不稳”他答不上来。这不是仿真做得不够细而是时域法本质上是“答案有了过程没有”——你只知道稳不稳不知道离临界有多远。而暂态能量函数法TEFTransient Energy Function恰好补上这个缺口它不靠逐步积分推时间而是构造一个李雅普诺夫函数把“稳定与否”转化成“扰动结束时的系统能量是否超过临界能量”这一条不等式判断。这份 105 页的 PPT 教案正是把这条从原理到应用的完整链路讲清楚了——从古典力学的能量视角出发到单机无穷大系统的直接法推演再到多机系统的 UEP、RUEP、等位能线这些进阶话题。它适合两类人一是正在学电力系统分析的研究生需要把教材里“直接法”那一章真正吃透二是做在线稳定评估的工程师想理解能量函数类方法到底怎么落地、边界在哪。PPT 里公式推导和物理图景结合得紧不是那种公式堆到底的讲义读起来负担不大。2. 先把物理图景立住能量、临界能量与稳定裕度的关系2.1 从“小球滚碗”理解暂态能量函数的三个关键点这份教案的开篇用了非常经典的力学类比——一个小球在碗壁上运动。无扰动时小球停在碗底也就是稳定平衡点S.E.P受到扰动后小球获得动能和势能开始沿碗壁往上滚如果碗壁足够高、摩擦力存在小球最终会耗散能量回到碗底如果扰动太大小球翻过碗沿就回不来了。这个“碗沿”对应的位置就是不稳定平衡点U.E.P碗沿高度对应的势能就是临界能量 Vcr。由此得到暂态能量函数法的三个关键点找平衡点、构造暂态能量函数 V、确定临界能量 Vcr。这三件事环环相扣任何一个出了问题后面的稳定判断都不成立。在单机无穷大系统里这三个关键点都有明确的解析解稳定平衡点 δs arcsin(Pe/Pm)不稳定平衡点 δu π − δs临界能量 Vcr V(δu)。而在多机系统里这三个问题每一个都会变得棘手——这正是教案里花费大量篇幅讲 AUEP、RUEP 的原因。2.2 为什么坐标变换是构造能量函数的前提条件教案里有一个容易忽视但极其重要的细节在构造暂态能量函数前必须做坐标变换。原系统的稳定平衡点不在坐标原点直接构造 V 函数会导致 V(0) ≠ 0这不满足李雅普诺夫函数的基本条件当 X → 0 时 V → 0。所以要将稳定平衡点 δs 作为新坐标原点定义新状态变量 x δ − δs。这个操作看起来只是平移坐标轴实际上它决定了能量函数能不能在零点处取到最小值是后续一切推导的前提。这个概念在多少教材里被一笔带过但在实际工程中如果你自己写程序构造能量函数时发现 V 在平衡点处不为零第一件事就检查坐标变换做没做对。2.3 临界能量的“故障后网络依赖”属性教案第 16 页强调了另一个关键结论临界能量 Vcr 与系统故障后的网络结构有关也就是网络吸收发电机发出动能的能力。计算时用故障后的运动方式Vcr 也用故障后的参数值。这句话对做过实际计算的人很有启发。初学容易把故障前、故障中、故障后的网络混在一起用但能量函数法的逻辑是故障切除后的系统决定了“能装下多少能量”而故障期间的轨迹决定了“系统实际获得了多少能量”。Vcr 只看故障后网络V(δ) 按故障后的系统方程计算二者必须同源。提示做单机无穷大系统算例时先画出故障前、故障中、故障后三条功角特性曲线再动手列能量表达式。三条曲线对应的 Pe 表达式不同这一步错了后面全错。3. 李氏定理与数学基础为什么能量函数法不是“拍脑袋”3.1 李雅普诺夫稳定性判据的工程解读教案第 78 页给出了李氏定理的完整表述。把数学语言翻译成工程语言对一个自治系统 Ẋ f(X)如果能找到一个标量函数 V(x)满足 V(0) 0、V(x) 0x ≠ 0且 V̇(x) ≤ 0那么系统在原点稳定如果 V̇(x) 0则渐近稳定V̇(x) 0 则不稳定。这里需要注意“稳定”和“渐近稳定”的差别。V̇ ≤ 0 只能说明系统能量不再增长但不保证回到原点——对应物理场景就是小球停在碗壁某个位置不动了如果存在无阻尼的极限环之类的情况V̇ 0 才说明能量持续耗散最终回到平衡点。工程上我们通常希望渐近稳定。暂态稳定分析关注的不是“平衡点附近的小扰动稳定”而是“大扰动后能否回到平衡点或稳定运行点”所以需要用到引力域稳定域的概念——只有初态点落在某个区域内系统才是渐近稳定的这个区域就是引力域。单机无穷大系统的引力域边界恰好由不稳定平衡点 δu 划定。3.2 线性定常系统的李函数构造实例Lyapunov 方程的由来教案第 910 页给出了一个完整的线性系统算例这一步是理解后续非线性推广的根基。对线性系统 Ẋ AX取二次型 V(x) XᵀPX 为候选李函数则V̇ ẊᵀPX XᵀPẊ (AX)ᵀPX XᵀP(AX) Xᵀ(AᵀP PA)X要让 V̇ 负定需要 AᵀP PA 为负定矩阵。工程上通常反过来任选一个正定矩阵 Q常用单位阵令 AᵀP PA −Q解出 P再校验 P 是否正定。这就是 Lyapunov 方程。这个推导的价值在于它给出了一个“可验证”的流程先选 V 的形式再通过求解方程验证 V̇ 的符号最后回到 V 的正定性检查。每一步都有明确的操作对象不靠猜。后续在暂态能量函数法中构造 V 时思路一脉相承——只不过非线性系统的 V 构造没有系统化的方法更加依赖对物理过程的把握。3.3 从线性到非线性能量函数构造为什么难线性系统有 Lyapunov 方程兜底而非线性系统比如电力系统的摆动方程没有通用的构造方法。教案给出的路径是“从物理意义出发”动能部分取 K ½Mω²势能部分通过对功角方程右侧积分得到。这种做法有明确的物理支撑——动能是速度的二次型势能是功角相关项在整个轨迹上的积分。但是多机系统的能量函数构造远没有这么顺畅。不同发电机之间相对运动、网络结构变化、负荷模型差异都会影响 V 是否严格满足李雅普诺夫条件。这也是为什么教案后续花了大量篇幅讲“多机系统的特殊问题”和“两种坐标系”。4. 单机无穷大系统的直接法推演从运动方程到临界切除时间4.1 状态方程建立与坐标原点移动单机无穷大系统的经典二阶摆动方程为M·d²δ/dt² D·dδ/dt Pm − Pe·sin(δ)其中 M 为惯性时间常数D 为阻尼系数Pm 为机械功率Pe·sin(δ) 为电磁功率。令状态变量 x₁ δ − δs功角偏移、x₂ ω角速度偏差则状态方程为dx₁/dt x₂ dx₂/dt (Pm − Pe·sin(δs x₁) − D·x₂) / M注意这里已经完成了坐标平移δs 是稳定平衡点对应的功角。这一步在教案第 18 页交代得很清楚。为什么要这么写因为后续能量函数要在“δs 处取极小值”如果不平移坐标V 在平衡点处不为零整个李雅普诺夫判据就失效了。4.2 动能与势能函数的构造过程对状态方程的第二个式子两边乘以 M·dδ再从 δs 积分到当前功角 δ得到动能和势能的表达式动能来自速度项V_k ½·M·ω²势能来自不平衡功率对功角的积分V_p ∫[δs→δ] (Pe·sin(δ) − Pm) dδ关于这个势能表达式需要多说两句。积分下限取稳定平衡点 δs所以 V_p(δs) 0这是坐标变换的直接收益。积分内部是被积函数 Pe·sin(δ) − Pm 关于功角的积分——它描述了系统偏离平衡点后“存储”的势能。当 δ 在 δs 与 δu 之间时被积函数为正势能随 δ 增大而增加到达 δu 时系统达到临界势能。而临界能量 Vcr V_p(δu)因为 δu 处 ω 0不稳定平衡点处速度为 0动能为 0所以临界能量只包含势能部分。把 δu π − δs 代入势能表达式即可求得 Vcr 的解析值。4.3 暂态稳定判据与临界切除时间的求法得到了系统实际能量 V(δ, ω) 和临界能量 Vcr判断就变成一条不等式V(δ, ω) Vcr → 稳定 V(δ, ω) Vcr → 不稳定 V(δ, ω) Vcr → 临界状态判别稳定裕度时定义裕度 (Vcr − V) / Vcr × 100%在故障切除时刻δ δc、ω ωc计算此时的 V(δc, ωc)与 Vcr 比较就能判断切除后系统是否稳定。临界切除时间 tcr 的求法是在故障持续期间逐步推进时间计算每个时刻的 V(δ, ω)找到第一个使 V(δ, ω) Vcr 的时刻。# 伪代码示例临界切除时间扫描计算流程 import numpy as np # 系统参数标幺值示例 M 8.0 # 惯性时间常数 D 0.1 # 阻尼系数 Pm 0.8 # 机械功率 Pe 1.2 # 电磁功率幅值 # 稳定平衡点与不稳定平衡点 delta_s np.arcsin(Pm / Pe) delta_u np.pi - delta_s # 临界能量: 对 (Pe*sin(delta) - Pm) 从 delta_s 积分到 delta_u from scipy.integrate import quad Vcr, _ quad(lambda d: Pe * np.sin(d) - Pm, delta_s, delta_u) # 故障期间逐步积分, 求 V(delta, omega) 首次超过 Vcr 的时刻 dt 0.001 t 0.0 delta delta_s omega 0.0 while t 1.0: # 故障期间运动方程 (假设恒机械功率, 电磁功率按故障中曲线) Pe_fault 0.2 # 故障中电磁功率幅值 d_delta omega d_omega (Pm - Pe_fault * np.sin(delta) - D * omega) / M delta d_delta * dt omega d_omega * dt # 构造实际能量 (用故障后网络计算势能) V_k 0.5 * M * omega**2 V_p, _ quad(lambda d: Pe * np.sin(d) - Pm, delta_s, delta) V_total V_k V_p if V_total Vcr: print(f临界切除时间约: {t:.4f}s) print(f切除时刻功角: {np.degrees(delta):.2f}°) break t dt这段代码的推导逻辑与教案完全一致故障期间用故障中的功角特性曲线推进轨迹计算能量时必须用故障后的网络参数势能积分中被积函数用 Pe·sin(δ) − Pm而不是故障期间的 Pe_fault·sin(δ) − Pm。这是最容易翻车的区别——轨迹推进看故障中方程能量评估看故障后方程两边不要混。注意上面代码是求解流程示意不是完整工程脚本。实际算例中故障期间 Pe_fault 可能随电气量变化而变化需要结合网络方程联立求解。4.4 等面积定则与能量函数法的等价关系教案第 11 页画的功角特性曲线abcd 包围面积等于 defg 包围面积这正是等面积定则的图像化表述。把等面积定则翻译成能量函数语言加速面积代表故障期间系统获得的动能增量减速面积代表故障后系统能吸收的最大动能。二者相等时对应临界切除状态此时 V(δc, ωc) Vcr。等面积定则和能量函数法本质上是同一个判据的两种表达区别只在视角等面积定则从“面积”看能量函数法从“数值”看。理解了这个等价关系你会发现“V(δ, ω) 与 Vcr 比较”和“加速面积与减速面积比较”是同一件事这对后续理解多机系统中的能量函数拆分解读例如发电机动能项、网络势能项的物理含义非常有帮助。5. 避坑手册能量函数法应用中的五个高频翻车点5.1 翻车点一把 Vcr 当常数忽视了故障后网络依赖现象同一台发电机同一个故障位置换了故障线路阻抗后直接用之前的 Vcr 判断稳定性得出错误结论。原因Vcr 是故障后网络的函数网络拓扑或参数变了δu 变了Vcr 也变。教案第 16 页明确写了“计算时用故障后的运动方式Vcr 也用故障后的值”这句话很多人扫一眼就过去了。解决每次计算 Vcr 前先确认故障后网络的等值参数已更新。在程序中把 Vcr 的计算放在“故障后网络参数已生成”之后而不是放在故障前的初始化阶段。5.2 翻车点二能量积分路径选错——用故障前轨迹评估故障后能量现象算出来的稳度裕度偏大临界切除时间偏高按结果配置保护动作时间实际仿真却失稳。原因V(δ, ω) 的势能部分是从 δs 沿实际轨迹积到当前 δ但实际计算中若错误地用故障前网络的 Pe−δ 曲线去积分得到的是“错误的系统存储能量”该系统根本不存在于故障后的物理过程中。解决严格区分三条 Pe−δ 曲线——故障前、故障中、故障后。实际获取能量用故障后的曲线轨迹推进用故障中的曲线初始平衡点用故障前的曲线求 δs。5.3 翻车点三多机问题硬套单机公式现象把单机无穷大系统的公式直接套到多机系统的某一台发电机上算出“稳定”结果与实际仿真不符。原因多机系统不存在唯一的“功角”发电机之间存在相对运动能量函数的构造需要先选惯性中心COI坐标系把各机功角转换为相对 COI 的角度。教案第 3336 页的多机系统特殊问题部分就是在讲这件事。解决多机系统计算前先求系统的 COI 坐标完成坐标系变换。判断稳定时不能只看单机能量要看系统总能量动能和势能的代数加总与临界能量比较。5.4 翻车点四忽略了阻尼项对能量耗散的影响现象构造能量函数时忽略阻尼算出的稳定裕度偏保守得到的临界切除时间明显小于时域仿真结果。原因阻尼项 D·ω² 在能量表达式中作为耗散项会持续消耗系统动能对稳定有正贡献。教案中的推导先将阻尼设为 0 做简化第 14 页但在对比实际系统时如果不把阻尼效应考虑进去误差会很大。解决工程评估中如果系统阻尼不可忽略建议在能量函数中加入阻尼耗散项做修正。更实用的做法是用能量函数法给初判、用时域仿真给终判两者结合。5.5 翻车点五把 AUEP 当成 RUEP 用现象多机系统分析中有人用“最近不稳定平衡点”AUEP法估算临界能量找出来的点并不是实际故障轨迹最可能穿越的边界点结果或保守或冒进结果不一致。原因AUEP 找的是“距离当前稳定平衡点最近的所有不稳定平衡点中能量最低的那个”它与故障轨迹的具体形状无关而 RUEP相关不稳定平衡点是“故障轨迹实际趋近的不稳定平衡点”它才真正对应了故障后系统的失稳模式。解决优先用 RUEP 法做故障相关性的临界能量评估。工程上常用 BCUBoundary of Stability region via Controlling UEP方法先定位相关不稳定平衡点再计算 Vcr。6. 多机系统进阶两种坐标系、AUEP/RUEP 的选择与验证技巧6.1 COI 坐标系与惯量中心参考多机系统里每台发电机都在摇摆没有“绝对功角”可言。工程上常用的做法是选择惯性中心COI做参考系。COI 的角度和角速度按下式定义δ_COI (Σ Mᵢ·δᵢ) / (Σ Mᵢ) ω_COI (Σ Mᵢ·ωᵢ) / (Σ Mᵢ)每台发电机相对 COI 的功角 θᵢ δᵢ − δ_COI角速度 ω̃ᵢ ωᵢ − ω_COI。然后在 COI 坐标系下构造能量函数V(θ, ω̃) ½·ΣMᵢ·ω̃ᵢ² V_p(θ)教案提到“两种坐标系”——绝对坐标系和 COI 坐标系。稳定平衡点的选取、能量函数的构造方式跟着变。初学多机系统先把 COI 坐标变换练熟再谈后续的 UEP 计算。多机系统的临界能量计算要比单机系统复杂得多教案给出的 AUEP 和 RUEP 法是两条不同的路径AUEP 法找到的是“最近不稳定平衡点”RUEP 法找的是“与故障轨迹相关的不稳定平衡点”。实际操作中AUEP 法偏保守RUEP 法更精确但计算量更大。教案第 83 页提到了一个实用做法扰动结束后系统实际穿越稳定边界的位置大概率集中在等位能线Potential Energy Boundary SurfacePEBS上沿着它去找相关的 UEP效率更高、定位也更稳。6.2 手写程序时的验证技巧纸上推公式是一回事把能量函数法写成程序是另一回事。我自己跑能量函数程序有一个固定的验证习惯先构造一个已知解析解的算例比如单机无穷大系统把程序算出的临界切除时间与时域仿真逐步扫描得到的结果对比误差控制在千分位以内再交给别人用。另外一个习惯是在程序里把动能项和势能项分开打印。多机系统里如果动能项在故障切除后下降而势能项上升说明系统正在把动能转化为势能是“向稳定走”的信号如果两者同向上升系统大概率在加速远离平衡点。这比只看总能量多了一个过程维度的判断依据。从那以后我每次做暂态稳定分析无论单机还是多机都强制自己先算一遍能量函数判据再回头按时域仿真校核一次两者对上了才敢写结论。这套流程虽然保守但能挡住大部分翻车事故。希望帮到你。本文还有配套的精品资源点击获取
简介暂态能量函数法是电力系统暂态稳定分析的重要直接法手段该PPT学习教案系统梳理了方法原理、数学描述与工程应用适合电气工程专业高年级本科生、研究生及从事电网稳定分析的工程技术人员使用。课件以古典力学能量概念切入逐步展开李雅普诺夫函数构造、临界能量判定、李氏定理与引力域等内容并延伸至单机无穷大系统直接法、多机系统特殊问题、两种坐标系、相关不稳定平衡点(RUEP)等专题脉络清晰、层层深入。资源包共1个pptx文件大小约1.23MB、全文共105页公式推导与示意图结合便于课堂授课或自学研读。目前已有161人浏览学习。学习后可掌握暂态能量函数的建模思路与稳定裕度评估方法理解故障轨迹及切除时刻对暂态稳定性的影响机制为电力系统暂态安全分析与控制实践打下扎实基础。1. 暂态能量函数法一套能算出稳定裕度的直接法而不是玄学做暂态稳定分析的人十有八九都有过这样的体会时域仿真曲线拉出来一大把功角摇摆看着是稳了但你要问他“稳定裕度是多少”“再切晚 0.1 秒还稳不稳”他答不上来。这不是仿真做得不够细而是时域法本质上是“答案有了过程没有”——你只知道稳不稳不知道离临界有多远。而暂态能量函数法TEFTransient Energy Function恰好补上这个缺口它不靠逐步积分推时间而是构造一个李雅普诺夫函数把“稳定与否”转化成“扰动结束时的系统能量是否超过临界能量”这一条不等式判断。这份 105 页的 PPT 教案正是把这条从原理到应用的完整链路讲清楚了——从古典力学的能量视角出发到单机无穷大系统的直接法推演再到多机系统的 UEP、RUEP、等位能线这些进阶话题。它适合两类人一是正在学电力系统分析的研究生需要把教材里“直接法”那一章真正吃透二是做在线稳定评估的工程师想理解能量函数类方法到底怎么落地、边界在哪。PPT 里公式推导和物理图景结合得紧不是那种公式堆到底的讲义读起来负担不大。2. 先把物理图景立住能量、临界能量与稳定裕度的关系2.1 从“小球滚碗”理解暂态能量函数的三个关键点这份教案的开篇用了非常经典的力学类比——一个小球在碗壁上运动。无扰动时小球停在碗底也就是稳定平衡点S.E.P受到扰动后小球获得动能和势能开始沿碗壁往上滚如果碗壁足够高、摩擦力存在小球最终会耗散能量回到碗底如果扰动太大小球翻过碗沿就回不来了。这个“碗沿”对应的位置就是不稳定平衡点U.E.P碗沿高度对应的势能就是临界能量 Vcr。由此得到暂态能量函数法的三个关键点找平衡点、构造暂态能量函数 V、确定临界能量 Vcr。这三件事环环相扣任何一个出了问题后面的稳定判断都不成立。在单机无穷大系统里这三个关键点都有明确的解析解稳定平衡点 δs arcsin(Pe/Pm)不稳定平衡点 δu π − δs临界能量 Vcr V(δu)。而在多机系统里这三个问题每一个都会变得棘手——这正是教案里花费大量篇幅讲 AUEP、RUEP 的原因。2.2 为什么坐标变换是构造能量函数的前提条件教案里有一个容易忽视但极其重要的细节在构造暂态能量函数前必须做坐标变换。原系统的稳定平衡点不在坐标原点直接构造 V 函数会导致 V(0) ≠ 0这不满足李雅普诺夫函数的基本条件当 X → 0 时 V → 0。所以要将稳定平衡点 δs 作为新坐标原点定义新状态变量 x δ − δs。这个操作看起来只是平移坐标轴实际上它决定了能量函数能不能在零点处取到最小值是后续一切推导的前提。这个概念在多少教材里被一笔带过但在实际工程中如果你自己写程序构造能量函数时发现 V 在平衡点处不为零第一件事就检查坐标变换做没做对。2.3 临界能量的“故障后网络依赖”属性教案第 16 页强调了另一个关键结论临界能量 Vcr 与系统故障后的网络结构有关也就是网络吸收发电机发出动能的能力。计算时用故障后的运动方式Vcr 也用故障后的参数值。这句话对做过实际计算的人很有启发。初学容易把故障前、故障中、故障后的网络混在一起用但能量函数法的逻辑是故障切除后的系统决定了“能装下多少能量”而故障期间的轨迹决定了“系统实际获得了多少能量”。Vcr 只看故障后网络V(δ) 按故障后的系统方程计算二者必须同源。提示做单机无穷大系统算例时先画出故障前、故障中、故障后三条功角特性曲线再动手列能量表达式。三条曲线对应的 Pe 表达式不同这一步错了后面全错。3. 李氏定理与数学基础为什么能量函数法不是“拍脑袋”3.1 李雅普诺夫稳定性判据的工程解读教案第 78 页给出了李氏定理的完整表述。把数学语言翻译成工程语言对一个自治系统 Ẋ f(X)如果能找到一个标量函数 V(x)满足 V(0) 0、V(x) 0x ≠ 0且 V̇(x) ≤ 0那么系统在原点稳定如果 V̇(x) 0则渐近稳定V̇(x) 0 则不稳定。这里需要注意“稳定”和“渐近稳定”的差别。V̇ ≤ 0 只能说明系统能量不再增长但不保证回到原点——对应物理场景就是小球停在碗壁某个位置不动了如果存在无阻尼的极限环之类的情况V̇ 0 才说明能量持续耗散最终回到平衡点。工程上我们通常希望渐近稳定。暂态稳定分析关注的不是“平衡点附近的小扰动稳定”而是“大扰动后能否回到平衡点或稳定运行点”所以需要用到引力域稳定域的概念——只有初态点落在某个区域内系统才是渐近稳定的这个区域就是引力域。单机无穷大系统的引力域边界恰好由不稳定平衡点 δu 划定。3.2 线性定常系统的李函数构造实例Lyapunov 方程的由来教案第 910 页给出了一个完整的线性系统算例这一步是理解后续非线性推广的根基。对线性系统 Ẋ AX取二次型 V(x) XᵀPX 为候选李函数则V̇ ẊᵀPX XᵀPẊ (AX)ᵀPX XᵀP(AX) Xᵀ(AᵀP PA)X要让 V̇ 负定需要 AᵀP PA 为负定矩阵。工程上通常反过来任选一个正定矩阵 Q常用单位阵令 AᵀP PA −Q解出 P再校验 P 是否正定。这就是 Lyapunov 方程。这个推导的价值在于它给出了一个“可验证”的流程先选 V 的形式再通过求解方程验证 V̇ 的符号最后回到 V 的正定性检查。每一步都有明确的操作对象不靠猜。后续在暂态能量函数法中构造 V 时思路一脉相承——只不过非线性系统的 V 构造没有系统化的方法更加依赖对物理过程的把握。3.3 从线性到非线性能量函数构造为什么难线性系统有 Lyapunov 方程兜底而非线性系统比如电力系统的摆动方程没有通用的构造方法。教案给出的路径是“从物理意义出发”动能部分取 K ½Mω²势能部分通过对功角方程右侧积分得到。这种做法有明确的物理支撑——动能是速度的二次型势能是功角相关项在整个轨迹上的积分。但是多机系统的能量函数构造远没有这么顺畅。不同发电机之间相对运动、网络结构变化、负荷模型差异都会影响 V 是否严格满足李雅普诺夫条件。这也是为什么教案后续花了大量篇幅讲“多机系统的特殊问题”和“两种坐标系”。4. 单机无穷大系统的直接法推演从运动方程到临界切除时间4.1 状态方程建立与坐标原点移动单机无穷大系统的经典二阶摆动方程为M·d²δ/dt² D·dδ/dt Pm − Pe·sin(δ)其中 M 为惯性时间常数D 为阻尼系数Pm 为机械功率Pe·sin(δ) 为电磁功率。令状态变量 x₁ δ − δs功角偏移、x₂ ω角速度偏差则状态方程为dx₁/dt x₂ dx₂/dt (Pm − Pe·sin(δs x₁) − D·x₂) / M注意这里已经完成了坐标平移δs 是稳定平衡点对应的功角。这一步在教案第 18 页交代得很清楚。为什么要这么写因为后续能量函数要在“δs 处取极小值”如果不平移坐标V 在平衡点处不为零整个李雅普诺夫判据就失效了。4.2 动能与势能函数的构造过程对状态方程的第二个式子两边乘以 M·dδ再从 δs 积分到当前功角 δ得到动能和势能的表达式动能来自速度项V_k ½·M·ω²势能来自不平衡功率对功角的积分V_p ∫[δs→δ] (Pe·sin(δ) − Pm) dδ关于这个势能表达式需要多说两句。积分下限取稳定平衡点 δs所以 V_p(δs) 0这是坐标变换的直接收益。积分内部是被积函数 Pe·sin(δ) − Pm 关于功角的积分——它描述了系统偏离平衡点后“存储”的势能。当 δ 在 δs 与 δu 之间时被积函数为正势能随 δ 增大而增加到达 δu 时系统达到临界势能。而临界能量 Vcr V_p(δu)因为 δu 处 ω 0不稳定平衡点处速度为 0动能为 0所以临界能量只包含势能部分。把 δu π − δs 代入势能表达式即可求得 Vcr 的解析值。4.3 暂态稳定判据与临界切除时间的求法得到了系统实际能量 V(δ, ω) 和临界能量 Vcr判断就变成一条不等式V(δ, ω) Vcr → 稳定 V(δ, ω) Vcr → 不稳定 V(δ, ω) Vcr → 临界状态判别稳定裕度时定义裕度 (Vcr − V) / Vcr × 100%在故障切除时刻δ δc、ω ωc计算此时的 V(δc, ωc)与 Vcr 比较就能判断切除后系统是否稳定。临界切除时间 tcr 的求法是在故障持续期间逐步推进时间计算每个时刻的 V(δ, ω)找到第一个使 V(δ, ω) Vcr 的时刻。# 伪代码示例临界切除时间扫描计算流程 import numpy as np # 系统参数标幺值示例 M 8.0 # 惯性时间常数 D 0.1 # 阻尼系数 Pm 0.8 # 机械功率 Pe 1.2 # 电磁功率幅值 # 稳定平衡点与不稳定平衡点 delta_s np.arcsin(Pm / Pe) delta_u np.pi - delta_s # 临界能量: 对 (Pe*sin(delta) - Pm) 从 delta_s 积分到 delta_u from scipy.integrate import quad Vcr, _ quad(lambda d: Pe * np.sin(d) - Pm, delta_s, delta_u) # 故障期间逐步积分, 求 V(delta, omega) 首次超过 Vcr 的时刻 dt 0.001 t 0.0 delta delta_s omega 0.0 while t 1.0: # 故障期间运动方程 (假设恒机械功率, 电磁功率按故障中曲线) Pe_fault 0.2 # 故障中电磁功率幅值 d_delta omega d_omega (Pm - Pe_fault * np.sin(delta) - D * omega) / M delta d_delta * dt omega d_omega * dt # 构造实际能量 (用故障后网络计算势能) V_k 0.5 * M * omega**2 V_p, _ quad(lambda d: Pe * np.sin(d) - Pm, delta_s, delta) V_total V_k V_p if V_total Vcr: print(f临界切除时间约: {t:.4f}s) print(f切除时刻功角: {np.degrees(delta):.2f}°) break t dt这段代码的推导逻辑与教案完全一致故障期间用故障中的功角特性曲线推进轨迹计算能量时必须用故障后的网络参数势能积分中被积函数用 Pe·sin(δ) − Pm而不是故障期间的 Pe_fault·sin(δ) − Pm。这是最容易翻车的区别——轨迹推进看故障中方程能量评估看故障后方程两边不要混。注意上面代码是求解流程示意不是完整工程脚本。实际算例中故障期间 Pe_fault 可能随电气量变化而变化需要结合网络方程联立求解。4.4 等面积定则与能量函数法的等价关系教案第 11 页画的功角特性曲线abcd 包围面积等于 defg 包围面积这正是等面积定则的图像化表述。把等面积定则翻译成能量函数语言加速面积代表故障期间系统获得的动能增量减速面积代表故障后系统能吸收的最大动能。二者相等时对应临界切除状态此时 V(δc, ωc) Vcr。等面积定则和能量函数法本质上是同一个判据的两种表达区别只在视角等面积定则从“面积”看能量函数法从“数值”看。理解了这个等价关系你会发现“V(δ, ω) 与 Vcr 比较”和“加速面积与减速面积比较”是同一件事这对后续理解多机系统中的能量函数拆分解读例如发电机动能项、网络势能项的物理含义非常有帮助。5. 避坑手册能量函数法应用中的五个高频翻车点5.1 翻车点一把 Vcr 当常数忽视了故障后网络依赖现象同一台发电机同一个故障位置换了故障线路阻抗后直接用之前的 Vcr 判断稳定性得出错误结论。原因Vcr 是故障后网络的函数网络拓扑或参数变了δu 变了Vcr 也变。教案第 16 页明确写了“计算时用故障后的运动方式Vcr 也用故障后的值”这句话很多人扫一眼就过去了。解决每次计算 Vcr 前先确认故障后网络的等值参数已更新。在程序中把 Vcr 的计算放在“故障后网络参数已生成”之后而不是放在故障前的初始化阶段。5.2 翻车点二能量积分路径选错——用故障前轨迹评估故障后能量现象算出来的稳度裕度偏大临界切除时间偏高按结果配置保护动作时间实际仿真却失稳。原因V(δ, ω) 的势能部分是从 δs 沿实际轨迹积到当前 δ但实际计算中若错误地用故障前网络的 Pe−δ 曲线去积分得到的是“错误的系统存储能量”该系统根本不存在于故障后的物理过程中。解决严格区分三条 Pe−δ 曲线——故障前、故障中、故障后。实际获取能量用故障后的曲线轨迹推进用故障中的曲线初始平衡点用故障前的曲线求 δs。5.3 翻车点三多机问题硬套单机公式现象把单机无穷大系统的公式直接套到多机系统的某一台发电机上算出“稳定”结果与实际仿真不符。原因多机系统不存在唯一的“功角”发电机之间存在相对运动能量函数的构造需要先选惯性中心COI坐标系把各机功角转换为相对 COI 的角度。教案第 3336 页的多机系统特殊问题部分就是在讲这件事。解决多机系统计算前先求系统的 COI 坐标完成坐标系变换。判断稳定时不能只看单机能量要看系统总能量动能和势能的代数加总与临界能量比较。5.4 翻车点四忽略了阻尼项对能量耗散的影响现象构造能量函数时忽略阻尼算出的稳定裕度偏保守得到的临界切除时间明显小于时域仿真结果。原因阻尼项 D·ω² 在能量表达式中作为耗散项会持续消耗系统动能对稳定有正贡献。教案中的推导先将阻尼设为 0 做简化第 14 页但在对比实际系统时如果不把阻尼效应考虑进去误差会很大。解决工程评估中如果系统阻尼不可忽略建议在能量函数中加入阻尼耗散项做修正。更实用的做法是用能量函数法给初判、用时域仿真给终判两者结合。5.5 翻车点五把 AUEP 当成 RUEP 用现象多机系统分析中有人用“最近不稳定平衡点”AUEP法估算临界能量找出来的点并不是实际故障轨迹最可能穿越的边界点结果或保守或冒进结果不一致。原因AUEP 找的是“距离当前稳定平衡点最近的所有不稳定平衡点中能量最低的那个”它与故障轨迹的具体形状无关而 RUEP相关不稳定平衡点是“故障轨迹实际趋近的不稳定平衡点”它才真正对应了故障后系统的失稳模式。解决优先用 RUEP 法做故障相关性的临界能量评估。工程上常用 BCUBoundary of Stability region via Controlling UEP方法先定位相关不稳定平衡点再计算 Vcr。6. 多机系统进阶两种坐标系、AUEP/RUEP 的选择与验证技巧6.1 COI 坐标系与惯量中心参考多机系统里每台发电机都在摇摆没有“绝对功角”可言。工程上常用的做法是选择惯性中心COI做参考系。COI 的角度和角速度按下式定义δ_COI (Σ Mᵢ·δᵢ) / (Σ Mᵢ) ω_COI (Σ Mᵢ·ωᵢ) / (Σ Mᵢ)每台发电机相对 COI 的功角 θᵢ δᵢ − δ_COI角速度 ω̃ᵢ ωᵢ − ω_COI。然后在 COI 坐标系下构造能量函数V(θ, ω̃) ½·ΣMᵢ·ω̃ᵢ² V_p(θ)教案提到“两种坐标系”——绝对坐标系和 COI 坐标系。稳定平衡点的选取、能量函数的构造方式跟着变。初学多机系统先把 COI 坐标变换练熟再谈后续的 UEP 计算。多机系统的临界能量计算要比单机系统复杂得多教案给出的 AUEP 和 RUEP 法是两条不同的路径AUEP 法找到的是“最近不稳定平衡点”RUEP 法找的是“与故障轨迹相关的不稳定平衡点”。实际操作中AUEP 法偏保守RUEP 法更精确但计算量更大。教案第 83 页提到了一个实用做法扰动结束后系统实际穿越稳定边界的位置大概率集中在等位能线Potential Energy Boundary SurfacePEBS上沿着它去找相关的 UEP效率更高、定位也更稳。6.2 手写程序时的验证技巧纸上推公式是一回事把能量函数法写成程序是另一回事。我自己跑能量函数程序有一个固定的验证习惯先构造一个已知解析解的算例比如单机无穷大系统把程序算出的临界切除时间与时域仿真逐步扫描得到的结果对比误差控制在千分位以内再交给别人用。另外一个习惯是在程序里把动能项和势能项分开打印。多机系统里如果动能项在故障切除后下降而势能项上升说明系统正在把动能转化为势能是“向稳定走”的信号如果两者同向上升系统大概率在加速远离平衡点。这比只看总能量多了一个过程维度的判断依据。从那以后我每次做暂态稳定分析无论单机还是多机都强制自己先算一遍能量函数判据再回头按时域仿真校核一次两者对上了才敢写结论。这套流程虽然保守但能挡住大部分翻车事故。希望帮到你。本文还有配套的精品资源点击获取
