DOI： 10.11817/j.issn.1672-7207.2021.05.012

单相BES-quasi-ZSI二倍频脉动的解析建模

魏祥林^{1, 2, 3}，王晗⁴，陈伟^{2, 3}，杨维满^{2, 3}，黄苏融¹

(1. 上海大学 机电工程与自动化学院，上海，200444；

2. 兰州理工大学 电气工程与信息工程学院，甘肃 兰州，730050；

3. 甘肃省新能源电力变换与控制工程研究中心，甘肃 兰州，730050；

4. 上海交通大学 电子信息与电气工程学院，上海，200240)

摘要：研究单相Ⅱ型电池储能型准Z源逆变器(battery-energy stored quasi-Z-source inverter，BES-qZSI)的二倍频脉动产生和传播的物理机理，运用状态空间平均法求解其二倍频脉动量的精确解析表达式；在分析其二倍频脉动量对qZS网络阻抗参数敏感性的基础上，提出一种基于主导阻抗参数抑制二倍频脉动的qZS网络不对称阻抗参数的设计方法，建立基于Matlab/Simulink的BES-qZSI-II电路仿真模型，并将仿真结果与计算结果进行比较。研究结果表明：本文所建解析模型正确。

关键词：准Z源逆变器；电池储能；二倍频脉动；解析模型；不对称阻抗参数

中图分类号：TM464            文献标志码：A            开放科学(资源服务)标识码(OSID)

文章编号：1672-7207（2021）05-1502-14

Analytic modeling of double-line frequency ripple in single-phase BES-quasi-ZSI

WEI　Xianglin^{1, 2, 3}, WANG　Han⁴, CHEN　Wei^{2, 3}, YANG　Weiman^{2, 3}, HUANG　Surong¹

(1. School of Mechatronic Engineering and Automation, Shanghai University, Shanghai 200444, China;

2. College of Electrical and Information Engineering, Lanzhou University of Technology, Lanzhou 730050, China;

3. Gansu Renewable Energy Power Conversion and Control Engineering Technology Research Center, Lanzhou 730050, China;

4. School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai 200240, China)

Abstract: The generating and propagating mechanism of the double-line frequency(2 ω) ripple of the battery-energy stored quasi-Z-source inverter(BES-qZSI-II) was studied. The analytical model of the 2 ω ripple was established by applying the state space averaging method. Based on analyzing the sensitivity of impedance parameters on the 2 ω ripple, a design procedure for the asymmetric impedance parameters of single-phase BES-qZS-II inverter was proposed, which was based on the idea that the 2 ω ripple is mainly suppressed by the dominant impedance parameter accordingly. A circuit simulation model was established based on Matlab/Simulink, and the simulation results were compared with those of calculation. The results show that the proposed analytical model to calculate the 2 ω ripple in single-phase BES-qZSI-II inverter is correct.

Key words: quasi-Z-source inverter; battery-energy stored; double-line frequency ripple; analytic model; asymmetrical impedance parameter

准Z源逆变器(quasi-Z-source inverter, qZSI)可实现DC-DC与DC-AC功能的单级集成，具有效率高、成本低、可靠性好、输入电流连续、qZS网络电容器电压等级低、直流电源与直流母线共地等优点，被视为分布式发电系统的新生代电能变换器，尤其适用于光伏发电系统^[1-3]。但是，基于qZSI的光伏发电系统没有集成储能装置，当光照强度不足时不能提供持续稳定的电能。因此，为了补偿光伏板与负荷之间的不平衡功率，平缓光伏功率随机波动对负荷的冲击，提高光伏发电系统的鲁棒性，扩大其正常运行范围，带储能装置的qZSI光伏发电系统成为众多学者研究的热点^[4-6]。

文献[7]首次提出了无外加充放电电路，直接将储能电池并接在电容器C₂两端的储能型qZSI拓扑(第I类储能型qZSI，简称BES-qZSI-I)，并通过理论分析、仿真和实验验证了电池充放电过程与机理。文献[8]和[9]指出，断续工作模式(discontinuous conduction mode, DCM)会限制BES-qZSI-I拓扑中储能电池的放电能力，提出将储能电池与电容器C₁并联的储能型qZSI拓扑(第II类储能型qZSI，简称BES-qZSI-II)，证明了在逆变输出相同功率的前提下相对于BES-qZSI-I而言，该拓扑的储能电池具有更宽的放电功率范围，并提出了一种基于储能电池电流闭环控制实现恒定直流母线峰值电压的控制策略。文献[10]指出BES-qZSI-I和BES-qZSI-II存在同样的缺点，即电容器电压受电池电压箝位，限制了系统设计和直流电压控制的灵活性，进而提出一种通过双向DC-DC变换器将储能电池与电容器C₁并联的BES-qZSI拓扑，其电池电压变化范围可以更宽，且电网可向电池反向充电；但该拓扑因结构和控制复杂，未能成为BES-qZSI的主流拓扑。

