光伏储能虚拟同步发电机并网仿真模型：VSG控制与参数整定

做新能源并网仿真的朋友，应该没少被一个问题折磨过：光伏出力一波动，电网频率就跟着抖；逆变器响应倒是快，可快得过头，电网需要的那点“惯性”它反而给不了。这几年大家都在讨论虚拟同步发电机（VSG），本质上就是让逆变器在控制策略上“装”成一台同步发电机的样子，把转子的惯量和阻尼特性补回来。这套光伏储能虚拟同步发电机并网仿真模型，就是用 Simulink 把这件事完整搭了出来——光伏侧做 MPPT，储能侧稳住母线电压，逆变器侧跑 VSG 控制，最后并入电网。如果你们课题组正好在研究新能源并网惯量支撑、VSG 参数整定，或者做毕设需要一套能改参数、能看波形的完整模型，这篇内容应该能帮你省下不少折腾的时间。

1. 为什么分布式电源需要“模拟”同步电机

先聊一个底层问题：电网过去一百年靠什么稳定运行？答案是同步发电机自带的转子。汽轮机或者水轮机拖着转子转，转子上有实实在在的转动惯量，负载突然增加时，转子转速会平滑地下降一点，这个过程是渐变的、有阻尼的，给了调度系统足够的反应时间。而光伏和储能经过逆变器并网时，直流侧的能量变换是电力电子开关在微秒级完成的，没有任何物理旋转部件，自然也就没有转子惯性。结果就是，微网或者高比例新能源电网一旦功率失配，频率变化率（RoCoF）可以飙得非常高，继电保护装置还没来得及动作，频率就已经越过边界了。

虚拟同步发电机的思路很直接：既然物理上没有转子，那就在控制算法里“虚拟”一个。VSG 控制算法实时观测逆变器的输出功率和电网频率，然后按照同步发电机的转子运动方程去计算一个虚拟转速和虚拟功角，再把这个虚拟转子角度当成逆变器输出电压的相位指令。于是，逆变器呈现出来的外特性就类似于一台同步发电机——它有了等效惯量 J，有了等效阻尼 D，甚至还能模拟同步电机的励磁调节特性。这就是为什么 VSG 被认为是解决高比例新能源并网稳定性问题的热门方案。

很多初学者会问：下垂控制（Droop Control）不也能让逆变器参与调频调压吗？为什么非得用 VSG？我的理解是：下垂控制本质上是静态特性，有功-频率曲线是一条直线，功率一变频率立刻跟着变，缺少了时间尺度上的“过渡过程”。VSG 的下垂特性是嵌入在二阶转子方程里的，频率变化是带惯性的渐变过程，这就在暂态阶段把频率变化率压住了。两者在稳态时的外特性可能比较接近，但暂态响应完全不同。电网最关心的恰恰是暂态过程。这套模型选择 VSG 而不是下垂，就是这个原因。

这套光伏储能 VSG 并网仿真模型的结构分成三层：直流侧是光伏阵列加 Boost 电路做 MPPT，储能电池通过双向 DC/DC 变换器接在直流母线上，负责把母线电压稳住；逆变器侧采用三相桥，经过 LCL 滤波器并入电网或者接入本地负荷；控制层核心是 VSG 算法，负责生成电压幅值和相位指令，再由电压电流双闭环去跟踪。整个链路在 Simulink 里分层搭建，每一层都能单独拆出来观察波形和调参，既适合做原理验证，也适合做参数敏感性分析。

2. VSG 控制算法：核心模块与数学模型拆解

VSG 控制算法是整个模型的心脏。这一节把有功-频率环、无功-电压环、虚拟阻抗、电压电流双闭环逐个拆开讲，每个环节背后的物理意义和控制逻辑都需要理解清楚，因为后续所有参数整定都建立在这些公式之上。

2.1 有功-频率环：把转子运动方程装进控制器

VSG 的核心是二阶转子运动方程：

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

