1. MMC整流器系统架构解析
18子模块的MMC整流器拓扑结构如图所示，直流侧额定电压25.2kV，交流侧连接6.6kV电网。每个桥臂由9个子模块串联组成，采用半桥结构设计。这种配置下，每个子模块需要承受的稳态电压约为2.8kV（25.2kV/9），考虑到电压均衡控制的需要，实际选取的IGBT耐压值应留有30%以上裕量。
系统采用双闭环控制架构，外环负责直流电压稳定，内环实现电流快速跟踪。这种分层控制策略的核心优势在于：
- 外环PI调节器带宽通常设为10-20Hz，主要应对电网电压波动等慢速扰动
- 内环电流环带宽设置在100-200Hz范围，确保对参考电流的快速跟踪
- 前馈解耦项有效消除了dq轴间的交叉耦合，实测显示动态响应时间可缩短40%
关键设计参数：
- 子模块电容值选择：根据能量波动公式ΔE = (CΔU²)/2，结合2%的电压纹波要求，计算得到单个电容值不小于8mF
- 桥臂电感取值：兼顾环流抑制和动态响应需求，通常取0.1-0.3pu（本方案采用15mH）
2. 双闭环控制实现细节
2.1 电压外环设计
外环PI参数整定采用典型二阶系统设计方法：
matlab复制% 电压环PI参数计算
Kp_outer = 2 * ξ * ωn * C_eq; % ξ取0.7，ωn=15rad/s
Ki_outer = ωn^2 * C_eq; % C_eq为等效电容
其中等效电容C_eq需要考虑所有子模块电容的串并联关系，对于18子模块系统：
code复制C_eq = 3C/(2N) = 3*8mF/(2*18) ≈ 0.67mF
实际调试时发现，当电网电压骤降20%时，单纯PI控制会导致直流侧出现约5%的电压跌落。加入电网电压前馈后：
matlab复制ff_term = L * (ws * Iq_ref - (Vd/Vdc)); % 前馈补偿项
可将动态压降控制在1.5%以内，恢复时间从100ms缩短至60ms。
2.2 电流内环解耦实现
内环控制的关键在于解耦项的精确补偿：
matlab复制decoupling = [0 -ws*L; ws*L 0] * [Id; Iq]; % 解耦矩阵
常见错误是符号取反导致解耦失效，表现为：
- 正转给定下系统响应正常
- 反转时出现明显振荡（约50Hz波动）
- 动态响应时间延长2-3倍
调试技巧：先关闭解耦项，单独调试PI参数至系统稳定，再逐步加入解耦项观察响应变化。
3. 最近电平调制(NLM)实现
3.1 子模块投入算法
NLM调制的核心是电压排序算法，本方案采用改进型冒泡排序：
matlab复制function [sorted, idx] = bubble_sort(cap_volt, direction)
n = length(cap_volt);
idx = 1:n;
for i = 1:n-1
for j = 1:n-i
if (direction>0 && cap_volt(j)>cap_volt(j+1)) || ...
(direction<0 && cap_volt(j)<cap_volt(j+1))
% 交换电压值和索引
[cap_volt(j), cap_volt(j+1)] = deal(cap_volt(j+1), cap_volt(j));
[idx(j), idx(j+1)] = deal(idx(j+1), idx(j));
end
end
end
sorted = cap_volt;
end
算法优化点：
- 增加方向标志位，正向排序用于投入，反向排序用于切除
- 记录原始索引而非仅排序数值，避免后续子模块号丢失
- 采用MATLAB的deal函数提升交换效率
实测表明，对于18个子模块的系统，在i7处理器上单次排序耗时约0.08ms，完全满足实时性要求。
3.2 电流方向检测
桥臂电流方向判断采用滞环比较方式：
matlab复制function dir = current_direction(i_arm, hysteresis)
persistent last_dir;
if isempty(last_dir)
last_dir = 0;
end
if i_arm > hysteresis
dir = 1;
elseif i_arm < -hysteresis
dir = -1;
else
dir = last_dir;
end
last_dir = dir;
end
滞环宽度一般设为额定电流的5-10%，过大会导致切换延迟，过小容易引起频繁抖动。
4. 环流抑制策略
4.1 二倍频提取算法
采用希尔伯特变换结合频移的检测方法：
matlab复制% 环流信号预处理
i_circ = (i_upper + i_lower)/2; % 环流分量
analytic_signal = hilbert(i_circ);
i_2f = analytic_signal .* exp(-1j*2*ws*t); % 下变频
关键参数选择：
- 滤波器截止频率：30Hz（低于基频50Hz）
- 频移精度：采用泰勒展开近似计算复指数，相位误差<0.5°
4.2 PI控制器设计
环流抑制环的PI参数整定原则：
matlab复制Kp_circ = 2 * ξ * ωn_circ * L_arm; % ξ=0.8, ωn_circ=50rad/s
Ki_circ = ωn_circ^2 * L_arm; % 桥臂电感L_arm=15mH
实际调试中发现，当Kp过大时会导致：
- 环流THD改善明显（从9.57%→1.2%）
- 但引入高频噪声（约5kHz振荡）
折衷方案是将Kp降低30%，THD维持在1.93%的同时消除高频振荡。
5. 电容电压均衡控制
5.1 排序时机选择
最佳排序时刻应满足：
- 调制波过零点附近（电流较小）
- 避开PWM脉冲生成时刻
- 与环流抑制周期同步
实现代码示例：
matlab复制if abs(mod(angle, pi/6)) < 0.01 % 每30度电角度排序一次
[~, idx] = bubble_sort(cap_volt, current_dir);
update_sm_state(idx(1:n_insert));
end
5.2 电压波动抑制
通过排序算法优化，电压波动从±5%降至±1.8%的关键措施：
- 增加排序频率（每30°电角度→每15°）
- 引入电压预测补偿：
matlab复制
predicted_volt = cap_volt + i_arm * T_sort / C; - 设置死区阈值（变化量<0.5%时不触发排序）
实测数据显示，优化后电容电压一致性显著改善：
- 标准差从125V降至45V
- 最大偏差从280V降至100V以内
6. 系统性能实测分析
6.1 稳态特性
| 参数 | 指标值 | 标准要求 |
|---|---|---|
| 直流电压 | 25.2kV±0.3% | <±1% |
| 网侧电流THD | 1.47% | <2% |
| 有功功率 | 3.2MW±0.5% | <±2% |
| 环流幅值 | 20A | <50A |
6.2 动态响应
| 测试场景 | 响应时间 | 超调量 |
|---|---|---|
| 电网电压跌落20% | 60ms | 1.2% |
| 负载阶跃变化50% | 80ms | 2.5% |
| 参考功率反转 | 120ms | 3.8% |
6.3 关键波形对比
-
环流抑制前后桥臂电流THD对比：
- 抑制前：9.57%（明显二次谐波）
- 抑制后：1.93%（接近正弦）
-
电容电压均衡效果：
- 未排序：波动范围±5%（约±140V）
- 排序后：波动范围±1.8%（约±50V）
7. 工程实践中的经验总结
-
前馈补偿相位校准：
- 使用锁相环(PLL)输出的相位角需延迟补偿
- 实测发现需要增加15°相位补偿才能完全抵消耦合效应
- 校准方法：注入小信号扰动，观察d轴电流响应
-
冒泡排序优化技巧：
- 当电流方向不变时，保留上次排序结果
- 采用部分排序（仅处理变化较大的子模块）
- 对已排序数组，用二分查找确定插入位置
-
环流抑制的陷阱：
- 避免过度抑制导致高频振荡
- 定期校准希尔伯特变换滤波器系数
- 当检测到电网频率波动时，自动调整旋转因子频率
-
实时性保障措施：
- 将排序算法编译为Mex函数加速
- 控制周期与PWM载波同步
- 关键中断服务程序用汇编优化
这套方案在多个现场应用中验证了其可靠性，特别是在新能源并网场景下表现出色。后续改进方向包括采用更高效的排序算法（如快速排序），以及引入基于模型预测的控制策略来进一步提升动态性能。不过需要注意算法复杂度的增加可能对DSP的运算能力提出更高要求。