基于VSG的T型三电平逆变器双机孤岛离网功率均分仿真

搞微电网仿真的人多半都有过这种经历单机 VSG 调得服服帖帖换到双机并联孤岛离网一合闸就出幺蛾子。环流过大、频率振荡、功率分不开波形乱七八糟。这个项目标题看着长其实核心就三件事虚拟同步发电机 VSG 控制、T 型三电平逆变器、两台机孤岛离网功率均分最后用 Simulink 把它跑通。这篇文章我不打算讲教科书式的理论铺垫而是按我实际搭模型、调参数、看波形的顺序把整个系统的控制链路、拓扑选择、均分机制、仿真实现和调试坑位全部拆开来讲。适合正在做微电网并联仿真、准备用 VSG 做离网控制、或者被三电平建模折腾过的朋友参考。1. 孤岛离网场景下逆变器为什么要假装成一台同步机1.1 传统下垂控制能分功率但撑不住频率冲击微电网离网运行时失去了大电网这个无限大母线撑腰所有逆变器必须自己建立电压和频率。最容易想到的控制就是下垂控制有功-频率下垂、无功-电压下垂。下垂控制的逻辑很直观输出功率大了就把频率和电压往下调一点相当于模拟同步机的一次调频特性。但这里有个先天问题下垂控制本质是比例控制没有惯量。它对应的是一个静态方程频率可以瞬间跳变到目标值。当负载突变、电机启动或者大功率设备投入时频率变化率 dω/dt 可以飙得很高这会让对频率敏感的负载出问题也会让并联系统在动态过程中振荡。你可以把下垂控制的逆变器看成一个没有飞轮的发电机风一吹转速就乱。VSG 干的事就是给逆变器装一个虚拟飞轮。它在控制算法里植入同步发电机的转子运动方程用软件模拟出一个物理惯量 J 和一个阻尼系数 Dp。这样一来逆变器在动态过程中就有了类似同步机的摇摆特性频率变化不再是瞬间完成而是按转子方程一步步走。1.2 VSG 的机械模拟惯量 J 和阻尼 D 从哪里来VSG 最核心的方程就两个转子运动方程和相位方程J * dω/dt Pref / ω0 - Pe / ω0 - Dp * (ω - ω0) dθ/dt ω其中 Pref 是有功功率参考值Pe 是实测输出有功功率ω0 是额定角频率J 是虚拟惯量Dp 是阻尼系数。你把右边三项翻译成人话就是机械转矩驱动转子加速减去电磁转矩阻碍转子加速减去阻尼转矩抵消频率偏差。功率参考值 Pref 和实测功率 Pe 的差值产生不平衡转矩推动虚拟转子加速或减速而阻尼项 Dp 的作用类似同步发电机的阻尼绕组频率跑偏时它往回拉。θ 则由 ω 积分得到它就是逆变器输出正弦参考电压的相位。在 Simulink 里实现这一步很简单用两个积分器搭一个二阶系统把 Pref、Pe、当前 ω 引进去输出 ω 和 θ。难的是 J 和 Dp 怎么取值这个放到后面参数设计部分详细讲。1.3 两台 VSG 并联都不做主、功率为什么能自动分单机 VSG 跑通之后你自然想知道两台 VSG 并联彼此没有通信线负载是怎么分到两台机头上的答案是通过公共母线频率这个无形的通信信号。设想一个 10kW 的负载突然接到公共母线上。在负载投入瞬间两台的输出功率暂时不够母线电压被拉低、频率有下降趋势。每台 VSG 都检测到自身输出功率 Pe 上升于是转子运动方程里不平衡转矩变大ω 开始下降。注意两台机共享同一条母线母线上只有一个频率。因此两台机的 ω 最终会被锁在同一个值上。稳态时 dω/dt 0转子运动方程变成Pe Pref - Dp * (ω - ω0)这就是 VSG 自带的频率下垂关系。只要两台机的 Pref 一样、Dp 一样稳态时它们输出的 Pe 就一样。有功功率自然均分不需要任何通信。这个特性是 VSG 在多机并联场景下最值钱的地方。你在波形里看到的频率变化过程其实就是虚拟转子在摇摆摇摆结束后功率达到新的平衡点。2. T 型三电平拓扑相比二极管钳位和两电平这台机器怎么选2.1 三种常用逆变器拓扑的关键差异标题里明确写了 T 型三电平这不是随便选的。做仿真时很多人图省事用两电平但如果在做中高压、大容量微电网项目三电平的波形质量和器件应力优势会直接决定方案是否可行。拓扑单个开关管电压应力输出电平数钳位方式主要短板两电平全直流母线电压2无谐波大、dv/dt 高、开关损耗大二极管钳位 NPC 三电平半直流母线电压3二极管钳位钳位二极管多、内外管损耗不均T 型三电平半直流母线电压3双向开关钳位中点电位需额外处理两电平的问题很好理解开关管要扛全母线电压开关频率一高损耗和尖峰都很痛苦。二极管钳位NPC是最早成熟的三电平方案但它每个相多了一串钳位二极管而且在实际运行中内管和外管损耗差异大散热设计不好做。T 型三电平的本质是两电平拓扑 双向钳位开关它用最少的器件达到三电平效果损耗分布更均匀在中低压场合比如 690V 以下非常有优势。2.2 T 型桥臂的电平输出逻辑和中性点电流T 型桥臂每个相有四个 IGBT结构可以理解成T1 和 T2 组成一个上下半桥T1 接母线正极T2 接母线负极T3 和 T4 组成一个双向钳位支路把输出端接到直流侧中点。输出逻辑一句话就能说清输出正电平 Vdc/2上半桥导通输出端接母线正极输出零电平 0中点双向开关导通输出端接直流中点输出负电平 -Vdc/2下半桥导通输出端接母线负极。注意 T3 和 T4 是背靠背反接的双向开关不是串联分压。两管同时导通时电流可以从输出端流进中点也可以从中点流出到输出端方向不同走不同的管和续流二极管。三电平相比两电平多出来的一个零电平让输出电压波形更接近正弦线电压谐波明显下降这是它最重要的收益。但代价是引入了中点电位问题当输出为零电平时负载电流流经直流中点如果上下电容充放电不平衡中点电位就会漂移导致某一半电容电压过高威胁器件安全。2.3 Simulink 建模的两种路线标准模块还是 IGBT 封装Simulink 里做 T 型三电平版本不同差别很大。新版本的 Simscape Electrical 元件库里有三电平变换器模块但要注意很多自带模块默认是二极管钳位 NPC 型不一定是你要的 T 型。如果你用的版本没有现成 T 型模块我的建议是直接手动搭。手动搭一个相桥臂其实不复杂四个 IGBT/Dio de 块按照 T 型拓扑连起来然后复制三份。直流侧用两个电容串联中点引出给双向钳位用。这种做法的好处有三个一是拓扑结构透明你可以随时测量任意开关管的电压电流二是换流路径完全可控NPC 模块里的内部钳位逻辑不会干扰你三是对照论文调参数时模型和理论能严格对得上。还有一个折中方案用 Universal Bridge 模块桥臂数选 3功率开关选 IGBT然后加一个二极管中点钳位。这种近似在只关心控制算法、不关心器件损耗时可以跑但波形细节上会有差别尤其是零电平换流时刻的电流路径。做双机并联功率均分这种控制层面的研究我更推荐手动搭建花半天时间一劳永逸。3. 双机功率均分的控制架构VSG 环、电压电流环、虚拟阻抗3.1 整个控制链路从功率计算到 PWM 的三段结构VSG 不是孤立的算法它必须嵌进逆变器的完整控制链里才起作用。我在仿真里用的控制结构是经典的三段式功率计算 - VSG 核心输出 θ 和 E - 电压外环 - 电流内环 - 调制 - PWM功率计算用瞬时功率法三相电压电流经过 Clarke 变换得到 αβ 分量然后算有功 P 和无功 Q。注意功率计算不能直接用瞬时值进 VSG因为瞬时功率里含有大量二倍频分量会让转子方程振荡。一般在功率后面加一个低通滤波器截止频率取基波频率的 1/10 到 1/5我用的是 100Hz 左右的惯性环节。VSG 核心输出两个量参考相位 θ 和电压幅值 E。θ 由转子方程积分得到E 则通过无功-电压调节得到。得到 θ 和 E 之后三相参考电压就是标准的正弦波ua E·sin(θ)ub E·sin(θ-120°)uc E·sin(θ120°)。这里有个容易忽略的地方VSG 的 E 对电压的控制是单向的它只能通过无功功率调节电压幅值但逆变器要真正输出这个电压还需要电压电流双环来跟踪。电压外环把参考电压和实际电压做差通过 PI 调节器给出电流内环的参考值电流内环再把电流参考和实测电流做差输出调制电压。双环的作用就是把 VSG 给出的理想电压指令转化成实际的 PWM 占空比。3.2 有功均分靠频率自同步无功均分靠什么有功均分的机制在前面已经说了靠的是公共母线频率和转子方程的下垂关系稳态时只要 Pref 和 Dp 对称P 自然分成一半。无功均分就没这么容易了。仔细想一下原因无功功率没有一个母线无功频率这样的全局同步信号。每台 VSG 只能测量自己输出端的无功功率然后通过 Q-V 下垂关系调整本地电压幅值。问题在于两台逆变器的线路阻抗不一样时同样的电压幅值差在不同阻抗上产生的无功功率差是完全不同的。这就是无功均分比有功均分难的根本原因。最常见的做法是给无功-电压控制加一个积分项或者用标准的 Q-V 下垂E E0 Kq * (Qref - Q)其中 Kq 是无功下垂系数Qref 是无功参考值。理论上Kq 相同、Qref 相同的两台机稳态时 Q 也应该均分。但实际仿真里你会发现因为线路阻抗差异和电压测量误差Q 总会有偏差。要解决这个问题就得引入虚拟阻抗。3.3 线路阻抗不一致时虚拟阻抗怎么补虚拟阻抗的思想是我虽然改变不了物理线路但可以在控制算法里故意加一段阻抗让逆变器的等效输出阻抗变成标准值。实现方式是在电压参考指令上减去一个阻抗压降Uref* E - Zv * ioZv 是虚拟阻抗io 是输出电流。比如设置 Zv Rv jXv其中 Rv 取 1Ω、Xv 对应 1mH 电感的感抗那么无论物理线路阻抗是多少从 VSG 角度看每台机的外特性都趋向一致。虚拟阻抗对无功均分的作用特别明显。本身低压微电网线路是偏阻性的而传统下垂是基于感性线路设计的控制量和实际物理量对不上。加了虚拟阻抗后你可以把等效输出阻抗掰成感性Q-V 下垂才真正打准。仿真里调这个参数非常直观一开始两台机无功差 30%把虚拟电感加上去差值能缩到 5% 以内。3.4 J、D 值设计从阻尼比到参数整定VSG 参数整定是我这次仿真花时间最多的地方。先给结论J 和 Dp 不是随便选的它们决定了 VSG 的振荡频率和动态响应。把转子运动方程做小信号线性化后可以得到一个二阶系统。特征方程近似为J * s² Dp * s Pe_sync / ω0 0其中 Pe_sync 是同步功率系数稳态时近似等于额定功率和线路阻抗决定的量。二阶系统的自然振荡角频率 ωn 和阻尼比 ζ 分别为ωn sqrt(Pe_sync / (J * ω0)) ζ Dp / (2 * sqrt(J * Pe_sync / ω0))也就是说J 越大振荡越慢系统惯性越大Dp 越大阻尼越强振荡衰减越快。但 Dp 同时也会加大稳态频率偏差因为前面说过稳态时 Pe Pref - Dp·(ω - ω0)Dp 越大同样的功率偏差对应更大的频率偏移。实际整定时我喜欢用一组工程经验值起步先按惯性时间常数 H 选 J关系是 J 2H·S/ω0²。对于 10kW、50Hz 的系统ω0 314 rad/s取 H 2s则 J ≈ 0.4 kg·m²。Dp 的初值可以取 J 的 3 到 5 倍跑出来看频率波形再微调。记住一个原则先保证 Dp 足够大让系统不振荡再调 J 拉惯量最后回头微调 Dp 控制频率稳态偏差。4. 从 0 到 1 搭出双机并联仿真模型结构、参数与步骤4.1 主电路搭接直流源、T 型桥、滤波与公共母线我用的是两套完全对称的主电路每台逆变器一个独立直流源直流电压 Udc 700V配合两个串联电容模拟直流中点。逆变器输出经过 LC 滤波器接到一条公共母线上公共母线再接本地负载。滤波参数需要算一下不能拍脑袋。LC 滤波器的截止频率一般取开关频率的 1/10 左右。开关频率选了 10kHzLC 截止频率取 1kHz 附近。取 L 3mHC 20μF算出来截止频率约为 1.03kHz符合要求。这个参数下10kHz 的开关纹波基本被滤掉而 50Hz 基波分量几乎无损通过。公共母线的连线阻抗我故意设得不对称1 号机到母线阻抗是 0.1 j0.314Ω2 号机到母线阻抗是 0.2 j0.628Ω。在 50Hz 下这两个阻抗分别对应约 1mH 和 2mH 的感抗。不对称阻抗是专门用来考验无功均分能力的如果控制算法只靠下垂不做虚拟阻抗补偿结果一定会出偏差。4.2 控制算法分模块实现VSG、双环、PWM 生成的步骤在 Simulink 里我按子系统分块搭每个子系统的职责非常单一。控制部分从测量模块开始。测量三相电压、三相电流后先做 Clarke 变换abc - αβ再算瞬时功率得到 P 和 Q。P 和 Q 各自过一个一阶低通滤波滤掉二倍频分量再送进 VSG 模块。VSG 模块内部就是转子运动方程的两个积分器和无功-电压调节器。有功路径输出 θ无功路径输出 E合并生成三相参考电压。电压电流双环放在同一个子系统里。电流内环用 PI 调节器加前馈解耦补偿电感上的耦合分量电压外环也是 PI输出电流内环的 dq 轴参考值。双环输出的是 dq 轴调制电压经过反 Park 变换得到三相调制波。PWM 部分用双载波同相层叠的方式。每相有两个同相三角载波调制波为正时和上三角载波比较生成正电平脉冲调制波为负时和下三角载波比较生成负电平脉冲。两个比较信号组合成四路开关管的门极信号。仿真里还要留死区时间默认设 2μs。4.3 关键参数表照抄这份参数能少走很多弯路参数数值说明直流电压 Udc700V对应线电压 380V 等级的常见直流母线额定功率 S10kVA单台容量滤波电感 L3mHLC 截止频率约 1kHz滤波电容 C20μF与 L 配合构成滤波器开关频率 fs10kHz双载波同相层叠VSG 虚拟惯量 J0.4 kg·m²对应 H≈2s阻尼系数 Dp50兼顾阻尼和稳态频率偏差无功下垂系数 Kq0.02无功每变化 1kVar电压幅值调整约 2%虚拟电阻 Rv0.5Ω抑制环流虚拟电感 Lv1mH补偿感性下垂特性1号线路阻抗0.1 j0.314Ω不对称线路阻抗测试条件2号线路阻抗0.2 j0.628Ω特意加大差异仿真解算器ode23tb变步长刚性问题适配电力电子开关模型最大步长50μs保证 PWM 开关时刻不丢4.4 并联合闸与负载投切的仿真实操双机并联不能一上来就两台同时启动那样合闸冲击电流会大得离谱。我用的顺序是先让 1 号机空载建立电压等母线电压稳定后再把 2 号机并入。2 号机并入前要做预同步检测公共母线电压的幅值、相位、频率调整 2 号机的 E 和 θ使两者差小于设定阈值我习惯控制幅值差小于 5V、相位差小于 5°然后才发合闸脉冲。在 Simulink 里可以用断路器模块加一个 Step 信号控制合闸时间。负载投切场景我做了三个第一个是 0s 启动 1 号机、0.5s 并入 2 号机、1s 时投入 8kW 2kVar 负载第二个是负载再加一倍第三个是切除一半负载看动态恢复。三个场景组合起来既能看功率均分稳态精度又能看 VSG 的动态摇摆特性。5. 仿真结果判读与常见调试问题5.1 并联合闸瞬间的母线电压建立过程第一眼要看的波形是公共母线电压。理想过程应该是1 号机空载启动后有短暂的电压建立过程E 从零按无功调节器爬升到额定值大约两三个工频周期后电压稳定在 380V 左右。随后 2 号机预同步完成合闸瞬间电流冲击不大母线电压基本连续。如果你看到合闸点母线电压出现明显凹陷或者尖峰说明预同步没做好尤其是相位差偏大。相位差大意味着两台机的电压矢量夹角大会在闭合回路里产生很大的环流。这时候别急着改控制参数先把预同步的判定阈值调严相位差压到 2° 以内再看。另一个值得注意的现象是合闸后的一小段时间内两台机之间存在低频功率波动。这是正常现象是 VSG 虚拟转子在摇摆。正常情况下波动幅度不断衰减几个周波后消失。5.2 突加负载时的功率均分波形怎么看投入 8kW 2kVar 负载后观察两台机的有功和无功波形。有功均分应该很快建立两台机最终稳定在 4kW 附近无功也应该接近 1kVar 附近。关键是看动态过程负载投入瞬间两台机的功率都会有一个尖峰随后靠 VSG 的摇摆特性调节回到平衡。有功均分如果稳态时偏差超过 2%优先检查 Pref 和 Dp 是否对称再看看功率滤波器的截止频率是否一致。无功均分如果偏差大那基本就是虚拟阻抗没起作用或者 Kq 选得不对称。我这次仿真在没加虚拟阻抗前三台机无功差 25% 左右加上 Rv0.5Ω、Lv1mH 后偏差缩到了 3% 以内效果立竿见影。还要看频率波形。负载投入瞬间频率应该有一个小幅度下降然后慢慢恢复到接近 50Hz。如果频率出现长时间等幅振荡那就是阻尼不够把 Dp 调大如果频率跌下去回不来稳态偏差过大说明 Dp 太大适当减小。5.3 最常踩的坑振荡发散、环流超标、中点漂移调试过程中我遇到过几个典型问题记录一下排查思路。第一个是双机并联后系统高频振荡发散。这个问题通常和虚拟阻抗有关。虚拟阻抗本质上是给电压参考加了一个随电流变化的修正量如果 Zv 里感性成分太大会和控制环路的电感补偿形成正反馈。解决办法是先只加 Rv不加 Lv等系统稳定后再慢慢增加 Lv每次增加都要看阻尼是否恶化。第二个是合闸瞬间环流超标峰值甚至超过额定电流。这个基本不用怀疑控制算法就是预同步没做到位。检查相位差是不是真的压下来了以及断路器合闸动作时刻和预同步完成时刻之间的延迟。仿真里断路器命令到真正合闉有导通延迟预同步判定通过后要预留这段时间再触发合闸。第三个是直流中点电位漂移导致 PWM 波形畸形。T 型三电平的中性点电流在零电平输出时流过直流电容不均等会导致两个电容电压分化。我在直流侧加了两路均压电阻同时在调制波里叠加了很小的零序分量辅助平衡。对于控制层面研究均压电阻够用如果做器件级损耗分析再考虑更精细的中点电位主动控制。另外还有一个很影响仿真体验的坑Simulink 解算器的选择。三电平 PWM 加双机并联系统刚性很强普通 ode45 可能跑得极慢甚至卡住。我用 ode23tb 配合 50μs 最大步长仿真速度能接受波形也不丢细节。变步长模式下如果最大步长取得太大PWM 的开关时刻被跳过波形上会出现奇怪的毛刺注意甄别。最后再说一句调试心得。VSG 双机并联和单机完全不是一回事单机调得好不代表并联能用。每次调参数我习惯只动一个变量并且把负载投切场景固定下来对比。J 和 Dp 的调法就是先定 Dp 消除振荡再调 J 改善惯性切忌两个一起动不然波形变了你根本不知道是哪个参数引起的。这套流程走下来双机均分、频率稳定、负载切换这几个目标都能达到剩下的就是在这个骨架上去做更深入的研究了。