其中，( J ) 是虚拟转动惯量，单位是 ( kg \cdot m^2 )，它决定了对频率变化的惯性响应强度；( D ) 是虚拟阻尼系数，它模拟了同步电机的阻尼绕组效应，作用是抑制振荡；( \omega ) 是 VSG 虚拟角速度，( \omega_g ) 是电网同步角速度；( P_m ) 是虚拟机械功率（由有功参考和调频环节给出），( P_e ) 是逆变器当前实际输出的有功功率。

在 Simulink 里实现这个方程不难：测量逆变器的输出电压和电流，计算出瞬时功率，经过低通滤波得到 ( P_e )；把 ( P_m ) 和 ( P_e ) 的差值乘以 ( 1/J ) 后做积分，得到虚拟角速度 ( \omega )；( \omega ) 再积分一次得到虚拟功角 ( \theta )，这个 ( \theta ) 就是 VSG 输出电压相位的参考值。为了频率静态误差，我通常还会在功率参考前并联一个调频器，也就是附加的功频下垂系数 ( k_f )，让 VSG 在一次调频过程中自动调整功率指令。

实现时有几个容易踩坑的地方：

• 单位问题：公式里的功率是标幺值还是实际值，积分器初始值、限幅值都要对应。我习惯全模型都用标幺值做控制，最后乘以基准值再进调制环节。
• 频率基准：( \omega_g ) 不能直接取 314.159 rad/s，如果做离网实验，应该取 VSG 自身输出电压频率的测量值作为反馈，否则负载突变时频率恢复会很慢。
• 积分限幅：虚拟角速度 ( \omega ) 的积分器必须设置限幅，通常限制在额定频率的 ±5% 范围内，否则仿真初始阶段可能出现瞬时频率飙升，数值上没有任何物理意义。

2.2 无功-电压环：模拟励磁调节特性

同步发电压机的励磁系统会维持机端电压恒定，无功功率增加时机端电压下降，励磁系统会自动增加励磁电流把电压顶回来。模拟这个特性，VSG 采用了无功-电压下垂控制：

[ E = E_0 + k_q (Q_{ref} - Q_e) ]

其中 ( E_0 ) 是空载电动势参考值，( Q_{ref} ) 是无功功率参考，( Q_e ) 是当前输出无功，( k_q ) 是无功下垂系数。为了模拟励磁绕组的时间延迟，通常会在 ( E ) 后面加一个一阶惯性环节，时间常数在 0.01~0.05 s 数量级。

这里要注意：( k_q ) 的取值直接影响电压静态特性。( k_q ) 太大，电压跌落严重；( k_q ) 太小，无功分配不均。工程上一般取额定无功从零增加到满载时电压跌落不超过 5%，也就是：

[ k_q = \frac{\Delta E}{\Delta Q} = \frac{0.05 E_{rated}}{Q_{rated}} ]

如果这套模型之后要做多机并联仿真，无功下垂系数就是各台 VSG 均分无功功率的关键参数。而且，计算出的 ( E ) 经过惯性环节之后，就是后续双闭环控制的电压外环参考指令。

2.3 虚拟阻抗与电压电流双闭环

VSG 输出的 ( E \angle \theta ) 并不直接作为逆变器的调制波，通常还要经过一个虚拟阻抗环节，得到实际的电压参考值：

[ v_{ref} = E \angle \theta - Z_v(s) \cdot i_o ]

为什么要加虚拟阻抗？因为实际逆变器的滤波阻抗非常小（尤其在采用了 LCL 滤波器之后），而同步发电机的同步电抗很大（通常 1~2 p.u.），两者外特性差异巨大。如果不加虚拟阻抗，VSG 控制的有功环和无功环会存在较强的耦合——有功变大时电压也跟着变，这就会削弱 VSG 的暂态稳定性。加入虚拟阻抗之后，可以把逆变器的等效输出阻抗往感性方向塑造，让有功-频率和无功-电压的耦合被有效隔离。

