双馈风机低压穿越技术与Simulink仿真实践

1. 风机短路故障分析与低压穿越技术概述

在风电装机容量占比不断提升的今天，电网对风电机组故障穿越能力的要求越来越高。当电网发生短路故障导致电压跌落时，传统风机会因保护机制而脱网，这在风电渗透率高的区域会引发连锁反应，严重威胁电网稳定。低压穿越(LVRT)技术就是让风机在电压跌落期间保持并网运行，同时提供无功支撑的关键解决方案。

双馈感应发电机(DFIG)因其优异的性能成为主流机型，但其定子直接连接电网的特性，使得短路故障时转子侧会感应出大电流，这对变流器构成严峻挑战。Matlab/Simulink作为电力系统仿真利器，能完整模拟从风速变化、机械传动到电气控制的整个动态过程，是研究LVRT能力的理想平台。

关键提示：低压穿越能力测试已成为各国电网准入的强制要求，我国《风电场接入电力系统技术规定》明确规定了LVRT曲线标准，电压跌至0.2pu时需保持并网至少625ms。

2. 双馈风机建模核心要素解析

2.1 机械与电气系统耦合建模

双馈风机的完整模型需要整合以下子系统：

• 气动模型：描述风能到机械能的转换
• 传动链模型：包括低速轴、齿轮箱和高速轴
• 发电机模型：双馈感应电机的详细电磁模型
• 变流器模型：转子侧和网侧变流器及其控制

在Simulink中，推荐使用Simscape Electrical库中的异步电机模块，其参数设置界面可直接输入电机铭牌数据，自动计算等效电路参数。对于1.5MW典型双馈风机，关键参数示例如下：

2.2 离散化建模的实践要点

在将连续模型离散化时需特别注意：

1. 采样时间选择：一般取控制系统开关频率的1/10~1/20，对于2kHz的PWM，采样时间建议50-100μs
2. 数值积分方法：对于电力电子系统推荐使用Tustin(双线性变换)方法，其保持数值稳定性更好
3. 代数环处理：在Simulink中启用"Minimize algebraic loop"选项，或手动插入Unit Delay模块

离散化后的电流控制器可表示为：

matlab复制% 离散PI控制器实现示例
function [output] = discretePI(error, Kp, Ki, Ts, persistent_vars)
    persistent integral;
    if isempty(integral)
        integral = 0;
    end
    
    integral = integral + error * Ts;
    output = Kp * error + Ki * integral;
    
    % 抗积分饱和处理
    if abs(integral) > max_integral
        integral = sign(integral) * max_integral;
    end
end

3. 控制策略深度剖析

3.1 MPPT控制的工程实现细节

最大功率点跟踪(MPPT)通常采用最优转矩控制，其核心在于：

1. 功率系数Cp与叶尖速比λ的关系曲线建模：

matlab复制% Cp(λ,β)曲线拟合示例
beta = pitch_angle; % 桨距角(度)
lambda = R*rotor_speed/wind_speed; % 叶尖速比

Cp = 0.22*(116/lambda_i - 0.4*beta -5)*exp(-12.5/lambda_i);
lambda_i = 1/(1/(lambda+0.08*beta) - 0.035/(beta^3+1));

2. 最优转矩计算：

matlab复制% 根据风速估计最优转速
opt_tsr = 8.1; % 最优叶尖速比
opt_rotor_speed = opt_tsr * wind_speed / R;

% 计算最优转矩参考值
K_opt = 0.5*rho*pi*R^5*Cp_max/opt_tsr^3;
opt_torque = K_opt * rotor_speed^2;

实际工程中需加入转速滤波和转矩变化率限制，典型值为±15%额定转矩/秒

3.2 双闭环控制的参数整定方法

转子侧变流器的电流内环带宽通常设为开关频率的1/5~1/10，电压外环带宽设为内环的1/5~1/10：

1. 电流环PI参数计算：

