基于比例谐振和谐波补偿控制技术的单相逆变并网研究

张　淼，苏协飞

(广东工业大学 自动化学院, 广东 广州 510006)



基于比例谐振和谐波补偿控制技术的单相逆变并网研究

张淼，苏协飞

(广东工业大学 自动化学院, 广东 广州 510006)

电网电压背景谐波的干扰在并网系统中是不可避免的，这会给并网电流带来大量的低阶谐波.提高开关频率能够抑制低阶谐波，却带来大量的高阶谐波.针对基于比例谐振控制(PR)和谐波补偿(HC)控制器的LCL型单相并网逆变器，详细分析了电流控制器参数对系统性能的影响，通过对并网电流的稳态误差、相位裕度和幅值裕度的分析得出满足要求的控制器各参数选择范围，结合系统的稳定性分析来选择控制器参数.仿真结果证明所给出的控制策略和参数设计方法是可行的.

并网逆变器;比例谐振控制;谐波补偿;LCL滤波器

受到环保意识增强和技术发展迅速等因素的影响，分布式能源的应用变得越来越广泛，光伏逆变并网技术也愈来愈被重视[1].在并网逆变器中，并网处常用的滤波器一般为L型滤波器或LCL滤波器.L型滤波器结构简单，但体积较大并且对高频谐波的抑制能力有限；LCL滤波器具有三阶低通滤波特性，不仅成本低而且对高频谐波有很好的衰减效果[2-3]，在逆变并网中被广泛使用.但LCL滤波器的频率响应存在明显的谐振，文献[4-5]分析了LCL滤波器各种谐振阻尼方案的运行机理，其中基于滤波电容电流反馈的有源阻尼方案不仅控制系统结构简单，而且阻尼效果较好.

除了高频谐波，公共并网点(PCC)附近的非线性负载会使PCC处的电网电压存在背景谐波，不仅影响锁相环性能还会给并网电流带来低阶谐波.

为消除电网电压背景谐波的影响，一般会采用基于静止坐标系的比例积分(PI)控制及电网电压前馈的控制策略[6].这种方法结构简单，且对于抑制电网干扰有不错的效果，但PI控制在静止坐标系中不能有效处理稳态误差，而且对于以低开关频率运行的大功率设备来说，电网前馈对抑制低阶谐波的作用十分有限.另一种方法是基于旋转坐标系的PI+ HC控制[7]，这种方法克服了前一种方法的缺点，但信号解耦、两次dq坐标变换和谐波补偿大大增加了系统的计算量.

在旋转坐标系中，PI控制可转换为静止坐标系中的PR控制[8-9]，若PR控制与HC控制相结合，则能够对指定频率处提供足够大的增益，达到无稳态误差控制，并有效抑制低阶谐波[10].该控制策略结构相对简单，容易实现，文献[11-14]已分析了PR+HC的控制策略在逆变并网系统的运用.

本文在建立单相并网逆变器模型基础上，通过滤波器电容电流的反馈来抑制LCL滤波器的谐振，采用PR+HC的控制策略来提高系统的性能，通过对HC给系统带来的影响和系统稳定性的分析来设计控制器的参数.

1　系统模型

图1为并网逆变器结构图，为简化分析，设直流母线电压Udc为常数，UAB为逆变桥的输出电压，Ug为电网电压.忽略器件的寄生电阻，L1，L2和C组成LCL滤波器.LCL滤波器的谐振频率为

图1　LCL并网逆变器

图2为系统的控制框图，其中Kpwm为逆变比例系数；Gpr(s)为PR控制器，Ghc(s)为谐波补偿控制器；Kc为电容电流Ic的反馈系数.

由图2，可得出系统的开环传递函数

(1)

其中Gc(s)=Gpr(s)+Ghc(s).

并网电流Ig(s)可表示为

图 2系统控制框图

Fig.2System block diagram

(2)

其中，

(3)

2　控制策略

2.1PR+HC控制器

(4)

(5)

