MMC变换器仿真与调制算法实现详解

1. MMC基础与仿真框架搭建

模块化多电平变换器（Modular Multilevel Converter, MMC）作为高压直流输电领域的革命性拓扑，其核心价值在于通过子模块（Sub-Module, SM）的级联组合实现高压输出，同时保持器件应力在合理范围。我在参与某±350kV柔性直流输电项目时，深刻体会到MMC仿真对于实际工程的前导验证价值。

1.1 MMC基本工作原理

MMC的典型三相结构包含六个桥臂，每个桥臂由N个子模块与一个桥臂电抗器串联构成。子模块通常采用半H桥结构，其核心功能状态如下：

• 投入状态（Inserted）：上管导通，子模块输出电压为正端电容电压
• 切除状态（Bypassed）：下管导通，子模块输出电压为零
• 闭锁状态（Blocked）：所有IGBT关断，用于故障保护

电容电压平衡是MMC稳定运行的关键。以10个子模块为例，理想情况下每个电容电压应维持在Vdc/N（假设Vdc=20kV，则每个电容2kV）。实际运行中由于器件特性差异，电压波动需控制在±5%以内。

1.2 MATLAB仿真环境配置

推荐使用MATLAB/Simulink R2021a及以上版本进行MMC仿真，关键工具箱包括：

• Simscape Power Systems（必备）
• Control System Toolbox（用于环流抑制）
• Parallel Computing Toolbox（加速仿真）

基础仿真框架搭建步骤：

1. 创建子系统封装每个子模块，包含IGBT、反并联二极管、电容等元件
2. 构建arm级模型，集成N个子模块和桥臂电感
3. 组合三个相单元，连接直流母线和中点
4. 添加测量模块（电压/电流传感器）
5. 搭建控制子系统（调制算法+闭环控制）

注意：仿真步长建议设为1e-6秒以下，使用ode23tb求解器兼顾精度与速度。首次运行时建议先采用简化模型（如N=4），验证通过后再扩展至实际规模。

2. 调制算法实现与对比

2.1 最近电平逼近调制（NLC）

NLC的本质是通过四舍五入实现电压量化，其MATLAB实现核心代码如下：

matlab复制function n = NLC(u_ref, N)
    % u_ref: 归一化参考电压 [-1,1]
    % N: 每臂子模块数
    n = round((u_ref + 1) * N / 2);  % 映射到[0,N]区间后取整
end

实测数据表明，当N=20时：

• THD：8.7%（无额外滤波）
• 开关损耗：仅为PWM方案的35%
• 实时性：执行时间<10μs（i7-11800H处理器）

但NLC存在明显局限：当N<10时，低次谐波分量急剧增加。例如N=4时，5次谐波可达15%，必须配合桥臂电抗器使用。

2.2 载波层叠PWM（CPS-PWM）

传统CPS-PWM采用垂直堆叠的三角载波，实现要点：

matlab复制carrier = zeros(N, samples);
for k = 1:N
    carrier(k,:) = (k-0.5)/N - sawtooth(2*pi*fc*t, 0.5)/N;
end

关键参数选择原则：

• 载波频率fc：通常取1-2kHz，过高会导致开关损耗剧增
• 调制比m：建议运行在0.8以下避免过调制
• 载波同步：必须确保三相载波相位差120°

2.3 载波移相PWM（PS-PWM）

PS-PWM通过相位交错提升等效开关频率，其实现方式：

matlab复制phase_shift = 2*pi/(2*N);  % 注意是2N相位间隔
for k = 1:N
    carrier(k,:) = sawtooth(2*pi*fc*t + (k-1)*phase_shift, 0.5);
end

实测对比数据（N=10, fc=1kHz）：

经验：对于实时性要求高的场合（如RT-LAB），建议采用CPS-PWM；对电能质量要求严格的场合，PS-PWM更优。

3. 环流抑制技术实现

3.1 二倍频环流机理

环流（Circulating Current）主要包含：

• 直流分量：由功率传输需求决定
• 二倍频分量：由相间电压不平衡引起
• 高频分量：开关过程导致