然而单相BES-qZSI具有普通单相逆变器的固有特性——逆变输出功率含有二倍频脉动量^[11]。二倍频脉动功率从交流侧传播到直流侧，会在BES-qZSI的直流母线电压、qZS网络的电压和电流、光伏板电压和储能电池电压中激励起二倍频脉动分量^[12]，从而增加系统损耗、缩短储能电池和光伏板的寿命^[13]，以及增大输出电压的谐波含量^[14]，影响光伏板的最大功率跟踪(maximum power point tracing, MPPT)的精度，因此必须加以限制。单相qZSI二倍频脉动的抑制方法可分为被动抑制法和主动抑制法两大类^[12]。主动抑制法通过附加有源滤波电路(APF)^[15-17]、采用先进控制策略^[18-19]、修正PWM技术^[20-22]等手段抑制二倍频脉动，但其控制实现复杂。被动抑制法基于二倍频脉动与阻抗参数的关系，通过加大qZS网络的电容和电感抑制二倍频脉动，实现简单，其实现的关键是建立二倍频脉动量与qZS网络阻抗参数的关系表达式；此外，阻抗网络(Z网络)根据结构与参数的不同又可分为对称阻抗网络和不对称阻抗网络。针对对称阻抗网络，文献[20]推导了单相Z源逆变器电容电压和电感电流的低频脉动与阻抗参数的关系，但其解析方程复杂，不利于阻抗参数设计；文献[3]用一个五阶的交流信号模型分析qZSI的二倍频脉动，但二倍频脉动量之间存在相互耦合，难以清楚描述二倍频脉动与阻抗参数之间的关系。文献[23]和[24]分别建立了单相ZSI和qZSI的大信号平均值模型，并提出抑制二倍频脉动的阻抗参数设计方法，但是把光伏板建模成恒压源不能真实地反映其对变换器的动态作用；文献[25]考虑了光伏板的动力学特性及其端电容的滤波和稳压作用，建立了qZSI二倍频脉动的综合模型用于设计阻抗参数，但此模型是在阻抗参数对称的条件下建立的模型，且未考虑储能电池的影响。对称阻抗网络有利于简化qZSI的分析，降低其动态模型的阶数，但又限制了阻抗参数选择的自由度，造成某些阻抗参数过大^[26]。

不对称阻抗网络是基于参数位置不同则对具体脉动抑制能力不同的原则，对各个阻抗网络参数进行独立设计而成的；相比于对称阻抗网络，其对电容和电感的设计更加集约化合理化，从而可以缩小阻抗网络体积、增加功率密度和提高系统寿命^[26]，尤其适合于光伏模块集成型单相qZSI。针对采用不对称阻抗网络的单相ZSI/qZSI，也有众多学者开展了研究。文献[27-29]研究了阻抗网络不对称性带来的控制问题，但未涉及二倍频脉动分析的解析模型和阻抗参数的设计方法。文献[3]提出利用不同阻抗参数抑制不同脉动(高频/低频，电压/电流)的概念，并提出一种qZS网络不对称阻抗参数设计方法，但未给出二倍频脉动的解析模型，且系统并没有接入储能电池。文献[26]建立了单相BES-qZSI-I拓扑的二倍频脉动的解析模型，并提出一种不对称阻抗网络参数的设计方法；为了使BES-qZSI-I光伏系统能够24 h不间断向负载供电，文献[13]对该拓扑及其控制策略进行改进，并提出一种同时兼顾白天运行拓扑与黑夜运行拓扑的阻抗参数设计方法，但所述的解析模型和参数设计方法仅针对单相BES-qZSI-I拓扑。

综上所述，目前，针对单相BES-qZSI-II拓扑开展其二倍频脉动解析模型和不对称阻抗网络参数设计方法的研究较少，因此，本文作者针对不对称阻抗网络的单相BES-qZSI-II拓扑，首先分析其二倍频脉动的发生原因和传播机理，建立其二倍频脉动量的解析模型；其次分析其光伏板电压、储能电池电压和直流母线电压中二倍频脉动量对阻抗网络参数(C₁，C₂，L₁，L₂)及输入电容器电容(C_in)的敏感性，从抑制二倍频脉动和高频脉动的角度提出了一种单相BES-qZSI不对称阻抗网络参数的设计方法；最后，基于Matlab/Simulink的电路仿真实验验证了本文所建二倍频脉动解析模型的正确性。

1  单相BES-qZSI-Ⅱ系统DC侧动态模型

1.1　系统结构

