风储VSG这几年在新能源并网圈子里讨论度很高。我去年基于虚拟同步发电机做了一版风储并网系统的Simulink仿真把风电机组、储能电池、双向DC/DC、网侧变流器和VSG控制串成了一条完整链路解决了单纯风机并网不提供惯性支撑、抗扰动能力差的问题。很多做微电网或风电的同学来问我风机本身有变流器为什么还要专门加储能再用VSG控制这个问题恰好是整套方案的核心。风电出力随机波动常规变流器又把自己和电网频率“隔离”了电网频率变化时风机只能当旁观者而VSG的本质是让变流器对外呈现出同步发电机的暂态特性——有虚拟转子、有惯量、有阻尼。储能在这里扮演“虚拟原动机”的角色补上功率缺口变流器则扮演“虚拟同步机”在并网点把风机和储能整理成一个可控的、带惯量支撑的电源。这篇内容适合正在搭新能源并网仿真、调VSG参数或者毕设卡在频率支撑这一块的朋友我会从原理拆到模型再讲参数整定和调参踩坑力求可以直接照着复现。1. 风储VSG的整体思路与原理拆解1.1 为什么风机侧必须加储能先算清功率缺口先说结论没有储能VSG控制做得再漂亮也是空架子。风电并网最大的麻烦是风功率不可控。风速变化时风轮捕获的机械功率直接变化但电网不会因为风速降了就对风机宽容它只关心并网点的功率是否平衡、频率是否稳定。传统直驱或双馈风机通过背靠背变流器并网转子机械侧和电网电气侧在控制上解耦风机基本不响应电网频率变化。换句话说风机本质上是“输出功率跟着风速走”的电流源不是“电压支撑型”的同步电源。加入储能以后能量平衡关系变成P_wind P_batt P_grid。风速快速下调时电池立刻放电填补功率缺口风速超过限值时电池吸收多余功率。储能还要维持直流母线电压稳定因为VSG输出的有功等于电网吸收的有功瞬时差额如果没人承担母线电容电压就会暴跌。我这边配置储能容量一般按风机额定功率的15%~20%来选支撑时间2~10秒不是问题具体看所需惯量支撑时长。容量太大浪费太小则VSG的虚拟惯量效果发挥不出来。功率缺口可以量化算一下。假设一台50kW风机从额定风速降到0.8倍额定风速按理想立方关系捕获功率大约降到50×0.8³25.6kW缺口接近24.4kW。如果要求支撑3秒需要电池放出的能量约为24.4kW×3s73.2kJ考虑DOD限制和DC/DC效率电池有效容量至少得按几十到一百Wh级别留裕度。这样算完你才知道储能容量和VSG支撑能力之间是什么关系不会随便拍脑袋选数。1.2 虚拟同步发电机核心方程J和D到底在模拟什么虚拟同步发电机并不是一个新电机它只是一套控制算法外加一个功率变换器。它把同步发电机的转子运动方程写进控制环最核心的有功-频率关系是J dω/dt P_ref - P_e - D(ω - ω_0)其中J是虚拟惯量ω是变流器输出电压的角频率P_ref是参考有功P_e是实测并网有功D是虚拟阻尼ω_0是额定角频率。这个式子跟同步发电机的摇摆方程同构。电网频率掉下去ω小于ω_0但P_ref还是那么大公式右边变成正dω/dt为正说明控制器会让频率按一定速率回升而不是瞬间跟着电网频率躺平。J越大频率变化率越慢惯量感越强D越大振荡衰减越快但稳态时会出现一个频差偏差。无功-电压部分模拟同步机的励磁调节常用形式是E E_0 k_q(Q_ref - Q_e)E是变流器内电势幅值Q_e是并网无功。电网电压跌落时Q输出自动增大变流器就有类似同步机励磁调节的外特性。相位θ再通过对ω积分得到。这样变流器输出的就是一个幅值由无功环决定、相位由有功环决定的电压源并网侧表现为低阻抗电压源而不是电流源。有人会问为什么不直接用传统的PQ下垂控制下垂控制相当于一个带阻感输出阻抗的电压源虽然有静态调差但没有J这一项扰动瞬间的频率变化率直接由电网决定。VSG引入了二阶动态多了一个可调的“虚拟转子惯性”这让变流器能主动扛一下频率变化给一次调频或储能动作争取时间。这是风储VSG和普通下垂控制最本质的区别。1.3 为什么用Simulink而不是PSCAD或Matlab纯M脚本Simulink在我这里属于控制算法验证的第一选择。模型库里有Wind Turbine、Asynchronous Machine、Battery、Universal Bridge、Discrete PWM Generator这些现成件主电路和控制系统在一个环境下跑调试时能直接看到中间变量。PSCAD在电磁暂态细节上确实强但要建风机储能VSG完整闭环尤其是做控制系统参数扫描Simulink更顺手。后面要做C代码生成Simulink Embedded Coder可以直接把VSG控制子系统生成成C代码扔进DSP里跑这一点比PSCAD方便得多。Simulink另一个好处是“虚实结合”。主电路放Simscape Electrical控制部分用离散模块后面做RCP或半实物仿真时改动很小。我之前试过用Matlab脚本批量跑不同风速曲线和参数组合Simulink配合并行计算一个工况矩阵能在一个上午扫完。做毕设或者写小论文这种方式输出图表的效率特别高。2. 系统架构与关键模块设计2.1 系统拓扑与能量流从风轮、电池到电网我用的拓扑是比较经典的风储一体化并网风电机组输出经不控整流或背靠背变流器送到直流母线储能电池通过双向DC/DC挂到同一母线上母线再经过三相两电平逆变器和LCL滤波器并入电网。双馈风机的转子侧本身有交流变流器我这里说的风储共用直流母线是一种简化的机组级方案适合做VSG研究。能量流是闭环的风轮机械功率→发电机功率→直流母线→VSG逆变器→电网。风速下降时前级功率减少但VSG按指令继续输出功率差额由电池补上风速上升时电池吸收。母线电容只是一个很小的能量缓冲真正扛功率差额的是电池。所以储能控制带宽一定要快否则暂态过程母线会出大坑。我在做拓扑时最看重两个接口一个是直流母线电压的稳定性一个是并网逆变器PWM的相位连续性。这两条线一旦打通后续调VSG参数基本不会崩。2.2 风电侧建模MPPT查表还是全阶模型这里有两种选择。第一种用Simulink里的Wind Turbine模块加Asynchronous Machine再做背靠背变流器最接近工程但模型复杂度高转速环、桨距角、转子侧和网侧PWM全都要调试。第二种也是我更推荐的用MPPT功率源来等效。所谓风储VSG重点在“风储”外部特性和并网支撑不在风机自身的电磁暂态细节所以我会建一条MPPT功率曲线随风速变化输出P_ref。MPPT曲线按公式P 0.5ρAC_p(λ,β)v³计算C_p可查表简化。仿真中用一个一阶惯性环节缓冲一下模拟风轮响应避免风速阶跃时机电功率突变太生硬。很多模型直接让风电功率瞬间跳变结果VSG和储能还没反应过来母线已经断了。我自己习惯在风电功率输出后面加一个时间常数0.1~0.2秒的低通滤波效果很接近真实风轮惯性。准备一个小技巧风速发生器用随机风速叠加阶跃这样比单纯阶跃更有说服力波形看起来也更有“现场感”。2.3 储能侧建模电池加双向DC/DC稳住母线电池用Thevenin等效模型就够了理想电压源加内阻加RC网络Simulink里的Battery模块能直接干这活。SOC按SOC初始SOC减∫(I/Q_n)dt得到初始SOC我一般设在50%~60%给双向波动留足空间。双向DC/DC用Buck/Boost实现控制目标就是稳住直流母线电压。外环读取母线电压与参考值400V或600V比较输出电流指令内环控制电感电流。这里要特别注意方向问题。电池放电和充电时电流方向不同控制器符号处理错了会出现“越稳压越崩”的情况。我踩过这个坑后来在电流内环里加了一个符号判断或者直接用带符号的参考电流问题就消失了。储能功率指令不直接给固定数而是由直流电压误差自动产生好处是VSG功率变化时储能能自动响应不用额外加协调算法这对仿真初期调试非常友好。2.4 并网侧两电平VSG变流器与LCL滤波逆变器用三相两电平桥调制波由VSG生成PWM用三角载波比较或SVPWMSimulink里Universal Bridge加Discrete PWM Generator很方便。LCL滤波器参数按并网逆变器惯例选L1约0.3~1.5mHC约10~47μFL2约0.1~0.5mH具体看电网阻抗和开关频率。为保证高频谐振衰减我会加阻尼电阻或有源阻尼仿真时如果没加阻尼波形会出现明显振荡。LCL设计其实不用一开始就调到最优先用一个L滤波器让系统跑通再换成LCL优化谐波这是降风险的思路。VSG控制器输出的调制波频率跟随虚拟角频率相位由dq坐标下的功率计算决定。我这里没有额外做dq旋转坐标系变换的详细展开但提醒一下如果选择在dq坐标系里做电压电流双闭环Park变换的角度一定要用VSG的相位角而不是电网锁相环的角度否则整个闭环是乱的。很多新手在这里栽跟头。3. 仿真模型搭建实操要点3.1 模块化分层搭建别把模型堆成一个毛线球我不建议把整个模型放在一个页面上也不要一直改手算系数。我一般会建立一个大项目包含三块init脚本、主模型、后处理脚本。主模型里分四个子系统风电子系统、储能子系统、VSG控制子系统、并网测量与断路器子系统。子系统之间用Goto/From传信号结果用To Workspace保存最后统一在Matlab脚本里画图。调参最好都做成变量在init脚本里赋值。比如Pnom、f0、J、D、Kp、Ki、Kq这些全都定义为工作区变量。重新仿真只改脚本不会出现多层模型里到处找参数的情况。总线信号用Bus Creator打包Scope用Data Inspector查看。模型做到这个程度后面做参数敏感性分析时你只要写一个for循环改变量再反复调sim(model)就能跑完一组曲线。如果模型是毛线球结构这种扫描根本没法做。3.2 虚拟惯量J和阻尼D怎么选估算公式加仿真校验这个应该是问得最多的。先把小信号推导给出来把VSG的有功-频率环小信号化忽略调速动态可以得到J Δδ D Δδ K_s Δδ ΔP_m其中K_s是同步功率系数约等于额定有功与电压角度的关系。因此自然振荡角频率ω_nsqrt(K_s/J)阻尼比ζD/(2sqrt(JK_s))。直观结论就是J决定惯量和振荡周期D决定阻尼和收敛速度。工程整定我按三步走第一步根据额定容量决定J初值。如果采用标幺值体系J_pu一般对应虚拟惯性时间常数H取2~4秒比较合理。折算到有名值J2H S_n/ω_n²。举例50kW系统H2sω_n2π×50≈314.16rad/sJ≈2×2×50000/314.16²≈2.03kg·m²控制里直接写2就可以。第二步取阻尼比ζ≈0.7~1.1利用ζD/(2sqrt(JK_s))估算D。K_s不方便估算时就先给D一个标幺值30~50然后增大到波形无超调。我习惯固定J从D10扫到D100看频率响应曲线。第三步再根据最大频率变化率rocof要求回调J。rocof过大就增大J恢复太慢就减小J。参数初始化代码可以写成这样% 风储VSG参数初始化示例 Pnom 50e3; % 额定功率 50kW f0 50; % 额定频率 Hz J 2.0; % 虚拟惯量单位可视为kg·m^2/rad D 50; % 虚拟阻尼初始标幺值 Kq 0.02; % 无功下垂系数电压偏差修正这里J和D的单位在不同实现方式下会变我建议仿真时全部用标幺值防止单位换算错误。标幺值下J通常落在0.5~5之间D落在10~100之间范围不算大扫描很快。3.3 求解器、步长与PWM离散化的坑Simulink里Simscape Electrical模型建议用Powergui配置为离散仿真定步长Ts1e-5秒左右。PWM开关频率10kHz时每个周期有100个采样点够用。如果求快变步长ode23tb加离散PWM也可以但离散PWM要求模型是时序逻辑变步长可能引起抖动我不推荐新手这么干。母线电容和LCL高频环节会让ode45步长变得很小仿真速度反而更慢。控制环的采样周期可以设为1e-5秒功率计算里的滤波时间常数设1e-3秒左右避免高频分量进入VSG控制。测量三相电压电流时要注意把测量模块的采样率设成Ts的整数倍否则输出波形会有混叠。Simulink的Scope显示默认会二次采样看起来毛刺很重先检查这个再怀疑控制有问题。3.4 四个必测工况风速突变、负载阶跃、无功支撑、电网频率扰动我把测试工况整理成四类每次调参前先跑一遍比对着单一工况调要可靠得多。风速突变工况在0.5秒时风速从12m/s降到9m/s观察并网有功、母线电压、SOC补偿情况。本地负载阶跃工况在0.8秒时并网点投入一部分有功负载模拟负荷扰动。无功阶跃工况投入无功负载验证无功-电压环。电网频率扰动工况在三相电压源里设置频率偏移验证VSG惯量响应。每个工况都要看两类指标暂态最低点或最高点以及恢复时间。把这些指标记录成表格放在init脚本里做成参数结构体跑完保存到mat文件后面画图直接用。我在实际调试中发现只跑风速阶跃容易忽略电压支撑问题只跑有功阶跃又容易忽略无功环路。只有四个工况一起跑才能暴露J和D参数在不同扰动下的矛盾取舍。4. 典型仿真结果与波形分析4.1 风速下降工况频率支撑与有功响应风速下降后风机输出从50kW降到约25kWVSG在一开始还按原参考输出差额由电池放电提供直流母线短暂跌落但被储能电流内环拉回。并网有功开始下降但下降速率较缓没有出现功率阶跃。这就是虚拟惯量在频率响应上的表现。从波形上看频率最低点出现在扰动后约0.2秒跌了0.2Hz左右随后约0.5秒恢复。如果没有VSG采用普通下垂控制频率变化率明显更陡。我看惯量效果主要看两个指标一是rocof也就是df/dt二是频率最低点。J大时rocof小最低点不明显但恢复更慢D大时振荡次数明显减少。这个工况下储能SOC的变化曲线能直接证明电池参与了支撑如果SOC曲线纹波很大且直流母线电压有明显波动说明储能电流环带宽不够。4.2 负载阶跃工况端电压与无功响应并网点有功负载增加后VSG输出有功增加储能短暂放电无功环同时动作端电压略有下降但被无功补偿拉回。如果无功下垂系数k_q设置太小电压偏差大如果太大无功分配振荡。我在仿真里一般取每1%电压对应10%~20%额定无功的下垂系数之后微调。这个工况最容易看出无功环和电压环的配合。很多模型里VSG只写了有功环无功环直接接一个固定电压源那样并网电压支撑效果就看不见了。风储VSG的价值在于它既支撑频率也支撑电压所以无功-电压环一定要参与闭环。调试时如果电压跌落瞬间无功出力没有快速上升先检查Q_e反馈的符号和滤波时间常数。4.3 直流母线动态与SOC变化我那次仿真里1.5秒内SOC从50%下降到49.2%放电电量约与功率缺口积分一致。这说明电池确实扛住了暂态功率而不是靠母线电容硬撑。母线电压在扰动后能回到参考值附近电压波动幅度控制在2%以内储能双向DC/DC的双闭环是稳定器。如果SOC曲线变化方向和功率缺口不一致说明储能控制环路的参考方向反了。如果SOC变化很小但母线电压跌落很大可能是双向DC/DC的电流限幅太低电池放不出那么多功率。遇到这种情况我一般先把电流限幅调到额定电流的1.5倍再观察波形不要一上来就改PI参数。4.4 J、D参数敏感性对比小波、大波和恰到好处参数敏感性分析建议做成表格不断扫描J和D。我常用三组对比小J大D、大J小D、中J中D。结果一般如下参数组合JD频率最低点恢复时间波形特征小J大D0.580最低点明显接近0.3Hz短约0.3s响应快阻尼好大J小D5.010最低点不明显约0.1Hz长约1s以上振荡多次恢复慢中J中D2.040折中约0.2Hz约0.5s振荡一次收敛适中结论很直观J增大相当于给系统加重抑制了频率变化率但牺牲恢复D增大相当于增加刹车减少振荡但带来稳态频差。两者需要折中。更有价值的做法是把J从0.5到5按步长0.5扫一遍绘制一组rocof和恢复时间曲线这样写报告时能直接说明整定依据比只给一组参数有说服力得多。5. 常见问题与排查技巧实录5.1 代数环与仿真中途卡死Simulink最烦的报错就是代数环。最常见的是VSG功率测量到控制到PWM到主电路之间形成瞬时代数依赖。解决方法是加Memory或Unit Delay打破循环。我会在反馈功率计算后面加一个一阶低通滤波或Sample Hold效果类似。仿真一开始就check代数环不要等跑不过才找。有些模型代数环隐藏得很深比如在积分器初始条件设置里引入了瞬时反馈这种要仔细看每个红色高亮路径。5.2 算法发散先把闭环拆开单测现象是仿真到0.3秒报Input port contains an Inf or NaN或者功率波形直接飞。原因通常是符号反了、J太小、D太小。我踩过一次把无功环符号写反结果电压发散到几千伏。排查顺序先把开环系统跑通再断开VSG闭环输出给一个固定频率和幅值的调制波看主电路是否正常然后闭合并调制。参数方面初值J给大一些比如2D给大一些比如50跑稳了再优化。不要一上来就追求极端参数那样发散后很难分辨是算法问题还是参数问题。5.3 并网瞬间冲击锁相环初始相位对不上并网瞬间断路器闭合如果相位不一致会有很大冲击。Simulink里三相电压源初始相角默认可能是0但PLL输出需要时间来锁定。处理办法让仿真先跑0.2秒让PLL锁定0.2秒后再闭合断路器或者直接把PLL初始相位设成与电压源一致。如果做更严谨的预同步用幅值、频率、相位三个判据满足后合闸。仿真中预同步的表达式可以在断路器控制逻辑里写比如当|U_mag-U_ref|1%|f_pll-f_0|0.01Hz|θ_pll-θ_g|0.01rad后给合闸脉冲。这样做的好处是后续半实物验证时策略可以直接沿用。5.4 SOC越限与直流母线失稳储能SOC越限会导致母线失稳。加入SOC保护当SOC小于20%或大于90%时限制电池电流指令同时把VSG功率参考值限幅到风机和储能可提供的范围。如果支持不了优先保证本地负载可以切换成离网模式。实际仿真里电池容量设置太小也会引发类似问题直接增大电池Ah容量比调控制器省事。我经常遇到有人在SOC保护逻辑里加了滞回却在母线电压参考值上没做修正导致保护动作后母线电压无法恢复。这里建议保护动作同时把DC/DC电流指令平滑过渡不要直接阶跃到零。5.5 波形毛刺太多测量环节先加滤波PWM开关带来的毛刺很常见别急着改控制。测量功率用Mean或RMS模块或一阶低通滤波截止频率设100~200Hz。三相瞬时电压波形看开关谐波的话把LCL高频谐振点调离开关频率的整数倍频带。波形记录用Scope的Sample Time不要低于Ts否则数据量巨大。我用Data Inspector导出数据再在Matlab里画图比截图清晰得多。还有一个经验顺序是先看“滤波后的功率”再看“瞬时电压电流”。如果在滤波后的功率波形上还能看到明显振荡才是控制问题瞬时波形毛刺大多属于测量或开关谐波不必过分纠结。最后说我个人调试这套风储VSG仿真的最大感受真正决定成败的不是风电模型有多精细而是储能不能在暂态中跟住VSG的功率差额。我一开始纠结于风机怎么建模花了不少时间后来把直流母线电压环和电池双向DC/DC调顺再回来调VSG的J和D进度飞快。建议后来者先把能量流打通再研究参数敏感性。仿真最后可以做一轮参数扫描把J、D、k_q全部设成外部变量一键输出一组曲线写报告或答辩时说服力直接拉满。这套流程我后来在好几套微电网仿真里重复使用都跑得挺稳希望也能帮你少走点弯路。