基于类数字陷波器的网侧电流反馈LCL型并网逆变器谐振的有源阻尼方法

摘要

本发明公开了一种基于类数字陷波器的网侧电流反馈LCL型并网逆变器谐振的有源阻尼方法，用于解决现有网侧电流反馈LCL型并网逆变器谐振的阻尼方法鲁棒性差的技术问题。技术方案是采用单环级联型阻尼方式，级联在主电流控制器之后，不需要额外的电压或电流传感器，不构成内环，简化了控制分析和设计，降低了成本。阻尼的方式是通过改变开环环路增益的相位，使得被修改的环路增益的相频提前穿越负180度，同时被修改的环路增益的相频超前部分被削减，整体相频被限制在负180度和负540度之间，使其能够阻尼从低频到奈奎斯特频率之间的宽范围的谐振频率，通用性良好，避免重复设计，同时面对LCL参数变化和网侧阻抗变动，其鲁棒性良好。

权利要求

1.一种基于类数字陷波器的网侧电流反馈LCL型并网逆变器谐振的有源阻尼方法，其特征在于包括以下步骤：步骤一、在离散时间域，类数字陷波器GNL的传递函数为式(1)，其相频 为式(2)，在数字控制器中采用直接数字II型实现；

式中，K＝(1+p0)2/(1+a1+a2)，p0,a1和a2是三个待定实系数，决定类数字陷波器的零极点，z是离散时间域z变换的算子，ω是频率变量，T是采样周期；

步骤二、预先定义位于临界频率fs/6之前的负180度相位穿越频率ωx，式中fs是采样频率；通过向原始环路增益OLTF中加入类数字陷波器，使得被修改的环路增益MOLTF的相位在ωx处穿越负180度，得到表达式(3)∠MOLTF＝GNLGcKPWMGp(z)z-1＝-π            (3)式中，∠MOLTF是被修改的环路增益的相频，Gc是电流控制器传递函数，KPWM是逆变器等效模型，Gp(z)是LCL滤波器低频模型，z-1是一个周期的控制延时；

步骤三、确定p0关于ωx的约束；参见式(1)，类数字陷波器有两个重实数极点，令其中一个极点在ωx处独立提供相位滞后，令另一个极点的相位滞后效果被两个零点抵消，表示为相位关系得到式(4)，将式(2)和式(4)代入式(3)，得到p0关于ωx的约束为式(5)；式中，∠表示求传递函数的相频；

∠(z+p0)＝∠(z2+a1z+a2)             (4)

步骤四、得到a2关于a1和p0的约束；化简式(4)，得到式(6)，要使得解有意义，式(6)分母不为0，因此进一步求解得到a2关于a1和p0的约束关系为式(7)；

a2＝a1p0+2p0cos(ωxT)+1＝a1p0+m              (7)式中，m＝2p0cos(ωxT)+1；

步骤五、求解a1；令环路增益的相位在两个频率点处相等，这两个频率点是类数字陷波器的相频过零点频率ωz即2πfz和临界频率fs/6；根据ωz的定义得到式(8)，化简得到式(9)；另外，根据LCL型并网逆变器的可控性要求，0<ωzT<π，则sin(ωzT)≠0，由此化简式(9)得到式(10)；同时，根据环路增益的相位在两个频率点处相等的条件得到式(11)；类数字陷波器在临界频率fs/6处的相频为式(12)；综合(10)-(12)三个方程得到一个可行的a1；

式中，k＝1，使得类数字陷波器相频位于(–π/2,π/2)之间；

至此，只要给定负180度相位穿越频率ωx，由式(5)、式(7)和式(10)-式(12)分别算出三个参数p0,a1和a2，够获得用于网侧电流反馈LCL型并网逆变器谐振阻尼的类数字陷波器：

说明书

基于类数字陷波器的网侧电流反馈LCL型并网逆变器谐振的

有源阻尼方法

技术领域

