火电一次调频阻尼不足？VSG虚拟同步发电机并联补偿Simulink仿真全解析

做电力系统仿真的朋友应该都有体会，传统火电单机参与一次调频，最怕的不是动作慢，而是阻尼不足带来的频率振荡。很多论文里拿火电调速器和汽轮机模型跑出来，扰动了之后频率曲线要么超调大、要么来回荡好几拍才平稳，物理上就是系统等效阻尼太低。VSG（虚拟同步发电机）这个思路出来之后，行业内很多人把它当成新能源并网的“伪装器”，但真正把它和火电联合、专门解决一次调频阻尼优化问题的模型，目前公开资料里并不算多。

我这段时间正好把一个“传统火电阻尼不足 + VSG并联补阻尼 + 一次调频特性优化”的Simulink仿真平台完整搭了一遍，跑了负荷阶跃、参数扫描、盲调优化几轮，踩了不少坑，也把可发论文的切入点和创新点理清楚了。这篇就按一个可复现、可对比、可发paper的标准，把模型结构、公式推导、参数整定、Simulink搭建步骤和问题排查一次性讲透。不管你是做研究生课题、电网仿真，还是刚入门VSG建模，这份记录应该都能让你少走一个月的弯路。

1. 问题根源：火电一次调频的阻尼从哪来，又为什么不够用

先别急着搭模型。要把“VSG增加火电阻尼”这件事讲明白，先得搞清楚火电自己那点阻尼到底是怎么来的、为什么总是捉襟见肘。这一节我用工程视角把物理过程拆开说，后面仿真分析也要靠这些概念。

1.1 一次调频的基本动作过程

电网频率跌了之后，火电机组要做的事情可以拆成三步。第一步，转速/频率测量环节发现频率偏离额定值；第二步，调速器（Governor）根据调差系数算出一个功率增量指令，典型的下垂特性表达式是：

[ \Delta P_g = -\frac{1}{R}\Delta f ]

其中R就是调差系数，火电机组一般取0.04~0.05，意思是频率偏差1%（0.5Hz），机组要调整20%~25%额定功率。第三步，汽轮机阀门动作，把高压缸、再热器、中低压缸里的蒸汽流量调整上去，机械功率增加，转子加速，频率回归稳态。整个过程在外特性上，就是一个带延时、带惯性的功率-频率闭环。

1.2 “阻尼不足”这顶帽子扣在谁头上

传统火电系统真正的薄弱点，不在调速器，而在汽轮机和电网的机电耦合环节。汽轮机有一个特别要命的东西叫“再热时间常数”，高压缸做功之后蒸汽还要进再热器再加热一次才到中低压缸，这个通道的时间常数一般在6~12秒，比调速器动作时间（0.2~0.5秒）慢一个数量级。结果就是：调速器指令已经下去了，汽轮机功率却像“半睡半醒”一样慢慢爬，系统在频率动态过程中长期处于一种“激励大、响应慢”的欠阻尼状态。

再叠加不同机组之间调速器参数差异、励磁系统动态、负荷电压特性反馈等因素，整个系统的等效阻尼比经常只有0.1~0.2，甚至更低。反映在仿真波形上就是频率超调大、衰减慢，重一点还会出现区间低频振荡。

1.3 阻尼不足表现出的具体症状

仿真里判断阻尼不足，主要看三个特征。一是频率偏差曲线出现明显超调，比如负荷突增5%，稳态频差可能是0.15Hz，但动态过程中频率最低点能把0.2Hz都穿过去，额外多跌了0.05Hz；二是恢复过程震荡，频率在两三个周期内有明显的来回摆动；三是功率分配混乱，几台机组在动态过程中你推我搡，有功出力和指令对不上。

这三个症状在火电+VSG联合模型里，正好可以当作“优化前”和“优化后”的对比指标。很多论文把这个过程写成“传统火电阻尼不足导致一次调频动态品质恶化，VSG的虚拟阻尼可以附加阻尼转矩改善系统动态特性”，逻辑就是这么来的。

