风储VSG——基于虚拟同步发电机的风储并网系统Simulink仿真这名字听着像又一个毕业论文选题。但真上手做起来你会发现它解决的问题特别具体风电场里一台台变流器越来越多传统火电机组逐渐退出后电网频率跌落时谁站出来提供惯量支撑VSG的核心思路说起来也不玄乎——把同步发电机的转子运动方程写成控制算法塞进变流器的控制回路里让电力电子设备的行为表现得像一台自带飞轮的发电机。这篇文章是我实际搭建这个仿真的完整复盘从系统架构到Simulink里的模块摆放从参数整定到波形判读再到我踩过的几个坑。如果你正在做风储并网、微电网、新能源变流器控制相关的仿真或者只是想把VSG从论文公式变成能跑的模型这篇文章应该能帮你少走不少弯路。1. 为什么风储并网需要虚拟同步发电机从“电网缺惯量”说起1.1 同步发电机凭什么稳频要理解VSG先得理解传统同步发电机在电网里到底干了什么。同步发电机的转子是一大坨旋转的铁芯和绕组转起来以后本身就带着巨大的旋转动能。当电网负荷突然增加、系统频率有下降趋势时转子会自然地减速把自己的一部分动能释放成电磁功率送给电网这个过程是物理机制自动完成的不需要任何控制器下发指令。用转子运动方程来描述就是J*(dω/dt) Tm - Te - D*(ω - ω₀)J是转动惯量Tm是机械转矩Te是电磁转矩D是阻尼系数。这个方程里最关键的是J那一项——有了惯量频率变化率就被抑制住了电网频率不会瞬间跌落。频率从50Hz跌到49.8Hz需要一定的时间这就给一次调频、AGC这些上层控制争取了宝贵的响应时间。1.2 风电和储能天生“缺惯量”风电变流器和光伏逆变器则完全是另一回事。它们本质上是一个受控电流源或受控电压源控制器的响应速度在毫秒级功率输出完全跟着指令走。问题是——它们的功率输出和电网频率之间没有任何物理上的“惯性耦合”。风速突然变小风机输出功率就掉了电网频率可能因此下跌但变流器本身并不知道电网正在面临压力。新能源渗透率越高电网里可用的旋转惯量就越少频率稳定性问题就会暴露出来。储能可以帮忙但储能变流器和风电变流器一样是电力电子设备它也没有天然的旋转惯量。这时候就轮到VSG出场了——让变流器在控制算法层面“模拟”同步发电机的转子行为。你把转子运动方程写进控制环让变流器的输出功率根据电网频率的变化自动作出反应频率降了就多输出一点功率频率升了就少输出一点。从电网侧看这套设备就像一台虚拟的同步发电机惯量、阻尼、调频调压特性一应俱全。风电储能的组合做VSG有个天然优势VSG在做惯量支撑和一次调频时需要短时间内改变输出功率这部分能量从哪来如果只有风电就得让风机减载运行留出备用容量这本身牺牲了发电量。但如果有储能就可以让储能去吸收或释放波动功率风电尽量保持在最大功率跟踪MPPT状态两者在直流母线上汇合由储能侧变流器统一执行VSG逻辑。这也是风储联合系统做惯量支撑最常用的工程方案。2. 风储VSG系统的物理接线与能量脉络直流母线汇集方案2.1 风力机与储能单元的并网位置做仿真之前先把系统框架在脑子里立起来。风储VSG并网系统在结构上一般分三层风机和储能单元作为一次能源层两个变流器作为中间执行层电网作为最终的接入对象。风机部分可以采用永磁同步发电机PMSG配全功率变流器也可以用双馈异步发电机DFIG。从仿真角度我个人建议用PMSG因为它的控制思路更直接机侧变流器负责MPPT转速控制网侧变流器管直流母线电压。风机的机械功率可以通过风速、叶尖速比λ和风能利用系数Cp(λ,β)计算出来在Simulink里可以搭一个简单的查表模型也可以直接用风电机组封装模块。储能单元则是电池加双向DC/DC变换器挂在直流母线上。双向DC/DC的作用有两个一是维持直流母线电压稳定二是按VSG控制器的需求快速吞吐功率。DC/DC的控制一般用双闭环外环是直流母线电压PI控制内环是电感电流控制。当VSG需要多输出功率时直流母线电压会短暂下降DC/DC外环立刻响应让电池释放能量把母线电压拉回来。2.2 直流母线方案的三个关键支路整个系统的能量脉络可以拆成三条支路来看。第一条是风电支路捕获的风能经过机侧整流器送进直流母线在MPPT模式下它只负责尽可能多地把风能转成电能不参与频率调节。第二条是储能支路电池通过双向DC/DC接到直流母线它的角色是“能量缓冲池”风电功率波动大的时候它来吸收VSG调频需要能量的时候它来释放。第三条是网侧变流器支路直流母线上的电功率经过三相全桥逆变器通过L或LCL滤波器送到电网上。VSG控制算法就作用在这条支路上——它直接决定逆变器输出的电压幅值、相位和频率。这里有个容易被忽略的点网侧变流器在传统PQ控制模式下用的是电网电压锁相角做坐标变换但在VSG模式下坐标变换的角度应该来自VSG自己的虚拟转子相位θ。这意味着系统在控制层面拥有了一个“虚拟的旋转坐标系”而不是被动跟随电网。这也是VSG和传统矢量控制之间最本质的控制架构差异。3. 在Simulink里搭VSG核心控制器逐块分解3.1 主电路模块怎么选Simscape还是受控源Simulink里搭建主电路有两种主流做法一种是直接用Simscape Electrical中的三相全桥IGBT模块、LCL滤波器、电网电源和变压器这种模型最贴近物理实际开关管频率、滤波特性都能真实反映另一种是简化受控源模型用受控电压源或受控电流源代替逆变桥这种方式仿真速度快适合先验证控制算法。做完整的风储VSG系统仿真我建议主电路用Simscape Electrical但把控制算法做成独立的控制子系统。开启Simscape环境只需要在命令行敲powerlib然后从Simscape Electrical Specialized Power Systems库里拖元件。逆变桥用Universal Bridge模块选IGBT/Diode作为开关器件把直流母线电压接到输入端输出端接LCL滤波器。这里注意直流母线侧的电源不能直接用理想电压源——因为你需要看到功率在母线之间的动态交换更合理的是用电池模型挂在母线上或者直接用一个带内阻的直流电压源加一个可控电流源来模拟风电/储能的功率注入。滤波参数我给一组可以直接落地的参考值额定功率20kW电网线电压380V/50Hz直流母线电压750VLCL滤波器取L₁1mH、C30μF、L₂0.5mH阻尼电阻Rd1Ω。载波频率取10kHz。这套参数下开关纹波抑制效果比较理想电流谐波能压到5%以内。3.2 VSG转子运动方程的Simulink表达VSG控制器是整个系统的灵魂它在Simulink里的实现并不复杂。有功-频率控制对应的转子运动方程离散化后长这样omega[n] omega[n-1] dt * (P_ref - P_e - D*(omega[n-1] - omega₀)) / (J * omega₀) theta[n] theta[n-1] dt * omega[n]在Simulink里最简单的做法是用一个MATLAB Function模块输入端接P_ref有功参考值、P_e实测输出有功、omega上一时刻角频率、theta上一时刻相位输出端给新的omega和theta。系统仿真用的离散步长如果是1e-5秒这个离散化误差足够小。我喜欢用这种方式而不是积分模块搭因为代码里可以顺手加限幅、滞回和初始化逻辑调试起来一目了然。这个函数的核心就是在做一件事算功率差除以虚拟惯量J积分得到角频率偏差再积分得到相位。两个积分环节叠加就把同步发电机的“转子运动”给模拟出来了。看一眼代码function [theta, omega] vsg_rotor(Pref, Pe, D, J, omega0, dt, omega_pre, theta_pre) domega (Pref - Pe) - D * (omega_pre - omega0); domega domega / (J * omega0); omega omega_pre domega * dt; theta theta_pre omega * dt; if theta 2*pi theta theta - 2*pi; end end无功-电压控制相对简单——无功误差进一个PI控制器输出作为电压幅值修正量叠加上去得到虚拟电势幅值E。有了E和θ三相参考电压就出来了Ea E*sin(theta) Eb E*sin(theta - 2*pi/3) Ec E*sin(theta 2*pi/3)这三路信号就是后续电压电流双闭环的参考值。3.3 电压电流双闭环与PWM的衔接VSG输出的是内电势参考电压但这个参考电压不能直接给PWM——还需要经过电压外环和电流内环的跟踪控制才能让逆变器实际输出逼近期望值。电压环输入是“参考电压-输出电压”的误差通过PI控制器输出电流环参考值电流环输入是“电流参考-电感电流”的误差再经过PI控制器输出调制波。这里有个容易被新手卡住的环节参考电压的三相正弦信号到底遇到什么坐标系控制。我的习惯是在VSG生成的θ坐标系下做dq变换。abc三相参考电压经过Clark-Park变换变成ud_ref、uq_ref电容电压实测值和电感电流实测值也变换到同一坐标系下然后做dq解耦控制。这样PI控制器只需要处理直流量控制带宽好设计。调制波生成之后接PWM Generator模块载波频率写10000输出脉冲连接到Universal Bridge的触发端。整个控制链路就是VSG转子方程A/Q控制器 → 参考电压 → 电压电流双闭环 → PWM → 逆变桥 → LCL滤波 → 并网点。4. 参数整定的工程手感先给初值再看波形调优4.1 转动惯量与阻尼系数的匹配逻辑参数整定是整个仿真里最耗时、也最考验经验的部分。先说转动惯量J。如果严格按物理公式换算设额定容量S20kVA交流侧额定角频率ω₀2π50≈314.16rad/s取惯量时间常数H2s那么J2H*S/ω₀²≈0.81 kg·m²。这个值可以作为仿真初值。但你完全不需要把这个值死记硬背——仿真平台里J的单位和你模型的功率基准一致时直接用J1左右起步基本都能跑出像样的动态曲线。重点是理解J的影响趋势J越大虚拟转子的“质量”越大频率变化率越小但频率回稳速度变慢容易出现低频慢振荡J太小频率动态看上去很“脆”负载稍微一动频率就剧烈抖动。调的时候记住这个方向就不会乱。阻尼系数D的直观作用是提供振荡衰减。D偏小的时候有功功率阶跃响应会出现明显的超调和持续振荡波形像弹簧减震器失效的悬架D偏大系统响应变得迟钝频率偏差恢复得很慢。我的经验是先用D10这样的值起步观察暂态波形然后按超调量来决定加减。惯量和阻尼之间是要配合着调的J和D共同决定了系统阻尼比如果超调大就把D加大如果响应慢就把J减小或者按阻尼比公式反着推。4.2 下垂系数和PI参数的初值参考VSG里还有两组下垂系数有功-频率下垂K_f和无功-电压下垂K_q。K_f的作用是做一次调频频率偏移造成的功率增量ΔPK_f*(f₀-f)。工程上比较直觉的做法是让频率偏差0.5Hz对应满功率变化即K_fP_rated/0.540000W/Hz。K_q同理无功电压下垂可以取额定无功容量对应的电压偏差在5%以内K_qQ_rated/ΔU比如Q_rated10kVar、允许压降20VK_q500Var/V。这两个系数给了初值之后一次调频调压的幅度基本就落在合理区间剩下的是微调。电压电流双闭环的PI参数也有理论路子可走。电流内环的带宽设为载波频率的十分之一左右比如fc10kHz时电流环带宽ω_ci≈1000rad/s此时Kp≈L₁ω_ci1mH10001Ki≈Rω_ciR是电感电阻和开关管压降折算电阻约0.1ΩKi≈100。电压外环带宽比电流环低一个数量级约200rad/sKp≈Cω_cv30μF*2000.006Ki≈(1/Road)*ω_cv。这些初值算出来之后边看波形边微调一般半天就能收敛。5. 仿真工况与波形判定拿什么指标衡量这套系统好不好5.1 工况一并网稳态与额定输出模型搭好、参数给上初值之后我习惯先跑一个最简单的工况系统先并网有功参考给12kW无功参考给0kVar风速恒定在10m/s让仿真跑1秒左右等所有波形进入稳态。这一步的核心是验证控制链路有没有大的原则性错误。稳态要看的波形包括三相并网电压电流是否正弦对称、电流谐波是否被压住、直流母线电压是否稳定在750V附近、电池SOC是否平稳变化。如果电流波形出现了明显的低频振荡或直流偏置八成是坐标变换角度错位。如果直流母线电压出现二倍频波动甚至崩溃多半是逆变桥和DC/DC之间的功率匹配有问题。5.2 工况二负荷突变的频率响应稳态能稳住之后加压测试。典型设计在2秒时让并网点负荷从10kW突增到20kW或者让功率参考值阶跃上升。这个时候就要紧盯VSG的输出有功P_e和虚拟角频率ω。理想响应是功率阶跃瞬间P_e快速跟随ω先有一个小的过渡跌落然后在一两秒内恢复到接近50Hz的稳态值。这一段波形就是VSG的“名片”——没有惯量的传统PQ控制不会有这种动态过渡它会瞬间完成功率调整但也不会向电网提供任何频率支撑信号。你换不同编号的J和D参数重跑这个工况波形差异非常直观这也是给导师展示VSG价值最方便的实验。5.3 工况三风速波动下储能的平抑作用第三个必需工况是风速突变。在3秒时把风速从10m/s阶跃降到8m/s风电输出功率明显下滑。对比两种场景储能不参与控制时电网侧的功率输出会跟着风电一起掉频率出现偏差储能参与VSG协调控制时电池会立刻补上功率缺口网侧输出功率基本维持平稳频率波动压得很小。判定标准就是看网侧变流器的有功输出曲线和电池功率曲线——电池功率应该对风电功率的下降形成镜像补偿。这一步做完整个“风储VSG”才是真正意义上的风储联合而不是风电和储能各跑各的。另外这张波形图也是论文里最“出片”的图之一建议务必截好。6. 实际调试中踩过的坑与绕行建议6.1 仿真步长与求解器选择直接决定成败先讲一个最坑的很多人第一次搭含IGBT高频开关的Simulink模型一仿真就卡死、发散或者报“Simulation stopped because the model contains singularities”。这些问题的根源基本都在步长。10kHz的PWM如果你用固定步长求解器步长至少要到1e-6秒仿真实体时间长1秒就意味着百万步计算抱怨慢之前先看看步长设置。我的惯用组合是Powergui选择“Discrete”采样时间设为1e-5秒求解器用离散固定步长如果用的是Simscape Electrical开关级模型建议改成“Continuous”加变步长ode23tb求解器把最大步长限制在1e-4秒。跑一次20kW系统1秒仿真时间开关级模型大约需要几分钟这正常。如果你想把一次调频的几十秒过程完整跑完时间太长的话可以把功率等级缩小或者控制算法先跑平均模型等逻辑验证完再切回开关模型。6.2 PLL锁相角和VSG相位“打架”的问题我最早把VSG控制接到并网系统上时犯过一个特别隐蔽的错误坐标变换角度我习惯性用了PLL锁相环输出的电网相角结果功率反馈、电流环全乱套了。原因很简单——VSG的本质是自己建立相位方向而PLL锁的是电网相位。如果PLL角度和VSG的虚拟转子相位存在一个不断变化的偏差整个dq坐标系就是错乱的控制性能自然一团糟。正确的做法是在VSG模式下所有坐标变换的角度都用VSG内部生成的θ。PLL在这个架构里只作为并网预同步过程使用或者干脆在并网稳定后不参与坐标变换。这个坑值得单独记下来因为报错不会提示你波形乱了你也很难一下想到是角度源头的问题。6.3 从平均模型过渡到开关模型控制逻辑验证提速最后一个建议和模型调试顺序有关。不要一上来就搭完整的IGBT开关模型——调试一圈下来PC风扇能把人逼疯。我现在的习惯是先搭平均模型逆变桥用受控电压源代替VSG控制器的转子方程和参考电压生成逻辑保持不变先跑通功率控制逻辑把Pi参数、下垂系数都调到满意再替换成Universal BridgePWM Generator的开关级模型回头只验证开关纹波和谐波指标。这样既保证了最终模型的“含金量”也把调试迭代周期压缩了一大截。储能电池的SOC初始值、初始相位、电容初值这些细节也不要忽略。电池模块的SOC没设置会让仿真一步都跑不动电容初始电压不给也会在前几个微秒造成大电流冲击。建议在模型里给所有积分环节加状态初始化省得每次重启都是“从零开始”。做风储VSG这个项目最深的体会是VSG并不是什么高不可攀的算法它就是把一道控制方程——转子运动方程——用工程手段落到变流器上。但越简单的东西越讲究细节坐标变换的角度取自哪里、惯量和阻尼怎么配合、储能响应够不够快、滤波器谐振点压没压住谐振峰。每一个看似小到可以忽略的取舍最后都会在你截图的波形上体现得清清楚楚。按这套架构搭完跑通你对“并网变流器控制”的理解会往上迈一个台阶。