[0001] 本发明涉及一种网侧电流反馈LCL型并网逆变器谐振的阻尼方法，特别涉及一种基于类数字陷波器的网侧电流反馈LCL型并网逆变器谐振的有源阻尼方法。

背景技术

[0002] 参照图10。文献1“Design and Analysis of Robust Active Damping for LCL Filters Using Digital Notch Filters,IEEE Transactions on Power Electronics,vol.32,no.3,pp.2360–2375,Mar.2017.”公开了一种基于数字陷波器的网侧电流反馈LCL型并网逆变器谐振尖峰的阻尼方法，该方法是单环网侧电流控制，不需要额外的电压或电流传感器；通过数字陷波器引入了一定的相位滞后，使得开环环路增益在谐振频率之前穿越负180度，根据奈奎斯特稳定性判据，避免了幅频谐振尖峰对应相频负180度，从而保证了系统稳定性。OLTF是原始环路增益，Gnotch是数字陷波器，MOLTF是被修改的环路增益。对于低谐振频率fr＝1.37kHz和高谐振频率fr＝2.98kHz，数字陷波器保证了系统的稳定性。对于高于但接近陷波频率的谐振频率fr＝2.11kHz，由于数字陷波器在陷波频率处引入正180度相位超前，使得谐振频率fr＝2.11kHz时的环路增益存在有效负穿越，系统不稳定。因此该方法只能阻尼较低或较高的谐振频率，不能阻尼高于但接近陷波频率的谐振频率，这就造成通用性差和不必要的重复设计，同时也意味着对于器件老化等原因造成的参数变动的鲁棒性差。

[0003] 从图10可以看出，背景技术公开的数字陷波器的在陷波频率处有固定的正180度的超前相位，该超前相位只和陷波频率有关，不能任意调整，当谐振频率出现在陷波频率附近时，如谐振频率fr＝2.11kHz及陷波频率为1.95kHz时，数字陷波器提供的正180度的超前相位将造成有效负穿越，幅频谐振尖峰对应负180度，因此，数字陷波器无法有效阻尼高于且接近陷波频率的谐振频率。

发明内容

[0004] 为了克服现有网侧电流反馈LCL型并网逆变器谐振的阻尼方法鲁棒性差的不足，本发明提供一种基于类数字陷波器的网侧电流反馈LCL型并网逆变器谐振的有源阻尼方法。该方法采用单环级联型阻尼方式，级联在主电流控制器之后，不需要额外的电压或电流传感器，不构成内环，简化了控制分析和设计，降低了成本。阻尼的方式是通过改变开环环路增益的相位，使得被修改的环路增益的相频提前穿越负180度，同时被修改的环路增益的相频超前部分被削减，整体相频被限制在负180度和负540度之间，使其能够阻尼从低频到奈奎斯特频率之间的宽范围的谐振频率，通用性良好，避免重复设计，同时面对LCL参数变化和网侧阻抗变动，其鲁棒性良好。

[0005] 本发明解决其技术问题所采用的技术方案：一种基于类数字陷波器的网侧电流反馈LCL型并网逆变器谐振的有源阻尼方法，其特点是包括以下步骤：

[0006] 步骤一、在离散时间域，类数字陷波器GNL的传递函数为式(1)，其相频 为式(2)，在数字控制器中采用直接数字II型实现。

[0007]

[0008]

[0009] 式中，K＝(1+p0)2/(1+a1+a2)，p0,a1和a2是三个待定实系数，决定类数字陷波器的零极点，z是离散时间域z变换的算子，ω是频率变量，T是采样周期。

[0010] 步骤二、预先定义位于临界频率fs/6之前的负180度相位穿越频率ωx，式中fs是采样频率。通过向原始环路增益OLTF中加入类数字陷波器，使得被修改的环路增益MOLTF的相位在ωx处穿越负180度，得到表达式(3)

[0011] ∠MOLTF＝GNLGcKPWMGp(z)z-1＝-π           (3)