注意：一次调频是“有差调节”，频率不能完全回到50Hz，这是正常的。如果你仿真里频率长期不回50，不是阻尼问题，是缺少二次调频环节。后面我会专门讲这个坑。

2. VSG补阻尼的思路：不是替换，而是并联增强

搞明白火电缺什么之后，VSG的角色就清晰了。它做的不是替代火电的调速器，而是像给汽车加了一套“电子稳定程序”，在火电机械系统响应慢、阻尼不足的频率区间里，主动提供一个可调的虚拟阻尼转矩。这也是这套模型能出成果的核心逻辑。

2.1 VSG的核心方程

VSG的本质，是用变流器的控制算法去模拟同步发电机的转子运动方程。最简化的标幺形式是：

[ J\frac{d\omega}{dt} = P_m - P_e - D(\omega - \omega_n) ]

其中J是虚拟惯量，D是虚拟阻尼系数，P_m是机械功率参考（由下垂控制器给出），P_e是VSG实际输出的电磁功率。这个方程最有价值的地方在于：惯量和阻尼不再是物理机械决定的固定值，而是两个可独立整定的控制参数。想要惯量大就调大J，想要阻尼强就调大D，二者互不干扰，这在真实火电机组上是不可能做到的。

VSG的有功-频率外环还带着一路下垂控制：

[ P_m = P_{ref} + K_p(f_n - f) ]

K_p和火电的调差系数1/R功能一样，都是按频率偏差调节功率。所以VSG并联进系统后，它既参与了一次调频的有功分配，又通过D项往系统里注入了实实在在的阻尼转矩，这就是“增加火电阻尼”的机制。

2.2 为什么VSG能补上火电缺的那部分阻尼

火力发电的阻尼来自转子本身的机械摩擦、涡轮叶片阻尼、调速器动态附加的阻尼，实际数值很难精确控制。VSG的D项则是直接把“频率偏差比例项”作为功率信号注入电流环，电气上响应速度是毫秒级，比汽轮机快了一两个数量级。等效到系统层面，相当于给火电转子运动方程额外并联了一个“快速阻尼绕组”。

这个思路在经典电机学里能找到影子——同步电机的阻尼绕组。只不过VSG用控制算法把这个阻尼做成了可调参数，而且不依赖任何物理绕组。所以在仿真里你会看到很有意思的现象：火电参数完全不变，只把VSG并联进去，系统频率振荡的衰减速度明显加快，超调量下降，调节时间缩短。这就是论文里那句“VSG对传统火电一次调频具有阻尼补偿作用”的最直观证据。

2.3 两种方案对比：替换式VSG vs 火电+VSG联合

我在搭模型之前也犹豫过，是直接建一个纯VSG并网模型，还是做“火电+VSG”联合系统。后来对比下来，差别非常明显。

纯VSG模型适合研究变流器自身控制策略，但对“传统火电”这个关键词几乎没有贡献，论文投出去容易被审稿人质疑“没有贴合电网实际”。火电+VSG联合模型则有两条线可以深挖：一条是VSG参数（J、D）如何影响火电主导的调频动态，另一条是火电调速器与VSG下垂控制器之间的参数协调。这两条线都是既有研究比较薄弱但实际需求很强的方向，模型稀缺性也高，这正是标题里“珍稀模型，可发paper”的含义。

2.4 可发paper的切入点

基于这个联合模型，我梳理了几个适合快速出文章的方向。第一个是“阻尼补偿系数整定”，用解析法或智能算法对D和J进行优化，给出不同火电渗透率下的推荐参数表；第二个是“多台VSG与火电的阻尼分配策略”，研究VSG出力比例不同时系统阻尼比的变化规律；第三个是“频率二次跌落抑制”，传统VSG在频率恢复阶段容易造成功率反调，结合火电的慢动态做协调控制，是审稿人比较买账的创新点。

