1. 双馈风力发电机（DFIG）系统概述
双馈感应发电机（Doubly-Fed Induction Generator, DFIG）作为现代风力发电系统的核心部件，因其独特的转子侧变频器结构和优异的变速恒频运行特性，已成为2-5MW功率范围风电机组的主流技术方案。与传统同步发电机相比，DFIG系统通过仅处理约30%额定功率的转子侧变流器即可实现宽转速范围（±30%同步转速）内的最大功率点跟踪（MPPT），大幅降低了变流器成本和系统损耗。
在实际风电场运行中，DFIG需要应对三大核心挑战：一是风速随机性导致的机械转矩波动，二是电网电压跌落时的暂态过电流保护，三是满足各国电网规范要求的低电压穿越（LVRT）能力。这些需求直接推动了矢量控制、直接功率控制等现代控制策略在DFIG系统中的应用演进。
2. DFIG数学模型构建原理
2.1 三相静止坐标系下的基本方程
DFIG的完整电磁模型包含6组电压方程和6组磁链方程。以定子侧为例，其电压方程可表示为：
code复制u_{abc,s} = R_s i_{abc,s} + dψ_{abc,s}/dt
其中下标s表示定子侧，R_s为定子电阻，i和ψ分别代表电流和磁链。类似的方程也存在于转子侧。这种三相模型虽然物理意义明确，但存在时变耦合项，直接用于控制器设计极为困难。
2.2 dq同步旋转坐标系变换
通过Park变换将三相量投影到与同步转速ω_s旋转的dq坐标系，得到解耦的电压方程：
code复制u_{ds} = R_s i_{ds} - ω_s ψ_qs + dψ_ds/dt
u_{qs} = R_s i_{qs} + ω_s ψ_ds + dψ_qs/dt
此时磁链方程简化为：
code复制ψ_ds = L_s i_ds + L_m i_dr
ψ_qs = L_s i_qs + L_m i_qr
其中L_s为定子全电感，L_m为互感。这种形式揭示了DFIG的电磁本质——通过转子电流i_dr、i_qr可独立控制定子侧的磁链和转矩。
2.3 功率传输方程
DFIG的有功功率P和无功功率Q在dq坐标系下表示为：
code复制P = 3/2 (u_ds i_ds + u_qs i_qs)
Q = 3/2 (u_qs i_ds - u_ds i_qs)
当定子电压定向于q轴（u_ds=0）时，控制策略可进一步简化。这种定向方式成为后续矢量控制的基础。
3. Simulink建模关键技术实现
3.1 电气子系统建模细节
在Simulink中构建DFIG本体模型时，需特别注意：
-
参数归一化处理：将实际值转换为标幺值（pu）可提高数值稳定性。例如额定功率S_base=2MW，定子电压U_base=690V，则阻抗基值Z_base=U_base²/S_base≈0.238Ω。
-
磁链计算模块：采用积分器实现ψ=∫(u-Ri)dt时，必须添加高通滤波器消除直流偏置，截止频率通常设为0.1Hz。实测表明，未处理的直流偏置会导致稳态误差超过5%。
-
坐标系变换模块：Park变换需要精确的转子位置角θ_r。实践中采用编码器信号经PLL处理后获得，角度误差应控制在±0.1°以内。
3.2 机械传动链建模
风力机与发电机之间的传动系统采用两质量块模型：
code复制J_t dω_t/dt = T_aero - K_s(θ_t-θ_r) - D_s(ω_t-ω_r)
J_g dω_r/dt = K_s(θ_t-θ_r) + D_s(ω_t-ω_r) - T_em
其中J_t和J_g分别为风轮和发电机惯量，K_s和D_s为轴系刚度和阻尼系数。在2MW机组典型参数中，K_s≈2.8×10⁶ N·m/rad，D_s≈3×10⁴ N·m·s/rad。该模型能准确反映0.1-10Hz范围内的扭振特性。
3.3 变流器与控制系统实现
转子侧变流器（RSC）和网侧变流器（GSC）均采用两电平电压源型拓扑，开关频率通常为2-4kHz。关键实现要点：
-
RSC电流环设计：
- 内环带宽取开关频率的1/5～1/10，如4kHz开关频率对应400-800Hz带宽
- PI参数计算：K_p = Lσ·ω_c，K_i = R·ω_c（Lσ为漏感，ω_c为截止角频率）
- 实测表明，d轴电流响应时间应<5ms，超调量<3%
-
GSC控制策略：
- 直流母线电压控制环带宽设为20-50Hz
- 采用电网电压定向（VOC）时，需锁相环（PLL）动态响应时间<20ms
- 典型参数下，直流电压波动应控制在±5%额定值以内
4. 先进控制策略仿真分析
4.1 传统矢量控制（VC）实现
基于定子磁链定向的VC策略框图包含：
- 外环：有功-无功功率控制器（带宽10-20Hz）
- 中环：转子电流解耦控制器（带宽100-200Hz）
- 内环：PWM调制环节
典型阶跃响应指标：
- 功率指令跟踪时间：80-120ms
- 额定工况下效率：94-96%
- LVRT期间直流过电压：<1.15pu
4.2 直接功率控制（DPC）优化
采用滞环比较器的经典DPC存在开关频率不固定的缺陷。改进方案包括：
-
空间矢量调制（SVM-DPC）：
- 开关频率固定为2kHz
- 功率纹波降低40%以上
- 需设计基于Lyapunov函数的电压矢量选择算法
-
模型预测控制（MPC-DPC）：
- 采样周期50-100μs
- 代价函数包含功率误差和开关损耗项
- 仿真显示THD可降至3%以下
4.3 低电压穿越（LVRT）策略验证
根据GB/T 19963-2021标准，在0.9pu电压跌落时：
-
硬件保护：
- 撬棒电路（Crowbar）触发阈值：1.2pu转子电流
- 动作时间：<2ms
- 能耗电阻功率：20%额定功率
-
软件策略：
- 正负序分离控制
- 无功电流优先注入（需提供0.2pu以上无功支撑）
- 仿真显示可维持0.625s不脱网
5. 仿真案例与结果分析
5.1 风速阶跃响应测试
设置风速从8m/s阶跃至12m/s：
- 叶尖速比λ从7.2调整至8.1
- 功率系数Cp保持在0.48最大值附近
- 仿真显示MPPT跟踪误差<1.5%
5.2 电网对称跌落测试
模拟0.3pu电压跌落持续0.5s：
- 转子电流峰值限制在1.8pu（未保护时可达3pu）
- 直流母线电压最高1.12pu
- 无功支撑电流达0.25pu
5.3 参数敏感性分析
关键发现：
- 转子电阻误差>10%会导致转矩控制偏差>5%
- 互感L_m误差影响磁链观测精度，需在线辨识
- 电流传感器相位偏差>1°会引发明显振荡
6. 工程实践中的关键问题
6.1 实时仿真挑战
当仿真步长>50μs时可能出现：
- 变流器开关细节丢失
- 电流环振荡（Nyquist频率限制）
- 解决方案：采用FPGA硬件在环（HIL）
6.2 模型降阶技巧
保留关键动态的简化方法：
- 忽略定子暂态（时间常数约100ms）
- 用一阶惯性环节代替PLL（带宽50Hz）
- 简化传动链为单质量块（仅研究>1Hz动态时）
6.3 代码生成优化
从Simulink到DSP的自动代码生成需注意：
- 定点数格式选择（电流环用Q12，速度环用Q20）
- 避免代数环（加入单位延迟模块）
- 关键中断服务程序（ISR）周期≤100μs