[0012] 式中，∠MOLTF是被修改的环路增益的相频，Gc是电流控制器传递函数，KPWM是逆变器等效模型，Gp(z)是LCL滤波器低频模型，z-1是一个周期的控制延时。

[0013] 步骤三、确定p0关于ωx的约束。参见式(1)，类数字陷波器有两个重实数极点，令其中一个极点在ωx处独立提供相位滞后，令另一个极点的相位滞后效果被两个零点抵消，表示为相位关系得到式(4)，将式(2)和式(4)代入式(3)，得到p0关于ωx的约束为式(5)。式中，∠表示求传递函数的相频。

[0014] ∠(z+p0)＝∠(z2+a1z+a2)@ωx             (4)

[0015]

[0016] 步骤四、得到a2关于a1和p0的约束。化简式(4)，得到式(6)，要使得解有意义，式(6)分母不为0，因此进一步求解得到a2关于a1和p0的约束关系为式(7)。

[0017]

[0018] a2＝a1p0+2p0cos(ωxT)+1＝a1p0+m           (7)

[0019] 式中，m＝2p0cos(ωxT)+1。

[0020] 步骤五、求解a1。令环路增益的相位在两个频率点处相等，这两个频率点是类数字陷波器的相频过零点频率ωz即2πfz和临界频率fs/6。根据ωz的定义得到式(8)，化简得到式(9)。另外，根据LCL型并网逆变器的可控性要求，0<ωzT<π，则sin(ωzT)≠0，由此化简式(9)得到式(10)。同时，根据环路增益的相位在两个频率点处相等的条件得到式(11)。类数字陷波器在临界频率fs/6处的相频为式(12)。综合(10)-(12)三个方程得到一个可行的a1。

[0021]

[0022]

[0023]

[0024]

[0025]

[0026] 式中，k＝1，使得类数字陷波器相频位于(–π/2,π/2)之间。

[0027] 至此，只要给定负180度相位穿越频率ωx，由式(5)、式(7)和式(10)-式(12)分别算出三个参数p0,a1和a2，够获得用于网侧电流反馈LCL型并网逆变器谐振阻尼的类数字陷波器：

[0028]

[0029] 本发明的有益效果是：该方法采用单环级联型阻尼方式，级联在主电流控制器之后，不需要额外的电压或电流传感器，不构成内环，简化了控制分析和设计，降低了成本。阻尼的方式是通过改变开环环路增益的相位，使得被修改的环路增益的相频提前穿越负180度，同时被修改的环路增益的相频超前部分被削减，整体相频被限制在负180度和负540度之间，使其能够阻尼从低频到奈奎斯特频率之间的宽范围的谐振频率，通用性良好，避免重复设计，同时面对LCL参数变化和网侧阻抗变动，其鲁棒性良好。

[0030] 下面结合附图和具体实施方式对本发明作详细说明。

附图说明

[0031] 图1是本发明方法中类数字陷波器的相频频率特性；

[0032] 图2是本发明方法中类数字滤波器在数字控制器中直接数字II型实现；

[0033] 图3是本发明方法实施例网侧电流反馈LCL型并网逆变器的单环控制结构，其中i2是网侧电流，是电流参考值，VO是逆变器桥臂输出电压；

[0034] 图4是本发明方法实施例的Bode图，其中GNL是类数字陷波器，OLTF是原始环路增益，MOLTF是被修改的环路增益；

[0035] 图5是本发明方法实施例被修改的环路增益的相位、频率和谐振频率的关系图；

[0036] 图6是本发明方法实施例逆变器侧电感L1和网侧阻抗Lg变动时闭环零极点图；

[0037] 图7是本发明方法实施例滤波电容C变动时闭环零极点图；

[0038] 图8是本发明方法实施例当滤波电容C变动时的仿真结果；

[0039] 图9是本发明方法实施例的动态性能仿真结果；

[0040] 图10是背景技术方法数字陷波器对不同谐振频率的有效和无效阻尼情况。

