基于简化模型的光伏组件最大功率跟踪仿真

基于简化模型的光伏组件最大功率跟踪仿真张亢a ，胥良a ，朱显辉a ，师楠b ，董芮麟a（黑龙江科技大学a.电气与控制工程学院；b.工程训练与基础实验中心，哈尔滨150022）摘要：为简化变化工况下光伏组件最大功率跟踪的计算，详细推导了变化工况下光伏组件输出特性超越方程的简化求解过程，并以该简化光伏组件输出模型为基础，结合典型的Boost 升压电路，利用扰动观察法构建了基于Matlab 平台的最大功率跟踪模型，给出了辐照度和温度变化条件下光伏组件最大功率跟踪的仿真结果，并与理论计算结果进行了对比分析，结果表明变化工况下所建最大功率跟踪模型的平均相对误差在1%以内，证明了所建模型的有效性和准确性。

关键词：光伏组件；Matlab 模型；最大功率跟踪；扰动观察法中图分类号：TM470文献标志码：A文章编号：1673-2928（2019）02-0017-04收稿日期：2018-11-27作者简介：张亢（1995-），女，河南驻马店人，硕士研究生，主要研究方向为新能源发电。

DOI:10.19329/ki.1673-2928.2019.02.0052019年3月第18卷第2期（总第98期）安阳工学院学报Journal of Anyang Institute of TechnologyMar,2019Vol.18No.2（Gen.No.98）0前言随着石化能源的日益枯竭和气候变暖的逐渐加剧，寻找一种可再生的清洁能源已成为全球能源供给领域的共识，能将储量巨大的太阳能转化为电能的光伏发电装置得到了广泛的重视。

光伏组件是将光能转换为电能的核心设备，因材料和制造工艺的影响，当前其光电转换效率仅为20%左右[1]。

在较低的光电转换效率下，使光伏组件的实际工作点维持在最大功率点附近，是改善光伏输出效率的有效手段，对其进行研究具有较好的理论和现实意义。

国内外专家学者提出了不同的最大功率跟踪策略，具有代表性的主要有：恒定电压法、电导增量法、扰动观察法等[2-3]。

恒定电压法是一种简化的最大功率跟踪方法[4-5]，该方法控制简单且易实现[6]，但忽略了辐照度和温度对光伏组件输出特性的影响，在工况变化时不能可靠地跟踪最大功率。

电导增量法是基于最大功率点处的电导为0的跟踪策略，控制效果好，稳定度高，但对控制系统的性能和传感器的精度要求也较高[2,7]。

扰动观察法通过施加扰动量来判断输出功率的调整方向，其结构简单，被测参数少，也比较容易实现[8]，是目前研究较多、也较为常用的最大功率点跟踪方法[5]。

此外，新型智能算法如粒子群优化算法、神经网络算法、最优梯度优化算法等也在光伏最大功率跟踪中受到了一定关注，但其需要预先给定阈值和权重等参数，很难在工程实践中得到广泛应用[9]。

本文采用扰动观察法作为光伏组件输出特性的最大功率跟踪策略，基于简化的光伏组件输出特性方程，结合典型的Boost 升压电路，在Matlab 仿真平台上搭建由光伏组件模块、MPPT 模块、脉宽调制模块以及升压电路组成的系统模型，给出光伏组件可靠的最大功率跟踪结果，并讨论辐照度和温度对光伏组件输出最大功率点的影响规律，以期为相关的研究提供借鉴。

1光伏组件数学模型光伏组件单指数模型的等效拓扑如图1所示。

图1中I ph 为光生电流（A ）；R S 为串联电阻（m Ω）；R sh 为并联电阻（k Ω）；I d 为流过二极管的电流（A ）；I sh 为流过并联电阻R sh 的电流（A ）；I 为组件输出电流（A ）；V 为组件输出电压（V ）。

由基尔霍夫电流定律可知I =I ph -I d -I sh （1）I d 可计算为I d =I 0{exp [](q (V +IR S ))/(NKT )-}1（2）式（2）中：I 0为流过二极管的暗电流；q 为电子电荷（1.6×10-19C ）；N 为二极管品质因子；K 为玻尔兹曼常数（1.38×10-23J /K ）；T 为组件的绝对温度（K ）。

基于欧姆定律可知，I sh 满足I sh =(V +IR S )/R sh （3）将式（2）、式（3）代入式（1）得I =I ph -I 0{exp [](q (V +IR S ))/(NKT )-}1-V +IR S )/R sh （4）安阳工学院学报2019年式（4）能较为完善地描述光伏组件的工作原理，已被广泛应用于光伏组件理论分析中。

但式（4）为包含五个未知参数I ph 、I 0、N 、R s 、R sh 的超越方程，而且厂家给定的技术手册并不提供上述各参数的具体取值，导致其求解过程比较烦琐，阻碍了在工程实际中的推广应用。

为简化光伏组件输出特性的计算过程，力图在厂商给定短路电流I d 、开路电压V k 、最大功率点电流I m 、最大功率点电压V m 等参数条件下，计算光伏组件的输出特性曲线。

文献[10]给出了以下近似：串联电阻R s 的阻值一般较小，远小于二极管正向导通电阻，故假设I d =I sh ，并联电阻R sh 的阻值较大，因此忽略（V +IR s ）/R sh 项。