虚拟阻抗可以取纯感性（( L_v )），也可以取阻感性叠加。参数整定时，虚拟电抗一般取 0.05~0.15 p.u.，太小起不到解耦效果，太大则会导致电压跌落过多。需要注意的是，虚拟阻抗在算法上是一个微分/差分运算，放在 Simulink 里要小心高频噪声被放大，通常需要在电流采样回路上加滤波器，或者在虚拟阻抗环节本身加一个小时间常数的低通。

接下来是电压电流双闭环，这是让逆变器实际输出电压跟踪参考波形的关键。电压外环控制 VSG 输出端电压，电流内环控制滤波电感电流，两者都采用 PI 控制器。带宽设计原则是：电压环带宽控制在 100~300 Hz，电流环带宽控制在 500~1000 Hz，也就是相差 3~5 倍，这样才能保证内环响应足够快但又不给外环引入振荡。

双闭环实现上有一个常被忽略的细节：电流内环的输出是调制波，但如果直接在 dq 坐标系下用 PI 控制器，d 轴和 q 轴之间会存在耦合电压项 ( \omega L i_d )、( \omega L i_q )。常用的做法是加入解耦项：在控制器输出端把这两项补偿掉。Simulink 模型里如果不做解耦，仿真波形在负载突变后会出现 dq 轴电流交叉振荡，看起来像参数没调好，实际上是耦合项没有处理。

3. 光伏储能直流侧：从 MPPT 到母线电压稳定

VSG 逆变器要输出有功功率，前提是直流侧母线电压得稳得住。光伏阵列的输出天然波动，光照一变化直流母线电压就会跟着起伏，储能系统在这里起到“压舱石”的作用。

3.1 光伏阵列与 Boost 电路 MPPT

光伏阵列部分在 Simulink 里可以使用光伏电池模型，按实际组件的开路电压、短路电流、最大功率点参数设置串并联数量，搭出与需求容量匹配的光伏阵列。为了让阵列始终工作在最大功率点，我一般采用扰动观察法（P&O）：

对比当前功率与上一时刻功率，如果功率增加就维持扰动方向，如果功率减小就反转扰动方向。扰动步长需要权衡：步长太大，MPPT 稳态输出波动大；步长太小，光照突变时跟踪速度跟不上。Simulink 仿真里，建议先把步长设在 0.5%~1% 的参考电压范围，跑一次全阴影工况，观察功率曲线是否平滑。

Boost 电路负责把光伏阵列的低电压抬升到直流母线电压等级，开关管的 PWM 信号由 MPPT 控制器生成。其中有个关键问题：Boost 电感的大小直接影响电流纹波和 MPPT 跟踪精度。一般来说电感电流纹波控制在 10% 以内，就可以用：

[ L = \frac{V_{pv} \cdot D}{f_{sw} \cdot \Delta I_L} ]

估算出最小值，其中 ( D ) 是占空比，( f_{sw} ) 是开关频率，( \Delta I_L ) 是允许的电流纹波峰值。

3.2 储能双向 DC/DC 与母线电压控制

储能电池通过双向 DC/DC 变换器并联到直流母线。当光伏输出功率大于 VSG 发出的有功时，储能吸收多余功率，变换器工作在 Buck 模式；当光伏输出不足时，储能释放功率，变换器工作在 Boost 模式。母线电压控制采用电压外环电流内环的双闭环结构，电压外环的输出作为电流参考，内环控制电池电流跟踪参考值。

实际仿真里最常见的错误是储能直接并入直流母线而没有加双向变流器控制，后果是光伏功率波动全部倾泻到 VSG 逆变器上，母线电压出现剧烈波动，VSG 的虚拟惯量再大也顶不住。储能侧控制器的带宽可以设置得比 VSG 的惯量响应带宽高一些（比如 50~100 Hz），这样能保证母线电压快速稳定。此外，储能电池的 SOC 状态在 Simulink 模型中可以简化成恒压源串联内阻，但要观察电池动态响应，建议用 SOC 模块的充电/放电模型。

3.3 光伏储能功率协调逻辑

整个系统的功率平衡关系很简单但很重要：

