双馈风机高电压穿越(HVRT)控制策略与工程实践

1. 风电系统高电压穿越（HVRT）的必要性

双馈感应发电机（DFIG）作为现代风电系统的主流机型，其高电压穿越能力已成为并网硬性要求。当电网电压骤升至1.3倍额定值时，传统控制策略会导致机组在100毫秒内脱网——这相当于让整个风电场在电网最需要支撑时"撂挑子"。国标GB/T 19963-2021明确规定：机组必须能在1.3 p.u.电压下坚持500毫秒以上，同时提供至少0.2 p.u.的感性无功支撑。

实际工程案例：某200MW风场曾因缺乏HVRT功能，在电网暂态过电压时集体脱网，导致区域电网崩溃，直接经济损失超3000万元。

2. DFIG高电压故障的物理本质

2.1 电磁耦合效应

当电网电压Vs突然升高时，定子磁链ψs=∫Vsdt会产生阶跃增长。根据转子电压方程vr=-Lm/Ls·dψs/dt-jωsψs，将同时产生：

• 瞬态分量（dψs/dt项）：引发高达2kV的电压尖峰
• 稳态分量（ωsψs项）：持续抬升转子电压基准

2.2 能量失衡危机

定子输出功率Ps≈Vs²/Rs随电压平方增长，而风轮机械功率Pm受风速限制无法突变。功率差额ΔP=Ps-Pm会通过变流器涌入直流母线，其动态过程为：

C·dVdc/dt = Pgrid - Pgen = ΔP

以1.5MW机组为例，1.3p.u.电压下ΔP可达200kW，若不及时处理，10mF的直流电容电压将在10ms内飙升超过1500V。

3. 三层防御体系构建

3.1 主动控制层设计

转子侧变流器(RSC)控制

matlab复制function [idr_ref, iqr_ref] = RSC_HVRT_Control(Vs_pu)
    % 参数设定
    k1 = 0.5;  % 降载系数
    k2 = 0.3;  % 无功补偿系数
    
    % 动态指令生成
    iqr_ref = -k1 * (Vs_pu - 1.0);  % 负iq减少有功输出
    idr_ref = +k2 * (Vs_pu - 1.0);  % 正id提供感性无功
end

关键参数整定原则：

• k1取值需使ΔP≤Chopper最大耗散功率
• k2取值需满足≥0.2p.u.无功要求

网侧变流器(GSC)控制

采用双闭环结构：

1. 外环电压控制：PI输出idg_ref
2. 外环无功控制：强制iqg_ref=+0.2p.u.
3. 内环电流控制：带宽≥500Hz

3.2 被动保护层实现

3.3 状态机协调逻辑

matlab复制function [mode] = HVRT_StateMachine(Vs_pu, Vdc, Vr_peak)
    persistent timer_start
    
    % 事件检测
    if Vs_pu >= 1.1 && isempty(timer_start)
        timer_start = tic;
        mode = 'HVRT_Mode';
    elseif Vs_pu < 1.05 || toc(timer_start) > 0.5
        mode = 'Normal_Mode';
        timer_start = [];
    else
        mode = 'HVRT_Mode';
    end
end

4. Simulink建模实操指南

4.1 主电路搭建步骤

1. 从Simscape Electrical库拖拽：

   • DFIG模块（设置Rs=2.6mΩ, Lls=0.087mH）
   • 背靠背变流器（IGBT参数Ron=1mΩ）
   • 直流母线电容（C=10mF）

2. 保护电路实现：

   matlab复制% Crowbar子系统
   if Vr_peak > 1200
       thyristor_gate = 1;
   else
       thyristor_gate = 0;
   end
   

4.2 关键仿真设置

• 求解器：ode23tb（适合电力电子系统）
• 步长：1μs（必须捕捉电压尖峰）
• 仿真时长：3s（包含扰动前后各1s）

5. 工程调试避坑指南

5.1 Crowbar误动作处理

现象：转子电压未达阈值却频繁触发
解决方法：

• 检查电压测量滤波时间常数（推荐10-50μs）
• 调整触发延时（建议5ms防抖）

5.2 GSC振荡问题

现象：无功电流出现10-20Hz波动
优化方案：

• 增加电流环阻尼比至0.8
• 在iqg_ref路径添加一阶惯性环节（T=0.01s）

5.3 参数敏感度分析

6. 进阶优化方向

6.1 预测控制算法

基于Luenberger观测器提前1-2个周期预判电压突变：

matlab复制function Vs_pred = Voltage_Predictor(Vs_history)
    % 二阶外推算法
    coeff = [1.5 -0.5];  % 预测系数
    Vs_pred = coeff * Vs_history(end-1:end);
end

6.2 风场协同控制

各机组HVRT策略参数优化：

• 中央控制器分配无功指令
• 考虑集电线阻抗影响

实测数据表明，协同控制可使风场整体无功响应速度提升40%