Simulink环境下的Active AFD主动反孤岛模型仿真

一、理论基础

l一个性能完善的光伏并网发电系统,需要各种保护措施保证用户的人身安全,同时防止设备因意外而造成的损坏。由于光伏发电系统和电网并联工作,因此光伏发电系统需能及时检测出电网故障并切断其与电网的连接。如果不能及时发现电网故障,就会出现光伏发电系统仍向局部电网供电的情形,从而使本地负载仍处于供电状态,造成设备损坏和人员伤亡。这种现象被称为孤岛效应。

目前孤岛检测主要分为被动式检测和主动式检测两种。被动式检测是利用电网断电瞬间,逆变器输出功率与局部负载功率不平衡造成的逆变器输出端电压值和频率变化作为孤岛效应检测的依据。其具有检测方法简单,对系统运行无干扰等优点。但是如果在电网停电的瞬间,逆变器输出功率与局部负载功率达到平衡,该方法将失去作用。

主动式检测是在系统工作中,对逆变器输出电流、频率或相位施加一定的扰动信号,并对其进行检测。如果电网正常,因电网的巨大平衡作用,逆变器输出不受扰动信号的影响;一旦电网出现故障,这些扰动量就会在逆变器输出端逐步累计,直至超出规定范围,从而反映出电网故障。与被动式检测方法相比,主动式检测方法具有精度高,检测盲区小的优点。但当局部电网存在多个分布式能源系统时,主动式检测效果下降,严重时甚至无效。

被动式检测主要有电压频率检测、相位检测、频率检测三种方法,在实际中均有一定的应用。但是由于被动式孤岛检测方法对逆变器输出功率与负载功率是否匹配有较高的要求,因此存在较大的检测盲区。所以在此不做详细的描述。

其中,主动频率偏移(Active Frequency Drift.AFD)孤岛检测方法的思想是通过锁相环(PLL)确定公共点电压的频率和相位.调整光伏逆变器输出电流的给定频率。使电流频率比公共点电压频率略高或略低。若电流半波已完成而电压过零点未到。则强制电流给定为零。直到电压过零点触发到来。电流才开始下一个半波。


AFD的稳态频率应使负载相角满足如下关系

基于simulink的Active anti-islanding-AFD主动反孤岛模型仿真_simulink

式中tz为电流过零点超前(或滞后)电压过零点的时间间隔。 cf为AFD的截断系数,T为电压周期。

二、案例背景

1.问题描述

首先,这个电力系统,我们做如下的设置,即:

230V/50Hz,单相(single phase or 1 phase)

   然后,仿真的具体要求为:

负载:小型电网 single phase 1kW grid-connected PV systems

根据 IEEE 1547 standard 对反孤岛的要求

·Voltage: 240V ± 10%,

·frequency:50Hz ± 1%

·Quality Factor, Qf:  Qf < 2.5 (Qf = 0 , 0.02,0.04, … , 2.4)

·Required Islanding detection time: tdetect<0.2s  and isolated to the load

·THD: <5%

2.思路流程

整个反孤岛模型,采用simulink进行建模,在matlab2010b中,建立仿真模型,其中AFD算法则采用S函数进行设计,并封装为AFD模块供simulink调用。



三、部分MATLAB程序

AFD的S函数部分程序如下:

登录后复制

