电力系统近似常数与近似公式
1. 引言
电力系统横跨多个数量级：电压等级从 10 kV 到 1000 kV，输送容量从 kW 到 GW，网络规模从单个变电站到跨越大陆的互联电网。然而，与这种宏观上的巨大差异形成对比的是，系统中许多重要的参数在数量级上却相当保守——它们被底层物理定律（电磁场方程、热力学定律）和工程约束（绝缘强度、材料极限）锁定在一个相对狭窄的区间内，并没有随规模同步扩张。这些”近似常数”是加深电力系统分析直觉的宝贵线索，也可以作为筛查数据合理性的辅助判据——当某一参数明显偏出这些区间时，往往值得重新核查。
电力系统的运行分析也发展出了一套超出一般教科书范畴的近似公式体系。这些公式建立在若干简化假设之上：把多机系统折算成单机等值，把复杂负荷线性化成频率阻尼，把分散机组的一次调频等效成一阶环节，把长线路用波阻抗和自然功率概括，把暂态稳定计算简化成线性化冲击量。这样做的代价是结果必然有偏差；但它的价值同样非常明确：它能把仿真结果背后的物理本质暴露出来，帮助分析人员在详细仿真之前先建立对系统的概念性直观理解。换言之，近似公式并不是为了替代详细仿真，而是让人在看到仿真结果时能够理解”为什么是这个量级、为什么朝这个方向变化”。
本文将这两个层次统一纳入一个框架，按如下思路组织：第 2 节约定全文统一的符号体系，避免后续公式推导中出现歧义；第 3 节梳理电力系统设备参数中的近似常数，理解为什么这些参数在宽广的工程范围内保持稳定；第 4—8 节在此基础上推导并解释频率动态、机电振荡、电压稳定、线路无功和暂态稳定等核心场景中最有工程解释力的近似公式；第 9 节给出一组综合算例，展示这些公式在工程估计中的实际用法；第 10 节讨论近似公式的边界与误差来源。
2. 符号约定
为避免单位混乱，先作以下约定。
标幺值与有名值：本文频率动态部分默认采用统一基准下的标幺化量，即功率缺额、频率偏差、负荷频率系数、一次调频系数等均已按同一基准值归一化。若改用有名值（MW、Hz），则必须同步修改各系数单位；例如若 $\Delta f$ 以 Hz 表示、$\Delta P$ 以 MW 表示，那么 $k_D$、$k_G$ 应以 MW/Hz 表示，不能把标幺公式和有名公式直接混用。
本文主要符号如下：
| 符号 | 含义 |
|---|---|
| $H$ | 单台发电机惯性常数（s），定义为转子额定动能与额定容量之比 $H = E_k / S_N$ |
| $T_s$ | 系统等值惯性时间常数（s），用于频率动态，对应全网折算惯性，$T_s = 2H$ |
| $T_j$ | 单台发电机惯性时间常数（s），用于单机振荡与暂态稳定分析，$T_j = 2H$ |
| $T_G$ | 一次调频综合时间常数（s） |
| $k_D$ | 负荷频率调节系数（标幺值） |
| $k_G$ | 发电机组功频静态特性系数（标幺值） |
| $k_s$ | $k_D + k_G$，系统综合调频系数 |
| $K_s$ | 同步转矩系数（标幺值），贯穿机电振荡与暂态稳定分析 |
| $E’_q$ | 发电机暂态电势（标幺值） |
| $U$、$U_g$ | 母线电压、送端电压（标幺值） |
| $X_\Sigma$ | 等值电抗（标幺值） |
| $Z_c$ | 线路波阻抗（Ω） |
| $P_N$ | 线路自然功率（MW） |
| $\omega_0$ | 系统额定角频率（rad/s），50 Hz 系统中 $\omega_0 = 2\pi \times 50 \approx 314$ rad/s |
注意 $T_s$ 与 $T_j$ 的含义不同：$T_j$ 针对单台机组，$T_s$ 是全网折算值；在单机无穷大系统的分析中两者才会对应同一台机组，在多机等值场景中 $T_s$ 应理解为全网惯性的加权合成。
3. 电力设备的近似常数
掌握符号体系之后，在进入动态公式推导之前，有必要先建立对设备参数本身的直觉。电力系统设备参数中存在若干近似常数，它们在宽广的容量和电压范围内保持相对稳定，其背后都有清晰的物理约束机制。这些常数是后续近似公式得以成立的重要前提，尤其是第 3.6 节关于发电机惯性时间常数 $H$ 的讨论，将直接为第 4、5、8 节的近似公式提供数量级基准。
3.1 单根导线载流量的”恒定性”
不论电压等级如何，组成架空线路的单根子导线额定载流量始终维持在 400 A 到 1000 A 左右。
这一现象的根本原因在于导线的热平衡机制：导线的稳态载流量取决于焦耳热与表面散热的平衡。最大允许电流 $I_{\max}$ 与导线半径 $r$ 的关系服从
\[I_{\max} \propto r^{1.5}\]由于趋肤效应以及粗导线在制造、展放上的机械困难，单根导线的物理直径被严格限制，从而锁死了单根导线的电流上限。
工程实践恰好印证了这一规律：10 kV 配电架空线路使用单根导线，载流量约 400 A；而 500 kV 骨干网架为输送大功率，采用四分裂导线（如每根子导线 577 A），平均单根导线的载流量依然在这个区间内。电压等级的提高，是通过增加分裂数来增大整体载流量，而非提高单根导线本身的极限。
3.2 架空线路的波速与波阻抗
对于架空输电线路，波速的理论公式为
\[v = \frac{1}{\sqrt{LC}} = \frac{1}{\sqrt{\mu_0 \varepsilon_0}} \approx 3 \times 10^5 \text{ km/s}\]这正是真空中光速的定义，其中 $L$ 是单位长度电感（H/m），$C$ 是单位长度电容（F/m）。考虑趋肤效应及大地回流损耗，实际波速会略低，但依然非常接近光速。
同样值得关注的是波阻抗（也称浪涌阻抗）
\[Z_c = \sqrt{\frac{L}{C}}\]波阻抗在不同电压等级之间的变化不大（虽然不能称为严格恒定）。这一特性有直接的工程意义：线路的自然功率基本只与运行电压有关，而与线路的具体参数关系较弱，从而为超高压线路的无功平衡分析提供了简便的估算路径——这一点将在第 7 节详细展开。
3.3 线路单位长度正序电抗
在 50 Hz 频率下，三相架空线路的正序电抗普遍在 0.2 到 0.4 Ω/km，无论电压等级是多少。这一规律的理论解释来自几何均距公式：
\[x_1 \approx 0.1445 \log_{10}\!\left(\frac{D_{eq}}{D_s}\right) + 0.0157 \quad (\Omega/\text{km})\]随着电压等级升高，相间绝缘要求导致等效几何均距（$D_{eq}$）变大；但为了抑制电晕，必须同步采用分裂导线，导致自几何均距（$D_s$）也相应变大。分子分母同步放大，对数值变化有限，因此 $x_1$ 在宽广的电压范围内保持稳定。
3.4 典型线路参数汇总
下表给出不同电压等级的典型线路参数（来源于电力规划设计总院《电力系统设计手册》），传播速度和波阻抗 $Z_c$ 由 $X_1$ 和 $C_1$ 换算得出。
| 电压等级 | 线路型号 | $R_1$（Ω/km） | $X_1$（Ω/km） | $C_1$（μF/km） | 传播速度（km/s） | $Z_c$（Ω） |
|---|---|---|---|---|---|---|
| 1000 kV | 8×LGJ-500 | 0.0076 | 0.264 | 0.0139 | 292600 | 246 |
| 500 kV | 4×LGJQ-300 | 0.027 | 0.280 | 0.0128 | 296000 | 264 |
| 220 kV | 2×LGJJ-240 | 0.065 | 0.303 | 0.0119 | 295400 | 285 |
| 110 kV | LGJ-150 | 0.63 | 0.442 | 0.00808 | 296700 | 417 |
从表中可以直观地看出：$X_1$ 在前三个电压等级几乎相同，传播速度在所有等级都接近光速，波阻抗随电压等级升高而缓慢下降（高压线路自然功率更大的原因之一，但主要影响因素还是电压等级不同）。
3.5 变压器短路阻抗百分比的一致性
变压器的短路阻抗百分比（$U_k\%$）具有跨越容量的一致性：从 10 kVA 的柱上变压器到 1000 MVA 的特高压主变压器，$U_k\%$ 通常落在 4% 到 15% 的狭窄带内（特高压主变略高，约 18%—24%）。
这一区间由两股对立的约束力共同决定。下限约束（机械强度）：短路时绕组受到的电磁力 $F \propto (1/U_k\%)^2$，$U_k\%$ 过小会导致短路电磁力呈平方级暴增，瞬间损坏变压器。上限约束（电网运行）：$U_k\%$ 过大会导致正常带负荷时的电压调整率过高，使得变压器二次侧电压随负荷大幅波动，难以满足供电质量要求。两种约束的交汇，将 $U_k\%$ 锁定在相对固定的区间内，并不随容量的增大而发生数量级的跳变。
3.6 旋转电机的两个近似常数
发电机机端电压的”物理天花板”
发电机容量从 10 MW 跃升至 1000 MW，机端额定电压仅从 10.5 kV 缓慢提升至 27 kV 左右，存在明显的物理上限。原因在于定子槽内空间极为宝贵：若大幅提高机端电压，所需的绝缘层厚度会急剧增加，而绝缘层导热极差，还会严重挤占定子铜导线的截面积，导致电机有效输出功率反而下降。
惯性时间常数 $H$ 的稳定性
各种容量的同步发电机，其惯性时间常数 $H$ 普遍维持在 2 到 6 秒的范围内。$H$ 恒定的根本原因，在于发电机的功率密度与转动能量密度受到同一组物理约束的锁定：
- 功率密度的上限由定子槽内的电磁负荷决定——导体电流密度受铜的导热极限约束，气隙磁通密度受硅钢片饱和约束，使得每单位体积所能输出的电功率大致恒定；
- 转动能量密度的上限由转子材料强度决定——转子表面线速度受钢材屈服强度约束，这锁定了单位质量转子所能储存的最大动能。
由于两者都被材料的物理极限锁定，且随容量增大时转子体积的增长同时贡献于动能和功率，比值 $H$ 在宽广的容量范围内自然保持稳定。
这一结论将在后文中反复用到：第 4 节的频率动态分析中，全网等值惯性时间常数 $T_s$ 正是由各机组 $H$ 加权合成，大型互联电网一般在 8—15 s 量级；第 5 节和第 8 节分析机电振荡与暂态稳定时，单机 $T_j$ 的典型区间（4—10 s）同样直接影响振荡频率和临界切除时间的数量级估计。
4. 频率动态近似：从一阶模型到含一次调频的二阶模型
上一节建立了对设备参数的直觉，尤其是发电机惯性时间常数 $H$ 决定了全网等值惯性 $T_s$ 的数量级。本节从摆动方程出发，以 $T_s$ 为核心参数，推导事故后频率响应过程中最有工程价值的四个量：初始频率变化率、稳态频差、最低点出现时刻，以及频率最低点。
4.1 一阶模型：只考虑惯性与负荷频率特性
频率动态最基本的出发点是摆动方程。若把系统等效为单机带集中负荷，并在事故后最初阶段把机械功率看成暂时不变，写成功率形式为
\[M\frac{d\omega}{dt} = P_m - P_L\]若负荷具有频率特性，可写作 $P_L = P_0(f/f_0)^{K_L}$，在额定频率附近线性化，得
\[\Delta P_L \approx D\,\Delta f, \qquad D = \frac{K_L P_0}{f_0}\]将物理惯性 $M$ 转换为系统等值惯性时间常数 $T_s$，事故造成初始有功缺额 $\Delta P_{OL0}$ 时，得到一阶频率方程
\[T_s\frac{d\Delta f}{dt} + k_D\Delta f = -\Delta P_{OL0}\]令 $T_f = T_s/k_D$，则频率偏差的一阶解为
\[\Delta f(t) = -\frac{\Delta P_{OL0}}{k_D}\left(1 - e^{-t/T_f}\right)\]这个结果揭示了三件重要事实：第一，若系统只有负荷频率特性而没有一次调频响应，则频率跌落是单调的，不会出现”先跌到最低点再回升”的过程；第二，时间常数 $T_f = T_s/k_D$ 越大，跌落越慢；第三，稳态频差为
\[\Delta f_\infty = -\frac{\Delta P_{OL0}}{k_D}\]这意味着仅靠负荷频率特性，频率只能停在偏离额定值的新平衡点——这正是一次调频和 AGC 存在的根本原因，也引出了下一节需要考虑的二阶模型。
4.2 二阶模型：计入一次调频
真正对运行人员最有价值的，是把一次调频考虑进去之后的二阶近似模型：
\[T_s\frac{d\Delta f}{dt} = -\Delta P_{OL}\] \[T_G\frac{d\Delta P_G}{dt} + \Delta P_G = -k_G\Delta f\] \[\Delta P_D = k_D\Delta f, \qquad \Delta P_{OL} = \Delta P_D - \Delta P_G + \Delta P_{OL0}\]这个模型的结构非常清楚：惯性决定”刚出事时系统先以多快的速度掉下去”，负荷频率特性和一次调频共同决定”之后能不能把下跌拉住、拉到什么程度”。
将方程联立消去中间量，得关于 $\Delta f$ 的二阶常系数微分方程：
\[T_sT_G\frac{d^2\Delta f}{dt^2} + \left(T_s + k_DT_G\right)\frac{d\Delta f}{dt} + \left(k_D + k_G\right)\Delta f + \Delta P_{OL0} = 0\]令 $k_s = k_D + k_G$，$T_f = T_s/k_D$，并定义
\[\alpha = \frac{1}{2}\left(\frac{1}{T_G} + \frac{1}{T_f}\right), \qquad \Omega = \sqrt{\frac{k_s}{T_sT_G} - \alpha^2}\]当 $k_s/(T_sT_G) > \alpha^2$ 时，系统属于欠阻尼情形，频率响应可写成
\[\Delta f(t) = -\frac{\Delta P_{OL0}}{k_s}\left[1 - 2A_m e^{-\alpha t}\cos(\Omega t + \varphi)\right]\]其中
\[A_m = \frac{1}{2\Omega T_s}\sqrt{k_sk_G}, \qquad \varphi = \arctan\!\left(\frac{1}{\Omega}\left(\frac{k_s}{T_s} - \alpha\right)\right)\]这个表达式包含了频率动态中最重要的三个特征：整体负偏移（尺度由 $\Delta P_{OL0}/k_s$ 决定）、指数衰减率（主要由 $\alpha$ 决定），以及摆动与回升的节奏（主要由 $\Omega$ 决定）。
4.3 四个最有工程价值的量
（1）初始频率变化率
在事故刚发生的瞬间，一次调频和负荷频率响应都来不及充分建立，因此
\[\boxed{\left.\frac{d\Delta f}{dt}\right|_{0^+} = -\frac{\Delta P_{OL0}}{T_s}}\]初始下降率只取决于功率缺额和系统惯性。惯性越大，系统在事故起点越”钝”；缺额越大，跌落起点越”陡”。结合第 3.6 节的结论，在高比例新能源系统中传统同步机被逆变器替代，等值惯性下降，频率脆弱性会显著增大。
（2）稳态频差
令 $t \to \infty$，指数项衰减为零，则
\[\boxed{\Delta f_\infty = -\frac{\Delta P_{OL0}}{k_s} = -\frac{\Delta P_{OL0}}{k_D + k_G}}\]负荷频率特性和一次调频静态特性在稳态上并列贡献。要真正恢复到原频率，还需要二次调频（AGC）参与。
（3）最低点出现时刻
欠阻尼情形下，令频率响应的导数为零，可求出频率最低点出现时刻 $t_m$。写成最安全的形式是
\[t_m = \frac{1}{\Omega}\operatorname{atan2}\!\left(2T_sT_G\Omega,\; k_DT_G - T_s\right)\]这里特意使用 $\operatorname{atan2}$ 而非普通反正切，是因为分母 $k_DT_G - T_s$ 在实际参数下可能为负，若直接用 $\arctan$ 主值，极易把象限选错，将最低点时刻算成负值或错误分支。
（4）频率最低点
将 $t_m$ 代入频率响应式，并利用最低点处导数为零的条件，可整理为
\[\boxed{\Delta f_{\min} = -\frac{\Delta P_{OL0}}{k_s}\left(1 + \sqrt{\frac{k_GT_G}{T_s}}\;e^{-\alpha t_m}\right)}\]这是欠阻尼且存在明确最低点的前提下成立的表达式。若系统参数使二阶响应变成过阻尼，则曲线会单调下滑至稳态，不再出现”先下后回”的最低点，上述公式不应机械套用。最低频率的有名值为 $f_{\min} = f_0(1 + \Delta f_{\min})$。
4.4 参数影响的物理解读
对上述公式作一个系统性的参数敏感性梳理，有助于在工程中快速定位问题的根因：
- $T_s$ 增大：直接减小初始频率下降率，但对稳态频差和最低点几乎没有直接影响。真正决定”频率跌落有多陡”的，是惯性而不是一次调频参数。
- $k_G$ 增大：显著减小稳态频差，同时通常能抬高最低点。但 $k_G$ 并不决定调频动作有多快，动作快慢更多体现在 $T_G$ 上。
- $T_G$ 越大：一次调频建立越慢，最低点通常越低、出现时间也越晚。”稳态差不多，但最低点很深”的问题往往不在于静态调差不够，而在于调频建立得太慢。
- $k_D$ 的作用常被低估：负荷频率特性不是机组控制器，但对事故后的实际频率过程有持续的”自然阻尼”作用，是无需控制指令的第一层缓冲机制。
5. 机电振荡频率的近似公式
上一节分析的频率动态，关注的是全网层面、秒级至分钟级的有功-频率耦合过程，核心参量是全网等值惯性 $T_s$。当视角切换到单台机组或机群与系统之间的相对摆动，则进入了机电振荡的分析领域——其时间尺度更短（通常 0.5—5 s），关注的核心量变为单机惯性时间常数 $T_j$ 和同步转矩系数 $K_s$。
5.1 基本推导
机电振荡最经典的近似来自线性化摆动方程。设发电机与系统之间的功角为 $\delta$，机组惯性时间常数为 $T_j$，在运行点 $\delta_0$ 附近作小扰动，忽略阻尼转矩和机械功率扰动，则摆动方程为
\[\frac{T_j}{\omega_0}\frac{d^2\Delta\delta}{dt^2} = -\Delta P_e\]采用暂态电势 $E’q$、系统侧电压 $U$ 和等值电抗 $X\Sigma$，电磁功率满足功角关系
\[P_e = \frac{E'_qU}{X_\Sigma}\sin\delta\]在运行点 $\delta_0$ 处线性化，得
\[\Delta P_e = K_s\,\Delta\delta, \qquad K_s = \frac{E'_qU}{X_\Sigma}\cos\delta_0\]$K_s$ 称为同步转矩系数，反映功角偏离平衡点时电磁功率的恢复能力。代入摆动方程，得到标准无阻尼振荡方程，其固有振荡频率为
\[\boxed{f_n = \frac{1}{2\pi}\sqrt{\frac{\omega_0 K_s}{T_j}}}\]同步转矩系数也可改写为
\[K_s = P_0\cot\delta_0\]这个写法很有启发性：在相同出力下，功角越大，$\cot\delta_0$ 越小，同步”刚度”越弱，振荡频率也就越低——这正是重载、远距离输电更容易出现低频振荡的微观机制。结合第 3.6 节，$T_j$ 的典型值在 4—10 s，代入典型的 $K_s \approx 1$—$2$ pu，可以验证上式给出的振荡频率确实落在 1—2 Hz 的局部模式范围内。
5.2 一个经常被误解的问题
在一些非正式表述中，常会听到”低频振荡是调速器或励磁调节器引起的”。严格来说，这种说法并不准确。振荡频率本身主要由同步刚度和惯性决定，即由 $K_s$ 和 $T_j$ 决定；调速器、励磁器、PSS 等控制装置更主要影响的是阻尼和模态耦合，从而影响”振不振得起来、振幅会不会被放大、衰减快不快”，而不是从根本上决定一个本来不存在的振荡频率。
更严谨的说法应当是：扰动激发系统本来就存在的机电模态；控制器参数改变的主要是这些模态的阻尼比和参与因子，有时也会通过改变运行点或等值刚度间接改变频率，但其首要作用并不是”凭空产生振荡频率”。
值得注意的是，第 5 节中的 $K_s$ 与第 8 节暂态稳定分析中的 $K_s^{\text{post}}$ 在形式上完全相同，但后者取的是故障后网络的参数。这并非偶然：线性化摆动方程在小扰动和冲击法近似中共享同一物理结构。
6. 电压相关的近似公式
频率和功角的近似分析主要关注有功功率和旋转机械的动态。而电力系统中还有一个独立但同等重要的维度：母线电压与无功功率之间的平衡。本节给出三个最常用的电压近似公式，从静态电压稳定极限到母线电压快算，再到电容器投切效果估算，由宏观到微观，逐层递进。
6.1 受端负荷最大有功极限
对一个典型的二节点系统：电源电压为 $U_g$，经总电抗 $X_\Sigma$ 向受端负荷供电，忽略电阻，若负荷功率因数角为 $\varphi$，则静态电压稳定极限对应的最大可送有功为
\[\boxed{P_{L,\max} = \frac{U_g^2}{2X_\Sigma}\frac{\cos\varphi}{1 + \sin\varphi}}\]当功率因数为 1 时，公式退化为 $P_{L,\max} = U_g^2/(2X_\Sigma)$。这个结果揭示了两个要点：线路电抗越大，电压稳定极限越低；负荷功率因数越滞后，极限也越低。它和第 8 节暂态稳定中的功角极限不是一回事，但两者都体现出一个共同事实：系统”刚度”下降时，可承受的传输负荷也随之下降。
6.2 受端最低电压
在同样的假设下，负荷达到静态电压稳定极限时，受端电压达到最小值：
\[\boxed{V_{\min} = \frac{1}{\sqrt{2 + 2\sin\varphi}}}\]其中电压已按送端基准归一化。几个典型值非常值得记忆：
- 当 $\cos\varphi = 1$ 时，$V_{\min} = 1/\sqrt{2} \approx 0.707$
- 当 $\cos\varphi = 0.95$ 时，$V_{\min} \approx 0.617$
- 当 $\cos\varphi = 0.707$（即 $P = Q$）时，$V_{\min} \approx 0.541$
运行人员可从这些数值直接获得直觉：负荷功率因数越差，受端电压一旦跌破某个阈值，系统实际上已经非常接近静态电压崩溃边缘。
6.3 用 $P$、$Q$ 快速估算母线电压及电容器投切效果
把计算点母线对外部系统等值成戴维南电源（系统等值电势 $E$，系统等值阻抗 $jX_i$），若已知该母线从系统侧收到的有功、无功分别为 $P$、$Q$，消去相位差后可得
\[P^2 + \left(Q + \frac{U^2}{X_i}\right)^2 = \frac{E^2U^2}{X_i^2}\]令 $S = E^2/X_i$（理解为该母线看向系统的短路容量），令 $y = (U/E)^2$，则上式化为
\[y^2 + \left(\frac{2Q}{S} - 1\right)y + \frac{P^2 + Q^2}{S^2} = 0\]取接近运行点的物理解：
\[\left(\frac{U}{E}\right)^2 = \frac{1}{2}\left[1 - \frac{2Q}{S} + \sqrt{\left(1 - \frac{2Q}{S}\right)^2 - \frac{4(P^2 + Q^2)}{S^2}}\right]\]当系统较强、$P,Q \ll S$ 时，可进一步展开为
\[\frac{U}{E} \approx 1 - \left(\frac{Q}{S} + \frac{P^2}{2S^2}\right)\]这个式子的工程含义是：母线电压对无功最敏感，对有功通常只表现为二阶效应。因此运行中一旦看到母线电压波动，第一反应通常应当先检查无功平衡，而不是先怀疑有功潮流本身。
由上式作微分，当研究对象是母线投切电容器时，通常 $\Delta P \approx 0$，于是立刻得到最实用的快算式
\[\boxed{\Delta U\,(\text{pu}) \approx \frac{\Delta Q_c}{S}}\]其中 $\Delta Q_c$ 为电容器注入母线的容性无功（约定投入为正），$S$ 为该母线的短路容量。这一公式揭示了一个极其重要的运行事实：同样一组电容器，在强系统母线上的电压效果很有限，而在弱系统母线上的电压效果会明显放大。 这也是为什么相同容量的并联电容器装在不同短路容量的母线上，投切后电压跳变量会差很多。
7. 线路自然功率与无功损耗的近似公式
第 3.2 节指出，波阻抗 $Z_c$ 在不同电压等级之间的变化不大。这一特性的直接推论，是线路自然功率可以用一个简洁的公式估算，并由此衍生出长线路无功行为的一个极为实用的判据。这与第 6 节的电压-无功分析形成互补：第 6 节关注的是母线节点的无功平衡，本节关注的是传输线路本身的无功特性。
7.1 自然功率
对无损长线路，线路在额定电压下的自然功率为
\[\boxed{P_N = \frac{U^2}{Z_c}}\]这里的物理意义非常深刻：当线路传输的有功恰好等于自然功率时，线路电感消耗的无功和线路电容发出的无功大体平衡，线路对外表现出近似”无功自平衡”的状态。
7.2 无功损耗估算
若线路额定电压下单位长度的充电功率为 $Q_N$，线路长度为 $l$，而实际传输有功为 $P$，则线路无功功率的近似损耗（或盈余）为
\[\boxed{\Delta Q_L = \left[\left(\frac{P}{P_N}\right)^2 - 1\right]Q_N l}\]这个公式有两个非常重要的判断：
- 若 $P < P_N$，则 $\Delta Q_L < 0$，说明线路总体上发无功（容性占优）；
- 若 $P > P_N$，则 $\Delta Q_L > 0$，说明线路总体上吸无功（感性占优）。
自然功率因此不只是一个”功率值”，更是一个分界点：它把长线路的无功行为分成了净容性区和净感性区，直接决定了电压控制策略的方向。这一判断与第 6.3 节的电容器投切效果估算可以联合使用：线路的无功净值给出系统侧需要补偿的方向，电容器投切公式则给出补偿量对应的电压变化。
8. 暂态稳定极限的近似估算
前面几节分别从频率（第 4 节）、功角小扰动（第 5 节）和电压-无功（第 6、7 节）三个维度建立了近似分析框架。本节回到功角稳定问题，但关注的不再是线性化的小扰动振荡，而是大扰动——故障期间能量冲击与系统在故障切除后恢复稳定之间的对抗。这是暂态稳定分析的核心，其准确评估需要完整的时域仿真，但工程快估同样有价值。
8.1 冲击法的工程思想
严格的暂态稳定分析采用等面积法或完整时域仿真。对于单机无穷大系统，等面积法写成
\[\int_{\delta_0}^{\delta_c}(P_m - P_{e,f})\,d\delta = \int_{\delta_c}^{\delta_{\max}}(P_{e,\text{post}} - P_m)\,d\delta\]在工程快估中，冲击法的思路是：故障相当于给转子”踹了一脚”，这一脚的力度（加速功率 $\Delta P_a$）乘以作用时间（$\Delta t$），决定了转子积累了多少”初速度”；系统的弹性（同步转矩系数 $K_s^{\text{post}}$）和惯性（$T_j$）共同决定这个初速度最终会弹出多大的功率摆幅。这与第 5 节小扰动振荡的物理图像一脉相承，只是从线性化的自由振荡推广到了大扰动下的强迫响应。
8.2 推导过程
第一步：故障切除时转子积累的角速度偏差
故障持续期间加速功率近似为常数 $\Delta P_a$，在故障持续时间 $\Delta t$ 内积分，得切除瞬间转子角速度偏差：
\[\Delta\omega_{cl} = \frac{\omega_0\,\Delta P_a}{T_j}\,\Delta t\]第二步：故障切除后的自由振荡
故障切除后，转子在故障后网络的新平衡点 $\delta_s = \arcsin(P_m/P_{\max}^{\text{post}})$ 附近自由振荡，其振荡角频率为
\[\omega_n = \sqrt{\frac{\omega_0 K_s^{\text{post}}}{T_j}}, \qquad K_s^{\text{post}} = P_{\max}^{\text{post}}\cos\delta_s\]第三步：从角度振幅换算到功率振幅
功率振荡的单边振幅（最大超调量）为
\[A = K_s^{\text{post}} \cdot \frac{\Delta\omega_{cl}}{\omega_n} = \omega_n \cdot \Delta P_a \cdot \Delta t\]写成频率形式 $\omega_n = 2\pi f_d$，即得
\[\boxed{A = 2\pi f_d \cdot \Delta P_a \cdot \Delta t}\]重要提示：此处 $f_d$ 应取故障后网络的参数（$K_s^{\text{post}}$ 和 $T_j$）计算，不能沿用故障前运行点。
8.3 稳定条件与暂稳极限
第一摆中，功率最高点不得超过故障后网络所能传输的最大功率 $P_{\max}^{\text{post}}$，否则系统失步。由此得暂稳极限的估算值：
\[\boxed{P_{st} \approx P_{\max}^{\text{post}} - D_p}\]其中 $D_p = 2\pi f_d \cdot \Delta P_a \cdot \Delta t$ 采用单边振幅定义。某些资料将 $D_p$ 定义为双边总摆幅，则公式改写为 $P_{st} \approx P_{\max}^{\text{post}} - D_p/2$，使用前必须先确认 $D_p$ 的定义，这是冲击法手算最常见的错误来源。
8.4 临界切除角与临界切除时间
冲击法回答”给定切除时间，摆幅有多大”。继电保护人员更关心的另一面是：在多长时间内必须切除故障，系统才不会失步？ 对故障期间 $P_{e,f} \approx 0$ 的三相金属性短路，等面积条件给出临界切除角
\[\delta_{cr} = \arccos\!\left[\frac{P_m}{P_{\max}}(\pi - 2\delta_0) - \cos\delta_0\right]\]由匀加速积分关系得临界切除时间
\[\boxed{t_{cr} \approx \sqrt{\frac{2T_j(\delta_{cr} - \delta_0)}{\omega_0\,P_m}}}\]其中 $\delta_{cr} - \delta_0$ 以弧度计，$P_{\max}$ 应取故障后网络的最大传输功率。
冲击法与临界切除时间公式形成双层结构：先用冲击法快速判断风险量级，再用 $t_{cr}$ 公式校核继电保护动作时间是否充裕，两者从不同方向刻画同一物理约束，互为补充。
9. 综合算例
以下给出一组综合估算算例，演示各近似公式在”调度快估场景”中的实际用法。各小题共用同一系统背景，目的在于展示计算流程和物理解读，而不是进行精确仿真。
9.1 事故后的频率最低点估算
设某受端系统以 50 Hz 运行，系统基准容量取 10 GW。一次大扰动造成有功缺额 800 MW，即 $\Delta P_{OL0} = 0.08$ pu。等值参数取
\[T_s = 8\text{ s}, \quad T_G = 5\text{ s}, \quad k_D = 1.2, \quad k_G = 4.0\]则 $k_s = k_D + k_G = 5.2$，$T_f = T_s/k_D = 6.667$ s。进一步有
\[\alpha = \frac{1}{2}\left(\frac{1}{5} + \frac{1}{6.667}\right) = 0.175\text{ s}^{-1}, \qquad \Omega = \sqrt{\frac{5.2}{8 \times 5} - 0.175^2} \approx 0.315\text{ s}^{-1}\]初始频率变化率
\[\left.\frac{df}{dt}\right|_{0^+} = -\frac{0.08}{8} \times 50 = -0.50\text{ Hz/s}\]稳态频差
\[\Delta f_\infty = -\frac{0.08}{5.2} \approx -0.769\text{ Hz}, \quad f_{ss} \approx 49.231\text{ Hz}\]最低点时刻
\[t_m = \frac{1}{0.315}\operatorname{atan2}(2 \times 8 \times 5 \times 0.315,\;1.2 \times 5 - 8) \approx 5.23\text{ s}\]频率最低点
\[\Delta f_{\min} = -\frac{0.08}{5.2}\left(1 + \sqrt{\frac{4.0 \times 5}{8}}\,e^{-0.175 \times 5.23}\right) \approx -0.02512\text{ pu}\]故 $f_{\min} = 50 \times (1 - 0.02512) \approx 48.744$ Hz。
作为对照，若完全不考虑一次调频（纯一阶模型），则在同一时刻频率约为 48.185 Hz，最低点将深约 0.56 Hz。这说明一次调频对最低点的抬升不只是改善稳态，而是显著改变了事故后前几秒最危险的那一段过程。
9.2 机组机电振荡频率估算
设某送端机组经网络等值电抗 $X_\Sigma = 0.55$ 与系统相连，$E’_q = 1.12$ pu，$U = 1.0$ pu，初始有功出力 $P_0 = 0.8$ pu，$T_j = 9$ s。
先求初始功角：$\delta_0 = \arcsin(0.8 \times 0.55 / 1.12) \approx 23.1°$。
同步转矩系数为
\[K_s = \frac{1.12}{0.55}\cos 23.1° \approx 1.873\]故小扰动振荡频率为
\[f_n = \frac{1}{2\pi}\sqrt{\frac{314 \times 1.873}{9}} \approx 1.29\text{ Hz}\]这个数值落在典型局部机电振荡频率范围内。若线路切除使 $X_\Sigma$ 增大，或机组出力升高使 $\delta_0$ 增大，则 $K_s$ 会下降，振荡频率也会随之降低——这正是输电断面潮流增大时振荡风险升高的微观机制。
9.3 静态电压稳定极限估算
设某受端负荷由 $U_g = 1.0$ pu 的电源经总电抗 $X_\Sigma = 0.32$ 供电，负荷功率因数 $\cos\varphi = 0.95$（$\sin\varphi \approx 0.312$）。
最大可送有功为
\[P_{L,\max} = \frac{1.0^2}{2 \times 0.32} \cdot \frac{0.95}{1 + 0.312} \approx 1.131\text{ pu} = 113.1\text{ MW（基准 100 MVA）}\]达到极限时的受端最低电压为 $V_{\min} = 1/\sqrt{2 + 2 \times 0.312} \approx 0.617$ pu。这说明在功率因数为 0.95 的情况下，受端电压一旦跌到 0.62 pu 左右，系统实际上已经非常接近静态电压崩溃边缘。
9.4 母线投切电容器时的电压变化估算
设某 220 kV 母线看向系统的短路容量约为 $S_{sc} = 5000$ MVA，现投入一组并联电容器 $\Delta Q_c = 60$ Mvar。
\[\Delta U\,(\text{pu}) \approx \frac{60}{5000} = 0.012\]母线电压约上升 1.2%，对应 220 kV 线电压变化约为 $\Delta U_{\text{kV}} \approx 2.64$ kV。切除同一组电容器则电压变化方向相反。
9.5 500 kV 长线路的无功估算
设某 500 kV 长线路自然功率 $P_N = 1000$ MW，单位长度充电功率 $Q_N = 1.2$ Mvar/km，线路长度 $l = 200$ km，充电功率基数为 $Q_N l = 240$ Mvar。
轻载（700 MW）时：$\Delta Q_L = [(0.7)^2 - 1] \times 240 = -122.4$ Mvar（发无功）。重载（1400 MW）时：$\Delta Q_L = [(1.4)^2 - 1] \times 240 = +230.4$ Mvar（吸无功）。判断的分界点正是自然功率 1000 MW。这一结论与 9.4 算例联合使用，即可判断该线路在不同潮流方式下是否需要投入或切除无功补偿装置。
9.6 暂态稳定冲击法快估
设某送端系统故障期间近似加速功率 $\Delta P_a = 0.9$ pu，故障切除时间 $\Delta t = 0.12$ s，故障后网络最大传输功率 $P_{\max}^{\text{post}} = 1.65$ pu，$T_j = 9$ s。
故障后稳定平衡角：$\delta_s = \arcsin(0.9/1.65) \approx 33.1°$，$K_s^{\text{post}} = 1.65\cos33.1° \approx 1.383$。
故障后振荡频率：$f_d = (1/2\pi)\sqrt{314 \times 1.383/9} \approx 1.106$ Hz。
单边振幅：$D_p = 2\pi \times 1.106 \times 0.9 \times 0.12 \approx 0.750$ pu。
暂稳极限估算值：$P_{st} \approx 1.65 - 0.750 = 0.900$ pu。
进一步用临界切除时间公式做第二层校核（取 $\delta_0 = 0.578$ rad）：
\[\delta_{cr} = \arccos\!\left[\frac{0.9}{1.65}(\pi - 2 \times 0.578) - \cos33.1°\right] \approx 75.9° = 1.325\text{ rad}\] \[t_{cr} = \sqrt{\frac{2 \times 9 \times (1.325 - 0.578)}{314 \times 0.9}} \approx 0.218\text{ s}\]与实际切除时间 $\Delta t = 0.12$ s 相比，裕量约为 45%。等面积法对同一场景的精确校核给出暂稳极限约为 1.239 pu，冲击法偏保守约 27%——偏保守的来源在于，冲击法用线性化的简谐振荡振幅近似代替了非线性功角曲线上的等面积积分。两种快估方法合用，形成互补的双层结构：冲击法给出”第一摆摆幅有多危险”的直观量，临界切除时间公式给出”继电保护动作时限还差多远”的工程判据。
10. 近似公式的边界与误差来源
近似公式之所以有用，恰恰是因为它主动丢掉了很多细节。知道它在哪些地方会失效，和知道它本身一样重要。
频率动态：单机等值模型适用于事故后的前若干秒到十几秒，尤其适合分析频率下降率、频率最低点和一次调频的主导作用。时间再拉长，就会逐渐受到锅炉、AGC、联络线控制、负荷恢复特性、直流控制等中长期环节影响，误差会明显放大。
机电振荡频率：公式给出的是模态频率的主导近似，并不能直接给出阻尼比和振幅。若振荡问题的核心是”为什么振幅放大了”“为什么互联系统模态阻尼变差了”，就必须进入特征值分析、参与因子分析和控制器参数分析。
静态电压稳定极限：公式本质上是二节点等值结果。若系统存在复杂无功补偿、分接头动作、电阻不可忽略、多节点电压耦合显著，则这个公式只能用来建立直觉，不宜替代潮流—连续潮流—QV/PV 曲线分析。
线路自然功率：最适合超高压长线路。对中短线路、含串补、含显著电阻或存在复杂并联补偿的线路，直接使用自然功率公式会明显粗糙。
暂稳公式：冲击法是最粗的近似，适合快速判断第一摆摆幅和传输水平是否处于危险区域；临界切除时间公式稍严一些，更直接回答继电保护时限是否充裕的问题。两者都不能代替正式的暂态稳定仿真，但合用可以大幅减少”不必要的详细仿真次数”，同时避免对明显危险场景的漏判。
11. 结语
本文从两个层次建立了对电力系统的近似认知。
第一个层次是近似常数：单根导线载流量锁定在 400—1000 A，波速接近光速，正序电抗在 0.2—0.4 Ω/km，变压器短路阻抗百分比在 4%—15%，发电机惯性时间常数在 2—6 s。这些常数的稳定性并不是巧合，而是物理定律和材料极限共同作用的结果。
第二个层次是近似公式：频率动态的二阶模型、机电振荡频率公式、静态电压稳定极限、母线电压快算与电容器投切估算、线路自然功率与无功损耗，以及暂态稳定的冲击法与临界切除时间估算。这些公式的价值不在于替代精确模型，而在于把仿真结果背后的物理机制暴露出来，让人知道系统的动态大致会朝哪个方向走、关键参数的敏感方向是什么。
对真正有经验的分析人员，往往先用近似公式判断量级、方向和敏感参数，再用详细模型确认边界、修正偏差和落地控制。完全离开近似公式而只依赖数值仿真，反而很容易失去对系统的直观感觉。
参考文献
[1] 周荣光. 电力系统理论精析[M]. 电子工业出版社, 2014.
[2] 王梅义. 大电网事故分析与技术应用[M]. 中国电力出版社, 2008.
[3] Kundur P. 电力系统稳定与控制[M]. 中国电力出版社, 2002.
[4] 电力规划设计总院. 电力系统设计手册[M]. 中国电力出版社.