CN111130366B - 三相vienna整流器的无网压传感器的预测控制方法 - Google Patents

三相vienna整流器的无网压传感器的预测控制方法 Download PDF

Info

Publication number
CN111130366B
CN111130366B CN202010023572.7A CN202010023572A CN111130366B CN 111130366 B CN111130366 B CN 111130366B CN 202010023572 A CN202010023572 A CN 202010023572A CN 111130366 B CN111130366 B CN 111130366B
Authority
CN
China
Prior art keywords
voltage
phase
grid
formula
vienna rectifier
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010023572.7A
Other languages
English (en)
Other versions
CN111130366A (zh
Inventor
李山
苏新柱
郭强
肖蕙蕙
陈岚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University of Technology
Original Assignee
Chongqing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University of Technology filed Critical Chongqing University of Technology
Priority to CN202010023572.7A priority Critical patent/CN111130366B/zh
Publication of CN111130366A publication Critical patent/CN111130366A/zh
Application granted granted Critical
Publication of CN111130366B publication Critical patent/CN111130366B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M7/00Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
    • H02M7/02Conversion of ac power input into dc power output without possibility of reversal
    • H02M7/04Conversion of ac power input into dc power output without possibility of reversal by static converters
    • H02M7/12Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
    • H02M7/21Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
    • H02M7/217Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only
    • H02M7/219Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only in a bridge configuration
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/16Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30Reactive power compensation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Rectifiers (AREA)

Abstract

本发明公开了一种三相VIENNA整流器的无网压传感器的预测控制方法,包括如下步骤:S1、对三相VIENNA整流器在三相静止坐标系中的数学模型进行Clarke变换,建立两相静止坐标系下三相VIENNA整流器的数学模型;S2、采用虚拟磁链算法估算电网电压;S3、采用二阶低通滤波器取代步骤S2中的积分项;S4、以功率跟踪误差最小确定k+1时刻iα和iβ的预测值
Figure DDA0002361661430000011
Figure DDA0002361661430000012
求得交流侧电压矢量Vα、Vβ,得到空间矢量调制信号。本发明具有无需网压传感器和锁相环,即可对实际电网电压准确快速地追踪,保持电网电压电流同相运行等优点。

Description

三相VIENNA整流器的无网压传感器的预测控制方法
技术领域
本发明属于电力电子变换技术领域,特别的涉及一种三相VIENNA整流器的无网压传感器的预测控制方法。
背景技术
VIENNA整流器是一种三电平拓扑,与两电平整流器相比,该型整流器谐波含量更低、效率更高。与传统三电平整流器相比,该型整流器所需要的开关器件更少,开关管所承受的应力低,也不需要设置死区时间。具有功率开关承受电压应力低、功率开关数量少、功率密度大、输入电流谐波含量少、对各种类型的PFC控制方法均具有良好的适应性的优点,因此在功率因数校正技术研究中得到了广泛的应用。
目前,三相VIENNA整流器控制策略主要分为电流控制和功率控制两大类。其中电流控制又可分为直接电流控制和间接电流控制,前者以网侧电流作为反馈和被控量形成电流闭环,后者无电流环,但两者都需采集电网电压形成电压环。而功率控制需提取电网电压基波幅值、相位等信息来控制网侧有功、无功功率,从而实现直流侧输出电压恒定、单位功率因数运行等控制目标。由此可见,现有的各控制策略都依赖于电网电压信息的正确提取。而在实际应用中,一般通过安装电网电压传感器来获得电网电压信号,而安装电网电压传感器既会增加系统成本和硬件复杂度、降低系统可靠性,还会引起噪声、分辨率受限等问题。
发明内容
针对上述现有技术的不足,本发明所要解决的技术问题是:如何提供一种无需网压传感器和锁相环,即可对实际电网电压准确快速地追踪,保持电网电压电流同相运行的三相VIENNA整流器的无网压传感器的预测控制方法。
为了解决上述技术问题,本发明采用了如下的技术方案:
一种三相VIENNA整流器的无网压传感器的预测控制方法,其特征在于,包括如下步骤:
S1、对三相VIENNA整流器在三相静止坐标系中的数学模型进行Clarke变换,建立两相静止坐标系下三相VIENNA整流器的数学模型:
Figure BDA0002361661410000011
其中,三相VIENNA整流器在三相静止坐标系中的数学模型为:
Figure BDA0002361661410000012
式中:Ea、Eb、Ec分别为三相静止坐标系下A、B、C相的电网电压;ia、ib、ic分别为三相静止坐标系下A、B、C相的电感电流;L为滤波电感;R为滤波电感和线路的等效电阻;Vao、Vbo、Vco分别表示A、B、C相交流侧与直流侧电容中性点o之间的电压;VoN表示o与电源中性点N之间的电压;Eα、Eβ、iα和iβ分别为两相静止坐标系下的电网电压和交流侧电感电压;Vα和Vβ为两相静止坐标系下交流侧电压分量;
S2、采用虚拟磁链算法估算电网电压:
Figure BDA0002361661410000021
S3、采用二阶低通滤波器取代步骤S2中的积分项,其传递函数如下:
Figure BDA0002361661410000022
式中:ωc为截止频率;
S4、在一个采样周期内对电流微分量进行离散化:
Figure BDA0002361661410000023
将上式代入两相静止坐标系下三相VIENNA整流器的数学模型,得到:
Figure BDA0002361661410000024
式中:Ts为三相VIENNA整流器的采样周期,且远小于电网基波周期;RTs乘积项较小,可忽略;
Figure BDA0002361661410000025
Figure BDA0002361661410000026
分别为k+1时刻iα和iβ的预测值,以功率跟踪误差最小确定
Figure BDA0002361661410000027
Figure BDA0002361661410000028
求得交流侧电压矢量Vα、Vβ,得到空间矢量调制信号;其中,
Figure BDA0002361661410000029
Figure BDA00023616614100000210
为:
Figure BDA00023616614100000211
式中:Pref和Qref为瞬时功率参考值;P和Q为瞬时功率;||E||=ω||ψ||,||ψ||为对ψ求模,ψ=∫Edt为电网虚拟磁链。
进一步的,还包括如下步骤:
S5、采用二步预测法补偿信号延迟,得到在两相静止坐标系下的期望交流侧电压矢量预测值
Figure BDA00023616614100000212
Figure BDA00023616614100000213
Figure BDA00023616614100000214
式中:
Figure BDA00023616614100000215
Figure BDA00023616614100000216
分别为k+2时刻iα和iβ的预测值:
Figure BDA0002361661410000031
式中:Q*、P*、E*
Figure BDA0002361661410000032
Figure BDA0002361661410000033
分别为Q、P、E、Eα和Eβ在k+1时刻的预测值。
进一步的,还包括如下步骤:
S6、以下一周期中点位置的电网电压预测值E(k+1/2)作为电网电压预测值E*,确定期望交流侧电压矢量Vα、Vβ,具体步骤为:
根据三阶拉格朗日插值法推导下一周期中点位置电压:
Figure BDA0002361661410000034
式中:x0=k-2、x1=k-1、x2=k、x=k+1/2;
代入后化简得到:
Figure BDA0002361661410000035
具体的,所述步骤S4中,k+1时刻iα和iβ的预测值
Figure BDA0002361661410000036
Figure BDA0002361661410000037
采用如下步骤得到:
先根据瞬时功率理论得到:
Figure BDA0002361661410000038
对上式进行求导并离散化,得到采样周期内整流器的瞬时功率预测值满足:
Figure BDA0002361661410000039
以满足功率跟踪误差最小为目标去预测电流,可定义目标函数:
F=ΔP2+ΔQ2
其中ΔP、ΔQ为瞬时功率参考值Pref、Qref与预测值P*、Q*之差;
根据极值定理,跟踪误差要达到最小,需满足
Figure BDA00023616614100000310
Figure BDA00023616614100000311
则有:
Figure BDA00023616614100000312
具体的,所述步骤S5中,k+2时刻iα和iβ的预测值
Figure BDA00023616614100000313
Figure BDA00023616614100000314
采用如下步骤得到:
对瞬时功率表达式进行二阶求导,并离散化得到:
Figure BDA0002361661410000041
确定二步预测优化后的目标函数:
F'=(Pref-P**)2+(Qref-Q**)2
式中:Q**和P**分别为Q和P在k+2时刻的预测值;
根据极值定理,跟踪误差要达到最小,需满足
Figure BDA0002361661410000042
则有:
Figure BDA0002361661410000043
具体的,所述步骤S1中,还包括如下:
S11、对三相VIENNA整流器进行建模前,假定:
1)电网平衡;
2)所有开关管和二极管均为理想器件;
3)开关频率远大于电网基波频率;
S12、定义A、B、C相的开关函数si,i为a、b和c;将开关函数Si分解成Sio、Sip、Sin三个开关,当开关导通时其值为1,关断为0,且满足:
Figure BDA0002361661410000044
即当Si导通时,Sip=0、Sin=0;当Si关断且ii>0,Sip=1、Sin=0;当Si关断且ii<0,Sip=0、Sin=1;其中,Si表示i相交流侧电压的电位状态,ii表示i相电流;
交流侧电压与直流侧电容中性点电压关系式:
Figure BDA0002361661410000045
式中:Vc1和Vc2分别为直流侧上下滤波电容上的电压;
由于电网平衡,即三相交流电压对称,则:
Figure BDA0002361661410000046
将上式与三相VIENNA整流器在三相静止坐标系中的数学模型合并后得到:
Figure BDA0002361661410000051
将上式与交流侧电压与直流侧电容中性点电压关系式相加得到:
Figure BDA0002361661410000052
具体的,所述步骤S2中,采用虚拟磁链算法估算电网电压前,采用如下步骤:
根据电网电压矢量超前虚拟电网磁链矢量π/2,且幅值增长ω倍的特性,可得:
Figure BDA0002361661410000053
其中,||ψ||为对ψ求模,ψ=∫Edt,进而推导出电网虚拟磁链与交流侧电压关系式:
Figure BDA0002361661410000054
对两相静止坐标系下三相VIENNA整流器的数学模型的两边同时积分,去掉微分项得:
Figure BDA0002361661410000055
将上两式结合可推算出基于虚拟磁链法估算电网电压的公式:
Figure BDA0002361661410000056
式中,ω为电网电压基波角频率;
综上所述,本发明具有如下优点:
1、无需网压传感器、锁相环,即可对实际电网电压准确快速地追踪,保持电网电压电流同相运行。
2、适用于三相VIENNA整流器其它拓扑结构,具有良好的输入、输出及动态特性。
附图说明
图1为三相VIENNA整流器的拓扑图。
图2为aβ坐标系下整流器稳态矢量图。
图3为电网电压观测器的结构框图。
图4为系统结构框图。
图5为系统稳态下电网电压实际值与估算值波形。
图6为额定负载运行、负载突增及负载突减三个时段的A相电网电压与输入电流波形。
图7为额定负载运行时A相电流谐波畸变率。
图8为直流侧电压波形。
图9为直流侧电容电压波形。
图10为电网瞬时有功功率、无功功率和功率因数波形。
具体实施方式
下面结合实施例对本发明作进一步的详细说明。
一、三值逻辑开关函数等效法推导VIENNA整流器在两相静止坐标系下的数学模型
如图1所示,为三相VIENNA整流器的拓扑图,图中Ea、Eb、Ec分别为A、B、C相电网电压,L为滤波电感,R为滤波电感和线路的等效电阻,Va、Vb、Vc分别为A、B、C相交流侧电压,C1、C2分别为直流侧上下滤波电容,Rload为阻性负载。
三相VIENNA整流器在三相静止坐标系中的数学模型为:
Figure BDA0002361661410000061
其中Vao、Vbo、Vco分别表示a、b、c与o之间的电压,VoN表示o与N之间的电压。为了便于推导Va、Vb、Vc表达式,假定所有开关管和二极管均为理想元件,电网平衡,开关频率远大于电网基波频率。此时,三相桥臂开关管可等效为一个三值逻辑开关函数:
定义A、B、C相的开关函数si,i为a、b和c;将开关函数Si分解成Sio、Sip、Sin三个开关,当开关导通时其值为1,关断为0,且满足:
Figure BDA0002361661410000062
其中,Si表示i相交流侧电压的电位状态,ii表示i相电流;即当Si导通时,Sip=0、Sin=0;当Si关断且ii>0,Sip=1、Sin=0;当Si关断且ii<0,Sip=0、Sin=1;各开关量之间的关系如下表1所示。
表1各开关量之间关系
Table 1 The relation between the switch quantities
S<sub>i</sub> i<sub>i</sub> s<sub>i</sub> s<sub>in</sub> s<sub>in</sub> s<sub>in</sub>
导通 电位o 1 0 0
关断 大于0 电位p 0 1 0
关断 小于0 电位n 0 0 1
交流侧电压与直流侧电容中性点电压关系式:
Figure BDA0002361661410000071
三相交流电压对称时,有:
Figure BDA0002361661410000072
合并式(1)和(4)得:
Figure BDA0002361661410000073
其中,Sip、Sin(i=a,b,c)可根据表1中Si、ii之间的关系得到。以a相为例,当Sa导通时,Sap=0、San=0;当Sa关断且ia>0,Sap=1、San=0;当Sa关断且ia<0,Sap=0、San=1。
Figure BDA0002361661410000074
将式(3)和(5)相加得:
Figure BDA0002361661410000075
式(1)进行Clarke变换,可得整流器在两相静止坐标系中的数学模型:
Figure BDA0002361661410000076
如果不考虑式(8)中的微分项,单从算式本身看来,应该是可以用来进行电网电压估算的。但在实际应用中,如果根据式(8)进行电网电压估算得到无网压传感器控制策略,其中的微分计算,会导致较高的噪声灵敏度。故下文提出虚拟磁链模型去掉微分项。
二、虚拟磁链算法估算电网电压
对式(8)两边同时积分,去掉微分项得:
Figure BDA0002361661410000081
根据磁链公式有电网虚拟磁链ψ=∫Edt。电网电压E、电网虚拟磁链ψ及交流侧电压V在αβ坐标系中的关系如图2所示。
图中ω为电网电压基波角频率。根据电网电压矢量超前虚拟电网磁链矢量π/2,且幅值增长ω倍的特性,并结合图2中各变量矢量位置关系,可得:
Figure BDA0002361661410000082
其中,||x||表示对x求模。由式(10)推导出电网虚拟磁链与交流侧电压关系式:
Figure BDA0002361661410000083
结合式(9)、式(11)可推算出基于虚拟磁链法估算电网电压的公式:
Figure BDA0002361661410000084
式中交流侧电压分量Vα、Vβ源于三相电压Va、Vb、Vc的Clarke变换。可见,基于虚拟磁链法的电网电压估算式中存在积分项,而积分环节的初值和积分偏置问题会导致整流器冲击电流过大甚至无法启动。故下面提出用二阶低通滤波器取代积分项来改进此电压估算策略。
三、基于二阶低通滤波器取代积分项改进的电压观测器
若对一个参数施以积分运算,其幅值会衰减为原来的ω倍、相角滞后π/2,因此,式(12)中积分项整体ω∫Vdt的输出幅值仍为V,只是相角会滞后π/2。
由一阶低通滤波器特性可知:输入信号频率ω=ωc(截止频率)时,输出信号的幅度为
Figure BDA0002361661410000085
相位偏移-π/4。故可利用两个增益K为
Figure BDA0002361661410000086
的一阶低通滤波器取代纯积分项,其传递函数如下:
Figure BDA0002361661410000091
根据式(12)和(13)设计一种基于二阶低通滤波器改进的电网电压观测器,从而避免纯积分运算,其结构如图3所示。
四、预测电流控制的基本原理
由于三相VIENNA整流器采样周期Ts远小于电网基波周期,故在一个采样周期内对电流微分量进行离散化得:
Figure BDA0002361661410000092
将式(14)代入式(8)得:
Figure BDA0002361661410000093
其中,RTs乘积项较小,故可忽略;i*代表网侧电流k+1时刻预测值。说明书中标有*号的变量均表示k+1时刻预测值,标有**号表示k+2时刻预测值。
根据瞬时功率理论有:
Figure BDA0002361661410000094
对式(16)求导并离散化,得到采样周期内整流器的瞬时功率预测值:
Figure BDA0002361661410000095
以满足功率跟踪误差最小为目标去预测电流,可定义目标函数F=ΔP2+ΔQ2,其中ΔP、ΔQ为瞬时功率参考值(Pref、Qref)与预测值(P*、Q*)之差。
根据极值定理,跟踪误差要达到最小,需满足
Figure BDA0002361661410000096
成立,则有:
Figure BDA0002361661410000097
由式(18)得到预测电流值i*,代入式(15)求得交流侧电压矢量Vα、Vβ后,即可得到空间矢量调制信号。
五、二步预测优化
在实际运用中,由于系统的信号采集和处理均无法瞬时完成,这将导致期望开关信号滞后一个周期。为此,本文利用二步预测法补偿信号延迟,通过算得V**代替V*导入空间矢量调制,以补偿实际应用中的滞后。
二步预测优化后的目标函数:
F'=(Pref-P**)2+(Qref-Q**)2。 (19)
对式(16)二阶求导,并离散化得:
Figure BDA0002361661410000101
根据
Figure BDA0002361661410000102
则有:
Figure BDA0002361661410000103
可求得在两相静止坐标系下的期望交流侧电压矢量预测值V*
Figure BDA0002361661410000104
六、电网电压预测值的改进
式(22)中电网电压预测值E*并不能直接得到,传统预测控制直接令E*≈E处理,但在电压较高而开关频率不高的情况下,这会产生不可忽略的误差。通常用下一周期中点位置值比直接用当前值E代替更为精准。故本文根据三阶拉格朗日插值法推导下一周期中点位置电压。
拉格朗日插值公式:
Figure BDA0002361661410000105
式中x0=k-2、x1=k-1、x2=k、x=k+1/2。
代入化简后可得:
Figure BDA0002361661410000106
再用E(k+1)≈E(k+1/2)代入式(22)可更为精准地得到期望交流侧电压矢量。
具体实施时,如图4所示,图中虚线框内左侧部分的工作流程如下:
1、将整流器三相输入电流ia、ib、ic经Clarke坐标变换,得到两相静止坐标系下αβ轴分量iα、iβ
2、根据三相交流侧电压的电位状态Sabc,结合整流器三相输入电流ia、ib、ic和直流侧电容电压Vc1、Vc2通过式(7)求得交流侧三相电压Va、Vb、Vc,再经Clarke变换即式(6)得到交流测电压在两相静止坐标系下的αβ轴分量Vα、Vβ
3、将得到的交流测电压αβ轴分量Vα、Vβ与交流侧电流αβ轴分量iα、iβ结合,通过设计两个增益为
Figure BDA0002361661410000107
的一阶低通滤波器取代虚拟磁链法推导的电网电压估算式(12)中的纯积分项,即一种基于二阶低通滤波器改进的电网电压观测器(图3),得到电网电压Eα、Eβ
4、将得到的电网电压Eα、Eβ结合交流侧电流αβ轴分量iα、iβ,根据瞬时功率理论,经瞬时功率计算式(16),得到整流器瞬时功率P、Q。
图中虚线框内右侧部分的工作流程为:
5、将直流侧负载电压Vdc与其电压参考值进行比较输出的偏差信号送入PI控制器,控制器生成控制信号,将得到的控制信号与负载电压Vdc相乘,得到有功功率参考值Pref
公式表述:Pref=[(Vdcref-Vdc)(Kp+Ki/s)]Vdc
6、将得到的有功功率参考值Pref和无功功率参考值Qref(=0)结合得到的整流器瞬时功率P、Q和电网电压Eα、Eβ,经过电网电流预测计算(19)得到网侧电流k+1拍预测值i*在两相静止坐标系下αβ轴分量iα *、iβ *
7、将所述αβ轴分量经瞬时功率预测计算式(17),得到采样周期内整流器瞬时功率k+1拍预测值P*、Q*
8、根据得到的瞬时功率预测值结合电网电压预测值,经二步预测优化策略计算网测电流k+2拍预测值,即通过式(22),得到网测电流k+2拍预测值i**的αβ轴分量iα **、iβ **
9、将得到的网测电流k+2拍预测值结合电网电压预测值,经交流侧电压矢量式(23),得到在两相静止坐标系下交流侧电压矢量预测值Vα *、Vβ *
10、根据得到的交流侧电压矢量预测值经过空间矢量脉宽调制,得到三相交流侧电压的电位状态Sabc
七、实验结果及分析
为验证本文所提无网压传感器改进预测控制方法的正确性,在Matlab/Simulink仿真平台中进行仿真。主电路参数如表2所示。
表2 VIENNA整流器参数指标
Table 2 Parameters of VIENNA rectifier
参数 数值 参数 数值
电网电压/V 220 网侧电感/<u>mH</u> 12
直流电压/V 700 直流侧电容/<u>mF</u> 4.4
开关频率/kHz 20 额定负载/Ω 24.5
仿真时长t=1.5s。系统运行至t=1s时,负载由原值Rload=24.5Ω突变为49Ω;运行至t=1.5s时,由49Ω突变为36.75Ω。
图5为系统稳态下电网电压实际值与估算值波形,可以看出实际值与估算值波形几乎保持一致,验证了本文提出的三相VIENNA整流器电网电压估算算法的正确性。图6为额定负载运行、负载突增及负载突减三个时段的A相电网电压与输入电流波形,可以看出整个过程中电网电压与相电流保持同相,功率因数几乎为1。
图7为在额定负载下,A相电流总谐波畸变率,其值为4.68%,说明系统具有良好的输入特性。由图8可知,系统启动及负载突变过程中,直流电压均能较快调节至给定值,稳定后静态误差几乎为0,功率突减一倍时超调量为13.14%,说明系统具有良好的动态特性。结合图8、图9可知,中点电位抑制效果明显,直流电压纹波仅为0.5V,具有良好的输出特性。
图10为整个过程,整流器输入有功功率、无功功率及功率因数仿真波形。该图说明VIENNA整流器发生负载突增、突减时,瞬时有功功率能快速准确地跟踪给定值,而无功功率稳定在零值附近,整个运行过程中功率因数始终大于0.95,实现了系统的功率因数校正。
由此可见,本申请基于三相VIENNA整流器提出了一种无网压传感器改进预测控制策略。该策略基于虚拟磁链理论推导出三相VIENNA整流器电网电压估算算法,并利用二阶低通滤波器改进算法中纯积分环节,以解决积分初值和偏置问题;控制部分采用基于瞬时功率理论的预测电流控制方法,并针对实际系统中采样和处理会引起信号滞后问题,提出二步预测电流法加以改进,以得到更为精准的交流侧电压矢量。最后对该策略进行了仿真验证,仿真结果表明,本申请提出的基于三相VIENNA整流器的无网压传感器改进预测控制策略具有以下特点:
1)无需网压传感器、锁相环,即可对实际电网电压准确快速地追踪,保持电网电压电流同相运行。
2)适用于三相VIENNA整流器其它拓扑结构,具有良好的输入、输出、动态及稳态特性。
以上所述仅为本发明的较佳实施例而已,并不以本发明为限制,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种三相VIENNA整流器的无网压传感器的预测控制方法,其特征在于,包括如下步骤:
S1、对三相VIENNA整流器在三相静止坐标系中的数学模型进行Clarke变换,建立两相静止坐标系下三相VIENNA整流器的数学模型:
Figure FDA0002665690380000011
其中,三相VIENNA整流器在三相静止坐标系中的数学模型为:
Figure FDA0002665690380000012
式中:Ea、Eb、Ec分别为三相静止坐标系下A、B、C相的电网电压;ia、ib、ic分别为三相静止坐标系下A、B、C相的电感电流;L为滤波电感;R为滤波电感和线路的等效电阻;Vao、Vbo、Vco分别表示A、B、C相交流侧与直流侧电容中性点o之间的电压;VoN表示o与电源中性点N之间的电压;Eα、Eβ、iα和iβ分别为两相静止坐标系下的电网电压和交流侧电感电压;Vα和Vβ为两相静止坐标系下交流侧电压分量;
S2、采用虚拟磁链算法估算电网电压:
Figure FDA0002665690380000013
S3、采用二阶低通滤波器取代步骤S2中的积分项,其传递函数如下:
Figure FDA0002665690380000014
式中:ωc为截止频率;
S4、在一个采样周期内对电流微分量进行离散化:
Figure FDA0002665690380000015
将上式代入两相静止坐标系下三相VIENNA整流器的数学模型,得到:
Figure FDA0002665690380000016
式中:Ts为三相VIENNA整流器的采样周期,且远小于电网基波周期;RTs乘积项较小,可忽略;
Figure FDA0002665690380000017
Figure FDA0002665690380000018
分别为k+1时刻iα和iβ的预测值,以功率跟踪误差最小确定
Figure FDA0002665690380000019
Figure FDA00026656903800000110
求得交流侧电压矢量Vα、Vβ,得到空间矢量调制信号;其中,
Figure FDA00026656903800000111
Figure FDA00026656903800000112
为:
Figure FDA0002665690380000021
式中:Pref和Qref为瞬时功率参考值;P和Q为瞬时功率;||E||=ω||ψ||,||ψ||为对ψ求模,ψ=∫Edt为电网虚拟磁链;
S5、采用二步预测法补偿信号延迟,得到在两相静止坐标系下的期望交流侧电压矢量预测值
Figure FDA0002665690380000022
Figure FDA0002665690380000023
Figure FDA0002665690380000024
式中:
Figure FDA0002665690380000025
Figure FDA0002665690380000026
分别为k+2时刻iα和iβ的预测值:
Figure FDA0002665690380000027
式中:Q*、P*、E*
Figure FDA00026656903800000214
Figure FDA00026656903800000215
分别为Q、P、E、Eα和Eβ在k+1时刻的预测值;
S6、以下一周期中点位置的电网电压预测值E(k+1/2)作为电网电压预测值E*,确定期望交流侧电压矢量Vα、Vβ,具体步骤为:
根据三阶拉格朗日插值法推导下一周期中点位置电压:
Figure FDA0002665690380000028
式中:x0=k-2、x1=k-1、x2=k、x=k+1/2;
代入后化简得到:
Figure FDA0002665690380000029
2.如权利要求1所述的三相VIENNA整流器的无网压传感器的预测控制方法,其特征在于,所述步骤S4中,k+1时刻iα和iβ的预测值
Figure FDA00026656903800000210
Figure FDA00026656903800000211
采用如下步骤得到:
先根据瞬时功率理论得到:
Figure FDA00026656903800000212
对上式进行求导并离散化,得到采样周期内整流器的瞬时功率预测值满足:
Figure FDA00026656903800000213
以满足功率跟踪误差最小为目标去预测电流,可定义目标函数:
F=ΔP2+ΔQ2
其中,ΔP、ΔQ为瞬时功率参考值Pref、Qref与预测值P*、Q*之差;
根据极值定理,跟踪误差要达到最小,需满足
Figure FDA0002665690380000031
Figure FDA0002665690380000032
则有:
Figure FDA0002665690380000033
3.如权利要求2所述的三相VIENNA整流器的无网压传感器的预测控制方法,其特征在于,所述步骤S5中,k+2时刻iα和iβ的预测值
Figure FDA0002665690380000034
Figure FDA0002665690380000035
采用如下步骤得到:
对瞬时功率表达式进行二阶求导,并离散化得到:
Figure FDA0002665690380000036
确定二步预测优化后的目标函数:
F'=(Pref-P**)2+(Qref-Q**)2
式中:Q**和P**分别为Q和P在k+2时刻的预测值;
根据极值定理,跟踪误差要达到最小,需满足
Figure FDA0002665690380000037
则有:
Figure FDA0002665690380000038
4.如权利要求1所述的三相VIENNA整流器的无网压传感器的预测控制方法,其特征在于,所述步骤S1中,还包括如下:
S11、对三相VIENNA整流器进行建模前,假定:
1)电网平衡;
2)所有开关管和二极管均为理想器件;
3)开关频率远大于电网基波频率;
S12、定义A、B、C相的开关函数si,i为a、b和c;将开关函数Si分解成Sio、Sip、Sin三个开关,当开关导通时其值为1,关断为0,且满足:
Figure FDA0002665690380000041
即当Si导通时,Sip=0、Sin=0;当Si关断且ii>0时,Sip=1、Sin=0;当Si关断且ii<0时,Sip=0、Sin=1;其中,Si表示i相交流侧电压的电位状态,ii表示i相电流;
交流侧电压与直流侧电容中性点电压关系式:
Figure FDA0002665690380000042
式中:Vc1和Vc2分别为直流侧上下滤波电容上的电压;
由于电网平衡,即三相交流电压对称,则:
Figure FDA0002665690380000043
将上式与三相VIENNA整流器在三相静止坐标系中的数学模型合并后得到:
Figure FDA0002665690380000044
将上式与交流侧电压与直流侧电容中性点电压关系式相加得到:
Figure FDA0002665690380000045
5.如权利要求1所述的三相VIENNA整流器的无网压传感器的预测控制方法,其特征在于,所述步骤S2中,采用虚拟磁链算法估算电网电压前,采用如下步骤:
根据电网电压矢量超前虚拟电网磁链矢量π/2,且幅值增长ω倍的特性,可得:
Figure FDA0002665690380000046
其中,||ψ||为对ψ求模,ψ=∫Edt,ω为电网电压基波角频率,进而推导出电网虚拟磁链与交流侧电压关系式:
Figure FDA0002665690380000051
对两相静止坐标系下三相VIENNA整流器的数学模型的两边同时积分,去掉微分项得:
Figure FDA0002665690380000052
将上两式结合可推算出基于虚拟磁链法估算电网电压的公式:
Figure FDA0002665690380000053
CN202010023572.7A 2020-01-09 2020-01-09 三相vienna整流器的无网压传感器的预测控制方法 Active CN111130366B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010023572.7A CN111130366B (zh) 2020-01-09 2020-01-09 三相vienna整流器的无网压传感器的预测控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010023572.7A CN111130366B (zh) 2020-01-09 2020-01-09 三相vienna整流器的无网压传感器的预测控制方法

