搞风储联合并网仿真最绕不开的一个词就是虚拟同步发电机VSG。我在Simulink里把整套“风电储能VSG控制”的系统完整跑通了这段经历值得好好记录一下。网上聊VSG的帖子很多但大多数只给控制框图不给参数或者给了模型参数但一跑就发散根本没有复现价值。这篇文章就从原理、系统结构、模型搭建到调试排错把这条完整链路讲透给后面做新能源并网仿真的同学当个参考。这套仿真解决的核心问题是风力发电本身没有转动惯量大规模接入电网后会拉低系统惯量水平导致频率和电压支撑能力下降。通过储能配合风机再用VSG控制让整个风储单元对外呈现“像一台同步发电机”的特性就能在并网逆变器里模拟出惯性和阻尼。适合电力电子、新能源并网方向的研究生、工程师还有刚开始接触VSG控制的入门者你可以直接照着我这套思路搭建自己的模型。1. 虚拟同步发电机到底在解决什么问题1.1 风电机组惯量不足的短板传统火力发电、水力发电机组转子本身就是个巨大的能量存储器。当系统负荷突变或者某台机组跳闸时同步发电机的转子会依靠自身机械惯性在短时间内释放或吸收动能阻止频率快速跌落或上升。这个“抗扰”能力就是电网稳定运行的基础也是电力系统分析里常说的转动惯量。风电的情况完全不同。直驱型风机PMSG通过全功率变流器并网转速和电网频率没有物理耦合双馈型风机DFIG虽然定子直接连电网但转子侧变流器也把转子转速和电网频率“解耦”了。这意味着风机对外表现出的等效惯量非常小几乎可以忽略。当风电占比超过一定水平后电网里“转不动”的和“不认频率”的机组越来越多一次调频能力下降频率变化率ROCOF明显增大严重时会触碰低频减载动作值。储能装置可以补救这个问题但储能本身是直流设备要通过逆变器才能和电网交互。如果逆变器只是简单地把功率充放出去那它本质上还是个“跟屁虫”需要依靠电网的相位信息才能工作电网出问题它也帮不上忙。要让储能风机真正具备同步发电机的特性就得从控制算法层面做文章。1.2 VSG控制的核心思想与数学模型虚拟同步发电机的思路很直接——在逆变器的控制算法里内嵌一条同步发电机的转子运动方程让逆变器的输出特性在外部看起来和一台同步发电机等效。同步发电机最关键的数学模型是转子运动方程J * dω/dt Pm - Pe - D * (ω - ω0)其中 J 是转动惯量Pm 是机械功率Pe 是电磁功率D 是阻尼系数ω 是角速度。VSG控制通过检测逆变器输出功率 Pe和有功功率指令 Pm对应同步机的原动机输入做差经过运动方程计算出虚拟角速度 ω再对 ω 积分得到虚拟转子角度 θ最后用这个 θ 作为PWM调制波的相位。这么做的好处是当频率跌落时因为 J 的存在输出频率不会瞬间掉下去而是有一个“缓慢下降”的过渡过程。这个过渡过程就是惯量支撑。D 则提供阻尼决定频率到达新稳态所需的调节时间。只要调好 J 和 D逆变器的频率响应就和同步发电机非常接近。无功—电压外特性同样可以模拟。同步发电机的励磁系统有自动电压调节器AVR维持机端电压在给定值附近。VSG中通过无功功率下垂环节生成虚拟励磁电动势 E再根据电压偏差修正E就能复现同步机的无功调压特性。1.3 风储联合VSG的必要性为什么一定要加储能因为VSG提供惯量支撑时需要实实在在的有功功率变化。当电网频率跌落VSG会通过运动方程自动增加输出有功这多出来的功率从哪来风机本体拿不出来——风机的功率由风速决定你没法让风速瞬间变大。唯一的来源就是储能储能系统在此时瞬间放电补齐惯量支撑所需的有功。所以“风储VSG”本质上是让风机提供能量来源、储能提供功率缓冲、逆变器提供并网接口三者在控制层协同工作。没有储能风机单独做VSG只能做个“空壳”响应一旦要求它真的输出惯量支撑功率直流母线电压就会剧烈跌落系统直接保护停机。这套仿真的价值就在于它把风机模型、储能模型、VSG控制算法、并网接口放在同一个环境下联合运行你可以直观看到风波动时、负荷突变时、并网浪涌时各模块之间的功率流动和暂态响应。2. 风储并网系统整体结构与方案选型2.1 系统拓扑选型与能量流我搭建的系统拓扑如下永磁直驱风机PMSG通过机侧整流器和Boost升压电路输出到直流母线蓄电池储能通过双向DC/DC变换器也接到直流母线直流母线再经过三相两电平逆变器和LCL滤波器升压并网。交流侧通过升压变压器接入10kV电网。这个拓扑最大的优点是能量可以双向流动。风速大时风机输出功率大于并网指令储能吸收多余能量风速低时储能放电补足并网功率需要惯量支撑时储能快速释放功率。直流母线电压由储能变换器稳定在1100V风机侧Boost只负责最大功率跟踪逆变器侧按照VSG算法输出PWM波三个控制目标解耦清晰便于分模块调试。和“风机直接通过背靠背变流器并网”的方案比多了一个储能支路和公共直流母线成本更高但灵活性也更强。对于光伏储能系统这个结构同样适用只需要把风机侧模型换成光伏阵列即可储能双向DC/DC和VSG逆变器部分完全复用。2.2 储能系统建模细节储能电池选择锂离子电池等效电路模型国内仿真主流选择用MATLAB/Simulink自带的Battery库配置Nominal voltage 750VRated capacity 200Ah。重点说一下双向DC/DC控制这是个容易出错的地方。储能侧变换器采用电压外环电流内环双闭环控制。外环采集直流母线电压与1100V给定量做差输入PI控制器输出为电池充放电电流参考值内环PI控制器根据电流偏差输出占空比。电池放电时变换器工作在Boost模式把750V提升到1100V充电时工作在Buck模式把1100V降到750V给电池充电。这里有一个关键点电压外环PI参数不能太激进。储能系统的响应速度决定了整个风储系统的动态性能如果外环PI太快直流母线电压容易振荡逆变器功率输出也会跟着抖动太慢则无法快速支撑直流母线电压。我实测整定为Kp2.5、Ki30时动态性能较理想母线电压波动能控制在±3%以内。2.3 并网逆变器容量与LCL滤波参数匹配逆变器容量按风场容量的1.5倍配置选3MW。风机风机额定2MW储能最大放电功率1MW两者叠加后最大功率约3MW。适配这个容量直流母线电压1100V交流侧线电压有效值690VPWM调制比最高约0.92留有约8%的电压裕量。LCL滤波器参数是这样确定的逆变器侧电感取0.4mH网侧电感取0.15mH滤波电容取60μF这样谐振频率约950Hz远高于基波频率50Hz且低于开关频率5kHz的1/4能有效抑制高频开关纹波同时不会引发低频谐振。阻尼电阻在滤波电容支路串联了0.5Ω用来补偿LCL谐振峰的负向幅值实测效果不错。电容支路串阻尼电阻会带来额外损耗但在仿真阶段这个方法效率最高、代码最简单。实际工程上更推荐用有源阻尼通过对电容电流或电压的反馈算法实现同样的阻尼效果但Simulink仿真阶段先用无源阻尼把系统跑通后续再优化比较合理。3. Simulink建模实操过程3.1 主电路模型搭建先搭风机模型。我用的PMSG是Simulink里内置的永磁同步电机模型极对数为8定子电阻2.875Ωd轴电感8.5mHq轴电感8.5mH永磁磁链1.8Wb。机侧整流器用Simscape Electrical里的Universal BridgeMOSFET二极管选项PWM驱动频率5kHz控制方式采用双闭环磁场定向控制FOC转速外环电流内环。风机气动模型单独写了一个S-Function输入风速、转速输出机械转矩包含最大功率跟踪曲线。这里踩过一个坑Simulink的S-Function对输入变量的数据类型要求严格风速信号在Simulink里默认是double类型而S-Function内部如果声明了single类型会报错所以一开始就统一用double类型写接口省得后面排查数据类型报错。储能电池双向DC/DC部分按照第二节参数搭建MOSFET和二极管选2kV/3kA规格电感取2mH开关频率和风机侧整流器都用5kHz。直流母线电容用两个并联的10000μF电解电容模型等效ESR设为1mΩ这样母线电压在功率突变时能保持相对平稳。LCL滤波器和升压变压器用Simscape元件搭建变压器漏抗设为0.06pu仿真精度会比较稳定。3.2 VSG控制算法实现VSG控制是整个模型的核心我在Simulink里用模块程序混合的方式实现。虚拟转子部分直接用积分模块搭建是可行但调参麻烦的我更推荐用MATLAB Function写运动方程。核心代码如下function [theta, w_vsg, Pm_comp] vsg_core(Pe_ref, Pe_meas, w_grid, J, D, isEnable) % 虚拟同步发电机核心运动方程 % 输入有功功率指令Pe_ref实际输出有功Pe_meas电网角速度w_grid % 输出虚拟转子角度theta虚拟角速度w_vsg补偿功率Pm_comp persistent theta_prev w_vsg_prev; if isempty(theta_prev) theta_prev 0; w_vsg_prev 2*pi*50; end if isEnable % 有功-频率下垂控制生成机械功率补偿 droop (w_vsg_prev - w_grid) * 10000; % 下垂系数10kW/(rad/s) Pm_comp Pe_ref - droop; % 转子运动方程使用欧拉法 dtheta w_vsg_prev; dw (Pm_comp - Pe_meas - D * (w_vsg_prev - w_grid)) / J; w_vsg w_vsg_prev dw * 0.0001; theta theta_prev dtheta * 0.0001; theta_prev theta; w_vsg_prev w_vsg; else theta theta_prev; w_vsg w_vsg_prev; Pm_comp 0; end end注意几个细节采样周期固定在100μs即和PWM开关周期同步这样可以避免离散模块与连续模块之间的混叠问题积分步长过大会导致高频振荡所以算法时间步长和模型主步长保持一致。无功控制部分我用的下垂特性Qref设定为0根据电压幅值偏差修正虚拟电动势EE E0 Kq * (Uref - Umeas)其中Kq取200V/kV即电压每差1kV虚拟电动势修正200V。电压检测用Park变换后d轴分量经低通滤波器平滑滤波时间常数0.02s。3.3 风电出力与储能协调逻辑协调策略直接用功率平衡表实现比我见到的有些论文里写的模糊控制更简单、复现性也更好。整个系统的控制优先级是并网有功指令由VSG算法给出风机侧始终按MPPT最大功率输出储能的功率指令由母线电压偏差决定。具体实现是这样的风机侧Boost变换器采用三大控制环——转速外环、电流内环、电压前馈。内置MPPT查表模块输入转速和风速输出功率参考值。储能DC/DC变换器采集母线电压与1.1kV参考值比较经过PI控制器给出电池功率指令。逆变器的有功功率参考值由VSG核心计算由虚拟同步机特性自动调节输出。用Simulink的Bus Selector和Bus Creator把所有功率指令汇总到一个信号总线里监控时直接用Scope查总线信号。需要注意的是风机和储能的启动时序。仿真开始前0.2s电网已经稳定风机逆变器闭锁储能DC/DC预充电到母线电压0.2秒解锁风机Boost开始MPPT0.5秒解锁逆变器且VSG使能端切换为1到1秒完成并网。这个时序设计非常重要同时解锁会让各模块初始电流突变数值算法容易发散。3.4 求解器与仿真步长配置求解器选择对仿真能否收敛影响很大这是我的体会最深的一处。VSG控制里包含连续积分环节和电力电子开关器件Simulink里默认的变步长ode45虽然能跑但速度很慢而且这种混合系统在开关跳变时刻经常出现误差超限。我最后选的是ode23tb变步长求解器相对容差1e-4最大步长1e-5s。如果非要定步长仿真推荐用ode3Bogacki-Shampine算法步长设5微秒。这样每个PWM开关周期内会计算200个点对于5kHz开关频率和50Hz基波的混叠问题处理比较好。仿真时间设置启动时序需要1秒加上稳态观察0.5秒所以仿真结束时间设2.5秒足够覆盖整个动态过程。注意初始步长不要设太大我遇到过因为初始步长过大而立即发散的情况设为max步长的1/10比较稳妥。4. 常见问题与排查技巧实录4.1 模型发散问题的系统性解决发散是仿真过程中最常见的问题几乎每个做VSG仿真的人都会遇到。我整理了四个最典型的原因以及对应的排除方法第一个是代数环Algebraic Loop问题。当我用Stateflow写VSG控制又希望输入输出在同一个仿真步内闭合时Simulink会提示检测到代数环求解器在每一步里都要迭代解方程某一步迭代不收敛就发散。解决方案是在反馈回路上加一个Unit Delay模块把反馈信号延迟一个采样周期。虽然引入了微小延迟但等效于实际控制器的执行延迟对动态特性影响不大。第二个是PI参数积分饱和。虚拟同步机启动瞬间功率误差很大PI控制器积分项迅速饱和导致输出PWM占空比长时间钳位在限幅值逆变器输出电流波形明显畸变。解决办法是给PI控制器后端增加限幅模块Saturation并采用反计算抗饱和结构——当输出达到限幅值时将实际输出和计算输出的差值反馈回积分输入端自动削弱积分量。第三个是滤波器参数与开关频率不匹配引起的谐波发散。LCL滤波器谐振频率如果与开关频率及其整数倍靠得太近谐波增益极高反馈信号产生振荡功率计算值抖动很大。检查方法是看频谱把输出电流接入Powergui FFT Analysis工具如果特定频率分量幅值异常突出就是滤波参数问题。我最后的配套参数上文已给出基本规避了这个问题。第四个是母线电容ESR过小理想电容导致的数值震荡。母线电压在开关切换瞬间理论上有一个没法通过的尖峰数值求解器为了追逐这个尖峰会不断缩小步长导致仿真速度极慢甚至直接中断。给母线电容并联ESR电阻后问题立刻缓解这就是仿真模型里必须加寄生参数的原因。4.2 频率响应异常与参数关联有段时间我调出来的VSG频率响应特别奇怪电网频率跌落0.1Hz后VSG输出主动功率先冲高再回落但没有惯量支撑的“缓降”过程直接快速恢复到指令值更像传统的PQ控制模式。查了代码发现是J值设置过小导致的。J值为5时运动方程中dw (Pm-Pe)/J的加速度值特别大虚拟角速度很快就会跟到新平衡点根本来不及体现惯量特性。把J值调到85单位kg·m²按系统功率基准3MW折算约14秒等效惯量时间常数后频率下降的过渡过程明显拉长顶部出现一个明显的“功率支撑平台”这就是我要的效果。阻尼系数D影响的是频率到达新的稳态值的时间。D值太小系统振荡次数多D值太大响应迟钝。经过实测比较D在25左右时系统响应速度和振荡抑制效果刚好平衡。需要注意的是惯量时间常数和阻尼系数的标幺化换算容易出错建议将所有变量都折算到系统的额定功率和额定角速度基准上可视化调节J、D时更直观。4.3 仿真运行速度慢整系统模型大约包含800多个状态变量用变步长求解器跑2.5秒仿真实际耗时约40分钟这对参数调试来说效率太低。我总结出了几个有效的提速方法一是尽量用离散模块替代连续模块。VSG运动方程用MATLAB Function按固定步长实现而不是用1/s连续积分模块。这能让Simulink把大块连续状态子系统的求解复杂度降下来。二是减少Scope显示的数据点采样率。高采样率Scope会拖慢仿真速度不监控的信号不要用Scope接用To Workspace模块以200Hz采样率录制关键波形后续用MATLAB脚本绘图。三是把不关心的子系统设置为“无操作”或直接用饱和模块替代。比如风机的详细PMSG模型在验证VSG响应阶段不是重点用理想电流源替代就行跑通了再换回详细模型。四是善用Simulink Accelerator模式。把代码编译成C代码执行提速效果在复杂模型上非常明显我个人实测是3到5倍。4.4 并网瞬间电流冲击问题第一次并列成功时系统会有明显冲击电网电流跃变幅值比较大还伴有高频振荡。最有效的处理方式是限制VSG的初始角度和电网角度差。正常锁相环发出的theta是从0开始的但VSG核心输出的theta在模型初始化阶段也是0并网瞬间如果两侧相位不一致就会产生冲击。我在VSG使能开关切换前加了一段预同步逻辑检测虚拟角度与电网角度的差值如果超过一定阈值就用PI控制器把这个差值先调回零附近再开逆变器这样并网冲击大幅减小。这个步骤在单台机仿真时可能不明显但在多机并联时会直接决定系统是否成功并联。5. 仿真结果分析与调参心得5.1 风功率波动场景下的响应仿真里风速设定在8m/s到12m/s之间随机波动叠加了湍流风模型风机的机侧整流器通过MPPT算法跟踪最大功率输出功率在0.9MW到1.8MW之间波动。储能双向DC/DC将波动部分的功率吸收并平滑直流母线电压稳定时间约0.3s电压偏差不超过3%。从波形上看逆变器出口的有功功率基本稳定在1.2MW设定值附近功率波动显著降低。储能电池的充放电电流曲线反映出蓄放功率方向随着风速变化反转在Simulink的Scope里能看到完整的转移过程。5.2 频率跌落时VSG的惯量支撑表现在2秒时设置一个接地短路故障模拟电网发生故障导致频率跌落整个风储系统的频率变化曲线呈现出明显的“先缓后急”特征。故障发生后0.5秒内VSG输出的有功功率从1.2MW上升到1.5MW支撑电网频率下滑由于储能放电能力限制功率峰值维持在1.5MW左右故障结束时VTG在转子运动方程作用下缓慢恢复整个支撑过程约0.8秒。对比相同条件下“无VSG控制”的逆变器频率跌落幅度从0.24Hz缩小到0.11Hz频率变化率也明显降低。这个对比结果直接验证了VSG的惯量支撑效果也说明了为什么新能源并网导则对具备惯量支持的机组有明确的需求。5.3 调参过程中的核心经验最后分享几个我在多次仿真调试中总结出来的关键经验J和D的整定顺序应该遵循“先D后J”。阻尼系数决定系统的稳定性先把阻尼调到临界值附近确保系统不振荡然后逐渐增大J值来增强惯量支撑效果。如果J过大或者D过小系统在扰动后会持续振荡这种情况下先降J再增D不要同时调两个参数。储能系统的功率响应速度必须快于VSG的功率指令变化速率。储能DC/DC的电流内环带宽只要足够高至少要支持500Hz以上的指令跟踪才能保证功率指令下发后电池侧能够快速响应。否则VSG发出支撑功率指令储能响应跟不上直流母线电压就会崩这个现象在仿真前期非常容易遇到。三相对称电压锁相环节带宽可以稍微放宽。实际PLL带宽取值在10Hz到20Hz之间比较合适带宽过低时锁相环跟踪变慢VSG的初始角度预同步时间会明显拉长带宽过高时噪声滤除能力下降角度波动会被放大。最后提醒一下仿真只是验证控制策略的第一步最终还是要靠半实物仿真平台如RT-LAB或者dSPACE跑实时代码才真正接近工程实际。不过能把Simulink模型调通等于完成了从理论到算法验证的关键一步后面的路就好走了。
