1. 项目背景与研究价值做并网逆变器仿真的多半都被三相电网不平衡折磨过。电网电压一出现塌陷直流侧二次纹波、并网电流畸变、功率振荡一个接一个冒出来。这时候如果再叠加三电平逆变器的中点电位漂移问题传统的两电平单dq控制基本是撑不住的。我这次就把三电平并网逆变器T型和NPC型都做了对比放在不平衡三相电网环境里用Simulink完整搭了一套仿真模型从主电路、正负序分离、电流控制器到三电平调制全部跑通今天把这套建模思路和实操过程整理出来。这篇文章不是什么高深的新算法更多是一份可以直接上手的工程整理。适合正在做三电平并网逆变器课题、需要快速搭出仿真验证平台的同学也适合想搞明白T型和NPC型到底怎么选、不平衡工况下负序电流怎么抑制的工程师。你跟着后面的章节走至少能复现出一个能稳定运行、波形可看、参数可调的仿真平台。1.1 为什么会盯上“不平衡电网三电平”这个组合三相不平衡不是实验室里人为制造出来的极端工况。实际并网现场里单相故障切除、大功率单相负荷投切、山区线路阻抗不对称都会让PCC点电压出现负序分量。负序分量一旦进来传统同步旋转坐标系下控制出来的电流就不是三相对称的正弦波了最直观的表现是三相电流幅值不一致、带明显畸变直流母线电压上出现100Hz纹波。两电平拓扑在低压小功率场景里很成熟但到了需要更高电压等级、更大功率密度的场合两电平有天然短板开关管要承受完整母线电压输出电平只有两种带来的dv/dt和开关损耗都比较大。三电平逆变器在直流母线电压相同的前提下每只开关管承受的电压应力更低输出相电压有三个台阶等效开关频率更高并网电流谐波更小滤波电感也可以适当缩小。真正把“不平衡电网”和“三电平”放在一起研究是因为问题会叠加。电网不平衡本身就需要正负序分离控制而三电平拓扑又有自己的中点电位平衡问题和窄脉冲问题。如果只把两电平那套控制直接搬到三电平上很容易出现一边抑制负序电流、一边中点电压疯狂漂移的情况。所以这个题目值得单独建一套模型去验证而不是简单替换拓扑。1.2 三电平拓扑选型T型与NPC型怎么权衡T型和NPC型是三电平逆变器里最常见的两种实现方式。两者每相都是四个开关管加四个续流二极管都能输出正、零、负三种电平但具体的钳位思路完全不同。NPC型中性点钳位型每相额外用了两个钳位二极管把输出端钳到直流母线中点。四个开关管是串联在母线上下的正常稳态下每个管子只承受一半母线电压所以更适合母线电压比较高的场合。缺点是钳位二极管在换流过程中会有反向恢复损耗杂散电感控制不好时会带来额外电压尖峰。T型T-NPC的每相结构更像一个“半桥加双向开关”外侧两个开关管负责把输出端接到正母线或负母线中间一对反向串联的开关管负责把输出端接到中点。它的优点是少了两只钳位二极管导通路径更短电流流经的器件数量少损耗特性在中低压大电流场合更占优。缺点也很明显外侧开关管要承受完整的直流母线电压对器件耐压要求比NPC型高。我当时在Simulink里分别搭了两种模型。官方库里有现成的三电平NPC桥臂但T型不一定有对应模块我干脆用单个IGBT和二极管搭了自定义桥臂。这样做的额外好处是每个器件的开关状态、关断电压、导通电流都能单独拉出来看对理解拓扑本质很有帮助。下面这张对比表是选型时最常用到的对比项NPC型三电平T型三电平每相开关管数44每相钳位二极管数20外侧开关阻断电压约半母线电压约完整母线电压导通路径较长经过钳位二极管较短器件较少中点电压控制必须必须常用场景中高压大功率低压大电流、光伏/储能PCS选型没有绝对的对错。如果你的研究侧重点是高压大容量NPC型更合适如果关注效率、损耗和体积T型往往表现更好。后面所有控制逻辑在两种拓扑下是可以共用的差异主要在开关信号和死区配置上。2. 控制策略与核心原理2.1 不平衡电网下绕不开的两个核心难题第一个难题是负序分量带来的电流畸变。电网电压不平衡时PCC电压可以分解成正序、负序、零序三组分量。并网逆变器一般通过变压器隔离或星形接法连接零序通路基本被切断所以重点看正序和负序。你把负序分量放到正序同步旋转dq坐标系里观察会发现它是一个以两倍电网频率旋转的交流量。传统PI控制器对直流误差能够无静差调节但对100Hz交流误差无能为力最终结果就是电流环输出里混入低频谐波并网电流不再对称。第二个难题是瞬时功率的二次振荡。用瞬时功率理论展开有功功率和无功功率里会出现电网电压正序分量与电流负序分量、电压负序分量与电流正序分量交叉相乘的项频率都是100Hz。这会导致直流母线电容上出现明显的二次纹波对电解电容寿命和下一级变换器都有影响。如果中点电位平衡没做好二次纹波还会进一步放大三电平输出波形质量迅速恶化。这里要明确一个控制目标你不可能同时把“输出电流完全三相对称”和“瞬时功率完全无波动”都做到除非给电网注入额外补偿装置。单台并网逆变器在不平衡电网里通常只能二选一。我的模型选择的是负序电流抑制策略也就是让逆变器输出三相对称正弦电流负序电流参考值设为零代价是直流侧仍会保留一部分二次功率脉动。这是工程上最常用、最容易验证的方案后续想改成功率波动抑制只要重新计算负序电流参考值就行。2.2 正负序分离从αβ到dq的完整思路负序抑制的前提是能准确把正序、负序分量从采样电压里分离出来。我用的办法是DSOGI双二阶广义积分器也叫SOGI正交信号发生器。核心思想是先在αβ静止坐标系下对每个轴上的电压量做一个带通滤波同时产生两个正交信号一个和原始信号同相位另一个滞后90度。有了两套正交信号之后正序和负序的αβ分量可以用简单矩阵运算算出来。正序分量公式是vα 0.5 * (vα - qvβ) vβ 0.5 * (qvα vβ)负序分量公式是vα- 0.5 * (vα qvβ) vβ- 0.5 * (-qvα vβ)这里的vα、vβ是SOGI滤波后的输出qvα、qvβ是它们的正交滞后信号。DSOGI的好处是结构对称、实现简单而且能滤掉一部分高次谐波比纯数学陷波器更稳定。关键要调好两个参数谐振角频率固定在2pi50阻尼比取0.707左右。阻尼比太小分离出来的信号会有明显振荡太大会让动态响应变慢建议在0.6到1.0之间试。正序分量分离出来后送进同步坐标系锁相环得到电网正序角度θ。负序dq变换用的角度则是-θ。在Simulink里可以拉一个标准的PLL3ph模块也可以自己搭Park变换和PI频率环。我实际用的是自己搭的结构因为三个相电压里如果负序分量太强标准PLL模块的动态会受影响自己搭反而更容易观察中间量。2.3 电流环和功率外环参数怎么算电流内环采用正序dq和负序dq两套控制器每一套就是两个独立的PI调节器。系统干路可以简化成一个电阻加电感的被控对象忽略电网电压前馈后的传函形式为G(s) 1 / (Ls R)按典型I型系统设计把PI零点与被控对象极点对消再设定期望电流环带宽ωcPI参数可以按下面两个公式估算Kp L * ωc Ki R * ωc我用的主电路参数是直流母线电压800V并网电感3mH线路等效电阻0.1Ω。电流环带宽不敢取太高取ωc 2 * pi * 800 ≈ 5026 rad/s算出来Kp约等于15.08Ki约等于502.65。这个参数在Semulink的离散PI控制器里可以直接填。需要注意的是实际仿真里R通常很小Ki如果太小会导致稳态误差消除慢我会在此基础上微调零点位置把PI零点保持在L/R附近。电流环里还有三个工程细节必须处理加入电网电压前馈把PCC点电压叠加到PI输出上让电流环只面对电感电流的动态加入交叉解耦项也就是正序dq下的±ωLiq/ωLid可以明显降低动态过程中dq轴之间的耦合加入输出限幅和抗积分饱和防止启动瞬间PI积分冲过头导致超调。功率外环就简单多了。采用正序电压定向也就是把正序d轴跟电网正序电压矢量对齐。此时q轴正序电压为0参考电流的正序d轴分量可以这样求id_ref P_ref / (1.5 * ed) iq_ref -Q_ref / (1.5 * ed)我给的参考有功是10kWed是电网相电压峰值约311V所以id_ref大约21.4A无功参考设0。负序dq电流参考值全部设0这样控制器就会有意识地把负序电流压到接近零。3. Simulink仿真模型搭建全程实录3.1 整体模型结构与主电路参数我在Simulink里的模型整体分为五块不平衡电网电源、三电平主电路、采样测量、控制算法、PWM调制与开关信号。主电路用的是Simscape Electrical Specialized Power Systems元件控制部分用普通Simulink模块中间通过电压电流测量模块和受控电压源/开关信号连接。主电路参数我给一套可以直接跑的配置参数名数值说明直流母线电压 Vdc800V理想直流源加分裂电容分裂电容 C1/C23mF / 3mF每只电容电压目标400V并网电感 L3mH滤除开关纹波线路等效电阻 R0.1Ω模拟线路和电感内阻开关频率 fsw10kHz三电平等效频率更高控制采样周期100us与PWM周期同步仿真步长2us保证载波脉冲不丢失直流源不能简单用一个固定800V电源跨接在正负母线上那样中点电位没有通路电容电压完全无法观察。正确做法是两个额定电压400V的电容串联再把正母线、中点、负母线分别引到IGBT桥臂。直流源可以用理想DC Voltage Source电压设为总电压两个分裂电容并联在母线上也可以直接用两个独立电压源分别接在正-中和中-负之间但后一种方式电容充放电动态会不真实我更推荐前一种。T型桥臂搭建有个容易出错的地方中间那对双向开关必须用两个共用门极信号的IGBT反向串联。每个IGBT都要自带反并联续流二极管两个管子背靠背串联中间连接点作为输出端。门极给同一路PWM信号靠电流方向自然换流。NPC型桥臂则要注意钳位二极管的极性上钳位二极管阴极接输出端、阳极接母线中点下钳位二极管阳极接输出端、阴极接母线中点接反了模型会出现奇异矩阵。3.2 不平衡电网模拟怎样才算是真正的“不平衡”很多新手做不平衡仿真直接把三相电压源的幅值随机改一下结果自己都不知道到底加了多少负序分量。这样后面分析负序抑制效果时完全没有参照。我建议用Simulink里的Three-Phase Programmable Voltage Source模块在这个模块里可以按时间节点设定每一相的幅值、相位和频率。最简单的做法是设一个Base值三相平衡电压设为311V、50Hz然后在0.4秒设置事件把A相电压幅值改成0.6 * 311B相和C相保持不变。这样就制造了一个典型的单相电压跌落负序分量比例很清楚。如果想做两相跌落就把A、B两相同时改成0.6倍幅值负序分量会更大对控制器冲击也更明显。为了确认不平衡程度我建议在测量模块后面再接一个Three-Phase Sequence Analyzer它会直接输出正序、负序和零序分量的幅值。三相不平衡度的定义是负序幅值除以正序幅值用这个模块的数值可以直接算。我用A相电压跌落到60%时实测负序电压约为58V左右不平衡度大概在18%量级这个工况已经远严于国标对公共连接点2%不平衡度的要求属于很苛刻的测试条件了。3.3 DSOGI-PLL与正负序dq变换的实现正负序分离部分我单独封装了一个子系统。输入是三相电压v_a、v_b、v_c先做Clarke变换得到vα和vβ然后各通一个SOGI模块。SOGI模块内部就是一个二阶滤波结构输出两个信号滤波后的同相信号以及滞后90度的正交信号。在Simulink里可以直接用Continuous Transfer Fcn或Discrete Transfer Fcn搭这两个传函。离散化注意使用Tustin法采样时间和控制周期保持一致。我最初用连续模块跑离散求解器老是出现小步长收敛问题后来全部改成离散模块仿真稳定性好了很多。PLL角度生成我放在正序αβ分量之后用一个简单的Park变换把vα、vβ旋转到dq坐标然后q轴误差进PI调节器输出角频率积分成角度θ。PI频率环带宽我设在100Hz左右太高会把负序分量里的残余扰动放大太低在电压突变时锁不住相位。用DSOGI之后本来就有一定滤波延迟PLL带宽不用太激进。正序电流和负序电流分别采样三相电流Clarke变换后用θ和-θ做Park变换。注意正序dq里的q轴量和负序dq里的q轴量含义完全不同千万别共用同一个θ。这里如果变量名没起好后面查问题时非常痛苦我习惯用后缀_pos和_neg区分所有信号比如id_pos、id_neg。3.4 三电平调制与中点电位平衡三电平的PWM调制我采用的是载波层叠PWM也就是上下两个三角载波一个从0到1一个从-1到0三相调制波放在同一坐标系下比较。输出逻辑不是简单的“大于载波输出1”而是根据所在区域判断输出P、O、N三种状态。Simulink里可以用Relational Operator和逻辑组合实现但我更推荐直接用MATLAB Function模块写三电平PWM判断函数代码短且不容易出错。下面是我初始化脚本里和调制相关的一部分方便你直接理解参数来源Vdc 800; L 3e-3; R 0.1; fsw 10e3; % 电流环PI参数 wc 2 * pi * 800; Kp_i L * wc; % 约 15.08 Ki_i R * wc; % 约 502.65 % SOGI参数 w0 2 * pi * 50; zeta 0.707; k_sogi 2 * zeta; % 中点平衡比例系数 Kpn 0.05;中点电位平衡是三电平绕不开的问题。我的做法是在三相调制波上叠加一个零序电压校正量让冗余开关状态的占空比发生变化进而控制中点电流方向。具体校正量近似按这个思路算v_offset Kpn * (Vdc2 - Vdc1) * sign(相电流)方向非常关键我第一版把符号写反了结果中点电压在两个方向上发散震荡一直找不到原因。调试时可以先设Kpn很小比如0.01观察中点电压是往平衡方向收敛还是往更偏的方向走确认方向后再逐步增大。Kpn太大会导致调制波过调制输出电流反而畸变通常不超过0.1。死区时间也必须加。我在每个IGBT开关信号回路里串了一个On-Delay模块死区设3us。别小看这3us如果忽略死区仿真波形在过零点会看到奇怪的台阶分析THD时还会得到和实际样机严重不符的结果。4. 仿真结果分析与调参经验4.1 不平衡前后波形对比先学会看门道仿真时间我设到1秒0.4秒前是平衡工况0.4秒后A相电压跌到60%。这时先看两项核心波形第一个是并网三相电流波形。没有负序抑制时A相电流幅值会明显变大或变小BC两相也不对称波形里能看出100Hz包络。加上负序抑制后三相电流幅值基本一致波形接近正弦说明负序电流被压下去了。这个结果是方案有效的最直接证据。第二个是直流分裂电容电压Udc1和Udc2。理想情况下两只电容电压稳定在400V附近实际不平衡工况下会有波动。如果中点电位平衡环节有效两个电压之间的差值平均值应该基本为零波动幅值控制在几伏到十几伏以内。如果看到两只电容电压像跷跷板一样来回偏说明中点平衡方向反了或者Kpn参数太小。对电流质量做定量分析用powergui里的FFT Analysis工具选A相并网电流基波设为50Hz采样10个周期最大频率显示到2kHz。负序抑制开启后我测到的THD在3%到5%之间主要谐波集中在10kHz开关频率附近关闭负序抑制后THD会直接翻倍甚至更高。记住FFT窗口要避开0.4s突变点取0.6s到0.8s比较稳定。4.2 关键参数对结果的影响调参顺序与推荐区间调参一定要按顺序来不然容易产生“按下葫芦浮起瓢”的错觉。我的经验是先调主电路和PWM再调电流环最后调不平衡工况相关的SOGI、PLL和中点平衡。电流环带宽是第一个要盯的参数。带宽太低电流跟踪慢负序抑制效果会打折扣带宽太高采样延迟和数字控制延迟会引入明显超调。用直流母线800V、电感3mH这个配置带宽800Hz到1kHz比较合适对应Kp大约15到19Ki大约500到630。如果你的电感更大Kp要跟着增大否则动态响应会变慢。SOGI阻尼比建议固定0.707除非你发现分离出来的负序分量波动很大。阻尼比调小会让滤波更尖锐但也更容易在暂态时振荡。中点平衡比例系数Kpn从0.02起步每次增加0.01观察中点电压波动和THD变化。死区时间按开关频率反比放大10kHz下取2到4us都合理死区太大反而会引入3、5次谐波。还有个容易被忽略的参数是PLL带宽。它决定了对相位跳变的响应速度。不平衡电压跌落的瞬间相位会有跳变PLL带宽太低会看到电流相位慢慢摆回来带宽太高则在稳态时出现相位抖动。100Hz左右的带宽是我在这个模型里试出来最稳的中值。4.3 我在现场踩过的三个坑第一个坑就是前面说的中点平衡方向搞反。当时波形里电容电压始终处于发散状态我一度以为是SOGI分离出来的正负极序有问题查了半天发现只是Kpn前面的符号错了。这个教训说明调新参数前一定要先做方向性验证哪怕算法是抄来的也要亲手确认反馈极性。第二个坑是PWM载波相位和死区导致开关脉冲丢失。用Simulink的PWM Generator模块时如果不小心把载波频率设成和控制频率不一样的数值会出现电流波形低频包络看着像不稳定实际上是调制比错位。我后来自己用MATLAB Function生成开关逻辑把载波周期固定为控制周期的两倍问题当场消失。第三个坑是仿真步长太大导致结果虚高。最初我把仿真步长设成20usPWM载波10kHz下每个载波周期只有5个采样点开关时刻量化误差严重FFT分析显示谐波分布乱七八糟。改成2us后波形才稳定。Simulink里连续电力电子模型用离散求解器时仿真步长至少要比载波周期小20倍这是个很实用的经验数值。5. 常见问题与排查技巧实录5.1 仿真发散问题先排除这些根因很多Simulink电力电子模型一运行就报“singular matrix”或者“simulation stops due to numerical error”多半不是控制算法问题而是主电路结构有问题。最容易出现的包括现象常见原因排查方法仿真一开始就发散分裂电容接法不对或桥臂存在悬空状态检查直流源、中点、桥臂连接关系T型桥臂双向开关失效两个IGBT门极信号不同步确认共用同一门极信号NPC桥臂报奇异矩阵钳位二极管极性接反核对上下钳位二极管阴极/阳极方向控制闭环后发散反馈极性错误或电流环无穷大断开闭环先开环测试PWM状态我的习惯是先把PWM占空比固定成常数比如0.5不带闭环看主电路是否正常输出三电平阶梯波。主电路正常后再接入电流环最后才接PLL和正负序分离。这个顺序能帮你把问题范围快速缩小。5.2 负序抑制效果不理想的排查步骤如果你按代码搭好了模型却发现三相电流还是不对称先别急着调PI参数按下面顺序检查第一步看DSOGI输出。把正序αβ分量和负序αβ分量拉到Scope里检查正序分量是否幅值稳定、负序分量是否存在。如果负序分量明显跟着电网电压变化说明分离本身是正常的。第二步看负序dq电流的给定值和反馈值。负序d、q参考都是0如果反馈id_neg和iq_neg不为零说明电流环没有把误差压下去。这时候检查负序PI控制器是否饱和、抗积分饱和是否卡住了输出。第三步看PWM的调制波是否同时叠加了正序和负序两路信号。有些模型正负序电流环输出的调制度在反变换到αβ时符号搞错负序补偿反而变成了正反馈。把叠加后的三相调制波和原始正弦调制波对比就能看出异常。如果这些都没问题那大概率是锁相环角度有偏差。负序dq变换使用-θ而-θ和θ不是简单取反就行反Park变换里的角度也要一起动。在我自己的模型里负序控制环节的所有Park反变换角度全部用-θ这个细节容易漏。5.3 模型跑太慢怎么加速Simscape电力电子模型天然比普通Simulink模型慢这是正常的。但如果你为了等项目跑完一等就是一下午有几个加速技巧很有效第一用离散电力电子模型而不是连续模型。Specialized Power Systems里把powergui切到Discrete并选择合适步长。离散模型和连续模型的精度差别在这个频段内并不大但仿真速度差距非常明显。第二把观察方便的Scope改成To Workspace模块或者只保存需要分析的数据。一个高采样率的Scope即使没打开也会拖慢仿真速度。我通常只在关键观察点加Scope多余示波器全部删掉。第三如果只是为了验证稳态控制策略可以用平均模型或简化开关模型。把每个桥臂等效成比例放大器开关频率的细节没有了但负序抑制、中点平衡这类控制效果仍然能看。等你确认控制逻辑无误后再切回详细开关模型做最终波形分析。最后再分享一个我个人比较受用的点把整个模型分成控制和主电路两个子系统中间用Goto/From或者总线对象连接这样调一个模块时不需要动其他模块排查问题定位也快。我后来还在这个模型基础上试过加入LCL滤波器和谐振阻尼控制结构变化不大但并网电流质量又上了一个台阶。如果你也在做三电平并网方向这个模型完全可以继续往LCL、弱电网适应性、多机并联这些方向延伸。先跑通再优化永远是仿真研究最靠谱的路径。