不平衡电网下VSG控制与PR电流环优化方案

1. 不平衡电网下VSG控制的技术背景与挑战

随着新能源发电占比突破40%临界点，电网运行特性正在发生根本性变化。以光伏、风电为代表的分布式电源大规模接入，导致电网电压不平衡问题日益突出。根据IEEE 1547-2018标准统计，电压不平衡度超过2%的工况出现频率较传统电网提升了3倍以上。这种不平衡主要表现为三相电压幅值不对称和相位偏移，其根源在于：

• 单相负载的不均衡分布
• 不对称短路故障
• 分布式电源出力波动

传统并网逆变器多采用PQ控制策略，在理想电网条件下表现良好。但实际工程案例表明，当电网电压出现20%以上的不平衡时，传统控制方案会导致：

1. 并网电流谐波含量激增（THD>5%）
2. 有功/无功功率出现2倍频波动
3. 系统动态响应速度下降40%以上

虚拟同步发电机(VSG)技术通过模拟同步电机的转子运动方程，为系统提供虚拟惯量和阻尼特性。其核心优势在于：

• 通过机械方程模拟实现频率支撑
• 提供类似同步机的无功-电压调节特性
• 增强系统抗扰动能力

然而，常规VSG控制在三相不平衡工况下暴露明显缺陷：

1. 电流内环采用PI控制时，对负序分量和谐波抑制能力有限
2. 功率环动态响应速度与电流质量存在矛盾
3. 锁相环在电压跌落时易产生相位跳变

2. 系统整体控制架构设计

2.1 分层控制结构解析

本方案采用四级分层控制架构，各层级功能明确且相互协同：

code复制[功率计算层] → [VSG控制层] → [电压环] → [电流环]

2.1.1 功率计算模块实现

基于瞬时功率理论构建三相功率计算模型，具体实现步骤：

1. 采集三相电压v_abc和电流i_abc
2. 进行Clarke变换得到αβ坐标系分量：

   code复制v_α = (2v_a - v_b - v_c)/3
   v_β = (v_b - v_c)/√3
   

3. 采用双二阶广义积分器(DSOGI)进行正负序分离：
   • 设置积分器截止频率ω_c=10rad/s
   • 正序提取时间常数τ=20ms
4. 计算瞬时功率：

   code复制p = v_α·i_α + v_β·i_β
   q = v_β·i_α - v_α·i_β
   

关键参数选择：DSOGI的Q值设为1.414，可在动态响应与滤波效果间取得平衡。

2.1.2 改进型锁相环设计

采用双同步坐标系解耦锁相环(DDSRF-PLL)，其创新点在于：

• 正序dq坐标系实现快速锁相（带宽50Hz）
• 负序dq坐标系专门抑制电压跌落干扰
• 交叉解耦补偿消除正负序耦合效应

实测性能指标：

• 相位跟踪误差<0.2°（不平衡度25%时）
• 响应时间<20ms（阶跃扰动下）
• 谐波抑制比>40dB

2.2 VSG核心算法实现

2.2.1 转子运动方程离散化

将同步机转子方程离散为差分形式：

code复制Δω = (P_ref - P_meas - D·ω)/(J·s)
θ = ∫ωdt

其中关键参数设计：

• 虚拟惯量J=0.2kg·m²（对应时间常数H=3s）
• 阻尼系数D=15N·m·s/rad

2.2.2 功率环动态调节

采用自适应阻尼系数策略：

code复制D = D_base + k·|Δω|

式中k=0.5为调节系数，当频率偏差Δω增大时自动增强阻尼效果。

3. PR控制在电流环的应用细节

3.1 PR控制器原理剖析

比例谐振(PR)控制器的传递函数为：

code复制G_PR(s) = k_p + 2k_r·ω_c·s/(s²+2ω_c·s+ω_0²)

与PI控制相比具有以下优势：

1. 在基频ω_0处提供无穷大增益
2. 对特定频率谐波具有选择性抑制
3. 相位滞后小（<5°@50Hz）

3.2 参数整定方法

通过频域分析法确定最优参数：

1. 谐振增益k_r：根据谐波抑制需求设定
   • 基波k_r=50
   • 3次谐波k_r=20
   • 5次谐波k_r=10
2. 截止频率ω_c：影响带宽
   • 取ω_c=10rad/s（约1.6Hz）
3. 比例系数k_p：决定动态响应
   • 按相位裕度>45°要求，取k_p=0.5

3.3 数字实现注意事项

在离散化过程中需采用Tustin变换保持频率响应特性：

code复制s = (2/T)·(z-1)/(z+1)

采样周期T=100μs时，需注意：

• 避免数值溢出（采用32位浮点运算）
• 谐振频率微调补偿（49.8-50.2Hz可调）
• 抗积分饱和处理（clamp限幅）

4. Simulink建模关键技巧

4.1 模型搭建要点

1. 主电路部分：

   • 采用Average Model提高仿真速度
   • 线路阻抗设置为0.1+j0.5Ω
   • 滤波电感L=2mH（ESR=0.05Ω）

2. 控制部分：

   • 使用MATLAB Function模块实现DSOGI
   • 采用Triggered Subsystem实现多速率控制
   • 添加Rate Transition模块处理跨时钟域信号

4.2 故障模拟设置

电压跌落工况配置：

• 跌落相别：C相
• 跌落幅度：40%（187V/311V）
• 持续时间：1s（1-2s区间）
• 恢复斜率：100V/s

4.3 仿真步长选择

采用变步长求解器ode23tb，设置：

• 最大步长：50μs
• 相对容差：1e-4
• 绝对容差：1e-6

5. 实测性能与问题排查

5.1 典型运行数据对比

5.2 常见问题解决方案

问题1：启动时电流冲击

• 原因：初始相位未同步
• 对策：添加软启动模块，在0.5s内线性增加功率指令

问题2：高频振荡

• 现象：开关频率附近出现谐振
• 解决方法：
  1. 增加RC阻尼支路（R=10Ω, C=0.1μF）
  2. 降低电流环带宽至800Hz

问题3：稳态误差

• 排查步骤：
  1. 检查PR控制器谐振频率是否偏移
  2. 验证PLL锁相精度
  3. 测量实际线路阻抗与模型匹配度

6. 工程应用建议

1. 硬件选型指南：

   • DSP建议选用TI C2000系列（如TMS320F28379D）
   • 电流传感器带宽需>50kHz（如LEM LAH-50P）
   • ADC采样率不低于100kSPS

2. 参数现场调试流程：

   • 先调电压环（带宽50Hz）
   • 再调电流环（带宽1kHz）
   • 最后优化VSG惯性时间常数

3. 电磁兼容设计：

   • 交流侧加装共模扼流圈
   • 直流母线布置薄膜电容（2.2μF/kW）
   • 控制板采用四层板设计

在实际项目中，我们发现在电网阻抗较大（>1Ω）的场合，需要适当增大滤波电感值至3mH，同时调整PR控制器的谐振增益k_r提高至80，这样可以有效抑制因线路压降导致的电流畸变。另外，建议在功率指令通道添加一阶惯性环节（时间常数100ms），避免功率突变对机械轴的冲击。