基于BES-qZSI-II的单相光伏发电系统的结构如图1所示。qZS网络由2个电抗器(其电感分别为L₁和L₂)、2个电容器(其电容分别为C₁和C₂)以及双向有源开关S₅组成。光伏板连接到BES-qZSI-II的直流输入端口1(1a-1b)，负载/电网连接到BES-qZSI-II的H桥交流端口2(2a-2b)，储能电池连接到BES-qZSI-II的直流端口3(3a-3b)，从而构成一个三端口系统。双向有源开关S₅取代了传统qZS网络的二极管，其驱动信号与H桥直通信号互补，可避免轻载或低功率因数负载下可能出现的断续现象，同时实现交流端口和直流端口之间能量的双向流通^[30]。与文献[8]和[9]中将光伏板等效为恒压源不同的是，输入电容器可对光伏板电压进行滤波和稳定，使光伏板同时具备电压源和电流源的特性。

图1　基于BES-qZSI-II的单相光伏发电系统

Fig. 1　Single-phase PV system based on BES-qZSI-II

光伏板端电压v_in为

(1)

式中：V_pv为光伏板开路电压；i_pv为光伏板输出电流；R_s为光伏板的内阻。

储能电池端电压表达式为

(2)

式中：为储能电池的开路电压；i_b为储能电池输出电流；R_b为储能电池的内阻。

1.2　系统DC侧动态模型

为了建立单相BES-qZSI-II直流侧的动态模型，将H桥逆变电路等效为电流源i_PN，则BES-qZSI-II的等效电路如图2所示。

图2　单相BES-qZSI-II的等效电路

Fig. 2　Equivalent circuit of BES-qZSI-II

当单相BES-qZSI-II的逆变桥工作于直通状态时，直流母线电压v_PN为0，其等效电路如图2(a)所示，此时S₅关断，其反并联二极管承受反压而截止，逆变器实现了boost升压功能，对应动态方程为

(3)

式中：为微分算子；和分别为流过2个电抗器的电流；为电容器2的端电压。

当单相BES-qZSI-II的逆变桥工作在非直通状态时，其等效电路如图2(b)所示，此时S₅导通，对应动态方程为

(4)

式中：i_PN为直流母线电流。

利用状态变量平均法，由式(3)和式(4)可得BES-qZSI-II直流侧的动态模型为

(5)

式中：“< >”是开关周期平均算子，表示变量在1个开关周期内的平均值；D为直通占空比，在1个开关周期内恒定不变。

2  二倍频脉动的产生原因与传播机理

2.1　二倍频脉动的产生原因

2.1.1　交流输出功率的二倍频脉动　

假设单相BES-qZSI-II输出的交流基频电压v_a和基频电流i_a分别为

(6)

式中：V_a和I_a分别为基频电压v_a和基频电流i_a的幅值；ω为输出交流电的基频角频率；φ为功率因数角。

单相BES-qZSI-II输出的交流瞬时功率p_o为

(7)

式中：第1项为交流输出功率的平均值，第2项为以二倍基频变化的交流输出功率，它叠加在周期平均功率之上引起输出功率脉动，因此称之为交流输出功率的二倍频脉动分量，可用表示：

(8)

2.1.2　直流母线电流的二倍频脉动　

在图2(a)所示BES-qZSI-II的直通状态，直流母线电压v_PN=0，直流侧与交流侧之间没有功率传递；在图2(b)所示的非直通状态，直流母线电压v_PN达到峰值V_PN，功率由直流侧传递到交流侧。因此，在1个开关周期内，BES-qZSI-II的交流侧和直流侧之间满足功率平衡关系：

(9)

根据PWM调制原理，直流母线电压峰值V_PN与交流输出电压v_a之间满足

(10)

式中：为调制函数，M为PWM调制度。

将式(7)和式(10)代入式(9)，可得直流母线电流为

(11)

式(11)表明，直流母线电流i_PN包含直流分量I_PN和二倍频脉动分量：

(12)

(13)

直流母线峰值电压平均值V_PN、光伏板电压平均值V_in与直通占空比D之间存在如下关系：

(14)

联列式(8)、式(10)、式(13)和式(14)，求解得到

(15)

式(15)表明，单相BES-qZSI-II交流输出功率的二倍频脉动分量是导致直流母线电流出现二倍频脉动量的根本原因。

2.2　二倍频脉动的传播机理

由微分方程组(5)，可得BES-qZSI-II直流侧的状态空间方程为

(16)

式中：是对角矩阵；状态矢量；输入矢量。

从物理意义上讲，直流母线电流i_PN作为qZS网络的输入量，其直流分量I_PN和二倍频分量会分别在qZS网络的电压和电流中激励起各自的响应分量。因此，稳态时状态变量中必然包含有直流分量和二倍频分量，即状态变量可写为

