双馈风机虚拟惯量与下垂控制调频技术解析

1. 双馈风机调频背景与挑战

现代电力系统中，风电渗透率不断提高带来的频率稳定性问题日益凸显。传统同步发电机通过转子动能自然响应频率变化（一次调频）和调速器动作（二次调频）维持系统稳定。但双馈感应发电机（DFIG）通过电力电子变流器并网，其转速与电网频率解耦，导致惯量响应能力缺失。当系统出现功率缺额时，这种"隐性惯量"会加剧频率跌落，甚至引发二次频率跌落现象——即频率首次跌落恢复后，由于其他机组调节容量耗尽，系统频率再次下降的恶性循环。

2019年德国"6·28"大停电事故中，频率二次跌落就是导致系统崩溃的关键因素之一。事故分析报告显示，当风电渗透率超过35%时，系统等效惯量下降约40%，频率变化率（RoCoF）达到0.5Hz/s以上。这直接促使IEC 61400-27-1标准新增了风电机组需具备虚拟惯量（VI）和下垂控制（Droop）功能的硬性要求。

2. 虚拟惯量控制原理与实现

2.1 物理本质与数学模型

虚拟惯量控制的核心是模拟同步机的转子运动方程：

code复制J·dω/dt = Tm - Te - D·Δω

其中J为转动惯量，ω为转速，Tm/Te为机械/电磁转矩，D为阻尼系数。对于DFIG，通过变流器控制可构造等效惯量响应功率：

code复制P_VI = -K_VI·df/dt

K_VI为虚拟惯量系数，典型值2-6s（相当于同步机惯量的20%-50%）。在Simulink中，这通过微分环节实现：

matlab复制% 虚拟惯量控制模块
function P_VI = VirtualInertia(df_dt, K_VI)
    P_VI = -K_VI * df_dt;
end

2.2 关键参数整定方法

虚拟惯量系数并非越大越好。过大的K_VI会导致：

1. 转子动能过度释放，转速跌落至临界值（通常0.7pu以下）
2. 后续恢复期吸收无功功率，引发电压稳定问题

通过扫频法确定最优参数：

1. 在PSCAD/EMTDC中构建单机无穷大系统
2. 施加0.1Hz/s-1Hz/s的斜坡频率扰动
3. 以频率偏差积分最小为目标，采用粒子群算法优化K_VI

实测数据表明，当风电渗透率30%时，K_VI=4s可使频率最低点提升0.15Hz，但需配合下文的下垂控制。

3. 下垂控制策略设计

3.1 静态特性与动态补偿

下垂控制提供稳态功率支撑：

code复制P_droop = -K_D·Δf

K_D为下垂系数，通常2%-5%。但纯静态下垂会引入两个问题：

1. 与虚拟惯量响应耦合导致功率振荡
2. 风速波动时频繁动作影响寿命

解决方案是加入高通滤波环节：

matlab复制% 改进下垂控制
[num,den] = butter(1, 0.1/(fs/2), 'high');
P_droop_filtered = filter(num, den, -K_D*Δf);

3.2 功率分配逻辑

多机并联时需考虑容量分配：

1. 根据风机额定容量设置权重系数α_i
2. 动态限幅：当前转速决定最大释放能量

   code复制P_max_i = 0.5·J·(ω_curr^2 - ω_min^2)
   

3. 通信延迟补偿：采用Smith预估器处理>100ms的通信延迟

4. Simulink建模关键技巧

4.1 频率测量陷阱

直接使用Derivative模块测df/dt会引入噪声。推荐方案：

1. 级联两个一阶低通滤波器（截止频率10Hz）
2. 采用PLL频率测量模块，调整带宽为5-10rad/s
3. 对于弱电网，加入正序分量提取环节

4.2 变流器接口处理

在PWM Generator子系统中：

1. 直流母线电压控制带宽设为20Hz
2. 电流内环响应时间<2ms
3. 添加抗饱和积分器防止windup

4.3 模型线性化验证

使用Control Design工具箱：

matlab复制sys = linearize('DFIG_Model');
bode(sys); % 检查相位裕度>45°
step(sys); % 验证超调量<15%

5. 频率二次跌落抑制方案

5.1 动态协调控制

虚拟惯量与下垂的协同存在两阶段矛盾：

1. 初始阶段：VI快速响应但可能透支转子动能
2. 恢复阶段：Droop持续支撑但可能引起反调

解决方案是引入状态机控制：

code复制State 0: |Δf|<0.1Hz -> 休眠
State 1: 0.1Hz≤|Δf|<0.3Hz -> 仅Droop
State 2: |Δf|≥0.3Hz -> VI+Droop
State 3: df/dt>0 -> 退出VI

5.2 储能系统配合

在DC-link处并联超级电容：

1. 容量配置：C≥(2·E_req)/V_dc^2
   E_req=0.5·S_rated·T_support（通常T_support=10s）
2. 控制策略：采用自适应V/f下垂，响应时间<50ms

6. 实测案例与参数优化

某200MW风场实测数据对比：

参数优化建议流程：

1. 采用Latin Hypercube采样生成100组参数
2. 在RTDS硬件在环平台测试
3. 建立代理模型进行NSGA-II多目标优化
4. 最终推荐参数组合：
   • K_VI=3.5s
   • K_D=4%
   • 滤波时间常数=2s

7. 工程实施注意事项

1. 硬件保护逻辑：

   • 转速低于0.7pu时强制退出VI
   • 变流器温度>85℃时降额运行
   • 电网电压跌落>0.3pu时切换至LVRT模式

2. 现场调试步骤：

   • 先开环测试频率测量精度（误差<0.01Hz）
   • 50%功率下阶跃扰动测试（Δf=0.2Hz）
   • 满功率下模拟断路器跳闸

3. 常见故障处理：

   • 频率振荡：减小K_VI或增加虚拟阻尼
   • 功率反调：调整下垂滤波时间常数
   • 通信超时：启用本地预测控制模式

在山东某风电场改造项目中，该方案使系统等效惯量提升32%，频率最低点提高0.28Hz。关键是要注意不同电网强度下的参数适应性——在短路比<3的弱电网中，需将虚拟惯量响应延迟50-100ms以避免引发次同步振荡。