VSG阻抗建模与仿真：从控制原理到弱电网稳定性分析实战

简介：本资源是一套面向电力电子与新能源并网方向研究生及科研工程师的VSG逆变器序阻抗建模MATLAB仿真工具包，聚焦弱电网下虚拟同步发电机并网稳定性分析这一核心问题，特别适用于阻抗建模、扫频法验证与小信号稳定性研究等场景。压缩包共9个文件（6个.m脚本、2个.txt说明文档、1个.slx Simulink模型），总大小仅66KB，轻量紧凑；其中包含完整可运行的VSG功率环仿真模型、正负序阻抗计算主程序（ZpZncalc.m）、扫频激励与频响提取脚本（FreqScan_test.m、FFTread.m等），以及关键参数配置与运行说明文档，所有代码均附详细中文注释。已有144人学习下载，用户可直接复现文献《虚拟同步发电机接入弱电网的序阻抗建模与稳定性分析》中的建模流程与扫频结果，在中高频段获得与理论高度一致的阻抗特性曲线，快速开展VSG并网谐振风险评估与控制器优化设计。

1. 项目缘起：从“并网稳定”到“阻抗建模”的思考

最近在调试一个虚拟同步发电机（VSG）控制的储能逆变器项目时，遇到了一个挺典型的问题：在弱电网环境下，系统偶尔会出现频率或电压的轻微振荡，虽然没到失稳的地步，但总让人心里不踏实。排查了一圈硬件和基础控制参数，都没发现明显问题。后来和团队里的老工程师讨论，他提了一句：“你测过它在正序和负序扰动下的阻抗特性吗？VSG的阻抗和传统电流源型逆变器可不太一样。” 这句话点醒了我。我们通常关注VSG的下垂控制、惯性模拟这些“功”层面的特性，却容易忽略其在“阻抗”这个“防”层面的表现。而电网的扰动，尤其是三相不平衡这类含有负序分量的扰动，恰恰会与逆变器的阻抗特性发生交互，影响并网稳定性。

于是，我决定放下手头的调试，回头补补课，系统地研究一下VSG的正负序阻抗建模与仿真。这个事不能只停留在理论推导，必须得在MATLAB/Simulink里把它“跑”出来，直观地看到阻抗曲线，分析其与电网阻抗的交互关系，这样才能对我们实际产品的参数整定和稳定性评估有直接的指导意义。这篇文章，我就把自己搭建模型、进行阻抗扫描仿真的完整过程、核心原理以及踩过的几个坑，详细地记录下来。无论你是正在研究VSG的学生，还是像我一样遇到类似工程问题的工程师，希望这份“实战笔记”都能给你提供一个清晰的参考路径。

2. VSG控制核心与阻抗特性的内在关联

要理解VSG的阻抗，必须先吃透它的控制核心。很多人把VSG简单地理解成“加了惯性和下垂的逆变器”，这没错，但要想建模，就得拆开看它的信号流。VSG的核心思想是模拟同步发电机的转子运动方程和励磁调节，其功率环（外环）通常由以下部分构成：

2.1 有功-频率环：模拟转子惯性这部分的核心方程是大家熟悉的摇摆方程：J * dω/dt = Pm - Pe - Dp*(ω - ωg)其中，J是虚拟转动惯量，Dp是阻尼系数，Pm是机械功率（对应参考有功功率Pref），Pe是电磁功率（对应实际输出有功功率P）。这个方程决定了VSG的频率动态响应。当电网频率变化时，VSG不会像电流源型逆变器那样立刻跟随，而是会通过储存或释放“虚拟动能”来提供惯性支撑，响应速度受J和Dp影响。这一点，直接影响了其在频域上的阻抗特性，尤其是在工频附近的表现。

2.2 无功-电压环：模拟励磁调节电压环通常采用下垂控制：E = E0 + Dq*(Qref - Q)这里，E是VSG输出的空载电动势幅值，E0是额定电压，Dq是无功-电压下垂系数，Qref和Q分别是无功参考值和实际值。这个环决定了VSG的电压调节特性。当并网点电压跌落时，VSG会增大无功输出以支撑电压。