function [sys,x0,str,ts]=mdlInitializeSizes%定义全局变量global f_i;            global f_vo;global f_v_hb;global f_v_lb;global theta_i;global theta_vo;global isIslanding;global num0;global num1;global step;global step1;global k;%以下指标根据IEEE1547来取值f_i        = 50; %电流频频率f_vo       = 50; %电压频率theta_i    = 0;  %电流相位theta_vo   = 0;  %电压相位isIslanding= 0;  %非孤岛效应num0       = 0;num1       = 0;k          = 10;step1      = 0.5;f_v_hb     = 50;f_v_lb     = 50;%获得系统默认的系统参数变量sizes    sizes                = simsizes;%连续状态的个数sizes.NumContStates  = 0;%离散状态的个数sizes.NumDiscStates  = 0;%输出变量的个数sizes.NumOutputs     = 2;%输入变量的个数sizes.NumInputs      = 2;%布尔变量,表示有无直接馈入。1表示有,0便是没有sizes.DirFeedthrough = 1;%采样时间个数,s函数支持多采样系统sizes.NumSampleTimes = 1;  %将结构体sizes赋值给syssys                  = simsizes(sizes);%初始变量状态x0                   = [];%系统保留值,必须为空str                  = [];%采样周期变量ts                   = [1e-4 0]; function sys=mdlOutputs(t,x,u)%定义全局变量global f_i;            global f_vo;global f_v_hb;global f_v_lb;global theta_i;global theta_vo;global isIslanding;global num0;global num1;global step;global step1;global k;%锁相环输出theta_vo=u(2);%电压相位为0时更新频率if abs(theta_vo) < 0.04    %数字锁相环输出的电压频率  f_vo=u(1);    end%判断是否不是孤岛效应if(isIslanding==0)      %并网电压相位是否过零      if abs(theta_vo)<0.04        %判断并网电压频率是否越界,根据1547规定,1%的调整,所以就是正负0.5Hz          if (f_vo>50.5) || (f_vo<49.5)                 sys         = [0 0];                %如果超出则为孤岛检测到了                isIslanding = 1;          else                       f_vo = u(1);              if(f_v_lb < f_vo)                  num0=0;              end              if(f_vo   < f_v_hb)                 num1=0;                 end              if (f_vo>=50)                    num1=num1+1;                    if(num1>0)                      step=k*step1;                    else                      step=step1;                    end                       f_i    = f_vo+step;                     f_v_hb = f_i;              else                    num0=num0+1;                    if(num0>0)                      step=k*step1;                    else                      step=step1;                    end                        f_i    = f_vo-step;                    f_v_lb = f_i;              end                theta_i=theta_vo;          end    else        if(pi-theta_i<0.03) && (pi-theta_vo>0.03)             theta_i=pi;        elseif(2*pi-theta_i<0.03) && (2*pi-theta_vo>0.03)             theta_i=2*pi;                            elseif (pi-theta_vo<0.03) && (pi-theta_i>0.03)            theta_i=pi;        elseif (2*pi-theta_vo<0.03) && (2*pi-theta_i>0.03)            theta_i=2*pi;        else            theta_i=theta_i+2*pi*f_i*1e-4;        end    end    %非孤岛输出正弦波    sys(1)=sin(theta_i);else    %孤岛输出零    sys(1)=0;endsys(2)=f_vo;1.2.3.4.5.6.7.8.9.10.11.12.13.14.15.16.17.18.19.20.21.22.23.24.25.26.27.28.29.30.31.32.33.34.35.36.37.38.39.40.41.42.43.44.45.46.47.48.49.50.51.52.53.54.55.56.57.58.59.60.61.62.63.64.65.66.67.68.69.70.71.72.73.74.75.76.77.78.79.80.81.82.83.84.85.86.87.88.89.90.91.92.93.94.95.96.97.98.99.100.101.102.103.104.105.106.107.108.109.110.111.112.113.114.115.116.117.118.119.120.121.122.123.124.125.126.127.128.129.130.131.132.133.134.135.

AFD的simulink整体模型如下:

基于simulink的Active anti-islanding-AFD主动反孤岛模型仿真_AFD_02

四、仿真结论分析

   运行simulink仿真的结果进行显示,具体的仿真结果如下所示:

基于simulink的Active anti-islanding-AFD主动反孤岛模型仿真_matlab_03

基于simulink的Active anti-islanding-AFD主动反孤岛模型仿真_主动反孤岛_04

基于simulink的Active anti-islanding-AFD主动反孤岛模型仿真_主动反孤岛_05

      从前面的仿真结果可知,采用AFD这种主动反孤岛的方式,系统能够较快的检测到孤岛,但是同时AFD对系统的波形影响较大。

免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删

QR Code
微信扫一扫,欢迎咨询~

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 155-2731-8020
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

手机不正确

公司不为空