在新能源并网领域，光伏逆变器与虚拟同步发电机的稳定性分析一直是研究热点。当电网强度较弱时，这两种并网接口的阻抗特性差异会显著影响系统的稳定裕度。本文将通过Simulink仿真，深入对比分析光伏并网逆变器与虚拟同步发电机的正负序阻抗建模方法，并探讨其在弱电网下的稳定性差异。无论你是电力电子方向的研究生，还是从事新能源并网开发的工程师，都能从本文获得一套完整的建模、仿真与稳定性分析流程。

1. 背景与核心概念

1.1 弱电网稳定性挑战

随着风电、光伏等分布式电源大规模接入，电力系统呈现出“弱电网”特性。所谓弱电网，通常指电网短路比较低、等效阻抗较大的情况。在这种场景下，并网逆变器与电网之间的动态交互变得复杂，容易引发宽频带振荡等稳定性问题。因此，精确建立并网接口的阻抗模型，是分析并预判系统稳定性的关键。

1.2 正负序阻抗建模的意义

在不对称运行或故障条件下，电网电压和电流会包含正序、负序和零序分量。对于三相三线制系统，零序分量通常不存在。正负序阻抗模型能够分别描述并网设备在正序和负序扰动下的动态响应，比传统的单输入单输出阻抗模型更能全面反映其在非理想电网条件下的特性。通过对比光伏逆变器（通常采用电网跟随控制）和虚拟同步发电机（模拟同步机惯性与阻尼）的序阻抗，可以深入理解不同控制策略对弱电网稳定性的影响机制。

1.3 技术路线概述

本文将遵循“理论建模-仿真实现-对比分析”的技术路线。首先，分别推导光伏并网逆变器与虚拟同步发电机的正负序阻抗数学模型。接着，在Simulink中搭建详细的仿真模型进行验证。最后，通过频域扫描（如阻抗比判据）和时域仿真，对比两者在典型弱电网参数下的稳定性表现。

2. 仿真环境准备与模型说明

2.1 所需软件与工具

本文仿真基于以下环境，读者可根据自身情况调整版本。

• 核心软件: MATLAB/Simulink (本文示例基于 R2021a， 建议使用 R2018b 及以上版本)。
• 必要工具箱: Simscape Electrical (原名 SimPowerSystems)， 用于搭建电力电子与电网模型。
• 辅助分析: MATLAB 脚本用于数据处理和频域分析（如绘制奈奎斯特图）。

2.2 示例模型总体结构

我们将创建两个独立的Simulink模型：PV_Inverter_Model.slx和VSG_Model.slx。它们共享相同的电网和负载条件，以便公平对比。模型主要包含以下部分：

1. 主电路: 直流电压源、三相逆变桥、LCL滤波器、变压器、电网等效阻抗（用于模拟弱电网）和无穷大电网。
2. 控制系统:
   • 光伏逆变器: 基于双闭环（电流内环、电压外环）的电网跟随控制，包含锁相环。
   • 虚拟同步发电机: 包含有功-频率下垂控制、无功-电压下垂控制，以及模拟转子运动方程的虚拟惯量环节。
3. 测量与扰动注入模块: 用于测量端口电压、电流，并在特定频率点注入正序或负序小信号扰动，以计算阻抗。

2.3 关键参数设置（示例）

以下参数为示例值，实际分析中需根据具体案例调整。

% 文件：init_parameters.m % 系统基准值 f_base = 50; % Hz， 基波频率 V_base = 220; % V， 相电压有效值（线电压380V） S_base = 10e3; % VA， 容量基准 % 计算基准阻抗 Z_base = V_base^2 / S_base; % LCL滤波器参数 L1 = 1.5e-3; % H， 逆变器侧电感 L2 = 0.5e-3; % H， 电网侧电感 Cf = 30e-6; % F， 滤波电容 Rd = 1; % Ohm， 阻尼电阻 % 弱电网等效阻抗（模拟电网强度变化） % 强电网时Zg较小，弱电网时Zg较大 Lg_weak = 5e-3; % H， 弱电网等效电感 Rg_weak = 0.1; % Ohm， 弱电网等效电阻 Zg_weak = Rg_weak + 1j*2*pi*f_base*Lg_weak; % 控制参数示例（光伏逆变器） % 电流环PI参数 Kp_i = 0.5; Ki_i = 50; % 电压环PI参数 Kp_v = 0.1; Ki_v = 5; % 锁相环参数 Kp_pll = 50; Ki_pll = 500; % 控制参数示例（虚拟同步发电机） % 有功-频率下垂系数 Dp = 5; % Hz/pu % 无功-电压下垂系数 Dq = 0.05; % pu/pu % 虚拟惯量J和阻尼系数D J = 0.2; % kg·m^2 (标幺值) D = 10; % N·m·s/rad (标幺值)

