三相MMC简化模型搭建实战：基于Sysplorer的建模与调试全攻略

搞MMC也有一阵子了后台经常有人问我“想入门MMC能不能用仿真软件搭个三相模型跑起来”这里面的坑我基本都踩过。今天这篇就用Sysplorer作为工具讲清楚怎么从零搭一个能跑的三相MMC简化模型把核心思路、参数选取、调试过程一次说透。不管是刚接触柔性直流的学生还是想转做电力电子仿真工程师的朋友照着这个思路搭一遍基本就能把MMC的运行机理摸个大概。1. 先搞懂MMC的底层逻辑子模块、桥臂和相单元模块化多电平变流器MMC这个名字拆开看就很好理解“模块化”说的是它由成百上千个相同的子模块堆叠而成“多电平”则是这些子模块协同工作后在交流侧输出的阶梯波电压。在动手建模之前这三点必须搞清楚子模块长什么样、桥臂怎么组织、三相系统怎么靠六个桥臂撑起来。1.1 子模块MMC的最小积木最常用的半桥子模块包含一个直流电容、两个IGBT开关管带反并联二极管和一个旁路晶闸管可选。工作时两个IGBT的开关组合决定电容是被“投入”到桥臂中还是被“旁路”。上管开通、下管关断电容接入桥臂子模块输出电压等于电容电压。上管关断、下管开通电容被旁路子模块输出电压接近0。上下管均关断此时处于闭锁状态电流只能走二极管这种状态一般只在故障或预充电时出现。所以子模块本质上就是一个可控的电压源只不过电压极性固定单极性幅值等于电容当前电压。整个MMC的输出电压调节靠的就是在任意时刻决定“桥臂上投入多少个子模块”然后用多个“0或Vc”的小台阶拼出期望电压。1.2 桥臂上桥臂与下桥臂的串联组合所谓“桥臂”就是一组子模块串联起来的支路。一个相单元由上桥臂和下桥臂组成两个桥臂的连接点就是交流输出端。直流母线侧分别连接到上桥臂顶端和下桥臂底端。假设每个桥臂有N个子模块上下桥臂各投入一部分子模块那么交流侧输出电压取决于上下桥臂投入数量的差值。同时直流侧电压需要支撑住所以上下桥臂投入的子模块数之和必须始终等于N。这里有个初学者很容易绕晕的点子模块电容电压不是恒定的。因为电流会不断给电容充放电所以电容电压一直在波动只不过设计时希望波动控制在额定值的一定比例内通常5%~10%。在仿真里如果看到电容电压波形有小幅锯齿波动那是正常现象别急着当故障处理。1.3 从单相到三相六个桥臂组成一个完整MMC三相MMC就是三个相单元并联在直流母线上每个相单元有上下两个桥臂。所以一共是六组桥臂每组桥臂有N个子模块。如果每个桥臂子模块数为N那么交流侧线电压的电平数可以达到2N1相电压的电平数是N1。这里有一个常见的搜索混淆点很多人搜“mmc排队模型中的等待的顾客总数是怎么计算的”那个MMC其实是排队论里的M/M/C模型C代表服务台数量跟咱们电力电子的模块化多电平变流器完全不是一回事。后面我会专门说这个。2. 为什么要把三相模型简化简化模型到底简化了啥既然MMC的子模块数量动辄几十上百个直接在仿真里搭详细模型会非常痛苦。比如一个桥臂20个子模块三相一共120个子模块每个子模块两个IGBT器件仿真步长要压到微秒级以下才能看清开关细节跑一个几十毫秒的工况可能需要几个小时。所以除非要做器件级损耗分析否则入门阶段一定要学会“简化”。2.1 详细模型的痛点开关频率高电磁暂态仿真步长必须小于开关周期否则波形失真。大量电力电子器件模型导致节点导纳矩阵规模爆炸每步求解都很慢。数值刚性严重遗忘电阻、寄生电容如果不设置好就发散。控制算法还没有调试通的时候用详细模型等于拿大炮打蚊子耗时且难以定位问题。2.2 三种常见的简化路线简化方式核心思路适用场景开关函数模型用理想开关代替IGBT忽略开关过程但保留子模块离散特性验证调制算法、子模块均压平均值模型把整个桥臂等效为一个受控电压源忽略单个子模块的开关动作系统级控制策略验证、功率传输特性分析戴维南等效模型用等效电阻串联电压源替代子模块串仿真精度高且速度快大规模MMC电磁暂态仿真初学者搭三相简化模型首推平均值模型。原因很简单我们想研究的是MMC在系统层面的响应比如直流电流、交流电压、功率调节、环流特性而不是某个IGBT的导通压降。平均值模型用受控源来模拟桥臂对外呈现的电压和电容储能效应仿真速度可以提升几个数量级。2.3 简化模型到底保留了哪些物理特性保留了子模块电容的储能效应所以能正确反映电容电压波动。保留了桥臂电感的续流和限流作用所以环流依然存在。保留了交直流侧功率平衡关系所以有功传输、无功调节的行为和详细模型一致。去掉了开关动作细节和器件损耗所以不需要关注死区、开关频率、谐波频谱的高频段。换句话说简化模型让你能用“宏观视角”看MMC的工作而不是陷在开关逻辑里。等你用简化模型把控制算法调顺了再回头换详细模型去做电磁兼容、绝缘协调、损耗评估那时候效率会高得多。3. Sysplorer建模实操从子模块到三相系统Sysplorer是一个支持多物理域建模的系统级仿真环境建模方式偏向物理拓扑连接很适合搭MMC这类功率变换电路。我用的版本搭建流程大概如下不同版本界面可能略有差异但核心逻辑是共通的。3.1 新建工程与元件库准备打开Sysplorer以后先新建一个“电路/电气”类型的工程。元件库一般包含基础被动元件电阻、电容、电感、受控源、电压表、电流表以及一些数学运算模块。如果找不到某个元件可以尝试在Modelica标准库中查找比如Modelica.Electrical.Analog下的理想受控电压源。本案例用到的核心元件清单直流电压源2个模拟直流母线正负端或者1个双极性直流源受控电压源6个对应六个桥臂的平均值模型电容6个对应六个桥臂的等效电容用来模拟子模块电容储能电感6个桥臂电感三相交流负载或无穷大电网用三相电压源代替接地参考点3.2 搭建单个桥臂的平均值模型平均值模型的核心是把“N个子模块串联”等效成“一个受控电压源 一个等效电容”。等效电容值等于单个子模块电容除以N因为串联后容值会缩小。具体的模型结构分两步实现第一步桥臂对外输出电压等于“该桥臂当前投入的子模块数 × 子模块电容电压”。在平均值模型中这个电压可以用一个受控电压源来表示控制量来自调制环节计算出的“投入子模块数量参考值”n_ref和当前电容电压vc。第二步电容电压怎么变化桥臂电流流过等效电容产生充放电。用公式表达就是dv_c/dt i_arm / C_equiv其中i_arm是流经桥臂的电流C_equiv C_sm / N。在仿真里用积分器对桥臂电流除以等效电容求积分就能得到电容电压。注意初始值要设成额定电容电压否则启动瞬间会有很大的暂态冲击。受控电压源的正负极性要严格和桥臂电流方向配合否则能量方向反了电容电压会越积越高或者一泄到底。我第一次搭的时候就因为参考方向没注意电容电压直接冲到5倍额定值仿真直接挂掉。3.3 搭建三相MMC整体拓扑三相MMC的电气拓扑连接如下每一相有上桥臂和下桥臂上桥臂顶端接直流正极下桥臂底端接直流负极。上桥臂和下桥臂的连接点作为三相交流输出端。六个桥臂的电感分别位于桥臂内部通常放在子模块串和输出端之间用于抑制环流和限制故障电流。在Sysplorer里连接电路时要注意直流侧必须提供一个参考地。通常把直流母线中点接地或者直接在下桥臂底端接参考地。如果忘记接地仿真器会提示孤立节点或者奇异矩阵这在我调试时出现过很多次。交流侧我用了一个三相电压源来模拟无穷大电网或者用三相RL负载也行。如果只是想看波形接一个纯电阻负载最容易。不过真实MMC一般通过变压器连接交流电网入门阶段直接用理想电压源更干净避免变压器参数引入额外变量。3.4 控制部分先开环跑通再说闭环新手最容易犯的错误是一上来就写闭环控制。我建议第一步先用最简单的开环调制也就是人为给定一个正弦调制波让上下桥臂的投入子模块数量按正弦规律变化。上下桥臂的投入数逻辑如下上桥臂投入数n_up (N/2) * (1 - m*sin(ωt))下桥臂投入数n_down (N/2) * (1 m*sin(ωt))其中m是调制比一般取0.8~0.9N是每个桥臂的子模块数。这样上下桥臂投入数之和恒等于N直流侧电压稳定交流侧输出就能得到一个接近正弦的阶梯波。在Sysplorer里需要生成正弦函数、加常数、乘系数这个可以用数学表达式模块实现。由于是平均值模型不需要真正的PWM调制器直接用连续的投入数参考值施加给受控电压源就行。开环跑通以后波形应该能看到交流侧电压为正弦桥臂电流中带有少量二倍频环流分量。这时候再考虑加环流抑制或者闭环功率控制否则后续控制参数调起来头疼。3.5 关键参数计算过程参数设计是整个模型能否跑起来的关键。我以一个10kV直流母线、1MW输出功率的简化三相MMC为例给出完整计算过程。直流母线电压Vdc 10kV桥臂子模块数取N 10则每个子模块额定电容电压Vc_nom Vdc / N 1000V注意这里没有考虑冗余实际工程中会取Vdc/(N-冗余数)但简化模型不需要管这个。交流侧相电压幅值调制比m 0.85则相电压峰值V_ph_peak m * Vdc / 2 0.85 * 5000 4250V对应相电压有效值约为3005V线电压有效值约为5205V。如果要接一个10kV交流电网需要配变压器但入门模型可以用理想电源直接给定幅值。桥臂电感工程上桥臂电抗标幺值一般在0.1~0.2 pu之间。取0.15 pu基准阻抗Z_base Vdc^2 / P 10000^2 / 1e6 100Ω桥臂电感L_arm 0.15 * Z_base / (2π * 50) 0.15 * 100 / 314.16 ≈ 0.0478H取47.8mH。子模块电容MMC的设计中电容储能与额定功率的比值一般取40kJ/MWh左右不对是kJ/MVA。通常取电容储能W 40kJ / MVA左右的水平。以1MW计算总储能 W_total 40kJ。因为每个桥臂有10个子模块六个桥臂共60个子模块每个子模块储能W_sm 40000 / 60 ≈ 666.7J子模块电容C_sm 2 * W_sm / Vc_nom^2 2 * 666.7 / (1000^2) 1.3334e-3 F即1333μF。这个值看上去比实际工程中的电容要大一些是因为我们取了较小的直流电压和较高的储能系数。实际以你的仿真需求为准。等效电容每个桥臂等效电容C_equiv C_sm / N 1333 / 10 133.3μF仿真步长平均值模型没有高频开关步长可以放宽到10~50微秒。如果是详细模型肯定要压到1微秒以下。最终参数表参数数值直流母线电压10kV桥臂子模块数10子模块电容电压1000V桥臂电感47.8mH子模块电容1333μF桥臂等效电容133.3μF调制比0.85载波频率若使用PWM2000Hz3.6 在Sysplorer中连接信号并运行将六个桥臂的受控电压源接入拓扑后需要在每个桥臂的等效电容积分器上设置初始值。上桥臂和下桥臂的电容初始值都设为1000V否则启动瞬间会有一个很长的电容充电过程波形看起来像是震荡衰减容易让人误判模型有问题。运行仿真前至少检查这几项所有受控电压源的控制输入是否来自正确的信号线。桥臂电流变量名是否一致积分模块有没有代数环。直流电压源的正负极方向和电容极性是否一致。是否有孤立节点或悬空节点。第一次运行时建议把仿真时长设短一点比如0.5秒步长50微秒。如果波形是发散的先暂停查看是哪个节点电压或电流越界。我习惯把每个桥臂电流和电容电压都加到观测窗口这样定位问题快得多。4. 调试之路仿真发散、波形畸变和环流问题让模型“跑起来”容易让模型“跑得对”难。我把自己调试过程中遇到的高频问题整理了一下基本都是平均值模型阶段一定会遇到的。4.1 一运行就发散的几个原因和排查顺序如果仿真刚开始几个毫秒就报错或数值爆炸按照这个顺序排查第一步检查积分器初始值。电容电压初始值为0而直流侧直接加10kV的话相当于瞬间给电容短路充电电流巨大仿真必崩。把等效电容初始值设为Vdc/N。第二步检查受控电压源的极性。上桥臂的受控电压源应该从上往下“降压”下桥臂从下往上“升压”。如果极性反了等效电容会被正向充电到超过额定值然后恶性循环。第三步检查是否有代数环。平均值模型里桥臂电压依赖电流电流又依赖电压如果信号连接方式不对仿真器可能无法求解。解决办法是在控制链路上加一个小惯性环节比如1/(Ts1)或者用Sysplorer的延迟模块断开瞬时反馈。第四步检查步长。平均值模型不是一定能用大步长如果电路中存在很小时时间常数比如很小的电容、很小的电感步长还是得相应缩小。我把50μs改成10μs以后很多莫名发散都消失了。4.2 交流侧电压波形畸变如果波形能跑但交流电压不是平滑正弦而是有明显畸变重点检查调制比和参考信号。调制比超过1时输出电压幅值会被直流电压限制产生削顶波形平顶。把调制比降到0.9以下。正弦参考信号如果带有直流偏置会导致上下桥臂投入数不对称交流侧出现直流分量。用数学模块减去平均值。桥臂电感太小电流纹波会增大在交流侧表现为谐波增加尤其低次谐波。适当增加电感。另外还有一个容易被忽略的点平均值模型输出的电压是瞬时值如果受控电压源的延迟太大会造成相位滞后看上去像波形畸变。此时检查模型里是否有不必要的滤波环节。4.3 环流问题二倍频负序环流MMC的环流是它的“专利烦恼”。环流主要包含二倍频负序分量它在三相之间流动不会流到交流侧和直流侧但会增大桥臂电流应力。在平均值模型中只要是用两个独立受控电压源模拟上下桥臂环流通道天然存在所以一定能观察到环流。环流偏大的表现是桥臂电流波形不再像正弦而是叠加了一个明显的二倍频波动。抑制环流最简单的方法是“加环流阻抗”也就是增大桥臂电感。入门阶段可以在上下桥臂之间或者三相中点之间加一个谐振支路但这会影响模型复杂度。更实用的是加环流抑制控制器检测三相桥臂电流的二倍频分量通过一个比例谐振控制器生成补偿电压叠加到调制波上。这个属于进阶内容但在简化模型上调通环流抑制之后再迁移到详细模型会非常顺畅。我调的环流抑制器参数一般是比例系数Kp取0.5~2谐振系数Kr取50~200谐振频率设为100Hz二倍基频。注意先开环看环流极性如果反馈接反环流不但不降反而会震荡起来。4.4 常见问题速查表现象可能原因解决方案仿真开始瞬间崩溃电容初始值为0设置电容初始电压为Vdc/N电容电压持续上升受控源极性反了交换受控电压源正负端交流波形削顶调制比太大m降至0.9以下桥臂电流低频震荡代数环或延迟过大加入小惯性环节缩小步长输出功率达不到预期调制比或电容电压偏低检查直流电压是否等于N*Vc_nom直流侧电流脉动大桥臂电感太小增加L_arm至0.15~0.2pu三相电压不对称三相参考信号相位不正确检查ABC三相相位差120°5. 快问快答单相MMC、PSCAD和“MMC排队模型”的区别模型跑通之后很多人接下来会有几个衍生问题这里一并说明白。5.1 单相MMC和三相MMC建模有什么不同单相MMC只有两个桥臂拓扑比三相简单很多适合用来验证子模块均压、调制策略等基础概念。把三相模型中取一个相单元断开另外两相就能得到一个单相MMC模型。但单相模型看不到环流在三相之间的流动路径因为环流本质上是三相之间循环流动的电流。如果你想研究二倍频负序环流必须搭三相模型。入门时先搭单相能快速理解子模块电压和输出波形的关系但想深入MMC运行特性三相简化模型是值得直接上的。5.2 为什么我选了Sysplorer而不是PSCADPSCAD/EMTDC在电力系统电磁暂态仿真中确实很经典但上手门槛不低而且网上的盗版问题不提也罢。Sysplorer这种系统级建模工具的好处是元件连接和信号控制更灵活尤其是和控制系统协同仿真的时候不需要像PSCAD那样在自定义模型上写大量Fortran代码。对于MMC这种“强耦合电气控制系统”的对象用Sysplorer搭平均值模型非常顺手改参数、加控制都更直观。当然如果之后要做绝缘配合、雷击过电压、断路器开断这类电磁暂态分析PSCAD依然是扎实的选择。工具没有绝对优劣关键是哪个能把你的问题快速跑清楚。我的建议是入门和算法验证用Sysplorer工程认证和电磁暂态分析再转PSCAD两者搭配效率最高。5.3 搜“MMC排队模型”到底是在搜什么这个属于搜索关键词的经典混淆。如果你搜“mmc排队模型中的等待的顾客总数怎么计算的”你看到的MMC其实是排队论中的“M/M/C”模型其中M表示到达间隔和服务时间服从负指数分布C表示服务台数量。那个等待顾客总数的计算是基于生灭过程和Little定律跟咱们的模块化多电平变流器没有任何关系。在写技术文章、做文档标题时建议把“MMC变流器”写全避免歧义。如果你真的是在查排队论的M/M/C模型记住基础公式等待顾客总数 (ρ^(C1)) / (C! * (1-ρ)) * P0其中 ρ λ/(Cμ)P0是系统空闲概率。这个和MMC变流器的电容电压公式一样都属于“各有各的MMC别搞混”的情况。5.4 简化模型跑通以后下一步怎么加详细功能三相简化模型稳定运行后可以按这个顺序逐步增加复杂度加入子模块电容电压排序均压算法观察电容电压均衡效果。换成最近电平逼近NLM或载波移相调制把平均值模型分解成离散子模块。加入二倍频环流抑制控制器观察桥臂电流波形改善。加入功率外环、电流内环让MMC能够控制有功和无功。接入变压器与交流电网模型做并网工况仿真。每加一步都在上一版模型基础上迭代不要一口气全改完否则出了bug都找不到源头。最后分享一个我自己踩过很多次的坑平均值模型虽然跑得快但它的数值特性比详细模型“硬”很多特别容易在积分环节产生数值漂移。所以我在所有电容电压积分器后面都加了一个很小的泄漏电阻或者用带复位功能的积分模块每隔一段时间强制校准一次初始值。这样做之后长时间仿真漂移的问题基本消失了。如果你在做MMC仿真的过程中也遇到电容电压慢慢爬升或者缓慢下降的情况不妨先试试这个土办法比调控制参数还管用。