(17)

式中：，，，分别为，，，的直流分量；带“~”的变量表示相应的二倍频分量。

3  二倍频脉动的解析建模

3.1　直流分量的解析表达式

稳态时，电容电压和电感电流的直流分量在1个开关周期内的平均变化率为0，因此，令式(5)中状态变量的微分项为0，即可得到

(18)

其中：I_b和I_pv分别为储能电池和光伏板的平均电流。

式(1)和式(2)中直流分量关系式为

(19)

联立式(18)和式(19)，可求得各状态变量的直流分量为

(20)

式(20)表明，qZS网络的稳态运行点不受阻抗参数的影响；即使考虑电抗器和电容器的等效串联电阻，其对稳态运行点的影响也很小，可以忽略。

3.2　二倍频脉动量的解析表达式

将各状态变量的直流分量从式(1)、式(2)和式(5)中分离，则稳态时各二倍频脉动量之间存在如下关系：

(21)

  (22)

令二倍频脉动量的表达式为

(23)

式中：“^”表示二倍频脉动量的幅值；φ_k(k=1，2，3，4，5)为二倍频脉动量的初相角。

联立式(21)~(23)，求得电容电压和电感电流的二倍频脉动分量的解析表达式

(24)

式中：

式(24)描述了直流母线电流的二倍频脉动分量与直流侧各电压和电流二倍频脉动量之间的关系；在一定的情况下，式(24)也反映了各二倍频脉动量与阻抗参数的关系。

4  二倍频脉动对阻抗参数的敏感性分析

根据式(21)，限制光伏板和储能电池的电压脉动量可限制其电流脉动量。因此，光伏板端电压v_in、储能电池端电压和直流母线电压v_PN中的二倍频脉动量必须被限制在工程容许范围之内。

4.1　二倍频电压脉动率的定义

为了定量地描述二倍频脉动量，引入二倍频电压脉动率，其定义分别为

(25)

式中：△v_PN，△v_in和分别为直流母线电压v_PN、光伏板端电压v_in和储能电池端电压的二倍频脉动率。

4.2　二倍频电压脉动率对阻抗参数的敏感性分析

假设单相BES-qZSI-II的稳态工作点为：V_pv=70 V，R_s=1 Ω，V_B=102 V，R_b=1 Ω，D=0.25，M=0.7，=33.3 V，阻感负载为1.2 mH+10 Ω，研究△v_PN，△v_in和△对qZS网络阻抗参数的敏感性。

4.2.1　△v_PN对阻抗参数的敏感性　

当L₁=1 mH，L₂=1 mH，C₂=3 000 μF时，直流母线电压的二倍频脉动率△v_PN随C₁和C_in的变化规律如图3(a)所示；当L₁=1 mH，L₂=1 mH，C_in=1 000 μF时，△v_PN随C₁和C₂的变化规律如图3(b)所示；当C₁=3 000 μF，L₂=1 mH，C_in=1 000 μF时，△v_PN随L₁和C₂的变化规律如图3(c)所示；当L₂=1 mH，C₂=1 000 μF，C_in=1 000 μF时，△v_PN随L₂和C₂的变化规律如图3(d)所示。

图3　二倍频脉动率△v_PN随阻抗参数的变化

Fig. 3　Relationship of △v_PN and impedance parameters

由图3(b)可见：随着电容C₂的增大，△v_PN显著减小；从图3(d)也可见：当C₂比较小时，△v_PN也会随着L₂的增加而缓慢减小；C₁和L₁对抑制△v_PN的作用不明显；C_in的增加对抑制△v_PN起负面作用，如图3(a)所示。

4.2.2　△v_in对阻抗参数的敏感性　

当L₂=1 mH，C₁=2 000 μF，C₂=3 000 μF时，光伏板端电压的二倍频脉动率△v_in随L₁和C_in变化规律如图4(a)所示；当L₂=2 mH，C₁=3 000 μF，C_in=1 000 μF时，△v_in随L₁和C₂变化规律如图4(b)所示；当C₂=3 000 μF，L₂=1 mH，C_in=1 000 μF时，△v_in随L₁和C₁变化规律如图4(c)所示；当L₁=1 mH，C₂=3 000 μF，C_in=1 000 μF时，△v_in随L₂和C₁变化规律如图4(d)所示。

图4　二倍频脉动率△v_in随阻抗参数的变化

Fig. 4　Relationship of △v_in and impedance parameters

由图4可知：C₁和L₁抑制二倍频脉动△v_in的能力几乎相等；C₂和L₂的增加对抑制△v_in起反作用；随着C_in的增加对抑制△v_in的变化是非单调的。

图5　二倍频脉动率△随阻抗参数的变化

Fig. 5　Relationship of △and impedance parameters