3. 核心原理与阻抗建模方法

3.1 阻抗测量原理：扰动注入法

在仿真中获取阻抗频率特性最直接的方法是扰动注入法。其基本原理是在并网接口点（PCC点）注入一个幅值很小、频率可变的电压或电流扰动信号，同时测量该点的电压和电流响应。阻抗即为电压响应与电流响应的比值。

• 正序阻抗测量: 注入正序三相平衡的小信号电压扰动v_p^~，测量正序电流响应i_p^~，则正序阻抗Z_p = v_p^~ / i_p^~。
• 负序阻抗测量: 注入负序三相平衡的小信号电压扰动v_n^~，测量负序电流响应i_n^~，则负序阻抗Z_n = v_n^~ / i_n^~。

在Simulink中，可以通过受控电压源串联在PCC点来实现扰动注入。

3.2 光伏并网逆变器阻抗模型推导

光伏并网逆变器通常采用电网电压定向的矢量控制。在dq同步旋转坐标系下，其小信号模型可以线性化。考虑锁相环的动态，正序阻抗模型在dq域是一个2x2的传递函数矩阵。经过推导和简化（在基频附近），其正序阻抗幅值在低频段呈感性，主要受电流环带宽和滤波器影响；在高频段，主要由LCL滤波器的容性支路决定。负序阻抗模型与正序类似，但频率响应关于零频率对称。其数学模型复杂，本文不展开全部推导，仿真中将通过测量直接获取。

3.3 虚拟同步发电机阻抗模型特点

虚拟同步发电机的核心是模仿同步发电机的转子运动方程和励磁调节器。其输出阻抗特性与传统同步机有相似之处：

• 正序阻抗: 在低频段，由于虚拟惯量J和阻尼D的存在，其阻抗特性表现为一个二阶系统，能够为系统提供惯性支撑和阻尼，有助于抑制低频振荡。
• 负序阻抗: VSG的控制通常基于正序分量，对负序分量的调节能力有限，其负序阻抗通常较小，且特性与滤波器参数关系更密切。 VSG的阻抗模型同样复杂，涉及功率计算、下垂控制、虚拟转子等多个环节的耦合。仿真测量是分析其阻抗特性的有效手段。

4. Simulink仿真实现与阻抗扫描

4.1 光伏逆变器模型搭建

1. 主电路搭建:

   • 使用DC Voltage Source模拟光伏阵列输出的直流电压（如800V）。
   • 使用Universal Bridge模块作为三相逆变器，器件选择IGBT/Diodes。
   • 按照L1-Cf-(L2+Rd)的结构搭建LCL滤波器。
   • 使用Three-Phase Series RLC Branch模拟弱电网阻抗Zg_weak。
   • 使用Three-Phase Programmable Voltage Source作为理想电网电压源。

2. 控制系统搭建:

   • 锁相环: 使用PLL (3ph)模块或自行搭建基于dq变换的PLL。
   • 坐标变换: 使用abc to dq0和dq0 to abc模块。
   • 双闭环控制:
     • 外环（电压/功率）: 根据直流电压指令或功率指令生成d轴电流参考值id_ref。q轴电流参考值iq_ref通常设为0以实现单位功率因数。
     • 内环（电流）: 采用PI调节器对id和iq进行跟踪控制，输出dq轴电压调制信号。
   • PWM生成: 将dq轴电压调制信号反变换到abc坐标系，与三角载波比较生成PWM脉冲驱动逆变桥。

   % 文件：PV_Controller.slx 中的部分初始化脚本（嵌入在Model Properties -> InitFcn中） % 设置PI控制器参数 Kp_i = 0.5; Ki_i = 50; Kp_v = 0.1; Ki_v = 5; % 设置PLL参数 Kp_pll = 50; Ki_pll = 500; % 设置载波频率 f_sw = 10e3; % 10kHz

4.2 虚拟同步发电机模型搭建

