1. 从“同步机缺位”说起为什么风电控制要引入VSG我最初接触虚拟同步机是在做微电网项目时当时被一个反直觉的现象震住了一台永磁直驱风电仿真模型在并网点用传统PQ控制跑得好好的可一旦电网频率波动或者进入弱网工况输出功率就开始抖波形和变形虫似的。后来才搞清楚问题不在风电本体而在控制策略——传统PQ控制本质上是一个“跟随型”控制它假设电网是无限大刚性电网风电机组只需要老老实实把功率送出去就行。但在风电渗透率越来越高的今天这个假设越来越不成立。电网频率的支撑、惯量的响应以前靠火电、水电这些同步发电机的转子机械惯量来扛现在换成电力电子换流器换流器本身没有任何物理旋转质量电网频率一波动它只能干瞪眼。这就是业内常说的“同步机缺位”。虚拟同步机VSGVirtual Synchronous Generator核心思路就一句话在换流器的控制算法里用数学方程“造”一台虚拟的同步发电机。你给它一个虚拟转动惯量它就有惯性给它一个虚拟励磁它就能调电压给它下垂特性它就能参与一次调频。控制程序替硬件补上了物理上不存在的转子。这个模型用Matlab/Simulink来做是目前最主流的验证路径。前期我搭过基于双馈风机的VSG也搭过光储VSG但永磁直驱风电配上VSG我认为是最值得吃透的组合——因为永磁直驱机组是全功率变流拓扑机侧和网侧完全解耦网侧换流器控制自由度很高非常适合承载VSG控制算法。如果这篇文章能给正在做风电并网控制仿真的朋友一些参考少走两三个弯路那这篇就值了。2. 永磁直驱风电仿真模型的基础架构风机、电机与换流器2.1 风力机模型Cp曲线怎么设MPPT怎么追踪很多新手上来就关心VSG控制怎么写结果把风力机模型当成“一个转速源”糊弄过去后面直流母线电压、功率波形全跟着错。风力机部分是整个模型的“能源源头”必须认真建。Simulink里建风力机模型常规做法是用标幺值方法。风速v、风轮转速ωr、桨距角β三个输入经过气动模型算出机械转矩Tm。气动模型的核心是叶尖速比λ和风能利用系数Cp公式λ ωr * R / vR是风轮半径Cp是λ和β的非线性函数常见经验公式为Cp 0.5176 * (116/λi - 0.4β - 5) * exp(-21/λi) 0.0068λ其中1/λi 1/(λ 0.08*β) - 0.035/(β^3 1)。这个公式来自Matlab官方风力机demo直接抄就行。搭建时注意把Cp算出来之后和λ、β存进查表模块方便后面观察最大功率点。MPPT最大功率追踪我用的是最优叶尖速比法已知最优叶尖速比λ_opt和当前风速v通过ωr* λ_opt * v / R获得参考转速再交给机侧换流器去跟踪。这个法子实现简单在风速变化缓慢的仿真场景下非常好用。阵风或湍流风工况下有误差但做VSG研究时重点是网侧控制策略风况扰动不宜太复杂否则干扰变量太多。2.2 PMSG与机侧换流器的建模要点永磁同步电机PMSG在Simulink里可以直接用库里的Permanent Magnet Synchronous Machine模块但我个人建议改用dq坐标方程自己写嵌入函数尤其是做VSG这种偏控制方向的研究。原因在于库里的电机模型电气参数齐全但你很难看到磁链、反电动势的内部关系自己用电压方程建模反而能加深对整个功率链路的理解。PMSG在d-q旋转坐标系下的定子电压方程如下ud Rsid Lddid/dt - ωeLqiquq Rsiq Lqdiq/dt ωe*(Ld*id ψf)其中ψf是永磁磁链ωe是电角速度。电磁转矩Te 1.5p[ψf*iq (Ld-Lq)idiq]。机侧换流器用id0控制即可这是表贴式永磁电机最常见的控制策略。转矩和iq成正比控制简单配合MPPT给的转速参考外环转速PI输出iq参考内环电流PI输出ud、uq再做坐标变换、SVPWM完事。经验提醒机侧换流器仿真中的直流侧电容不要省略它的作用类似于“功率缓冲池”。机侧送来多少功率、网侧送走多少功率不平衡时全靠电容电压的升降来体现。电容值取小了直流母线纹波大取大了动态响应慢。一般在仿真里先按经验公式C P_N / (2ω_gV_dc*ΔV_dc)估算再微调。2.3 直流母线与网侧换流器的角色分工整个永磁直驱风电模型可以简化成三段风轮PMSG产生功率→ 机侧换流器整流MPPT控制→ 直流母线 → 网侧换流器逆变VSG控制→ LCL/LC滤波 → 电网。机侧只管一件事让风机工作在最大功率点。网侧只管另一件事把直流功率平滑地变成交流功率送上网同时用VSG算法主动参与电网频率/电压支撑。这两者的解耦核心就是直流母线电压。机侧转速环维持MPPT直流电压的稳定任务交给网侧VSG。在VSG控制里有功功率环本身就有维持直流电压稳定的能力——因为直流功率和交流有功功率平衡了直流电压自然稳定。这也是VSG相对PQ控制的一个天然优势PQ控制下必须额外加一个直流电压外环但VSG的有功-频率环可以通过频率的“暂态摆动”来缓冲不平衡功率直流电压的波动反而更小。这一点我在仿真里对比过多次后面讲波形结论时再展开。3. 虚拟同步机控制的Matlab/Simulink实现细节3.1 有功-频率控制环路转动惯量与阻尼的实现这是整个VSG的灵魂所在。同步发电机的转子运动方程是J * dω/dt Tm - Te - D * (ω - ω0)翻译成电力电子语言就是Tm对应机械输入功率即风电机组输出的直流功率Pm/ωTe对应电磁输出功率Pe/ωD(ω-ω0)对应阻尼项工程实现时更多直接用功率形式Pm - Pe - D * (ω - ω0) J * ω0 * dω/dtSimulink里搭这个环路的顺序是采样网侧瞬时有功功率Pe用瞬时功率公式P udid uqiq在dq轴下计算减去参考有功功率Pref对应MPPT的出力得到不平衡功率ΔPΔP减去阻尼项D*(ω-ω0)之后乘上1/(J*ω0)再积分就得到虚拟角频率ω。到这里有一个关键操作积分器初值必须设置为ω0工频50Hz对应的角频率314.15rad/s否则并网瞬间频率从0爬升冲击电流大到直接触发保护。再往下虚拟频率ω积分得到虚拟相位θ这个θ就是后面坐标变换和SVPWM用的相位角。注意这里完全不需要锁相环PLL——VSG自己生成了相位这是和PQ控制的本质差异。电网扰动时VSG的相位由转子运动方程决定能主动响应频率变化而不是被动跟踪。3.2 无功-电压控制环路与虚拟励磁调节同步发电机另一大特性是励磁调节无功输出影响机端电压机端电压偏差反过来调节励磁。工程化到VSG里就是一条下垂控制虚拟励磁调节E E0 ΔE E0 (Qref - Qe) * Kq kq * (Vref - Vm)解释一下E是虚拟电动势幅值E0是空载电动势Qref、Qe是参考无功和实时无功Kq是无功-电压下垂系数Vref是额定电压参考Vm是并网点电压幅值。有些文献把后面这项称为“电压调节器”实际上就是模拟AVR自动电压励磁调节器的静态特性。下垂系数Kq怎么定有一个工程惯例无功功率从0变到额定无功Q_N时电压允许偏差为额定电压的2%~5%所以Kq ≈ ΔV / ΔQ 0.03 * V_N / Q_N。Simulink建模时就按这个数量级取值然后用试错法微调。值得留意的是VSG的无功环和电网之间是存在耦合的你输出的无功多了并网点电压升高电压升高了无功又会下降。这个负反馈是稳定性的来源。仿真中一旦发现无功持续发散或者振荡优先检查Kq是不是取得太大而不是去调PI参数。3.3 虚拟阻抗和电压电流双闭环怎么接VSG输出的E和θ是合成电压参考的基本要素但E和θ直接作为调制波给SVPWM是不行的——必须经过电压电流双闭环整形同时引入虚拟阻抗来改善阻尼特性。思路如下VSG内部的电压参考源是E∠θ先经过虚拟阻抗得到电压环参考值Vref_dq E∠θ - (Rv jωLv) * I_dq其中I_dq是滤波电感电流Rv、Lv就是虚拟电阻和虚拟电感。虚拟阻抗的作用相当于在换流器和电网之间加了一个“看不见的阻抗”可以抑制换流器与电网之间的功率振荡。弱网环境下这个值非常关键我一般取Rv在0.05~0.2Ω模拟阻尼绕组Lv取0.5~3mH数量级上与滤波电感相当不宜过大否则电压跌落很明显。电压环输出作为电流环参考再经过电流环PI得到调制电压进SVPWM。两个PI的传递函数和参数整定方法下一节专门讲。整个控制链路在Simulink里搭出来后应该是非常清晰的一条主线功率环算幅值相位 → 虚拟阻抗整形 → 电压环 → 电流环 → 调制。4. 参数整定的实际经验从衰减振荡到正常波形4.1 转动惯量J和阻尼系数D的匹配原则刚接触VSG的人最容易在这里翻车J和D是凭空设的参数不知道怎么给初值。给得太小VSG和普通PQ控制没区别惯量支撑形同虚设给得太大功率响应慢甚至出现低频振荡。我常用的初值计算方法是把虚拟同步机等效成一个二阶系统来看。忽略阻尼时的自然振荡角频率大约为ωn sqrt(Pe0 * ω0 / (J * ωn_pll))简化处理时常用ωn ≈ sqrt(PN * ω0 / J)其中Pe0是额定功率ω0是工频角频率。工程上要求振荡频率在0.5~2Hz区间太靠近电网工频会造成谐振太小则动态响应太慢。反解J的范围J Pe0 * ω0 / ωn²以额定功率2MW的风机为例ω0314.15rad/s希望振荡频率在1Hz左右ωn6.28rad/s那么J的数量级是2e6 * 314 / 39.4 ≈ 1.6e7 kg·m²。这里算出来的数值看起来很大很正常因为你的功率基准是兆瓦级的。对应在标幺值体系里J*实际上就只有几秒的时间常数。阻尼系数D的整定逻辑是阻尼比ζ D/(2√(Jω0/ωn))工业上喜欢ζ在0.7~1之间兼顾响应速度和超调量。反推D 2ζ√(Jω0PN/ω0) ≈ 2ζ√(J*PN)。按上面J算出来D在1e5这个量级。仿真时的调参口诀我记下来了发现功率响应超调大、振荡多先增大D发现响应太慢、像个迟钝的胖子先减小J如果两者都不行再检查虚拟阻抗是不是太小、电压环带宽是不是不够。这个顺序基本上覆盖我遇到过的90%的振荡问题。4.2 电压环、电流环PI参数从哪开始调双闭环的整定套路在风电控制里已经非常成熟和普通并网逆变器几乎一样。原则是内环快、外环慢电流环带宽是电压环的5~10倍电压环带宽是功率环的5倍以上。电流环整定方法工程上常用“零点对消法”假设L为滤波电感、R为寄生电阻电流环PI参数取Kp_i L * ω_ciω_ci是期望的电流环带宽Ki_i R * ω_ci期望带宽一般取1~2kHz。仿真模型里滤波电感L0.5mHR0.01Ω取ω_ci2π*1500≈9400rad/s得到Kp_i≈4.7Ki_i≈94。这个作为Initial Guess基本能跑出像样的电流波形。电压环带宽取电流环的1/10左右150Hz左右按类似方法推。开始时可以先给一组保守值Kp_v1、Ki_v50看电压波形再微调。核心观察对象是电压环输出的d轴调制电压是否平稳——如果高频抖动明显多半是内环太快或外环相位裕度不足。4.3 预同步与并网切换的处理VSG要并网必须先解决一个“相位打架”问题VSG自己生成相位θ电网也有自己的相位θg。如果两者不一致就合闸会产生巨大的冲击电流。解决办法是加一个预同步环节并网前检测VSG端电压与电网电压的幅值差、频率差、相位差通过一个缓慢的PI调节器去微调VSG的ω和E直到差值收敛到零附近再发并网指令。Simulink里实现预同步有两种方式一是用一个Boolean信号控制“并网接触器”模块导通前一直让预同步环工作二是把预同步的输出叠加在VSG的ω和E上等误差小于阈值后切掉。我推荐第二种更平滑。一个我踩过的具体坑预同步环的PI参数不能快。我把Kp调到0.5觉得“响应快一点不好吗”结果预同步过程中相位来回穿越始终收敛不了。后来改成Kp0.02、Ki0.1相位差缓慢平滑地逼近零合闸冲击电流被压到了额定电流的5%以内。预同步这个环节本来就是“慢工出细活”别求快。5. 仿真调试中的五个常见坑5.1 启动瞬间冲击电流仿真模型启动时直流母线电容初始电压为零VSG控制尚未稳定很容易在第一个周期就打出骇人的冲击电流。解决这个问题的标准操作分三步先让机侧换流器工作不给网侧发脉冲等直流母线电压建立到额定值比如1100V网侧VSG控制启动但先给一个较小的功率参考值比如0.1pu让电流缓慢爬升并网预同步完成后再切换到大功率参考。很多人忽略第三步直接Pref一步给到1pu结果功率环和电流环互相“抢戏”波形至少振荡5到6个周期才稳定下来。正确的做法是用Ramp信号斜坡信号把Pref从0平滑抬升到期望值斜坡上升时间取0.2~0.5秒仿真里的冲击电流能压得很漂亮。5.2 PLL在VSG中的角色和隐患前面提到VSG不需要PLL来生成功率环的相位但真实模型中你依然绕不开PLL——因为采样回来的三相电压电流需要经过坐标变换才能进dq轴计算而变换角度本身要来自于某个基准。VSG的主相位θ来自转子运动方程但这个θ对应的是换流器内部虚拟电动势的相位不是电网电压的相位。在并网状态下要获取电网电压的dq分量你仍然需要一个PLL或等效的电压定向算法。这就带来一个隐患如果PLL用不好反而会“污染”VSG。我之前遇到过PLL带宽设太高电网电压稍有畸变PLL输出的角度抖动叠加在VSG相位上后电流波形出现了明显的毛刺。后来把PLL带宽降低自然频率10Hz左右、阻尼比0.7问题迎刃而解。经验是VSG体系里PLL是辅助角色负责测量不负责主导能慢则慢。5.3 仿真步长与求解器设置Matlab/Simulink里跑这种高频电力电子模型求解器设置直接影响仿真速度和波形精度。我踩过最狠的坑是为了仿真速度快把步长设成5e-5秒SVPWM的调制效果完全失真波形丑得没法跟人说“这是VSG”。后来统一改成离散求解器discrete步长选2e-6秒或更小功率电路和控制电路全部用离散模块不要混用连续和离散模块否则仿真速度极慢还会报代数环问题脉宽调制频率10kHz采样频率与步长匹配。如果模型太大、仿真时间太长可以把步长放宽到5e-6秒但一定要先对比几个关键波形电流、功率、直流电压确认误差可接受再大量跑工况。做参数扫描扫J、扫D时可以用parfor配合Simulink的batchsim那效率提升立竿见影。5.4 功率波动与直流母线电压的耦合VSG在惯量支撑过程中会短时吸取或释放功率这个功率差额全部体现在直流母线电压的变化上。直流电容若太小电网频率跌落时VSG还没来得及“出力”直流电压就先崩了。这实际上是物理极限不是控制能完全解决的——大扰动工况下必须依赖储能或者直流侧电容的容量。相应地仿真中观察直流电压的跌落幅度能反过来判断虚拟惯量是否给得合适如果惯量支撑瞬间直流电压跌落超过10%说明这个J已经超过当前直流电容能承受的物理极限光调J和D也挽回不了必须增大电容或加储能。这算是我从“波形对不齐”到理解系统本质约束的一个顿悟时刻分享给各位。5.5 离散化误差在SVPWM中的放大SVPWM本质上需要一个精确的电压矢量作用时间计算离散化步长过大时Ton、Toff的计算精度下降低速或轻载工况下会直接体现为电流谐波增大。在永磁直驱风电场景里风速低时风机出力小、电流很小谐波占比就显得很高波形一团糟。解决建议SVPWM模块自己写函数不要用离散化较粗的库模块同时在计算占空比时用当前步长作为时间基准而不是用固定运行的载波周期。另外可以输出电流的THD指标定期观察Matlab的powergui里有FFT分析工具确保THD在5%以内否则控制参数再漂亮也会被谐波拖累。6. 结果怎么读典型波形与工况验证6.1 有功阶跃虚拟惯量的平滑效应模型搭好并稳定运行后第一个必做工况是有功参考阶跃比如风速突变导致MPPT输出的Pref从0.5pu阶跃到0.8pu观察并网有功功率、VSG虚拟频率和电流波形。传统PQ控制的有功阶跃响应是一条“硬直线”功率立即跳到新值几乎无惯性但电流冲击大、直流电压会有一个明显的凹陷VSG响应则是平滑过渡有功功率以45度斜坡状爬升虚拟频率有一个微小的瞬时跌落再恢复——这正是惯性作用的可视化表现。惯量越大爬坡越缓频率暂态偏差越小但进入稳态的时间越长。这个波形的对比一定要做因为它同时也是你验证J和D参数正确性的直接手段。6.2 频率扰动一次调频响应与惯量支撑第二个关键工况是电网频率扰动在电网侧用可编程电压源将频率从50Hz跌落到49.8Hz-0.4%持续1秒后恢复。记录VSG输出有功功率和虚拟频率的变化。正常结果应该是频率跌落的瞬间VSG因惯性效应输出一个短暂的功率突增这就是“惯量支撑”随后因为有功-频率下垂特性功率稳定在一个高于扰动的值这就是“一次调频”。扰动结束后有功功率应平滑回落。如果惯量支撑不明显说明J偏小如果调频稳态功率偏差不符合下垂计算的期望值检查Pe的计算是否正确尤其要注意dq坐标变换角度是否对齐。这一步做完VSG的核心价值就看得清清楚楚了。6.3 VSG与PQ控制的电网适应性对比最后一个建议做的验证把VSG模型里的控制切换到传统PQ控制保留同一套主电路参数在同样的弱网条件电网阻抗较大、短路比SCR2甚至更低下对比两者表现。PQ控制在这种工况下大概率会出现功率振荡甚至失稳因为它的PLL在弱网下动态变差而VSG因为有虚拟惯量和虚拟阻抗耐受性明显更好。这张对比图放在论文里或者技术汇报里分量完全不一样——它把“惯量支撑”从抽象概念变成了肉眼可见的仿真曲线。一点额外的体会实际做下来VSG仿真最难的不是算法本身而是要让整套系统在“像样”的动态里稳定运行。每一次参数改动波形为什么变、物理上对应什么机制这比最终那张漂亮的波形图更有价值。永磁直驱风电这个平台确实很适合做VSG研究——全功率变流结构让控制灵活度大MPPT的出力波动反而成了检验VSG动态性能的天然扰动源。参数和模型的具体细节如果你也在调这套系统欢迎拿你的波形出来交流参数整套盲调那种苦跑过仿真的人都懂。
永磁直驱风电虚拟同步机VSG控制：Matlab/Simulink仿真实现与参数整定
1. 从“同步机缺位”说起为什么风电控制要引入VSG我最初接触虚拟同步机是在做微电网项目时当时被一个反直觉的现象震住了一台永磁直驱风电仿真模型在并网点用传统PQ控制跑得好好的可一旦电网频率波动或者进入弱网工况输出功率就开始抖波形和变形虫似的。后来才搞清楚问题不在风电本体而在控制策略——传统PQ控制本质上是一个“跟随型”控制它假设电网是无限大刚性电网风电机组只需要老老实实把功率送出去就行。但在风电渗透率越来越高的今天这个假设越来越不成立。电网频率的支撑、惯量的响应以前靠火电、水电这些同步发电机的转子机械惯量来扛现在换成电力电子换流器换流器本身没有任何物理旋转质量电网频率一波动它只能干瞪眼。这就是业内常说的“同步机缺位”。虚拟同步机VSGVirtual Synchronous Generator核心思路就一句话在换流器的控制算法里用数学方程“造”一台虚拟的同步发电机。你给它一个虚拟转动惯量它就有惯性给它一个虚拟励磁它就能调电压给它下垂特性它就能参与一次调频。控制程序替硬件补上了物理上不存在的转子。这个模型用Matlab/Simulink来做是目前最主流的验证路径。前期我搭过基于双馈风机的VSG也搭过光储VSG但永磁直驱风电配上VSG我认为是最值得吃透的组合——因为永磁直驱机组是全功率变流拓扑机侧和网侧完全解耦网侧换流器控制自由度很高非常适合承载VSG控制算法。如果这篇文章能给正在做风电并网控制仿真的朋友一些参考少走两三个弯路那这篇就值了。2. 永磁直驱风电仿真模型的基础架构风机、电机与换流器2.1 风力机模型Cp曲线怎么设MPPT怎么追踪很多新手上来就关心VSG控制怎么写结果把风力机模型当成“一个转速源”糊弄过去后面直流母线电压、功率波形全跟着错。风力机部分是整个模型的“能源源头”必须认真建。Simulink里建风力机模型常规做法是用标幺值方法。风速v、风轮转速ωr、桨距角β三个输入经过气动模型算出机械转矩Tm。气动模型的核心是叶尖速比λ和风能利用系数Cp公式λ ωr * R / vR是风轮半径Cp是λ和β的非线性函数常见经验公式为Cp 0.5176 * (116/λi - 0.4β - 5) * exp(-21/λi) 0.0068λ其中1/λi 1/(λ 0.08*β) - 0.035/(β^3 1)。这个公式来自Matlab官方风力机demo直接抄就行。搭建时注意把Cp算出来之后和λ、β存进查表模块方便后面观察最大功率点。MPPT最大功率追踪我用的是最优叶尖速比法已知最优叶尖速比λ_opt和当前风速v通过ωr* λ_opt * v / R获得参考转速再交给机侧换流器去跟踪。这个法子实现简单在风速变化缓慢的仿真场景下非常好用。阵风或湍流风工况下有误差但做VSG研究时重点是网侧控制策略风况扰动不宜太复杂否则干扰变量太多。2.2 PMSG与机侧换流器的建模要点永磁同步电机PMSG在Simulink里可以直接用库里的Permanent Magnet Synchronous Machine模块但我个人建议改用dq坐标方程自己写嵌入函数尤其是做VSG这种偏控制方向的研究。原因在于库里的电机模型电气参数齐全但你很难看到磁链、反电动势的内部关系自己用电压方程建模反而能加深对整个功率链路的理解。PMSG在d-q旋转坐标系下的定子电压方程如下ud Rsid Lddid/dt - ωeLqiquq Rsiq Lqdiq/dt ωe*(Ld*id ψf)其中ψf是永磁磁链ωe是电角速度。电磁转矩Te 1.5p[ψf*iq (Ld-Lq)idiq]。机侧换流器用id0控制即可这是表贴式永磁电机最常见的控制策略。转矩和iq成正比控制简单配合MPPT给的转速参考外环转速PI输出iq参考内环电流PI输出ud、uq再做坐标变换、SVPWM完事。经验提醒机侧换流器仿真中的直流侧电容不要省略它的作用类似于“功率缓冲池”。机侧送来多少功率、网侧送走多少功率不平衡时全靠电容电压的升降来体现。电容值取小了直流母线纹波大取大了动态响应慢。一般在仿真里先按经验公式C P_N / (2ω_gV_dc*ΔV_dc)估算再微调。2.3 直流母线与网侧换流器的角色分工整个永磁直驱风电模型可以简化成三段风轮PMSG产生功率→ 机侧换流器整流MPPT控制→ 直流母线 → 网侧换流器逆变VSG控制→ LCL/LC滤波 → 电网。机侧只管一件事让风机工作在最大功率点。网侧只管另一件事把直流功率平滑地变成交流功率送上网同时用VSG算法主动参与电网频率/电压支撑。这两者的解耦核心就是直流母线电压。机侧转速环维持MPPT直流电压的稳定任务交给网侧VSG。在VSG控制里有功功率环本身就有维持直流电压稳定的能力——因为直流功率和交流有功功率平衡了直流电压自然稳定。这也是VSG相对PQ控制的一个天然优势PQ控制下必须额外加一个直流电压外环但VSG的有功-频率环可以通过频率的“暂态摆动”来缓冲不平衡功率直流电压的波动反而更小。这一点我在仿真里对比过多次后面讲波形结论时再展开。3. 虚拟同步机控制的Matlab/Simulink实现细节3.1 有功-频率控制环路转动惯量与阻尼的实现这是整个VSG的灵魂所在。同步发电机的转子运动方程是J * dω/dt Tm - Te - D * (ω - ω0)翻译成电力电子语言就是Tm对应机械输入功率即风电机组输出的直流功率Pm/ωTe对应电磁输出功率Pe/ωD(ω-ω0)对应阻尼项工程实现时更多直接用功率形式Pm - Pe - D * (ω - ω0) J * ω0 * dω/dtSimulink里搭这个环路的顺序是采样网侧瞬时有功功率Pe用瞬时功率公式P udid uqiq在dq轴下计算减去参考有功功率Pref对应MPPT的出力得到不平衡功率ΔPΔP减去阻尼项D*(ω-ω0)之后乘上1/(J*ω0)再积分就得到虚拟角频率ω。到这里有一个关键操作积分器初值必须设置为ω0工频50Hz对应的角频率314.15rad/s否则并网瞬间频率从0爬升冲击电流大到直接触发保护。再往下虚拟频率ω积分得到虚拟相位θ这个θ就是后面坐标变换和SVPWM用的相位角。注意这里完全不需要锁相环PLL——VSG自己生成了相位这是和PQ控制的本质差异。电网扰动时VSG的相位由转子运动方程决定能主动响应频率变化而不是被动跟踪。3.2 无功-电压控制环路与虚拟励磁调节同步发电机另一大特性是励磁调节无功输出影响机端电压机端电压偏差反过来调节励磁。工程化到VSG里就是一条下垂控制虚拟励磁调节E E0 ΔE E0 (Qref - Qe) * Kq kq * (Vref - Vm)解释一下E是虚拟电动势幅值E0是空载电动势Qref、Qe是参考无功和实时无功Kq是无功-电压下垂系数Vref是额定电压参考Vm是并网点电压幅值。有些文献把后面这项称为“电压调节器”实际上就是模拟AVR自动电压励磁调节器的静态特性。下垂系数Kq怎么定有一个工程惯例无功功率从0变到额定无功Q_N时电压允许偏差为额定电压的2%~5%所以Kq ≈ ΔV / ΔQ 0.03 * V_N / Q_N。Simulink建模时就按这个数量级取值然后用试错法微调。值得留意的是VSG的无功环和电网之间是存在耦合的你输出的无功多了并网点电压升高电压升高了无功又会下降。这个负反馈是稳定性的来源。仿真中一旦发现无功持续发散或者振荡优先检查Kq是不是取得太大而不是去调PI参数。3.3 虚拟阻抗和电压电流双闭环怎么接VSG输出的E和θ是合成电压参考的基本要素但E和θ直接作为调制波给SVPWM是不行的——必须经过电压电流双闭环整形同时引入虚拟阻抗来改善阻尼特性。思路如下VSG内部的电压参考源是E∠θ先经过虚拟阻抗得到电压环参考值Vref_dq E∠θ - (Rv jωLv) * I_dq其中I_dq是滤波电感电流Rv、Lv就是虚拟电阻和虚拟电感。虚拟阻抗的作用相当于在换流器和电网之间加了一个“看不见的阻抗”可以抑制换流器与电网之间的功率振荡。弱网环境下这个值非常关键我一般取Rv在0.05~0.2Ω模拟阻尼绕组Lv取0.5~3mH数量级上与滤波电感相当不宜过大否则电压跌落很明显。电压环输出作为电流环参考再经过电流环PI得到调制电压进SVPWM。两个PI的传递函数和参数整定方法下一节专门讲。整个控制链路在Simulink里搭出来后应该是非常清晰的一条主线功率环算幅值相位 → 虚拟阻抗整形 → 电压环 → 电流环 → 调制。4. 参数整定的实际经验从衰减振荡到正常波形4.1 转动惯量J和阻尼系数D的匹配原则刚接触VSG的人最容易在这里翻车J和D是凭空设的参数不知道怎么给初值。给得太小VSG和普通PQ控制没区别惯量支撑形同虚设给得太大功率响应慢甚至出现低频振荡。我常用的初值计算方法是把虚拟同步机等效成一个二阶系统来看。忽略阻尼时的自然振荡角频率大约为ωn sqrt(Pe0 * ω0 / (J * ωn_pll))简化处理时常用ωn ≈ sqrt(PN * ω0 / J)其中Pe0是额定功率ω0是工频角频率。工程上要求振荡频率在0.5~2Hz区间太靠近电网工频会造成谐振太小则动态响应太慢。反解J的范围J Pe0 * ω0 / ωn²以额定功率2MW的风机为例ω0314.15rad/s希望振荡频率在1Hz左右ωn6.28rad/s那么J的数量级是2e6 * 314 / 39.4 ≈ 1.6e7 kg·m²。这里算出来的数值看起来很大很正常因为你的功率基准是兆瓦级的。对应在标幺值体系里J*实际上就只有几秒的时间常数。阻尼系数D的整定逻辑是阻尼比ζ D/(2√(Jω0/ωn))工业上喜欢ζ在0.7~1之间兼顾响应速度和超调量。反推D 2ζ√(Jω0PN/ω0) ≈ 2ζ√(J*PN)。按上面J算出来D在1e5这个量级。仿真时的调参口诀我记下来了发现功率响应超调大、振荡多先增大D发现响应太慢、像个迟钝的胖子先减小J如果两者都不行再检查虚拟阻抗是不是太小、电压环带宽是不是不够。这个顺序基本上覆盖我遇到过的90%的振荡问题。4.2 电压环、电流环PI参数从哪开始调双闭环的整定套路在风电控制里已经非常成熟和普通并网逆变器几乎一样。原则是内环快、外环慢电流环带宽是电压环的5~10倍电压环带宽是功率环的5倍以上。电流环整定方法工程上常用“零点对消法”假设L为滤波电感、R为寄生电阻电流环PI参数取Kp_i L * ω_ciω_ci是期望的电流环带宽Ki_i R * ω_ci期望带宽一般取1~2kHz。仿真模型里滤波电感L0.5mHR0.01Ω取ω_ci2π*1500≈9400rad/s得到Kp_i≈4.7Ki_i≈94。这个作为Initial Guess基本能跑出像样的电流波形。电压环带宽取电流环的1/10左右150Hz左右按类似方法推。开始时可以先给一组保守值Kp_v1、Ki_v50看电压波形再微调。核心观察对象是电压环输出的d轴调制电压是否平稳——如果高频抖动明显多半是内环太快或外环相位裕度不足。4.3 预同步与并网切换的处理VSG要并网必须先解决一个“相位打架”问题VSG自己生成相位θ电网也有自己的相位θg。如果两者不一致就合闸会产生巨大的冲击电流。解决办法是加一个预同步环节并网前检测VSG端电压与电网电压的幅值差、频率差、相位差通过一个缓慢的PI调节器去微调VSG的ω和E直到差值收敛到零附近再发并网指令。Simulink里实现预同步有两种方式一是用一个Boolean信号控制“并网接触器”模块导通前一直让预同步环工作二是把预同步的输出叠加在VSG的ω和E上等误差小于阈值后切掉。我推荐第二种更平滑。一个我踩过的具体坑预同步环的PI参数不能快。我把Kp调到0.5觉得“响应快一点不好吗”结果预同步过程中相位来回穿越始终收敛不了。后来改成Kp0.02、Ki0.1相位差缓慢平滑地逼近零合闸冲击电流被压到了额定电流的5%以内。预同步这个环节本来就是“慢工出细活”别求快。5. 仿真调试中的五个常见坑5.1 启动瞬间冲击电流仿真模型启动时直流母线电容初始电压为零VSG控制尚未稳定很容易在第一个周期就打出骇人的冲击电流。解决这个问题的标准操作分三步先让机侧换流器工作不给网侧发脉冲等直流母线电压建立到额定值比如1100V网侧VSG控制启动但先给一个较小的功率参考值比如0.1pu让电流缓慢爬升并网预同步完成后再切换到大功率参考。很多人忽略第三步直接Pref一步给到1pu结果功率环和电流环互相“抢戏”波形至少振荡5到6个周期才稳定下来。正确的做法是用Ramp信号斜坡信号把Pref从0平滑抬升到期望值斜坡上升时间取0.2~0.5秒仿真里的冲击电流能压得很漂亮。5.2 PLL在VSG中的角色和隐患前面提到VSG不需要PLL来生成功率环的相位但真实模型中你依然绕不开PLL——因为采样回来的三相电压电流需要经过坐标变换才能进dq轴计算而变换角度本身要来自于某个基准。VSG的主相位θ来自转子运动方程但这个θ对应的是换流器内部虚拟电动势的相位不是电网电压的相位。在并网状态下要获取电网电压的dq分量你仍然需要一个PLL或等效的电压定向算法。这就带来一个隐患如果PLL用不好反而会“污染”VSG。我之前遇到过PLL带宽设太高电网电压稍有畸变PLL输出的角度抖动叠加在VSG相位上后电流波形出现了明显的毛刺。后来把PLL带宽降低自然频率10Hz左右、阻尼比0.7问题迎刃而解。经验是VSG体系里PLL是辅助角色负责测量不负责主导能慢则慢。5.3 仿真步长与求解器设置Matlab/Simulink里跑这种高频电力电子模型求解器设置直接影响仿真速度和波形精度。我踩过最狠的坑是为了仿真速度快把步长设成5e-5秒SVPWM的调制效果完全失真波形丑得没法跟人说“这是VSG”。后来统一改成离散求解器discrete步长选2e-6秒或更小功率电路和控制电路全部用离散模块不要混用连续和离散模块否则仿真速度极慢还会报代数环问题脉宽调制频率10kHz采样频率与步长匹配。如果模型太大、仿真时间太长可以把步长放宽到5e-6秒但一定要先对比几个关键波形电流、功率、直流电压确认误差可接受再大量跑工况。做参数扫描扫J、扫D时可以用parfor配合Simulink的batchsim那效率提升立竿见影。5.4 功率波动与直流母线电压的耦合VSG在惯量支撑过程中会短时吸取或释放功率这个功率差额全部体现在直流母线电压的变化上。直流电容若太小电网频率跌落时VSG还没来得及“出力”直流电压就先崩了。这实际上是物理极限不是控制能完全解决的——大扰动工况下必须依赖储能或者直流侧电容的容量。相应地仿真中观察直流电压的跌落幅度能反过来判断虚拟惯量是否给得合适如果惯量支撑瞬间直流电压跌落超过10%说明这个J已经超过当前直流电容能承受的物理极限光调J和D也挽回不了必须增大电容或加储能。这算是我从“波形对不齐”到理解系统本质约束的一个顿悟时刻分享给各位。5.5 离散化误差在SVPWM中的放大SVPWM本质上需要一个精确的电压矢量作用时间计算离散化步长过大时Ton、Toff的计算精度下降低速或轻载工况下会直接体现为电流谐波增大。在永磁直驱风电场景里风速低时风机出力小、电流很小谐波占比就显得很高波形一团糟。解决建议SVPWM模块自己写函数不要用离散化较粗的库模块同时在计算占空比时用当前步长作为时间基准而不是用固定运行的载波周期。另外可以输出电流的THD指标定期观察Matlab的powergui里有FFT分析工具确保THD在5%以内否则控制参数再漂亮也会被谐波拖累。6. 结果怎么读典型波形与工况验证6.1 有功阶跃虚拟惯量的平滑效应模型搭好并稳定运行后第一个必做工况是有功参考阶跃比如风速突变导致MPPT输出的Pref从0.5pu阶跃到0.8pu观察并网有功功率、VSG虚拟频率和电流波形。传统PQ控制的有功阶跃响应是一条“硬直线”功率立即跳到新值几乎无惯性但电流冲击大、直流电压会有一个明显的凹陷VSG响应则是平滑过渡有功功率以45度斜坡状爬升虚拟频率有一个微小的瞬时跌落再恢复——这正是惯性作用的可视化表现。惯量越大爬坡越缓频率暂态偏差越小但进入稳态的时间越长。这个波形的对比一定要做因为它同时也是你验证J和D参数正确性的直接手段。6.2 频率扰动一次调频响应与惯量支撑第二个关键工况是电网频率扰动在电网侧用可编程电压源将频率从50Hz跌落到49.8Hz-0.4%持续1秒后恢复。记录VSG输出有功功率和虚拟频率的变化。正常结果应该是频率跌落的瞬间VSG因惯性效应输出一个短暂的功率突增这就是“惯量支撑”随后因为有功-频率下垂特性功率稳定在一个高于扰动的值这就是“一次调频”。扰动结束后有功功率应平滑回落。如果惯量支撑不明显说明J偏小如果调频稳态功率偏差不符合下垂计算的期望值检查Pe的计算是否正确尤其要注意dq坐标变换角度是否对齐。这一步做完VSG的核心价值就看得清清楚楚了。6.3 VSG与PQ控制的电网适应性对比最后一个建议做的验证把VSG模型里的控制切换到传统PQ控制保留同一套主电路参数在同样的弱网条件电网阻抗较大、短路比SCR2甚至更低下对比两者表现。PQ控制在这种工况下大概率会出现功率振荡甚至失稳因为它的PLL在弱网下动态变差而VSG因为有虚拟惯量和虚拟阻抗耐受性明显更好。这张对比图放在论文里或者技术汇报里分量完全不一样——它把“惯量支撑”从抽象概念变成了肉眼可见的仿真曲线。一点额外的体会实际做下来VSG仿真最难的不是算法本身而是要让整套系统在“像样”的动态里稳定运行。每一次参数改动波形为什么变、物理上对应什么机制这比最终那张漂亮的波形图更有价值。永磁直驱风电这个平台确实很适合做VSG研究——全功率变流结构让控制灵活度大MPPT的出力波动反而成了检验VSG动态性能的天然扰动源。参数和模型的具体细节如果你也在调这套系统欢迎拿你的波形出来交流参数整套盲调那种苦跑过仿真的人都懂。
