基于Simulink的11电平MMC并网控制模型搭建与仿真调试

在Simulink里把11电平三相MMC逆变器并网控制模型完整跑通，我前前后后折腾了一个多月。第一版模型搭得很快，波形出来也“像那么回事”，但一查相电压波形，根本不是11级阶梯，子模块电容电压乱跳，桥臂电流里二倍频分量大得夸张。后来发现，问题不在Simulink本身，而在自己对MMC的调制、均压、环流抑制和并网控制这几个环节的理解太“线性”了——总以为把两电平那套dq控制直接搬过来就行，实际上MMC多出来的不只是电平数，而是一整套能量分配逻辑。

这篇文章就按我自己从零搭模型的实际顺序来写：先说明为什么选11电平这个规格，再逐块拆解子模块均压与调制、环流抑制、并网控制策略，最后把Simulink建模细节和调试经验完整交代出来。适合正在做储能并网、光伏中高压并网，或者准备用仿真验证MMC控制算法的人。即使你之前只玩过两电平三相VSC，沿着这条路径也能把MMC的坑基本摸一遍。

1. 为什么是“11电平”：MMC选型背后的工程逻辑

1.1 中高压并网场景下，两电平和三电平的尴尬位置

先聊一个听起来很基础但值得说透的问题：为什么要用MMC做并网逆变器？如果把逆变器放在10kV甚至35kV的中高压并网场合，两电平VSC会遇到两个硬约束。第一个是器件耐压，常用的IGBT耐压等级就那么多，两电平拓扑里每个开关器件要承受完整的直流母线电压，直流母线到了20kV以上，要么选器件串联，要么换更高耐压的器件，串联均压本身就是一个非常棘手的问题。第二个是波形质量，两电平开关频率做得再高，输出电压还是在两个电平之间跳变，dv/dt大、共模干扰重，并网滤波器体积跟着膨胀，这在储能站、光伏电站这种需要“安静、长寿”的场合并不友好。

三电平NPC能缓解一部分问题，但它的器件结构和控制复杂度是跳跃式上升的，中点电位平衡问题在动态工况下够调一阵子，而且电平数到这里基本到头了，想往更多电平走，NPC的钳位二极管和电容数量完全不可控。

MMC的思路则完全不一样：它不追求让单个器件承受高压，而是把高压拆成很多个低压子模块串联，每个子模块就是一个半桥加一个电容，承受的电压只有直流母线的N分之一。电平数上去之后，输出电压谐波大幅下降，单个器件开关频率也可以压得很低，开关损耗和滤波器体积都跟着降下来。这也是MMC在柔性直流输电、中高压储能PCS、STATCOM这些场景里能站住脚的底层原因。

1.2 11电平到底意味着多大的模型规模

先把“电平数”这件事拧清楚。对三相MMC来说，每相有上下两个桥臂，每个桥臂由N个子模块串联。输出相电压相对于直流母线中点，可以取N+1个电平。所以标题里的11电平，翻译过来就是每桥臂N=10个子模块，上桥臂10个、下桥臂10个，单相一共20个，三相一共60个半桥子模块。每个半桥子模块又含2个IGBT加反并联二极管、1个直流电容，所以整机里半导体器件数量是120个IGBT加120个二极管——数字听起来唬人，但每个器件只承受Vdc/10的电压，比如直流母线20kV时每个子模块电容只有2kV左右，器件选型反而宽松。

一个容易混淆的点是：这里说的11电平是相电压对直流中点的电平数，如果看线电压，因为两个相电压相减，台阶数会更多。和实时热搜里经常出现的三电平逆变器放在一起对比，就很直观。

实际做仿真时，器件数量的增加直接影响的是计算耗时。60个子模块如果全部用电力电子开关器件级建模，Simulink每个步长都要解一大堆开关节点，跑一次几百毫秒的工况能急死人。所以本文后面的建模方案会特意避开“器件级”建模，这个后面展开。