4.2.3　△对阻抗参数的敏感性　

当L₁=1 mH，L₂=1 mH，C_in=1 000 μF时，qZS网络电容器端电压的二倍频脉动率△随C₁和C₂变化规律如图5(a)所示；当L₂=1 mH，C₂= 3 000 μF，C_in=1 000 μF时，△随L₁和C₁变化规律如图5(b)所示；当C₂=3 000 μF，L₁=1 mH，L₂=1 mH时，△随C₁和C_in变化规律如图5(c)所示；当L₁=1 mH，C₂=3 000 μF，C_in=1 000 μF时，△随C₁和L₂变化规律如图5(d)所示。

由图5可知：L₁，L₂和C₂的增大会使二倍频脉动率△增大；抑制△只能依靠C₁和C_in的增加，而C₁的抑制能力比C_in的更强。

综上可知，二倍频电压脉动率△v_PN，△v_in和△与阻抗参数之间的关系如表1所示。

表1　二倍频电压脉动率与阻抗参数之间的关系

Table 1　Relationship of 2 ω ripple and impedance parameters

4.3　qZS网络阻抗参数的设计方法

单相BES-qZSI-II的阻抗网络中除二倍频电压电流脉动之外，还有高频电压电流脉动。通常抑制高频脉动所需的阻抗参数远小于抑制二倍频脉动的最小值，因此，在设计qZS网络参数时，首先应根据抑制高频脉动的需要确定阻抗参数的初始值，然后再对某些参数做等步长递增以解决二倍频脉动的抑制问题。文献[3]，[25]和[26]对高频脉动与阻抗参数的关系已经进行了深入研究，它们应该满足

(26)

(27)

式中：T_s为开关周期；k_sh为每个开关周期的直通次数；γ_i为电流的允许高频脉动率；γ_υ为电压的允许高频脉动率。

由表1可见，抑制某个二倍频电压脉动率，可以增大的阻抗参数不只1个，其中能起到显著抑制作用的参数被称之为主导参数；另一方面，某个主导参数的增大又可能会增加其他二倍频脉动率。因此，通过阻抗参数的合理设计来抑制单相qZS网络的二倍频电压脉动，不但要合理选择主导参数，而且还存在参数整定顺序的优化问题，以避免后续参数的整定对前已整定二倍频脉动的恶化。

基于表1，对阻抗网络参数抑制不同脉动的功能做如下约定：C_in和L₂主要用于抑制高频电压脉动和电流脉动，C₂主要用于抑制二倍频电压脉动率△v_PN，L₁主要用于抑制二倍频电压脉动率△v_in，C₁主要用于抑制二倍频脉动率△。基于上述功能约定，单相BES-qZSI阻抗网络的参数设计流程如图6所示。图6中，，和为期望的二倍频电压脉动率，和分别为电容和电感的单步增量。

图6　单相BES-qZSI-II阻抗网络的参数设计流程

Fig. 6　Parameters design procedure for single-phase BES-qZSI-II network

5  仿真验证

利用MATLAB/Simulink搭建了单相BES-qZSI-II光伏发电系统的电路模型，将电路仿真结果与解析计算结果进行比较，以验证解析模型的正确性。由于二倍频脉动量与阻抗参数的关系是主电路的自身特性，为了避免控制器性能对此特性的影响，仿真验证中采用开环且直通比恒定的控制方式。

文献[25]和[26]通过电路仿真结果验证解析模型的正确性，并进行了实验验证，研究表明：由于功率器件管压降和杂散损耗的存在，多个电压的实验测试值较解析计算值和电路仿真值都偏小，而基于理想器件的电路仿真结果则与解析计算结果几乎完全一致。因此，本文仅基于电路仿真结果的解析模型验证是合理而且充分的。

单相BES-qZSI-II仿真电路的参数为：M=0.7，D=0.25，V_pv=61 V，V_B=90 V，R_s=0.1 Ω，R_b=0.2 Ω，负载阻抗为10 Ω+1.2 mH，开关频率10 kHz。仿真验证选取的3个算例对应的不对称阻抗参数如表2所示。

表2　验证算例对应的阻抗参数

Table 2　Impedance parameters for validation examples

表3　算例1的二倍频脉动幅值仿真与计算结果比较

Table 3　Comparison of simulation results and calculation results of 2ω ripple for example 1

表4　算例2的二倍频脉动幅值仿真与计算结果比较

Table 4　Comparison of simulation results and calculation results of 2ω ripple for example 2

表5　算例3的二倍频脉动幅值仿真与计算结果比较

Table 5　Comparison of simulation results and calculation results of 2ω ripple for example 3

3个算例的二倍频脉动幅值仿真值与计算值的比较分别如表3~5所示。

