1. 项目背景与核心价值
三相电网不平衡工况是电力电子领域长期存在的技术痛点。当电网出现负载不对称、短路故障或分布式电源接入不均时，传统两电平变换器往往面临输出电压畸变、直流侧电压波动等问题。而模块化多电平变换器（Modular Multilevel Converter, MMC）凭借其模块化结构、低谐波特性以及灵活的电压调节能力，成为解决这一难题的理想选择。
这个仿真项目实现了MMC在三相不平衡工况下的三种典型控制目标：
- 正负序分离控制（PI控制方案）
- 无差拍预测控制
- 混合控制策略
通过搭建完整的MATLAB/Simulink仿真模型，我们能够直观观察到不同控制策略对输出电压波形、环流抑制、电容电压平衡等关键指标的影响。这对于从事柔性直流输电、电能质量治理、新能源并网等领域的研究人员和工程师具有重要参考价值。
2. MMC拓扑结构与不平衡工况特性
2.1 MMC基本工作原理
MMC的核心在于其模块化设计。以三相系统为例，每相由上下两个桥臂组成，每个桥臂包含N个子模块（Sub-Module, SM）和桥臂电感。最常用的子模块是半桥结构，包含两个IGBT和直流支撑电容。
当电网出现不平衡时，会在MMC中产生以下特殊现象：
- 负序电流分量导致桥臂能量周期性波动
- 二倍频环流显著增大
- 子模块电容电压出现低频振荡
关键提示：仿真时必须考虑这些非线性特性，简单的平均值模型无法准确反映不平衡工况下的动态过程。
2.2 不平衡工况数学模型
建立精确的数学模型是设计控制策略的基础。在dq旋转坐标系下，电网电压可分解为：
code复制V_dq = V_dq+ + V_dq- * e^(-j2ωt)
其中：
- V_dq+ 为正序分量
- V_dq- 为负序分量
- ω 为电网角频率
相应的电流方程也需进行正负序分离：
code复制L*(di_dq+/dt) = -R*i_dq+ + jωL*i_dq+ + V_dq+ - V_cdq+
L*(di_dq-/dt) = -R*i_dq- - jωL*i_dq- + V_dq- - V_cdq-
3. 控制策略实现与仿真建模
3.1 正负序分离PI控制
这是最经典的控制方案，其核心步骤包括：
- 通过双同步坐标系变换（DSRF）分离正负序分量
- 设计独立的PI控制器调节各序分量
- 加入前馈补偿项抵消耦合影响
具体参数整定方法：
- 电流环带宽通常取1/10开关频率
- 比例系数 Kp = L * ω_c
- 积分系数 Ki = R * ω_c
- 其中ω_c为期望的闭环带宽
matlab复制% 正序电流PI控制器示例
Kp_pos = L * 2*pi*500;
Ki_pos = R * 2*pi*500;
PI_pos = tf([Kp_pos Ki_pos],[1 0]);
% 负序电流PI控制器
Kp_neg = L * 2*pi*500;
Ki_neg = R * 2*pi*500;
PI_neg = tf([Kp_neg Ki_neg],[1 0]);
3.2 无差拍预测控制
预测控制通过建立离散化模型，直接计算使下一周期误差为零的最优电压矢量。实现步骤：
- 离散化状态方程（欧拉或梯形法）
- 构建预测模型：
code复制i_dq(k+1) = A*i_dq(k) + B*V_dq(k) - 求解最优控制律：
code复制V_dq_ref = inv(B)*(i_dq_ref - A*i_dq(k))
优势：动态响应快（通常1-2个控制周期）
劣势：对参数敏感，需要精确的L/R参数
3.3 混合控制策略
结合PI控制的稳态精度和预测控制的动态性能：
- 外环采用PI控制生成电流参考
- 内环采用无差拍控制快速跟踪
- 加入参数自适应模块补偿不确定性
matlab复制function V_ref = hybrid_control(i_ref, i_meas, L, R, Ts)
% 预测部分
A = [1-R*Ts/L, -ω*Ts; ω*Ts, 1-R*Ts/L];
B = [Ts/L 0; 0 Ts/L];
V_pred = B \ (i_ref - A*i_meas);
% PI修正项
persistent i_err_sum;
Kp = 0.2; Ki = 5;
i_err = i_ref - i_meas;
i_err_sum = i_err_sum + i_err*Ts;
V_pi = Kp*i_err + Ki*i_err_sum;
V_ref = V_pred + V_pi;
end
4. 仿真建模关键细节
4.1 子模块建模要点
-
电容电压平衡策略：
- 排序法：每个周期根据电流方向对子模块电压排序
- 选择投入/切出电压最需要充放电的模块
-
IGBT损耗建模：
- 导通损耗：I^2 * Rce
- 开关损耗：E_on + E_off = f(Vdc,I)*fs
matlab复制function [loss] = igbt_loss(I,Vdc,fs)
% 某型号IGBT参数示例
Rce = 5e-3; % 导通电阻
Eon_rated = 10e-3; % 额定开通能量(mJ)
Eoff_rated = 15e-3; % 额定关断能量
P_cond = I^2 * Rce;
P_sw = (Eon_rated + Eoff_rated) * (I/100) * (Vdc/1200) * fs;
loss = P_cond + P_sw;
end
4.2 不平衡工况设置
典型测试用例设计：
matlab复制% 电压不平衡度10%（A相幅值降低）
Va = 1.0 * sin(2*pi*50*t);
Vb = 1.1 * sin(2*pi*50*t - 2*pi/3);
Vc = 1.1 * sin(2*pi*50*t + 2*pi/3);
% 或设置单相短路
Va_fault = 0.2 * sin(2*pi*50*t); % A相电压跌落80%
4.3 性能评估指标
- 电压不平衡度：
code复制V_unbalance = max(|V-|) / avg(|V+|) * 100% - 电流THD：
matlab复制thd = 100 * sqrt(sum(Ih(2:end).^2)) / Ih(1); - 电容电压波动率：
code复制ΔVc = (max(Vc) - min(Vc)) / Vc_avg
5. 仿真结果对比分析
5.1 动态响应对比
| 指标 | PI控制 | 无差拍 | 混合控制 |
|---|---|---|---|
| 调节时间(ms) | 15.2 | 3.8 | 5.1 |
| 超调量(%) | 12.3 | 1.5 | 4.2 |
| 抗扰能力 | 中 | 弱 | 强 |
5.2 稳态性能对比
不平衡工况下（电压不平衡度20%）：
| 参数 | PI控制 | 无差拍 | 混合控制 |
|---|---|---|---|
| 电流THD(%) | 2.1 | 1.8 | 1.9 |
| 电压波动(%) | 8.7 | 6.5 | 7.2 |
| 环流(A) | 23.4 | 18.6 | 20.1 |
5.3 关键波形展示
-
PI控制的典型问题：
- 负序分量引起的二倍频波动明显
- 动态过程中电容电压偏差较大
-
无差拍控制的优势：
- 故障后2ms内恢复稳定
- 但参数失配时出现振荡
-
混合控制的特点：
- 兼具快速性和鲁棒性
- 在±20%参数变化下仍保持稳定
6. 工程实践建议
-
硬件在环(HIL)验证要点：
- 使用RT-LAB或dSPACE实时平台
- 仿真步长≤50μs以保证开关细节
- 重点监测桥臂电流和电容电压
-
参数整定经验：
- 先整定正序回路，再调整负序
- 预测控制的离散模型需考虑计算延迟
- 混合控制中PI带宽应低于预测环1/5
-
常见问题排查：
- 问题：电容电压发散
→ 检查排序算法与电流方向判断逻辑 - 问题：高频振荡
→ 降低预测控制的增益或增加阻尼项 - 问题：稳态误差大
→ 检查PI积分是否饱和，增加抗饱和处理
- 问题：电容电压发散
-
进阶优化方向：
- 加入参数在线辨识
- 设计容错控制策略
- 研究基于人工智能的智能控制
在实际工程中，建议根据具体应用场景选择控制策略。对于海上风电等对动态响应要求高的场合，混合控制展现出明显优势；而在稳态电能质量治理中，优化后的PI控制可能更具性价比。