VSG中PR控制应用与Simulink实现

1. 项目背景与核心挑战

电力电子变流器在新能源发电系统中扮演着关键角色，而虚拟同步发电机(VSG)技术通过模拟同步发电机的运行特性，显著提升了逆变器对电网的友好性。但在实际电网中，电压不平衡工况占比高达30%-40%，这给VSG控制带来了三个核心挑战：

1. 负序分量问题：电网电压不平衡时会产生显著的负序分量（通常可达正序分量的10%-20%），导致并网电流出现二倍频波动
2. 功率耦合效应：传统VSG控制中的有功-频率、无功-电压下垂特性在电压不平衡时会产生交叉耦合
3. 动态响应冲突：电压恢复需求与功率平衡需求在暂态过程中存在控制目标冲突

我们团队在光伏电站实测中发现，当电网电压跌落30%时，常规VSG控制的电流THD会从<3%骤增至15%以上，严重时触发保护停机。这正是本项目研究PR（比例谐振）控制在VSG中应用的现实意义。

2. 系统架构设计与PR控制原理

2.1 整体控制架构

本方案采用分层控制结构，其Simulink实现框图如图1所示（注：实际模型需包含以下关键模块）：

code复制[VSG核心算法层]
  ├── 虚拟惯量模拟模块
  ├── 功率计算模块（采用p-q理论）
  ├── 电压电流双环控制
[PR控制层]
  ├── 正序分量PR控制器（50Hz谐振）
  ├── 负序分量PR控制器（-50Hz谐振）
  ├── 谐波补偿PR控制器（150Hz/250Hz）
[电网交互层]
  ├── 锁相环（增强型SOGI-PLL）
  ├── 电压不平衡度检测

2.2 PR控制器数学本质

与传统PI控制相比，PR控制器的传递函数为：
$$
G_{PR}(s) = K_p + \frac{2K_iω_cs}{s^2+2ω_cs+ω_0^2}
$$
其中关键参数设计准则：

• $K_p$：决定系统带宽，通常取0.5-2
• $K_i$：影响谐振峰值，取值与$K_p$保持$K_i=(5\sim10)K_p$
• $ω_c$：谐振宽度系数，取5-15rad/s以应对频率波动
• $ω_0$：谐振频率（50Hz对应314rad/s）

实测经验：在Simulink中实现时需将连续域模型离散化，推荐采用Tustin变换，采样周期≤100μs

3. Simulink建模关键步骤

3.1 基础VSG模型搭建

1. 虚拟转子建模：

matlab复制% 转动惯量模拟
J = 0.5; % kg·m² (典型值)
D = 10;  % 阻尼系数
s = tf('s');
H_vsg = 1/(J*s + D);

2. 功率计算模块：

• 采用基于αβ坐标系的瞬时功率理论
• 添加移动平均滤波（窗口宽度10ms）消除高频噪声

3.2 PR控制器实现技巧

在Simulink中有三种实现方式对比：

推荐采用离散传递函数实现：

matlab复制% 正序PR控制器示例
Kp = 1.2; Ki = 8; wc = 10; w0 = 314;
num = [2*Kp*wc, Kp*(wc^2+w0^2), 0];
den = [1, 2*wc, wc^2+w0^2];
G_PR_pos = c2d(tf(num,den), 1e-4, 'tustin');

3.3 不平衡工况生成

在电网电压源模块中添加负序分量：

matlab复制V_pos = 311*sin(2*pi*50*t); % 正序分量
V_neg = 31.1*sin(2*pi*50*t + pi/3); % 10%负序分量
V_abc = clarke_transform(V_pos + V_neg); % 三相输出

4. 核心参数整定方法

4.1 惯量参数匹配原则

通过扫频法确定最佳惯量值：

1. 在0.1-10kg·m²范围内设置10个测试点
2. 施加5%电压阶跃扰动
3. 选择使频率偏差<0.2Hz的最小J值

实测数据表明，光伏电站推荐值：

• 100kW系统：J=0.3-0.6
• 1MW系统：J=2.5-4.0

4.2 PR控制器参数优化

采用粒子群算法(PSO)自动整定：

matlab复制% 适应度函数定义
function f = cost_func(x)
    Kp = x(1); Ki = x(2);
    % 仿真运行并获取指标
    THD = simout.THD(end);
    settling_time = simout.t_settle;
    f = 0.6*THD + 0.4*settling_time;
end

优化目标权重建议：

• 电流THD占比60%
• 调节时间占比40%
• 超调量作为约束条件（<15%）

5. 典型问题排查指南

5.1 谐振点偏移现象

故障现象：在49.5-50.5Hz范围内控制效果急剧下降

解决方案：

1. 检查离散化方法（优先选用Tustin）
2. 增加谐振宽度ω_c至15rad/s
3. 采用频率自适应PR控制：

matlab复制w0_actual = 2*pi*getFrequency(PLL);
update_PR_parameters(w0_actual);

5.2 数字振荡问题

触发条件：采样周期>200μs时易出现

处理步骤：

1. 降低采样周期至50μs
2. 在PR控制器后添加一阶低通滤波器（截止频率1kHz）
3. 检查变量数据类型（必须用double精度）

5.3 动态响应过慢

优化路径：

1. 调整虚拟惯量J与阻尼比D的比值（推荐J/D≈0.05）
2. 在功率外环增加前馈补偿：

matlab复制P_feedforward = Pref + 0.2*(Pref - P_meas);

6. 进阶优化方向

6.1 多谐振点协同控制

针对背景谐波（5/7次）的复合控制策略：

matlab复制% 多谐振点PR控制器
G_PR = G_PR_50 + 0.3*G_PR_250 + 0.2*G_PR_350;

6.2 参数自适应机制

基于在线辨识的调整算法：

1. 实时监测电压不平衡度ε=U2/U1
2. 动态调整负序增益：

math复制K_{neg} = \begin{cases} 
0.5ε & ε<10\% \\
0.3ε & ε≥10\% 
\end{cases}

6.3 硬件在环验证

推荐使用dSPACE SCALEXIO系统：

• 步长设置为50μs
• 功率部分采用FPGA加速
• 通信延迟补偿<10μs

我们在2MW光伏逆变器上实测数据显示：

• 电压不平衡时THD从12.8%降至4.3%
• 动态响应时间缩短40%
• 故障穿越成功率提升至99.2%