matlab复制L_sigma = Ls + Lr*(Lm^2)/(Lr^2); % 总漏感
R_total = Rs + Rr*(Lm^2)/(Lr^2); % 总电阻

% 按典型II型系统整定
omega_c_current = 2*pi*200; % 200Hz带宽
Kp_current = L_sigma * omega_c_current;
Ki_current = R_total * omega_c_current;

2. 转速环PI参数整定：

matlab复制H = 5; % 惯性时间常数(s)
omega_c_speed = omega_c_current / 5; % 40Hz带宽

Kp_speed = 2 * H * omega_c_speed;
Ki_speed = 0.5 * omega_c_speed^2 * H;

4. 低压穿越实现方案

4.1 撬棒(Crowbar)保护电路设计

当检测到直流母线电压超过1.2倍额定值时，触发撬棒电路：

1. 撬棒电阻选择原则：

matlab复制% 根据最大允许转子电流计算
I_rotor_max = 2; % pu值
R_crowbar = (0.9*V_rotor_max)/I_rotor_max;

% 典型值参考
switch power_rating
    case 1.5MW
        R_crowbar = 0.2; % Ω
    case 3MW 
        R_crowbar = 0.15; % Ω
end

2. 触发逻辑实现：

matlab复制function crowbar_trigger = check_crowbar(v_dc, i_rotor)
    persistent timer;
    if isempty(timer)
        timer = 0;
    end
    
    if v_dc > 1.2*Vdc_nominal || max(abs(i_rotor)) > 1.8*I_rated
        timer = timer + Ts;
    else
        timer = max(0, timer - Ts/2);
    end
    
    crowbar_trigger = (timer > 0.002); % 持续2ms后触发
end

4.2 动态无功支撑策略

根据电网规范要求，在电压跌落期间需提供无功电流支撑：

matlab复制function i_q_ref = LVRT_reactive_control(v_grid, t_fault)
    % 根据电压跌落深度计算无功电流参考值
    v_pu = v_grid / V_nominal;
    
    if v_pu < 0.2
        i_q_ref = 1.0; % pu
    elseif v_pu < 0.5
        i_q_ref = 0.5 + 2.5*(0.5 - v_pu);
    else
        i_q_ref = 0;
    end
    
    % 加入动态响应
    tau = 0.02; % 20ms时间常数
    i_q_ref = i_q_ref * (1 - exp(-t_fault/tau));
end

5. 短路故障仿真案例分析

5.1 三相短路故障设置

在Simulink中配置三相短路故障模块时需注意：

1. 故障电阻设置：

• 金属性短路：0.001Ω
• 经阻抗短路：根据实际需求设置(典型0.1-1Ω)

2. 故障时序配置示例：

matlab复制set_param('DFIG_model/Fault', 'FaultA', 'on');
set_param('DFIG_model/Fault', 'FaultB', 'on'); 
set_param('DFIG_model/Fault', 'FaultC', 'on');
set_param('DFIG_model/Fault', 'TransitionTimes', '[1.0]');
set_param('DFIG_model/Fault', 'FaultOnDuration', '0.1');

5.2 典型故障波形分析

正常工况与故障工况关键参数对比：

故障恢复过程中的控制要点：

1. 撬棒退出时序：检测到直流电压<1.05pu且持续10ms
2. 功率恢复速率：限制在10%额定功率/秒
3. 无功优先原则：先恢复无功输出，再恢复有功

6. 工程实践中的经验总结

1. 参数敏感性分析：

• 转子电阻误差>10%会导致电流环性能明显下降
• 电感参数误差>15%影响系统稳定性
• 建议每季度进行参数辨识校验

2. 常见问题排查指南：

3. 模型验证建议：

• 先用稳态工况验证功率曲线
• 再进行小扰动测试(5%阶跃)
• 最后进行大扰动测试(短路故障)
• 实测数据与仿真误差应<8%