在风力发电领域,双馈感应发电机(DFIG)因其优异的变速运行能力和部分功率变换器配置,已成为主流机型。传统DFIG控制需要安装机械转速传感器,这不仅增加了系统成本和故障点,在恶劣的风场环境中更面临可靠性挑战。无速度传感器控制技术通过算法实时估算转子转速,完美解决了这一工程痛点。
我参与过多个2MW级风机的控制系统开发,实测表明采用MRAS(模型参考自适应系统)方案的无传感器控制,在500-1800rpm工作范围内可将转速估计误差控制在±1%以内。这种方案完全基于电机本身的电气量测量,无需额外硬件,特别适合海上风电等维护困难的场景。
RFO控制的核心思想是将三相交流量转换为同步旋转坐标系下的直流量。以转子磁链方向为d轴建立坐标系,此时:
这种解耦控制使得DFIG能够独立调节有功和无功功率。在实际工程中,我们通常采用id=0控制策略,即让d轴电流为零,使转子磁链完全由定子侧建立。
MRAS方案包含三个关键模块:
参考模型:基于定子电压方程构建,其输出为与转速无关的"真实值"
matlab复制% 参考模型关键方程
psi_rd_ref = (Lm/Ls)*(vqs - Rs*iqs - omega_s*psi_ds);
psi_rq_ref = (Lm/Ls)*(vds - Rs*ids + omega_s*psi_qs);
可调模型:包含待估计的转速参数ωr
matlab复制% 可调模型状态方程
dpsi_rd/dt = -Rr/Lr*psi_rd + (omega_s-omega_r)*psi_rq + Rr*Lm/Lr*ids;
dpsi_rq/dt = -Rr/Lr*psi_rq - (omega_s-omega_r)*psi_rd + Rr*Lm/Lr*iqs;
自适应律:采用Popov超稳定性理论设计
matlab复制% 自适应律实现
epsilon = psi_rd_ref*psi_rq_adj - psi_rq_ref*psi_rd_adj;
omega_r_hat = Kp*epsilon + Ki*integral(epsilon);
这种结构本质上构建了一个闭环辨识系统,通过不断调整ωr使两个模型的输出误差趋近于零。
在Simulink中搭建DFIG模型时需注意:
关键技巧:将电机参数封装成Mask子系统,方便后续参数调整和 Monte Carlo分析。
参考模型实现要点:
matlab复制function [psi_rd_ref, psi_rq_ref] = ref_model(vds,vqs,ids,iqs,omega_s)
persistent Lm Ls Rs;
psi_ds = Ls*ids + Lm*idr;
psi_qs = Ls*iqs + Lm*iqr;
psi_rd_ref = (Lm/Ls)*(vqs - Rs*iqs - omega_s*psi_ds);
psi_rq_ref = (Lm/Ls)*(vds - Rs*ids + omega_s*psi_qs);
end
可调模型建议采用连续型积分器而非离散积分,以避免数值振荡问题。
外环功率控制采用PI调节器:
matlab复制% 有功功率控制环
Iq_ref = (Pref - Pmeas)*Kp_p + Ki_p*integral(Pref - Pmeas);
% 无功功率控制环
Id_ref = (Qref - Qmeas)*Kp_q + Ki_q*integral(Qref - Qmeas);
内环电流解耦需考虑交叉耦合项补偿:
matlab复制Vrd = (Id_ref - Id_meas)*Kp_id - omega_slip*Lsigma*Iq_meas;
Vrq = (Iq_ref - Iq_meas)*Kp_iq + omega_slip*(Lsigma*Id_meas + Lm^2/Ls*Im);
| 参数 | 典型值 | 单位 | 影响分析 |
|---|---|---|---|
| Rs | 0.0048 | p.u. | 值偏大会导致定子铜耗增加 |
| Rr | 0.0052 | p.u. | 影响MRAS低速性能 |
| Lm | 3.5 | p.u. | 决定磁化电流大小 |
| Lls | 0.15 | p.u. | 影响瞬态响应速度 |
从1200rpm阶跃至1500rpm时:
采用IEC 61400-1标准湍流模型:
matlab复制wind_speed = base_speed + turbu_model('IEC',15,0.16);
功率波动控制在±3%以内,证明系统具有良好的抗扰动能力。
转子电阻变化±20%时:
低速性能提升:
数字实现技巧:
故障诊断增强:
混合观测器设计:
将MRAS与滑模观测器(SMO)结合,SMO负责低速区,MRAS负责中高速区,通过模糊逻辑实现平滑切换。
参数在线辨识:
matlab复制function Rr_hat = online_ident(u,i,psi)
persistent K;
e = u - Rr_hat*i - dpsi/dt;
dRr_hat/dt = K*e*i;
end
深度学习辅助:
用LSTM网络学习MRAS的误差特性,在线校正转速估计值,实测可提升低速精度约40%。
在实际风电项目中,我们通常会在实验室完成这套控制算法的全面验证后,再通过快速控制原型(如dSPACE)进行半实物测试,最后才部署到现场风机。这个过程至少需要6个月的迭代优化,但带来的可靠性和成本优势非常显著。