一种光伏发电系统最大功率点智能跟踪方法
【技术领域】
[0001] 本发明涉及新能源开发及应用技术领域，特别涉及一种光伏发电系统最大功率点 智能跟踪方法，该方法可实现外界环境变化时光伏发电系统最大功率点的精确跟踪。
【背景技术】
[0002] 太阳能已成为目前应用最广泛的可利用清洁能源，光伏太阳能技术在诸多太阳能 应用技术中脱颖而出。然而，光伏太阳能电池在发电的过程中具有极其强烈的非线性，并且 极易受到多种环境因素的影响，很难稳定运行于系统最大功率点的位置。因此，如何使光伏 发电系统能够在多变的环境下准确的运行于最大功率点的位置成为了能源有效利用的关 键，对光伏发电系统进行最大功率跟踪（MPPT)显得非常重要。
[0003] 为了提高光伏发电系统的效率，获得尽可能多的电能，国内外对光伏发电系统的 MPPT控制算法进行了广泛的研究，提出了多种MPPT控制算法，主要有：恒定电压法、扰动观 测法、电导增量法。然而，传统的最大功率点跟踪方法总存在着诸如响应速度慢，跟踪精度 不够，系统工作点在最大功率点附近振荡等不足之处，在实际应用中很难取得很好地效果， 对光伏发电系统所带来的优势不明显，无法实现高精度快速跟踪最大功率点的目的，很难 进行大规模推广。近年来随着科学技术的不断发展，智能优化方法步入人们的视野，并不断 成为各个技术领域应用的核心技术手段，利用智能优化的方法解决工程问题已成为了工程
技术领域必然的发展趋势，因此智能优化方法在光伏发电系统的最大功率点跟踪领域的应 用成为了可能，并且成为该领域方向的研究热点。
【发明内容】
[0004] 为了解决现有技术中存在的问题，本发明的目的在于提出一种光伏发电系统最大 功率点智能跟踪方法。当外界环境变化时，其可实现对光伏发电系统最大功率点的精确跟 踪，避免传统跟踪方法在最大功率点附近振荡的缺点，提尚光伏发电系统的输出效率。
[0005] 为达到以上目的，本发明是采取如下技术方案予以实现的：
[0006] -种光伏发电系统最大功率点智能跟踪方法，包括下述步骤：
[0007] 步骤一，构建考虑太阳辐射变化和温度影响时的光伏发电系统中光伏电池阵列数 学模型如下式：
[0008] ΓΤ
[0009] 式中，U为光伏电池阵列电压，I为对应的输出电流，，Ψ2 = (Um/Uoc-l)/ln(l-Im/Isc),Δυ=-β·ΔT-Rs ·ΔΙ,ΔΙ=α·S/Sref ·ΔΤ+(S/Sref-1) ·Isc,ΛΤ= ?；=Τa+t。*S;Sraf为太阳辐射参考值，取lkW/m2,Traf为温度参考值，取25°C; 4为短路电流，1为开路电压，In、Un分别为最大功率点电流和电压，α为电流随温度变化 系数，β为电压随温度变化系数，以上参数都在31^、1'1^条件下获取；Rs为光伏模块串联电 阻；S为总太阳辐射，T。为太阳电池温度，Ta为环境温度，t。为电池模块温度系数； W-MJU^MJ
[0010]步骤二，构建跟踪目标函数(j--VP|(C'丨。聲-w I丨,/(屮』/ ) 给定搜索分段数η= 1〇,电压间隔λ=iWη，初始电压值υ(1)= [λ, 3λ,…，u。。-λ]τ， 初始步长δ=λ/4,初始点Λ-,. =w，（/ = 1,2,….//oor(v/2))，方向因子μ= [1，-1]，速度 因子C2,c2= 2,k= 1，全局比较因子ξ=λ/1〇〇,终止条件ε= 10 3;
[0011] 步骤三，如果F(Χι+μiδ) <F(Xl)，令λ·>λ:.+ ，转步骤四；否则如果 FUi+μ；；δ) <FUJ，令< =x,转步骤四；否则，令x>x(，转步骤四；其中，〖为本 轮比较搜索中的优异值；
[0012] 步骤四，如果) <F({/m)，转步骤五；否则转步骤六；其中，?/广为第k次迭代 后的电压值；
[0013] 步骤五，如果δ>ξ，令17丨叫^4 = 切*):)，f= 1 + 1，转步 骤三；否则令=<:，K=[，糾> -#.1) + ? (?1 -#')，14+1，转步骤三；其 中，?/$为tTdtiF(L/f)所对应的c/f,ITf-为第k+Ι次迭代后的电压值；
[0014] 步骤六，如果δ> ε，则令δ=δ/2，τ:=?/,(?、Lf+ll=(7广，灸=1 + 1，转步 骤三，否则，，其中，1为最大功率点对应的电压值。
[0015] 与现有技术相比，本发明具有以下有益的技术效果：
[0016] 1、本发明所述方法不需要求解目标函数的导数，所以在解决较为复杂的光伏发电 系统最大功率点优化模型时非常有效；
[0017] 2、本发明所述方法同时从多个初始值开始搜索，能快速找到最大功率点位置，动 态响应速度快；
[0018] 3、本发明所述方法在速度函数中考虑局部寻优及全局寻优的综合因素，不易陷入 局部最优解，减少了系统在最大功率点振荡的能量损失。
[0019] 为了更清楚的理解本发明，以下结合附图对本发明作进一步的详细说明。
【附图说明】
[0020] 图1为本发明方法的步骤流程图；
[0021] 图2为本发明方法实施例的总体框图；
[0022] 图3为本发明方法实施例的MATLAB/Simulink仿真控制模型图；
[0023] 图4为本发明方法实施例的最大功率点仿真计算结果图；
[0024] 图5为本发明方法实施例的最大功率点智能跟踪效果图；
【具体实施方式】
[0025] 如图1所示（其中图1中大写的Y是YES的缩写，代表满足所属条件的情况；大写 的N是N0的缩写，代表不满足所属条件的情况），一种光伏发电系统最大功率点智能跟踪方 法，包括下述步骤：
[0026] 步骤一，构建考虑太阳辐射变化和温度影响时的光伏发电系统中光伏电池阵列数 学模型如下式：
[0027]
[0028] 式中，U为光伏电池阵列电压，I为对应的输出电流，ψ =jj1f，Ψ2 = (IVUc-D/lnd-L/IJ，Δυ=-β·AT-RS·ΔΙ，ΔΙ=a.S/Sref ·AT+(S/Sref-1) .Isc，ΛΤ= ?；=Τa+t。*S;Sraf为太阳辐射参考值，取lkW/m2,Traf为温度参考值，取25°C; Is。为短路电流，U。。为开路电压，I 分别为最大功率点电流和电压，α为电流随温度变化 系数，β为电压随温度变化系数，以上参数都在31^、1'1^条件下获取；Rs为光伏模块串联电 阻；S为总太阳辐射，T。为太阳电池温度，Ta为环境温度，t。为电池模块温度系数； W--AV. u-m
[0029] 步骤二，构建跟踪目标函数F 平如队一切+泣―私ψ! /(平凡):， 给定搜索分段数η= 1〇,电压间隔』==^./"，初始电压值υ(1)= [λ，3λ，···，υ^λ]τ4??始步长3 = /4，初始点\=&',1')，（/ = 1,2, - .//〇〃/'(///2))，方向因子4 = [1，-1]，速度因 子C2,c2= 2,k= 1，全局比较因子ξ=λ/1〇〇,终止条件ε= 10 3;
[0030] 步骤三，如果F(Χι+μiδ) <F(Xl)，令磚+MJ，转步骤四；否则如果 Ρ(Χι+μ2δ) <F(Xl)，令转步骤四；否则，令<=七，转步骤四；其中，《为本轮 比较搜索中的优异值；
[0031] 步骤四，如果转步骤五；否则转步骤六；其中，t/f1为第k次迭代 后的电压值；
[0032] 步骤五，如果δ>ξ，令 # ⑴=.<，+l1+q(W+"-1)，Λ4 + 1，转步 骤三；否则令#+11 =心X=#+1'丨-Cl(f/糾丨-+ ，/c4+:l啭 中，L/J)为minF(C/f>)所对应的l/fSC/f+1>为第k+1次迭代后的电压值；
[0033] 步骤六，如果δ> ε，则令δ=δ/2, <=?/广，?/〉?+1??，?=?'+1，转步 骤三，否则，，其中，1]_为最大功率点对应的电压值。
[0034] 实施例：
[0035] 本实施例可以验证本发明的有效性和准确性。图2为实施例的总体框图。选取 的光伏电池阵列在光辐射度为lkW/m2,温度为25°C的标准测试条件下基本参数为：1" = 5· 51A，、Vn= 18. 17V，、Voc= 21. 6V、，Isc= 6· 06A、，α= 7. 3X10 3A/°C，、β= 0· 108V/°C。 基于MATLAB/Simulink工具箱，根据本发明方法构建的光伏电池仿真模型如图3所示。
[0036] 假设随着天气的变化，光辐射度从1000W/m2逐渐降为800W/m2乃至500W/m2,根据 仿真计算得到的理论最大功率点的结果如图4所示，而根据本发明方法得到的跟踪效果如 图5所示。结合图4及图5,可以看出，基于本发明方法，当光照强度发生变化时，该系统能 够快速地做出调节，使其工作在最大功率点附近，稳定性较好，而且很好的完善了恒压法当 环境温度一定时光照强度发生较大改变从而使系统无法进行最大功率点跟踪的现象。
[0037] 本实施例没有详细叙述的部件、工艺及字母表示属本行业的公知部件、和常用手 段及常识，这里不一一叙述。
【主权项】
1. 一种光伏发电系统最大功率点智能跟踪方法，其特征在于：包括下述步骤： 步骤一，构建考虑太阳辐射变化和温度影响时的光伏发电系统中光伏电池阵列数学模 型如下式：式中，U为光伏电池阵列电压，I为对应的输出电流，=VU，！；= T a+t。· S ;Sraf为太阳辐射参考值，取lkW/m 2, Traf为温度参考值，取25°C ;I sc为短路电流，U。。为开路电压，In、Un分别为最大功率点电流和电压，α为电流随温度变化系 数，β为电压随温度变化系数，以上参数都在S"f、Traf条件下获取；1为光伏模块串联电阻； S为总太阳辐射，T。为太阳电池温度，T a为环境温度，t。为电池模块温度系数； 步骤二，构建跟踪目标函数给定 搜索分段数η = 10,电压间隔λ =ιν/η，初始电压值u(1)=[人，3人，~，1]。。-人]\初始 步长δ = λ/4,初始点$ = [/广，（i = 1，2,…floor ( η/2))，方向因子μ = [1，-1]，速度 因子C 2, c 2= 2, k = 1，全局比较因子ξ = λ/1〇〇,终止条件ε = 10 3; 步骤三，如果F(Xl+ μ丨δ ) < F(Xl)，令Λ-;乂V,.十以八，转步骤四；否则如果F(Xl+ μ 2 δ ) <F(Xl)，令-<=七+ /^，转步骤四；否则，令x丨=x,，转步骤四；其中，《为本轮比较搜索中 的优异值；步骤四，如果，:(<):<尸〇7严），转步骤五；否则转步骤六；其中，为第k次迭代后的 电压值； 步骤五，如果 1，转步骤三； 否则令转步骤三；其中， C为min 、所对应的C/f >，为第k+Ι次迭代后的电压值； 步骤六，如果S > ε，则令δ = δ/2，χ:=?7广，04-+1!=0'1^ = 1^1，转步骤三， 否则，=〇，其中，1]_为最大功率点对应的电压值。
【专利摘要】本发明公开了一种光伏发电系统最大功率点智能化跟踪方法，首先，根据光伏发电系统中光伏电池阵列数学模型构建了目标函数；其次，通过从多个电压区间进行智能搜索的方式，快速获得多个最优解；最后，通过在速度函数中引入全局比较因子ξ，实现了局部寻优及全局寻优的兼顾。本发明避免传统跟踪方法跟踪精度不高、容易在最大功率点附近振荡的缺点，可实现外界环境变化时光伏发电系统最大功率点的精确跟踪，动态响应速度快，不易陷入局部最优解，减少了系统在最大功率点振荡时的能量损失，提高光伏发电系统输出效率。
【IPC分类】G05F1/67
【公开号】CN105334901
【申请号】CN201510812726
【发明人】智勇, 郑伟, 拜润卿, 梁福波, 李养俊, 陈仕彬, 乾维江, 郝如海, 张彦凯, 高磊, 邢延东, 祁莹, 魏乔, 章云
【申请人】国网甘肃省电力公司电力科学研究院, 国网甘肃省电力公司, 国家电网公司, 西安电子科技大学
【公开日】2016年2月17日
【申请日】2015年11月20日