LCL逆变器型直驱风电机的稳定性

2021-09-26苏勋文安鹏宇王浠再卢宸赐

黑龙江科技大学学报 2021年5期

苏勋文，安鹏宇，王浠再，卢宸赐

(黑龙江科技大学电气与控制工程学院，哈尔滨 150022)

0 引言

随着风能的利用率不断提高，风力发电以其独特的优势得到迅速发展。伴随电网中新能源接入数量的增加，架设更多更长的输电线路会向电网中引入大量线路阻抗，使电网逐渐呈现出弱网特性[1-3]，风电场阻抗与电网阻抗的相互影响导致系统失稳，为解决这一问题，高本锋[4]、宋瑞华[5]和苏田宇[6]提出调整控制器参数法、附加阻尼法与配置机网侧阻尼控制等多种方法改善风电并网后系统的运行稳定性。Cespedes等[7]指出可通过更改附加的滤波装置实现系统的稳定运行，但采用单电感滤波器需要较大电感值的电感，其滤波器成本高、体积大且动态性能差，已经不符合大规模风电并网的需求[8-9]，通过LCL型滤波器能够选取较小的电感实现同样的滤波效果[10]。Wen等[11]给出了阻抗法的单电感型滤波器的并网逆变器阻抗模型，其滤波电感值较大。夏丹妮[12]指出了通过阻抗法的基于LC型并网逆变器的直驱风电机阻抗模型，其小信号模型未考虑电容电流支路。

为了进一步提高并网稳定性，减小直驱风电机并网成本，笔者以阻抗法的LCL逆变器型直驱风电机为研究对象，综合锁相环、电流内环与功率外环的特性，在两相旋转坐标系下推导其等效阻抗模型，通过广义奈奎斯特判据分析部分参数对稳定性的影响。

1 直驱风电机网侧小信号模型

1.1 直驱风电机基本原理

直驱风力发电机原理如图1所示。

图1 直驱风力发电机原理Fig. 1 Schematic diagram of PMSG

直驱风电机由风轮、发电机、变流器及其控制系统组成。在整个直驱风力发电机运行时，可以忽略风轮、轴系产生的机械振荡，将风轮、发电机和发电机侧变流器通过戴维宁定理等效为一个直流电压源，通过调整电压值的大小实现控制功率的大小[8,12]。

1.2 dq坐标系下的开环小信号模型

直驱风电机网侧逆变器的主电路结构如图2所示，其中，Udc为直流电压稳态值，Cdc为直流电压母线稳压电容，S1～S6为反并联自关断二极管，Lf1、Lf2和Cf分别为LCL滤波器的滤波电感与滤波电容，Rf1、Rf2和RC分别为LCL滤波器上滤波电感与滤波电容上的寄生电阻，Lg为网侧电感，Rg为网侧电感上的寄生电阻，ea、eb和ec分别为网侧三相电压。

图2 网侧逆变器的主电路拓扑Fig. 2 Topology of three-phase LCL grid-connected inverter

为了从并网点处得到整个系统的阻抗模型，需要先得到其传递函数，根据图2所示结构进行派克变换，得到在两相旋转坐标系(dq坐标系)下的小信号模型，如图3所示。

图3 qd两相旋转坐标系下的小信号模型Fig. 3 Small-signal model of LCL inverter in qd frame

图3中，Dsd、Dsq、dsd和dsq分别表示在两相旋转坐标系下主电路电压的稳态值与扰动值，udc为直流电压扰动值，isd、isdc和usd，isq、isqc和usq分别表示两相旋转坐标系下d轴与q轴并网点电流扰动值、电容支路电流扰动值和并网点电压扰动值。

由图3可以看到

对于开环输出阻抗Zout，可以令dsd、dsq和udc为0，则有

(1)

Lf=Lf1+Lf2，

Rf=Rf1+Rf2，

式中：GRL——LCL滤波器阻抗传递函数矩阵；

GR1——LCL滤波器逆变电路侧阻抗传递函数矩阵；

Lf——LCL滤波器等效电感；

Rf——LCL滤波器等效电阻。

对于开环主电路占空比到并网点电流的传递函数Gid，可以令usd、usq和udc为0，则

(2)

式中，Gdc——dq坐标系下直流电压稳态值的传递函数矩阵。

根据基尔霍夫电压定律，并网点处电压扰动关系为

(3)

式中：GR2——LCL滤波器网侧阻抗传递函数矩阵；

GC——LCL滤波器电容支路阻抗传递函数矩阵。

将式(3)分别代入式(1)和(2)，可以得到并网点电压到并网点电流的传递函数为

占空比到并网点电流的传递函数为

此时，开环情况下的主电路传递函数如图4所示。其中，us为主电路中的电压扰动值，is为主电路电流扰动值，ds为主电路占空比扰动值。

图4 开环主电路传递函数Fig. 4 Inverter power stage model with open loop

1.3 dq坐标系下的闭环模型

在并网逆变器的控制策略中，通过锁相环(PLL)获取控制系统与电网系统的相位信息，使并网逆变器与电网同步。而锁相环内部使用PI控制器进行控制，原理如图5所示。图中，ua、ub和uc分别为锁相环采样点处的三相电压，ud和uq分别为经过派克变换后的d轴和q轴电压，θPLL为主电路系统与控制系统间在扰动存在时的相位差。

图5 锁相环原理Fig. 5 Principle of PLL

主电路系统与控制系统间在扰动存在时的相位差θPLL可以定义为

(4)

式中：KpP——锁相环PI调节器比例系数；

KiP——锁相环PI调节器积分系数。

当电网侧出现扰动时，电网与控制系统就会出现相位差，两坐标系之间的相位关系矩阵为

(5)

由于扰动量导致的相位变化很小，因此，可以将cosθPLL≈1，sinθPLL≈θPLL，则式(5)为

主电路与控制电路的关系为

(6)

式中：ucd——控制电路中d轴扰动值；

ucq——控制电路中q轴扰动值；

Usd——主电路中d轴电压稳态值；

Usq——主电路中q轴电压稳态值。

联立式(4)和(6)，可得

θPLL=GPLLusq,

式(6)可以变为

因此，可以得到主电路电压、电流和占空比到控制电路电压、电流和占空比的传递函数矩阵

式中：Isd——主电路中d轴电流稳态值；

Isq——主电路中q轴电流稳态值；

Dsd——主电路中d轴占空比稳态值；

Dsq——主电路中q轴占空比稳态值；

GPLL——锁相环传递函数矩阵。

含锁相环的主电路控制模型如图6所示。图中，ic为控制电路电流扰动值，dc为控制电路占空比扰动值。

图6 含锁相环的主电路传递函数Fig. 6 Main system tramsfer function with PLL

此时的输出阻抗为

除了通过锁相环对直驱风电机系统的并网进行控制，还通过电流内环和功率外环两个环节实现系统的控制。在电流内环中，通过一个PI调节器实现对电流的控制，电流内环控制矩阵为

式中：KPi——电流内环PI调节器比例系数；

KIi——电流内环PI调节器积分系数。

在系统通过派克变换从三相静止坐标系到两相旋转坐标系的变换后，电流内环还会存在正负序电流的耦合，因此，还需要加入解耦矩阵：

含锁相环与电流内环的主电路控制模型如图7所示。图中，iref为电流参考值。由于该模型为小信号模型，为了尽可能减小使扰动值，故将电流参考值设为0。

图7 含锁相环与电流环的主电路传递函数Fig. 7 Main system oransfer function with PLL and current control

此时，输出阻抗为

式中，E——二阶单位矩阵。

为了实现对电流内环中的参考电流进行更加准确的控制，现选择在电流内环的控制器之前加入功率外环，确保逆变器产生正确的无功功率和有功功率。在控制电路的两相旋转坐标系中，有功功率与无功功率的表达式为

(7)

式中：P——有功功率；

Q——无功功率；

icd——控制电路中d轴电流稳态值；

icq——控制电路中q轴电流稳态值；

ucd——控制电路中d轴电压稳态值；

ucq——控制电路中q轴电压稳态值。

当系统稳定时，通过对式(7)进行线性化可以得到：

式中：GiPQ——电流到功率的传递函数矩阵；

GvPQ——电压到功率的传递函数矩阵；

Isd——主电路中d轴电流稳态值；

Isq——主电路中q轴电流稳态值；

Usp——主电路中d轴电压稳态值；

Usq——主电路中q轴电压稳态值。

功率外环的PI控制器传递函数矩阵为

式中：KpPQ——功率外环PI调节器比例系数；

KiPQ——功率外环PI调节器积分系数。

含锁相环、电流内环与功率外环的主电路控制模型如图8所示，Pref、Qref为有功功率与无功功率参考值，与电流参考值选取目的一致，故设为0。

图8 含锁相环电流环和功率环的主电路传递函数Fig. 8 Small-signal model with PLL, current and power control

此时的输出阻抗为

2 直驱风电机的稳定性

在传统的直驱风电机并网中，采用单个电感进行滤波。为了更加优异的并网效果，需要不断增加电感值的大小，随着电感大小的增加，电网强度便会随之下降，导致系统的动态性能同步降低。将通过对比单电感型滤波器和LCL型滤波器在直驱风电机并网中的实际效果。

2.1 LCL型滤波器谐振

根据图2所示的直驱风电机并网侧电路结构，其基尔霍夫电压方程为

(8)

式中：if1——滤波器逆变器侧电感电流值；

if2——滤波器网侧电感电流值；

uc——滤波器电容支路电压值；

uinv——逆变器输出电压值；

uR1——滤波器逆变器侧电感寄生电阻值；

uR2——滤波器网侧电感寄生电阻值；

eg——并网侧电压值。

为了简化计算，忽略滤波器寄生电阻，对式(8)进行拉普拉斯变换可得,逆变器输出电压到并网电流的传递函数为

(9)

由式(9)可知，传递函数HLCL的分母不含二次项，此时，系统易受到尖峰谐振的干扰，现通过在电容支路添加一个电阻RC，获得新的传递函数矩阵为

