📚 课程目录第4章 VSG有功-频率控制原理

4. VSG有功-频率控制原理

虚拟同步机（VSG）之所以能在电力电子变换器上“模拟”出同步发电机的频率支撑行为，核心在于它把同步机的转子运动方程搬进了控制器。理解这一节，关键不是记住公式，而是弄清三个物理量——虚拟惯量、虚拟阻尼、一次调频——各自在频率动态中扮演什么角色，以及如何把连续时间的转子方程变成数字控制器里能一步步执行的差分方程。

一、从同步机转子运动方程到VSG控制方程

同步发电机的转子运动方程（忽略机械损耗）为：

J·dω/dt = Tm − Te − D·(ω − ω0)

其中 J 为转动惯量，ω 为转子角频率，Tm 为机械转矩，Te 为电磁转矩，D 为阻尼系数，ω0 为额定角频率。VSG 把这套方程移植到逆变器控制中：用“虚拟惯量 J”替代真实转子惯量，用“虚拟阻尼 D”替代机械阻尼，用有功功率指令与实测功率之差构造“虚拟转矩”，从而让逆变器输出频率具有惯性响应和阻尼振荡的能力。

工程上更常用有功—频率形式表达。设极对数为1，机械角频率与电角频率一致，将转矩乘以 ω 近似为功率，可得常用的 VSG 有功-频率控制方程：

J·ω0·dω/dt = Pref − Pe − D·ω0·(ω − ω0) − Kf·(ω − ω0)

式中 Pref 为有功功率指令，Pe 为实测输出有功，Kf 为一次调频（下垂）系数。这里需要先明确符号约定：本文及图示均以 ω 高于 ω0 为正偏差，功率以逆变器输出为正方向。当频率升高时，VSG 应减少有功输出，因此下垂项和阻尼项都表现为“负反馈”。

二、虚拟惯量、虚拟阻尼与一次调频的分工

三个参数在频率动态中的职责完全不同，混淆它们会导致整定方向错误。

虚拟惯量 J：决定频率变化的“缓急”。J 越大，频率变化率 df/dt 越小，频率跌落或上升越平缓，但响应变慢，过大时可能引起功率振荡。它对应的是频率的微分行为。
虚拟阻尼 D：抑制频率偏差和功率振荡，相当于给系统增加“摩擦”。D 越大，频率偏差被压得越小，但过大会使动态变迟钝，并影响稳态功率分配。
一次调频系数 Kf：决定频率稳态偏差与功率增量的比例关系，即下垂特性。它对应的是频率的比例行为，决定稳态时频率偏离额定值多少、功率变化多少。

可以这样理解：J 管“过程快慢”，D 管“振荡强弱”，Kf 管“稳态落点”。三者共同作用，才能让 VSG 既有惯性支撑，又能参与一次调频。

三、离散化实现方法

数字控制器（DSP 或 MCU）无法直接求解连续微分方程，必须离散化。最常用的是前向欧拉法或后向欧拉法。以采样周期 Ts 对转子方程离散，采用前向欧拉：

ω[k+1] = ω[k] + (Ts / (J·ω0)) · [Pref − Pe[k] − D·ω0·(ω[k] − ω0) − Kf·(ω[k] − ω0)]

进一步整理为增量形式，便于定点实现：

Δω = (Ts / (J·ω0)) · (Pref − Pe[k] − (D·ω0 + Kf)·(ω[k] − ω0))
ω[k+1] = ω[k] + Δω
θ[k+1] = θ[k] + ω[k+1]·Ts

其中 θ 为虚拟转子角，用于生成调制波相位。实际实现时，通常把 (D·ω0 + Kf) 合并为一个等效阻尼系数，减少一次乘法。若采用后向欧拉，则需把 ω[k+1] 代入方程右侧，得到隐式差分，稳定性更好但计算稍复杂。

离散化时有三个工程要点必须注意：

采样周期 Ts 的选择：Ts 应远小于系统最小时间常数。若 J 很小而 Ts 偏大，前向欧拉可能发散。一般要求 Ts ≤ 0.1·J·ω0/(D·ω0+Kf) 量级。
功率测量滤波：Pe 通常来自低通滤波后的瞬时功率，滤波延迟会等效增加阻尼或降低相位裕度，整定 D 时需一并考虑。
限幅与抗积分饱和：频率偏差和功率指令都需限幅，防止故障或大扰动下虚拟转子角越限导致失步。

四、参数整定示例

以下为一组教学示例参数（非实测结果，仅用于说明整定思路）：额定频率 50 Hz，ω0 = 314.16 rad/s，逆变器额定有功 10 kW，采样周期 Ts = 100 μs。若希望频率跌落 0.5 Hz 时增发 5 kW 有功，则一次调频系数：

Kf = ΔP / Δω = 5000 / (2π·0.5) ≈ 1592 W·s/rad

虚拟惯量可按“等效惯性时间常数”选取。若希望 VSG 的等效惯性时间常数 H = 2 s，则：

J = 2·H·S / ω0² = 2·2·10000 / 314.16² ≈ 0.405 kg·m²

虚拟阻尼 D 通常按保证系统阻尼比 ζ 在 0.6～0.8 之间选取，可先取 D·ω0 与 Kf 同量级，再通过阶跃扰动观察频率超调和恢复时间微调。

整定顺序建议：先定 Kf（决定稳态下垂），再定 J（决定频率变化率），最后调 D（抑制振荡）。不要同时大范围改动三个参数，否则难以判断是哪个参数在起作用。

五、频率动态过程与结果判读

当负载突增时，VSG 的频率响应可分为三个阶段：首先是 J 主导的惯性阶段，频率缓慢下降；随后 D 与 Kf 共同作用，频率下降速度减缓并出现最低点；最后进入稳态，频率稳定在由 Kf 决定的下垂点。若观察到频率持续等幅振荡，说明 D 偏小；若频率恢复过慢且超调大，说明 J 或 D 偏大；若稳态频率偏差与预期不符，应检查 Kf 和功率测量标度。

上图把连续域方程、三个参数的分工、离散化实现和频率动态响应放在同一张图中对照。可以看到，J 决定曲线下降的斜率，D 和 Kf 共同决定最低点深度与恢复速度，而 Kf 最终决定稳态频率偏差 Δf_ss。离散化实现则把这一连续过程拆成每个采样周期执行一次的增量计算，使 VSG 能够在数字控制器中实时运行。

离散化实现时，若采用前向欧拉且 Ts 过大，可能出现数值发散；若采用后向欧拉，稳定性更好但需注意计算延迟。无论哪种方法，都应对 ω 和 θ 设置限幅，防止大扰动下虚拟转子角越限导致失步。

总结来说，VSG 有功-频率控制的核心就是：用 J 模拟转子惯性，用 D 模拟阻尼，用 Kf 实现一次调频，再通过离散化把连续转子方程变成控制器可执行的差分方程。整定时先定 Kf 确定稳态下垂，再定 J 确定频率变化率，最后调 D 抑制振荡，并通过频率动态曲线判断参数是否合适。