[ P_{pv} + P_{bat} = P_{vsg} + P_{dc_loss} ]

光伏功率变化时，储能迅速补偿差值，让 VSG 输入功率尽量平滑。仿真中我用一个阶梯光照信号做测试：光伏功率每秒跳变一次，观察 VSG 输出有功是否平缓跟随。从波形上清晰地看到，如果没有储能，光伏跳变时 VSG 输出功率会出现明显的冲击——虚拟惯量只能延缓这个过程，而不能凭空产生能量，这正好说明了储能和 VSG 是互补配合的关系。

4. Simulink 模型搭建与参数整定实操

这一节讲模型的整体搭建流程和参数整定思路。按下面的步骤操作，跑出一个稳定的 VSG 并网波形基本没有问题。

4.1 模型整体架构与模块分层

我习惯把 Simulink 模型分成三层：

• 顶层为功率电路层：光伏阵列、Boost 电路、储能电池、双向 DC/DC、三相逆变器、LCL 滤波器、电网模型。用断路器模块模拟并网/离网切换。
• 第二层为控制层：MPPT 控制、储能控制、VSG 算法、电压电流双闭环、SVPWM（或 SPWM）调制。
• 第三层为信号测量层：电压电流测量、功率计算、频率计算、以及把这些信号通过 Goto/From 标签送到控制层。

分层的好处是：信号可以干净地导出到 MATLAB 工作区，方便后续做批处理参数扫描。

逆变器主电路参数（以 50 kW 系统为例）参考如下：

LCL 滤波器设计的目标是把开关频率附近的纹波衰减到 20% 以下，同时避免谐振峰过高。设计完成后必须加阻尼电阻或者采用有源阻尼，否则在 Simulink 仿真里很容易出现高频谐振发散。

4.2 VSG 控制模块的实现方式

Simulink 里实现 VSG 控制有两种常见方式：

第一种是纯模块搭：用积分器、增益模块、加法器把转子方程搭出来。优点是模型结构直接对应公式，适合教学；缺点是模块数量多，参数修改时要一个个点进去，做批量仿真时很痛苦。

第二种是用 S-Function 或者 MATLAB Function 写代码实现 VSG 主控制算法。代码方式更紧凑，参数导出到 MATLAB 脚本里一键赋值，方便做参数扫描和优化。我在实际仿真中更推荐代码方式，核心算法代码并不复杂，转子方程加上各种限幅和逻辑也就几十行。一个典型的 MATLAB Function 接口包括：输入 ( P_e )、( Q_e )、( \omega_g )，输出 ( V_d )、( V_q )、( \theta )。内部状态包括虚拟角速度、虚拟功角和电动势幅值，用持久变量保存积分器状态。

4.3 虚拟惯量与阻尼参数的计算方法

VSG 参数整定是整个模型最需要讲清楚的部分，因为 J 和 D 直接决定了暂态响应特性。

虚拟惯量 ( J ) 的含义是：当功率变化 ( \Delta P ) 时，虚拟角速度的变化率是 ( \Delta P / J )。如果取 ( J ) 过小，VSG 失去惯量支撑的初衷，频率变化率会太大；如果 ( J ) 过大，有功功率会出现长周期振荡，而且恢复时间拖得很长，这会让调频性能变差。

工程上初始取值建议从额定功率和额定角速度出发：

[ J = \frac{P_{rated}}{\omega_{rated}^2} \cdot T_j ]

其中 ( T_j ) 是虚拟惯性时间常数，物理意义是 VSG 从额定功率变化到零功率时，转速从额定值跌落到零所需的时间。同步发电机的惯性时间常数通常在 4~12 秒，VSG 的 ( T_j ) 可以取小一点，比如 2~8 秒，在模型里先取 4 秒试跑。

虚拟阻尼系数 ( D ) 的作用是提供功率振荡的阻尼力矩，合适时能快速衰减暂态振荡。D 的初值可以按照同步电机的阻尼比概念来估算，设期望的阻尼比为 ( \zeta = 0.7 )，则：