2.3 内环控制：电流环或功率直接计算VSG的外环输出（频率ω和电压幅值E）经过坐标变换后，生成电压指令。这个电压指令通常有两种处理方式：

1. 直接作为调制波：这种方式结构简单，但抗扰性稍差，相当于电压源特性更明显。
2. 作为电流环的给定：外环输出的电压指令经过一个虚拟阻抗（通常是感抗）后，转化为电流参考值，再经过快速的电流内环跟踪。这是目前更主流、性能也更优的方式，因为它通过电流环限流，更安全，且能明确塑造输出阻抗。

我们重点讨论第二种方式。这里引入的“虚拟阻抗”Zv = Rv + jωLv是一个关键设计参数。它并不代表实际的物理阻抗，而是一个控制算法中的概念。它的作用主要有两个：一是帮助实现功率解耦（尤其是在dq坐标系下），二是主动塑造VSG的输出阻抗。是的，你可以通过设计Rv和Lv，来改变VSG在特定频段呈现出来的阻抗大小和相位，这是优化其与电网阻抗交互、提升稳定性的重要手段。

2.4 从控制框图到阻抗模型理解了上述结构，我们就可以建立从扰动（并网点电压的dΔV，qΔV）到响应（输出电流的dΔI，qΔI）的小信号模型。这个过程需要：

1. 在稳态工作点（某个特定的有功P0、无功Q0、电网电压Vg0下）对VSG的非线性方程进行线性化。
2. 考虑锁相环（PLL）的动态。对于VSG，其频率由自身控制产生，PLL通常只用于初始同步和角度微调，其动态影响与传统电流源逆变器相比可能不同，但建模时仍需纳入。
3. 将功率计算、虚拟阻抗、电流环等环节的传递函数全部考虑进去，最终得到一个多输入多输出（MIMO）的传递函数矩阵G(s)，使得[ΔId; ΔIq] = G(s) * [ΔVd; ΔVq]。
4. 这个G(s)的逆，在某种意义上就反映了VSG的输出导纳（阻抗的倒数）。但更常用的方法是，通过施加正序或负序的电压扰动，并测量电流响应，直接在频域计算阻抗Z(s) = ΔV(s) / ΔI(s)。

3. 在Simulink中构建可进行阻抗扫描的VSG仿真模型

理论推导是基础，但仿真才是验证和获得直观认识的利器。我们的目标是在Simulink中搭建一个VSG模型，并能够方便地注入扰动信号，进行频率扫描，自动计算出正序和负序阻抗。下面是我的建模步骤和关键点。

3.1 VSG主电路与控制模块搭建首先，建立一个典型的三相两电平电压源型逆变器模型，直流侧接一个稳定的直流电压源。控制部分按照第二节的描述搭建：

• 功率计算模块：采用基于瞬时功率理论的p-q计算或低通滤波后的有功无功，确保在扰动下也能正确计算。
• VSG核心算法模块：实现摇摆方程和电压下垂方程。这里注意，摇摆方程需要离散化（如采用前向欧拉法），并合理设置积分步长。
• 虚拟阻抗模块：在dq坐标系下实现Vd_ref' = Vd_ref - (Rv*Id + ω*Lv*Iq)，Vq_ref' = Vq_ref - (Rv*Iq - ω*Lv*Id)。这里的Vd_ref， Vq_ref是VSG算法输出的内电势分量。
• 电流内环模块：采用典型的PI控制器，输出调制电压。PI参数需要根据开关频率和滤波器参数精心设计，保证足够的带宽和稳定性裕度。
• 坐标变换与PLL模块：使用基于二阶广义积分器（SOGI）的锁相环或同步参考系锁相环（SRF-PLL）。对于阻抗扫描，PLL的带宽设置很重要，带宽过高可能会在阻抗曲线中引入额外的相位滞后。