具体实施方式

[0041] 参照图1-9。本发明基于类数字陷波器的网侧电流反馈LCL型并网逆变器谐振的有源阻尼方法具体步骤如下：

[0042] 1、类数字陷波器GNL传递函数的结构和设计过程。

[0043] 类数字陷波器的传递函数和相频特性为：

[0044]

[0045]

[0046] 式中，GNL是类数字陷波器的传递函数， 是类数字陷波器的传递函数的相频，K＝(1+p0)2/(1+a1+a2)，p0,a1和a2是三个待定实系数，决定类数字陷波器的零极点，z是离散时间域z变换的算子，ω是频率变量，T是采样周期。

[0047] 类数字陷波器GNL的设计过程包括以下步骤：

[0048] 步骤一：预先定义负180度相位穿越频率ωx，本实施例取1600πrad/s。在引入类数字陷波器后，根据穿越频率定义，有：

[0049] ∠MOLTF＝GNLGcKPWMGpz-1＝-π              (3)

[0050] 式中，∠MOLTF是被修改的环路增益的相频，Gc是电流控制器传递函数，KPWM是逆变器等效模型，Gp是LCL滤波器的低频模型，z-1是数字控制中一个周期的控制延时。

[0051] 步骤二：类数字陷波器有两个重实数极点，令其中一个极点在ωx处独立提供相位滞后，令另一个极点的相位滞后效果被两个零点抵消，表示为相位关系为式(4)，将式(2)和(4)代入式(3)，得到式(5)，其中ωx和p0是单调关系，由式(5)计算出p0等于–0.3355。

[0052] ∠(z+p0)＝∠(z2+a1z+a2)@ωx               (4)

[0053]

[0054] 式中，∠表示求传递函数的相频。

[0055] 步骤三：对式(4)进行化简，得到式(6)，要使得解有意义，式(6)分母不为0，据此进一步求解得到式(7)。式(7)表明p0被式(5)约束之后，a1和a2是线性关系。

[0056]

[0057] a2＝a1p0+2p0cos(ωxT)+1＝a1p0+m                (7)

[0058] 式中，m＝2p0cos(ωxT)+1，代入数值得m＝0.4120。

[0059] 步骤四：

[0060] 参照图1。经过上述计算，a2被a1线性表出，则只剩下a1待定。对 求a1的导数，得到式(8)，根据LCL型并网逆变器的可控性要求，0<ωT<π，因此在ω>ωx时， 是a1的减函数，表明随着a1的增加， 单调减少，即 的过零点频率ωz(2πfz)向右移动。为验证，取a1为-1.5到1，步进0.5，相应的a2由式(7)约束，得到一系列类数字陷波器的相频特性，随着a1增加，其过零点逐渐右移，这表明类数字陷波器的超前相位可以通过a1来调整。

[0061]

[0062] 再令被修改的环路增益MOLTF的相位在两个频率点处相等，这两个频率点是类数字陷波器的相频过零点频率ωz即2πfz和临界频率fs/6。这样，根据ωz的定义得到式(9)，化简得到式(10)。另外，根据LCL型并网逆变器的可控性要求，0<ωzT<π，则sin(ωzT)≠0，由此化简式(10)得到式(11)。同时，根据环路增益的相位在两个频率点处相等的条件得到式(12)。类数字陷波器在临界频率fs/6处的相频为式(13)。

[0063]

[0064]

[0065]

[0066]

[0067]

[0068] 式中，k＝1，使得类数字陷波器相频在(–π/2,π/2)之间。

[0069] 综合式(11)-式(13)，调用数值分析软件MATALB提供的fsolve.m程序解得一个可行的a1等于–0.2391，根据式(7)得到对应的a2等于0.4997，K等于0.3503。此时，得到离散时间域类数字陷波器GNL的传递函数为式(14)：

[0070]

[0071] 参见图2。类数字陷波器在数字控制器中的实现为：

[0072]

