光伏逆变器并网仿真复现：SMIB模型下PLL与dq电流控制全解析

去年年底我们团队做仿真工具国产化迁移第一批要复现的就是一套光伏逆变器并网模型。原来在 Simulink 里跑得挺好换到 MWORKS Sysplorer 之后才发现这根本不是“重新画一遍图”的事。PLL 锁相环的初始相位怎么给、dq 电流控制的前馈解耦怎么接、SMIB单机无穷大母线下电网电压怎么等效每一步没想透波形上立刻就能看出来。这篇把我搭这套仿真模型的完整过程记录下来从系统架构到具体参数从 PLL 整定公式到调试时踩过的坑。不需要你有很深的电力电子基础但至少要知道三相电压、PI 控制器大概是什么。读完一遍之后你至少能照着参数表把模型复现出来更重要的是知道为什么要这样搭。1. 为什么重做这套光伏并网仿真项目背景与 SMIB 的取舍1.1 先交代项目时间点2024 年初我们接到一个任务把团队沉淀了多年的仿真资产从国外商业软件往国产科学计算与建模平台迁移第一批候选就是 MWORKS Sysplorer。理由不复杂许可证成本逐年走高供应链受限公司层面要求核心研发流程有可替代方案。技术选型对比做了两三个月最后定下 Sysplorer看中的是它既支持 Modelica 多领域物理建模又有类 Simulink 的框图建模环境光伏逆变器这种“电路控制”混合系统正好能覆盖。真正动起手来才发现网上能找到的参考资料非常少尤其是“光伏逆变器并网 SMIB 仿真 PLL dq 电流控制”这种具体组合基本全靠啃软件自带的库模型和帮助文档。这篇文章算是我把过程整理出来希望能让后面的人少走弯路。1.2 SMIB 到底是什么意思SMIB 全称 Single Machine Infinite Bus单机无穷大母线系统。在电力系统分析里它表示你研究的对象一台发电机、一台逆变器本质上就是一个“电源机器”通过线路阻抗接入一个理想电网。理想电网的特点是电压幅值恒定、频率恒定、容量无限大不受研究对象任何影响。为什么光伏并网研究喜欢用 SMIB 当作初始场景因为它的边界条件最干净。当你要验证逆变器控制策略是不是正确时不需要考虑电网阻抗波动、电压跌落、谐波背景这些复杂扰动先把“一个逆变器接到一个理想电网上、输出功率可控”这件事做对、做稳再去叠加电网故障和弱网工况才有意义。平台迁移也一样先用 SMIB 这种基准算例验证模型移植前后结果一致后面扩展才有底气。1.3 这套模型适合谁复现如果你正在做下面几件事这篇文章基本就是给你写的一是刚接触光伏逆变器并网仿真想找一个结构清晰、参数完整、能跑的入门模型二是已经在 Simulink 里搭过类似模型正在评估迁移到 Sysplorer 的工作量三是已经搭好了模型但 PLL 经常锁不住、dq 电流纹波大、启动就过流需要一份排查思路参考。我尽量把每一步的原理和参数来源都写清楚这样即使版本不同、参数不同你也能举一反三。2. 系统级架构主电路长什么样控制链路怎么走2.1 主电路与光伏侧等效整体拓扑是光伏并网逆变器最常见的两级式结构前级是光伏阵列加 Boost 升压电路负责把光伏板输出电压升到直流母线电压同时跑 MPPT 最大功率点跟踪后级是三相桥式逆变器负责把直流母线电压逆变成三相交流经过滤波电感并入电网。在 SMIB 场景里前级光伏阵列和 Boost 我直接用恒压直流源替代了不是偷工减料而是刻意简化。原因有两个第一本次仿真的核心研究对象是 PLL 锁相和 dq 电流内环MPPT 动态响应比电流环慢一个数量级以上短时间内对并网控制的影响可以忽略第二把直流侧固定成恒定电压等于先验证“在理想直流条件下交流侧控制正不正确”如果这一步都做不好叠加了 MPPT 波动只会更难排查。等整套系统验证通过再把 Boost 和光伏模型接回来循序渐进。滤波环节用的是 L 型电感不是 LCL。LCL 滤波在高频谐波抑制上确实更好但它会引入谐振峰阻尼设计不好反而容易振荡。第一版迁移模型讲究稳妥L 单纯、直观、模型简单对验证控制策略来说足够。LCL 的谐振阻尼问题适合建模稳定之后再专门研究。2.2 双闭环控制链路控制结构是经典的电压外环加电流内环双闭环。先看信号流采样三相电网电压和三相并网电流电网电压送进 PLL锁出电网电压矢量的实时相位角 θ用这个 θ 对三相电流做 Clarke-Park 变换得到 dq 轴电流分量 id 和 iq直流母线电压外环 PI 根据电压误差给出 id 参考值iq 参考值通常设 0单位功率因数运行两个电流 PI 分别控制 id、iq 跟随参考值输出经过前馈解耦和电网电压前馈得到 dq 轴电压指令再经过反 Park 变换得到三相调制波送进 SVPWM 产生驱动脉冲。一句话概括PLL 解决“坐标系靠不靠得准”的问题dq 电流环解决“电流跟不跟得上”的问题直流电压外环解决“能量平不平衡”的问题。三个环节互相依赖调试顺序也基本按照这个依赖关系来。2.3 时间尺度差异为什么重要双闭环能成立的前提是内环比外环快得多。电流环面对的是电感这种小时间常数对象典型带宽几百到上千 rad/s响应时间在毫秒级直流母线电压环面对的是母线电容时间常数大得多响应时间在几十到上百毫秒。带宽设计上两者至少拉开 5 到 10 倍外环看到的内环近似成一个增益为 1 的快速随动环节控制设计才能解耦。这不是经验主义而是经典的“时间尺度分离”思路搞清楚了才明白为什么电流环 PI 参数看起来“很大”、电压环 PI 参数看起来“很小”。3. PLL 锁相环的建模与参数整定让坐标系真正对齐电网3.1 SRF-PLL 的信号流结构并网逆变器里最常用的锁相方案叫 SRF-PLL同步参考坐标系锁相环。名字听着复杂本质就是一个闭环角度跟踪器先把三相电网电压 Ua、Ub、Uc 做 Clarke 变换得到 αβ 两相电压再用当前估计的角度 θ_est 做 Park 变换得到 dq 轴电压 ed 和 eq。这里有个关键点当 θ_est 恰好等于真实电网电压矢量角度时eq 会等于 0ed 等于电网电压幅值。所以把 eq 当作误差信号送进 PI 调节器PI 输出是角频率修正量叠加到额定角频率 2π×50 上再经过积分器得到 θ_est这个 θ_est 又反过来驱动 Park 变换形成闭环。这个环路在物理上可以类比成一个自动追踪目标的雷达天线。目标位置是真值角度天线指向是估计角度eq 就是两者的偏差信号PI 不断调整转动速度直到偏差归零。一旦 eq 稳定在 0 附近就说明估计角度和电网真实相位对齐了。3.2 PLL 参数计算过程PLL 的 PI 参数不能瞎拍它有个清晰的小信号模型。设电网电压幅值为 Vm相位误差为 θ_e在小偏差条件下有 eq ≈ Vm × θ_e。从 θ_e 到最终角度输出的开环传递函数是 Vm × (kp ki/s) × 1/s即一个含一个零点、两个积分环节的系统。闭环特征方程为s² Vm × kp × s Vm × ki 0对照标准二阶系统 s² 2ζωn × s ωn² 0就能解出kp 2ζωn / Vmki ωn² / Vm以我的模型参数为例电网线电压有效值 380V相电压幅值 Vm 311V。锁相环闭环带宽取 50Hz也就是 ωn 2π×50 ≈ 314 rad/s阻尼比取典型值 0.707。代入计算kp 2×0.707×314 / 311 ≈ 1.43ki 314² / 311 ≈ 317实际取值我用的是 kp 1.5、ki 300非常接近。注意 PLL 带宽不能取太高锁相环本质是从畸变、不平衡的电网电压里提取基波正序相位带宽越高电网谐波和负序分量就越容易混进估计相位里带宽太低又会锁得慢、动态响应差。50Hz 对工频并网是一个折中偏稳的选择。3.3 为什么会出现“PLL 没有锁定”这是并网仿真里最经典的问题之一。PLL 没锁定多数不是参数错了而是初始状态没给对。仿真刚开始时如果 θ_est 的初始值和电网电压 A 相初始相位差得远PLL 要经历一个较大的瞬态追逐过程。在追逐过程中可能表现为频率振荡、dq 电压波形剧烈抖动严重时直接导致系统发散。处理办法有两个一是把 PLL 积分器的初始角度设置成电网 A 相电压的初始相位角让模型一开机就在“锁定附近”这是最推荐的做法二是等 PLL 稳定后再投入电流环和 PWM可以在控制逻辑上加一个使能信号。另外要检查 PLL 的量纲和基准。不同坐标变换采用的幅值不变还是功率不变约定会直接影响前向通道增益和 PI 参数计算结果。例如功率不变变换下 d 轴电压不再是 Vm 而是 1.5 倍或 √3/2 倍的关系没统一基准PI 参数按公式算出来自然锁不正。4. dq 电流内环解耦、前馈和 PI 参数的计算逻辑4.1 从三相方程到 dq 方程三相电感滤波器的电压方程写出来是三个互相耦合的一阶微分方程。这三个量都在以工频振荡直接设计 PI 控制器很别扭。把它变换到与电网电压同步旋转的 dq 坐标系后原来 50Hz 的正弦电流变成了直流分量 id 和 iqPI 控制器对直流参考值可以做到无静差跟踪这就是 dq 控制的核心价值。推导后的 dq 轴电压方程如下L × did/dt Vd - ed - R × id ωL × iq L × diq/dt Vq - eq - R × iq - ωL × id其中 Vd、Vq 是逆变器输出端电压在 dq 轴的分量ed、eq 是电网电压在 dq 轴的分量。注意到一个问题id 的方程里出现了 iqiq 的方程里出现了 id这两项 ωL × iq 和 ωL × id 就是交叉耦合项。如果不处理它d 轴电流的变化会扰动 q 轴q 轴也会反过来扰动 d 轴两个 PI 环会“打架”。4.2 前馈解耦的具体形式解决交叉耦合的办法是前馈补偿在 d 轴电压指令里加上减去 ωL × iq在 q 轴电压指令里加上 ωL × id正好把方程里的耦合项抵消掉。同时把电网电压 ed、eq 也前馈进去让电流环不用费力去对抗电网电压。最终控制律写成Vd_ref PId ed - ωL × iq Vq_ref PIq eq ωL × id这样做完之后id 和 iq 两个通道近似解耦各自独立成一个“RL 串联负载 PI 控制器”的单输入单输出系统。这也是为什么很多教材上直接说“dq 电流控制本质上是把三相耦合系统拆成两个独立的直流系统来控”。不解耦的后果不一定发散但动态响应明显变差电流调节过程会出现很长的振荡拖尾。4.3 PI 参数整定的计算推导电流环的 PI 参数完全可以算出来不必试凑。被控对象是电感加电阻的串联传递函数为 1/(LsR)。PI 控制器传递函数为 kp ki/s。开环传函是两者相乘。一个非常实用的做法是零极点对消令 kp L×ωcki R×ωc这样 PI 的零点正好抵消被控对象的极点闭环传函退化成纯一阶惯性环节带宽就是 ωc。代入模型参数L 3mHR 0.03Ω电流环闭环带宽取 1000 rad/s约 159Hz远小于 10kHz 开关频率合理。于是kp 3mH × 1000 3 ki 0.03Ω × 1000 30注意一个实操细节这里 R 是电感等效串联电阻和桥臂导通压降等效电阻的总和实际中 R 可能比 0.03Ω 更小或更大直接按公式算出的 ki 往往偏小导致稳态时积分环节补偿电网电压波动和前馈误差的能力不足。我见过不少人在这个基础上直接把 ki 放大十倍比如取 300也能稳定但代价是引入更大超调。我的建议是先按公式算再在 ±50% 范围内试凑每次只动一个参数配合波形观察。4.4 PWM 与直流母线电压约束电流环输出的是 dq 轴电压指令要经过反 Park 变换变成三相调制波再进 SVPWM 生成驱动信号。SVPWM 相比 SPWM 的优势是直流电压利用率高约 15%这一点在并网逆变器里很重要。以 700V 直流母线为例SVPWM 能输出的最大相电压幅值是 Vdc/√3 ≈ 404V而额定工况下我们需要的相电压幅值大约 311V有一定裕量来补偿电感压降和动态调节。但裕量不等于无限PI 输出超调大时照样会顶到限幅。调制波一旦进入饱和区输出电流就会畸变id、iq 出现明显的非线性纹波。所以电流 PI 的输出限幅值要设置成和 SVPWM 线性区范围匹配限幅范围留出 10% 到 20% 的余量防止动态过程中过量饱和。积分环节还要加抗饱和处理否则每次限幅期间积分器还在累积退饱和之后会出现很长时间的过冲。5. 在 Sysplorer 中落地模块拆分、参数表和仿真配置5.1 物理层与信号层拆分思路在 Sysplorer 里建模最重要的一件事是把模型按“物理域”和“信号域”分开组织。主电路属于物理域包括直流源、三相桥臂、滤波电感、线路阻抗和电网电压源建议用 Modelica 电气模型库搭这样能拿到物理单位一致的电压电流信号。控制回路属于信号域包括 PLL、坐标变换、PI 调节器、SVPWM 算法用框图算子搭信号流和 Simulink 风格一致。两域之间通过测量模块和受控源接口相连。这种拆分的好处是调试方便。物理域出问题波形上看电压电流大概率能看出是过流还是振荡控制域出问题直接看 dq 分量和调制波的波形定位快。混在一个域里的话信号量纲混乱排查成本直线上升。Sysplorer 具体模块名称在不同版本可能略有差异我自己用的版本里主电路相关的三相理想电压源、IGBT 桥臂、电感和测量元件都在电气库目录下找得到控制部分用求和、增益、积分、PID 等基础算子组合。如果你装的版本模块叫法不同对照功能找即可原理不变。5.2 关键参数表参数表整理如下所有数值都是我实际用过的可以直接当初始值使用参数数值含义电网线电压有效值380VSMIB 理想电网电网频率50Hz工频直流母线电压700V模拟 Boost 输出滤波电感 L3mHL 型滤波器等效电阻 R0.03Ω电感串联电阻等线路阻抗0.1Ω 0.1mH逆变器到无穷大母线开关频率10kHzSVPWM 载波频率电流环带宽1000 rad/s内环动态指标电流环 kp / ki3 / 30按公式计算所得PLL 带宽50Hz阻尼比 0.707PLL kp / ki1.5 / 300按公式计算取整直流电压环 kp / ki0.15 / 15外环带宽约 15Hz直流电压环参数我没有展开推导它面对的电容器对象模型简单带宽设置比电流环低一个量级先知道这个原则即可。5.3 搭建步骤与仿真求解器设置我的搭建顺序是先搭主电路不接控制用固定三相调制波验证 PWM 和功率电路能正常工作再单独搭 PLL只接电压测量观察 θ 是否按 50Hz 速度旋转、eq 是否趋近 0PLL 确认稳定之后才接电流环先给固定 id_ref 做电流跟踪测试最后再加直流电压外环闭合。每步都出波形、确认正常再进下一步。仿真求解器设置上PWM 是 10kHz 的开关信号固定步长必须能分辨开关边沿我用的步长是 1e-5 秒也就是 10kHz 一个周期正好 10 个计算步。再大容易在开关切换处出数值问题。仿真时间不用长0.5 秒足够看到启动暂态加稳态过程。如果模型编译报代数环错误优先检查 PLL 的相位反馈路径、电流环的前馈路径必要时在反馈采样环节加一个单位延迟断开代数环这在离散化控制里本来就是合理的做法。6. 调试实录启动冲击、PLL 未锁定和纹波问题的完整排查6.1 第一次启动就过流PLL 初始相位没对齐第一次闭环跑起来电流直接冲到了额定值的五倍以上眼看着波形从零飙到天上然后系统发散。第一时间怀疑是 PI 参数错了但算过好几遍没有问题。后来把电网电压 A 相和锁相环估计角度画在同一张图上真相立刻暴露PLL 积分器初始角度设的是 0但电网 A 相电压初相是 30 度开机瞬间坐标系就偏着电流环按错误的坐标系给了错误的电压指令自然产生巨大冲击。解决办法是读取电网电压初始相位写到 PLL 积分器的初始条件里。不同软件实现方式不同但思路相同让估计角度和真实角度从 t0 那一刻就基本对齐。改完之后启动过流峰值降到了额定值的 1.2 倍以内完全可接受。6.2 PLL 频率振荡和“未锁定”现象还有一次现象是 PLL 输出频率在 48Hz 到 52Hz 之间来回晃整个仿真像“打了摆子”。查了半天问题出在我把 PLL 的 PI 输出直接接到了积分器忘了加额定角频率前馈。正确的 PLL 结构里PI 应该只输出频率修正量2π×50 的额定角频率要独立加上去。少了这个前馈PI 得靠积分慢慢把频率从 0 推到 50Hz过程中还容易超调。这个坑对应很多人在仿真里遇到的“PLL 没有锁定”现象。检查要点就两个一是初始相位有没有对齐二是额定角频率前馈有没有加。把这两个地方修好锁相环基本不会出大问题。如果锁上之后还有轻微波动那大概率是电网电压幅值基准没统一导致整定的阻尼比和计算值不一致。6.3 电流有稳态误差和 2 倍频纹波电流环正常跟踪之后我又发现两个小毛病。一个是 id 稳态值比参考值低了大概 5%而且动态响应慢。原因是电网电压前馈项 ed、eq 没加电流环靠积分器硬扛电网电压这个“大扰动”。补上 ed、eq 前馈之后稳态误差基本消失动态响应也明显利索。另一个是 dq 电流波形上叠加了一个 100Hz 的纹波这种一般不是电流环的问题而是电网侧或者锁相环角度里带了二次谐波。PLL 输出角度只差哪怕 1 度的波动变换到 dq 域就会在电流上看到同频纹波。排查顺序很重要先看 ed、eq 是否干净再看 PLL 角度有没有二次谐波最后才是动电流环的参数。不要一上来就盲目调大 ki那样可能暂时压住纹波但会牺牲动态性能和稳定性。6.4 判断仿真结果是否正确的检查清单很多新手看仿真波形趋势差不多就认为对了其实有几个硬指标可以快速自检。第一PLL 锁相后的 ed 应该约等于电网电压幅值 311Veq 应该约等于 0θ 以 50Hz 速率线性增长频率稳定在 50Hz 附近。第二id 稳态应该精确跟踪参考值iq 稳态接近 0单位功率因数下计算的无功功率接近 0。第三三相电流波形正弦、对称THD 不出奇地高动态过程峰值电流不超过额定值的 1.5 到 2 倍。第四调制波幅值在 SVPWM 线性区范围内没有长期饱和。这四组波形全部合格这套 SMIB 并网模型才算真正跑通。最后分享一个我自己养成的小习惯把 PLL 参数、电流环参数、被控对象的电气参数全部列成一张表每次只改一个变量记录一组波形。调参数最怕同时动好几个旋钮出了问题根本不知道是哪个改动引起的。锁相环锁不上别急着怀疑 PI 公式先检查初始相位和基准统一性电流追踪不稳先确认解耦和电网电压前馈是否接对再动带宽。这套 SMIB 仿真做完你对并网逆变器控制的整体理解会比只看书深刻得多毕竟工程里的逻辑只有亲手点火试过才算真正长在自己身上。