关键技巧1：模型初始化VSG模型启动时需要稳定的初始状态。我的做法是：先让逆变器以V/f控制模式（即开环输出电压幅值和频率）启动，建立空载电压。然后，在并网开关闭合前，将VSG控制器的积分器初始值设置为当前状态（如频率初始化为电网频率，电压初始化为并网点空载电压），实现平滑切换。这能避免合闸冲击，也让阻抗扫描从一个合理的稳态工作点开始。

3.2 正负序扰动注入与测量电路这是阻抗扫描仿真的核心。我们希望在并网点（PCC）注入一个幅值很小、频率可变的电压扰动，同时测量扰动引起的电流响应。

• 扰动注入点：在电网阻抗（用一个RL串联电路模拟）之后、PCC点之前，串联一个受控电压源。这个电压源用于注入扰动信号。
• 扰动信号生成：我们需要生成正序和负序的扰动。对于一个频率为f_inj（相对于基波50Hz的偏移）的正序扰动，其三相表达式为：Va_inj = A * sin(2*pi*(50+f_inj)*t)Vb_inj = A * sin(2*pi*(50+f_inj)*t - 2*pi/3)Vc_inj = A * sin(2*pi*(50+f_inj)*t + 2*pi/3)负序扰动则将B相和C相序互换。 在Simulink中，可以用三个Sine Wave模块，通过设置不同的相位来生成。扰动幅值A要足够小（如额定电压的1%-5%），以确保系统响应在线性范围内。
• 测量环节：精确测量PCC点的三相电压和逆变器输出的三相电流。使用Simulink的“电压测量”和“电流测量”模块。务必注意测量点的位置，电压测量点应在扰动注入点之后（即PCC点），电流测量点应在逆变器输出侧。

3.3 阻抗计算与频扫自动化实现手动改变扰动频率、运行仿真、记录数据是非常低效的。我们需要利用MATLAB脚本来自动化这个过程。

1. 参数化模型：将扰动频率f_inj、扰动类型（正序/负序）等设置为Simulink模型的工作区变量（如f_inj，seq_type）。
2. 编写频扫脚本：在MATLAB中编写一个.m脚本文件。

   % 定义频率扫描范围，例如从1Hz到1000Hz，对数均匀取点 f_list = logspace(0, 3, 50); % 10^0 到 10^3 Hz，取50个点 Z_pos = zeros(length(f_list), 1); % 存储正序阻抗 Z_neg = zeros(length(f_list), 1); % 存储负序阻抗 % 设置仿真参数 sim_time = 1; % 仿真时间，要保证达到稳态 Amp_inj = 10; % 扰动幅值 (V) for i = 1:length(f_list) f_inj = f_list(i); % 1. 设置正序扰动并仿真 seq_type = 1; % 1代表正序 sim('VSG_Impedance_Model'); % 运行Simulink模型，模型名需对应 % 等待仿真稳定后，取最后几个周期的数据进行FFT分析 [Z_pos_mag, Z_pos_phase] = calculate_impedance_from_simout(t, Vpcc, Iinv, f_inj, 50); Z_pos(i) = Z_pos_mag * exp(1j*Z_pos_phase); % 2. 设置负序扰动并仿真 seq_type = -1; % -1代表负序 sim('VSG_Impedance_Model'); [Z_neg_mag, Z_neg_phase] = calculate_impedance_from_simout(t, Vpcc, Iinv, f_inj, 50); Z_neg(i) = Z_neg_mag * exp(1j*Z_neg_phase); end % 3. 绘制阻抗曲线（伯德图或奈奎斯特图） figure; subplot(2,1,1); semilogx(f_list, 20*log10(abs(Z_pos)), 'b-', 'LineWidth', 1.5); hold on; semilogx(f_list, 20*log10(abs(Z_neg)), 'r--', 'LineWidth', 1.5); ylabel('阻抗幅值 (dBΩ)'); grid on; legend('正序阻抗', '负序阻抗'); subplot(2,1,2); semilogx(f_list, angle(Z_pos)*180/pi, 'b-', 'LineWidth', 1.5); hold on; semilogx(f_list, angle(Z_neg)*180/pi, 'r--', 'LineWidth', 1.5); ylabel('阻抗相位 (度)'); xlabel('频率 (Hz)'); grid on;

