做风电并网仿真的人，大概率绕不开一件事：新能源占比高了以后，电网的惯量和阻尼被稀释得很厉害。风储VSG，也就是把虚拟同步发电机算法用在风储联合系统上，这几年在Simulink仿真里特别火。简单说，就是用储能配合风机，让逆变器在机电暂态尺度上模仿一台真实同步发电机的转动惯量和阻尼特性，从而让新能源场站具备主动支撑频率和电压的能力。这篇文章我就以一套典型的“风储+VSG”并网仿真模型为例，把从原理、建模、参数计算到波形分析的完整过程拆给你看。不论你是正在做毕业设计，还是刚接触新能源并网仿真的工程师，只要按着这个思路搭，至少能省出两三个星期的试错时间。
1. 为什么风储系统非要引入VSG？——虚拟同步发电机原理拆解
1.1 风电并网的“惯量危机”到底是怎么回事
传统电力系统的主力是火电和水电机组，它们都是旋转设备，转子本身带着巨大的转动惯量。当系统出现功率缺额、频率开始下降时，这些转子会靠自身动能释放能量，延缓频率跌落的速度，给一次调频和备用电源争取时间。风电和光伏接入后，电力电子变流器把风机和电网完全解耦，风机转子的动能根本“看不到”电网频率的变化——这就是所谓惯量缺失。用大白话说，电网像一辆大货车，同步机是自带刹车蓄能的重载卡车，而新能源并网逆变器就像一个没有方向盘助力和缓冲的小推车，反应快但压不住冲击。
逆变器如果没有额外控制，本质上就是一个“跟屁虫”：锁相环测电网相位，然后按功率指令往电网灌电流。电网频率稍微波动，它没有任何天然的惯性反作用力。风机数量一多，整个系统的等效惯量被稀释，频率跌起来又急又深，传统同步机的调速器还没反应过来，保护就动作了。这就是新能源并网考核里为什么反复强调“惯量支撑”的原因。模拟同步机不是要复刻它的转子，而是把它的“对外行为特征”用控制算法复现出来。
1.2 VSG的核心思路：把逆变器“伪装”成同步机
虚拟同步发电机（VSG）的核心思路，是在逆变器控制算法里加入一个同步电机数学模型。同步电机对外表现出来的两条最关键特性，一条是有功-频率调节特性，另一条是励磁-电压调节特性。VSG就通过这两条通道，分别模拟同步机的转子运动方程和励磁调节规则。最常用的转子运动方程是：
J * (dω/dt) = P_ref - P_e - D*(ω - ω0)
其中J是虚拟惯量，D是虚拟阻尼，P_ref是给定有功功率，P_e是逆变器实际输出有功功率，ω是虚拟角速度。这方程说明一件事：输出功率和指令不等时，角速度就会变化，而这个变化过程带着惯性、带着阻尼，不会瞬间跳变。
加上这一层之后，逆变器的输出频率不再由电网锁相环硬性决定，而是由这个“虚拟转子”的动能来决定。电网频率跌了，它会像真实同步机一样把储能或者直流母线电容里的能量释放出来，先顶一下频率；功率平衡恢复后，再慢慢回到额定值。这套机制让VSG做的波形非常“有电机味”，尤其在弱电网、孤岛以及黑启动场景里，它能表现出和常规跟网型逆变器完全不同的坚韧度。
1.3 储能在VSG架构里扮演的角色
单靠风机做VSG有点力不从心。风机出力随风速波动，很多时候根本没有足够的“腾挪空间”来提供惯量和调频功率。比如风速低了，风机自身的功率都不能满足并网额定出力，哪还有余量去顶频率？储能恰好可以弥补这一点：电池响应速度快，双向DC/DC可以快速充放电，在VSG需要瞬时功率支撑时立刻出力，等风功率恢复或者负载变化后再把电补回来。
所以在仿真架构里，储能的作用是“压舱石”和“缓冲池”。风机负责输出主功率，储能负责吸收功率波动、支撑直流母线电压、在VSG功率指令变化时提供瞬时能量。整个系统的控制链路就变成了：风机发电、储能调峰调频、VSG模拟同步机外特性。仿真中储能通常用锂电池加双向Buck/Boost实现，通过功率外环控制电池充放电，使直流母线电压稳定，同时满足VSG并网功率指令。
2. 风储VSG仿真系统整体架构与模型选型
2.1 拓扑选择：背靠背直流母线方案与交流母线方案
风储联合的结构形式很多，仿真里最常见的其实是共直流母线拓扑。风机经整流器或者机侧变流器接到直流母线，储能通过双向DC/DC也接到同一条直流母线，然后共用一台VSG逆变器并入电网。这种结构的好处是只有一个并网接口，VSG算法只需作用在那台逆变器上，控制链路清晰，调试也方便。
另一种是共交流母线拓扑，风机、储能都通过各自的变流器挂在交流母线上，再统一升压并网。这种方式更贴近实际风电场，但仿真复杂度会高很多，要处理多台变流器之间的功率均流、模式切换和环流抑制。如果你的研究目标是VSG本身的惯量支撑特性和并离网切换，强烈建议先用共直流母线方案，先把控制逻辑跑通、波形看明白，再向多机并联扩展。直接在初期就上交流母线多机并联，光调均流就能把人磨没了。
2.2 控制框架：有功-频率VSG、无功-电压下垂、储能功率分配
VSG控制框架大致分成三块：有功-频率控制、无功-电压控制、底层电压电流内环。有功-频率这一块就是前面说的转子运动方程，输入是功率指令与实际输出之差，输出是虚拟频率和功角。无功-电压部分采用Q-V下垂，模拟同步机的励磁系统，输出虚拟内电势幅值。最后把虚拟内电势幅值和相位合成为一个三相正弦参考电压，再交给电压电流双闭环，生成逆变器PWM驱动信号。
储能那边的控制相对独立。检测直流母线电压，与参考值做差，PI控制器输出电流参考，再经过电流内环控制Buck/Boost占空比，维持母线电压稳定。如果想做上层功率调度，比如风储联合平滑出力，可以在储能功率外环加上限幅和功率分配模块。注意一点：储能和VSG不能各调各的，否则直流母线电压和VSG功率指令会产生耦合振荡。我的做法是在VSG有功指令里加了一个低通滤波，同时让储能母线电压控制器的带宽低于VSG虚拟惯量带宽，这样两条链路不容易打架。
2.3 储能系统建模：锂电池、双向DC/DC与SOC估算
Simulink里锂电池建模建议直接用Simscape Electrical自带的Battery模块，不用自己写复杂的电化学方程组。但几个关键参数必须设置准确：额定电压、额定容量、SOC初值、内阻和充放电特性曲线。Battery模块的SOC输出可以直接连到仿真示波器里观看，方便排查异常。
双向DC/DC通常用半桥结构实现，放电时Boost升压给直流母线供电，充电时Buck降压给电池充电。控制上需要判断母线电压偏差方向来决定工作模式。如果嫌麻烦，可以用Simscape Electrical里的Two-Quadrant DC/DC Converter模板，里面已经封装好了可控开关和反向二极管，直接接控制信号就行。SOC估算用安时积分法，仿真时要注意容量单位换算：电池容量如果是Ah，要换算成As再积分。我见过不少人在这里翻车，SOC在几十秒内从70%跳满或者跳零，根源往往是单位没换算对。
3. Simulink仿真模型搭建实操：从零开始拼系统
3.1 软件版本与基础仿真设置
我这里用的是MATLAB R2021a以上版本，Simulink自带的Simscape Electrical工具箱基本够用。如果你装的是学生版，先确认License里有没有Simscape Electrical，否则后续很多电气元件都拖不出来。可以在MATLAB命令行敲 ver 命令检查工具箱列表。
仿真类型建议这么设：主电路部分用连续模式，控制部分用离散模式，采样频率设在10kHz左右。不要整个模型都用固定步长1e-6跑，那样速度慢到怀疑人生。求解器选ode23t或者ode15s，采用变步长，最大步长限定为开关周期的1/10左右，既保证稳定性也兼顾速度。我最初用全固定步长1e-5跑，一个10秒工况要半小时；后来改成混合步长方案，同样工况三分钟完事，波形也几乎一样。
3.2 风电机组模型的简化：要不要用PMSG/DFIG？
如果你研究的是风储VSG整体外特性，重点在惯量支撑和并网控制，那风机模型不必上完整的气动+PMSG+最大功率跟踪那一套。完整风机模型包含风速模型、桨距角控制、转矩控制、发电机和整流器，仿真步长会被拖到微秒级，极其折磨。
我的简化做法是：用一个可控直流源配合风速-功率查表模块，模拟风机输出功率随风速变化。风速信号可以用阶跃、斜坡或者随机风，查表得到机械功率，再经过一个一阶惯性环节（模拟风机转子惯性）接到直流母线。如果后面要研究风机内部动态对VSG的影响，再按需加入PMSG模型。记住：先简后繁，先把VSG调通，再往里塞细节，否则报错满天飞。
3.3 储能双向DC/DC模块：快速建出来
在Simscape Electrical里，最快的方式是用两个理想开关和两个二极管组装一个半桥变换器。上桥臂IGBT配合反并联二极管，下桥臂同样，电感接在中间点，一端连电池，另一端通过半桥接直流母线。控制信号用PWM发生器产生。
控制策略写起来并不复杂：测量直流母线电压，与参考值500V比较，误差经过PI控制器得到电流参考；再与电池电流比较，电流误差经过PI控制器得到占空比。放电模式下PWM作用于下管，充电模式下PWM作用于上管。模式切换用Stateflow比较方便，但简单模型里用一个条件判断模块也可以，判断母线电压偏差的符号来选择输出到哪路PWM。仿真中要小心两个桥臂不能同时导通，PWM信号要有死区，否则开关器件会直通短路。
3.4 VSG逆变器主电路与LCL滤波器搭建
逆变器直接用Universal Bridge模块，选IGBT/Diode，三相桥结构。直流母线正负极接储能电容，容量一般取2200μF到4700μF，根据系统容量和动态响应要求调整。输出侧接LCL滤波器，具体参数我放在后面一章详细给推导。注意给IGBT设置合适的开通关断延时和缓冲吸收电路，否则仿真中电压尖峰会直接吹爆模型。
LCL滤波器由逆变侧电感L1、网侧电感L2和滤波电容Cf构成。一个最容易忽略的点是滤波电容支路必须加阻尼——要么串阻尼电阻，要么做有源阻尼控制。不加阻尼的LCL在一个特定频率处会出现负谐振峰，逆变器等效阻抗在那里变成负值，电流波形会高频振荡。仿真中如果你看到并网电流叠加了一股1kHz左右的衰减振荡，十有八九就是LCL阻尼没做好。
3.5 VSG控制算法：S函数还是Simulink模块？
VSG核心算法用S-Function写最灵活，也可以用Simulink基础模块搭。我更推荐用MATLAB Function或者S-Function Builder，因为转子运动方程是一组微分方程，需要做积分运算，用纯模块搭初期值和延迟很容易出问题。核心代码很短，基本就是欧拉积分。
下面给个示意性的伪代码，实际使用时还要加限幅和低通滤波：
if t == 0 theta = 0; % 初始功角 omega = 2*pi*50; % 初始角速度 end % 虚拟转子运动方程 dtheta_dt = omega - omega0; domega_dt = (P_ref - Pe - D*(omega - omega0)) / J; theta = theta + dtheta_dt*Ts; omega = omega + domega_dt*Ts; % 励磁调节 E_ref = E0 + n*(Q_ref - Qe); % 三相参考电压 Va = E_ref*sin(theta); Vb = E_ref*sin(theta - 2*pi/3); Vc = E_ref*sin(theta + 2*pi/3);在实际模型里，P_ref和Q_ref来自上层调度或下垂计算，Pe和Qe来自功率测量模块。参考电压生成后要经过限幅和滤波，再送给电压电流双闭环。电压电流内环用PI控制器，带宽设为开关频率的十分之一到五分之一，PWM载波频率取10kHz左右。
3.6 并网开关、三相变压器与电网模拟器
并网开关用Three-Phase Breaker模块，控制信号可以用一个Step在固定时间合闸。合闸之前建议先让VSG做预同步，等相位和幅值都对齐了再合闸，否则冲击电流能冲到额定值的几倍。变压器用Three-Phase Transformer模块，电压等级比如0.69kV/10kV，连接组别常用Dyn11。这个连接组别千万别选错，否则参考相位和电网侧相位之间会差30度，合闸瞬间必炸仿真。
电网用Three-Phase Source模拟理想无穷大母线，内阻设小一点比如0.001Ω。想模拟弱电网场景，就把电网阻抗调大，或者用单独的RL串联阻抗网络代替源内阻。注意电网侧和逆变器侧必须有公共参考地，变压器中性点接地方式也要设对，否则Simulink会报“floating neutral”错误。
4. VSG关键参数计算与调参实战记录
4.1 虚拟惯量J和阻尼系数D怎么算
虚拟惯量J不是随便拍脑袋定的，它要对应同步机的惯性时间常数H。H的定义是转子储能除以额定容量，同步机一般在2到10秒之间。换算公式是：
J = 2HS_ref / (ω0^2)
比如额定容量S_ref=100kVA，额定角速度ω0=2π*50，取H=4秒，那么J撑死也就三千分之几，但代入实际数值算出来约等于8.1 kg·m²。很多人对这个单位没概念，确实，虚拟惯量不是真实物理转子，但它在数学上模拟的就是转子惯量的效果，所以数值级别可以参考真实同步机，再根据仿真波形微调。
阻尼系数D代表功率偏差对角频率变化的阻尼作用。D太大，系统会显得“拖沓”，频率恢复慢但超调小；D太小，频率振荡衰减慢，甚至持续波动。经验上可以先取D = J * 10~20倍，再看系统闭环极点。更严格的做法是把J和D代入小信号模型，画根轨迹，按期望阻尼比0.7~0.8来整定。工程上更多是仿真中直接试：先把J定住，从小到大调D，看频率曲线的超调量和收敛时间；再反过来微调J。我最终用的是J=8、D=120，频率波动两个周期收敛，效果很理想。
4.2 有功下垂系数m与无功下垂系数n的标幺化计算
有功下垂系数m = Δf / ΔP，表示有功功率每变化多少，频率允许偏移多少。频率偏差限制通常为0.5Hz，有功功率允许变化100%，那么标幺值下 m_pu = (0.5/50)/(100/100)=0.01。在Simulink有名值控制模型里，如果直接用实际频率误差乘以m，量纲就是Hz/W，很容易出错。所以建议把反馈量先归一化，P_ref和Pe都除以额定功率，频率差除以50Hz，再做下垂计算，这样可以避免单位混乱。
无功下垂系数n类似，n = ΔU / ΔQ。电压偏差允许5%，无功变化100%，那么n_pu = 0.05。无功-电压实际控制中还要考虑线路感抗导致的电压降落，所以在并网点要取电压有效值作为反馈。Simulink里测有效值需要动态测量模块，注意设定好低通滤波时间常数，通常20~50ms，否则电压波动反馈进去会加剧振荡。
4.3 LCL滤波参数设计经验公式
LCL滤波器的谐振频率fr需要避开控制带宽，也要避开开关频率及其边带，一般取开关频率的1/10到1/5。已知开关频率fsw=10kHz，取谐振频率fr≈1kHz。电感宽度由电流纹波限制决定，用公式：
L_total = Udc / (8 * Δi_max * fsw)
Δi_max一般取额定电流峰值的10%~20%。假设直流母线电压700V，额定电流峰值参量先算出来，把L_total分成L1和L2，通常L1:L2=1:0.5~0.7。电容值再用谐振频率公式反推：Cf = 1 / (L_total * (2π*fr)^2)。
举个例子：100kVA系统，线电压380V，额定电流约152A，峰值215A，纹波取30A，那么L_total = 700 / (83010000) ≈ 0.29mH，取300μH。L1取180μH，L2取120μH，Cf由谐振频率1kHz反算约846μF，这个值偏大，说明要把谐振频率抬高到2kHz左右。实际中还要考虑电容无功损耗，一般Cf不超过系统总容量的5%。所以LCL参数设计是个反复迭代的过程，仿真中如果发现并网电流高次谐波大，先调电容支路阻尼，再动电感值。
4.4 储能容量与SOC初始值设定技巧
储能容量要和VSG所需支撑时间匹配。比如系统在惯量支撑阶段需要额外提供10%额定功率、持续5秒，那需要的能量就是额定容量×10%×5s/3600s。以100kVA系统为例，10%功率是10kW，支撑5秒就是14Wh，对应50Ah/200V电池大约能支撑5分钟级别的短时功率支撑，远够用。所以仿真里电池额定容量不用选太大，否则SOC变化曲线会平淡得看不出控制效果。
初始SOC建议设在70%左右，同时保留向上和向下调节空间。如果初始SOC设成100%，充电过程就被限幅了，电池只能放不能充，模拟结果会失真；设成0%则放电能力不足。仿真时还要注意DC/DC的电流方向定义：Battery模块默认输出电流为正，放电时电流流出电池，SOC下降；充电时电流流入电池，SOC上升。如果发现SOC曲线方向反了，查一下电流测量模块的正负连接。
5. 仿真结果分析与波形解读（含典型工况演示）
5.1 风速阶跃下系统频率与有功响应
我做的第一个典型工况是风速在1秒时从10m/s突变到12m/s，风机输出功率从80kW爬向100kW。如果风机直接整流并网，直流母线电压会瞬间抬升，并网功率突变，频率几乎跟着电网走，没有任何支撑现象。加了VSG之后，P_ref从80kW平滑过渡，由于虚拟惯量J的存在，输出有功出现了几十毫秒的斜坡，频率偏差也被限制在允许范围内，这就是惯量支撑的直接体现。
这里要提醒一下：风机简化的功率接口不能太“硬”。风速阶跃本身变化剧烈，如果查表模块直接输出功率阶跃，VSG输入的P_ref就跟着突变，虚拟角速度会产生剧烈振荡。更自然的做法是在风机功率输出后面加一个时间常数5~10秒的一阶惯性环节，模拟风机转子储存的动能。这样风速阶跃后，风机输出功率是渐变过程，VSG的响应才有真实的物理感。
5.2 负载突变时VSG的频率恢复能力
本地负载在1.2秒增加了50kW。并网场景下，电网相当于无穷大母线，VSG输出频率由电网牵制，所以频率几乎不变，但有功功率振荡会体现在功率波上。离网场景下，频率偏移会更明显：负载增加瞬间，VSG输出有功跟不上，虚拟转子减速，频率降下来；随后储能快速输出功率，频率回升到新稳态。观察虚拟角速度ω和并网点频率波形，可以看到典型的二阶系统欠阻尼响应，有衰减振荡。
调参时最直观的方法是看频率超调量和恢复时间。如果系统欠阻尼，频率会振荡好几周才收敛，那就增大D；如果频率变化太慢、动态功率支撑滞后，那就减小J。总之要记住J管“惯性”、D管“阻尼”，两者共同决定VSG的动态响应形状。不要看到振荡就盲目降低所有增益，那会把响应调到又慢又软。
5.3 并网转孤岛运行模式切换表现
这是VSG最能体现价值的地方。并网状态下电网相当于无穷大母线，VSG功角随着电网频率走；当并网开关在2.0秒断开，系统必须无缝转孤岛运行，由风机、储能和本地负载组成独立微网。如果VSG控制里一直依赖PLL锁定大电网相位，脱网瞬间会失去参考，输出电压幅值和相位直接崩溃。
正确做法是让VSG的虚拟角速度直接积分得到相位，不依赖PLL。这样脱网后，VSG依靠自身的转子运动方程和储能功率输出维持电压和频率。仿真里把并网点开关设为2.0秒断开，观察三相电压电流波形。我实测下来，切换瞬间电压波形基本连续，没有冲击电流；频率从50.05Hz先降到49.85Hz，再经过0.6秒恢复到约50Hz。这个结果说明储能和VSG的控制切换逻辑是有效的。
6. 仿真中高频踩坑问题与排查技巧
6.1 仿真心步长太小导致速度慢、卡死
很多新手会把整个模型设置成固定步长1e-6，结果跑1秒仿真要半小时，最后直接失去耐心。实际上主电路开关频率10kHz，步长设为1e-5或5e-6完全够。IGBT开关瞬间的电压电流尖峰用内部小步长去解，不需要全局都采用极微小步长。推荐全局用变步长，最大步长限制在1/10开关周期；如果你需要观察开关高频分量，再单独在功率模块内部设定局部采样时间。
如果仿真卡死不收敛，优先检查是否存在代数环。最常见的情况是MATLAB Function的输出直接又作为输入反馈回去，Simulink会报“Algebraic Loop”错误。解决办法是加一个Memory模块或者单位延迟(ZOH)打断这个环。另一个常见原因是电网阻抗或变压器漏感太小，导致数值刚性问题，此时把最小步长调小一点，或者换用ode23t求解器，往往会有奇效。
6.2 电流波形震荡发散，先查寄生参数还是控制参数？
遇到电流波形发散，先稳住心态，按顺序排查。第一步查主电路：LCL参数是否合理，电容支路有没有加阻尼，直流母线电容是否足够大。第二步查控制参数：电流环PI带宽是否过高，通常建议电流环带宽是开关频率的1/10到1/5；电压环带宽再低一个数量级，否则内外环会打起来。第三步查延迟：数字控制里采样、计算、PWM更新都有延迟，如果不补偿，相当于在环路里引入了一个相位滞后的纯延迟环节，很容易振荡。
我在仿真里遇到过一个问题：电压电流双闭环里我用了同一个采样触发信号，导致电流环采样值和PWM占空比更新之间存在一拍延迟，高频振荡一直消不掉。后来在电流环输出加了个一拍延迟补偿，振荡立刻消失。这个细节在校验模型时很容易被忽略，但影响很大。
6.3 VSG和锁相环PLL到底什么关系
很多刚接触VSG的人会混淆概念，以为VSG也要靠PLL锁电网相位才能工作。实际上VSG属于构网型控制，它的相位由虚拟转子运动方程自己产生，不需要锁电网相位。而传统跟网型逆变器则必须用锁相环测电网相位，按照电网相位去控制输出电流。这个区别决定了VSG可以工作在弱电网甚至孤岛场景，跟网型逆变器在这些场景下稳定性会急剧下降。
但VSG也不是完全不用PLL。在并网预同步阶段，需要把VSG输出电压相位和电网相位对齐，这时一定要用PLL测量电网相位，通过前馈或积分校正让两者一致，再合闸。并网正常运行后，可以断开PLL前馈，让VSG自由运行。很多好的Demo仿真里会保留PLL用于实时监测，但核心控制回路不能锁死电网相位，否则VSG就退化成了带功率外环的跟网型控制，所有构网特性荡然无存。
6.4 储能SOC偏离预期，问题出在哪
SOC跳变异常一般有三个原因。第一，电池容量单位搞错：Battery模块里容量单位默认Ah，安时积分法里电流单位是A，积分时间单位是s，需要换算成As，否则积分结果大错特错。第二，初始SOC设置不合理：设成100%导致充电限幅，设成0%导致放电限幅，仿真会失真。第三，DC/DC控制方向反了：电池充放电电流符号颠倒，SOC曲线违反能量守恒，这是仿真里最隐蔽的错。
排查方法是加一个SOC观测示波器和电池电流波形，同时播放模型运行过程。如果SOC曲线上升，但电池电流显示放电，那肯定有反符号问题。把电流传感器的方向或者PI控制器输出极性调一下即可。另外，储能容量如果设置太充裕，10秒仿真里SOC基本一条直线，看不出变化，也容易让人误以为控制没生效。建议把容量设得相对紧张一些，让SOC动态曲线有可见的升降，这样验证控制器响应更方便。
最后分享一个我自己在做这个仿真时的体会：很多第一次接触VSG的人会把注意力全放在控制算法代码上，反复调J和D，却忽略了整个风储系统的功率平衡链路。实际上，如果不把风机的功率波动、储能的SOC约束和直流母线电压控制当作一个整体去建模，VSG的波形再漂亮也没法落地。仿真最大的价值其实是帮你把“风机-储能-逆变器-电网”这条能量链路的逻辑想清楚，参数只是表现。建议你先把简化模型跑稳定，再加风速场景和并离网切换，这比一上来就堆PMSG、六脉波变压器、完整气动模型高效得多。希望这份记录能让你少走一些弯路。
