基于虚拟同步发电机的微网逆变器均衡控制方法
技术领域
本发明涉及一种微网控制方法,尤其是一种基于虚拟同步发电机的微网逆变器均衡控制方法。
背景技术
近年来,随着新能源发电单元在电力系统的渗透率不断提升,与此同时传统集中式一次能源逐渐减少,系统的转动惯量逐渐减小,频率波动变大,且一次能源的间歇性特性更加剧了电网的频率波动,使得系统的频率稳定性问题日趋严峻。且随着分布式电源电网渗透率的提高,其在电网中的角色将发生变化,分布式电源将不再只是向电网提供电能,还应能具备一定的电网电压和频率支撑能力,以稳定电网运行。然而目前基于电流源并网控制方式的分布式电源并不具备电压支撑能力,且并入不同电网结构时其稳定性和动稳态性能有较大差异,电网适应性较差。
传统电力系统中,同步发电机组(Generator Set-Genset)的下垂特性以及转动惯量大等因素,在维持系统的电压和频率稳定方面起着关键作用。Genset平稳和调节系统频率的过程可以分为三个阶段:第一阶段为Genset的惯性稳频,即依靠Genset自身转动惯量抑制系统的快速频率波动;第二阶段为一次调频,即当频率波动量超出一定值通过改变原动机功率输入来调节频率;第三阶段为二次调频,即而当系统功率恢复平衡后,调整一次调频指令将频率控制在额定频率值,从而实现频率的无差控制。显然,分布式发电系统中以电力电子装置实现的分布式发电单元(以下简称DGU)若能模拟或者部分模拟Genset的上述特性,使其像Genset一样参与频率和电压的调节过程,就可以降低分布式电源对电网的不利影响,解决分布式电源大规模并网应用中的相关技术瓶颈问题。而能模拟或者部分模拟Genset频率电压控制特性的电力电子电源装置就被称为虚拟同步发电机(VirtualSynchronous Generator,VSG)。VSG需要运行在两种模式下,并网和孤岛并联运行。
VSG并网运行时,需要对电网的电压和频率稳定性进行一定的支撑,并在电网频率异常的情况下做出限功率运行,孤岛并联运行时,需要向负载提供较高的电能质量。另外VSG需要运行在并网和孤岛两种模式下,当发生模式转换时,VSG应具有无缝切换能力。
针对上述问题,国内外的专家学者们提出了一些方法,主要有:
题为“具有转动惯量和阻尼自趋优的超越虚拟同步发电机方法”的中国发明专利申请说明书(CN105186554A)给出了一种根据转动惯量和阻尼自我调节来对电网频率差进行补偿控制,然而阻尼系数对稳态状态下的下垂特性有影响,控制上存在耦合,不利于参数设计。
题为“应用于虚拟同步发电机的离并网控制方法及系统”的中国发明专利申请说明书(CN105207261A)公开的技术方案中,采用快速终端滑模控制策略对离/并网状态进行切换,可以使得电网电压与负载电压的误差值快速平稳趋近于零,但是没有给出电网电压和负载电压之间的频率和相位差控制问题。
题为“一种基于虚拟同步发电机的解耦控制方法及装置”的中国发明专利申请说明书(CN105915140A)给出了一种电压电流的dq解耦控制方法,可以实现解耦,但是解耦特性依赖于状态变量的估计特性,且控制方法复杂。
总之,现有VSG技术并网模式下阻尼特性与下垂特性不能实现解耦控制,离网并联模式下很难同时兼顾动态响应与稳态均流之间的关系,并离网切换存在过渡过程。
发明内容
本发明要解决的技术问题为克服上述各种技术方案的局限性,针对VSG技术并网模式下的解耦控制与离网模式下的均衡控制,以及并离网模式切换下的过渡过程等问题,提供一种基于虚拟同步发电机的微网逆变器均衡控制方法。
本发明的目的是这样实现的。本发明提供了一种基于虚拟同步发电机的微网逆变器均衡控制方法,主要步骤如下:
步骤1,采样及坐标变换;
所述采样包括采集以下数据:微网逆变器滤波电容电压uca,ucb,ucc,微网逆变器桥臂侧电感电流iLa,iLb,iLc,微网逆变器并网点电网电压ea,eb,ec;
所述坐标变换包括对以下数据进行坐标变换:对微网逆变器滤波电容电压uca,ucb,ucc和桥臂侧电感电流iLa,iLb,iLc分别进行单同步旋转坐标变换得到滤波电容电压的dq分量Ucd,Ucq和桥臂侧电感电流的dq分量ILd,ILq;
步骤2,根据步骤1中得到的滤波电容电压的dq分量Ucd,Ucq,通过通用的微分离散化方程计算滤波电容电流的dq分量Icd,Icq;根据步骤1得到的桥臂侧电感电流的dq分量ILd,ILq和滤波电容电流的dq分量Icd,Icq,经过输出电流计算方程得到输出电流的dq分量Iod,Ioq;经过有功功率计算方程和无功功率计算方程得到平均有功功率P和平均无功功率Q;对微网逆变器并网点电网电压ea,eb,ec经过锁相环环节得到电网角频率ωg;
步骤3,根据步骤2中得到的平均有功功率P、电网角频率ωg和微网逆变器给定的有功功率指令P0、微网逆变器给定有功功率指令P0时的额定角频率ω0,经过功角控制方程及电网频率逻辑判断环节得到虚拟同步发电机的角频率ω,对ω积分得到虚拟同步发电机的矢量角θ;
步骤4,根据步骤2中得到的平均无功功率Q和微网逆变器给定的无功功率指令Q0、电压指令U0,经过无功控制方程得到虚拟同步发电机的端电压U*;
步骤5,先根据步骤4中得到的端电压U*和步骤1中得到的滤波电容电压dq的分量Ucd,Ucq,通过电压控制方程得到电流指令信号再根据电流指令信号步骤1中的桥臂侧电感电流的dq分量ILd,ILq和步骤2得到的滤波电容电流的dq分量Icd,Icq,通过加权电流控制方程得到控制信号Ud,Uq;
步骤6,将步骤5中得到的控制信号Ud,Uq经过单同步旋转坐标反变换得到三相桥臂电压控制信号Ua,Ub,Uc,再根据Ua,Ub,Uc生成开关管的PWM控制信号。
优选地,步骤2中所述平均有功功率P和平均无功功率Q的计算步骤包括:
步骤2.1,计算滤波电容电流的dq分量Icd,Icq;
令滤波电容电压Ucd,Ucq的离散序列为Ucd(n),Ucq(n),滤波电容电流dq分量Icd,Icq的离散序列为Icd(n),Icq(n),则计算滤波电容电流的通用的微分离散化方程为:
其中,C为滤波电容,Ts为微网逆变器采样频率,K为离散序列点数,n,k为自然数,即n=0,1,2,3,4......,k=0,1,2,3,4......;
根据上述方程可以求得滤波电容电流Icd,Icq的离散序列为Icd(n),Icq(n),从而可得滤波电容电流的dq分量Icd,Icq;
步骤2.2,计算输出电流的dq分量Iod,Ioq;
根据步骤2.1得到的滤波电容电流的dq分量Icd,Icq,经过输出电流计算方程得到输出电流的dq分量Iod,Ioq,所述的输出电流计算方程为:
Iod=ILd-Icd
Ioq=ILq-Icq
步骤2.3,根据有功功率计算方程和无功功率计算方程计算平均有功功率P和平均无功功率Q;
有功功率计算方程为:
无功功率计算方程为:
其中,Qpq为功率计算方程品质因数,ωh为陷波器需要滤除的谐波角频率,s为拉普拉斯算子,τ为一阶低通滤波器的时间常数,h为待抑制的谐波次数。
优选地,步骤3中所述功角控制方程为:
其中,ω0为微网逆变器给定有功功率指令P0时的额定角频率,m为功角控制下垂系数,J为模拟同步发电机机组的虚拟转动惯量,s为拉普拉斯算子,D1为微网逆变器频率反馈系数,D2为电网频率反馈系数;
采用具有回差特性的电网频率逻辑判断环节为:
当ωg≥ωg1时,令D1=-D2<0,即:
当ωg<ωg2时,令ω0=ωg,D1=-D2<0,即:
其中,ωg1为电网频率低阈值,ωg2为电网频率高返回值,ωg1<ωg2。
优选地,步骤4中所述无功控制方程为:
U*=U0+nq(Q0-Q)
其中,U0为微网逆变器给定无功功率指令Q0时的额定输出电容电压、nq为无功-电压下垂系数。
优选地,步骤5中所述电压控制方程为:
其中,Kp为电压环比例控制系数、Ki为电压环积分控制系数、Kr为电压环谐振控制器比例系数,Qu为电压环准谐振调节器品质因数,ωh为陷波器需要滤除的谐波角频率,s为拉普拉斯算子,h为待抑制的谐波次数。
优选地,步骤5中所述电流控制方程为:
其中,Kpi为电流环比例控制系数,Kri电流环谐振控制器比例系数,w1为电感电流的权重系数,w2为电容电流的权重系数,Kf为电压前馈系数,Qi为电流环准谐振调节器品质因数,s为拉普拉斯算子。
采用本发明后,对于采用虚拟同步发电机技术的微网逆变器,具备了如下优点:
1、电网频率低的异常情况下可以实现限功率运行,且只需要改变外环功率指令给定方式,而需要改变电压双环控制器,减小了过渡过程的冲击;模式切换运行时只需改变外环功率指令给定方式,不需要改变电压双环控制器,减小了过渡过程的冲击,减小了无缝切换时间。
2、虚拟阻尼不影响稳态下垂均分特性,与下垂特性分离控制与设计,相互解耦,提高了系统性能
3、采用基于电容电流与桥臂侧电感电流加权控制方案,在输出电压动态性能与输出阻抗之间取得均衡,有利于总体性能的提高。
4、抑制直流分量,且并网运行模式下具有较低的电流谐波畸变率。
5、给出了一种通用微分离散化函数,可以根据不同的系统特性设计微分离散化函数,有利于提高微分离散化稳定性及灵活设计幅相特性。
附图说明
图1是本发明的基于虚拟同步发电机的微网逆变器拓扑结构。
图2是本发明的虚拟同步发电机的功率外环控制框图。
图3是本发明的虚拟同步发电机的电压电流双环控制框图。
图4是本发明的基于虚拟同步发电机的微网逆变器等效数学模型。
具体实施方式
图1是本发明的实施例中基于虚拟同步发电机的微网逆变器拓扑结构。包括直流源Udc、直流侧滤波电容Cdc、三相半桥逆变电路、LC滤波器,直流侧滤波电容Cdc并联在所述直流源Udc的两端,直流源Udc的两个电源输出端分别与三相全桥逆变电路的两个输入端相连,三相全桥逆变电路的三相输出端与LC滤波器的三相输入端一一对应相连,LC滤波器的三相输出端分别与Dyn11型变压器的三角型侧相连接,变压器星型侧与三相电网Ea、Eb、Ec相连,电网相电压有效值为E,Lg为三相电网感抗对应的电感,LC滤波器由桥臂侧电感L和滤波电容C组成。
下面结合附图对本发明的优选方式作进一步详细的描述。
具体的,本实施例中的参数如下:直流母线电压Udc为550V,输出交流线电压有效值为380V/50Hz,额定容量为100kW,微网逆变器桥臂侧电感为L=0.5mH,微网逆变器滤波电容为C=200μF。变压器为100kVA270/400V Dyn11型变压器,微网逆变器采样频率fs为10kHz,因而Ts=100μs。
参见图1、2、3和4,本发明提供的一种基于虚拟同步发电机的微网逆变器均衡控制方法,主要步骤如下:
步骤1,采样及坐标变换;
所述采样包括采集以下数据:微网逆变器滤波电容电压uca,ucb,ucc,微网逆变器桥臂侧电感电流iLa,iLb,iLc,微网逆变器并网点电网电压ea,eb,ec。
所述坐标变换包括对以下数据进行坐标变换:对微网逆变器滤波电容电压uca,ucb,ucc和桥臂侧电感电流iLa,iLb,iLc,分别进行单同步旋转坐标变换得到滤波电容电压dq的分量Ucd,Ucq和桥臂侧电感电流dq的分量ILd,ILq。
步骤2,根据步骤1中得到的滤波电容电压的dq分量Ucd,Ucq,通过通用的微分离散化方程计算滤波电容电流的dq分量Icd,Icq;根据步骤1得到的桥臂侧电感电流的dq分量ILd,ILq和滤波电容电流的dq分量Icd,Icq,经过输出电流计算方程得到输出电流的dq分量Iod,Ioq;经过有功功率计算方程和无功功率计算方程得到平均有功功率P和平均无功功率Q;对微网逆变器并网点电网电压ea,eb,ec经过锁相环环节得到电网角频率ωg。
步骤2.1,计算滤波电容电流的dq分量Icd,Icq;
令滤波电容电压Ucd,Ucq的离散序列为Ucd(n),Ucq(n),滤波电容电流Icd,Icq的离散序列为Icd(n),Icq(n),则计算滤波电容电流的通用的微分离散化方程为:
其中,C为滤波电容,Ts为微网逆变器采样频率,K为离散序列点数,n,k为自然数,即n=0,1,2,3,4......,k=0,1,2,3,4......。
根据上述方程可以求得滤波电容电流Icd,Icq的离散序列为Icd(n),Icq(n),从而可得滤波电容电流Icd,Icq。
通用离散化方程的参数选择综合考虑差分方程稳定性条件,微分的频率响应以及DSP计算量。在本实施例中,取N=7,K=2,kn=4,kn-1=2,kn-2=1,。
步骤2.2,计算输出电流的dq分量Iod,Ioq;
根据步骤2.1得到的滤波电容电流的dq分量Icd,Icq,经过输出电流计算方程得到输出电流的dq分量Iod,Ioq,所述的输出电流计算方程为:
Iod=ILd-Icd
Ioq=ILq-Icq
步骤2.3,根据有功功率计算方程和无功功率计算方程计算平均有功功率P和平均无功功率Q;
有功功率计算方程为:
无功功率计算方程为:
其中,Qpq为功率计算方程品质因数、ωh为陷波器需要滤除的谐波角频率、s为拉普拉斯算子、τ为一阶低通滤波器的时间常数,h为待抑制的谐波次数。
在本实施例中,考虑主要滤除的谐波次数为2次和3次谐波,因此选取h=2,3,此时ωh=628.3186rad/s,942.4779rad/s。一阶低通滤波器主要考虑滤除高次谐波,且不影响动态响应,一般取τ≤2e-3s,本例取值τ=1.5e-4s;品质因数Qpq主要考虑陷波器的滤波效果,在本例中,选取Qpq=0.5。
步骤3,根据步骤2中得到的平均有功功率P、电网角频率ωg和微网逆变器给定的有功功率指令P0、微网逆变器给定有功功率指令P0时的额定角频率ω0,经过功角控制方程及电网频率逻辑判断环节得到虚拟同步发电机的角频率ω,对ω积分得到虚拟同步发电机的矢量角θ。
其中,功角控制方程为:
其中,ω0为微网逆变器给定有功功率指令P0时的额定角频率,m为功角控制下垂系数,J为模拟同步发电机机组的虚拟转动惯量,s为拉普拉斯算子,D1为微网逆变器频率反馈系数,D2为电网频率反馈系数。
采用具有回差特性的电网频率逻辑判断环节为:
当ωg≥ωg1时,令D1=-D2<0,即:
当ωg<ωg2时,令ω0=ωg,D1=-D2<0,即:
其中,ωg1为电网频率低阈值,ωg2为电网频率高返回值,ωg1<ωg2。
功角控制方程表明了微网逆变器有功功率下垂曲线关系、虚拟惯量大小和阻尼大小。其中,虚拟惯量标明了系统频率的变化率,为了保证系统频率变化平稳,需要有较大的虚拟惯量;然而虚拟惯量相当于在系统中加入了一阶惯性环节,太大的虚拟惯量有可能导致系统的不稳定。因而参数选择需要折中处理。为保证系统稳定性,在本实施例中,惯性时间常数范围在τvirtual=Jω0m≤2e-3s;功角控制方程中的有功功率下垂曲线关系包括三个系数,功角控制下垂系数m表示下垂曲线的斜率,取值原则为100%的有功功率变化时,频率变化0.5Hz以内;给定有功功率指令P0和相对应的额定角频率ω0表示下垂曲线的位置关系,主要考虑微网逆变器输出有功功率为P0时,其输出频率大小。
在本实施例中,功角控制下垂系数取值为根据惯性时间常数取值原则取τvirtual=Jω0m=1.5e-3s,可得J=0.2kg·m2,为保证控制运行时能量不流向直流侧,给定有功功率指令取值为P0=1kW,此时对应的额定角频率取值为ω0=314.1593rad/s。
具有回差特性的电网频率逻辑判断环节表明了功率外环调节器频率补偿特性与限功率运行特性的切换逻辑关系。当电网频率大于一定值时,微网逆变器调节功率输出来稳定电网频率波动,当电网频率低于一定值超出微网逆变器功率调节范围时,微网逆变器应该限功率运行。根据运行的电网频率波动范围以及微网逆变器功率调节范围以及下垂系数选取原则来选取电网频率逻辑判断的逻辑值。在本实施例中,选择ωg1=49.2Hz,ωg2=49.5Hz。D1,D2表明了外环功率环的阻尼特性,根据上述方程基于虚拟同步发电机的微网逆变器数学模型如图4所示,进而可得有功功率传递函数为:
其中,为功角传递函数,E为电网相电压有效值,X为微网逆变器每相等效输出阻抗。在本实施例中,微网逆变器的等效输出阻抗为额定阻抗的5%,因而Ks等效为Ks≈20×100kW。
根据控制系统二阶振荡方程可得系统的阻尼为其中ζ>0,将m,J,ω0,Ks带入可得D1的取值范围为D1<40,在本实施例中,取ζ=0.7,则D1=-15640,D2=15640。
步骤4,根据步骤2中得到的平均无功功率Q和微网逆变器给定的无功功率指令Q0、电压指令U0,经过无功控制方程得到虚拟同步发电机的端电压U*。
无功控制方程为:
U*=U0+nq(Q0-Q)
其中,U0为微网逆变器给定无功功率指令Q0时的额定输出电容电压、nq为无功-电压下垂系数。
无功-电压下垂系数nq取值原则为100%的无功功率变化时,电压幅值变化在2%之内;给定无功功率指令Q0和相对应的额定输出电容电压U0表示下垂曲线的位置关系,主要考虑微网逆变器输出无功功率为Q0时,其输出电压大小。
在本实施例中,无功-电压下垂系数取值为给定无功功率指令Q0考虑系统输出无功功率为Q0=0,此时对应的额定输出电容电压U0=380V。
步骤5,先根据步骤4中得到的端电压U*和步骤1中得到的滤波电容电压dq的分量Ucd,Ucq,通过电压控制方程得到电流指令信号再根据电流指令信号步骤1中的桥臂侧电感电流的dq分量ILd,ILq和步骤2得到的滤波电容电流的dq分量Icd,Icq,通过加权电流控制方程得到控制信号Ud,Uq。其中,电压控制方程为:
其中,Kp为电压环比例控制系数、Ki为电压环积分控制系数、Kr为电压环谐振控制器比例系数,Qu为电压环准谐振调节器品质因数,ωh为陷波器需要滤除的谐波角频率,s为拉普拉斯算子,h为待抑制的谐波次数。
电压控制方程中的参数主要考虑控制系统的稳定性和动稳态性能;在本实施例中,取Kp=0.03,Ki=0.8,准谐振调节器主要考虑消除系统中的奇次谐波,取h=3,5,7,9,11,因而角频率分别等于ωh=942.5rad/s,1570.8rad/s,2199.1rad/s,2827.4rad/s,3455.8rad/s。
品质因数Qu主要考虑谐振调节器的增益和稳定性,在本例中,选取Qu=0.7;准谐振控制器比例系数综合考虑电压环的动稳态控制性能和系统稳定性,在本例中,选取Kr=100。
电流控制方程为:
其中,Kpi为电流环比例控制系数,Kri电流环谐振控制器比例系数,w1为电感电流的权重系数,w2为电容电流的权重系数,Kf为电压前馈系数,Qi为电流环准谐振调节器品质因数,s为拉普拉斯算子。
电流控制方程中的参数主要考虑控制系统的阻尼特性和直流分量抑制能力;在本实施例中,取Kpi=0.05,准谐振调节器主要考虑消除系统中的直流分量,取h=1,因而角频率等于ωh=314.15rad/s。
品质因数Qi主要考虑谐振调节器的增益和稳定性,在本例中,选取Qi=0.7;准谐振控制器比例系数综合考虑电流环的直流分量抑制能力和系统稳定性,在本例中,选取Kri=50。
电感电流与电容电流加权反馈控制环节主要考虑微网逆变器孤岛运行输出电压的动态响应与并联均流之间的均衡。在本实施例中,取w1=0.3,w2=0.7。
步骤6,将步骤5中得到的控制信号Ud,Uq经过单同步旋转坐标反变换得到三相桥臂电压控制信号Ua,Ub,Uc,再根据Ua,Ub,Uc生成开关管的PWM控制信号。