同步发电机与可再生能源并网的暂态稳定性分析

1. 同步发电机与可再生能源并网系统的暂态稳定性挑战

在新型电力系统建设背景下，同步发电机(Synchronous Generator, SG)与电网跟随型可再生能源(Grid-Following Renewable Energy Sources, GFLR)的协同运行已成为典型场景。这种混合发电模式虽然能够兼顾系统惯量支撑与清洁能源消纳，但也带来了前所未有的暂态稳定挑战。

传统电力系统中，同步发电机的暂态稳定性分析主要关注"首摆稳定性"(First-Swing Stability)，即故障清除后第一个振荡周期内系统保持同步运行的能力。然而，当系统中存在大量GFLR时，其特有的低电压穿越(Low Voltage Ride Through, LVRT)控制特性会与SG的转子动态产生复杂的交互作用，可能引发以下新现象：

1. 周期性LVRT切换：SG转子角度(δ)的摆动会导致并网点电压(Uw)波动，当Uw低于LVRT激活阈值(Uin)时，GFLR会进入LVRT模式；而当电压恢复超过退出阈值(Uout)时，又返回正常工作模式。这种周期性切换会持续影响SG的转子动态。

2. 能量累积效应：不同LVRT控制策略（圆形限幅与矩形限幅）会以不同方式改变GFLR的输出功率(Pw)，进而通过功率耦合项(αEsI_w^ds)影响SG的加速/减速能量平衡。某些情况下，这种交互会导致每个振荡周期都产生净能量输入，最终引发多摆失稳(Multi-Swing Instability)。

3. 控制模式依赖：研究发现，圆形限幅(Circular Limiting)策略由于允许更大的有功电流输出，更容易恶化SG的首摆稳定性；而矩形限幅(Rectangular Limiting)配合缓慢恢复策略，则可能因持续的功率扰动导致多摆失稳。

2. 系统动态交互机制解析

2.1 电压-角度耦合机理

SG转子角度(δ)与GFLR并网点电压(Uw)之间存在紧密的耦合关系，可通过以下电压方程描述：

U_w = √[(αE_s cosδ + (1-α)U_g)^2 + (αE_s sinδ)^2] - βI_w sinη

式中关键参数：

• α = Y_s/(Y_s + Y_g)：表征SG与电网连接的紧密程度
• β = 1/(Y_s + Y_g) + 1/Y_w：系统导纳参数
• η = φ_w - θ：GFLR电流相位与PLL角度差

当δ增大时，U_w会相应降低，这种反比关系是LVRT周期性触发的物理基础。通过求解U_w=U_in和U_w=U_out对应的临界角度δ_in和δ_out，可以预测LVRT的激活/退出时机。

2.2 能量动态分析框架

采用扩展的暂态能量函数(Transient Energy Function, TEF)方法分析系统稳定性。定义系统的Lyapunov函数为：

V(δ,Δω) = 1/2 T_J Δω^2 + ∫_{δ_s}^δ (P_e - P_m - ΔP_w) dδ

其中能量变化率满足： dV/dt = -DΔω^2 + ΔP_w Δω

ΔP_w的正负和时变特性决定了GFLR对SG的阻尼效果：

• 当ΔP_w < 0且Δω > 0（相位I、IV）：提供正阻尼
• 当ΔP_w < 0且Δω < 0（相位II、III）：产生负阻尼

这种阻尼特性的不对称性正是多摆失稳的根源。

3. LVRT控制策略的稳定性影响

3.1 圆形限幅策略分析

圆形限幅的电流指令满足： i_d^{ref} = min{P^{ref}/U_w, √(I_max^2 - i_q^2)}

主要特性：

1. 优先保证无功电流支持(i_q = K_q(0.9 - U_w))
2. 最大化有功电流输出以减少功率损失
3. 在δ上升和下降阶段产生的ΔP_w近似对称，净能量输入ΔV_wa≈0

稳定性风险：

• 当K_q较小时，为维持总有功输出，i_d会趋近I_max
• 导致ΔP_w > 0，在Δω > 0阶段产生负阻尼
• 特别不利于首摆稳定性

3.2 矩形限幅策略分析

矩形限幅直接约束有功电流： i_d^{ref} = i_d,lim^{ref}

配合的恢复策略通常为： i_d^* = v_d t （斜坡恢复）

稳定性特点：

1. 严格限制有功电流，首摆阶段ΔP_w较小
2. 缓慢恢复阶段产生时变ΔP_w(t)
3. 通过积分计算可得每个周期净能量输入： ΔV_wb = (C_2 - C_1)(δ_out - δ_N) > 0

