微电网孤岛与并网模式切换的Matlab仿真实现

1. 项目概述

微电网作为分布式能源系统的重要组成部分，其孤岛与并网运行模式的平滑切换一直是电力系统领域的研究热点。这个基于Matlab2016a平台的仿真项目，主要针对含燃气轮机和双馈感应发电机(DFIG)的微电网系统，研究其在不同运行模式下的控制策略实现。

在实际工程中，微电网需要具备两种基本能力：一是孤岛运行时能够自主维持电压频率稳定；二是并网运行时可以实现功率的优化分配。本项目通过Simulink搭建的仿真模型，验证了燃气轮机采用主从控制策略作为系统主电源，配合DFIG的最大功率跟踪控制，能够有效实现上述功能需求。

2. 系统架构设计

2.1 整体结构框图

系统由以下核心组件构成：

• 燃气轮机发电机组（主控单元）
• DFIG风力发电系统（从属单元）
• 本地负载
• 静态开关（用于孤岛/并网切换）
• 中央控制器

各组件通过交流母线连接，燃气轮机采用同步发电机模型，通过调速器和励磁系统实现电压/频率调节；DFIG系统包含背靠背变流器和最大功率点跟踪(MPPT)算法。

2.2 主从控制策略设计

燃气轮机作为主控单元，其控制策略包含两个层级：

1. 一次调频控制：通过下垂特性实现功率-频率调节

   • 频率下垂系数设定为4%
   • 有功功率参考值根据负载需求动态调整

2. 二次调频控制：采用PI调节器消除稳态误差

   • 比例系数Kp=0.8
   • 积分时间常数Ti=15s

DFIG作为从属单元，主要实现：

• 转子侧变流器控制：维持定子电压稳定
• 网侧变流器控制：实现直流母线电压稳定
• MPPT控制：通过最优转矩法追踪最大功率点

3. Simulink建模实现

3.1 燃气轮机模型搭建

在Matlab2016a的Simulink环境中，燃气轮机模型包含以下关键模块：

1. 原动机模型：

   matlab复制% 燃气轮机传递函数模型
   num = [1];
   den = [0.04 1];
   G_turbine = tf(num,den);
   

2. 同步发电机参数设置：

   • 额定功率：500kW
   • 额定电压：400V
   • Xd' = 0.15 pu
   • 惯性时间常数H=3.5s

3. 励磁系统采用IEEE AC1A标准模型：

   • 电压调节器增益Ka=200
   • 励磁机时间常数Te=0.8s

3.2 DFIG控制系统实现

DFIG的Simulink模型重点实现了：

1. 转子侧变流器控制：

   • 采用双闭环控制结构（外环电压、内环电流）
   • 开关频率：5kHz
   • LCL滤波器参数：L=2mH, C=50μF

2. MPPT算法实现：

   matlab复制% 最优转矩法核心算法
   function T_ref = MPPT_OT(omega,P_prev)
       K_opt = 0.48; % 最优转矩系数
       T_ref = K_opt * omega^2;
       if P_prev < P_current
           K_opt = K_opt * 0.98; 
       end
   end
   

4. 模式切换控制策略

4.1 孤岛到并网切换

切换过程分为三个阶段：

1. 预同步阶段：

   • 检测电网电压幅值、频率、相位
   • 调节燃气轮机输出使差值<2%
   • 同步检查时间窗口：5个周期

2. 闭锁阶段：

   • 禁用DFIG的MPPT功能
   • 设置燃气轮机为PQ控制模式
   • 持续时间：100ms

3. 并网执行：

   • 静态开关闭合指令
   • 延时50ms后恢复MPPT

4.2 并网到孤岛切换

关键处理流程：

1. 电网故障检测：

   • 电压跌落>15%持续60ms
   • 频率偏差>0.5Hz

2. 控制模式切换：

   • 燃气轮机转为V/f控制
   • DFIG切换为电压支持模式
   • 变流器限流值设为1.2pu

3. 负载再分配：

   • 优先保障重要负载
   • 卸载非关键负载（如有必要）

5. 仿真结果分析

5.1 稳态性能

在额定负载条件下：

• 电压偏差：<1.2%
• 频率偏差：<0.15Hz
• 总谐波畸变率(THD)：<3%

5.2 动态响应

负载突变（50%→100%）时：

• 频率最大跌落：0.8Hz
• 恢复时间：1.2s
• 电压暂降：4.5%

5.3 模式切换特性

孤岛→并网切换：

• 冲击电流：<1.5倍额定
• 相位差：<5°
• 完成时间：280ms

6. 关键参数调试经验

6.1 下垂系数优化

通过实测发现：

• 频率下垂系数在3%-5%时系统稳定性最佳
• 功率下垂系数建议取2%-3%
• 不同运行模式下需采用自适应调整：

  matlab复制if island_mode
      R_f = 0.04;
  else
      R_f = 0.03;
  end
  

6.2 PI参数整定技巧

对于燃气轮机调速器：

1. 先设Ti=∞，调整Kp至系统出现轻微振荡
2. 取该Kp值的60%作为最终值
3. 逐步减小Ti直至动态响应达标

实测推荐值：

• 转速环：Kp=1.2, Ti=0.5s
• 功率环：Kp=0.8, Ti=1.2s

7. 常见问题解决方案

7.1 仿真发散问题

可能原因及处理：

1. 代数环问题：

   • 在反馈回路中加入单位延迟模块
   • 步长改为变步长ode23t

2. 初始值不匹配：

   • 使用Powergui的Load Flow初始化
   • 检查电机初始滑差设置

7.2 切换过程振荡

典型解决方法：

1. 增加预同步时间窗口
2. 在控制回路中加入加速度限制：

   matlab复制function omega_ref = rate_limiter(omega_cmd)
       persistent prev_omega;
       max_rate = 0.1; % pu/s
       omega_ref = prev_omega + sign(omega_cmd-prev_omega)*min(abs(omega_cmd-prev_omega),max_rate*Ts);
       prev_omega = omega_ref;
   end
   

8. 模型优化建议

8.1 实时性优化

对于大型模型：

1. 将连续系统改为离散系统

   • 采样时间取开关周期的1/10
   • 使用Fixed-step求解器

2. 分模块编译：

   matlab复制% 将DFIG子系统编译为S-function
   rtwbuild('DFIG_Controller');
   

8.2 精度提升方法

1. 采用平均值模型替代开关细节
2. 对关键变量使用double精度：

   matlab复制set_param('model/Block','DataType','double');
   

3. 启用代数环检测选项

在实际工程应用中，这种控制策略已经证明可以在85%以上的工况下实现无缝切换。不过需要注意，当微电网中含有大量电力电子设备时，可能需要额外考虑虚拟同步机(VSG)控制策略来增强系统惯性。