CN111313467A - 基于参数联合设计的lcl型逆变器并网装置及控制方法 - Google Patents
基于参数联合设计的lcl型逆变器并网装置及控制方法 Download PDFInfo
- Publication number
- CN111313467A CN111313467A CN202010176277.5A CN202010176277A CN111313467A CN 111313467 A CN111313467 A CN 111313467A CN 202010176277 A CN202010176277 A CN 202010176277A CN 111313467 A CN111313467 A CN 111313467A
- Authority
- CN
- China
- Prior art keywords
- res
- coefficient
- sin
- cos
- omega
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000013016 damping Methods 0.000 claims abstract description 42
- 238000005070 sampling Methods 0.000 claims abstract description 16
- 230000001629 suppression Effects 0.000 claims abstract description 16
- 239000003990 capacitor Substances 0.000 claims abstract description 9
- 230000003068 static effect Effects 0.000 claims description 14
- 230000009466 transformation Effects 0.000 claims description 10
- 239000011347 resin Substances 0.000 claims description 9
- 238000010586 diagram Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/01—Arrangements for reducing harmonics or ripples
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Inverter Devices (AREA)
Abstract
本发明公开了一种基于参数联合设计的LCL型逆变器并网装置及控制方法。该装置包括三电平逆变器、数字处理控制模块和驱动电路,数字处理控制模块包括采样单元、闭环控制单元、有源阻尼单元和正弦脉宽调制单元。方法为:选取LCL滤波器电容容值、网侧电感和逆变器侧电感感值;通过奈奎斯特稳定性分析,并结合谐振抑制条件,确定PR控制器比例参数和有源阻尼反馈系数的可行范围;在所确定的可行域中使用拉格朗日乘子法选取合适的PR控制器比例参数和有源阻尼反馈系数,使控制***获得较大的带宽以及良好的动态性能。本发明具有硬件成本低、控制准确、适用范围广的特点,可以有效抑制LCL滤波器谐振频率的谐波分量,并降低入网电流的畸变率。
Description
技术领域
本发明属于电力电子变换技术领域,特别是一种基于参数联合设计的LCL型逆变器并网装置及控制方法。
背景技术
LCL滤波器由于其结构简单,高频滤波性能好,输出谐波含量少等优点,在新能源分布式并网发电场合得到了广泛的应用。但由于LCL滤波器的固有特性,LCL滤波器的谐振特性会大大降低输出侧的电能质量。目前对于LCL滤波器谐振频率的问题,主要有两种解决方法:(1)采用外加硬件阻尼电路抑制LCL谐振;(2)采用软件控制方法抑制LCL谐振。由于第一种方法会增加硬件成本,故一般使用后一种方法。在理想电网条件下,现有的有源阻尼控制方法相对成熟,如基于陷波器的有源阻尼控制方法、基于滤波电容电流反馈的有源阻尼控制方法、基于多状态量混合反馈的有源阻尼控制方法等。然而在实际情况中,电网并不是理想的基波电网,电网中会存在低频谐波,在非理想电网条件下,需要多个PR控制器并联控制,这需要***具有足够的开环带宽,这给有源阻尼的LCL滤波器谐振抑制控制带来了困难。
发明内容
本发明的目的在于提供一种适用于非理想电网条件下基于参数联合设计的LCL型逆变器并网装置及控制方法,以在非理想电网条件下实现LCL滤波器的谐振抑制。
实现本发明目的的技术解决方案为:一种基于参数联合设计的LCL型逆变器并网装置,包括三电平逆变器、数字处理控制模块和驱动电路,其中三电平逆变器为LCL型NPC三电平逆变器,数字处理控制模块包括采样单元、闭环控制单元、有源阻尼单元和正弦脉宽调制单元;
采样单元分别采集LCL滤波器网侧的三相电压信号、LCL滤波器网侧的三相电流信号,传送至闭环控制单元;
闭环控制单元根据采集到的信号,经Clarke变换将静止abc坐标系下的网侧电压、网侧电流变换到静止αβ坐标系下;将静止αβ坐标系下网侧电流的α、β轴分量iα、iβ输入有源阻尼单元;
闭环控制单元和有源阻尼单元得到的调制波信号相加,送至正弦脉宽调制单元,正弦脉宽调制单元的输出端经过驱动电路接入三电平逆变器每相桥臂的各个开关管。
进一步地,所述数字处理控制模块为TMS320F28377D与EPM1270T芯片。
一种基于参数联合设计的LCL型逆变器并网控制方法,包括以下步骤:
步骤1、在每个开关周期内,数字处理控制模块的采样单元分别采集LCL滤波器的网侧电压信号ea、eb、ec和网侧电流信号ia、ib、ic;
步骤2、闭环控制单元根据步骤1中采集到的信号,经Clarke变换将静止abc坐标系下的网侧电压、网侧电流变换到静止αβ坐标系下;
步骤3、计算z域下加入有源阻尼环及PR控制器后***的开环传递函数,利用奈奎斯特判据对***进行稳定性分析;
步骤4、对LCL滤波器谐振频率谐波抑制条件进行分析;
步骤5、在满足稳定性条件以及谐振抑制条件的范围内,使用拉格朗日乘子法,选取使***带宽最大的PR控制器的比例系数Kp和有源阻尼反馈系数H,以获得优化的闭环控制***;
步骤6、以电流正弦化为目标计算电流给定,并以网侧电流作为反馈量与所求得的电流给定量相减,经比例谐振调节器后与有源阻尼环输出相加,再经过Clarke反变换,输出三相调制波信号;
步骤7、将步骤6所得的三相调制信号经正弦脉宽调制单元生成脉宽调制信号,该脉宽调制信号通过驱动电路控制逆变器开关管的工作状态。
进一步地,步骤3所述对***进行稳定性分析,具体如下:
***开环传递函数Gop(s)为:
其中,L1为LCL滤波器逆变器侧电感感值,L2为LCL滤波器网侧电感感值,C为LCL滤波器电容容值,Kp为PR控制器的比例系数,H为有源阻尼反馈系数,ωres为LCL滤波器谐振频率,Ts为***采样时间,ωres的表达式如下:
将z=sin(ωTs)-jcos(ωTs)代入上式,整理得控制***开环频率响应方程为:
其中A(ω)为***开环频率响应方程中分子实部的系数、B(ω)为***开环频率响应方程中分子虚部的系数、C(ω)为***开环频率响应方程中分母实部的无H项系数、D(ω)为***开环频率响应方程中分母实部的含H项系数、E(ω)为***开环频率响应方程中分母虚部的无H项系数、F(ω)为***开环频率响应方程中分母虚部的含H项系数;
A(ω)=Ts 2[(ωresTs-sin(ωresTs))cos(ωTs)4-(2ωresTs cos(ωresTs)-2sin(ωresTs))cos(ωTs)3-(6sin(ωTs)2-1)(ωresTs-sin(ωresTs))cos(ωTs)2+6sin(ωTs)2(ωresTscos(ωresTs)-sin(ωresTs))cos(ωTs)+sin(ωTs)2(sin(ωTs)2-1)(ωresTs-sin(ωresTs))]
B(ω)=Ts 2sin(ωTs)[4(ωresTs-sin(ωresTs))cos(ωTs)3-6(ωresTs cos(ωresTs)-sin(ωresTs))cos(ωTs)2-(4sin(ωTs)2-2)(ωresTs-sin(ωresTs))cos(ωTs)+2sin(ωTs)2(ωresTs cos(ωresTs)-sin(ωresTs))]
D(ω)=L1L2{(ωresTs-sin(ωresTs))cos(ωTs)4-2(ωresTs+ωresTscos(ωresTs)-2sin(ωresTs))cos(ωTs)3-2[3(ωresTs-sin(ωresTs))sin(ωTs)2-2ωresTs cos(ωresTs)-ωresTs+3sin(ωresTs)]cos(ωTs)2+2(ωresTscos(ωresTs)+ωresTs-2sin(ωresTs))(3sin(ωTs)2-1)cos(ωTs)+(ωresTs-sin(ωresTs))sin(ωTs)4-2(2ωresTscos(ωresTs)+ωresTs-3sin(ωresTs))sin(ωTs)2+ωresTs-sin(ωresTs)}
F(ω)=L1L2sin(ωTs){4(ωresTs-sin(ωresTs))cos(ωTs)3-6(ωresTscos(ωresTs)+ωresTs-2sin(ωresTs))cos(ωTs)2+4[(sin(ωTs)2-3)sin(ωresTs)-(sin(ωTs)2-2cos(ωresTs)当-1)ωresTs]cos(ωTs)-2(sin(ωTs)+1)[2sin(ωresTs)-ωresTs(cos(ωresTs)-1)](sin(ωTs)-1)}
ωres≤ωs/6时:
Kp和H的取值范围如下:
其中A(ωh)为当ω=ωh时***开环频率响应方程中分子实部的系数、B(ωh)为当ω=ωh时***开环频率响应方程中分子虚部的系数、C(ωh)为当ω=ωh时***开环频率响应方程中分母实部的无H项系数、D(ωh)为当ω=ωh时***开环频率响应方程中分母实部的含H项系数、E(ωh)为当ω=ωh时***开环频率响应方程中分母虚部的无H项系数、F(ωh)为当ω=ωh时***开环频率响应方程中分母虚部的含H项系数;
当ωs/6<ωres<ωs/2时:
Kp和H的取值范围如下:
其中Kp为比例谐振控制器的比例系数,H为有源阻尼反馈系数,Ts为数字控制***采样周期,ωs为数字控制***采样角频率,ωh为***穿越频率,Top(z)为***开环传递函数。
进一步地,步骤4所述对LCL滤波器谐振频率谐波抑制条件进行分析,具体如下:
由于数字控制中零阶保持器的存在,LCL滤波器的实际的谐振频率ωres在使用有源阻尼控制环路后会发生偏移,偏移后的谐振频率为ωres',谐振抑制分析结果如下:
其中
a(ωres′)=-(D(ωres′)2+F(ωres′)2)·((A(ωres′)2+B(ωres′)2)·(D(ωres′)2+F(ωres′)2)′-(D(ωres′)2+F(ωres′)2)·(A(ωres′)2+B(ωres′)2)′)
b(ωres′)=4·(A(ωres′)2+B(ωres′)2)′·(C(ωres′)D(ωres′)3+C(ωres′)D(ωres′)F(ωres′)2+D(ωres′)2E(ωres′)F(ωres′)+E(ωres′)F(ωres′)3)-2·(A(ωres′)2+B(ωres′)2)·(C(ωres′)D(ωres′)3+C(ωres′)D(ωres′)F(ωres′)2+D(ωres′)2E(ωres′)F(ωres′)+E(ωres′)F(ωres′)3)′
c(ωres′)=2(A(ωres′)2+B(ωres′)2)′·(3C(ωres′)2D(ωres′)2+C(ωres′)2F(ωres′)2+D(ωres′)2E(ωres′)2+3E(ωres′)F(ωres′)+4C(ωres′)D(ωres′)E(ωres′)F(ωres′))-(A(ωres′)2+B(ωres′)2)·(3C(ωres′)2D(ωres′)2+C(ωres′)2F(ωres′)2+D(ωres′)2E(ωres′)2+3E(ωres′)F(ωres′)+4C(ωres′)D(ωres′)E(ωres′)F(ωres′))′
d(ωres′)=4·(A(ωres′)2+B(ωres′)2)′·(C(ωres′)3D(ωres′)+C(ωres′)D(ωres′)E(ωres′)2+C(ωres′)2E(ωres′)F(ωres′)+E(ωres′)3F(ωres′))-2·(A(ωres′)2+B(ωres′)2)·(C(ωres′)3D(ωres′)+C(ωres′)D(ωres′)E(ωres′)2+C(ωres′)2E(ωres′)F(ωres′)+E(ωres′)3F(ωres′))′
e(ωres′)=-(C(ωres′)2+E(ωres′)2)·((A(ωres′)2+B(ωres′)2)·(C(ωres′)2+E(ωres′)2)′-(C(ωres′)2+E(ωres′)2)·(A(ωres′)2+B(ωres′)2)′)
其中a(ωres')为当ω=ωres'时约束条件中H4项的系数,b(ωres')为当ω=ωres'时约束条件中H3项的系数,c(ωres')为当ω=ωres'时约束条件中H2项的系数,d(ωres')为当ω=ωres'时约束条件中H项的系数,e(ωres')为当ω=ωres'时约束条件中常数项的系数,A(ωres')为当ω=ωres'时***开环频率响应方程中分子实部的系数、B(ωres')为当ω=ωres'时***开环频率响应方程中分子虚部的系数、C(ωres')为当ω=ωres'时***开环频率响应方程中分母实部的无H项系数、D(ωres')为当ω=ωres'时***开环频率响应方程中分母实部的含H项系数、E(ωres')为当ω=ωres'时***开环频率响应方程中分母虚部的无H项系数、F(ωres')为当ω=ωres'时***开环频率响应方程中分母虚部的含H项系数。
进一步地,步骤5中Kp和H的选取如下:
max Kp
h2(Kp,H,ωh,ωres′)=a(ωres′)H4+b(ωres′)H3+c(ωres′)H2+d(ωres′)H+e(ωres′)=0
上式为KP和H的最终参数选取方程,其中gi(kp,H,ωh,ωres')且i=1,2,3表示不等式约束条件,hj(kp,H,ωh,ωres')表示等式约束条件且j=1,2。
本发明与现有技术相比,其显著优点在于:(1)通过网侧电流反馈有源阻尼,没有增加硬件成本,实现了LCL谐振抑制控制;(2)选取使***通带较大的控制参数,降低了输出电流的畸变率,提高了波形质量。
附图说明
图1是本发明基于参数联合设计的LCL型逆变器并网装置的结构示意图。
图2是LCL型NPC三电平逆变器并网***控制框图。
图3是NPC三电平并网逆变器的拓扑图。
图4是ωres≤ωs/6时加入有源阻尼之后的仿真波形图。
图5是ωs/6<ωres<ωs/2时加入有源阻尼之后的仿真波形图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
结合图1,本发明基于参数联合设计的LCL型逆变器并网装置,包括三电平逆变器、数字处理控制模块和驱动电路,其中三电平逆变器为LCL型NPC三电平逆变器,数字处理控制模块包括采样单元、闭环控制单元、有源阻尼单元和正弦脉宽调制单元;采样单元分别采集LCL滤波器网侧的三相电压信号、LCL滤波器网侧的三相电流信号,传送至闭环控制单元;闭环控制单元根据采集到的信号,经Clarke变换将静止abc坐标系下的网侧电压、网侧电流变换到静止αβ坐标系下;将αβ坐标系下网侧电流的α、β轴分量iα、iβ输入有源阻尼单元;闭环控制单元和有源阻尼单元得到的调制波信号相加,送至正弦脉宽调制单元,正弦脉宽调制单元的输出端经过驱动电路接入三电平逆变器每相桥臂的各个开关管。
作为一种具体示例,所述数字处理控制模块为TMS320F28377D与EPM1270T芯片。
本发明基于参数联合设计的LCL型逆变器并网装置的控制方法,包括以下步骤:
步骤1、在每个开关周期内,数字控制模块的采样单元分别采集LCL滤波器的网侧电压信号ea、eb、ec,网侧电流信号ia、ib、ic,直流侧上电容电压UC1和直流侧下电容电压UC2;
步骤2、闭环控制单元根据步骤1中采集到的信号,经Clarke变换将静止abc坐标系下的网侧电压、网侧电流变换到静止αβ坐标系下;
Clarke变换转换矩阵为Tabc/αβ:
经过此步骤,得到静止αβ坐标系下网侧电压的α、β轴分量eα、eβ和网侧电流的α、β轴分量iα、iβ;
步骤3、计算z域下加入有源阻尼环及PR控制器后***的开环传递函数,利用奈奎斯特判据对***进行稳定性分析。
LCL型NPC三电平逆变器并网***控制框图如图2所示,其中:
Gc(s)为电流控制器,其传递函数如下所示:
Kp为比例控制器增益,Kr为基波谐振控制器增益,ωi为基波谐振控制器角频率,Krh(h=3,5,7…)为各谐振控制器增益,h为谐波次数,ωih(h=3,5,7…)为各谐振控制器角频率,ωo为电网基波电压角频率。其中PR控制器可以实现对基波电流的无静差控制,HC控制器可以抑制并网电流中由电网背景谐波产生的谐波分量。
Gad(s)为有源阻尼环节,其传递函数为:
Gad(s)=HL1L2s2
其中,H为有源阻尼系数,L1为LCL滤波器逆变器侧电感感值,L2为LCL滤波器网侧电感感值。
GZOH(s)为零阶保持器,其传递函数为:
其中Ts为***采样时间。
G1(s)为LCL型逆变器的逆变器输出电压至入网电流的传递函数:
其中,C为LCL滤波器电容容值,ωres为LCL滤波器谐振频率,ωres的表达式如下:
***开环传递函数Gop(s)为:
其中,L1为LCL滤波器逆变器侧电感感值,L2为LCL滤波器网侧电感感值,C为LCL滤波器电容容值,Kp为PR控制器的比例系数,H为有源阻尼反馈系数,Ts为***采样时间。
将z=sin(ωTs)-jcos(ωTs)代入上式,整理可得控制***开环频率响应方程为:
其中A(ω)为***开环频率响应方程中分子实部的系数、B(ω)为***开环频率响应方程中分子虚部的系数、C(ω)为***开环频率响应方程中分母实部的无H项系数、D(ω)为***开环频率响应方程中分母实部的含H项系数、E(ω)为***开环频率响应方程中分母虚部的无H项系数、F(ω)为***开环频率响应方程中分母虚部的含H项系数。
A(ω)=Ts 2[(ωresTs-sin(ωresTs))cos(ωTs)4-(2ωresTs cos(ωresTs)-2sin(ωresTs))cos(ωTs)3-(6sin(ωTs)2-1)(ωresTs-sin(ωresTs))cos(ωTs)2+6sin(ωTs)2(ωresTscos(ωresTs)-sin(ωresTs))cos(ωTs)+sin(ωTs)2(sin(ωTs)2-1)(ωresTs-sin(ωresTs))]
B(ω)=Ts 2sin(ωTs)[4(ωresTs-sin(ωresTs))cos(ωTs)3-6(ωresTs cos(ωresTs)-sin(ωresTs))cos(ωTs)2-(4sin(ωTs)2-2)(ωresTs-sin(ωresTs))cos(ωTs)+2sin(ωTs)2(ωresTs cos(ωresTs)-sin(ωresTs))]
D(ω)=L1L2{(ωresTs-sin(ωresTs))cos(ωTs)4-2(ωresTs+ωresTs cos(ωresTs)-2sin(ωresTs))cos(ωTs)3-2[3(ωresTs-sin(ωresTs))sin(ωTs)2-2ωresTs cos(ωresTs)-ωresTs+3sin(ωresTs)]cos(ωTs)2+2(ωresTs cos(ωresTs)+ωresTs-2sin(ωresTs))(3sin(ωTs)2-1)cos(ωTs)+(ωresTs-sin(ωresTs))sin(ωTs)4-2(2ωresTs cos(ωresTs)+ωresTs-3sin(ωresTs))sin(ωTs)2+ωresTs-sin(ωresTs)}
F(ω)=L1L2sin(ωTs){4(ωresTs-sin(ωresTs))cos(ωTs)3-6(ωresTscos(ωresTs)+ωresTs-2sin(ωresTs))cos(ωTs)2+4[(sin(ωTs)2-3)sin(ωresTs)-(sin(ωTs)2-2cos(ωresTs)-1)ωresTs]cos(ωTs)-2(sin(ωTs)+1)[2sin(ωresTs)-ωresTs(cos(ωresTs)-1)](sin(ωTs)-1)}
(1)当ωres≤ωs/6时
当H<0时,不能满足稳定性条件;
当0≤H≤[2cos(ωresTs)-1]ωres 3Ts 2C/[sin(ωresTs)-ωresTs(2cos(ωresTs)-1)]时,为使***满足奈奎斯特稳定性判据,则需要满足|Gop(z=ejωhTs)|<1,即Kp和H需要满足如下条件:
其中A(ωh)为当ω=ωh时***开环频率响应方程中分子实部的系数、B(ωh)为当ω=ωh时***开环频率响应方程中分子虚部的系数、C(ωh)为当ω=ωh时***开环频率响应方程中分母实部的无H项系数、D(ωh)为当ω=ωh时***开环频率响应方程中分母实部的含H项系数、E(ωh)为当ω=ωh时***开环频率响应方程中分母虚部的无H项系数、F(ωh)为当ω=ωh时***开环频率响应方程中分母虚部的含H项系数。
当[2cos(ωresTs)-1]ωres 3Ts 2C/[sin(ωresTs)-ωresTs(2cos(ωresTs)-1)]<H时,不能满足稳定性条件;
(2)当ωs/6<ωres<ωs/2时
当H<[2cos(ωresTs)-1]ωres 3Ts 2C/[sin(ωresTs)-ωresTs(2cos(ωresTs)-1)]时,不能满足稳定性条件;
当[2cos(ωresTs)-1]ωres 3Ts 2C/[sin(ωresTs)-ωresTs(2cos(ωresTs)-1)]≤H≤0时,为使***满足奈奎斯特稳定性判据,则需要满足|Gop(z=ejωhTs)|<1,即Kp需要满足如下条件:
当0<H时,不能满足稳定性条件;
步骤4、对LCL滤波器谐振频率谐波抑制条件进行分析;
LCL滤波器谐振频率的谐波抑制条件为20lgAωres'≤-3dB,ωres′为加入有源阻尼后LCL滤波器的谐振频率,Aωres'的表达式为:
即Kp和H需要满足如下条件:
其中
a(ωres′)=-(D(ωres′)2+F(ωres′)2)·((A(ωres′)2+B(ωres′)2)·(D(ωres′)2+F(ωres′)2)′-(D(ωres′)2+F(ωres′)2)·(A(ωres′)2+B(ωres′)2)′)
b(ωres′)=4·(A(ωres′)2+B(ωres′)2)′·(C(ωres′)D(ωres′)3+C(ωres′)D(ωres′)F(ωres′)2+D(ωres′)2E(ωres′)F(ωres′)+E(ωres′)F(ωres′)3)-2·(A(ωres′)2+B(ωres′)2)·(C(ωres′)D(ωres′)3+C(ωres′)D(ωres′)F(ωres′)2+D(ωres′)2E(ωres′)F(ωres′)+E(ωres′)F(ωres′)3)′
c(ωres′)=2(A(ωres′)2+B(ωres′)2)′·(3C(ωres′)2D(ωres′)2+C(ωres′)2F(ωres′)2+D(ωres′)2E(ωres′)2+3E(ωres′)F(ωres′)+4C(ωres′)D(ωres′)E(ωres′)F(ωres′))-(A(ωres′)2+B(ωres′)2)·(3C(ωres′)2D(ωres′)2+C(ωres′)2F(ωres′)2+D(ωres′)2E(ωres′)2+3E(ωres′)F(ωres′)+4C(ωres′)D(ωres′)E(ωres′)F(ωres′))′
d(ωres′)=4·(A(ωres′)2+B(ωres′)2)′·(C(ωres′)3D(ωres′)+C(ωres′)D(ωres′)E(ωres′)2+C(ωres′)2E(ωres′)F(ωres′)+E(ωres′)3F(ωres′))-2·(A(ωres′)2+B(ωres′)2)·(C(ωres′)3D(ωres′)+C(ωres′)D(ωres′)E(ωres′)2+C(ωres′)2E(ωres′)F(ωres′)+E(ωres′)3F(ωres′))′
e(ωres′)=-(C(ωres′)2+E(ωres′)2)·((A(ωres′)2+B(ωres′)2)·(C(ωres′)2+E(ωres′)2)′-(C(ωres′)2+E(ωres′)2)·(A(ωres′)2+B(ωres′)2)′)
其中a(ωres')为当ω=ωres'时约束条件中H4项的系数,b(ωres')为当ω=ωres'时约束条件中H3项的系数,c(ωres')为当ω=ωres'时约束条件中H2项的系数,d(ωres')为当ω=ωres'时约束条件中H项的系数,e(ωres')为当ω=ωres'时约束条件中常数项的系数,A(ωres')为当ω=ωres'时***开环频率响应方程中分子实部的系数、B(ωres')为当ω=ωres'时***开环频率响应方程中分子虚部的系数、C(ωres')为当ω=ωres'时***开环频率响应方程中分母实部的无H项系数、D(ωres')为当ω=ωres'时***开环频率响应方程中分母实部的含H项系数、E(ωres')为当ω=ωres'时***开环频率响应方程中分母虚部的无H项系数、F(ωres')为当ω=ωres'时***开环频率响应方程中分母虚部的含H项系数。
步骤5、在满足稳定性条件以及谐振抑制条件的范围内,使用拉格朗日乘子法,选取使***带宽最大的PR控制器比例环节系数Kp和有源阻尼反馈系数H,以获得具有良好动态性能的闭环控制***,即求取满足如下条件的Kp和H:
max Kp
h2(Kp,H,ωh,ωres′)=a(ωres′)H4+b(ωres′)H3+c(ωres′)H2+d(ωres′)H+e(ωres′)=0
上式为KP和H的最终参数选取方程,其中gi(kp,H,ωh,ωres')(i=1,2,3)表示不等式约束条件,hj(kp,H,ωh,ωres')(j=1,2)表示等式约束条件。
步骤6、以电流正弦化为目标,计算电流给定,并以网侧电流作为反馈量与所求得的电流给定量相减,经比例谐振调节器后与有源阻尼环输出相加,再经过Clarke反变换,输出三相调制波信号;
步骤6.1、通过闭环控制单元和有源阻尼单元,得到静止αβ坐标系下的4路调制波信号uαh、uαpr、uβh、uβpr,将静止αβ坐标系下的α轴下的两个调制波信号uαh、uαpr相加得到:
uα=uαh+uαpr
β轴下的两个调制波信号uβh、uβpr相加得到:
uβ=uβh+uβpr
经过此步骤,得到静止αβ坐标系下的调制波信号uα、uβ;
步骤6.2、将静止αβ坐标系下uα、uβ转换到三相静止坐标系下,转换矩阵为Tαβ/abc:
经过此步骤,得到三相静止坐标系下的三相调制波信号ua、ub、uc;
步骤7、将步骤6所得的三相调制信号经正弦脉宽调制单元生成脉宽调制信号,该脉宽调制信号通过驱动电路控制逆变器开关管的工作状态,具体为:
将三相静止坐标系下的三相调制波信号ua、ub、uc送至正弦脉宽调制单元,生成脉宽调制信号,该脉宽调制信号通过驱动电路控制三电平逆变器开关管的工作状态,实现LCL滤波器谐振抑制的控制。
NPC三相三电平逆变器的调制规则如图3所示,以a相桥臂为例,在uaref的正半周,当uaref大于载波时,令Sa1、Sa2导通,a相桥臂输出高电平,当uaref小于载波时,令Sa2、Sa3导通,a相桥臂输出零电平;在uaref的负半周,当uaref小于载波时,令Sa3、Sa4导通,a相桥臂输出低电平,当uaref大于载波时,令Sa2、Sa3导通,a相桥臂输出零电平;b、c相桥臂的调制规则相同。
实施例1
本实施例利用MATLAB中的Simulink工具搭建了三电平逆变电路,直流电经直流母线电容后由三电平逆变电路逆变输出三相电压,经LCL滤波电路输出平稳的三相正弦电压。
仿真过程中的电气参数设置如表1:
表1
定义拉格朗日函数如下:
其中ωi是不等式约束的拉格朗日乘子,vj是等式约束的拉格朗日乘子,gi(Kp,H,ωh,ωres')表示不等式约束条件,hj(Kp,H,ωh,ωres')表示等式约束条件,gi(Kp,H,ωh,ωres')和hj(Kp,H,ωh,ωres')可以表示为:
h2(Kp,H,ωh,ωres′)=a(ωres′)H4+b(ωres′)H3+c(ωres′)H2+d(ωres′)H+e(ωres′)=0
将L(Kp,H,ωh,ωres',wi,vj)和L(Kp,H,ωh,ωres',wi,vj)分别定义为拉格朗日函数L(Kp,H,ωh,ωres',wi,vj)的一阶和二阶梯度,最大值的必要条件和充分条件为:
其中,d为拉格朗日函数的下降方向,G(Kp,H,ωh,ωres')为所有可行下降方向所形成的集合,G(Kp,H,ωh,ωres')的表达式为:
根据以上求解过程,电流控制器的比例增益在如下点达到最大值:
图4是ωres≤ωs/6时加入有源阻尼之后的仿真波形图,图5是ωs/6<ωres<ωs/2时加入有源阻尼之后的仿真波形图,可以看出,使用本发明中的一种基于参数联合设计的LCL型逆变器并网装置及方法,有效的抑制了网侧电流中的谐振频率次谐波,降低了电流的总谐波畸变率。
综上所述,本发明基于参数联合设计的LCL型逆变器并网控制方法,利用奈奎斯特判据对此***进行稳定性分析,计算抑制LCL滤波器谐振所需满足的条件,通过参数联合设计方法和拉格朗日乘子法得到可行域中使截止频率最大的参数作为***参数,将有源阻尼单元输出与闭环控制单元输出相加,再经Clarke反变换后得到三相调制波,经正弦脉宽调制单元生成正弦脉宽调制信号,该正弦脉宽调制信号通过驱动电路控制三电平逆变器各个开关管的工作状态,实现LCL滤波器谐振抑制的控制。本发明通过网侧电流反馈有源阻尼抑制了LCL谐振频率次谐波,减少了输出电流的谐波,提高了波形质量,有利于并网逆变器的并网,具有重大的工程应用价值。
Claims (6)
1.一种基于参数联合设计的LCL型逆变器并网装置,其特征在于,包括三电平逆变器、数字处理控制模块和驱动电路,其中三电平逆变器为LCL型NPC三电平逆变器,数字处理控制模块包括采样单元、闭环控制单元、有源阻尼单元和正弦脉宽调制单元;
采样单元分别采集LCL滤波器网侧的三相电压信号、LCL滤波器网侧的三相电流信号,传送至闭环控制单元;
闭环控制单元根据采集到的信号,经Clarke变换将静止abc坐标系下的网侧电压、网侧电流变换到静止αβ坐标系下;将静止αβ坐标系下网侧电流的α、β轴分量iα、iβ输入有源阻尼单元;
闭环控制单元和有源阻尼单元得到的调制波信号相加,送至正弦脉宽调制单元,正弦脉宽调制单元的输出端经过驱动电路接入三电平逆变器每相桥臂的各个开关管。
2.根据权利要求1所述的基于参数联合设计的LCL型逆变器并网装置,其特征在于,所述数字处理控制模块为TMS320F28377D与EPM1270T芯片。
3.一种基于参数联合设计的LCL型逆变器并网控制方法,其特征在于,包括以下步骤:
步骤1、在每个开关周期内,数字处理控制模块的采样单元分别采集LCL滤波器的网侧电压信号ea、eb、ec和网侧电流信号ia、ib、ic;
步骤2、闭环控制单元根据步骤1中采集到的信号,经Clarke变换将静止abc坐标系下的网侧电压、网侧电流变换到静止αβ坐标系下;
步骤3、计算z域下加入有源阻尼环及PR控制器后***的开环传递函数,利用奈奎斯特判据对***进行稳定性分析;
步骤4、对LCL滤波器谐振频率谐波抑制条件进行分析;
步骤5、在满足稳定性条件以及谐振抑制条件的范围内,使用拉格朗日乘子法,选取使***带宽最大的PR控制器的比例系数Kp和有源阻尼反馈系数H,以获得优化的闭环控制***;
步骤6、以电流正弦化为目标计算电流给定,并以网侧电流作为反馈量与所求得的电流给定量相减,经比例谐振调节器后与有源阻尼环输出相加,再经过Clarke反变换,输出三相调制波信号;
步骤7、将步骤6所得的三相调制信号经正弦脉宽调制单元生成脉宽调制信号,该脉宽调制信号通过驱动电路控制逆变器开关管的工作状态。
4.根据权利要求3所述的基于参数联合设计的LCL型逆变器并网控制方法,其特征在于,步骤3所述对***进行稳定性分析,具体如下:
***开环传递函数Gop(s)为:
其中,L1为LCL滤波器逆变器侧电感感值,L2为LCL滤波器网侧电感感值,C为LCL滤波器电容容值,Kp为PR控制器的比例系数,H为有源阻尼反馈系数,ωres为LCL滤波器谐振频率,Ts为***采样时间,ωres的表达式如下:
将z=sin(ωTs)-jcos(ωTs)代入上式,整理得控制***开环频率响应方程为:
其中A(ω)为***开环频率响应方程中分子实部的系数、B(ω)为***开环频率响应方程中分子虚部的系数、C(ω)为***开环频率响应方程中分母实部的无H项系数、D(ω)为***开环频率响应方程中分母实部的含H项系数、E(ω)为***开环频率响应方程中分母虚部的无H项系数、F(ω)为***开环频率响应方程中分母虚部的含H项系数;
A(ω)=Ts 2[(ωresTs-sin(ωresTs))cos(ωTs)4-(2ωresTscos(ωresTs)-2sin(ωresTs))cos(ωTs)3-(6sin(ωTs)2-1)(ωresTs-sin(ωresTs))cos(ωTs)2+6sin(ωTs)2(ωresTscos(ωresTs)-sin(ωresTs))cos(ωTs)+sin(ωTs)2(sin(ωTs)2-1)(ωresTs-sin(ωresTs))]
B(ω)=Ts 2sin(ωTs)[4(ωresTs-sin(ωresTs))cos(ωTs)3-6(ωresTscos(ωresTs)-sin(ωresTs))cos(ωTs)2-(4sin(ωTs)2-2)(ωresTs-sin(ωresTs))cos(ωTs)+2sin(ωTs)2(ωresTscos(ωresTs)-sin(ωresTs))]
D(ω)=L1L2{(ωresTs-sin(ωresTs))cos(ωTs)4-2(ωresTs+ωresTscos(ωresTs)-2sin(ωresTs))cos(ωTs)3-2[3(ωresTs-sin(ωresTs))sin(ωTs)2-2ωresTscos(ωresTs)-ωresTs+3sin(ωresTs)]cos(ωTs)2+2(ωresTscos(ωresTs)+ωresTs-2sin(ωresTs))(3sin(ωTs)2-1)cos(ωTs)+(ωresTs-sin(ωresTs))sin(ωTs)4-2(2ωresTscos(ωresTs)+ωresTs-3sin(ωresTs))sin(ωTs)2+ωresTs-sin(ωresTs)}
F(ω)=L1L2sin(ωTs){4(ωresTs-sin(ωresTs))cos(ωTs)3-6(ωresTscos(ωresTs)+ωresTs-2sin(ωresTs))cos(ωTs)2+4[(sin(ωTs)2-3)sin(ωresTs)-(sin(ωTs)2-2cos(ωresTs)-1)ωresTs]cos(ωTs)-2(sin(ωTs)+1)[2sin(ωresTs)-ωresTs(cos(ωresTs)-1)](sin(ωTs)-1)}
当ωres≤ωs/6时:
Kp和H的取值范围如下:
其中A(ωh)为当ω=ωh时***开环频率响应方程中分子实部的系数、B(ωh)为当ω=ωh时***开环频率响应方程中分子虚部的系数、C(ωh)为当ω=ωh时***开环频率响应方程中分母实部的无H项系数、D(ωh)为当ω=ωh时***开环频率响应方程中分母实部的含H项系数、E(ωh)为当ω=ωh时***开环频率响应方程中分母虚部的无H项系数、F(ωh)为当ω=ωh时***开环频率响应方程中分母虚部的含H项系数;
当ωs/6<ωres<ωs/2时:
Kp和H的取值范围如下:
其中Kp为比例谐振控制器的比例系数,H为有源阻尼反馈系数,Ts为数字控制***采样周期,ωs为数字控制***采样角频率,ωh为***穿越频率,Top(z)为***开环传递函数。
5.根据权利要求3所述的基于参数联合设计的LCL型逆变器并网控制方法,其特征在于,步骤4所述对LCL滤波器谐振频率谐波抑制条件进行分析,具体如下:
由于数字控制中零阶保持器的存在,LCL滤波器的实际的谐振频率ωres在使用有源阻尼控制环路后会发生偏移,偏移后的谐振频率为ωres',谐振抑制分析结果如下:
其中
a(ωres′)=-(D(ωres′)2+F(ωres′)2)·((A(ωres′)2+B(ωres′)2)·(D(ωres′)2+F(ωres′)2)′-(D(ωres′)2+F(ωres′)2)·(A(ωres′)2+B(ωres′)2)′)
b(ωres′)=4·(A(ωres′)2+B(ωres′)2)′·(C(ωres′)D(ωres′)3+C(ωres′)D(ωres′)F(ωres′)2+D(ωres′)2E(ωres′)F(ωres′)+E(ωres′)F(ωres′)3)-2·(A(ωres′)2+B(ωres′)2)·(C(ωres′)D(ωres′)3+C(ωres′)D(ωres′)F(ωres′)2+D(ωres′)2E(ωres′)F(ωres′)+E(ωres′)F(ωres′)3)′
c(ωres′)=2(A(ωres′)2+B(ωres′)2)′·(3C(ωres′)2D(ωres′)2+C(ωres′)2F(ωres′)2+D(ωres′)2E(ωres′)2+3E(ωres′)F(ωres′)+4C(ωres′)D(ωres′)E(ωres′)F(ωres′))-(A(ωres′)2+B(ωres′)2)·(3C(ωres′)2D(ωres′)2+C(ωres′)2F(ωres′)2+D(ωres′)2E(ωres′)2+3E(ωres′)F(ωres′)+4C(ωres′)D(ωres′)E(ωres′)F(ωres′))′
d(ωres′)=4·(A(ωres′)2+B(ωres′)2)′·(C(ωres′)3D(ωres′)+C(ωres′)D(ωres′)E(ωres′)2+C(ωres′)2E(ωres′)F(ωres′)+E(ωres′)3F(ωres′))-2·(A(ωres′)2+B(ωres′)2)·(C(ωres′)3D(ωres′)+C(ωres′)D(ωres′)E(ωres′)2+C(ωres′)2E(ωres′)F(ωres′)+E(ωres′)3F(ωres′))′
e(ωres′)=-(C(ωres′)2+E(ωres′)2)·((A(ωres′)2+B(ωres′)2)·(C(ωres′)2+E(ωres′)2)′-(C(ωres′)2+E(ωres′)2)·(A(ωres′)2+B(ωres′)2)′)
其中a(ωres')为当ω=ωres'时约束条件中H4项的系数,b(ωres')为当ω=ωres'时约束条件中H3项的系数,c(ωres')为当ω=ωres'时约束条件中H2项的系数,d(ωres')为当ω=ωres'时约束条件中H项的系数,e(ωres')为当ω=ωres'时约束条件中常数项的系数,A(ωres')为当ω=ωres'时***开环频率响应方程中分子实部的系数、B(ωres')为当ω=ωres'时***开环频率响应方程中分子虚部的系数、C(ωres')为当ω=ωres'时***开环频率响应方程中分母实部的无H项系数、D(ωres')为当ω=ωres'时***开环频率响应方程中分母实部的含H项系数、E(ωres')为当ω=ωres'时***开环频率响应方程中分母虚部的无H项系数、F(ωres')为当ω=ωres'时***开环频率响应方程中分母虚部的含H项系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010176277.5A CN111313467A (zh) | 2020-03-13 | 2020-03-13 | 基于参数联合设计的lcl型逆变器并网装置及控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010176277.5A CN111313467A (zh) | 2020-03-13 | 2020-03-13 | 基于参数联合设计的lcl型逆变器并网装置及控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111313467A true CN111313467A (zh) | 2020-06-19 |
Family
ID=71158727
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010176277.5A Pending CN111313467A (zh) | 2020-03-13 | 2020-03-13 | 基于参数联合设计的lcl型逆变器并网装置及控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111313467A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112003319A (zh) * | 2020-07-09 | 2020-11-27 | 西安理工大学 | 一种应用于双向并网变换器的双电流反馈控制方法 |
CN112532096A (zh) * | 2020-12-03 | 2021-03-19 | 南京理工大学 | 一种适应弱电网的lcl型逆变器并网装置及方法 |
CN113258615A (zh) * | 2021-06-08 | 2021-08-13 | 西南石油大学 | 并网逆变器频率自适应控制方法、装置、设备及存储介质 |
CN113629984A (zh) * | 2021-07-27 | 2021-11-09 | 浙江大学 | 一种基于双环电流控制策略的三相lcl型sapf参数设计方法 |
CN114094802A (zh) * | 2021-12-06 | 2022-02-25 | 南京理工大学 | 一种拓宽正阻尼区间的lcl型逆变器并网装置及方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109004649A (zh) * | 2018-06-29 | 2018-12-14 | 南京理工大学 | 一种基于有源阻尼的lcl滤波器谐振抑制装置及方法 |
-
2020
- 2020-03-13 CN CN202010176277.5A patent/CN111313467A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109004649A (zh) * | 2018-06-29 | 2018-12-14 | 南京理工大学 | 一种基于有源阻尼的lcl滤波器谐振抑制装置及方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112003319A (zh) * | 2020-07-09 | 2020-11-27 | 西安理工大学 | 一种应用于双向并网变换器的双电流反馈控制方法 |
CN112532096A (zh) * | 2020-12-03 | 2021-03-19 | 南京理工大学 | 一种适应弱电网的lcl型逆变器并网装置及方法 |
CN112532096B (zh) * | 2020-12-03 | 2022-09-23 | 南京理工大学 | 一种适应弱电网的lcl型逆变器并网装置及方法 |
CN113258615A (zh) * | 2021-06-08 | 2021-08-13 | 西南石油大学 | 并网逆变器频率自适应控制方法、装置、设备及存储介质 |
CN113258615B (zh) * | 2021-06-08 | 2021-10-08 | 西南石油大学 | 并网逆变器频率自适应控制方法、装置、设备及存储介质 |
CN113629984A (zh) * | 2021-07-27 | 2021-11-09 | 浙江大学 | 一种基于双环电流控制策略的三相lcl型sapf参数设计方法 |
CN113629984B (zh) * | 2021-07-27 | 2023-02-28 | 浙江大学 | 一种基于双环电流控制策略的三相lcl型sapf参数设计方法 |
CN114094802A (zh) * | 2021-12-06 | 2022-02-25 | 南京理工大学 | 一种拓宽正阻尼区间的lcl型逆变器并网装置及方法 |
CN114094802B (zh) * | 2021-12-06 | 2024-04-12 | 南京理工大学 | 一种拓宽正阻尼区间的lcl型逆变器并网装置及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109004649B (zh) | 一种基于有源阻尼的lcl滤波器谐振抑制装置及方法 | |
CN111313467A (zh) | 基于参数联合设计的lcl型逆变器并网装置及控制方法 | |
CN112532096B (zh) | 一种适应弱电网的lcl型逆变器并网装置及方法 | |
CN106787910A (zh) | 应用于并网逆变器电流控制的改进重复控制器设计方法 | |
CN110429603B (zh) | 六开关七电平有源电力滤波器及补偿方法 | |
CN112653342B (zh) | 一种静止坐标系下的复矢量电流环解耦控制装置及方法 | |
CN112636348B (zh) | 一种模块化三相电流型并网逆变器控制方法 | |
CN114583702B (zh) | Lcl型有源电力滤波器的并联电流控制***及控制方法 | |
CN112787491B (zh) | 一种三级式ac/dc电源的输入电流谐波抑制方法 | |
CN112165271A (zh) | 一种并网变流***及其模型预测控制方法 | |
CN110912150A (zh) | 一种基于虚拟阻抗的混合有源滤波器并网方法 | |
CN115276445A (zh) | 弱网下基于vsg的lcl并网逆变器谐振抑制及稳定性分析方法 | |
CN116683520A (zh) | 一种基于谐振扩展状态观测器的pcs自抗扰控制方法 | |
CN113067372A (zh) | 一种改善lcl滤波并网控制性能的有源阻尼方法及电路 | |
CN107528587B (zh) | 一种基于pir调节器的高精度快速宽频单相软锁相环 | |
CN116979535A (zh) | 一种用于有源电力滤波器的双重谐波电流检测方法 | |
Wei et al. | Performance analysis and improvement of output current controller for three-phase shunt active power filter | |
Pan et al. | LCL APF based on fractional-order fast repetitive control strategy | |
CN110707908A (zh) | 一种基于自适应电流谐波抑制的逆变器电流控制*** | |
CN210608909U (zh) | 一种基于自适应电流谐波抑制的逆变器电流控制*** | |
Tan et al. | Deadbeat plus fractional-order repetitive control for grid-tied inverter | |
Kuo et al. | FPGA control of a three-phase three-level T-type NPC grid-connected inverter with LCL filter | |
CN110880790B (zh) | 含lcl拓扑结构的并网功率变换器的控制方法 | |
CN112994496B (zh) | 一种任意有功工况下恒速率均压的模块化多电平变换装置 | |
CN111769563B (zh) | 一种lcl有源电力滤波器模型及其预测控制方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200619 |
|
RJ01 | Rejection of invention patent application after publication |