做风储并网仿真这行也有年头了最让我印象深刻的不是某条控制曲线调得多漂亮而是刚接触虚拟同步发电机VSG时的一次真实翻车模型结构跟下垂控制基本一样只是往控制环里多加了两个积分环节结果上电瞬间频率波形抖得像心电图直流母线直接掉到保护值。后来才明白VSG这层“虚拟转子”看着简单实际上把惯量、阻尼、锁相环交互、直流电容能量小时常数这些耦合问题全拉进来了。这篇就以“风储VSG-基于虚拟同步发电机的风储并网系统Simulink仿真”为主线把我在Simulink里从零搭建模型、整定参数、设计工况、排查振荡的完整过程写出来给正要复现这套模型的同行一条能直接落地的路线也把那些教科书里不讲、只有真正跑模型才会踩到的坑交代清楚。1. 为什么风储系统要引入虚拟同步发电机1.1 传统并网变流器缺了什么东西普通的新能源并网逆变器用PQ控制或者下垂控制本质上是把变流器当成一个受控电流源或压控电压源看到的是“功率指令-输出电流”的直接映射。这种控制方式带来的问题在强电网下不明显但在弱电网里非常致命系统一旦缺了一台同步机的转动惯量频率变化率就会变得很大扰动来了没人“撑住”电网。同步发电机为什么能稳不是因为它的调速器反应快而是因为它有庞大的转子转动惯量功率不平衡时转速不会突变这是物理学层面的惯性支撑。变流器没有转子但它可以通过控制算法在逻辑上“虚构”一个转子。这就是虚拟同步发电机名称的由来把同步发电机的二阶运动方程搬到变流器的控制系统里让变流器在外部特性上表现为有惯量、有阻尼、能参与一次调频和一次调压的电压源。1.2 风储VSG系统的分工逻辑风电本身是间歇性电源风机转子虽然有物理转动惯量但那部分动能和电网频率的解耦程度取决于风机变流器的控制策略常规情况下并不能直接给电网提供惯性响应。储能则不同它本质上是能量双向流动的功率缓冲池有能力在短时间尺度上释放或吸收有功功率。风储VSG就是把这两个单元组合起来风电负责提供稳态能量储能负责支撑暂态功率。仿真的重点于是落在两部分一是储能变流器实施的VSG控制是否能在频率跌落时快速释放有功二是风电侧功率波动会不会导致虚拟转子频率偏移超过允许范围。所以这套模型本质上是一个“能量调度惯量响应”的联合验证平台。提示商业软件里已经有现成的VSG模型块但强烈建议自己搭控制环。自己搭一遍才能真正理解J和D参数变化时功角曲线的形态为什么会变。直接拖现成模块出了问题往往束手无策。1.3 VSG和下垂控制之间到底差在哪下垂控制也参与一次调频它和VSG最核心的区别在于响应速度和动态过程。下垂控制的功率变化和频率偏差完全同步公式是P P0 kp(f0 - f)频率一变功率指令立刻跟上没有过渡。VSG则多了一个频率变化率路径等效成转子运动方程频率从正常状态到新稳态之间会有一段“爬坡-过冲-回落”的过程这个过程的形状由惯量J决定收敛速度由阻尼D决定。用通俗说法理解下垂控制是看到变化就动手VSG是看到变化先“转起来”再动手。这个延迟换来的好处是电网频率不会因为大量变流器同时响应而出现更剧烈的高频振荡。2. 风储并网Simulink模型的整体架构2.1 直流母线架构 vs 交流并联架构搭建风储VSG模型之前首先要确定拓扑。仿真里最常见、也最容易收敛的是直流母线架构风力发电机组经整流器接到直流母线储能电池通过双向DC/DC变换器也接到同一条直流母线然后由一台双向并网变流器统一逆变并网。VSG算法就放在这台并网变流器的控制器里。这个架构最大的好处是只需要一套VSG控制环风电和储能的功率耦合直接在直流母线上完成不需要在两个变流器之间做通信配合。建模时可以省掉风机的转速控制、桨距控制等复杂细节把风电场等效成一个“向直流母线注入功率的可控电流源”储能和DC/DC负责维持直流母线电压并网变流器负责VSG控制。交流并联架构则把风电场变流器和储能变流器在交流侧并联各自独立并网。这种拓扑更接近实际工程但两个变流器同时参与VSG控制时存在功率分配问题仿真调试难度成倍增加。新手起步阶段优先选直流母线架构跑通之后再往交流并联扩展。2.2 各模块的角色与接口定义我建议在动手搭模型之前先列一个接口定义表。Simulink模型里模块之间最怕信号名混乱、单位不一致、坐标变换基准不统一。表风储VSG仿真模型模块划分与接口定义模块功能关键接口单位/基准风电机组等效源模拟风电功率输入直流注入电流 I_wA储能电池提供功率缓冲端电压 U_bat内阻 R_batV / Ω双向DC/DC稳定直流母线电压直流母线电压 V_dcV并网变流器执行VSG控制逆变并网三相端电压、输出电流V / ALCL滤波器滤除开关谐波输出至电网的电压电流V / A电网模型提供理想电压源和线路阻抗三相电压源、R/L阻抗V / ΩVSG控制器模拟同步机转子运动方程有功功率、频率、参考电压kW / Hz / V2.3 搭建顺序建议按照我多次搭建和帮人排查的经验从零开始不要按“从风电往并网点搭”的顺序而应该按“电网往直流母线反着搭”的顺序。先搭电网、滤波器、并网变流器主电路用开环正弦波调制跑通整流/逆变基本功能。然后把VSG控制环接进去不接风电和储能直流侧先用一个理想恒压源供电。这时候只验证并网变流器能不能在设定功率下稳定输出。确认无误后再把直流电压源替换成“储能双向DC/DC风电等效电流源”的组合。这样每一步故障范围都非常小出现发散或者振荡第一时间就能定位到具体环节。3. VSG控制器在Simulink中的具体实现3.1 从同步发电机运动方程到控制框图虚拟同步发电机算法的核心是同步电机的转子运动方程J * dω/dt P_ref / ω_n - P_e / ω_n - D * (ω - ω_n)其中J是虚拟惯量单位kg·m²D是虚拟阻尼系数P_e是变流器实测输出的电磁功率ω_n是额定角频率ω是虚拟角频率。这个方程说的是当输出功率P_e不等于参考功率P_ref时虚拟角频率ω就会变化变化快慢由J决定变化过程中的阻尼强弱由D决定。仿真实现上只需要把这个方程用积分器搭出来再对虚拟角频率积分得到虚拟功角θ作为PWM参考波的相位基准。搭建时要注意方程的标幺化处理。如果用实际单位J的量级和D的量级差距很大控制器内部数值容易溢出。建议统一折算到标幺值下运算功率除以额定功率角频率除以额定角频率在模型里所有控制和测量量都用标幺值只在坐标变换和PWM生成环节还原成实际值。3.2 有功-频率环的搭建细节有功-频率环输入是参考功率和实测功率的差值经过运动方程后输出虚拟角频率和相位。这里最容易犯的一个错误是把参考功率定死为恒定值导致频率跟踪能力丢失。正确做法是叠加下垂环节让虚拟调速器能根据频率偏差修正功率指令公式为P_order P_ref kp * (f_n - f_measured)其中kp是有功-频率下垂系数等效于同步发电机的调差系数。仿真里测得的电网频率建议不要直接用PLL的输出去做下垂反馈PLL动态会影响整个VSG的阻尼特性。更稳妥的做法是用VSG自己内部的虚拟角频率ω做一次调频反馈也就是把运动方程里的(ω-ω_n)项同时作为下垂作用的输入。把运动方程的下垂项和阻尼项区分开D * (ω - ω_n)管的是阻尼振荡kp * (f_n - f)管的是稳态调频。两者作用不同D太小振荡kp太小频率静差大。很多模型出现问题根本原因是把这两个参数合二为一失稳后分不清是阻尼不够还是调频系数不合适。3.3 无功-电压环的搭建细节无功-电压环模拟的是同步发电机的励磁调节器经典实现是E E_ref kq * (Q_ref - Q)E是虚拟电动势幅值Q是实测无功功率kq是无功-电压下垂系数。为了提高模型的暂态响应精度推荐在此基础上增加一个一阶惯性环节模拟励磁绕组的时间常数E E_ref kq / (1 T * s) * (Q_ref - Q)T一般取0.01~0.05秒。这个惯性环节的物理意义是励磁电流的调节不可能瞬间完成它在抑制励磁调节与电网之间的高频互动上作用很大。Simulink里用“Transfer Fcn”块实现注意初始值要设置成额定励磁电动势E_ref避免启动瞬间电压出现跳变。3.4 虚拟阻抗与电压电流双闭环VSG输出的虚拟电动势E和相位θ组合起来就是一个三相参考电压源不能直接拿去产生PWM中间还需要经过“虚拟阻抗-电压环-电流环”三层校正。虚拟阻抗的作用是调整VSG的外特性让它的输出阻抗主要表现为感性从而和电网阻抗之间获得合理的功率传输特性。实际模型里可以省略虚拟阻抗计算直接用LC滤波器的滤波电容电压作为电压环反馈电流作为电流环反馈。电压外环用PI控制电流内环用P或PI控制。这里有一个关键经验电流内环的PI参数不能只按滤波器电感设计还需要考虑数字控制一拍延迟。延迟时间等于控制周期Ts如果Ts2e-5秒那电流环主导极点的响应速度不宜超过1/(20*Ts)否则相位裕度不够很容易在高频段激发振荡。3.5 PWM与测量链路Simulink里生成PWM信号最方便的方式是用内置的SVPWM发生器或者自己实现基于载波的SPWM。注意参考电压幅值不能超过直流母线电压除以√3超出会进入过调制区此时输出电压波形严重畸变电流谐波激增VSG的功率计算也会被污染。测量环节建议用Simscape Electrical自带的“Three-Phase V-I Measurement”块测完三相电压电流后经Clark变换到αβ坐标系再用Park变换到dq坐标系。坐标变换的初相位必须和电网电压初始相位对齐否则启动瞬间会有很大的冲击电流。我搭模型时习惯先让PLL跑到稳态再闭合功率环相当于给系统一个“并网预同步”过程避免初始状态不匹配导致的不收敛。4. 参数整定先算后调不会瞎试4.1 虚拟惯量J的计算方法虚拟惯量最直接的整定方法是从同步发电机的惯性时间常数H出发。H的物理含义是在额定功率下转子的动能可以维持额定转速运转多长时间典型同步机H在2到10秒之间。折算关系为J 2 * H * S_B / ω_n²假设仿真系统的基准容量S_B100kVA额定频率50Hzω_n314.16rad/s取H4s则J 2 * 4 * 100000 / (314.16²) ≈ 8.1 kg·m²这个值可以作为初始值。J太小虚拟转子反应太快失去提供惯量的意义J太大扰动后的频率震荡幅度大恢复时间拉长而且直流母线侧的瞬时功率波动也会更为剧烈。一般以仿真波形中频率变化率是否限制在0.5Hz/s以内为调整目标。4.2 阻尼系数D的整定方法阻尼系数D在运动方程里对应的物理量是阻尼转矩从二阶系统角度D影响阻尼比ζζ D / (2 * √(J * K))其中K为同步转矩系数由功角特性决定。实际模型中更直观的整定方法是按频率偏差计算假设允许频率偏差1Hz对应额定频率2%的偏差输出的阻尼功率达到额定功率则D S_B / (2π * Δf) 100000 / (2π * 1) ≈ 15915考虑到功角特性中的同步转矩贡献这个初值往往偏大实际仿真时可以从计算值的一半开始逐步增大。调参时观察频率衰减振荡的衰减包络衰减太慢说明D偏小频率收敛到新稳态太硬、没有惯性过渡感说明D偏大。4.3 电流环与滤波器的参数配合VSG外特性最终要通过电流内环去实现所以控制参数和主电路参数必须联动整定。先确定LC滤波器参数再设计电流环。滤波器谐振频率一般取开关频率的1/10到1/20。如果开关频率10kHzLC谐振频率设在1kHz左右比较合适。公式f_res 1 / (2π * √(L * C))取L0.5mHC20μF则f_res≈1592Hz满足要求。电流环PI参数的整定建议用工程法电流环带宽取开关频率的1/10即1kHz。这一带宽下PI参数可以按典型的“零极点对消期望带宽”方法计算Kp L * ω_bandwidthKi R * ω_bandwidth。比如L0.5mHR0.05Ωω_bandwidth2π*1000≈6283则Kp≈3.14Ki≈314。注意Simulink里的PI控制器有两种形式的表达一种是Kp Ki/s一种是Kp * (1 1/(Ti*s))接线前先确认参数含义否则“对应同一转速误差”时增益差了一个Ki倍数整定结果完全不可比。5. 典型工况设计与结果判断5.1 负荷突增工况验证惯量和阻尼推荐设计一个最简单但最有说明力的工况在t1.0秒时电网负载从50kW突增到80kW。记录以下曲线电网频率曲线看频率跌落速率和最低点。VSG有功输出曲线看储能释放功率的动态过程。虚拟角频率ω曲线看VM内转子速度的摆动。没有VSG的系统负荷突增时频率会陡降变化率大。接入VSG后频率变化率明显放缓最低频率点滞后于负荷突变点之后在阻尼作用下收敛到新的稳态。如果频率曲线出现持续的等幅或发散的振荡首先回查D参数是否过小其次是电流内环带宽是否被滤波器谐振峰的相位滞后“咬住”。5.2 风速阶跃工况验证风储协同第二个常用工况是风速阶跃。比如在t1.5秒风电注入直流母线的功率从60kW降至30kW。此时并网变流器输出的总功率如果仍维持在80kW差额50kW就要由储能系统补偿。观察直流母线电压波动和储能SOC变化。这个工况的关键观察点是直流母线电压在风电功率突降瞬间会有一个突然下坠双向DC/DC的电压环会快速抬升放电电流补偿。如果直流母线电压波动超过额定值的5%比如700V直流母线掉到665V以下说明储能侧DC/DC响应带宽不够需要把电压环PI调快或者增加直流母线电容值。5.3 与下垂控制的对比波形验证VSG价值最直观的方式是设置一个对比组同一套模型把VSG运动方程屏蔽掉直接采用有功-频率下垂控制然后在相同负荷突变工况下对比电网频率变化率。具体操作是把运动方程块替换为ω ω_n - kp * (P_out - P_ref)其他参数完全不变。对比结果会说明问题下垂控制的频率变化率在负荷突变瞬间是跳变的VSG则是斜线平滑过渡。这个对比实验做完基本就能直观理解为什么高比例新能源电网要追求“虚拟惯量”而非单纯的下垂响应。6. 调试中绕不开的坑6.1 启动瞬间发散新手遇到的绝大多数发散问题都出在初始条件。积分器初始值设置不当尤其是虚拟功角θ的初始值没有和电网电压相位对齐闭合瞬间会有巨大的冲击电流Simulink直接报“Nonlinear solver”或者波形飞到天文数字。解决方法先单独跑一个不带动量方程的开环工况读取PLL稳定后的相位值把这个值作为VSG虚拟角度积分器的初始条件。同时把电压环PI输出初始值设定为额定调制波电流环PI初始值设为0。Simulink的Integrator块可以直接设置Initial condition这一步不要偷懒直接关系到能否开机就稳定。6.2 离散步长与代数环模型用离散仿真时步长选择和控制周期保持一致。如果控制周期Ts2e-5秒而仿真步长也设为2e-5秒功耗计算经过滤波块后会有延迟。功率计算输出和运动方程输入之间如果有直接代数回路Simulink会给出“Algebraic Loop”警告严重时仿真速度极慢。解决代数环的标准做法是在功率计算后加一个Unit Delay或者Memory块人为产生一拍延迟。这一拍延迟在连续系统里看起来像是降低了一点点控制性能但在离散系统里它实际上对应了真实的数字控制延迟仿真结果反而更贴近实际工程实现。6.3 直流母线电压失控直流母线电压击穿或跌落常见原因是储能DC/DC初始状态和风电功率不匹配。搭建模型时判断逻辑一定要清晰风电向直流母线注入多少功率并网变流器就要向外输出多少功率。只要风电功率和逆变功率之间存在长期偏差差额一定会反映到直流母线电容上要么充电到过压要么放电到欠压。VSG本身只控制并网变流器“按虚拟转子特性输出功率”并不负责“平衡所有功率”储能系统才是平衡器。调VC控制作用时不要把目标放在让直流母线电压纹波为零应该允许它在一定范围内波动否则储能DC/DC的响应应力很大仿真步长需要压缩到极低算力负担成倍增加。6.4 SVPWM过调制导致电流畸变当系统启动、频率摆动剧烈时VSG输出的参考电压幅值可能短暂超出SVPWM的线性调制范围。过调制区电流谐波瞬间增大控制环看到谐波后进一步加大调节力度就会形成恶性循环。处理方式是在VSG电动势输出后加一个限幅环节限制调制波幅值在0.95倍线性调制范围内。限幅不会破坏VSG的物理特性只会让虚拟电动势在极端工况下暂时“饱和”符合同步发电机励磁限制器的思路。6.5 电网阻抗设置不当导致弱电网振荡许多仿真模型把电网直接设成理想电压源内阻为0。这种情况下VSG即使参数并不太合理也往往能跑出好看的波形一换成弱电网环境就问题暴露。建议从建模第一天就把电网阻抗设置成实际值假设并网点短路比为5基准容量100kVA、线电压380V时电网阻抗大约为Z_g U² / (SCR * S_B) 380² / (5 * 100000) ≈ 0.289Ω按X/R7分配电感和电阻则R≈0.04ΩX≈0.28ΩL≈0.89mH。把这一小段阻抗加到电网电压源前面仿真结果才具备参考价值。6.6 储能SOC到达边界导致功率掉链子长时仿真时还要关注储能SOC边界。如果仿真时间只有10秒SOC几乎不变但如果做分钟级长时仿真SOC会持续下降到边界此时双向DC/DC应该限制放电功率而VSG控制环并不会感知这个限制它依然会按照指令要功率两者冲突会导致模型发散。解决方法是把储能SOC监测信号接入VSG的功率参考值生成环节当SOC低于20%时逐步降低有功功率参考值并且限幅。这种逻辑在后来的实际项目中被证明是必须的提前加上可以省掉大量长时仿真的无谓排查。7. 基于模型进一步扩展的思路模型基本跑通之后可以沿着三个方向扩展。第一是控制策略层面在VSG运动方程中加入虚拟惯量自适应频率变化率大的时候自动增大J频率接近稳态时恢复默认值能明显提升弱电网下的暂态响应。第二是模式切换让储能变流器并网时运行VSG模式离网时自动切换为电压源模式带本地负荷这套逻辑在微电网孤岛场景很有用。第三是多机并联方向两台VSG同时并网时由于虚拟阻抗和PLL耦合会产生环流振荡此时需要引入下垂共享和无功均分逻辑。这三个方向都是当前学术界和工程界的热点但前提都是先把单机模型的基础打牢。我见过不少同行一上来就想做多机协同结果单机的参数都没摸清最后只能在互相矛盾的现象里打转。基础模型里的每一个参数都应该具备物理解释这是做仿真最有价值的部分。最后再分享一条实际的调节经验模型跑通之后不要急着追求复杂的工况先用负荷阶跃测试把D从0逐渐增大观察频率曲线的阻尼包络收敛过程。这个过程做过一次你对J和D的理解会远超看一百篇论文。J决定“振荡的节奏”D决定“振荡的收敛”记住这两个作用任何调节都不会跑偏。