4. 新能源并网小干扰稳定分析

新能源场站经变流器并网后，系统动态特性与同步机主导的传统电网有本质差异：大量电力电子设备在宽频段内呈现低惯量、弱阻尼特征，小干扰稳定问题不再局限于低频机电振荡，而是扩展到次同步、超同步乃至高频段。小干扰稳定分析的目标，是在平衡点附近将非线性系统线性化，通过状态矩阵的特征信息判断系统对微小扰动的响应是收敛还是发散，并进一步定位"是谁在振荡、谁该负责、该在哪里加阻尼"。

一、从非线性模型到状态矩阵

设系统动态方程为 dx/dt = f(x, u)，在稳态平衡点 x₀ 处做一阶泰勒展开，得到线性化模型 Δẋ = A·Δx，其中 A = ∂f/∂x |ₓ₀ 即状态矩阵（雅可比矩阵）。A 的维数等于系统独立状态变量个数，新能源并网系统中通常包括：同步机转子角与转速、励磁与调速环节状态、双馈/直驱风机的锁相环（PLL）状态、直流母线电压、电流内环与外环积分状态、虚拟同步机功角与频率状态等。工程上 A 的阶数可达数百甚至上千阶，因此实际分析依赖特征值分解而非手工推导。

对 A 求特征值 λᵢ = σᵢ + jωᵢ：实部 σᵢ 决定该模式的衰减/发散速率，虚部 ωᵢ 对应振荡角频率 fᵢ = ωᵢ/(2π)。判断准则为：σᵢ < 0 稳定，σᵢ = 0 临界，σᵢ > 0 失稳。阻尼比 ζᵢ = −σᵢ / √(σᵢ² + ωᵢ²) 是工程上更直观的指标，通常要求 ζᵢ ≥ 5%（弱阻尼预警线），低于 2%~3% 即视为高风险模式。

二、参与因子：把"模式"与"设备"对应起来

特征值只告诉我们"系统里存在哪些振荡模式"，但不告诉我们"哪个模式由哪些设备主导"。参与因子矩阵 P 解决了这一映射问题。设 A 的右特征向量矩阵为 V（列向量 vᵢ），左特征向量矩阵为 W（行向量 wᵢᵀ），满足 WᵀV = I，则参与因子定义为：

p_ki = |w_ki · v_ik| / (Σ_k |w_ki · v_ik|)

其中：
  v_ik —— 第 i 个模式对第 k 个状态变量的可观性（右特征向量分量）
  w_ki —— 第 i 个模式对第 k 个状态变量的可控性（左特征向量分量）
  p_ki —— 第 k 个状态变量在第 i 个模式中的参与程度

参与因子越大，说明该状态变量与该模式的耦合越强。工程判读时，把同一台设备（如同一台风机）的多个状态变量参与因子求和，即可得到"设备—模式"参与度，从而回答"这个振荡模式主要由哪几台机组、哪个控制环节引起"。

三、振荡模式识别与分类

新能源并网系统的振荡模式按频段大致可分为三类：低频段（0.1~2.5 Hz）以同步机功角振荡、虚拟同步机功率振荡为主；次/超同步段（数 Hz 至数十 Hz）常见于双馈风机与串补电网的次同步谐振、PLL 与电流环耦合振荡；高频段（数百 Hz 至 kHz）多与变流器开关频率、LCL 滤波器谐振、多机并联阻抗交互相关。识别模式时，除频率与阻尼比外，还应结合参与因子分布、模式形状（各机组相位关系）以及是否随某参数（如 PLL 带宽、电流环增益、短路比 SCR）单调变化来综合判断。