[0073] 式中，v[n]是离散中间变量，x[n]是离散输入变量，y[n]是离散输出变量，n是数字控制的第n个采样时刻。

[0074] 2、被修改的环路增益的频率特性。

[0075] 为了清楚直观地解释说明提出的方法，这里给出通用的网侧电流反馈控制结构，并利用Bode图分析环路增益的频率特性，采用的参数见表1。

[0076] 表1 参数表

[0077]

[0078]

[0079] 参照图3。网侧电流反馈LCL型并网逆变器的控制模型包括电流控制器Gc，一个周期延时z-1，逆变器模型KPWM，LCL滤波器模型Gp。类数字陷波器GNL级联在Gc之后，用于调整环路增益的相频，被修改的环路增益MOLTF的表达式为(16)。

[0080] MOLTF＝GcGNLz-1KPWMGp                  (16)

[0081] 其低频近似模型的相频表达式为：

[0082]

[0083] 式中，sign为符号函数，定义为：

[0084]

[0085] 3、Bode分析和仿真验证。

[0086] 下面利用Bode图和闭环零极点分析和验证环路增益频率特性的改善，以及对参数变动鲁棒性的提高。

[0087] 参照图1。当不考虑谐振角频率ωr的作用时，原始相频是一条光滑的曲线，在加入类数字陷波器后，由于a1的不同取值，相频特性在穿越角频率ωx之后会出现大于负180度的超前区域和小于负540度的滞后区域。当谐振角频率ωr出现在这些区域时，被修改的环路增益存在有效负穿越，将导致系统不稳定。因此，本发明令被修改的环路增益MOLTF的相位在两个频率点fz和fs/6相等来约束a1，消除了超前和滞后区域，保证了稳定性。

[0088] 参照图4。根据上述，令被修改的环路增益MOLTF的相位在两个频率点fz和fs/6相等，由式(24)和式(25)理论计算的fz为2.75kHz，相位为–238.4–180＝–418.4度，和Bode分析得到的fz＝2.74kHz，相位负419度一致。

[0089] 参照图5。当谐振角频率从ωr扫过2000πrad/s到1×104πrad/s时，被修改的环路增益MOLTF的相位被限制在负180度和负540度之间。这样，即使谐振角频率ωr因为器件老化或者网侧阻抗变化时，被修改的环路增益也不存在有效负穿越，提高了对参数变动的鲁棒性。

[0090] 参照图6和图7。进一步验证鲁棒性。LCL滤波器参数和网侧阻抗变动时，参数变化设定为：L1从50％的0.9mH到150％的2.7mH，Lg从0mH到10mH，C从50％的1.5uF到470％的14.1uF，系统闭环极点位于单位圆内，系统稳定，因此，鲁棒性良好。

[0091] 利用MATLAB/Simulink仿真软件进行仿真验证。

[0092] 参照图8。仿真验证类数字陷波器对参数变动的鲁棒性。初始电容C值为14.1uF，对应谐振频率1.27kHz，系统稳定运行；在0.1秒时，从C变为3uF，谐振频率变为2.98kHz，经过一段时间瞬态过程，系统仍进入稳定状态；在0.2秒时，C再次变为14.1uF，系统仍能稳定，和Bode分析一致；在0.27秒移除类数字陷波器阻尼，由于谐振频率小于临界频率1.67kHz，所以系统发散。

[0093] 参照图9。给定LCL的参数，仿真验证类陷波器的动态性能。当网侧阻抗Lg＝0mH时，系统带宽约400Hz，在参考信号从4A条变到10A以及反过来跳变时，动态过程较快；当网侧阻抗Lg＝10mH时，由于比例系数恒定，系统带宽下降，所以动态过程减慢，但是系统仍然保持稳定，验证了Bode分析。另外，如果谐振频率小于临界频率1.67kHz，移除类数字陷波器会使得系统发散，当谐振频率大于临界频率1.67kHz，移除类数字陷波器，系统仍然稳定。