3. 设计calculate_impedance_from_simout函数：这个函数是数据处理的核心。其思路是：
   • 对采集到的时域电压电流信号进行FFT，精确提取在扰动频率f_inj处的正序或负序分量。这里需要用到对称分量法。
   • 对于正序阻抗，计算Z_pos = FFT(V_abc @ f_inj)_positive / FFT(I_abc @ f_inj)_positive。
   • 对于负序阻抗，计算Z_neg = FFT(V_abc @ f_inj)_negative / FFT(I_abc @ f_inj)_negative。
   • 注意处理基波分量（50Hz）的影响，通常需要加窗（如汉宁窗）并选择适当的FFT长度以提高频率分辨率。

关键技巧2：确保线性响应与稳态阻抗扫描的前提是小信号线性化。因此，扰动幅值必须足够小。如何判断？可以做一个简单的测试：将扰动幅值减半，重新扫描某个频率点，如果计算出的阻抗值变化很小（<5%），则认为工作在线性区。另外，每次改变频率后，必须保证仿真时间足够长，让系统的瞬态响应完全衰减，只留下对注入扰动的稳态响应。对于低频段（如几Hz），可能需要更长的仿真时间。

4. 仿真结果分析与典型VSG阻抗特性解读

运行完自动化脚本，我们就能得到VSG在特定工作点（P0, Q0）下的正负序阻抗曲线。下图展示了一组典型的仿真结果（基于一个具体参数集的示意性曲线）。

4.1 正序阻抗的“感-容-感”变化仔细观察正序阻抗相位曲线，你可能会发现一个有趣的现象：从极低频的感性，到中频段可能变为容性，再到高频段又回到感性（由LCL滤波器决定）。这个“感-容-感”的变化是VSG区别于传统电流源逆变器（通常呈现为“容性”或“阻容性”）的关键。中频段的容性特性，源于电流内环的相位滞后和PLL的动态。这个容性区域如果与电网的感性阻抗在幅值上相交且相位差满足奈奎斯特稳定判据，就可能引发谐振失稳，这也是弱电网下VSG容易出问题的频段。

4.2 负序阻抗为何不同？负序阻抗与正序阻抗之所以不同，根本原因在于VSG控制是对称分量法下的“序解耦”并不完全。

1. 功率环的影响：VSG的有功、无功功率计算是基于正序分量的（或者说是基于正序旋转坐标系的）。负序电压/电流分量在功率计算中会产生二倍频波动，但经过控制环的低通滤波后，其对功率环的稳态影响很小。因此，功率环（摇摆方程和电压下垂）对负序扰动的“调节”作用很弱。
2. 电流内环与虚拟阻抗：电流内环在dq坐标系下是对称的，理论上对正负序电流的跟踪能力相同。虚拟阻抗Zv也是对称的。因此，在中高频段，负序阻抗主要由电流环和虚拟阻抗决定，可能与正序阻抗中高频部分相似。
3. PLL的影响：PLL主要跟踪正序电压相位。负序分量会导致PLL输出中出现二倍频波动，这个波动会通过坐标变换耦合到控制中，影响负序阻抗，尤其是在频率接近二倍基频（100Hz）的区域。

因此，负序阻抗在低频段更像一个被动的阻抗（主要由滤波器和虚拟阻抗决定），而在中高频段逐渐与正序阻抗趋同。负序阻抗特性对于评估VSG在电网电压不平衡条件下的运行能力至关重要，比如能否抑制负序电流、是否会引起不平衡工况下的振荡。

4.3 关键参数对阻抗曲线的影响通过改变VSG参数重复扫描，可以直观看到其影响：

• 增大虚拟惯量J：会显著提升极低频段（<2Hz）正序阻抗的幅值，使VSG在频率扰动时“更硬”，惯性响应更慢。
• 增大阻尼系数Dp：会降低极低频段阻抗的幅值，并使相位更偏向阻性，有助于抑制低频振荡。
• 增大虚拟电感Lv：会整体提升中低频段阻抗的幅值（使其更感性），这有助于增强对高频谐波的衰减，但可能降低系统在弱电网下的稳定性（因为感抗叠加）。
• 改变电流环带宽：电流环带宽越高，其中频段阻抗的“容性凹陷”区域越往高频移动，且凹陷程度可能变浅。

