基于蚁群优化与变步长扰动观察算法的光伏系统mppt方法
【专利摘要】本发明公开了一种基于蚁群优化与变步长扰动观察混合算法的局部遮阴光伏系统MPPT策略。所述方法在局部遮阴条件下光伏系统出现多峰值时解决搜索全局最大功率点问题。当适用情况出现时，使用蚁群算法跟踪光伏输出的前期阶段，避免扰动观察法导致系统陷入局部最大功率点的现象出现，同时提高收敛速度。当蚁群算法第二次迭代结束后，算法停止，开始使用变步长扰动观察法跟踪其后期阶段，利用该算法的鲁棒性，避免蚁群算法导致系统出现稳态震荡。最终将更加迅速、精确和稳定的跟踪到局部遮阴光伏系统的全局最大功率点。
【专利说明】
基于蚁群优化与变步长扰动观察算法的光伏系统MPPT方法
技术领域
[0001] 本发明属于光伏发电技术领域，具体设及一种能迅速、精确、稳定地捜索光伏系统 最大功率点的跟踪策略，尤其是基于蚁群优化与变步长扰动观察混合算法在光伏系统局部 遮阴情况下跟踪全局最大功率点的光伏系统MPPT方法。
【背景技术】
[0002] 目前，现有的光伏系统最大功率点跟踪算法分为两大种类，非智能算法与智能算 法。其中非智能算法包括扰动观察法、电导增量法等，运些方法虽然收敛速度较快、稳态震 荡较小，但当光伏系统出现局部阴影情况时，例如日出、日落W及光伏阵列周围的建筑物和 树木等形成的阴影，运些情况会大大降低系统效率，W上非智能算法很容易陷入局部极值 点。智能算法主要包括蚁群算法、粒子群优化算法、模拟退火算法等，具有全局而且高效的 优化性能、通用性强等优点，但系统结构复杂，计算规模较大，致使硬件成本过高，不适用于 工业大规模应用。
[0003] 例如，南昌大学申请的专利文献，公开号为CN104793691A，其公开了一种基于蚁群 算法的局部阴影下光伏阵列全局MPPT方法，技术方案可W概括为利用蚁群算法结合PI控制 器捜索光伏阵列工作于全局最大功率点处的电压，在最佳电压处采用固定小步长的扰动观 察法跟踪光伏阵列的全局最大功率点。由于该方案采用传统蚁群算法，迭代次数过多，计算 量较大，运样会导致收敛速度过慢，收敛时间过长；同时在后期阶段采用小的定步长作为扰 动量，运样同样会导致一个缓慢的收敛过程。其算法仿真图如图1所示。
[0004] 因此，在实际使用过程中，需要一个更加优异的最大功率点跟踪算法来提高太阳 能的利用率，与此同时能增加系统的收敛速度。
【发明内容】
[0005] 针对【背景技术】中提到的问题，本发明提供一种基于蚁群优化与变步长扰动观察混 合算法的局部遮阴光伏系统MPPT方法，用W利用蚁群优化算法的全局捜索能力和变步长扰 动观察算法的局部捜索能力，更加迅速、精确和稳定的跟踪到局部遮阴光伏系统的全局最 大功率点。
[0006] 本发明的目的是W下述方式实现的：
[0007] -种基于蚁群优化与变步长扰动观察算法的光伏系统MPPT方法，包括W下步骤： [000引步骤1、确定蚁群规模i及运动步长；
[0009 ]步骤2、确定蚁群初始位置；
[0010] 步骤3、采集光伏阵列输出电压化V和输出电流Ipv，计算输出功率P，每个妈蚁位置 所对应的输出功率被认为是该位置上的信息素 T ;
[0011] 步骤4、蚁群进行迭代计算，含有高信息素的妈蚁仍停留在原位，其他妈蚁根据公 式(1)调整自己的位置，其中5是蚁群由原位置移动到最大信息素妈蚁位置时的单位向量：
[001。
Cl)
[0013] 步骤5、当第一次迭代结束后，重复步骤3和步骤4完成第二次迭代，找到此时的"最 伴'妈蚁abest所对应的最大功率点为Pbest,妈蚁位置即其所对应的占空比为Dbest,蚁群算法 终止；
[0014] 步骤6、W蚁群算法两次迭代产生的最佳数据Pbest和化est作为初始数据，启动变步 长扰动观察算法，根据系统要求确定功率变化量允许最小值eP和电压变化量允许最小值 eU;
[0015] 步骤7、计算此时功率变化量A P的绝对值是否小于功率变化量允许最小值eP，若 是，转到步骤8，如否，转到步骤9;
[0016] 步骤8、计算此时电压变化量A U的绝对值是否大于电压变化量允许最小值eU，若 是，转到步骤9，如否，转到步骤11;
[0017]步骤9、根据公式(2)确定扰动步长A D:
[001 引 (2)
[0019] 其中a为变步长速度因子，dP= AP，加=AU;
[0020] 步骤10、根据dP是否为正数进行步长调节，若dP为正数，则根据公式（3)计算更新 的占空比;若dP为负数，则根据公式(4)计算更新的占空比：
[0021] D(k)=D(k-l)+AD (3)
[0022] D(k)=D(k-l)-AD (4)。
[0023] 步骤11、获得全局最大功率点。 / 1 V
[0024] 所述步骤1中，蚁群规模i定位范围为6-12,设置初始运动步长So为口％相勺 范围为50-70，第一次迭代运动步长Si为Soe^，第二次迭代运动步长52为5〇6 一2。
[0025] 所述步骤1中，蚁群规模i定位范围为9，设置初始运动步长So为0.06，第一次迭代 运动步长Si为0.06e^，第二次迭代运动步长S2为0.06e一2。
[0026] 所述步骤2中，将蚁群初始位置定位于0.1和0.9之间，区域内被均匀分割成i-1部 分。
[0027] 所述步骤9中，变步长速度因子a设定范围为0.001-0.003。
[0028] 相对于现有技术，本发明具有W下明显优点：
[0029] 本发明捜索速度快，费时少，捜索效率与输出P-U曲线复杂程度无关，具有较强的 适应性，有效避免陷入局部最优点。从仿真实验得到的数据，如图2所示，在相同溫度条件 下，设置两种不同光照强度模式1和模式2,分别W本发明与【背景技术】两种算法进行仿真。图 1为【背景技术】算法，图2为本发明算法。不难看出，该方案的收敛时间较短。
【附图说明】
[0030] 图1是现有的基于蚁群算法的局部阴影下光伏阵列全局MPPT方法的算法仿真图。 [0031 ]图2是本发明的混合算法的算法仿真图。
[0032]图3是本发明的方法的原理框图。
[0033] 图4是本发明的方法的步骤流程图。
[0034] 图5是本发明的应用示例图。
【具体实施方式】
[0035] 下面结合附图及应用实例对本发明进一步介绍。
[0036] 如图3和图4所示，一种基于蚁群优化与变步长扰动观察算法的光伏系统MPPT方 法，包括W下步骤：
[0037] 步骤1、确定蚁群规模i及运动步长;本实施例中，蚁群规模i数量为9,11为60,根据
得出初始运动步长So为0.06，第一次迭代运动步长Si为0. 〇6e^，第二次迭代 运动步长82为0.066一2。
[0038] 步骤2、确定蚁群初始位置;将蚁群初始位置定位于0.1和0.9之间，蚁群被定位于 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9。区域内被均匀分割成 8 部分。
[0039] 步骤3、采集光伏阵列输出电压Upv和输出电流Ipv,计算输出功率P，每个妈蚁位置 所对应的输出功率被认为是该位置上的信息素 T ;
[0040] 步骤4、蚁群进行迭代计算，含有高信息素的妈蚁仍停留在原位，其他妈蚁根据公 式(1)调整自己的位置，其中5是蚁群由原位置移动到最大信息素妈蚁位置时的单位向量： 幽]
…
[0042] 步骤5、当第一次迭代结束后，重复步骤3和步骤4完成第二次迭代，找到此时的"最 佳"妈蚁abest所对应的最大功率点为Pbest，妈蚁位置即其所对应的占空比为Dbest，蚁群算法 终止;本发明的蚁群算法只迭代两次，两次是经过多次试验后得出的最好的结果，收敛速度 快同时计算结果比较精确。
[0043] 步骤6、W蚁群算法两次迭代产生的最佳数据Pbest和化est作为初始数据，启动变步 长扰动观察算法，根据系统要求确定功率变化量允许最小值eP和电压变化量允许最小值 eU；
[0044] 步骤7、计算此时功率变化量A P的绝对值是否小于功率变化量允许最小值eP，若 是，转到步骤8，如否，转到步骤9;
[0045] 步骤8、计算此时电压变化量A U的绝对值是否大于电压变化量允许最小值eU，若 是，转到步骤9，如否，转到步骤11;
[0046] 巧骤9、巧据公式(2)确定扰动步长A D: 1234 (2) 2
[004引其中a为变步长速度因子，用于调整跟踪速度，设定范围为0.001-0.003,优选地， 可设定为 0.002;dP= AP，加=AU; 3 步骤10、根据dP是否为正数进行步长调节，若dP为正数，则根据公式（3)计算更新 的占空比;若dP为负数，则根据公式(4)计算更新的占空比： 4 D(k)=D(k-l)+AD (3)
[0051] D(k)=D(k-l)-AD (4)。
[0052] 步骤11、获得全局最大功率点。
[0053] 图5中使用Boost升压式变换器连接光伏阵列与负载，其主要优点在于对光伏阵列 的电磁干扰较小，驱动电路简单。Boost输出电压被错制在负载两端的电压，通过改变占空 比D就可W改变变换器输入电压，而Boost变化器输入电压即为光伏阵列的输出电压，因此 改变D就能改变光伏阵列工作点的电压，通过基于蚁群优化与变步长扰动观察混合算法，最 终能够将工作点稳定在全局最大功率点。所示算法位于MPPT控制器中。
[0054] W上所述的仅是本发明的优选实施方式，应当指出，对于本领域的技术人员来说， 在不脱离本发明整体构思前提下，还可W作出若干改变和改进，运些也应该视为本发明的 保护化围。
【主权项】
1. 一种基于蚁群优化与变步长扰动观察算法的光伏系统MPPT方法，其特征在于，包括 以下步骤： 步骤1、确定蚁群规模i及运动步长； 步骤2、确定蚁群初始位置； 步骤3、采集光伏阵列输出电压UPV和输出电流IPV，计算输出功率P，每个蚂蚁位置所对应 的输出功率被认为是该位置上的信息素 τ; 步骤4、蚁群进行迭代计算，含有高信息素的蚂蚁仍停留在原位，其他蚂蚁根据公式（1) 调整自己的位置，其中g是蚁群由原位置移动到最大信息素蚂蚁位置时的单位向量：步骤5、当第一次迭代结束后，重复步骤3和步骤4完成第二次迭代，找到此时的"最佳" 蚂蚁abe3St所对应的最大功率点为Pbe3St，蚂蚁位置即其所对应的占空比为Dbe3St，蚁群算法终 止； 步骤6、以蚁群算法两次迭代产生的最佳数据Pb(3St和Db(3St作为初始数据，启动变步长扰 动观察算法，根据系统要求确定功率变化量允许最小值eP和电压变化量允许最小值eU; 步骤7、计算此时功率变化量ΔΡ的绝对值是否小于功率变化量允许最小值eP，若是，转 到步骤8，如否，转到步骤9; 步骤8、计算此时电压变化量AU的绝对值是否大于电压变化量允许最小值eU，若是，转 到步骤9，如否，转到步骤11; 步骤9、根据公式(2)确定扰动步长Δ D:其中α为变步长速度因子，dP= AP，dU= AU; 步骤10、根据dP是否为正数进行步长调节，若dP为正数，则根据公式（3)计算更新的占 空比;若dP为负数，则根据公式(4)计算更新的占空比： D(k)=D(k-l)+AD (3) D(k)=D(k-l)-AD (4)。 步骤11、获得全局最大功率点。2. 根据权利要求1所述的基于蚁群优化与变步长扰动观察算法的光伏系统MPPT方法， (1 ^ 其特征在于，所述步骤1中，蚁群规模i定位范围为6-12,设置初始运动步长_：^ ~/% /+U ， η的范围为50-70，第一次迭代运动步长δ^δο。，第二次迭代运动步长δ2为δ〇θ气3. 根据权利要求2所述的基于蚁群优化与变步长扰动观察算法的光伏系统MPPT方法， 其特征在于，所述步骤1中，蚁群规模i定位范围为9,设置初始运动步长δ〇为0.06,第一次迭 代运动步长Si为0. ，第二次迭代运动步长δ2为0.06e_2。4. 根据权利要求1所述的基于蚁群优化与变步长扰动观察算法的光伏系统MPPT方法， 其特征在于，所述步骤2中，将蚁群初始位置定位于0.1和0.9之间，区域内被均匀分割成i-1 部分。5. 根据权利要求1所述的基于蚁群优化与变步长扰动观察算法的光伏系统MPPT方法， 其特征在于，所述步骤9中，变步长速度因子α设定范围为0.001-0.003。
【文档编号】G05F1/67GK105955394SQ201610482720
【公开日】2016年9月21日
【申请日】2016年6月24日
【发明人】苏海滨, 曹晓, 曹一晓, 常海松, 韩小鹏, 段刚强, 冯利, 郭鸿奇
【申请人】华北水利水电大学