Publications (2)

Publication Number Publication Date
CN111130366A CN111130366A (zh) 2020-05-08
CN111130366B true CN111130366B (zh) 2020-10-27

Family

ID=70487535

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010023572.7A Active CN111130366B (zh) 2020-01-09 2020-01-09 三相vienna整流器的无网压传感器的预测控制方法

Country Status (1)

Country Link
CN (1) CN111130366B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111555642A (zh) * 2020-05-18 2020-08-18 西安交通大学 一种预测直接功率控制方法
CN111628660B (zh) * 2020-06-05 2021-08-27 浙江鲲悟科技有限公司 压缩机变频调速系统及直流母线电压的控制方法
CN111628659A (zh) * 2020-06-05 2020-09-04 浙江鲲悟科技有限公司 压缩机变频调速系统及直流母线电压的控制方法
CN112701939B (zh) * 2020-12-07 2022-05-24 华南理工大学 Vienna整流器电流预测控制方法
CN112448604A (zh) * 2020-12-14 2021-03-05 上海空间电源研究所 一种整流器控制电路
CN114362554A (zh) * 2021-12-16 2022-04-15 西安理工大学 一种Vienna整流器的定频控制方法、系统、计算机设备及存储介质
CN116633126B (zh) * 2023-07-24 2023-10-17 成都希望森兰智能制造有限公司 一种无网压传感器的变频器老化系统功率因数控制方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102780222A (zh) * 2011-05-12 2012-11-14 Abb公司 用于控制在孤岛模式下的电网的方法和设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101771361A (zh) * 2010-02-11 2010-07-07 哈尔滨工业大学 无交流电压传感器并网逆变器的直接功率控制方法
US10651654B2 (en) * 2017-11-07 2020-05-12 State Grid Corporation Of China Model predictive controller for autonomous hybrid microgrids
GB2574645B (en) * 2018-06-14 2020-07-15 Zhong Qingchang Passive virtual synchronous machine with bounded frequency and virtual flux
CN109787491A (zh) * 2019-01-09 2019-05-21 三峡大学 基于虚拟磁链的三相Vienna整流器预测直接功率控制方法
CN110034690B (zh) * 2019-03-29 2020-10-23 南通大学 一种Vienna整流器模型预测虚拟磁链控制方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102780222A (zh) * 2011-05-12 2012-11-14 Abb公司 用于控制在孤岛模式下的电网的方法和设备