[ D = 2 \zeta \sqrt{\frac{J \cdot K}{1}} ]

其中 ( K ) 是功频下垂系数。不过说实话，仿真模型里直接扫描参数比手动计算有效率得多。我的做法是：把 J 和 D 设为模型工作区的变量，写一个循环，分别扫描 J 从 5 到 30（kg·m²）和 D 从 10 到 100 的范围，观察频率动态响应和功率超调量，选出综合效果最好的组合。这个方法比死磕公式速度更快，也更贴近实际调参流程。

4.4 并网预同步与启动流程

VSG 并网有一个特别容易出问题的场景，就是模型启动那一刻。如果 VSG 生成的电压幅值和相位与电网电压不一致，断路器闭合瞬间会产生巨大的冲击电流，仿真立刻发散。所以我在模型里加了预同步控制流程：

1. 检测电网电压幅值、频率和相位。
2. 把 VSG 输出电动势幅值指令设为电网电压幅值。
3. 把虚拟角速度参考值设为电网角速度。
4. 通过相位调节器（一个小 PI 控制器）逐步把 VSG 输出电压相位拉向电网相位。
5. 当相位差小于 0.1 度、幅值差小于 1% 时，闭合断路器。

这个流程在 Simulink 里实现需要加一个启停逻辑信号，通过断路器模块的外部控制端控制合闸。如果没有预同步，即使仿真软启动把功率斜坡给到额定值，合闸后的冲击电流也会让模型在前面几步就已经崩溃。

5. 仿真调试中的典型问题与排查经验

下面这组问题是我在实际复现和调试 VSG 模型时踩过的坑，每条都对应具体的现象和解决思路，做仿真遇到类似情况可以直接对着排查。

5.1 启动瞬间电流飙升

现象：模型一开始运行，逆变器电流直接冲到几千安培，仿真报错或波形出现不可收敛的尖峰。

原因分析：大概率是预同步没有做彻底，或者 VSG 的电动势幅值初始化和电网电压不一致。也有一种可能是断路器闭合瞬间 VSG 相位还在快速变化，PI 调节器没有跟踪到位。

解决方法：

• 先把预同步控制时长从 0.5 s 延长到 1.5 s，确保相位差降到极小值再合闸。
• 在合闸前的 VSG 控制里，把虚拟机械功率 ( P_m ) 初始值设成 0，不要让 VSG 在空载状态突然输出额定功率。
• 检查软件限幅：电动势幅值最好不要超过电网电压幅值的 1.2 倍。

5.2 VSG 有功功率低频振荡

现象：功率响应曲线出现明显振荡，振荡频率在 2 Hz 左右，持续很久不衰减。

原因分析：虚拟惯量 J 相对虚拟阻尼 D 偏大，系统阻尼比不足。或者虚拟阻抗取值不合适，导致有功环和无功环耦合，外环之间相互激发振荡。

解决方法：

• 按前面说的参数扫描方法，增大 D 或者减小 J，把阻尼比调到 0.7 左右。
• 检查虚拟阻抗是否取了纯阻性，如果是，把 ( R_v ) 调小，( L_v ) 调大，让阻抗角往感性偏移。
• 检查电压外环带宽是不是和 VSG 的机械动态带宽太接近。VSG 的惯量响应通常集中在 0.5~5 Hz，电压环带宽在 100 Hz 以上才安全。

5.3 锁相环与 VSG 的相位冲突

现象：模型并网稳态运行正常，但电网电压发生小幅跌落时，VSG 输出功率出现短时振荡甚至失步。

原因分析：很多模型里 VSG 和锁相环（PLL）是同时存在的——PLL 用来检测电网相位，VSG 内部又在生成自己的相位，二者天然存在动态耦合。并网模式下，如果把 PLL 测到的电网相位直接当成 VSG 的角度基准，在电网扰动时 PLL 的动态响应会和 VSG 的虚拟角速度积分发生相互作用，产生振荡。

解决方法：VSG 并网模式下实际上不需要让 PLL 直接控制输出相位，PLL 只用于检测频率和预同步阶段的对相。运行时让 VSG 的虚拟转子自己积分生成相位，只有频率测量需要从 PLL 取。这在结构上就避免了相位双闭环冲突。如果必须保留 PLL 做保护逻辑，把 PLL 带宽降低到 5 Hz 以下，可以减弱耦合效应。

5.4 仿真步长与离散化带来的数值不稳定

现象：模型在变步长求解器下波形正常，但固定步长跑一定时间后，高频部分出现数值振荡甚至发散。

原因分析：Simulink 在固定步长离散系统下，LCL 滤波器的谐振峰会在采样频率附近产生虚假谐振。或者 PWM 载波频率和仿真步长整数倍关系不匹配，导致 PWM 波形出现次谐波。

解决方法：

• 固定仿真步长必须至少满足奈奎斯特条件。10 kHz 开关频率下，仿真步长建议不大于 ( 1/(10 \times f_{sw}) )，也就是 1e-6 s。这是一回事吗？其实 ( 1/(10 \times 10000) = 1e-5 ) s，也就是 10 μs。如果条件允许，用 5 μs 更稳。
• LCL 滤波器加阻尼电阻，或者在数字实现里面用有源阻尼算法替换纯微分项。
• 如果只关心 VSG 外特性，可以采用平均值模型替代 PWM 细节，把三相桥等效成受控电压源，仿真步长就能放大到 50~100 μs，速度提升明显。

下面把这些问题整理成一张速查表，方便建模时对照：

5.5 模型调试中的“代数环”问题

如果你用纯模块搭 VSG 控制环，很可能会碰到代数环（Algebraic Loop）报错。代数环出现的原因通常是：某个模块的输出经过几层反馈网络又回到了输入端，而路径上没有任何记忆元件（积分器、延迟单元）打破这个循环。转子方程里的 ( P_e ) 测量和 ( \theta ) 输出之间如果中间没有状态量隔离，Simulink 解算代数环时会大幅降低仿真速度，严重时甚至无法求解。

解决的办法有两种：一是在反馈路径上放一个 Memory 模块或者一阶滤波器，打破代数环；二是把所有控制算法放进 MATLAB Function 里面，利用内部状态变量天然的记忆特性来避免代数环。我的经验是：只要控制链路稍复杂，尽量用 MATLAB Function 实现 VSG 核心算法，这是最省心的方案。

6. 这套模型能扩展到哪些方向

模型跑通之后，可以扩展的实验方向还不少。如果你们是科研用途，比较有价值的方向包括：

第一是 VSG 参数对系统稳定性的影响分析。利用这个仿真模型做扫参实验，观察不同虚拟惯量、阻尼系数下频率最大偏差和恢复时间的变化规律，画出三维稳定域图，这本身就是一个完整的系统分析成果。

第二是多机并联与功率分配。在现有模型上复制几套 VSG，通过线路阻抗连接到同一个交流母线，研究不同线路阻抗和虚拟阻抗下的功率均分效果。并联台数多了之后，传统下垂的灵敏度问题和 VSG 的全局稳定性问题都会浮现出来，非常能锻炼叠代分析能力。

第三是孤岛与并网模式切换。在电网模型前面加一个断路器，模拟电网断电和恢复的完整过程，观察 VSG 是否能在孤岛模式下带负载稳定运行，在重新并网时又能通过预同步恢复。这套控制逻辑和实验在真实装置上做成本高、风险大，但是用仿真模型来做却是又安全又直观。

从我个人这些年的仿真经验来看，这套模型的价值主要在于它把“新能源发电侧动态特性”和“电网支撑需求”放到了同一个时间尺度下观察。以前你看到光伏功率波动只会想到 MPPT 波形怎么调，加上 VSG 之后会意识到：真正决定并网品质的，是逆变器控制策略和直流侧能量管理之间的配合。把一个 VSG 仿真模型调稳、调透之后，你对电力电子变流器控制的认识会明显超过只做单机双闭环的阶段，至少看电网频率的视角会完全不同。