1. 主电路部分与光伏逆变器模型基本相同。
2. VSG核心控制算法:
   • 有功-频率控制:

     % 模拟转子运动方程 % dw/dt = (Pm - Pe - D*w) / (2H) ， 其中H=J*w0/2, w为角速度偏差 % 在Simulink中可用积分器、增益和加减法器实现

   • 无功-电压控制:

     % U_ref = U_set + Dq * (Q_set - Q) % 其中U_ref为输出电压幅值参考，Q为测量无功功率

   • 电压电流生成: 根据计算出的频率f（决定相位角theta）和电压幅值U_ref，生成三相电压参考信号vabc_ref。然后通过电流环跟踪这个电压参考，从而模拟同步机的内电势特性。

4.3 阻抗扫描仿真步骤

1. 搭建扰动注入回路: 在PCC点与电网阻抗之间，插入一个Controlled Voltage Source。其输入信号为一个小幅值的正弦扰动，频率f_dist从1Hz到2000Hz扫描。
2. 配置仿真求解器: 使用变步长ode23tb或固定步长离散求解器，步长设置需小于扰动信号最高频率周期的1/10。
3. 编写自动化扫描脚本:

   % 文件：impedance_scan.m clear; close all; clc; load_system('PV_Inverter_Model.slx'); % 或 'VSG_Model.slx' f_dist_list = logspace(0, log10(2000), 100); % 1Hz到2000Hz， 取100个对数点 Zp_measured = zeros(size(f_dist_list)); Zn_measured = zeros(size(f_dist_list)); for idx = 1:length(f_dist_list) f_dist = f_dist_list(idx); % 1. 设置模型工作区中的扰动频率参数 set_param('PV_Inverter_Model/扰动源', 'Frequency', num2str(f_dist)); % 2. 运行仿真 simOut = sim('PV_Inverter_Model', 'StopTime', '1'); % 仿真1秒，待稳态 % 3. 获取PCC点电压、电流数据 t = simOut.logsout.get('t').Values.Data; v_abc = simOut.logsout.get('v_pcc').Values.Data; i_abc = simOut.logsout.get('i_pcc').Values.Data; % 4. 提取稳态周期数据进行FFT，计算正负序分量 [Vp, Vn] = calc_seq_components(v_abc(end-2000:end， :), f_dist, 1/f_base); [Ip, In] = calc_seq_components(i_abc(end-2000:end， :), f_dist, 1/f_base); % 5. 计算阻抗 Zp_measured(idx) = Vp / Ip; Zn_measured(idx) = Vn / In; end % 绘制阻抗伯德图 figure; subplot(2,1,1); semilogx(f_dist_list, 20*log10(abs(Zp_measured)), 'b-', 'LineWidth', 1.5); hold on; semilogx(f_dist_list, 20*log10(abs(Zn_measured)), 'r--', 'LineWidth', 1.5); legend('正序阻抗|Zp|', '负序阻抗|Zn|'); ylabel('幅值 (dB)'); grid on; title('光伏逆变器序阻抗频率特性'); subplot(2,1,2); semilogx(f_dist_list, angle(Zp_measured)*180/pi, 'b-', 'LineWidth', 1.5); hold on; semilogx(f_dist_list, angle(Zn_measured)*180/pi, 'r--', 'LineWidth', 1.5); legend('正序阻抗相位', '负序阻抗相位'); xlabel('频率 (Hz)'); ylabel('相位 (度)'); grid on;

   注：calc_seq_components是一个需要自定义的函数，用于从三相时域信号中提取指定频率的正负序相量。

5. 稳定性对比分析与时域验证

5.1 基于阻抗比的稳定性判据（奈奎斯特判据）

对于并网系统，稳定性可以通过比较源侧阻抗Zs（电网阻抗）和负载侧阻抗Zl（逆变器阻抗）来判断。常用判据是阻抗比Zs/Zl的奈奎斯特曲线是否包围(-1, j0)点。更直观的方法是观察正序阻抗的幅相特性与电网阻抗的交点（即幅值相等处）的相位差。

• 稳定条件: 在|Zl| = |Zs|的频率点，相位差应大于-180度（即相位裕度为正）。
• 对比分析: 在同一弱电网阻抗Zg_weak下，分别绘制光伏逆变器和VSG的正序阻抗Zl_pv和Zl_vsg。比较两者与Zg_weak的幅相曲线，可以定性判断哪种拓扑在相同条件下相位裕度更大，即更稳定。

5.2 仿真结果对比示例