典型案例：

• 某实际系统中，当恢复速率v_d=2.5A/s时
• 连续3个振荡周期累积能量达70.57(kA·rad²/s)
• 最终导致SG在第三摆失去同步

4. 稳定性判据与边界计算

4.1 首摆稳定判据

保守的首摆临界能量： V_f^{max} = ∫_{δ_s}^{π-δ_s} [P_m - P_e + αE_s I_max] dδ

应用条件：

• 采用圆形限幅且K_q较小
• 计算故障清除时刻能量V_c
• 若V_c < V_f^{max}则首摆稳定

4.2 多摆稳定判据

考虑最恶劣场景的临界能量： V_m^{max} = V_m1^{max} - ΔV_wr

其中： V_m1^{max} = ∫_{δ_s}^{δ_in} [P_m - P_e] dδ ΔV_wr = αE_s I_max (δ_M - δ_N)

工程应用方法：

1. 计算故障清除能量V_c
2. 若V_c < min(V_f^{max}, V_m^{max})则全周期稳定
3. 否则需采取稳定控制措施

5. 反馈线性化控制器设计

5.1 控制原理

基于反馈线性化理论，设计附加电流控制器： I_add = [P_m - P_e^* - K_1Δω^* - K_2(δ^* - δ_d)]/(αE_s)

实现：

1. 精确对消非线性项(P_m - P_e^*)
2. 引入线性状态反馈提供可调阻尼
3. 通过Lyapunov函数证明全局渐近稳定

5.2 工程实现方案

实际控制律转换为： Δi_d^* = k_A [T_J s/(αE_s)] [k_1(ω_p - ω_g) + (θ - θ_s)]

实施要点：

1. 在LVRT退出后激活，持续3-5秒
2. k_1通常取2-3，k_A取1.5-2
3. 需限制输出范围避免饱和

5.3 控制效果验证

CHIL测试结果表明：

1. 首摆不稳定案例：
   • 无控制时E_-=0.75(kA·rad²/s)
   • 控制后E_+=0，成功抑制失稳
2. 多摆不稳定案例：
   • 无控制时周期性能量累积70.57
   • 控制后转为能量耗散-35.83
3. 动态响应时间<100ms

6. 工程实践建议

基于研究成果，提出以下工程实施建议：

1. LVRT参数整定原则：

   • 避免U_out - U_in < 0.02p.u.（防止高频切换）
   • 不建议单纯通过调整阈值解决稳定性问题

2. 控制器部署方案：

   • 优先在GFLR集群的中央控制器中实现
   • 需与SG的励磁控制协调（时间常数匹配）

3. 系统级保护配置：

   • 设置多摆失稳预测模块
   • 与现有失步保护配合（延时差异化）

4. 对GFM并网的启示：

   • 虚拟阻尼可抑制多摆振荡
   • 需注意故障后电流持续受限场景

7. 实测数据与案例分析

某实际风电汇集区域（SG容量4×630MVA，风电1.5GW）的仿真验证：

1. 圆形限幅场景：

   • K_q=1.5，故障清除时间120ms
   • 首摆临界能量1.76(p.u.)，实际1.82(p.u.)
   • 仿真结果：首摆失稳

2. 矩形限幅+慢恢复场景：

   • v_d=2A/s，U_out=0.95p.u.
   • 第三摆出现增幅振荡
   • 控制器投入后2.5秒平息

3. 控制器性能指标：

   • 响应时间：85ms
   • 超调量：<5%
   • 对风电机组出力影响：<3%

8. 技术延伸与未来方向

本研究的理论框架还可扩展到以下领域：

1. 混合发电系统：

   • SG与虚拟同步机(VSG)的交互
   • 多类型GFLR的集群效应

2. 新型控制策略：

   • 基于深度强化学习的自适应控制
   • 分布式协同控制架构

3. 稳定性评估工具开发：

   • 全周期暂态稳定快速评估算法
   • 数字孪生平台中的实时仿真

未来工作将重点关注：

1. 计及PLL动态的高阶模型
2. 故障后SG与GFLR的协同恢复策略
3. 面向高比例可再生能源的稳定标准修订

通过这项研究，我们建立了SG与GFLR交互稳定性的完整分析框架，提出的控制方法已在我国多个可再生能源基地得到应用验证。实践证明，这种基于非线性控制的解决方案比传统参数调整方法具有更好的适应性和鲁棒性。