提示：如果是硕士课题起步，建议选第一个方向，先把单机+VSG的阻尼补偿机理做透，再往多机扩展。模型复杂度和论文工作量都更可控。

3. Simulink建模与实操：从单机到联调的完整流程

这一节是整个项目的重头戏。我用的版本是Matlab R2022a，Simulink + Simscape Electrical。下面按模型拓扑、火电模块、VSG模块、负荷阶跃四步拆开讲，每一步都给出模块选择和参数设置。你不用完全复刻我的数值，但结构建议照搬。

3.1 模型整体拓扑

整个模型是一条“火电- VSG - 无穷大电网”的单线图，中间通过一个升压变压器和线路阻抗连接到理想电压源。负荷扰动用Step模块直接加在机端母线或者无穷大系统母线侧，观测点放在机端频率和联络线功率上。顶层模型建议按功能拆子系统，别全堆在一个窗口里，后面改参数和查错都方便：

• FirePlant_Sub：火电机组完整模型（调速器+汽轮机+同步机）
• VSG_Sub：虚拟同步发电机（控制+变流器+LC滤波）
• Grid_Sub：变压器+线路+无穷大电源+负荷开关
• Measure_Sub：PLL测频、功率计算、频率偏差计算

3.2 火电机组模型怎么搭

火电部分没有必要求助Simscape里那些复杂的详细模型，用传递函数就能把一次调频的动态特性表现得足够好，仿真速度还快。我用的是经典“Governor + Turbine + 同步机转子”结构。调速器传递函数取一阶惯性：

[ G_g(s) = \frac{1}{1 + T_g s} ]

Tg取0.2秒。汽轮机如果忽略再热，就是一阶惯性，但咱们要研究的就是再热时间常数带来的阻尼不足，所以必须把再热环节建出来。我用二阶模型：

[ G_t(s) = \frac{1 + F_{HP} T_{RH} s}{(1 + T_{CH} s)(1 + T_{RH} s)} ]

参数取Tch=0.2s（汽容时间常数）、Trh=7s（再热时间常数）、Fhp=0.3（高压缸功率占比）。这三个数非常关键，Trh越小汽轮机响应越快，阻尼特性就越好，论文里做“火电阻尼不足”场景时，可以把Trh拉到7~10s，效果很明显。

同步机部分不需要建电磁暂态级模型，用二阶转子运动方程即可。如果你想要更贴近并网波形，也可以用Simscape Electrical里的SM模型，但计算量会大很多，跑一次5%负荷阶跃可能要几分钟。我建议先把传递函数模型跑通，整体逻辑验证没问题之后再升级到详细电磁模型。

3.3 VSG模型怎么搭

VSG的Simulink实现方式有两条路线：一是用Simscape Electrical的封装变流器模块，二是在普通Simulink子系统里自己写控制方程，输出电流源接到电网。我强烈推荐后者，原因有三：模型轻、调试直观、参数完全透明。具体做法是在Subsystem里搭三层控制：

第一层是有功-频率下垂，输入频率偏差，输出机械功率参考值Pm，比例系数Kp由调差系数换算；第二层是转子运动方程，把Pm减去输出功率Pe，再减去阻尼项D*(ω-ωn)，经过1/J积分得到角速度，再积分得到相位；第三层是电流指令生成，根据有功功率指令和电压幅值生成dq轴电流参考值，接一个三相受控电流源。

转子运动方程这一块，我直接用了一个MATLAB Function模块，代码很简短：

function [w, theta] = vsg_swing(Pm, Pe, D, J, wn, dt) persistent dw theta_old; if isempty(dw) dw = 0; theta_old = 0; end ddw = (Pm - Pe - D * dw) / J; dw = dw + ddw * dt; w = wn + dw; theta = theta_old + w * dt; theta_old = theta; end

这里直观体现了VSG“惯量+阻尼”的双重效果：J决定频率变化的快慢，D决定振荡衰减的快慢。你在Scope里看到的VSG输出频率曲线，本质上就是这个方程的解。