其中二倍频分量危害最大，会导致：

• 电容电压波动增加300%
• 桥臂电流有效值上升50%
• 器件温升提高20K以上

3.2 负序环流抑制算法

基于旋转坐标系的抑制方案实现步骤：

1. 提取相间电压差v_diff = (va + vc - 2vb)/3
2. Clark变换到αβ坐标系：

   matlab复制alpha = (2*v_diff_a - v_diff_b - v_diff_c)/3;
   beta = (v_diff_b - v_diff_c)/sqrt(3);
   

3. Park变换到负序同步旋转坐标系（-2ωt）：

   matlab复制d = alpha.*cos(-2*omega*t) + beta.*sin(-2*omega*t);
   q = -alpha.*sin(-2*omega*t) + beta.*cos(-2*omega*t);
   

4. PI调节器输出抑制信号：

   matlab复制Kp = 0.5; Ki = 100;
   icir_ref = Kp*q + Ki*integral(q);
   

参数整定技巧：

• 先设置Ki=0，增大Kp至系统开始振荡
• 取振荡时Kp值的60%作为最终值
• 逐步增加Ki，观察电容电压波动不再明显改善时停止

3.3 抑制效果评估

某工程案例实测数据：

避坑指南：避免过度抑制导致系统刚度下降，建议控制环流在额定电流的10%-15%为宜。

4. 电容电压平衡策略

4.1 全局排序法优化

传统冒泡排序时间复杂度O(N²)，对于N>100的系统难以实时实现。改进方案：

matlab复制[~, idx] = sort(cap_voltage, 'descend');
active_SMs = idx(1:n_actual);

% 增加电压差滞环
delta_V = 0.02 * Vdc/N;  % 2%滞环宽度
for i = 1:N-1
    if abs(cap_voltage(i) - cap_voltage(i+1)) < delta_V
        % 保持原排序优先级
        active_SMs = previous_active_SMs;
        break;
    end
end

实测性能对比（N=200）：

4.2 分组均衡策略

对于大规模MMC（N>500），建议采用分级控制：

1. 第一级：桥臂间均衡（调节有功功率分配）
2. 第二级：桥臂内分组均衡（每10个SM为一组）
3. 第三级：组内单个SM均衡

实现代码框架：

matlab复制function update_SM_status()
    % 桥臂间均衡
    P_avg = mean(P_arm);
    for k = 1:6
        delta_P(k) = Kp_arm*(P_avg - P_arm(k));
    end
    
    % 分组均衡
    for m = 1:num_groups
        [~, idx] = sort(group_voltage(m,:));
        active_SMs(m) = idx(end);
    end
end

5. 双闭环控制设计

5.1 电压-电流解耦控制

dq坐标系下的解耦方程：

code复制v_d = (R + Ls)i_d - ωLi_q + e_d
v_q = (R + Ls)i_q + ωLi_d + e_q

MATLAB实现关键代码：

matlab复制% 电流内环
Kp_i = 0.5; Ki_i = 20;
v_d_ref = Kp_i*(id_ref - id) + Ki_i*integral(id_ref - id) - omega*L*iq;
v_q_ref = Kp_i*(iq_ref - iq) + Ki_i*integral(iq_ref - iq) + omega*L*id;

% 电压外环
Kp_v = 0.1; Ki_v = 5;
id_ref = Kp_v*(Vdc_ref - Vdc) + Ki_v*integral(Vdc_ref - Vdc);

5.2 参数辨识技巧

1. 电抗器参数辨识：

   • 施加10%额定电压阶跃
   • 记录电流上升曲线
   • 根据τ=L/R计算参数

2. 耦合系数校准：

   matlab复制% 测试步骤
   step_id = 0.1*rated_current;
   measure_oscillation = max(iq) - min(iq);
   L_actual = L_nominal * (step_id/measure_oscillation);
   

5.3 动态性能测试

某±200kV MMC测试数据：

调试心得：先整定电流环带宽为1/10开关频率，再设电压环带宽为电流环的1/5-1/10。动态测试时建议从20%额定工况逐步提升，避免直接满负荷冲击。