LCL型并网逆变器电流控制稳定性分析

成林俞++刘红++文丹丹

摘 要： 在LCL型光伏并网逆变器中，电流控制器的比例控制系数设计不合理易造成系统不稳定，甚至损坏逆变器。对采用逆变侧电感电流反馈的LCL型三相光伏并网逆变器进行数学建模，运用Routh?Hurwitz稳定性判据基于系统离散化模型进行稳定性分析，可求得比例控制系数的精确范围。分析逆变侧电感、滤波电容以及网侧电感的取值与比例控制系数的关系，探讨滤波电感和电容的取值对系统稳定性的影响。通过Matlab仿真及在50 kW光伏逆变器实验平台上进行测试，验证了所进行的稳定性分析的正确性。

关键词： 逆变器； 电流控制器； 稳定性判据； 比例控制

中图分类号： TN710?34； TM464 文献标识码： A 文章编号： 1004?373X（2015）16?0141?04

Stability analysis on current control of grid?connected inverter with LCL filter

CHENG Linyu1， LIU Hong2， WEN Dandan1

（1. Zhejiang Province Engineering Laboratory of Electrical Digital Design Technology， Wenzhou University， Wenzhou 325035， China；

2. Zhejiang Eifeisun Energy Technology Co.， Ltd.， Wenzhou 325032， China）

Abstract： The unreasonable design of proportional control coefficient of current controller is easy to lead to system instability or even damage of the photovoltaic grid?connected inverter with LCL filter. A mathematical model of three?phase photovoltaic grid inverter with LCL filter was established， which adopts inductive current feedback in inverter side. and then applies the Routh?Hurwitz stability criterion is used to analyze the stability of the inverter system on the basis of system discretization model， with which the precise range of proportional control coefficient can be derived. The relationship between values of the inverter side inductance， filtering capacitor and network side inductance， and proportional control coefficient is analyzed. The influence of filtering inductance and capacitance values on the system stability is discussed. The correctness of the stability analysis was demonstrated by the results of simulation with Matlab and test on a test?bed for 50 KW photovoltaic grid?connected inverter.

Keywords： inverter； current controller； stability criterion； proportional control

0 引 言

人类社会的快速发展使能源消耗量剧增，石油等不可再生能源的日益枯竭引起全球各界的担忧。石油、煤等化石能源燃烧排放出大量的二氧化碳、二氧化硫等气体引起温室效应、酸雨等环境问题，对人类的生存和发展构成严重的威胁。发展新能源已经成为全球能源变革的大趋势。太阳能是人类可利用的最丰富的清洁能源而得到世界各个国家的重视，我国先后实施“金太阳示范工程”和颁布一系列政策支持分布式光伏发电，大力推动了光伏产业的发展。

逆变器作为光伏并网系统的核心，其控制技术是整个并网系统的重点。电流谐波含量是逆变器并网技术规范中的一项重要指标，在并网逆变器中一般采用适当的输出滤波器抑制并网电流谐波含量。L滤波器的体积大且滤波效果不理想，LCL滤波器减小了体积且对高频电流谐波分量有更强的抑制能力而得到广泛的应用[1?4]。由于LCL滤波器是三阶系统，二阶谐振零极点多，若控制器设计不当容易造成系统的不稳定[5?8]。

在工程应用领域，并网逆变器常用的控制方法有比例积分控制（PI）、比例谐振控制（PR）等[9?10]。比例控制对逆变器控制系统的稳定起主导作用，积分控制和谐振控制是为消除静差而引入。为此，本文运用一种适合工程应用的离散化方法，对采用逆变侧电流反馈的LCL型三相并网逆变器的数学模型进行离散化分析，根据Routh?Hurwitz稳定性判据可以得到比例控制系数的精确范围。在此基础上，探讨电感和电容的取值变化与比例控制系数的关系，由此探讨电感和电容对系统稳定性的影响。endprint

1 系统建模及离散化

1.1 光伏逆变器模型

图1为三相LCL型并网逆变系统原理图。直流电压Udc由太阳能电池板产生，经全桥逆变电路后将直流电转换成交流电，通过控制逆变侧电感电流的频率和相位来跟踪电网电压的频率和相位，使输出电流保持为与电网电压同频同相的正弦信号，从而达到并网运行的目的。图中L1为逆变侧电感，C为滤波电容，L2为网侧电感。

图1 三相LCL型光伏并网逆变器原理图

1.2 系统建模及离散化

三相逆变器的调制电压[uM]设计为电网电压前馈[ug]加上电流给定[i*]与电流反馈[i]的误差乘以比例控制系数[kp]的组合，即：

[uM=ug+kP?i*-i] （1）

其中，电网电压前馈[ug]抵消电网电压的变化，因此比例控制的控制方程式为：

[u=kP?i*-i] （2）

为便于分析，假设三相电网平衡，电感、电容及开关管均为理想器件，则三相逆变器的每一相的工作状态一样。建立单相LCL滤波器的模型，可分析三相LCL滤波器。单相控制系统数学模型如图2所示。

图2 控制系统模型

这里取逆变侧电感电流[i1]，电容电流[iC]，及电容两端的电压[uC]为状态变量，得到下列方程：

[ddti1iCuC=00-1L100-1L1+1L201C0?i1iCuC+1L1uu0]

由于实际工程控制中一般是离散数字控制系统，因此需将上式离散化，即已知当前状态取值的条件下，求解出下一时刻状态的取值。这里利用Laplace变换求解，状态方程的离散化形式为：

[i1n+1iCn+1uCn+1=B1i1niCnuCn+B2?un]

采用比例控制，反馈电流从L1处采样，由式（2）可知：

[un=kP?i*n-i1n] （3）

可知加入比例控制后的离散状态方程为：

[i1n+1iCn+1uCn+1=Mi1niCnuCn+B2?i*n]

其中转移矩阵为：

于是，[|λ|<1]等价于[x]的实部小于零。

根据Routh?Hurwitz稳定性判据条件，求得比例控制系数的范围为：

[0

式中：[ωr=L1+L2CL1L2]。

2.2 仿真验证

仿真参数：逆变侧电感L1=160 μH，滤波电容C=450 μF，网侧电感L2=40 μH，采样周期T=100 μs。根据上述理论推导可知，求得比例控制系数的范围为：03.16取值kP=3.2时，逆变器系统谐振，电流波形发散，不能稳定控制。由此可知，该离散化方法能够精确计算出比例控制系统的稳定控制范围。

图3 kP=1时逆变器输出电流波形

2.3 实验验证

为进一步验证所进行的稳定性分析的正确性和有效性，在50 kW光伏逆变器平台上进行实验，根据该平台的参数计算出相应的比例控制系数。图6所示为kP取不同值时，逆变器的A相电流输出波形。由图可知，当kP=3.1时，系统稳定；当kP=2.05时，系统接近谐振，输出波形较差；当kP=2.1时，系统谐振。

图4 kP=3.16时逆变器输出电流波形

图5 kP=3.2时逆变器输出电流波形

图6 kP取不同值时逆变器输出电流波形

3 比例控制系数kP的影响因素

3.1 电容C的取值对kP的影响

设函数[fC=2L1+L2T+2L2ωrL1?tanωrT2]，由[0<ωrT<π]可知电容的取值范围：

[C>T2?L1+L2π2?L1?L2=31.7e-6]

本仿真的参数为：L1=（160e-6） H，L2=（40e-6） H，

T=（100e-6） s。电容C的取值范围为（32e-6）～（1e-3） F时，函数[fC]的曲线如图7所示。

图7 电容C的取值对比例控制系数的影响

分析以上仿真结果可知，当电容C的取值大于一定值的时候，此时[tanωrT2≈ωrT2]，电容C对[kP]的取值无影响，[kP]趋于稳定：

[fC=2L1+L2T+2L2ωrL1?tanωrT2≈2L1+L2T+2L2ωrL1?ωrT2=2L1+L2T+L2?TL1]

3.2 电感L1的取值对kP的影响

本仿真的参数为：C=（450e-6） F，L2=（40e-6） H，T=（100e-6） s：

设函数[fL1=2L1+L2T+2L2ωrL1?tanωrT2]。