运行阻抗扫描脚本后，可能得到如下典型结论（具体形状与参数紧密相关）：

1. 光伏逆变器: 其正序阻抗在低频段（如10-100Hz）可能呈现较大的负相位（容性特性），容易与感性的弱电网阻抗在幅值交点处产生较小的相位裕度，甚至小于-180度，从而引发低频振荡（如10-30Hz）。
2. 虚拟同步发电机: 由于其控制引入了虚拟惯量，其正序阻抗在低频段（如0.1-10Hz）表现出类似同步机的感性特性，相位滞后较少。因此，在与相同弱电网阻抗的交点处，通常能获得更大的相位裕度，稳定性更好。

5.3 时域仿真验证

为了验证频域分析的结论，需要在弱电网参数下进行时域阶跃响应或扰动测试。

1. 测试场景: 设置电网阻抗为Lg_weak和Rg_weak。系统稳定运行后，在t=0.5秒时，令光伏逆变器或VSG的有功功率参考值阶跃增加20%。
2. 观察变量: PCC点电压、输出电流、有功功率、频率（VSG）。
3. 预期结果:
   • 不稳定的系统（相位裕度很小或为负）会出现明显的功率或电流振荡，且振荡可能发散或持续。
   • 稳定的系统在经历一段动态调节后，会平稳到达新的稳态点。VSG由于有惯量支撑，功率变化会更平滑，频率跌落更小。

   % 在Simulink中设置阶跃信号模块，从0.8pu阶跃到1.0pu。 % 通过Scope观察PCC三相电流波形。 % 稳定案例：电流波形在扰动后经过几个周期的衰减恢复正弦。 % 不稳定案例：电流波形出现幅值不断增大的振荡。

6. 常见问题与排查思路

在搭建模型和仿真过程中，可能会遇到以下典型问题：

7. 工程实践建议与扩展方向

7.1 参数设计与优化原则

• 光伏逆变器: 在弱电网下，应遵循“低带宽、强鲁棒”的设计原则。优先保证锁相环和电流环在电网阻抗变化范围内的稳定性，而非追求最快的动态响应。可以采用基于电网阻抗在线估计的自适应控制。
• 虚拟同步发电机: 虚拟惯量J和阻尼D的选取需在惯性支撑效果与动态响应速度之间折衷。J越大，惯性越强，但频率恢复越慢。D影响振荡的衰减速度。需要通过时域仿真和根轨迹分析来优化。

7.2 仿真到实践的桥梁

1. 模型保真度: 本文仿真为理想开关模型。实际应用中，需考虑数字控制延迟（计算延迟、PWM更新延迟）、采样噪声、死区效应等非理想因素，这些都会影响高频段的阻抗特性。
2. 稳定性增强措施:
   • 有源阻尼: 在LCL滤波器控制中引入电容电流或逆变器侧电流反馈，主动抑制谐振峰。
   • 阻抗重塑: 通过在控制环路中增加特定的补偿环节，主动改变逆变器输出阻抗的形状，使其在关键频段避开与电网阻抗的 adverse interaction。
   • 宽频带同步技术: 采用基于二阶广义积分器的锁频环等，提升在背景谐波和不对称电压下的同步性能。

7.3 扩展研究方向

1. 多机并联稳定性: 研究多台光伏逆变器或VSG并联运行时，彼此之间的阻抗交互以及对电网稳定性的影响。
2. 序阻抗不对称性影响: 深入研究负序阻抗特性对电网电压不平衡条件下设备运行性能（如损耗、电流畸变）的影响。
3. 硬件在环验证: 将本文的Simulink控制模型编译下载到实际DSP或FPGA控制器中，通过硬件在环实验平台进行验证，更贴近工程实际。
4. 与储能结合: 研究将VSG控制应用于光储一体机或储能变流器，平抑功率波动，提供更强的电网支撑能力。

通过本文的梳理，你应该已经掌握了在Simulink中建立光伏逆变器与虚拟同步发电机序阻抗模型，并进行弱电网稳定性对比研究的基本方法。从理论推导到仿真实现，再到结果分析，这一完整流程是电力电子与电力系统交互研究的基础。关键在于理解不同控制策略如何塑造设备的“外在性格”（阻抗特性），并学会利用频域和时域工具去评估系统稳定性。建议读者动手复现本文的仿真模型，通过调整参数观察系统行为的变化，这是深化理解最有效的途径。