Also Published As

Publication number Publication date
CN111130366A (zh) 2020-05-08

Similar Documents

Publication Publication Date Title
CN111130366B (zh) 三相vienna整流器的无网压传感器的预测控制方法
CN107317490B (zh) 一种基于三相Vienna整流器的无差拍预测直接功率控制方法
CN109004669B (zh) 基于干扰观测器补偿的三相并网逆变器改进无源控制方法
CN110034690B (zh) 一种Vienna整流器模型预测虚拟磁链控制方法
CN107611971B (zh) 针对网压谐波畸变的网侧逆变器谐振全阶滑模控制方法
CN110165924B (zh) 一种单相pwm整流器的改进无差拍控制方法
CN110460250B (zh) 一种三相pwm整流器直接功率控制方法
CN110492494B (zh) 基于内模理论的电压源变换器直接功率控制方法
CN110011359B (zh) 一种有限集模型预测控制下的并网逆变器参数辨识方法
CN112653342B (zh) 一种静止坐标系下的复矢量电流环解耦控制装置及方法
CN114079399B (zh) 基于线性自抗扰控制的并网逆变器电流环控制系统和方法
CN112994482B (zh) 一种双pwm变频系统的模型预测功率动态补偿控制方法
CN107425743A (zh) 一种基于预测偏差反馈矫正的并网逆变器mpc方法
CN103595280A (zh) 电压不平衡下太阳能发电系统逆变器无差拍控制方法
CN102969913B (zh) 主动前端整流器模型预测控制参数不匹配补偿方法
CN108631624B (zh) 一种基于三维调制的级联h桥整流器及其控制方法
CN114301361B (zh) 一种基于母线电流控制的无电解电容永磁同步电机驱动系统控制方法
CN216959694U (zh) 基于线性自抗扰控制的并网逆变器电流环控制系统
CN112350600B (zh) 基于扰动估计的模块化多电平换流器功率解耦控制方法
CN115955130A (zh) 一种单相九电平pwm整流器的分数阶滑模控制方法
CN115566717A (zh) 一种基于改进型粒子群算法扰动观测器复合滑模的逆变器控制策略
CN111756261B (zh) 一种pwm整流器控制方法和装置
CN110707949B (zh) 定频pwm整流器的控制方法
CN110048466B (zh) 一种双馈风力发电系统网侧变流器模型预测控制方法
CN109391165B (zh) 一种基于模块化多电平变换器环流模型的干扰补偿方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant