DFIG低电压穿越保护与Simulink仿真实践-嵌云网-嵌入式AI开发资源站

DFIG低电压穿越保护与Simulink仿真实践

玫瑰好吃

1. 项目概述：DFIG低电压穿越保护的必要性

双馈感应发电机（DFIG）作为现代风力发电的主流机型，其最大技术挑战在于应对电网电压骤降时的稳定运行问题。当电网发生短路故障导致电压突然跌落至额定值的20%-80%时，传统DFIG系统会在毫秒级时间内因转子侧过电流而触发保护脱网。根据中国GB/T 19963-2011《风电场接入电力系统技术规定》，风电机组必须在不脱网的情况下持续运行625毫秒（对20%电压跌落工况），这对控制策略提出了严苛要求。

我在参与某2MW风电场改造项目时，曾遇到因LVRT（Low Voltage Ride Through）性能不达标导致全场脱网的事故。事后分析发现，转子侧变流器的瞬时过电流达到额定值的3.2倍，而原厂设计的Crowbar保护存在150ms的动作延迟。这个惨痛教训促使我深入研究基于Simulink的LVRT仿真方法，下面分享的建模过程已在实际项目中得到验证。

2. DFIG系统架构与LVRT原理

2.1 正常运行时的能量流动路径

典型DFIG系统的拓扑结构包含三个关键部分：

机械侧：风力机通过齿轮箱驱动发电机转子（转差率通常为±30%）
电气侧：定子直接连接电网，转子通过背靠背变流器实现有功/无功解耦控制
控制层：采用矢量控制策略，通过调节转子电流实现最大功率点跟踪（MPPT）

在稳态运行时，转子侧变流器（RSC）仅需处理转差功率（约30%额定功率），这是DFIG相比全功率变流系统的核心优势。

2.2 电压跌落时的暂态过程

当电网发生三相短路时，定子电压突然跌落会引发两个致命问题：

定子磁链暂态：根据法拉第定律，磁链不能突变，会在定子绕组中产生直流分量磁链，其衰减时间常数τs=Ls/Rs（约100-300ms）
转子感应过电压：旋转的转子绕组切割直流磁链，感应出滑差频率的暂态电动势，其幅值可达：
[
E_{r_transient} = \frac{L_m}{L_s} \cdot (1-s) \cdot \Delta V_s
]
其中Lm为互感，Ls为定子电感，s为转差率，ΔVs为电压跌落幅度

关键提示：当电压跌落80%时，转子开路电压可能瞬间达到额定值的3-5倍，远超IGBT模块的耐压水平（通常为1.5倍裕量）

3. Crowbar保护策略设计要点

3.1 硬件电路实现方案

常见的Crowbar拓扑有三种：

被动式：采用压敏电阻（MOV）并联在转子侧，成本低但响应不可控
主动式：晶闸管（SCR）串联电阻，通过门极触发实现快速投入（本文采用方案）
混合式：IGBT与二极管反并联，支持双向导通和快速切除

在Simulink中建模时，需要特别注意SCR的关断特性——只有当电流过零时才能自然关断，这直接影响保护退出时机。

3.2 投切逻辑设计规范

基于国标要求，我们的Stateflow状态机包含以下判定条件：

投入条件（OR逻辑）：

转子电流瞬时值 > 1.8 pu（持续5ms）
直流母线电压 > 1.25 pu
检测到电网电压跌落（dip < 0.8 pu）

切除条件（AND逻辑）：

电网电压恢复至0.9 pu以上
转子电流衰减至1.2 pu以下
Crowbar持续导通时间 ≥ 60ms（防止误切除）

mermaid复制stateDiagram
    [*] --> Normal
    Normal --> CrowbarActive: 触发投入条件
    CrowbarActive --> Normal: 满足切除条件
    CrowbarActive --> Fault: 持续500ms未切除
    Fault --> [*]

3.3 关键参数计算示例

以1.5MW机组为例（基准电压690V）：

Crowbar电阻选择：
- 目标将转子电流限制在1.5pu内
- 已知转子额定电阻Rr=0.01pu，暂态电抗X'=0.15pu
- 所需总阻抗Z=Vs_max/I_lim=3.0/1.5=2.0pu
- 计算得R_cb=√(2²-0.15²)-0.01≈1.98pu
- 实际值：1.98×(690²/1.5MΩ)=0.63Ω
动作时间验证：
- SCR触发延迟<100μs
- 电流检测滤波窗口=1/4周期（5ms@50Hz）
- 总响应时间<6ms，满足国标<10ms要求

4. Simulink建模全流程

4.1 主电路搭建技巧

使用Simscape Electrical库时，建议采用以下配置：

matlab复制% DFIG参数设置示例
DFIG.Rs = 0.023;       % 定子电阻(pu)
DFIG.Lls = 0.18;       % 定子漏感(pu)  
DFIG.Lm = 2.9;         % 互感(pu)
DFIG.J = 0.1;          % 转动惯量(kg·m²)

% 电网故障设置
Fault.StartTime = 1.0; % 故障起始时间(s)
Fault.Rf = 0.05;       % 短路阻抗(pu)

实用技巧：

为加快仿真速度，可将变流器简化为受控电压源
使用"Three-Phase VI Measurement"模块获取正序电压幅值
在转子回路添加"Three-Phase Breaker"模块模拟Crowbar动作

4.2 Stateflow逻辑实现

状态机的关键代码片段：

matlab复制function [cb_cmd] = CrowbarLogic(V_grid, I_rotor, t_cb)
    persistent state;
    
    % 投入条件判断
    if (max(abs(I_rotor)) > 1.8 || V_grid < 0.8) && isempty(state)
        state = "Active";
        t_cb = 0;
    end
    
    % 切除条件判断
    if strcmp(state, "Active")
        t_cb = t_cb + Ts;
        if V_grid > 0.9 && max(abs(I_rotor)) < 1.2 && t_cb >= 0.06
            state = [];
        elseif t_cb > 0.5
            state = "Fault";
        end
    end
    
    cb_cmd = strcmp(state, "Active");
end

4.3 典型故障场景测试

建议分阶段验证：

轻度跌落测试（0.7pu, 100ms）
- 检查Crowbar是否不误动
- 观察转子电流是否通过变流器调节抑制
深度跌落测试（0.2pu, 625ms）
- 验证Crowbar在5ms内动作
- 监测直流母线电压波动<1.3pu
重复跌落测试（0.5pu, 200ms间隔）
- 检查保护系统能否多次可靠动作
- 评估机组暂态转矩冲击(<2.5pu)

5. 仿真结果分析

5.1 关键波形对比

在80%电压跌落工况下（1.0-1.625s）：

参数	无Crowbar	有Crowbar
转子峰值电流	3.2 pu	1.6 pu
直流母线电压	1.45 pu	1.18 pu
恢复时间	脱网	350 ms
暂态转矩	2.8 pu	1.9 pu

实测发现：Crowbar电阻值增加20%可进一步降低转矩脉动，但会延长恢复时间，需根据电网要求权衡

5.2 工程优化建议

动态电阻调整：根据跌落深度自适应改变R_cb值
混合保护策略：在Crowbar投入初期配合桨距角快速调节
故障预测：利用PMU数据提前10ms预判电压跌落

6. 进阶研究方向

6.1 与储能系统协同控制

在直流母线侧增加超级电容储能：

故障期间提供短时功率支撑（100-200ms）
抑制Crowbar动作期间的直流电压波动
仿真时需添加双向DC-DC变换器模型

6.2 多机协调LVRT

对于风电场集群：

主控单元实时通信各机组状态
优先让转速较高的机组投入Crowbar
通过SVG动态补偿无功缺额

我在实际项目中验证过，这种策略可将全场脱网风险降低60%以上。