基于上述两点近似，光伏组件的I -V 方程可简化为I =I d (I -K I {exp []V /(K 2V k )-}1)（5）在最大功率点时，将V =V m ，I =I m 代入式（5）得I m =I d (I -K I {exp []V m /(K 2V k )-}1)（6）由于在常温条件下exp[V m /(K 2V k )]远远大于1，因此可忽略式中的“-1”项，解出K 1，即K 1=(1-I m /I d )exp []-V m /(K 2V k )（7）将开路条件I =0，V =V k 代入式（5），并将式（7）代入式（5）得[]1-(I m /I d )-1=exp []-V m /(K 2V k )·[]exp(1/K 2)-1（8）由于exp(1/K 2)远远大于1，忽略式中的“-1”项，解出K 2，即K 2=[](V m /V k )-1[]ln(1-I m /I d )-1（9）将k 1、k 2的表达式代入式（5）得到不含未知参数的光伏组件I -V 特性方程。

上述得到的光伏组件I -V 方程只适用于环境温度25℃，辐照度1000W/m 2标准测试环境，在工况变化时，需对上述方程进行如下修正[11]。

ìíîïïïïI d 1=I dref (1+a ΔT )G /G 0I m 1=I mref (1+a ΔT )G /G 0V k 1=V kref (1-c ΔT )ln(e +b ΔG )V m 1=Vmref (1-c ΔT )ln(e +b ΔG )（10）式中：I d 1，I m 1，V k 1，V m 1分别为光伏组件在不同辐照度和温度下的短路电流、最大工作点电流、开路电压、最大工作点电压；△T =T -T 0，△G =G -G 0，T 、G 分别为同一时刻光伏组件的实际温度和辐照度，T 0、G 0分别为标准测试条件下组件的参考温度（25℃）和参考辐照度（1000W/m 2）；e 为自然对数的底数，其值约为2.71828；a 、b 、c 为拟合系数，其典型值推荐为：a =0.0025/℃，b =0.0005/（W/m 2），c =0.00288/℃。

将式（7）、（9）、（10）代入方程（5）便得到任意工况下光伏组件输出的I -V 方程。

2系统仿真模型的建立利用前述的简化I -V 方程，以TDB125×125-72-P 单晶硅光伏组件为例搭建模型，单晶硅光伏组件参数详见文献[12]。

结合上述各模块和Boost 电路原理建立了光伏组件最大功率跟踪仿真系统如图2所示。

-++-++_+-图2光伏组件最大功率跟踪仿真系统图2仿真系统中的输入电压、电流通过零阶保持器进行保持之后，与记忆模块中保持的前一次采样信号做差，分别得到电压和功率的变化量，两者相乘的结果输入到Sign 函数中，当电压变化量与功率变化量的乘积为正时，Sign 函数输出1，否则Sign 函数输出-1，并以Sign 函数的单位输出与仿真步长的乘积判断最大功率跟踪的调整方法，步长设置为0.001。

脉宽调制模块中引入零阶保持器对MPPT 模块输出的参考电压进行保持，并将该值与三角波发生器产生的三角波进行比较，比较结果输入switch 转换开关。

之后输出幅值为1的理想PWM 控制波形控制Boost 电路的占空比，从而控制输出电压的大小。

引入Step 阶跃函数实现突变的辐照度和温度，仿真分别观察温度恒定（25℃）、辐照度突变时和辐照度恒定（1000W/m 2）、温度突变时光伏组件的输出功率。

温度恒定时，通过设定Step 函数使18第二期组件的辐照度在400（W/m 2）、600（W/m 2）、1000（W/m 2）变化；辐照度恒定时，同样通过设置三个Step 函数的参数使组件的温度在25℃、50℃、75℃变化。

3仿真结果与讨论辐照度和温度变化时光伏组件的输出功率分别如图3、图4所示。

图3、图4给出了不同工况下光伏组件的最大功率跟踪结果，在辐照度和温度分别变化时，组件输出功率能跟随相应的变化而变化。

时间（s）功率（W ）图3温度恒定时组件输出功率时间（s ）功率（W ）图4辐照度恒定时组件输出功率为进一步论证所建模型的跟踪精度，利用式（10）给出的辐照度和温度对光伏最大功率的影响规律，对不同工况下组件的最大输出功率的理论值进行计算，并将仿真结果与理论值进行比较，结果如表1所示。

需要澄清的一点是，本文采用简化光伏组件输出特性模型并不会对式（10）的理论计算结果产生影响。

因此，将式（10）的理论计算结果作为比较的参考值是合理的，可行的。

表1最大功率跟踪结果T 恒定G 恒定时间(s)0-3s 3-6s 6-9s 0-3s 3-6s 6-9s实测值(W)61.1596.31174.20173.93171.66168.47理论值(W)61.6796.75174.59174.59172.15168.13相对误差(%)0.840.450.220.380.280.20绝对误差(W)±0.99±0.75±1.2±1.2±1.1±1.6由表1可知，在温度不变，辐照度变化时，光伏组件的输出特性曲线发生变化，所构建简化函数的最大功率跟踪模型能够使组件工作在最大功率点附近；同样，在辐照度不变，温度变化时，所给模型也能达到预期目的。

并且，不同工况下最大功率跟踪模型的仿真结果与理论值的平均相对误差最大不超过1%，仿真结果的绝对误差不超过±2W ，表明所建模型在辐照度和温度变化时均能准确地跟踪其最大功率。