由表3和表4可见：在其他参数保持不变的情况下，当C₁由4 400 μF减小为2 000 μF时，二倍频脉动率△v_in和的变化最显著，上升率分别高达24%和8.4%，这与表1反映的C₁对△v_in和起主导抑制作用的结论是一致的；虽然和的二倍频脉动量和均有增大，但是直流母线电压的二倍频脉动量反而减小了，这是不对称阻抗网络中和的相位不同造成的。从表4和表5可知：当其他参数保持不变，L₁由2 mH减小为1 mH时，二倍频脉动率△v_in的上升率高达101%，而反而减小了，这也与表1反映出的L₁与△v_in和的关系是一致的。

图7　算例1对应电路仿真波形与解析计算波形的比较

Fig. 7　Waveform comparison between analytic calculation and circuit-based simulation for example 1

算例1得到的电路仿真波形与解析计算波形的对比如图7所示，图7(a)和图7(b)所示分别为阻抗负载的工频50 Hz的电压和电流的波形，图7(c)~(k)所示分别为直流侧各主要电压和电流的波形。从图7可见：直流侧各电压和电流的电路仿真波形和解析计算波形都呈现出二倍频(100 Hz)脉动的特点，而且仿真得到的二倍频脉动量的幅值和相位与解析计算结果几乎完全相等；解析计算波形相对于电路仿真波形的最显著区别在于其光滑性，这是由于二倍频脉动解析模型是忽略了开关过程而基于状态空间平均法建立的，因此解析计算波形不能反映开关频率脉动特性。比较图7(f)和图7(h)可以发现，i_pv和i_L1虽然具有相同的直流平均分量和二倍频脉动分量，但前者的高频脉动比后者的要小得多，这是因为光伏板端电容C_in吸收了中部分高频分量；图7(c)和图7(d)表明，电容电压和的直流平均电压不相等，而且它们的二倍频脉动量也是不相等，这是qZS阻抗参数的不对称性造成的；从图7(g)和图7(h)可知：电抗器电流i_L1和i_L2的直流平均值不相等，它们的二倍频脉动量的幅值和相位也是不相等，这既与不对称阻抗参数有关，也与储能电池带来的电路结构不对称性有关；从图7(i)可知：电池电流i_b的平均值为负，说明此时电池工作在充电状态；图7(k)中v_PN的仿真波形反映了直流母线电压的脉冲状特性，而且其峰值电压的直流分量和二倍频脉动分量与图7(j)中的2个电容电压之和(+)的仿真和计算结果一致。

从表3~5可知，解析计算所得二倍频脉动幅值与电路仿真值的相对偏差率仅为1%左右；从图7可见：解析计算波形中的二倍频脉动分量的幅值和相位与电路仿真结果几乎完全相等。因此，仿真结果证明本文所建BES-qZSI-II的二倍频脉动解析模型是正确的。

需要说明的是，尽管在解析模型建立过程中忽略了电力电子器件的导通压降，但是二倍频脉动量的实际值与解析计算值之间也几乎没有偏差，这是因为电力电子器件管压降影响的是直流平均电压，而二倍频脉动量叠加在直流平均电压上，因此，二倍频脉动量不受电力电子器件管压降的影响。

6  结论

1) 单相输出功率的二倍频功率分量是激励单相BES-qZSI-II产生二倍频电压脉动和电流脉动的根源。

2) 基于各阻抗元件的作用不同而采用不对称阻抗参数的qZS网络有利于实现其体积和成本的最优化。

3) 从抑制电压和电流脉动的角度设计qZS网络的阻抗参数，不但要正确选择抑制各脉动量的主导参数，而且要合理安排参数整定顺序。

4) 储能电池的并入和阻抗参数的不对称使得qZS网络的2个电容器直流电压、2个电抗器直流电流的平均值都不相等，而且它们的二倍频脉动量也不再对称。

5) 本文所建解析模型适用于对BES-qZSI-II的二倍频脉动量做出快速准确的评估。

参考文献：

[1] 李媛, 彭方正. Z源/准Z源逆变器在光伏并网系统中的电容电压恒压控制策略[J]. 电工技术学报, 2011, 26(5): 62-69.

LI Yuan, PENG Fangzheng. Constant capacitor voltage control strategy for Z-source/quasi-Z-source inverter in grid-connected photovoltaic systems[J]. Transactions of China Electrotechnical Society, 2011, 26(5): 62-69.

[2] ABU-RUB H, IQBAL A, MOIN AHMED S, et al. Quasi-Z-source inverter-based photovoltaic generation system with maximum power tracking control using ANFIS[J]. IEEE Transactions on Sustainable Energy, 2013, 4(1): 11-20.

