这几年做新能源并网仿真VSG虚拟同步发电机基本是绕不开的一个方向。风电、光伏这类电源和传统同步机组最大的区别就是电力电子接口几乎没有旋转惯量和阻尼电网频率波动一来它们既不会“顶一下”也不会主动抑制振荡。所以这两年大家都在研究怎么让变流器“表现得像一台同步机”而风储VSG就是很典型的一条技术路线——风机负责发电储能负责补功率VSG控制负责给整个系统提供虚拟惯量和一次调频能力。这篇文章我会完整拆一遍我自己搭的风储VSG并网Simulink仿真模型从系统架构、VSG控制原理、参数整定到Simulink里每个关键模块怎么搭、参数怎么设、遇到发散和振荡怎么排查全部写清楚。适合正在做新能源并网仿真、微电网控制或者毕业论文需要仿真支撑的同学参考内容偏工程实践不是教材理论复述。1. 风储VSG系统方案让风储机组“表现得像一台同步机”1.1 为什么要用VSG惯量问题的本质传统电力系统里同步发电机有实实在在的转子系统频率变化时转子的旋转动能会自动参与调节——频率下降转子减速释放动能频率上升转子加速吸收动能。这就是电力系统频率稳定的第一道防线。但风电经过背靠背变流器并网电机和电网之间完全是“电气隔离”的电网侧变流器本质上就是一个受控电流源或电压源它对系统频率的变化毫无物理响应。新能源渗透率越高系统惯量越低频率跌落速度越快这已经是业内共识的问题。储能加入之后情况就变了。储能系统能快速吞吐有功功率但如果没有合适的控制策略它依然只是一个“被动执行指令”的电源。VSG控制的本质就是把同步发电机的转子运动方程“写进”逆变器的控制算法里让逆变器外特性具备惯量、阻尼和一次调频能力。为什么要“虚拟”因为我们没有真实的转子但可以通过控制算法在功率响应上等效出一个转子。为什么要“同步发电机”因为电网的稳定运行框架是基于同步机的让变流器贴近同步机特性是最平滑、最安全的并网方式。1.2 整体拓扑风机、储能、VSG变流器怎么连我搭的模型采用共直流母线结构风机单元和储能单元都汇接到同一条直流母线上母线上再挂一台VSG逆变器经LCL滤波和并网开关接入三相电网。这种结构的优点是控制层次清晰——风机管发电、储能管直流母线电压、VSG管并网交流侧三个环节各自闭环调试时可以单独验证每个部分。风机单元风速输入到风力机模型计算捕获的机械功率带动永磁同步发电机发电机侧整流器负责最大功率跟踪MPPT。风机的直流侧输出功率就按当前风速下的最大功率给出。储能单元电池组通过双向DC/DC变换器接到直流母线控制目标是稳住直流母线电压。风机功率波动时储能自动吸收或释放差额功率。VSG逆变器三相两电平桥控制核心是VSG算法生成虚拟电动势和虚拟相位再通过电压电流双闭环和PWM输出实现并网功率控制和惯量支撑。为什么不用两台变流器分别控制风机和网侧核心是因为VSG本身需要“功率可调、能量可缓冲”风机本身很难做快速的功率双向调节而储能恰好补上这个短板。所以风储共直流母线再统一用一个VSG接口并网是工程上比较合理、仿真也容易收敛的方案。2. VSG核心原理转子运动方程与关键参数整定2.1 虚拟同步机的“运动方程”VSG控制的核心在功率环。同步发电机的转子运动方程可以写成J * dω/dt Tm − Te − Dp(ω − ω0)其中J是虚拟转动惯量ω是虚拟角速度Tm是虚拟机械转矩Te是电磁转矩Dp是阻尼系数ω0是额定角速度。把这个方程搬到功率域两边乘以ω可以得到更常用的形式J * ω0 * dω/dt P_ref − P_e − Dp * ω0 * (ω − ω0)这个式子的工程含义非常直接当参考功率P_ref和实际输出功率P_e出现偏差时虚拟角速度会变化而角速度的变化又通过积分生成虚拟相位θ进而改变输出电压相位最终调节输出功率。这个闭环就是VSG的“灵魂”——它天然具备一次调频能力频率跌落时输出功率自动增加频率升高时自动减少。无功—电压环则相对简单通常用下垂关系E E_ref n_q * (Q_ref − Q_e)输出的虚拟电动势幅值E结合虚拟相位θ就生成了三相参考电压u_a E * sin(θ) u_b E * sin(θ − 2π/3) u_c E * sin(θ 2π/3)这个三相参考电压会送给电压电流双闭环最终驱动PWM。2.2 惯量J、阻尼D、下垂系数的取值与计算参数整定是VSG仿真里最容易踩坑的地方我最初搭建时也经历过参数不当导致的持续振荡。先说结论J决定虚拟惯量大小J越大频率变化越平缓惯量支撑效果越明显但J过大会导致功率振荡、响应变慢。D起阻尼作用D越大振荡衰减越快但D过大会让功频下垂特性变硬相当于一次调频能力下降。参数整定的工程思路一般有两种从惯性时间常数出发或者从功频下垂关系出发。我自己的做法是先从惯性时间常数出发。以我模型里的一组典型参数为例参数数值说明系统额定容量 S100 kVA基准功率额定频率 f050 Hz电网额定频率额定角速度 ω0314.16 rad/s2πf0惯性时间常数 H2 s模拟同步机惯量水平虚拟转动惯量 J2HS/ω0² ≈ 4.05 kg·m²由于J 2HS/ω0²阻尼系数 D8.1 N·m·s/rad取J/D ≈ 0.5 s这里惯性时间常数H参考的是同步发电机经典取值范围。同步机组的H通常在2~10秒之间新能源VSG可以适当取小我一般取2~4秒。算J的公式就是J 2HS/ω0²所有单位都要换成标准单位S取Wω0取rad/s否则算出来量级会差很远。阻尼系数D的取法有两种思路一种是按J/D的时间常数来取这个时间常数一般在0.5~2秒我常用0.5秒所以D J/0.5 ≈ 8.1另一种是从功频下垂系数反推比如允许频率偏差0.2Hz时输出100%功率那么稳态下垂系数就是Dp * ω0 * 0.2*2π ≈ S推导出Dp ≈ 79.6这个值和时间常数法差得较多但两者在不同控制结构下对应。实际调试中我更推荐先把D取大一点让系统稳定再逐步减小D观察阻尼效果。J则相反先把J取小让系统响应快再增大J观察惯量支撑效果整个过程就像调两个互相制约的旋钮最终找到一组“响应快且不振荡”的参数。无功下垂系数n_q的整定更直观允许电压偏差5%时输出100%额定无功那么n_q 0.05 * E_n / Q_n。比如额定电压幅值311V、额定无功30kVarn_q就在0.5 V/Var这个量级附近具体跟容量有关需要用表格中的基准值换算。2.3 有功频率环和无功电压环的控制结构有功频率环的功能是把功率偏差转化为虚拟相位结构上是积分链P_ref和P_e相减经过J*ω0和时间常数处理生成dω/dt积分得到ω再积分得到θ。注意这里不直接对ω做限幅要对dω/dt和θ做处理防止启动冲击。无功电压环生成的是虚拟电动势幅值E它在下垂的基础上通常还加一个PI调节器来消除稳态误差。我的模型里用的是下垂加积分的形式外环输出E_ref限幅在0.9~1.1倍额定电压范围内。这个环节还要重点提一下预同步问题。直接并网的瞬间如果VSG输出电压和电网电压存在幅值差或相位差会产生很大的冲击电流严重时直接导致仿真发散。我在模型中加了预同步环节并网前用虚拟相位追踪电网相位幅值追踪电网幅值等两者误差足够小再合闸。这是仿真从“能转”到“能并网”的关键一步。3. Simulink仿真模型搭建从顶层到底层逐个模块拆解3.1 顶层结构与求解器设置Simulink模型顶层我按功能拆成五个模块区域风速与风力机模块、机侧整流与MPPT模块、储能双向DC/DC模块、VSG控制与逆变器模块、电网与负荷模块。每个模块用子系统封装输入输出接口留清楚。这样有两方面好处一是调试时能单独跑每个子系统方便定位问题二是参数集中投放在模型工作区不用每次改参数都进子系统里翻。求解器设置建议用变步长我用的求解器是ode23tb相对误差设为1e-4最大步长设为1e-4秒。不要用默认的变量步长设置高压开关动作、PWM脉冲这种快速暂态如果步长控制不好很容易出现莫名其妙的数值震荡。如果模型比较大可以把PWM的载波频率降低到5kHz左右采样时间继续控制在1e-5秒能明显加快仿真速度。3.2 风机模型与储能直流母线单元风机模型在VSG仿真里不需要特别精细重点在功率特性。风力机捕获的机械功率公式是P_m 0.5 * ρ * π * R² * C_p * v³其中ρ取1.225 kg/m³叶轮半径R取10m风速v是输入Cp是风能利用系数。Cp我直接用查表模块实现在最优叶尖速比下Cp取0.48这样不同风速下按MPPT曲线给出对应的机械功率就能模拟出发电功率随风速变化的特性。如果你不需要研究风机内部的电磁过程用这个简化模型完全够用而且能大幅降低仿真收敛难度。机侧整流器我用的三相不可控整流加Boost电路或者直接用一个受控电流源来等效风机输出功率两种做法我都试过。受控电流源的方式最稳定适合重点关注VSG行为的场景带电力电子开关的完整模型更真实但仿真时间长、收敛难度大。我的建议是先搭受控电流源版本把控制逻辑跑通再逐步替换成完整模型不要一上来就所有模块全开。储能单元的核心是双向DC/DC加上电池模型。电池我用的是Simscape Electrical里的Battery模块额定电压800V容量100Ah。双向DC/DC的控制目标锁定直流母线电压外环是电压PI内环是电感电流PI。母线电压给定值我取1100V这个值要高于逆变器交流侧电压峰值对应的直流侧最小电压否则逆变器无法正常工作。储能控制的最大难点在功率限幅。电池不可能无限吞吐功率所以我在电流内环之前加了一个基于SOC的功率限幅模块当SOC低于20%时限制放电功率高于90%时限制充电功率。很多初学者仿真发散其实是电池在极端工况下被强制输出不可能的功率造成的。加了这个限幅模型在各种工况下都能保持稳定。3.3 VSG控制算法在Simulink中的实现VSG控制算法是整个模型的核心我用Matlab Function模块实现核心方程再配合Simulink标准积分器这样的组合比纯模块连线更清晰排错也方便。Matlab Function核心代码如下function [theta, E, omega] vsg_core(P_ref, P_e, Q_ref, Q_e, omega, V_g) % 有功-频率环转子运动方程 domega_dt (P_ref - P_e)/(J*omega0) - Dp*(omega - omega0)/J; % 这里用了简化形式实际模型中通过积分器完成omega的积分 % 无功-电压环下垂积分 E E_ref n_q*(Q_ref - Q_e) k_iq*integral_q; % 虚拟相位 theta integral_omega; end实际搭建时为避免代数环我倾向于用积分器模块而非Matlab Function内部的积分操作。具体接法如下功率计算模块采集逆变器输出三相电压uabc和三相电流iabc通过Park变换转到dq坐标系用公式P 1.5*(v_di_d v_qi_q)、Q 1.5*(v_qi_d − v_di_q)计算瞬时功率。有功频率环P_ref和P_e做差经过增益模块1/(J*ω0)得到dω/dt送到Integrator模块输出ωω再与ω0做差乘Dp形成阻尼反馈ω直接送到第二个Integrator得到相位θ。无功电压环Q_ref和Q_e做差送PI控制器含限幅输出加上E_ref得到虚拟电动势E。参考电压生成用E和θ生成三相正弦参考电压公式就是前面的sin函数表达式。电压电流双闭环外环电压PI输出电流参考值内环电流PI输出调制波最后接PWM Generator载波频率10kHz。电压电流双闭环的参数整定有一个经验顺序先整定内环电流环再整定外环电压环。电流环PI的带宽要高于电压环3~5倍我常用电流环Kp2、Ki100电压环Kp0.2、Ki20这个量级。这个参数不是万能的但作为初始值起步非常有效。电感电容参数也要提一下。LCL滤波器参数我取逆变器侧电感L10.5mH网侧电感L20.2mH滤波电容Cf20μF阻尼电阻Rd2Ω。LCL滤波器的谐振频率要避开PWM载波频率和基波频率按这个参数算下来谐振频率在2kHz左右和载波10kHz拉开距离和基波50Hz更远留足了裕量。4. 调试、验证与经验总结4.1 三大类典型仿真问题与排查第一类问题是仿真直接发散。这种情况90%出在求解器设置和代数环上。求解器方面换成ode23tb或ode15s最大步长缩到1e-4基本能解决大半发散问题。代数环出现的原因是功率计算、控制算法和输出环节之间存在直接反馈回路而没有延迟。排查方法是在反馈通路上加一个Memory模块或Unit Delay模块断开代数环代价是引入一拍延迟相位裕度会略降但仿真稳定性大幅提升。第二类问题是低频功率振荡。这种振荡的频率通常在几赫兹到十几赫兹波形上就是有功功率和频率来回摆动振幅缓慢衰减或持续振荡。低频振荡根源几乎都在J和D的配合不当。J太大、D太小系统像“没有阻尼的秋千”越荡越厉害。我的调整策略是保持J不变逐步增大D看振荡是否收敛如果D已经很大还振荡再减小J。记住一个原则J决定振荡幅度D决定衰减速度。第三类问题是并网瞬间冲击电流过大。这个问题的根源是预同步没做好。排查时先看并网前VSG输出电压幅值和电网电压幅值是否一致相位是否趋于一致。我在模型里用了一个简单的预同步控制并网指令发出前0.2秒通过PLL锁电网相位用PI调节器把虚拟相位慢慢“拉”到电网相位附近幅值误差控制在5%以内再合闸。加了预同步之后冲击电流从几百安降到几十安效果非常明显。4.2 仿真结果怎么看频率响应、功率响应、电压响应模型跑通之后判断仿真是否正确不能只看“没发散”要观察三类波形。第一类是功率响应。设置一个阶跃第0.5秒并网第1秒把P_ref从20kW阶跃到40kW观察VSG输出的有功功率。正确响应是有功功率快速上升经过一个较小的超调后稳定在40kW整个调节过程持续约0.3~0.5秒。如果超调大且持续振荡说明阻尼不足或J偏大。第二类是频率响应。在电网侧加频率扰动比如第2秒电网频率从50Hz下降到49.8Hz。此时VSG输出的有功功率应该自动上升这就是虚拟同步机的一次调频响应。你要观察的是功率增加量是否和下垂系数对应以及频率过渡过程是否平滑。J越大功率变化越平滑但达到新稳态的时间也越长。第三类是直流母线电压响应。让风速从10m/s突变到8m/s风机输出功率瞬间减少此时储能应该迅速释放功率维持母线电压稳定。如果母线电压波动在±5%以内且能快速恢复说明储能控制带宽和VSG响应配合得不错。如果母线电压出现持续性振荡大概率是储能电流内环PI参数偏大。4.3 几个让模型更贴近工程的细节第一个细节是加限幅和抗积分饱和。VSG输出的虚拟电动势E要限幅调制波要限幅无功功率积分器要限幅储能电流环输出要限幅。所有PI控制器都要用带抗积分饱和的变体否则在限幅期间积分器会无限累积退出限幅后系统会出现严重的超调甚至发散。Simulink里直接选PID Controller模块把积分限幅选项打开就行。第二个细节是风机MPPT和VSG的协调。风机MPPT按最大功率输出但VSG在频率支撑时需要增加有功出力两者可能出现矛盾。工程上通常的做法是让风机按MPPT运行频率支撑功率由储能提供VSG只负责接口控制。我在模型中的处理方式是P_ref由上层调度给定不直接跟随风机的MPPT功率而是跟随“风机当前出力储能当前出力”这个总量这样能保证直流母线功率平衡。第三个细节是记录关键状态量用于分析。建议在模型中用Scope模块同时观察直流母线电压、三相电流、VSG虚拟角速度、输出有功功率、输出无功功率这5个信号。单看其中任何一个都很难发现系统问题但放在一起看问题原因会清楚很多。尤其是虚拟角速度和有功功率两个波形放在一起几乎能判断出所有控制参数的调节方向。第四个细节是参数的可视化配置。建议把风机容量、储能容量、直流母线电压、滤波器参数、VSG的J和D、PI参数都放在MATLAB工作区或数据字典中不用每次打开子系统翻模块参数。我在模型初始化脚本中用结构体统一管理参数这样换一组容量参数时只改脚本里的几个数就够了模型本身一行不用动。模型初始化脚本会在仿真开始时自动运行用set_param和base workspace变量传递参数方便又不会出错。踩过几次坑之后我最大的体会是VSG风储仿真70%的时间花在参数整定和模型调试上并不夸张核心公式其实就那么几个。但正是这几个公式对应的物理意义决定了你能不能把“虚拟同步”这个理念变成一套稳定运行的仿真系统。我搭这套模型时经历了从参数随手设导致的长时间振荡到系统地按惯性时间常数和阻尼时间常数整定后一次跑通的转变。如果你现在正被振荡和发散折腾别急着改一堆参数先弄清楚每个参数对应的物理意义再有针对性地调效率高很多。