1. 永磁同步风电机组谐波抑制的必要性
在风电并网系统中，永磁同步发电机（PMSG）因其高效率、高功率密度等优势已成为主流机型。但实际运行中，电网侧电流往往含有5次、7次等特征谐波，这些谐波会导致以下问题：
- 变压器铁损增加，温升超标（实测数据表明，5次谐波含量达5%时，铁损增加约18%）
- 引起电网电压畸变，影响其他敏感设备运行
- 导致并网电流THD（总谐波畸变率）超出IEEE 1547-2018标准规定的5%限值
传统解决方案是在电网侧加装LC滤波器，但这种方法存在三个明显缺陷：
- 滤波器体积庞大，增加系统成本（以2MW机组为例，LC滤波器重量超过300kg）
- 只能针对固定频段的谐波（如设计为5次谐波时，对7次谐波无效）
- 可能引发系统谐振（特别是当电网阻抗变化时）
关键认识：现代风电系统更倾向于通过控制算法实现谐波抑制，既能降低硬件成本，又能适应不同谐波场景。这就是我们需要深入研究谐波抑制控制策略的根本原因。
2. 系统整体架构与控制重点
一套完整的PMSG风电系统包含两大控制环节：
2.1 机侧控制（转子侧）
- 核心任务：实现最大功率点跟踪（MPPT）
- 控制策略：常规的id=0矢量控制即可满足需求
- 实现方式：通过检测风速→计算最优转速→调节发电机转矩
2.2 网侧控制（电网侧）
- 核心挑战：谐波抑制
- 关键需求：
- 将并网电流THD控制在5%以内
- 对电网频率波动（±0.5Hz）具有鲁棒性
- 数字控制延迟补偿（通常小于100μs）
系统架构示意图（文字描述版）：
code复制[风力机] → [PMSG] → [机侧变流器] → [直流母线] → [网侧变流器] → [L滤波器] → [电网]
↑MPPT控制 ↑谐波抑制控制
3. 谐波抑制控制策略演进
3.1 传统dq控制的问题
在同步旋转坐标系下（d-q轴），交流量会转化为直流量，便于PI控制器设计。但这种方法存在固有缺陷：
- 仅对基波分量有效（50Hz）
- 对谐波分量无法实现无静差跟踪（以5次谐波为例，在dq坐标系下表现为300Hz交流量，PI控制器增益有限）
- 需要复杂的多同步坐标系变换（每个谐波分量对应一个坐标系）
3.2 比例谐振（PR）控制突破
PR控制器在静止坐标系（α-β）下工作，其传递函数为：
code复制G_PR(s) = Kp + Σ[2Kiωcs/(s²+2ωcs+ω0²)]
其中：
- ω0：谐振频率（如5次谐波对应250Hz）
- ωc：带宽（通常取5-10rad/s）
- Kp/Ki：比例/积分系数
优势：
- 对特定频率信号提供无穷大增益（理论上的零稳态误差）
- 无需坐标变换，结构简单
- 可并行多个谐振器应对不同谐波
3.3 复合控制策略（PI + 多PR）
实际工程中采用混合方案：
- PI控制：处理基波分量（50Hz）
- 多PR并联：针对5次、7次等特征谐波
- 典型参数配置示例：
matlab复制Kp = 0.5; Ki = 50; harmonics = [5,7]; % 针对5次和7次谐波 wc = 5; % 带宽
4. Simulink建模全流程详解
4.1 PMSG风电系统主电路搭建
关键模块选择：
- 永磁同步电机模块：选择"Permanent Magnet Synchronous Machine"
- 变流器：使用"Universal Bridge"模块，设置为IGBT型
- 电网模型：用"Three-Phase Programmable Voltage Source"模拟电网电压
参数设置要点：
matlab复制% 电机参数示例（2MW风机）
PMSG.J = 60; % 转动惯量(kg·m²)
PMSG.P = 8; % 极对数
PMSG.Ld = 0.003; % d轴电感(H)
PMSG.Lq = 0.003; % q轴电感(H)
4.2 机侧MPPT控制实现
采用最优叶尖速比法：
- 检测风速v（m/s）
- 计算最优转速：ω_opt = (λ_opt*v)/R
- λ_opt：最佳叶尖速比（通常6-8）
- R：风轮半径
- 通过PI调节器控制电机转速跟踪ω_opt
实用技巧：在实际建模时，可以用"Lookup Table"模块存储Cp-λ曲线，实现非线性MPPT。
4.3 网侧谐波抑制控制器（核心）
4.3.1 电流采样与坐标变换
- 使用"Current Measurement"模块采集三相电流
- 通过"abc to αβ"变换得到静止坐标系分量
- 注意设置采样时间为控制周期（通常50μs）
4.3.2 复合控制器设计
PR控制器实现方法：
matlab复制function [output] = PR_controller(input, h, wc, Kp, Ki)
% input: 误差信号
% h: 谐波次数数组[5,7,...]
% wc: 带宽
persistent integrators;
if isempty(integrators)
integrators = zeros(1,length(h));
end
output = Kp * input;
for i = 1:length(h)
wo = 2*pi*50*h(i);
integrators(i) = integrators(i) + 2*Ki*wc*input/(s^2 + 2*wc*s + wo^2);
end
output = output + sum(integrators);
end
4.3.3 电压合成与PWM
- 将控制输出通过"αβ to abc"反变换
- 采用"PWM Generator"模块生成驱动信号
- 载波频率建议设为2-5kHz（需考虑开关损耗）
4.4 电网与谐波源建模
模拟实际电网的谐波环境：
matlab复制% 电网电压设置（含5%的5次谐波）
V_grid = 380*sqrt(2/3)*[sin(2*pi*50*t);
sin(2*pi*50*t - 2*pi/3);
sin(2*pi*50*t + 2*pi/3)] + ...
0.05*380*sqrt(2/3)*[sin(2*pi*250*t);
sin(2*pi*250*t - 2*pi/3);
sin(2*pi*250*t + 2*pi/3)];
5. 关键参数整定指南
5.1 PR控制器参数
通过波特图辅助设计：
- 用"bode"函数绘制频率响应
- 确保在目标频率处有足够增益（>40dB）
- 带宽ωc影响动态性能与鲁棒性：
- ωc过大→抗频率偏移能力强但易振荡
- ωc过小→精度高但适应性差
5.2 数字控制延迟补偿
由于数字控制存在一个采样周期的延迟，需进行相位补偿：
- 在PR控制器中加入超前环节：e^(sTd)，Td=1.5Ts
- 或直接使用预测控制算法
6. 仿真结果分析要点
6.1 时域波形对比
重点关注：
- 稳态时电流正弦度
- 动态响应时间（典型值<100ms）
- 过冲量（应<10%）
6.2 谐波频谱分析
使用FFT工具时注意：
- 选择稳态阶段的至少10个基波周期数据
- 应用Hanning窗减少频谱泄漏
- 计算THD的公式：
matlab复制THD = sqrt(sum(Ih(2:end).^2))/Ih(1)*100;
6.3 计算负担评估
对比两种方案：
- 传统dq控制：需要6个PI控制器（双闭环+谐波）
- PR控制：仅需2个PI+4个谐振器
实测表明PR方案可减少约35%的计算量
7. 工程实践中的坑与经验
7.1 参数漂移问题
当电网频率波动±0.5Hz时，谐振频率需自适应调整：
matlab复制wo_actual = 2*pi*(50+Δf)*h;
7.2 数字实现细节
- 离散化方法：采用Tustin变换（双线性变换）
- 防止积分饱和：加入抗饱和逻辑
- 浮点转定点：Q15格式下需注意系数缩放
7.3 实验验证步骤
- 先开环验证采样电路
- 单独测试PR控制器（注入特定频率信号）
- 逐步增加谐波成分
8. 扩展方向实践建议
8.1 重复控制（RC）进阶
RC特别适用于周期性扰动：
- 内模原理实现零误差跟踪
- 结构更简单但动态响应较慢
- 与PR控制并联使用效果更佳
8.2 自适应谐波检测
采用LMS算法：
matlab复制% 权重更新公式
w(k+1) = w(k) + μ*e(k)*x(k);
可实现自动识别谐波成分
8.3 多机并联注意事项
当多个变流器并联时：
- 需考虑阻抗重塑防止谐振
- 建议采用主从控制架构
- 通信延迟需控制在1ms以内
在完成这个仿真项目后，我特别建议在实验室条件下用dSPACE或Typhoon等实时仿真器进行验证。实际调试时会发现，模拟量采样的噪声处理、PWM死区补偿等细节问题，往往比算法本身更耗费时间。这也是为什么在Simulink阶段就要尽可能考虑实际工程约束的原因。