3.4 负荷扰动与仿真参数配置

负荷阶跃是验证一次调频特性的标准动作。我在母线上接了一个可变负荷，初始负荷700MW，在t=10s时突增50MW（相对基准容量1000MVA是0.05pu），仿真时长60s，足够观察完整动态过程。求解器如果跑传递函数模型，用ode45固定步长0.001s就够了；如果加了详细开关模型，建议改ode15s，不然仿真速度和收敛性都会遭罪。

我这里给你一组我最终采用的参数表，覆盖火电、VSG和系统三部分，可直接作为初值：

4. 参数整定：J和D到底怎么选，才叫“优化”

模型能跑起来之后，第二个核心问题就是参数怎么定。VSG的J和D不是拍脑袋给的，它们直接决定补偿效果。这一节我把我实际用过的整定流程和结果对比写清楚，包括工程扫参法和后续可用的智能优化思路。

4.1 D的整定：从阻尼比反推

VSG向系统注入的阻尼效果，可以折算成系统总阻尼比的变化。对于简化二阶系统，阻尼比近似公式是：

[ \xi \approx \frac{D}{2\sqrt{KJ}} ]

其中K是同步功率系数，标幺化后通常接近1。也就是说，阻尼比和D成正比、和根号J成反比。如果目标阻尼比是0.5，J标幺化后取0.2（对应H≈4s），那么D取80~100左右就能达到理想区间。这个值不是越大约好，D太大会让VSG输出功率对频率偏差过于敏感，和火电抢功率，甚至会激发电流环振荡。

我按D=40/80/160三档做了一轮仿真，负荷阶跃0.05pu，结果如下：

可以看到，D从40加到80，频率最低点和调节时间改善明显；从80加到160，继续改善但幅度变小，而VSG输出功率振荡幅度开始增大，说明已经接近上限。所以对于我这套模型，D取80~120最合理。

4.2 J的整定：惯性时间常数与频率变化率

J的作用是“托住”频率变化率RoCoF。J越大，频率落得越慢，但对阻尼振荡没有帮助，还会让恢复过程变得拖沓。折算成惯性时间常数来看，火电本身的H大概在4~6s，VSG取2~8s都能和系统匹配。我实测下来，J太小时频率曲线会出现高频抖动，J太大时频率最低点之后的爬升变得懒洋洋，18s都回不到稳态频差附近。最终我固定J对应的H=4s做基准，调整D作为主要优化变量。

4.3 优化方法：粒子群整定思路

手工扫参能找到可行解，但要说“优化”，审稿人通常会希望看到有目标函数、有搜索算法。我的做法是定义ITAE指标：

[ ITAE = \int_0^T t \cdot |\Delta f(t)| dt ]

这个指标越小，频率偏差收敛得越早越快。然后用粒子群优化（PSO）搜索J和D的最优组合。粒子群参数我用的是种群20个、迭代30次、惯性权重0.8降到0.2，上下界分别设J∈[1,8]s、D∈[30,200]。跑下来的最优组合在H=3.5s、D=105附近，和手工整定的结果基本一致。用PSO的意义在于：你可以画出目标函数随J/D变化的收敛曲线，放进论文里非常加分。

建议：初学阶段不要一上来就智能优化。先用扫参法把J-D平面的敏感度摸清楚，再上算法，否则优化结果出来你都不知道它合不合理。

4.4 不同参数组的仿真对比表格

论文里常规的“参数敏感性”表格，建议按下面这个思路做：固定D、变化J三档，固定J、变化D三档，再加一组“最优组”和“无VSG组”作对照。每组记录频率最低点、稳态频差、调节时间、RoCoF最大值四个指标。这五组数据一摆出来，整个“VSG改善一次调频特性”的结论就非常扎实了。我建议仿真数据全部用同一随机种子、同一求解器配置，确保对比可信。

5. 仿真结果分析与特征指标提取