式(4)为PR控制器的传递函数，其中Kp为比例系数；Kr为基波谐振项的积分增益；ωo=2fo，ωo为基波角频率，fo为基波频率；式(5)为HC控制器的传递函数，它是由若干个谐振项组成的，其中ξ为阻尼系数，决定着谐振项的带宽；Kn为谐波补偿谐振项的积分增益，n为谐波补偿次数；h为最高谐波补偿次数.

2.1.1比例系数Kp

在LCL滤波器中，电容C主要滤除高阶谐波，则在系统的低频段(小于fr/2), LCL滤波器可简化为L型滤波器[15-16]，式(1)可改写为

(6)

T(s)在截止频率fc的幅频为

|T(j2|.

(7)

若截止频率fc足够大，Gc(s)可简化为一个比例环节[17]，则|Gc(j2fc)|≈Kp，令|T(j2fc)|=1，得

(8)

2.1.2系统的稳态误差

在计算稳态误差时，式(6)依然有效.将式(6)代入式(2)，得

(9)

PR控制器在fo处提供了足够大的增益，并网电流Ig(s)与参考电流Iref(s)的相位误差被消除，但仍存在一定的幅值误差.由于|T(j2fo)|≫1，根据式(6)和式(9)，系统的幅值误差为

(10)

其中并网电流Ig，参考电流Iref和电网电压Ug为有效值.

HC控制器中各谐振项对带宽范围外的频段的影响很小，则忽略HC控制器对fo处的影响，|Gc(s)|在fo处可化简为|Gc(j2fo)|≈Kp+Kr.若希望系统的幅值误差η≤ηst，将ηst和|Gc(j2fo)|≈Kp+Kr代入式(10)，得

(11)

2.1.3系统的相位裕度

文献[18]指出，若系统的截止频率大于PR控制器的转折频率，则PR控制器可化简成一个拥有相同积分系数的PI控制器.同理，若系统的截止频率fc大于Gc(s)的转折频率fz，则Gc(s)可化简为

(12)

其中，Gc2(s)积分项的积分系数等于Gc(s)所有谐振项的系数之和.Gc(s)中各谐振项对其带宽以外频段的幅频增益影响不大，因此在截止频率fc或更高频段处Gc(s)弱化成一PI控制器.Gc(s)和Gc2(s)的波特图如图3所示，可以看出在Gc(s)的转折频率fz前后，Gc(s) 和Gc2(s)幅频曲线的斜率都为-20 dB/dec和0 dB/dec，相频曲线都从-90°上升到0°，当fc>fz时，Gc(s)≈Gc2(s).

图3　 Gc(s)与Gc2(s)的波特图(h=7)

根据式(1)和式(12)，系统的相位裕度为

(13)

若希望系统的相位裕度PM≥PMst，将PMst代入式(13)可得

Kx≤Kx_PM=

(14)

2.1.4谐振项积分系数Kr和Kn的约束范围

由式(11)和式(14)可得Kr和Kn的约束关系：

Kr_η≤Kr≤Kx_PM.

(15)

(16)

2.2有源阻尼系数Kc的约束范围

系统开环传递函数的相频曲线在fr处穿越-180°，系统在fr处的幅值裕度可表达为

GM=-20lg|T(j2fr)|.

(17)

由于fr比hfo要大很多，Gc(s)简化为Gc(s)≈Kp.假设Kx=Kr,若希望系统的幅值裕度GM≥GMst，把式(1)和式(8)代入式(17)，把式(8)和式(11)代入式(13)，分别得到Kc关于fc的边界值函数：

(18)

(19)

根据图3可知，Gc(s)的输出与电容电流反馈量的误差为PWM调制的调制波，为了避免调制波与载波多次交截，调制波的最大变化率应低于载波的变化率，因此要求[17]

(20)

3　实例

3.1系统参数

为验证分析的正确性，在MatlabSimulink中搭建一个1kW的模型，具体参数如表1所示.为保证系统稳定，这里取PMst=45°，GMst=3dB,ηst=0.8%.

表1　系统参数

3.2参数选择

3.2.1阻尼系数ξ和最高谐波补偿次数h

根据式(13)式建立PM关于ξ的函数PM(ξ)，则PM′(ξ)<0，PM关于ξ单调递减.文献[19]中取ξ=0.001，Ghc(s)拥有极小的带宽从而提升了谐波补偿的精确度，Ghc(s)对系统的负面影响可以忽略不计，补偿效果也更好.这里为了提高控制器对电网频率的适应性，以文献[20]为标准，电网频率允许偏差±0.5Hz，取ξ=0.01.

若Kn满足式(16)的约束条件，h对PM的影响可以忽略.在低阶谐波中，一般3次、5次和7次谐波为主流谐波，在这里取h=7.

3.2.2截止频率fc与比例系数Kp

根据由式(18)、(19)和(20)可得到Kc关于fc的取值范围，如图4的阴影部分.关于系统fc的限制，有如下两点：

图4　Kc和fc的取值范围

(1) fc>fz：这样可以保证Gc(s)化简的准确性，也不会因Gc(s)谐振项的相频在其谐振频率(nωo)有-180°的相移导致系统的PM过小.

(2) fc≤0.3fr：为避免对fr的谐振阻尼产生干扰，系统带宽不能超过fr.另外，受到采样和调制的延时影响，fc应限制在0.3fr之内[21].

较大的Kp让系统形成较大的带宽，让系统的瞬态响应更快.由式(8)可知Kp与fc近似成正比，fc尽量取大，取fc=0.3fr=1496Hz，得Kp=0.109.fc确定后，由图4可得到Kc的取值范围[Kc_GM, Kc_max].

3.2.3电容电流反馈系数Kc

根据式(13)建立PM关于Kc的函数PM(Kc)，可知PM′( Kc)<0，PM关于Kc单调递减.为了不让系统的PM过小，Kc不宜取大.取Kc=0.15.

3.2.4基波谐振项积分系数Kr

Kr越大，fo处的幅频越大，基波的追踪效果越好，但系统的PM也越小(证明方法与前文类似).另外，由前面分析可知加入HC后，系统的PM会减小.根据式(11)、(13)和(15)可得Kr的取值范围为[Kr_GM, Kx_PM)，取这个范围的中值，Kr= 64.304.

3.2.5谐波补偿谐振项积分系数Kn

根据式(2)，代表电网电压对并网电流的干扰或系统输出导纳的传递函数可表示为

(22)

根据式(16)可得出Kn的取值范围，由式(22)作出输出导纳的幅值|Gadm(s)|在补偿频率150Hz、250Hz和350Hz处关于Kn的曲线图，如图5所示.

图5|Gadm(s)|在150、250和350Hz处关于Kn的曲线图

Fig.5Curvesof|Gadm(s)|asthefunctionofKnat150, 250and350Hz

结合图5，关于Kn的取值有如下分析.

(1) 当Kn<6时， 随着Kn的值越大，各曲线中|Gadm(s)| 值越小，而|Gadm(s)|的值越小代表系统输出导纳越小，系统对低阶谐波的抑制效果越好；当12≥Kn≥6时，除了曲线3的|Gadm(s)|值还在递减外，曲线1和曲线2的|Gadm(s)|值基本稳定.

(2) 当Kn≥12时，曲线3的|Gadm(s)|值也稳定下来.随着Kn的值增大，三条曲线|Gadm(s)|的值基本不变,即Kn的值继续增大对相关频率处低阶谐波的抑制已经十分有限.

(3) 根据(13)式建立PM关于Kn的函数PM(Kn)，PM′( Kn)<0，PM关于Kn单调递减.为保证系统的PM值，Kn的值不宜取太大.

综合上述3点，取K1= K2=6， K3=12.

4　仿真结果

图6是根据上述参数得出的系统开环波特图.可以看出，在添加HC控制后，在谐振项的带宽内谐波补偿频率的幅频都拥有足够大的增益[22]，使系统产生极小的输出导纳去抑制低阶谐波.当控制器是单一PR控制时，系统GM=3.4dB，PM=60.8°；当控制器是PR+HC控制时，系统GM=3.2dB，PM=54.8°.系统的幅值裕度和相位裕度稍有下降，但对系统稳定性的影响不大.

图6　T(s)的波特图

在有背景谐波环境下对系统进行仿真，其中背景谐波环境为在基频(50Hz)电网电压中叠加15V(有效值)的3次谐波电压分量、10V(有效值)的5次谐波电压分量和7V(有效值)的7次谐波电压分量，叠加谐波分量后，电网电压THD为13.82%.图7是系统采用PR+HC控制策略的仿真结果图，并网电流关于采用PR+HC与单一PR控制策略的谐波率见表2.可以看出在电网电压背景谐波的干扰下，PR+HC的控制策略使并网电流的3次、5次和7次的谐波率与总谐波率明显下降.

图7　电网电压Ug，电流ig仿真图(Gc(s)=Gpr(s)+Ghc(s))

Gc(s)HD/%3次谐波5次谐波7次谐波THD/%PR2.102.131.823.90PR+HC0.820.360.132.08

5　结论

本文建立了单相LCL型并网逆变器在静止坐标系下数学模型，分析了PR+HC控制器中各参数对系统性能的影响，并基于系统的幅值裕度、相位裕度和稳态误差，给出了PR和HC控制器的参数设计方法.与单一的PR控制相比，仿真证明了控制器Gc(s)在添加了HC的情况下能有效地消除主流低阶谐波而不影响系统的稳定性.

[1]KOUROS，LEONJI，VINNIKOVD，etal．Grid-connectedphotovoltaicsystems：anoverviewofrecentresearchandemergingPVconverertechnology[J]．IEEEIndustrialElectronicsMagazine，2015，9(1)：47-61．

[2]LISERREM，BLAABJERGF，HANSENS．DesignandcontrolofanLCL-filtered-basedthree-phaseactiverectifier[J]．IEEEIndustrialElectronicsMagazine，2005，41(7) ：1281-1291．

[3]SOSAJ，ESCOBARG，VAZQUEZG，etal.ComparativeevaluationofLandLCLfiltersintransformerlessgridtiedconvertersforactivepowerinjection[C]∥2014IEEEInternationalAutumnMeetingonPowerElectronicsandComputing．Mexico:IEEE, 2014：1-6．

[4] 张兴，曹仁贤，等. 太阳能光伏逆变并网及其逆变控制[M]．北京：机械工业出版社，2010：160-176．

[5] 许津铭，谢少军，肖华锋．LCL滤波器有源阻尼控制机制研究[J]．中国电机工程学报，2012，32(13)：27-32．

XUJM，XIESJ，XIAOHF．ResearchoncontrolmechanismofactivedampingforLCLfilters[J]．ProceedingsoftheCSEE，2012，32(12)：27-33．

[6]WANGXH，RUANXB，LIUSW，etal．Fullfeedforwardofgridvoltageforgrid-connectedinverterwithLCLfiltertosuppresscurrentdistortionduetogridvoltageharmonics[J]．IEEETransactionsonPowerElectronics，2010， 25(2)：3119-3127．

[7]GOLESTANS，MONFAREDM，GUERREROJM．Anovelcontroltechniqueforsingle-phasegrid-connectedinverters[C] ∥2011InternationalConferenceonElectricalEngineeringandInformatics．Bandung:IEEE, 2011：1-6．

[8]TWININGE，HOLMESDG．Gridcurrentregulationofathree-phasevoltagesourceinverterwithanLCLinputfilter[J]．IEEETransactionsonPowerElectronics，2003，18(7)：888-895．

[9]ZMOODDN，HOLMESDG．StationaryframeCurrentregulationofPWMinverterswithzerosteady-stateerror[J]．IEEETransactionsonPowerElectronics，2003，18(7)：814-822．

[10]TEODORESCUR，BLAABJERGF，BORUPU，etal．Anewcontrolstructureforgrid-connectedLCLPVinverterwithzerosteady-stateerrorandselectiveharmoniccompensation[C]∥2004NineteenthAnnualIEEEAppliedPowerElectronicsConferenceandExposition．[S.l.]:IEEE, 2004，1：580-586．

[11]CASTILLAM,MIRETJ,MATASJ,etal.Controldesignguidelinesforsingle-phasegrid-connectedphotovoltaicinverterswithdampedresonantharmoniccompensators[J]．IEEETransactionsonIndustrialElectronics，2009，56(14)：4492-4501．

[12]XUJM,XIESJ,TANGT.Evaluationsofcurrentcontrolofinweakgridcaseforgrid-connectedLCL-filteredinverter[J]．IETPowerElectronics，2013，6(6)：227-234．

[13]XUJM,TANGT,XIESJ.Researchonlow-ordercurrentharmonicsrejectionsforgrid-connectedLCL-filteredinverters[J].IETPowerElectronics, 2014, 7(9): 1227-1234.

[14]LIB,ZHANGM,HUANGL，etal．Arobustmulti-resonantprregulatorusingdirectpoleplacementdesignstrategy[C]∥2013Twenty-EighthAnnualIEEEAppliedPowerElectronicsConferenceandExposition．LongBeach，CA，USA:IEEE,2013：960-966．

[15]DANNEHLJ，WESSELSC，FUCHSFW．Limitationsofvoltage-orientedPIcurrentcontrolofgrid-connectedPWMrectifierswithLCLFilters[J]．IEEETransactionsonIndustrialElectronics，2009，56(6)：380-388．

[16]TEODORESCUR，LISERREM，RODRíGUEZP.Gridconverterforphotovoltaicandwindpowersystems[M]．UnitedKingdom：Wiley-IEEEPress，2011：313-315．

[17]BAOCL，RUANXB，WANGXH，etal．Step-by-stepcontrollerdesignforLCL-typegrid-connectedinvertercapacitor-with-current-feedbackactive-damping[J]．IEEETransactionsonPowerElectronics，2014，29(7)：1239-1253．

[18]CHENXR，RUANXB，YANGDS，etal．Step-by-stepcontrollerdesignofvoltageclosed-loopcontrolforvirtualsynchronousgenerator[C]∥2015IEEEEnergyConversionCongressandExposition．Montreal，QC，2015：3760-3765．

[19]CASTILLAM,MIRETJ,MATASJ，etal．Controldesignguidelinesforsingle-phasegrid-connectedphotovoltaicinverterswithdampedresonantharmoniccompensators[J]．IEEETransactionsonIndustrialElectronics，2009，56(14)：4492-4501．

[20] 中国国网电力科学研究院，中国电力科学研究院．Q/GDW617-2011 光伏电站接入电网技术规定[S]．北京：国家电网公司，2011．

[21]TANGY，LOHPC，WANGP，etal．GeneralizeddesignofhighperformanceshuntactivepowerfilterwithoutputLCLfilter[J]．IEEETransactionsonIndustrialElectronics，2012，59(7)：1443-1452．

[22] 宋修奎， 唐雄民，章云，等.一种适合于单相电压型逆变电源的并联控制方案研究[J].广东工业大学学报，2015，32(2)：79-85.

SONGXK,TANGXM,ZHANGY,etal.Amethordofparallelcontrolforsingle-phasevoltagesourceinverters[J].JournalofGuangdongUniversityofTechnology, 2015,32(2):79-82.

Research on Proportional-Resonant Controller and Harmonic Compensation for Grid-Connected Inverter

Zhang Miao, Su Xie-fei

(School of Automation, Guangdong University of Technology, Guangzhou 510006,China)

The distortion of background harmonic voltage is inevitable in grid-connected system, which would result in serious low-order harmonics in the injected grid current. The low-order harmonics can be suppressed by increasing the switching frequency, while it also brings many high-order harmonics. Based on PR+HC-based current controller, the effects of the current controller parameters on steady-state error, phase margin and gain margin in LCL-type grid-connected inverter are investigated. A method is given to obtain the controller parameters’ region that meets the requirements, and the optimal parameters could be obtained on the basis of analysis in system stability. Simulation results verify the proposed control scheme and parameters design method.

grid-connected inverter; proportional-resonant control; harmonic compensation; LCL filter

2015- 12- 24

广东省自然科学基金资助项目(2015A030313487)

张淼(1968-)，男，教授，主要研究方向为电力电子与电力拖动及可再生能源发电控制技术.E-mail：Bezhangm@gdut.edu.cn

10.3969/j.issn.1007- 7162.2016.05.011

TM615

A

1007-7162(2016)04- 0059- 06