同理，可以得到单电感滤波器的逆变器输出电压到并网电流传递函数为

在Matlab/Simulink中对上述三种情况下的幅频特性分析结果如图9所示，单滤波电感滤波器的电感值为3 mH，LCL型滤波器电感值为0.3 mH。含电容支路电阻的LCL型滤波器电感值为0.3 mH。

从图9可以看出，三种滤波器具有相近的滤波效果，对于电容支路不含电阻的LCL型滤波器，在某一特定频率下会出现尖峰，造成系统失稳，在引入电容支路电阻后，尖峰得到抑制，保证系统能够在高频段下稳定运行。

图9 滤波器幅频特性波德图Fig. 9 Bode of filters amplitude frequency characteristics

2.2 阻抗模型

单电感型滤波器与LCL型滤波器进行对比如图10所示。

图10 直驱风电机输出阻抗模型的波德图Fig. 10 Bode of impedences with different filters

由于dq、qd分量的幅值与相位趋势与dd、qq轴相近，故仅对dd与qq轴，通过其波德图分析稳定性。由于在相同控制模式与控制参数下，单电感滤波器需要更大值的电感才能保证系统稳定，将LCL型滤波器的电感值调整至单电感型滤波器相同值，通过Matlab/Simulink搭建不同滤波器下的直驱风电机模型，由图10可以看到，在1到100 Hz内，无论d轴q轴的曲线都有较高程度的拟合；在100到1000 Hz内，dq分量和qq分量的负序范围都有明显增大，其中，电容支路不串接电阻的LCL型滤波器变化尤为明显，这种情况下系统极易失稳，电容支路串接电阻的LCL滤波器对整个系统的稳定效果更优，可以保证直驱风电机系统并网的稳定性需求。

3 仿真与结果分析

由于在风电并网中，电网显感性，呈现弱电网状态，控制器参数会对并网稳定性产生不可忽略的影响。为了探究不同控制器参数下的LCL型并网逆变器的直驱风电机系统的稳定性，通过广义奈奎斯特判据(GNC)进行论证。需要建立直驱风电机输出阻抗模型与网侧阻抗模型的回比矩阵N(s)，根据广义奈奎斯特判据对系统稳定性的判别依据，当回比矩阵N(s)的GNC曲线围绕(-1，j0)点时，系统失稳，出现振荡的风险加大；反之，系统稳定[12-13]。关系式为

(10)

式中：Zg(s)——电网等效阻抗矩阵；

ZPMSG(s)——直驱风电机等效阻抗矩阵；

Rg——电网等效电阻；

Lg——电网等效电感。

主电路仿真参数：PN为2 MW，Uabc为690 V，Lf1为0.8 mH，Lf2为0.8 mH，Cf为20 μF，Cdc为4 700 μF，Udc为1 200 V和Lg为0.6 mH。控制电路参数：KpP为40,KiP为180,KpPQ为0.5,KiPQ为2,KPi为0.6,KIi为35。

3.1 电网电感对直驱风电机稳定性的影响

由于电网阻抗与主电路的各项参数会产生耦合现象，电网等效电感值越大，弱网特性越明显[12,14]，对系统的稳定运行的状态影响越大。稳定运行下的回比矩阵广义奈奎斯特图如图11所示。电网电感分别为0.5、1.0与1.5 mH时的奈奎斯特图见图12。

图11 稳定时的直驱风电机奈奎斯特图Fig. 11 GNC plots in stable case

图12 不同电网电感的直驱风电机奈奎斯特图Fig. 12 GNC plots with different grid inductance

由图11、12可见，此时曲线没有包围稳定判别点(-1,j0)，文中均将采用局部放大图对稳定性进行分析。在直驱风电机系统并网时，线路上的寄生电阻较小，电网强弱主要受线路上的电感影响。通过图12可见，当电网强度变弱，即电网电感逐渐增大，由0.5 mH增大至1.5 mH时，红色实线逐渐包围稳定判别点。根据广义奈奎斯特稳定性判据，系统从稳定状态过渡到了不稳定状态，如果继续增加电网电感，系统将出现不可抑制的振荡。

3.2 控制器参数对直驱风电机稳定性的影响

通过将讨论锁相环比例系数与电流环比例系数对系统的影响。电网电感为0.5 mH，锁相环比例系数分别为40、80、120时的奈奎斯特曲线如图13所示。

图13 不同锁相环比例系数的直驱风电机奈奎斯特图Fig. 13 GNC plots with different proportionality coefficient of PLL

与不同滤波器的波德图反映出的结果近似，由于在两相旋转坐标系(dq坐标系)中推导的阻抗模型已计及解耦，因此dq与qd两个通道的耦合分量在不同锁相环比例系数下的波德图显示出的变化趋势与dd与qq分量近似，故在这里仅展示dd分量的幅值与相位和qq分量的幅值与相位。由图13可见，锁相环比例系数从40到120变化时对系统的稳定性影响不大。

为了保证结论的正确性，结合不同比例系数下的波德图进行判定，如图14所示。