参数整定完之后，就要把仿真结果整理成论文里能用的图表和结论了。这一节告诉你怎么读波形、提取哪些指标、怎么设计对照工况才能说服审稿人。

5.1 典型仿真波形怎么读

看频率偏差曲线，重点关注五个时间节点：扰动发生前稳态值、扰动瞬间的初始斜率（RoCoF）、频率最低点、恢复过程中的超调、最终稳态频差。无VSG组的特点是初始斜率不大但振荡明显，末端拖尾很长；加了VSG之后最直观的变化是振荡幅度减小、进入稳态的时间缩短。

看功率波形时，注意火电和VSG的分工：火电承担主要的稳态功率增量（因为容量大），VSG在动态早期承担较大部分增量（因为响应快），随后逐步退出、把功率让给火电。这个“快出慢退”的过程就是VSG补偿火电慢动态的微观机制，也是论文分析的重点段落。

5.2 需要提取的特征指标

论文尽量不要只贴波形，评审看的是量化指标。我整理了一个指标清单：

• 最大频率偏差（Frequency Nadir），单位Hz
• 稳态频率偏差，单位Hz
• 最大频率变化率（RoCoF），单位Hz/s
• 调节时间（频率进入±2%误差带的时间）
• 频率振荡次数
• VSG动态有功支撑比例（动态峰值/稳态增量）
• 系统等效阻尼比估计

前五个指标直接从Scope或者导出的数据矩阵里用MATLAB脚本计算，后两个需要做一点简单的数据处理。阻尼比可以用对数衰减法估计：从振荡波形里取相邻两个波峰值，用公式：

[ \xi = \frac{-\ln(\frac{A_2}{A_1})}{\sqrt{(2\pi)^2 + \ln^2(\frac{A_2}{A_1})}} ]

这个数值能直接反映“火电阻尼不足”和“VSG补阻尼”前后的系统动态品质变化，是整篇论文最有说服力的一列数据。

5.3 边界工况与鲁棒性验证

做科研不能只看一个工况。我建议再补三组边界测试：一是负荷扰动幅值加大到0.1pu，验证大扰动下VSG阻尼是否仍然有效；二是再热时间常数从7s改成10s，模拟不同火电机组老化/类型差异；三是VSG容量从100MW改成50MW和150MW，测试VSG接入比例的鲁棒性。这几组工况跑完之后，你可以得到一张“VSG容量-D最优值”的关系表，这种数据在投稿时特别受欢迎，因为审稿人自己都可以拿去复现。

6. 常见问题与排查技巧实录

这一节是一路踩坑踩出来的，优先级最高，建议先保存。Simulink里模型跑不出来，大部分坑其实都不在控制算法本身，而在仿真配置和环节交互上。

6.1 仿真发散/不收敛怎么办

最常见的现象是波形直接飞上天，数值冲到了1e10。这种情况十有八九是代数环问题：VSG输出功率和电流源之间形成了瞬时反馈，求解器每个步长都要迭代却总不收敛。解决办法是给VSG的有功测量通道加一个低通滤波器，时间常数取0.01~0.02s，切断代数环。另外，检查一下受控电流源的相位是否用了解算得到的theta，如果theta直接反哺到功率计算里，也可能形成正反馈。

还有一招很好用：把变步长求解器切换成固定步长。我通常设0.001s，两万步以内跑完60s仿真，速度完全可以接受。固定步长下模型不容易发散，找问题也容易得多。

6.2 频率一直在漂移、不回到50Hz

这是新手最容易误判的一个“问题”。如果只有一次调频，稳态频差本来就存在。我用0.04调差系数的火电单独承担50MW扰动时，稳态频差大约0.2Hz；加入VSG后，由于VSG调差系数0.03更小，两个调差系数并联等效值变小，稳态频差可以压到0.12Hz左右。这个变化是合理的，不是Bug。如果你的目标是频率无差，那就得叠加二次调频，也就是在VSG或火电调速器前加积分环节，但这属于AGC的范畴，一次调频论文里一般不掺和。