5. 基于阻抗模型的稳定性分析与工程应用启示

得到阻抗曲线不是终点，而是为了应用。最重要的应用之一就是基于阻抗比的稳定性判据。

5.1 奈奎斯特稳定判据与阻抗比曲线对于并网逆变器系统，可以从源（电网阻抗Zg(s)）和负载（逆变器阻抗Zinv(s)）的角度建模。稳定性判据可以简化为：观察环路增益L(s) = Zinv(s) / Zg(s)的奈奎斯特曲线是否包围(-1, j0)点。更工程化的方法是绘制阻抗比Zg(s) / Zinv(s)的奈奎斯特曲线，看其是否包围(-1, j0)点。或者，更直观地，在伯德图上比较Zinv和Zg的幅频特性，并检查相位差。

• 稳定条件：在所有频率下，当|Zinv| > |Zg|（即逆变器阻抗更大）时，系统稳定。如果存在某个频率点f_c使得|Zinv(f_c)| = |Zg(f_c)|，则需检查该频率点处的相位差。若Phase(Zinv) - Phase(Zg)的绝对值小于180°，则稳定；否则不稳定。

在我们的仿真中，Zinv就是测得的VSG正序阻抗。我们可以假设或测量一个典型的电网阻抗Zg（通常是感性的，Zg = jωLg）。将两者绘制在同一张伯德图上进行对比。

5.2 弱电网下的稳定性风险识别假设电网短路比（SCR）较低，即Lg较大，Zg的幅值曲线会整体上移。这时，它很可能与VSG阻抗曲线在中频段（比如几十到几百Hz）相交。如果相交点恰好位于VSG阻抗呈现容性相位的区域，而Zg是感性相位，两者的相位差可能接近180°，这就满足了振荡的条件。仿真结果可以清晰地预测这个风险点。

5.3 指导控制器参数优化阻抗模型为我们提供了量化优化控制器参数的依据。例如，如果我们发现系统在50Hz附近存在稳定性风险，可以：

1. 调整虚拟阻抗Zv：增加虚拟电阻Rv可以在中频段为阻抗引入正的电阻分量，提高阻尼，压制谐振峰。这是最直接有效的手段之一。
2. 优化电流环参数：适当降低电流环的带宽或调整其相位补偿，可以改变中频段阻抗“容性凹陷”的位置和深度，使其避开与电网阻抗的危险交互区域。
3. 引入有源阻尼策略：在控制中引入基于电容电流或并网电流反馈的阻尼项，这相当于在阻抗模型中并联了一个虚拟电阻，可以有效重塑高频段阻抗。

5.4 对产品设计与测试的启示

1. 设计阶段：在确定VSG控制参数（J， Dp， Lv， Rv， 电流环PI）时，不应只关注空载波形和阶跃响应，必须进行阻抗扫描仿真，评估其在典型弱电网阻抗谱下的稳定性边界。可以形成一套参数设计的“阻抗约束”准则。
2. 测试验证：在实验室环境下，可以使用电网模拟器或阻抗扫描仪，实际测量VSG样机的输出阻抗，与仿真结果进行对比校验。这是验证模型准确性和产品可靠性的关键一步。
3. 故障穿越：电网电压跌落时，VSG的阻抗特性会因工作点突变而动态变化。可以基于阻抗模型，分析在不同跌落深度下，系统是否仍能保持稳定，从而优化故障穿越控制策略。

通过这一整套从建模、仿真到分析、应用的过程，我们就能将抽象的“阻抗”概念，转化为具体指导VSG控制器设计、优化和稳定性评估的强力工具。它让我们不再盲目试凑参数，而是有的放矢地塑造逆变器的并网特性，使其更智能、更稳健地融入未来的高比例新能源电网。

本文还有配套的精品资源，点击获取