1.3 为什么先从11电平上手，而不是21电平、51电平

我的看法是，11电平是“刚好能覆盖MMC全部核心问题”的最小规格。电平数再少，比如5电平，子模块电容电压波动的问题也有，但环流和均压之间的耦合效应不够明显，很容易让人忽略；电平数太高，比如21电平甚至更高，模型规模成倍上涨，仿真时间实在不划算，而且控制算法本身还没有本质变化。更现实的是，11电平的相电压阶梯已经足够逼近正弦，你可以在示波器里很直观地看到“阶梯波叠加PWM”的形态，这对理解MMC工作原理非常友好。等11电平玩明白了，把N改成20、30只是改参数的事，新增的问题主要是开关状态组合爆炸和排序均压的计算负担，而不是控制原理变了。

2. 子模块均压与调制策略：仿真中第一个绕不过的坎

2.1 桥臂功率脉动让电容电压天然“不稳”

MMC子模块电容的运行方式和两电平直流母线电容很不一样。两电平的母线电容主要吸收的是直流-交流功率差，频率高但幅度可控；MMC的每个子模块电容则要持续被桥臂电流充放电，而且桥臂电流本身是“直流分量+工频分量+环流分量”的混合体。

可以这样理解：子模块的开关函数近似为占空比信号，其中既有直流偏置也有基波调制分量，而桥臂电流同样含有直流和基波分量。两个含基波分量的信号相乘，就会产生二倍频的瞬时功率分量。也就是说，电容上的能量本来就在以电网两倍频率的节奏往复振荡。这是物理特性决定的，不可能完全消除，只能通过参数设计把电压波动幅度压在允许范围内，再通过均压控制让每个子模块的波动步调一致、不发散。

仿真中常见的错误是把电容电压波动直接归结为“均压没做好”，其实不是。就算均压做得完美，每个电容电压也都有±5%左右的周期性波动，区别在于波动是均匀还是个别模块特别大。这个认知能帮你少走很多弯路。

2.2 PS-PWM与排序均压：Simulink里的搭配选择

调制策略决定了怎么把子模块的投入组合出来，而均压策略决定了在每个控制周期让具体哪个子模块投入。这两个问题在MMC里是绑在一起解决的。

我的模型用的是载波移相PWM（PS-PWM），也就是每个桥臂的N个子模块各配一路三角载波，载波相位依次错开2π/N，共享同一路调制波。这样单个子模块的开关频率可以很低，但桥臂输出电压的等效开关频率是单个载波频率的N倍，波形谐波性能很好。载波移相PWM在Simulink里实现很直接：用N个Repeating Sequence生成三角波，或者用Clock做载波计数，相位依次加1/N个周期，然后和调制波做比较。

均压部分，我推荐排序法。思路很简单：每个控制周期采集桥臂电流的方向，判断当前是在给子模块电容充电还是放电；然后把本桥臂N个子模块的电容电压排序，充电时优先投入电压最低的几个，放电时优先投入电压最高的几个。这样电压高的少充多放，电压低的多充少放，自然收敛。逻辑不复杂，但要注意排序法和载波移相配合时，子模块的开关状态可能频繁翻转，实际工程上会在排序结果里加滞环或者限制开关频率，仿真里为了提速可以先不管这些。

和其他调制策略放在一起看，选型逻辑会更清楚：

2.3 先算对参数再开搭：电容和桥臂电感的数量级

很多仿真模型发散，首要原因不是控制参数，而是主电路参数拍脑袋定的。这里给出一套10kV/10MVA并网算例的参数推导过程，照这个数量级去搭，配合控制参数调整会顺利很多。

先确定规模：直流母线电压Vdc取20kV，每桥臂N=10，子模块电容额定电压就是Vc=Vdc/N=2000V。系统额定容量Sn=10MVA。