[3] ZHOU Yan, LIU Liming, LI Hui. A high-performance photovoltaic module-integrated converter(MIC) based on cascaded quasi-Z-source inverters(qZSI) using eGaN FETs[J]. IEEE Transactions on Power Electronics, 2013, 28(6): 2727-2738.

[4] LIU Yushan, GE Baoming, ABU-RUB H, et al. Modelling and controller design of quasi-Z-source inverter with battery-based photovoltaic power system[J]. IET Power Electronics, 2014, 7(7): 1665-1674.

[5] 李佳琪, 陈健, 张文, 等. 高渗透率光伏配电网中电池储能系统综合运行控制策略[J]. 电工技术学报, 2019, 34(2): 437-446.

LI Jiaqi, CHEN Jian, ZHANG Wen, et al. Integrated control strategy for battery energy storage systems in distribution networks with high photovoltaic penetration[J]. Transactions of China Electrotechnical Society, 2019, 34(2): 437-446.

[6] LI Xiangjun, HUI Dong, LAI Xiaokang. Battery energy storage station(BESS)-based smoothing control of photovoltaic(PV) and wind power generation fluctuations[J]. IEEE Transactions on Sustainable Energy, 2013, 4(2): 464-473.

[7] CINTRON-RIVERA J G, LI Yuan, JIANG Shuai, et al. Quasi-Z-Source inverter with energy storage for Photovoltaic power generation systems[C]// 2011 Twenty-Sixth Annual IEEE Applied Power Electronics Conference and Exposition(APEC). Fort Worth, TX, USA: IEEE, 2011: 401-406.

[8] GE Baoming, ABU-RUB H, PENG Fangzheng, et al. An energy-stored quasi-Z-source inverter for application to photovoltaic power system[J]. IEEE Transactions on Industrial Electronics, 2013, 60(10): 4468-4481.

[9] GE Baoming, PENG Fangzheng, ABU-RUB H, et al. Novel energy stored single-stage photovoltaic power system with constant DC-link peak voltage[J]. IEEE Transactions on Sustainable Energy, 2014, 5(1): 28-36.

[10] RASIN Z, AHSANULLAH K, RAHMAN M F. Design and simulation of quasi-Z-source grid-connected PV inverter with bidirectional power flow for battery storage management[C]// 39th Annual Conference of the IEEE Industrial Electronics Society. Vienna, Austria: IEEE, 2013: 1589-1594.

[11] ZHANG Haiyu, LI Xiao, GE Baoming, et al. Capacitance, dc voltage utilization, and current stress: comparison of double-line frequency ripple power decoupling for single-phase systems[J]. IEEE Industrial Electronics Magazine, 2017, 11(3): 37-49.

[12] LIANG Zipeng, HU Sideng, HE Xiangning. Analysis and suppression strategy for the double-line frequency pulsation in single-phase quasi-Z-source converter[J]. IEEE Transactions on Power Electronics, 2019, 34(12): 12567-12576.

[13] LIANG W, LIU Y, GE B, et al. Night operation, analysis, and control of single-phase quasi-Z-source photovoltaic power system[J]．IET Renewable Power Generation, 2019,13(15): 2817-2829．

[14] 董帅, 张千帆, 王睿, 等. SVPWM控制时双向Z源逆变器电容电压纹波分析[J]. 电工技术学报, 2017, 32(24): 107-114.

DONG Shuai, ZHANG Qianfan, WANG Rui, et al. Analysis of capacitor voltage ripple for Bi-directional Z-source inverters based on SVPWM[J]. Transactions of China Electrotechnical Society, 2017, 32(24): 107-114.

[15] 高志强, 纪延超, 孙勇, 等. Z源直流母线电压波动的抑 制[J]. 哈尔滨理工大学学报, 2011, 16(4): 86-89.

GAO Zhiqiang, JI Yanchao, SUN Yong, et al. Suppression of voltage fluctuation on DC link voltage of Z-source[J]. Journal of Harbin University of Science and Technology, 2011, 16(4): 86-89.

[16] GE Baoming, LIU Yushan, ABU-RUB H, et al. An active filter method to eliminate DC-side low-frequency power for a single-phase quasi-Z-source inverter[J]. IEEE Transactions on Industrial Electronics, 2016, 63(8): 4838-4848.

[17] 李佳逸, 雷勇, 李媛, 等. 单相准Z源逆变器二倍频纹波的抑制及应力分析[J]. 电力系统保护与控制, 2018, 46(9): 25-33.

LI Jiayi, LEI Yong, LI Yuan, et al. Suppression and stress analysis of double-line frequency ripple of single-phase quasi Z source inverter[J]. Power System Protection and Control, 2018, 46(9): 25-33.