小干扰稳定分析流程：从状态矩阵到弱阻尼环节定位 ① 线性化与特征值 Δẋ = A·Δx A = ∂f/∂x |ₓ₀（雅可比） λᵢ = σᵢ + jωᵢ fᵢ = ωᵢ / 2π ζᵢ = −σᵢ / √(σᵢ²+ωᵢ²) 判据：σᵢ<0 稳定；ζᵢ<5% 弱阻尼 ② 参与因子分析 p_ki = |w_ki · v_ik| / Σ|w_ki · v_ik| v_ik：可观性（右特征向量） w_ki：可控性（左特征向量） 按设备聚合： P_dev = Σ p_ki（同设备状态） 输出：模式 ↔ 设备/控制环节映射 ③ 模式识别与定位 低频 0.1~2.5 Hz：功角/VSG 次/超同步：PLL、电流环 高频：LCL、开关频率 辅助判据： · 参与因子分布 · 模式形状（相位关系） · 参数灵敏度 ∂ζ/∂K ④ 弱阻尼环节定位与改善措施 PLL 锁相环 带宽过高 → 弱阻尼 参与因子：θ_pll 高 措施：降低带宽 或加 PLL 阻尼补偿 电流内环 增益过高 → 高频振荡 参与因子：i_d, i_q 高 措施：降低 Kp/Ki 或加有源阻尼 直流母线电压环 带宽与 PLL 耦合 参与因子：U_dc 高 措施：错开带宽 加虚拟惯量/阻尼 电网强度（SCR） SCR 低 → 阻抗交互强 参与因子：多机耦合 措施：加 SVG/调相机 或降低外环带宽 参数调整后重新计算特征值，验证 ζ 是否提升至 ≥5%

四、弱阻尼环节定位的工程做法

定位弱阻尼环节，核心思路是"特征值—参与因子—参数灵敏度"三步联动。以某直驱风电场经弱电网并网为例（以下为教学示例参数，非实测）：额定容量 200 MW，由 40 台 5 MW 机组组成，并网点短路比 SCR = 2.0，PLL 带宽设为 50 Hz，电流内环 Kp = 0.8、Ki = 20。对状态矩阵求特征值，发现一对共轭特征值 λ = −0.35 ± j 62.8，对应振荡频率约 10 Hz，阻尼比 ζ ≈ 0.56%，属于典型弱阻尼次同步模式。

进一步计算参与因子，发现该模式下 PLL 的锁相角状态 θ_pll 与 q 轴电流状态 i_q 的参与因子分别为 0.31 和 0.27，远高于其他状态；按机组聚合后，40 台机组的 PLL 参与度合计占该模式的 68%。这说明该 10 Hz 振荡主要由 PLL 动态主导。做参数灵敏度扫描：将 PLL 带宽从 50 Hz 降至 20 Hz，该模式阻尼比从 0.56% 提升至 3.8%；再在 PLL 中引入阻尼补偿项，阻尼比可进一步提升至 6.2%，满足 ζ ≥ 5% 的工程要求。这一过程就是"识别模式 → 定位环节 → 调整参数 → 复算验证"的完整闭环。

工程判读要点：参与因子高不等于一定是"故障源"，还需结合模式形状判断相位关系。若多台机组在同一模式下参与因子均高且相位接近，说明是全场共性环节（如 PLL 带宽统一偏高）；若仅个别机组参与因子突出，则可能是该机组参数异常或运行工况偏离。

五、分析中的常见误区

  • 只看特征值实部，忽略阻尼比。实部绝对值小但频率高的模式，阻尼比可能极低，实际振荡衰减很慢，必须用 ζ 判断。
  • 参与因子只看单个状态，不做设备聚合。单状态参与因子高可能只是坐标选择问题，按设备聚合后才能对应到真实物理对象。
  • 忽略运行点变化。新能源出力、电网短路比、机组投退都会改变 A 矩阵，弱阻尼模式可能只在特定工况下出现，需做多工况扫描。
  • 把线性化结论直接外推到大扰动。小干扰稳定分析只适用于平衡点附近的小扰动，大扰动失稳需结合时域仿真判断。
← 上一章 📖 返回目录 下一章 →