子模块电容用一个常用的储能常数法估算：C约等于2·Sn·H/(3·N·Vc²)。H是储能时间常数，工程上按40 kJ/MVA左右取值。代入得C≈2×10×10⁶×0.04/(3×10×2000²)≈6.67mF，考虑裕量我取了8mF。这个值决定了电容电压波动幅度，8mF对应约±5%的波动，仿真波形基本符合预期。如果想看电容电压波动对控制的影响，可以试着把电容减半到4mF，波动会明显加大，控制器更容易表现出耦合问题。

桥臂电感按标幺值设计。取桥臂电感0.1pu，交流侧基值Lbase=Vg²/(ω·Sn)=100×10⁶/(314×10×10⁶)≈31.8mH，所以L_arm≈3.18mH。折算到交流侧的等效相电感约为L_arm/2≈1.6mH，这个值后面做电流内环会用到。桥臂电感不要取得太小，否则环流抑制压力很大；也不要超过0.15pu，否则动态响应变差，而且电感上的基波压降会侵蚀可用调制范围。

3. 环流怎么来的，怎么在模型中压下去

3.1 二倍频环流的形成机理与波形识别

环流是MMC一个非常“有辨识度”的问题。它本质上是在直流母线正负极之间、通过三相桥臂形成内部回路流动的电流，不经过外部电网，因此在外部端口是看不到的。你只能在桥臂电流里发现异常：桥臂电流除了直流分量和基波分量之外，还有明显的二倍频成分，频率在50Hz系统里就是100Hz。

产生机理还是回到功率波动。每相桥臂的瞬时功率是脉动的，三相加起来大致恒定，但单相内部总有能量从直流侧流进桥臂电容再流出的过程。由于三相桥臂的功率脉动存在相位差，这些内部能量就顺着三相桥臂和直流母线之间循环流动，形成二倍频环流，而且三相合起来呈现负序特性的旋转对称关系。环流幅值会随着传输功率增大而增大，严重时能把桥臂电流的峰值顶高百分之二三十，器件应力和损耗都跟着涨。

仿真里怎么识别？把上、下桥臂电流相加再除以2，减去理论上的直流电流分量（约等于直流母线电流的三分之一），剩下的就是环流。如果这个量里能看到清晰稳定的100Hz正弦波，幅度远超噪声水平，环流抑制就该出手了。

3.2 三种环流抑制方案的仿真实现对比

抑制环流常见有三条路，我在仿真里都试过，各自特点比较鲜明。

第一种是单纯加大桥臂电感。电感对二倍频环流有天然的阻抗作用，效果直观，但代价也直观：电感体积、成本、基波压降全上去了，而且电感过大还会引发和子模块电容之间的谐振风险。所以这个方法只适合作为辅助，不能作为主要手段。

第二种是在二倍频旋转坐标系下做PI控制，也就是经典的CCSC方法。思路是把三相环流分量做一个坐标变换，用2ω角速度旋转的坐标系把它变成dq直流量，然后用PI控制器把d、q轴误差压到零，输出再反变换成三相加到桥臂调制波上。这个方法物理意义清楚，PI参数调起来有规律，但坐标变换和正负序方向容易搞混，我在这上面吃过亏。

第三种是用静止坐标系下的PR谐振控制器，谐振频率对准100Hz。控制器形式一般是kP+kR·s/(s²+ω₀²)，在100Hz处提供高增益，把环流的二倍频分量压下去，对其他频段影响很小。仿真里这个方案最省事，不用变换坐标，直接检测环流加减到调制波上就行。我的模型最终用的就是PR方案，谐振增益取比例增益的10倍左右，中心频率按电网频率实时自适应会更好，固定100Hz也无伤大雅。

3.3 环流抑制和均压策略的耦合处理

