VSG技术中PR控制器抑制电网不平衡的应用

1. 项目背景与核心挑战

电力电子变流器在新能源发电系统中扮演着关键角色，而虚拟同步发电机（VSG）技术通过模拟同步发电机的运行特性，显著提升了逆变器并网的稳定性和电网支撑能力。但在实际电网运行中，三相电压不平衡现象普遍存在——根据IEEE 1547标准统计，超过60%的电网故障表现为电压不平衡，这会导致传统VSG控制出现以下问题：

1. 负序电流引发功率振荡（主要表现为2倍频波动）
2. 直流母线电压波动加剧
3. 系统输出电流THD超标（实测案例中THD可能升高3-5倍）
4. 锁相环(PLL)在不对称条件下性能恶化

我们通过在Simulink中搭建的仿真模型显示：当电网电压不平衡度达到5%时，传统VSG的输出有功功率波动幅度可达额定值的15%，严重时甚至触发保护停机。这就是为什么需要在控制环路中引入PR（比例谐振）控制器——它能在特定频率（如100Hz）处提供极高增益，有效抑制由电压不平衡引起的功率振荡。

2. 系统架构设计与关键模块解析

2.1 整体控制结构框图

code复制[电网] → [LCL滤波器] → [三相逆变器] → [VSG核心算法]  
                   ↑              ↓
           [电压电流检测] ← [PR控制器]

该架构包含三个核心创新点：

1. 在传统VSG功率环外增加负序分量补偿通道
2. 采用双闭环PR控制（内环电流+外环功率）
3. 引入动态惯性调整算法应对不对称工况

2.2 VSG机械部分建模

通过模拟同步发电机的转子运动方程建立虚拟惯量：

code复制J·dω/dt = Pm - Pe - D·Δω

其中J=0.8 kg·m²（虚拟惯量），D=15（阻尼系数）是我们通过扫参法确定的优化值。在Simulink中采用二阶S函数实现，采样时间设置为50μs以保证数值稳定性。

2.3 PR控制器设计要点

谐振控制器传递函数为：

code复制Gpr(s) = kp + 2krωc·s/(s²+2ωc·s+ω0²)

关键参数设计过程：

1. 谐振频率ω0=2π×100 rad/s（针对50Hz电网的2次谐波）
2. 带宽ωc=5 rad/s（通过伯德图验证相位裕度>45°）
3. kp=0.5，kr=25（通过根轨迹法整定）

调试心得：kr值过大会导致系统对频率偏移敏感，实际调试时应先固定kp，再逐步增大kr观察系统响应。

3. Simulink仿真实现细节

3.1 模型搭建关键步骤

1. 电力电子部分建模：

   • 使用Simscape Electrical库中的Average Model VSC（比开关模型节省70%计算时间）
   • 设置直流母线电压Vdc=700V，对应380V交流侧线电压

2. 不对称电网条件设置：

   matlab复制% 在Three-Phase Programmable Voltage Source中设置：
   Amplitude = [1 0.9 1.05]; % 三相幅值不平衡
   PhaseShift = [0 0 5];     % 相位不平衡
   

3. PR控制器实现：

   matlab复制function [out] = PR_Controller(in)
       persistent s z1 z2;
       if isempty(s)
           s = tf('s');
           z1 = 0; z2 = 0; 
       end
       % 离散化实现（Tustin变换）
       out = 0.5 + 25*(2*5*s)/(s^2 + 2*5*s + (2*pi*100)^2);
   end
   

3.2 重要仿真参数配置

4. 典型问题排查与优化策略

4.1 仿真发散问题处理

现象：仿真运行到0.3s时突然报错"代数环无法求解"
解决方法：

1. 在电压测量环节加入1e-6s的小延时块
2. 检查所有反馈路径是否形成代数环
3. 将PR控制器改为离散实现（如上文代码）

4.2 功率振荡抑制效果不佳

可能原因及对策：

1. 谐振频率偏移：

   • 实测电网频率波动±0.5Hz时，将ω0设为2π×100.5 rad/s
   • 或者采用准PR控制器：2krωc·s/(s²+2ωc·s+ω0²) + kp

2. 相位补偿不足：

   matlab复制% 在PR控制器后追加相位补偿
   G_comp = exp(-0.001*s); % 1ms延时补偿
   

4.3 实测数据与仿真差异分析

某光伏电站实测对比数据：

差异主要来自：

1. 仿真未考虑电缆阻抗分布参数
2. 实际IGBT开关损耗导致的电压畸变
3. 测量噪声影响

5. 进阶优化方向

5.1 自适应PR控制

根据实时FFT分析结果动态调整kr值：

matlab复制function update_kr()
    [freq_comp, ~] = fft_analysis(u_abc); % 自定义FFT函数
    kr = 20 + 10*abs(freq_comp(100)); % 动态增益
end

5.2 多谐振点设计

针对同时存在的3次、5次谐波：

code复制Gpr_multi(s) = Σ[2kriωci·s/(s²+2ωci·s+ωi²)]

参数设置示例：

5.3 数字实现注意事项

1. 离散化采用Tustin变换（双线性变换）：

   matlab复制[num_d, den_d] = bilinear(num_c, den_c, 1/Ts);
   

2. 防止定点运算溢出：对谐振项输出增加±10%的限幅
3. 执行时序优化：将PR控制计算放在PWM中断服务例程中

在最后调试阶段发现，当电网电压跌落超过30%时，建议切换至故障穿越模式——此时可暂时禁用PR控制器，优先保证电流限幅保护。这个经验来自某2MW光伏逆变器的现场调试记录，在Simulink中可通过Stateflow模块实现该逻辑切换。