一种改进型分区域光伏MPPT控制策略的研究_李月英.pdf

2023.1Vol.47No.1研 究 与 设 计收稿日期：2022-06-21作者简介：李月英(1982)，女，河南省人，硕士，副教授，主要研究方向为电工理论与新技术、控制理论与控制工程。一种改进型分区域光伏MPPT控制策略的研究李月英，左明鑫(郑州科技学院 电子与电气工程学院，河南 郑州 450064)摘要：针对传统最大功率跟踪控制算法难以解决跟踪速度和稳态震荡的矛盾，根据光伏组件输出电压、电流和功率的关系，把光伏功率-电压(P-U)曲线分为线性区、最大功率区和急速变化区三部分，提出了分区域变步长复合最大功率跟踪控制策略；在线性区和急速变化区通过定步长算法快速追踪到最大跟踪点附近，进入最大功率区后通过牛顿差值法拟合曲线，并计算曲线极值点即为最大功率点。Simulink仿真和实验结果表明：改进型控制策略不仅具有较快的跟踪速度，而且有效地降低了稳态震荡。关键词：跟踪速度；稳态震荡；分区域；牛顿插值；变步长中图分类号：TM 615文献标识码：A文章编号：1002-087 X(2023)01-0112-05DOI:10.3969/j.issn.1002-087X.2023.01.025Research on an improved subregional photovoltaic MPPT controlstrategyLI Yueying，ZUO Mingxin(School of Electronic and Electrical Engineering,Zhengzhou University of Science and Technology,Zhengzhou Henan 450064,China)Abstract：According to the traditional maximum power tracking control algorithm,it is difficult to solve thecontradiction between tracking speed and steady-state oscillation.According to the relationship between the outputvoltage,current and power of the photovoltaic module,the PV power-voltage(P-U)curve is divided into three parts:linear area,maximum power area and rapid change area.A compound maximum power tracking control strategywith variable step size in different regions is proposed.In linear area and rapid change area,the maximum trackingpoint is quickly tracked by fixed step algorithm,and the curve is fitted by Newton difference method after enteringthe maximum power area,and the extreme point of the curve is calculated as the maximum power point.Simulinksimulation and experimental results show that the improved control strategy not only has a faster tracking speed,butalso effectively reduces the steady-state oscillation.Key words:tracking velocity;steady-state oscillation;subregional;Newton interpolation;variable step size太阳能凭借其取之不尽用之不竭且绿色无污染的特点，成为新能源开发的重要方向之一，而光伏最大功率点追踪(MPPT)技术就是有效利用太阳能的重要途径。目前，各种 MPPT控制技术层出不穷。文献1-2采用的定步长算法存在跟踪精度和稳态精度难以同时兼顾的矛盾，即步长过大，容易造成稳态震荡，步长过小，追踪速度会变慢；文献3采用自适应面积差法避免了传统算法动静态的矛盾，不但可以实现快速跟踪，而且稳态精度高；文献4-5中提出了改进型电导增量法(INC)，可以实时调整跟踪步长，同样具有良好的跟踪速度和稳态精度；文献6采用的变步长INC算法克服了动静态控制不可调和的矛盾；文献7采用的群粒子算法虽然能取得不错的效果，但是这些算法实际中难以实现。针对传统算法的问题，本文根据光伏组件输出电压、电流和功率的关系，把P-U曲线分为线性区、最大功率区和急速变化区三部分，三个区域均采用不同的步长，在离最大功率点(MPP)较远的时候采用大步长以快速接近MPP点，在 MPP点附近，则采用小步长避免震荡。1 光伏组件工程模型和工作特性曲线1.1 光伏组件工程模型光伏组件的工程模型如图1所示。光伏组件的电流-电压关系如式(1)所示。I=Iph-I0exp q(I Rs+U)AKT-1-(I Rs+U)Rsh(1)式中：Iph、I分别为光生电流、输出电流；Rsh、Rs分别为旁漏电阻、串联电阻；U 为负载两端电压；K 为玻尔兹曼常数，K=1.381023J/K；q为电子电荷量，q=1.61019C；T、A 分别为光伏电池温度和二极管特性因子。图1光伏组件电路模型1122023.1Vol.47No.1研 究 与 设 计标准情况下开路电压(Uoc)、短路电流(Isc)、最大功率点对应电压(Um)、最大功率点对应电流(Im)，可以从光伏电池厂家处获取，简化式(1)，可得式(2)式(4)。I=Isc-1-C1exp(UC2Uoc)-1(2)C1=(1-ImIsc)exp(-UmC2Uoc)(3)C2=(UmUoc-1)ln(1-ImIsc)1(4)1.2 光伏组件工作特性曲线标准情况(1 000 W/m2，25)下，得到 I-U、P-U、dP/dU-U和dI/dU-U的光伏组件输出特性曲线，如图2所示。从图 2(c)可以看出，在 P-U曲线左侧，P-U曲线几乎为一条斜率恒定的曲线，接近 MPP 时，曲线表现为上凹特性，在MPP附近曲线近似为开口向下的抛物线，在MPP的右半侧曲线很陡峭，变化较急剧。因此，根据曲线输出特性，把 P-U曲线分为三个部分：线性区、最大功率区和急速变化区。2 改进型MPPT控制策略2.1 光伏曲线区域的判断由图2可以看出，通过判断实际工作点的位置，选择不同的步长。这里根据输出功率对输出电压的倒数来选择步长，由dP/dU=I+UdI/dU得：(1)当|I|UdIdU，即满足IU|dIdU时，dPdU I，则工作在线性区；(2)当|I|UdIdU，即满足IU|dIdU时，dPdU UdIdU，则工作在急速变化区；(3)如不满足条件(1)和(2)，则说明工作在最大功率区。2.2 不同区域的步长选择(1)线性区步长选择由图 2(c)可知，MPP右侧曲线斜率绝对值是左侧曲线斜率绝对值的4倍左右，由于线性区功率变化相对缓慢，因此左侧要想快速追踪 MPP，步长应设置为右侧的 4倍。设急速变化区的步长为U，则线性区步长为4U，但是如果线性区的步长4U太大，会导致在线性区末端处直接跨越最大功率区到急速变化区，因此线性区的步长4U不能大于最大功率区的跨度。(2)急速变化区步长选择由前面的分析可知，急速变化区曲线较为陡峭，因此选择小步长为 U，防止在追踪过程中使实际工作点在线性区和急速变化区之间来回徘徊，造成能量损失和功率波动。(3)最大功率区步长选择由图 2(c)可知，在最大功率区范围内，P-U曲线近似为抛物线，在这个区域范围内抛物线上选取 3个采样点，根据这 3个采样点，利用牛顿差值法得到曲线函数多项式表达式，进一步通过多项式表达式求出极值点，即为最大功率点。选取的3个采样点分别为A1x1，f(x1)、A2x2，f(x2)和A3x3，f(x3)，可写出基于牛顿插值法的多项式公式：L(x)=l1(x)f(x1)+l2(x)f(x2)+l3(x)f(x3)(5)其 中：l1(x)=(x-x2)(x-x3)(x1-x2)(x1-x3)；l2(x)=(x-x1)(x-x3)(x2-x1)(x2-x3)；l3(x)=(x-x1)(x-x2)(x3-x1)(x3-x2)图2光伏组件输出特性曲线1132023.1Vol.47No.1研 究 与 设 计式(6)可进一步写出其对应的一般式：L(x)=ax2+bx+c(6)其中：a=f(x1)(x1-x2)(x1-x3)+f(x2)(x2-x1)(x2-x3)+f(x3)(x3-x1)(x3-x2)(7)b=f(x1)(x2+x3)(x1-x2)(x1-x3)+f(x2)(x1+x3)(x2-x1)(x2-x3)+f(x3)(x1+x2)(x3-x1)(x2-x3)(8)c=f(x1)x2x3(x1-x2)(x1-x3)+f(x2)x1x3(x2-x1)(x2-x3)+f(x3)x1x2(x3-x1)(x3-x2)(9)式(6)求导可得式(7)：L(x)=2ax+b=0(10)可得最大功率点对应的电压为xmpp：xmpp=-b2a(11)由式(6)和式(11)可得，当 xmpp=b/2a，所拟合近似抛物曲线的极大值点为 Lmppb/2a，(4acb2)/4a，此点就是 MPP的理论位置。因此，在进入最大功率区域内部时，只需要选择三个点，通过拟合曲线计算，就可以直接得到给定电压 xmpp，从而精确定位出MPP。但是这三个点的选取不是随意三点，因为如果三点全部位于 MPP的一侧(例如全部在左侧)，那么拟合的曲线将计算不出MPP对应位置。所以，选取的三点要分布于 MPP的两侧，就要加上三点的选取规则：被选取的三个点需要满足情况1或情况2。情况1：(U1U2U3)&(P2P1，P2P3)(12)情况2：(U1U2U3)&(P2P1，P2P3)(13)三个点只要满足情况1或情况2，就能利用牛顿插值来拟合曲线，从而完成计算。2.3 改进型MPPT控制策略工作原理实现改进型 MPPT控制工作原理如图 3所示。由图 3可知，首先初始化电压、电流和功率及其变化量的大小；在线性区和急速变化区的扰动方向由扰动法来确定，如果PU0，则实际工作点在 MPP左侧，保持原扰动方向继续扰动；如果 PU0，则实际工作点在 MPP 右侧，需要向原扰动方向的反向进行扰动；设定功率调整最小值U和电压调整最小值P，当功率和电压变化很小时，可认为已经追到MPP，此时实际工作点保持原地不动；当功率和电压变化超过设置阈值图3改进型MPPT控制流程1142023.1Vol.47No.1研 究 与 设 计时，判断实际工作点所在区域，实行区域分步长控制。(1)当|dI/dU|I/U|时，判定为急速变化区，采用步长 step实现追踪，由于急速变化区曲线较为陡峭，因此步长step不能太大。(3)当条件(1)和(2)都不满足时，说明此时进入了最大功率区，输出此时电压Uk，以Uk为基准，在0.93Uk和1.04Uk范围内，选取满足(U1U2U3)&(P2P1，P2P3)或(U1U2U3)&(P2P1，P2P2)的三个采样点：(U1，P1)、(U2，P2)和(U3，P3)，通过这个三个采样点利用牛顿插值法可以拟合出一个具有峰值的曲线。结合式(7)和式(8)，得到变量a和b的值，并且一次性定位到极值点 Lmppb/2a，(4acb2)/4a(即 MPP 处)。由于外界环境实时变化，在 Lmppb/2a，(4acb2)/4a保持周期时间T，不断重复图3中的步骤。3 Simulink仿真和DSP实验平台验证3.1 Simulink仿真结果分析在 Matlab R2017a/Simulink中搭建 MPPT仿真平台，如图4 所示。光伏电池参数：Uoc=24.1 V，Isc=14.2 A，Um=18.92 V，Im=12.24 A。标准状况下(1 000 W/m2，25)，光伏电池输出功率为 231.7 W；600 W/m2、25 时，输出功率为 152.8 W。Boost 电路参数设置：前级电容 C=150 F，后级电容 C1=550F，电感 L=15 mH，负载 Rload=10。Simulink 仿真时间设置为1 s，功率调整最小值P=1 W，电压调整最小值U=0.5 V。在 Simulink平台上，保持温度为 25 不变，通过设置光照强度(S)变化来观察输出功率波形；光照强度 S设置为：00.4 s，S=1 000 W/m2；0.41 s，S=600 W/m2。图5、图6和图7分别为大步长(Step=0.6)扰动观察法、小步长(Step=0.01)扰动观察法和改进型MPPT算法控制条件下的Simulink仿真输出功率波形。由图5可以看出，当采用大步长扰动法时，启动阶段到达新稳态用时 0.06 s，跟踪速度较快，稳态震荡较大，当 t=0.4 s时，光照强度由 1 000 W/m2突变为 600 W/m2，从突变开始到稳定后，用时约为0.041 s，跟踪速度相对较快，反应迅速。为了对比稳态震荡，统一在 0.80.9 s 时间段内，对波形局部放大，此时功率在 138.4149.7 W之间震荡，并没有真正跟踪到MPP，因为步长大，导致工作点越过MPP，反复在其左右侧来回跳动。由图6可以看出，当采用小步长扰动法时，启动阶段到达新稳态用时0.13 s，跟踪速度较慢，但稳态震荡较小，当 t=0.4s时，光照强度由 1 000 W/m2突变为600 W/m2，从突变开始到稳定后，用时约为 0.07 s，跟踪速度相对较慢。在 0.80.9 s时间段内，对图 6 中的波形局部放大，此时功率在 147.3151.2W之间波动，震荡现象不明显，基本追踪到了MPP。由图7可以看出，当采用改进型MPPT算法时，启动阶段，稳定后用时仅为0.054 s，启动较快。从图7中波形可以看出，稳定后波形震荡较小，当t=0.4 s时，光照强度由1 000 W/m2突变为600 W/m2，从突变开始到稳定后，用时仅为0.038 s，反应迅速。同样，在0.80.9 s时间段内，功率在149.7152.7 W，与小步长扰动法有着几乎相同的稳态效果。3.2 基于TMS320F28379S平台验证基于 TMS320F28379S 平台搭建了基于 Boost 电路的 10kW 光伏系统。Boost参数如下：前级电容 C=1 mF，后级电容C1=1.3 mF，L=2.2 mH，Boost功率管开关频率f=40 kHz。为了便于比较，设置曲线1(Uoc=334 V，曲线适配比FF=0.75，Umpp=250 V，Pmpp=10 kW)和曲线2(Uoc=680 V，FF=0.75，Umpp=510 V，Pmpp=10 kW)，在直流源上把曲线1切换到曲线2。图8所示为三种不同控制条件下的输出结果。图8(a)：固定大步长U=10 V，曲线 1切换到曲线 2，从曲线1对应的稳态到达曲线2对应的新稳态用时约为0.61 s，追踪相对速度较快，两种曲线下对应的稳态，输出电压和电流有一定的震荡，振幅较大。图 8(b)：固定小步长 U=5 V，曲线 1切换到曲线 2，从曲线 1对应的稳态到达曲线 2对应的新稳态用时约为 1.2 s，追踪速度较慢，反应迟钝，两种曲线下对应的稳态，输出电压和电流震荡小，稳态精度相对较高。图4Simulink仿真平台图5固定大步长扰动观察法控制条件下的输出功率波形图6固定小步长扰动观察法控制条件下的输出功率波形图7改进型MPPT控制条件下的输出功率波形1152023.1Vol.47No.1研 究 与 设 计图8(c)：在改进型控制算法条件下，由曲线 1切换到曲线2后，从曲线1对应的稳态到达曲线2对应的新稳态用时仅仅为0.48 s，追踪速度快，有着和固定大步长扰动法相近的动态特性；在保留大步长良好的动态特性基础上，稳态时输出电压和电流波形稳态精度较高。由仿真和实验结果可知：(1)仿真平台，通过判断|dI/dU|和|I/U|的关系，来确定实际工作点所在的区域。开始启动阶段，固定大步长扰动法、固定小步长扰动法和改进型 MPPT 算法分别用时 0.06、0.13 和0.054 s，固定大步长扰动法和改进型MPPT算法具有相近的动态特性，跟踪速度比小步长扰动法提升了一倍左右。稳态方面，改进型算法稳态震荡范围(149.7152.7 W)，幅度为3 W，远远小于固定大步长扰动法稳态震荡范围(138.4149.7 W)，幅度为11.3 W。(2)基于 TMS320F28379S实验平台，在曲线突然切换时，固定大步长扰动法、固定小步长扰动法和改进型 MPPT算法分别用时0.61、1.2和0.48 s，见表1。稳态方面，改进型算法稳态震荡范围(9 8009 980 W)，幅度为 180 W，远远小于固定大步长扰动法稳态震荡范围(9 3159 910 W)，幅度为595 W；固定小步长扰动法稳态震荡范围(9 6789 969 W)，震荡幅度为291 W，但追踪速度较慢。通过实验对比可知改进算法兼顾了跟踪速度和稳态精度。4 结论仿真和实验表明：本文所提出的改进型 MPPT策略通过分区域分步长控制，在线性区和急速变化区，根据曲线特性，选取步长快速，进入最大功率区，保证了跟踪速度，即动态特性；在最大功率区，基于三个采样点，利用牛顿插值法拟合曲线一步计算出最大功率点，从而避免了稳态震荡；解决了传统控制算法对目标控制动静态难以兼顾的矛盾。参考文献：1MOHD-ZAINURI M A A，MOHD-RADZI M A，SOH A C，et alDevelopmentofadaptiveperturbandobserve-fuzzycontrolmaximum power point tracking for photovoltaic boost DC-DCconverterJ.IET Renewable Power Generation，2014，8(2)：183-1942ABDELSALAM A K，MASSOUD A M，AHMED S，et alHigh-performance adaptive perturb and observe MPPT technique forphotovoltaic-based microgridsJ.IEEE Transactions on PowerElectronics，2011，26(4)：1010-10213李鑫，张红，王瑞圳，等基于改进型面积差法 MPPT策略研究J.控制工程，2019，6(6)：1177-11824TEY K S，MEKHILEF S Modified incremental conductanceMPPT algorithm to mitigate inaccurate responses under fast-changing solar irradiation levelJ.Solar Energy，2014(1)：333-3425侯文宝，张刚，田国华一种改进自适应电导增量法在光伏系统的应用J.电气传动，2017(9)：64-696刘洋，张红基于改变比例因子的变步长最大功率跟踪算法J.电源技术，2017，41(5)：777-779,7937MANICKAM C，RAMAN G P，RAMAN G R，et alFireworksenriched P&O algorithm for GMPPT and detection of partialshadinginPVsystemsJ.IEEETransactionsonPowerElectronics，2017，32(6)：4432-4443图8三种不同控制条件下的输出结果表 1 实验平台三种算法对比 算法 固定大步长 固定小步长 改进型 追踪时间/s 0.61 1.20 0.48 稳态电流波动范围/A 17.721.5 19.620.7 19.420.4 稳态电压波动范围/V 495527 508518 505521 折算功率波动范围/W 9 3159 910 9 6789 969 9 8009 980 116