这里要提醒一个容易忽略的细节：均压排序和环流抑制不是完全独立的。排序法频繁改变子模块投入状态，本身就会给环流注入扰动；而环流抑制器输出的是共模性质电压，作用在上下桥臂上时又会轻微改变调制比和子模块充放电分布。两者同时打开以后，如果参数都调得很激进，可能会出现“均压把环流压住，环流又让均压发抖”的奇怪振荡。

我的处理方法是分步调试。先把均压排序单独开起来，观察电容电压收敛；再把环流抑制单独开起来，看100Hz分量是否被压下去；最后两个同时开，观察波形是否保持稳定。如果出现耦合振荡，优先把环流抑制的谐振增益降一点，通常比动均压参数更有效。这个顺序看着笨，但确实能少掉很多头发。

4. 并网侧控制：锁相、dq电流环与功率外环的搭法

4.1 SRF-PLL锁相环节的搭建与参数

并网控制的第一步是拿到电网电压的幅值、频率和相位，最常用的是同步旋转坐标系锁相环（SRF-PLL）。原理不复杂：把三相电网电压先做Clark变换到αβ，再用PLL输出的角度θ做Park变换得到dq分量。如果θ对准了电网电压矢量，q轴分量应该等于0；只要有偏差，q轴就会有一个误差信号，用PI控制器把这个误差拉回零，PI输出就是角频率偏差，积分后得到θ。

仿真搭建时要注意带宽整定。PLL带宽太低，锁相过程慢，并网瞬间动态差；PLL带宽太高，电网电压里的背景谐波会污染相位，dq变换出来的角度都在抖动。我一般按带宽10~30Hz、阻尼比0.707来设计PI参数，网压畸变严重时往低取。Simulink里各步骤都是现成模块：Clark可以写一个3×2的增益矩阵，Park用Fcn或者查表做三角函数，PI限幅要做，积分器初值可以设电网额定角频率，这样启动阶段收敛快一点。

4.2 dq轴电流内环：对象、PI参数与解耦前馈

电流内环的被控对象不能直接套两电平的公式。MMC的桥臂电感上下并联，折算到交流侧的等效电感约为L_arm/2，再考虑交流侧如果加了滤波电感就一起加上。在这个基础上，dq轴电流方程里存在交叉耦合项ωL·iq和ωL·id，必须用前馈补偿项把这个耦合拆开，不然d轴指令一变，q轴电流会跟着晃。

电流环PI的整定我喜欢按带宽法来。给定电流环带宽f_bw，大概取1/5到1/10的等效开关频率。载波频率1.5kHz时，等效开关频率接近15kHz，但实际控制周期有限，带宽取300~500Hz比较稳妥。然后根据等效电感L_eq和等效电阻R_eq，算出kp≈L_eq×2πf_bw，ki≈kp×R_eq/L_eq，或者直接按时间常数来。实际上这个公式出来是一个比较合理的初值，仿真里微调10%~20%就能出不错的波形。

Simulink里实现时，控制器用离散化比较省事，采样周期取100μs左右即可。PI块建议用Tustin法离散，避免前向欧拉在大增益时产生极限环。

4.3 直流电压外环与无功通道

电流内环给了快速的有功/无功电流控制能力，但电流指令谁给出？有功电流指令一般来自直流电压外环。MMC的直流侧等效电容是几十个子模块电容串联的结果，动态特性比两电平直流母线慢很多，所以直流电压外环的带宽必须显著低于电流环。我习惯把外环带宽取电流环的1/5以下，PI参数按同一套带宽法算，出来的效果是直流电压缓慢而稳定地跟随给定，不出现振荡。

无功电流指令在仿真里可以直接设0，实现单位功率因数并网。如果想看动态，也可以给一个无功阶跃，观察并网电流的无功分量变化。这里需要注意符号约定：常见的电网电压定向坐标系里，有功功率P=1.5·ud·id，无功功率Q=-1.5·ud·iq，iq方向定义不同会差一个负号，调试前先把坐标变换方向确认清楚。

4.4 控制输出如何折算到上下桥臂调制波

