4. 小信号稳定性与特征值分析
低频振荡的本质，是系统在某个稳态运行点附近受到小扰动后，状态量之间通过发电机转子运动、励磁调节和网络功率交换形成的弱阻尼能量往复。要判断这种往复是衰减、等幅还是发散，最直接的方法是把系统在运行点处线性化，得到状态空间模型，再用特征值、阻尼比、振荡频率和参与因子来定量描述。本节围绕这条主线，说明建模、线性化、特征值计算和结果判读的完整过程。
一、状态空间建模：把动态系统写成统一形式
小信号分析的第一步是建立状态空间模型。对电力系统而言，状态变量通常取发电机转子角 δ、转速偏差 Δω、暂态电势 E′q、励磁系统状态量等；代数变量包括节点电压幅值和相角。设状态向量为 x，代数向量为 y，输入为 u，则系统可写成微分—代数方程：
ẋ = f(x, y, u)
0 = g(x, y, u)其中 f 描述转子运动、励磁和调速等动态，g 描述网络功率平衡。稳态运行点 (x₀, y₀) 满足 f(x₀, y₀, u₀)=0 且 g(x₀, y₀, u₀)=0。只有先求出潮流解和初值，线性化才有意义。
二、线性化：在运行点附近取一阶近似
在运行点处对 f 和 g 求偏导，得到增量方程：
Δẋ = A Δx + B Δy + E Δu
0 = C Δx + D Δy + F Δu若 D 可逆，消去 Δy，得到标准状态方程 Δẋ = Asys Δx + Bsys Δu，其中 Asys = A − B D⁻¹ C。Asys 就是小信号分析的核心矩阵。它的维数等于状态变量个数，通常每台发电机贡献 2～4 阶，系统规模越大，矩阵越庞大。
三、特征值、阻尼比与振荡频率
求解 Asys 的特征方程 det(λI − Asys) = 0，得到特征值 λ。对于一对共轭复根 λ = σ ± jω，可进一步计算：
- 阻尼比 ζ = −σ / √(σ² + ω²)，ζ 越大衰减越快；
- 振荡频率 f = ω / (2π)，单位 Hz；
- 衰减时间常数 τ = 1 / |σ|，反映振荡包络衰减到 1/e 所需时间。
工程上通常关注 0.1～2.5 Hz 的机电振荡模式。若某模式 ζ 小于 3%～5%，就认为阻尼偏弱，存在低频振荡风险；若 σ 为正，则系统小信号不稳定，振荡会持续发散。
四、参与因子：定位振荡中的关键机组
特征值只说明“有没有振荡”，参与因子则回答“谁在振荡”。对第 i 个特征值 λi，其左、右特征向量分别为 vi 和 wi，参与因子定义为：
p_ki = |v_ki · w_ki| / (v_iᵀ w_i)pki 越大，说明第 k 个状态量在该模式中参与越深。若某台机组的 Δω 参与因子最大，它就是该振荡模式的主导机组；若多台机组参与因子接近，则表现为区域间振荡。参与因子还能帮助判断励磁、调速等控制环节对阻尼的贡献。
五、计算示例与结果判读
假设某两机系统在某一运行点下，线性化后得到一对共轭特征值 λ = −0.35 ± j4.82。按公式计算：
ζ = 0.35 / √(0.35² + 4.82²) ≈ 0.0724 → 7.24%
f = 4.82 / (2π) ≈ 0.767 Hz
τ = 1 / 0.35 ≈ 2.86 s该模式频率约 0.77 Hz，属于典型区域间低频振荡；阻尼比 7.24% 高于 5% 的警戒线，说明在当前运行方式下阻尼尚可，但若负荷加重或联络线潮流增大，ζ 可能继续下降。此时应结合参与因子查看主导机组，再评估是否投入电力系统稳定器（PSS）或调整励磁参数。
若另一模式特征值为 λ = +0.12 ± j3.60，则 σ > 0，系统小信号不稳定，振荡会发散。此时必须优先处理：检查励磁系统增益是否过高、PSS 是否退出、网络拓扑是否削弱了同步力矩，再重新计算特征值验证。
六、工程应用要点
- 运行点决定结论：同一系统在不同负荷、不同开机方式下，特征值可能完全不同，必须针对关键运行方式分别计算。
- 关注弱阻尼模式：优先筛选 ζ 最小的若干模式，而不是只看绝对频率。
- 参与因子指导措施：参与因子大的机组加装 PSS 效果最明显；参与因子分散的模式则需要多机协调。
- 线性化与非线性仿真结合：特征值给出稳定性和阻尼的定量指标，时域仿真用于验证大扰动下的实际表现。