如果频率没有稳定在某个定值，而是一直爬坡，那就要检查负荷是否用常数功率模型，以及输出功率测量是否被低通滤波器严重滞后。我之前就遇到过RC滤波时间常数取得太大（0.5s），导致功率反馈量永远跟不上调度值，系统就一直在“慢爬”。

6.3 VSG与火电之间出现功率分配不均

VSG加入之后，动态过程中VSG功率峰值可能冲到额定容量的两倍。这个在控制逻辑上是正常的，因为J和D产生的“瞬时分担”就是靠高功率输出实现的。但如果稳态时VSG还在抢火电的功率份额，就要检查下垂系数Kp和调差系数R是否匹配。两个电源并联，有功分配按各自调差系数反比，也就是：

[ \frac{\Delta P_{VSG}}{\Delta P_{thermal}} = \frac{R_{thermal}}{R_{VSG}} ]

如果两者的调差系数差得太多，稳态功率分配就会失衡。我建议火电R取0.04，VSG的Kp换算之后的等效R取0.03~0.05之间，两者保持同一数量级，功率分配才合理。

6.4 常见问题速查表

7. 论文框架建议与后续扩展

模型跑通了，参数整定了，波形也出来了，最后一步就是把成果整理成能投稿的文章。我把这套工作安排进了一篇论文的骨架里，供你参考。

7.1 论文结构怎么组织

摘要部分直接点明“传统火电因再热时间常数等因素导致一次调频阻尼不足，本文提出火电- VSG联合系统，通过虚拟阻尼补偿改善频率动态特性”，摘要里要给出量化结果，比如“频率超调量降低40%，调节时间缩短35%”。引言部分按“火电调频现状→阻尼不足问题→VSG技术优势→本文工作”四段式推进，参考文献多引近三年的VSG和一次调频相关文章，审稿人很在意这点。

模型部分重点写两件事：一是火电的再热式汽轮机模型和调速器模型，给出传递函数；二是VSG的转子运动方程和下垂控制，突出J和D两个可调参数。仿真部分用第4节的参数表和第5节的指标表格，对比“无VSG”“VSG固定参数”“VSG优化参数”三组曲线，再给一组参数敏感性和鲁棒性分析。创新点不用多，2~3个就够，其中一个必须是“火电阻尼补偿”这个原创性较强的结论。

7.2 后续扩展的几个实战方向

这套模型可以直接往三个方向扩展。一是从单机扩展到多机系统，研究VSG接入位置和容量分布对全网一次调频阻尼的影响，这已经是电网规划级的课题了；二是把VSG的固定D升级为自适应阻尼，根据频率偏差和变化率实时调整D值，这个方向工作量不大、创新点明确，适合快速出文章；三是把火电锅炉动态考虑进去，建立“锅炉-汽轮机-调速器-VSG”完整的源网协调模型，研究分钟级频率稳定问题，这个深度比默认的一次调频模型高一个档次。

我在实际跑模型的时候，最大的一个体会是：VSG的阻尼补偿效果并不像很多人想的那样“参数越大越好”，它和火电的再热时间常数、调差系数强耦合，存在一个合理区间。找到这个区间，并把它从机理上解释清楚，就是论文的核心贡献所在。所以不要急着堆算法，先把物理过程用仿真数据讲圆，审稿人那一关自然就过了。

最后再分享一个小技巧：整套模型跑完之后，把所有Scope数据用export到MATLAB工作区，写一个自动计算ITAE、频率最低点、阻尼比的脚本，参数扫描的时候一键出表。我当时把整个扫参过程自动化之后，一天能跑完三四十组参数，这个效率在写论文的时候特别顶用。如果你在搭建过程中遇到我上面没提到的问题，也欢迎按这个思路继续往下排——先把离散和代数环清干净，再谈参数优化，这个顺序永远不会错。