控制器算出来的是三相交流侧的相电压参考值v_phase。MMC里这个参考值不能直接当调制波用，因为上下桥臂各有各的直流偏置。我用的分配关系是：

• 上桥臂每个子模块的调制波：v_up_j = (Vdc/2 - v_phase_ref + v_cir_comp) / N
• 下桥臂每个子模块的调制波：v_down_j = (Vdc/2 + v_phase_ref + v_cir_comp) / N

这里Vdc/2/N是子模块的直流电压偏置，v_cir_comp是环流抑制器输出的共模补偿电压，上下桥臂加的是同一个量，不会影响交流输出电压。得到的调制波要除以子模块电容额定电压Vc做归一化，落在0~1区间内参与PWM比较。这个接口是最容易出错的地方，我见过不少模型在这里符号反了，输出相电压直接少一半或者多一堆直流偏置，波形怎么调都调不对。

5. Simulink模型从零搭起的实操记录

5.1 模型架构与交付物的层次

一个能跑顺的MMC并网模型，结构上一定要分层，不然调试时你会迷失在几十个子模块的连线里。我的顶层模型分成五个子系统：三相电网与并网接口、MMC主电路、调制与均压控制、并网控制器（PLL+外环+内环）、环流抑制器。

网侧用一个理想三相电压源串线路阻抗即可，并网点加一个断路器，模拟0.2s左右合闸并网。主电路部分不推荐用Simscape自带的IGBT器件模型去搭60个子模块，仿真代价太大；我采用的是理想开关函数模型，也就是每个子模块等效成一个受控电压源，输出电压等于开关信号乘以电容电压，电容电压则通过对桥臂电流做积分获得。这样既保留了子模块电容储能、电压波动和环流这些核心物理特性，又避开了器件级求解的负担，仿真速度可以接受。

每个桥臂封装成一个子系统，里面有10个子模块子系统和桥臂电感。子模块建议用Mask封装，把C、Vc0作为参数暴露出来，这样后续改电平数也就是改参数加复制模块的事。

5.2 子模块理想开关函数模型的落地

子模块理想模型的数学关系非常简洁：

• 开关函数s=1表示子模块投入，电容接入桥臂；s=0表示旁路，电容被短路。
• 子模块端口输出电压v_sm = s × vc
• 电容电流ic = s × i_arm
• 电容电压vc = 1/C ∫ ic dt，初始值设为Vdc/N

在Simulink里，这个模型可以写成离散化的Matlab Function，也可以用积分器加逻辑运算搭。核心逻辑如下：

function [v_sm, vc_next] = sm_ideal(s, i_arm, vc, C, Ts) % 理想半桥子模块模型 % s: 开关信号，1为投入，0为旁路 % i_arm: 桥臂电流，从Simscape电流传感器获取 % vc: 上一时刻电容电压，初始外部设为Vdc/N ic = s * i_arm; vc_next = vc + (1/C) * ic * Ts; % 前向欧拉离散积分 v_sm = s * vc; end

这里有个关键坑：如果直接用Controlled Voltage Source的输入输出和电流测量形成闭环，很容易产生Simulink代数环，模型跑不动或者结果明显不对。解决方法是把开关信号s通过一个Memory或者Unit Delay延迟一个仿真步长再送入子模块，相当于PWM在执行上有一拍延迟。这个延迟在工程里本来就有，仿真中接受它即可。

5.3 启动、步长、代数环与数值发散

仿真参数这块，我的配置是定步长离散求解器，步长取1e-5秒，载波频率1.5kHz时一个载波周期有约66个仿真步长，波形解析度足够。定步长主要理由是PWM比较信号在每个步长的确定性更好，变步长求解器遇到开关时刻会频繁细化步长，仿真慢且容易触发数值噪音。