[18] GE Baoming, LIU Yushan, ABU-RUB H, et al. Current ripple damping control to minimize impedance network for single-phase quasi-Z source inverter system[J]. IEEE Transactions on Industrial Informatics, 2016, 12(3): 1043-1054.

[19] 李志勇, 王思君, 陈有根. Quasi-Z源逆变器纹波自注入APF及其控制策略[J]. 电机与控制学报, 2019, 23(10): 77-84.

LI Zhiyong, WANG Sijun, CHEN Yougen. A ripple self-injection active power filter control scheme for Quasi-Z-source inverter[J]. Electric Machines and Control, 2019, 23(10): 77-84.

[20] YU Yifan, ZHANG Qianfan, LIANG Bin, et al. Single-phase Z-source inverter: analysis and low-frequency harmonics elimination pulse width modulation[C]// 2011 IEEE Energy Conversion Congress and Exposition. Phoenix, AZ, USA: IEEE, 2011: 2260-2267.

[21] LIU Yushan, GE Baoming, ABU-RUB H, et al. Hybrid pulsewidth modulated single-phase quasi-Z-source grid-Tie photovoltaic power system[J]. IEEE Transactions on Industrial Informatics, 2016, 12(2): 621-632.

[22] 黄守道, 张阳, 罗德荣. Z源逆变器在风电并网系统中的电容电压纹波抑制策略[J]. 电工技术学报, 2015, 30(S2): 135-142.

HUANG Shoudao, ZHANG Yang, LUO Derong. Ripple suppression for the capacitor voltage of the Z-source grid-connected wind generation system[J]. Transactions of China Electrotechnical Society, 2015, 30(S2): 135-142.

[23] YU Yifan, ZHANG Qianfan, LIU Xiaofei, et al. DC-link voltage ripple analysis and impedance network design of single-phase Z-source inverter[C]// Proceedings of the 14th European Conference on Power Electronics and Applications. Birmingham, UK: IEEE, 2011: 1-10.

[24] SUN Dongsen, GE Baoming, YAN Xingyu, et al. Modeling, impedance design, and efficiency analysis of quasi-Z-source module in cascaded multilevel photovoltaic power system[J]. IEEE Transactions on Industrial Electronics, 2014, 61(11): 6108-6117.

[25] LIU Yushan, GE Baoming, ABU-RUB H, et al. Comprehensive modeling of single-phase quasi-Z-source photovoltaic inverter to investigate low-frequency voltage and current ripple[J]. IEEE Transactions on Industrial Electronics, 2015, 62(7): 4194-4202.

[26] LIANG Weihua, LIU Yushan, GE Baoming, et al. Double-line-frequency ripple model, analysis, and impedance design for energy-stored single-phase quasi-Z-source photovoltaic system[J]. IEEE Transactions on Industrial Electronics, 2018, 65(4): 3198-3209.

[27] THELUKUNTLA C S, VEERACHARY M. Analysis of single phase asymmetric Z-source inverter[C]// 2015 International Conference on Energy, Power and Environment: Towards Sustainable Growth(ICEPE). Shillong, India: IEEE, 2015: 1-6.

[28] HU Sideng, LIANG Zipeng, HE Xiangning. Ultracapacitor-battery hybrid energy storage system based on the asymmetric bidirectional Z-source topology for EV[J]. IEEE Transactions on Power Electronics, 2016, 31(11): 7489-7498.

[29] KAYIRANGA T, LI Hongbo, LIN Xinchun, et al. Abnormal operation state analysis and control of asymmetric impedance network-based quasi-Z-source PV inverter(AIN-qZSI)[J]. IEEE Transactions on Power Electronics, 2016, 31(11): 7642-7650.

[30] GUO Feng, FU Lixing, LIN C H, et al. Development of an 85-kW bidirectional quasi-Z-source inverter with DC-link feed-forward compensation for electric vehicle applications[J]. IEEE Transactions on Power Electronics, 2013, 28(12): 5477-5488.

(编辑  赵俊)

收稿日期： 2020 -05 -28; 修回日期： 2020 -08 -24

基金项目(Foundation item)：国家自然科学基金资助项目(51867016) (Project(51867016) supported by the National Natural Science Foundation of China)

通信作者：魏祥林，博士研究生，讲师，从事新能源变换与控制、电机及其控制等研究；E-mail：wxl2004@lut.edu.cn

引用格式： 魏祥林, 王晗, 陈伟, 等. 单相BES-quasi-ZSI二倍频脉动的解析建模[J]. 中南大学学报(自然科学版), 2021, 52(5): 1502-1515.

Citation: WEI Xianglin, WANG Han, CHEN Wei, et al. Analytic modeling of double-line frequency ripple in single-phase BES-quasi-ZSI[J]. Journal of Central South University(Science and Technology), 2021, 52(5): 1502-1515.