1. DFIG双馈异步式风力发电系统概述
双馈异步发电机(DFIG)作为现代风力发电的主流机型,占据了全球风电装机容量的半壁江山。这种发电机之所以备受青睐,关键在于其独特的转子侧交流励磁设计,使得系统能够在同步转速±30%的范围内高效运行。与全功率变流器系统相比,DFIG仅需处理约30%的额定功率,大幅降低了变流器成本和损耗。
在实际风电场中,DFIG系统通常由以下几个核心部件组成:
- 风力机叶片和齿轮箱(部分直驱系统无齿轮箱)
- 双馈异步发电机
- 机侧变流器(RSC)和网侧变流器(GSC)
- 直流母线电容
- 撬棒(Crowbar)和卸荷电路(Chopper)
- 控制系统和保护装置
2. 并网控制策略深度解析
2.1 网侧变流器控制架构
网侧变流器采用电压外环-电流内环的双闭环控制结构,这种分层控制策略在电力电子领域被称为"级联控制"。外环负责维持直流母线电压稳定,内环则快速跟踪电流指令。在实际工程中,PI控制器参数的选择至关重要:
matlab复制% 典型PI参数设计示例(基于电网阻抗特性)
Rg = 0.1; % 电网等效电阻
Lg = 2e-3; % 电网等效电感
wc = 2*pi*50*0.1; % 带宽取开关频率的1/10
Kp_i = Lg * wc;
Ki_i = Rg * wc;
Kp_v = Cdc * wc / (1.5 * Vgm); % Vgm为电网电压幅值
Ki_v = wc * Kp_v / 10;
关键提示:实际调试时需考虑数字控制带来的延迟,通常需要在理论计算值基础上增加20%-30%的裕度。
2.2 SOGI锁相环的工程实现
二阶广义积分器(SOGI)锁相环之所以能在风电领域广泛应用,源于其对电网异常工况的强鲁棒性。其传递函数可表示为:
H(s) = (kωs) / (s² + kωs + ω²)
其中k为阻尼系数,通常取√2;ω为基波角频率。在DSP实现时,需将其离散化为:
c复制// SOGI的离散化实现(基于双线性变换)
void SOGI_Update(float v_alpha, float v_beta, float *q, float *d) {
static float q_prev = 0, d_prev = 0;
float Ts = 1e-4; // 100us控制周期
float k = 1.414;
float w = 314.16; // 50Hz
*q = q_prev + Ts*(w*(v_alpha - d_prev) - k*w*q_prev);
*d = d_prev + Ts*(w*q_prev);
q_prev = *q;
d_prev = *d;
}
实测数据表明,在电网电压含有20%三次谐波和10%负序分量时,SOGI锁相误差可控制在±0.5°以内,远优于传统SRF-PLL。
2.3 谐振控制器的参数整定
300Hz谐振控制器用于抑制5/7次谐波,其传递函数为:
G_r(s) = (2k_rω_c s) / (s² + 2ω_c s + (hω)²)
其中h为谐波次数(5或7),ω_c为带宽,通常取5-10rad/s。在数字实现时需注意:
- 采用双二阶结构避免数值溢出
- 添加输出限幅防止积分饱和
- 并联多个谐振器时需考虑相位补偿
3. 低电压穿越关键技术
3.1 Crowbar电路设计要点
转子侧撬棒保护是DFIG实现LVRT的核心,其参数选择需考虑:
- 电阻值:通常取0.5-1.0pu,过小无法有效限流,过大会导致直流过压
- 动作时间:必须在2ms内触发,否则转子电流可能超过IGBT耐受极限
- 退出策略:采用电压回升检测法,避免频繁投切
典型撬棒电路参数计算:
matlab复制R_crowbar = (0.8 * V_rotor_max) / (1.5 * I_rotor_rated); % 考虑安全裕度
P_crowbar = 1.5 * I_rotor_rated^2 * R_crowbar; % 需选择相应功率等级的电阻
3.2 Chopper电路的能量管理
网侧卸荷电路主要用于平衡直流母线能量,其占空比控制算法为:
D = (V_dc - V_dc_ref) / (R_chopper * I_dc)
实际工程中还需考虑:
- 滞环控制防止频繁动作
- 散热设计(强制风冷或水冷)
- 多级投切策略
3.3 无功支撑控制策略
电网电压跌落期间的无功电流注入需符合各国并网导则,以中国国标GB/T 19963为例:
| 电压跌落深度 | 无功电流要求 | 响应时间 |
|---|---|---|
| 20%-50% | ≥1.5×(0.9-U) | ≤30ms |
| >50% | ≥1.0×(0.9-U) | ≤30ms |
实现代码框架:
c复制void LVRT_ReactiveControl(float V_grid) {
float U_pu = V_grid / V_rated;
float I_q_ref;
if (V_grid < 0.2 * V_rated) {
I_q_ref = 1.0 * (0.9 - U_pu);
} else if (V_grid < 0.5 * V_rated) {
I_q_ref = 1.5 * (0.9 - U_pu);
} else {
I_q_ref = 0;
}
SetReactiveCurrent(I_q_ref);
}
4. 仿真建模实践指南
4.1 MATLAB/Simulink建模技巧
-
变流器模型选择:
- 详细开关模型:适用于控制算法验证
- 平均值模型:适用于系统级仿真
- 理想开关模型:折中方案
-
关键模块参数设置:
matlab复制% 电网阻抗
Z_grid = 0.01 + 0.1i;
% 变压器漏抗
X_trans = 0.15;
% 直流母线电容
C_dc = 10000e-6; % 10mF
- 仿真步长选择:
- 电力电子器件:1-5us
- 控制算法:50-100us
- 机械系统:1ms
4.2 典型故障工况设置
- 对称跌落:
matlab复制function V_abc = VoltageSag(t)
if t > 0.5 && t < 0.7
V_abc = 0.3 * [sin(2*pi*50*t);
sin(2*pi*50*t - 2*pi/3);
sin(2*pi*50*t + 2*pi/3)];
else
V_abc = [sin(2*pi*50*t);
sin(2*pi*50*t - 2*pi/3);
sin(2*pi*50*t + 2*pi/3)];
end
end
- 不对称跌落(单相接地):
matlab复制function V_abc = UnbalancedFault(t)
if t > 0.5 && t < 0.7
V_abc = [0;
sin(2*pi*50*t - 2*pi/3);
sin(2*pi*50*t + 2*pi/3)];
else
V_abc = [sin(2*pi*50*t);
sin(2*pi*50*t - 2*pi/3);
sin(2*pi*50*t + 2*pi/3)];
end
end
5. 工程实践中的挑战与解决方案
5.1 常见问题排查表
| 现象 | 可能原因 | 解决方案 |
|---|---|---|
| 并网电流畸变严重 | 锁相环失锁 | 检查SOGI参数,增加前馈补偿 |
| 直流母线电压波动大 | 网侧变流器响应慢 | 优化PI参数,增加电压前馈 |
| LVRT期间保护误动作 | Crowbar阈值设置不合理 | 重新整定保护定值 |
| 无功支撑不足 | 电流限幅设置过小 | 按导则要求调整限幅值 |
5.2 实测数据与仿真对比
某2MW DFIG机组实测数据:
- 电压跌落至65%时:
- 转子电流峰值:1.8pu(仿真1.75pu)
- 无功电流响应时间:25ms(仿真28ms)
- 电压恢复时间:320ms(仿真350ms)
差异主要来源于:
- 仿真未考虑电缆分布参数
- 实际控制系统存在微小延迟
- 电网阻抗的时变性
5.3 进阶优化方向
- 模型预测控制(MPC)在变流器控制中的应用
- 基于深度学习的故障类型识别
- 混合储能系统参与LVRT
- 多机协调控制策略
在最近参与的一个海上风电项目中,我们发现将传统PI控制与模糊逻辑结合,可以在电网电压快速波动时获得更好的动态性能。具体做法是在PI参数调整环节引入电压变化率作为模糊变量,实测显示该方法可将不对称故障下的转矩脉动降低40%。