启动问题比参数更容易翻车。子模块电容初始电压如果设成0，模型一开始就是一个低压电容串联进高压回路，充电冲击电流能冲到额定值的几十倍，波形完全失真。正确做法是把所有子模块电容初始电压设为Vdc/N，让模型在启动瞬间就处在正确的直流偏置点附近；并网断路器再延迟几十毫秒合闸，进一步避免电网和主电路之间的初始相位冲击。

数值发散还有一个常见来源：控制器输出限幅没做好。调制度一旦超过1，PWM进入过调制区，输出波形出现大量低次谐波，电流激增，整个模型发散。调制波上限一定要钳位在0.98左右，留一点死区，别等到波形飞了再查。

5.4 波形验证清单

模型跑通之后，我用一张检查清单逐项核对，避免“看起来在动”就以为正确。这里给出正常和异常的特征对照，调试时非常实用。

6. 从仿真到实物之前，还有哪些坑要留个心眼

6.1 标幺值换算与样机参数映射

仿真模型用的是10kV/10MWA这个量级，做完之后想往实验室样机迁移，很多人直接按百分比缩小参数，最后发现环流抑制效果和电容电压波动完全对不上。正确做法是按标幺值换算：桥臂电感保持0.1pu左右，根据样机电压和容量用L=L_pu·V²/(ωS)反算实际亨利值；子模块电容保持相同的储能时间常数和电压波动率反算实际法拉值。电压等级变了，但标幺值体系不变，控制参数和波形形态才能在工程上延续。

6.2 器件数量、开关频率与冗余设计

11电平MMC用了120个IGBT，器件数量多，故障概率也随之上去了。工程上通常会把每桥臂子模块数做N+1甚至N+2，多出来的模块平时不参与调制，坏了之后通过旁路开关顶上去。仿真里也可以模拟这个场景：人为把一个子模块的开关信号改成旁路，观察均压和环流会不会突变，这是验证容错控制逻辑最直接的方式。

单个子模块的开关频率低不代表散热压力小。MMC的损耗分布和两电平不一样，导通关断损耗相对低，但电容的纹波电流和电感的铜损可能更突出。做热设计时别把注意力全放在IGBT上，子模块电容的纹波电流耐受能力往往是短板。

6.3 我踩过的三个典型坑

第一个坑是载波移相的相位没真正错开。我用了一个载波源复制出N路信号，忘了给每一路加相位偏移，结果所有子模块的PWM比较结果几乎同步变化，相电压波形直接退化成两电平。排查方式很粗暴：把N路载波信号画在同一张图上，看相位是否均匀铺满一个周期。载波移相的“移相”二字，落实在模型里就是每个子模块三角载波相位依次增加1/N周期。

第二个坑是排序均压的采样频率不够。我开始把均压计算放在和控制系统一样慢的100μs采样周期上，10个子模块的排序在一拍内完成，看似没问题，但实际效果很差，电容电压有高频分量逃过了采样。后来把均压的执行频率提到和PWM载波一致的周期，效果立刻改善。排序这种事，时间分辨率决定上限。

第三个坑是PR环流抑制器的离散化。直接在连续域设计好参数放到离散模型里，谐振点会漂移，100Hz的增益实际上落到90Hz或者110Hz，环流抑制效果打折，甚至引入新的低频振荡。解决办法是用Tustin变换加预畸变，或者直接在离散域设计谐振器，让中心频率精确对准100Hz。这个细节在文档里很少写，但仿真结果差之毫厘失之千里。

最后说一点个人体会。模型的每个数字都对应物理意义：载波相位错开，就意味着桥臂各子模块轮流分担电压；排序均压，就是让能量在时间和子模块之间重新分配；环流抑制，则是把三相之间来回倒腾的无功能量按回去。把这些想明白之后，看波形就不再是“像不像正弦”这么简单，你能从每一个纹波里读出它从哪来、过不过分。建议你把环流抑制这一路单独摘掉，跑一次对比，亲眼看100Hz分量在桥臂电流里涨起来的那一下，对MMC的理解会扎实很多。