梁智超，李夢(mèng)達(dá)，鄭旭彬，姚林萍

(上海電機(jī)學(xué)院電氣學(xué)院，上海 201306)

太陽(yáng)能是世界上最豐富的清潔能源之一，研究太陽(yáng)能光伏發(fā)電對(duì)促進(jìn)“碳達(dá)峰，碳中和”[1]具有積極意義。近年來(lái)，通過(guò)最大功率點(diǎn)跟蹤(maximum power point tracking，MPPT)的控制來(lái)提高光伏電池陣列的太陽(yáng)能轉(zhuǎn)換效率備受關(guān)注，其主要挑戰(zhàn)在于如何在非均勻光照下的多峰值P-V 曲線(xiàn)實(shí)現(xiàn)快速精確的跟蹤。目前已有很多學(xué)者研究不同的MPPT 算法，文獻(xiàn)[2]提到的擾動(dòng)觀察法和電導(dǎo)增量法在在步長(zhǎng)選擇上容易出現(xiàn)問(wèn)題，對(duì)于局部陰影或者外界環(huán)境劇烈變化下的多峰值的P-V 曲線(xiàn)容易陷入局部功率極值點(diǎn)的情況；文獻(xiàn)[3]提到的粒子群算法能夠解決局部陰影情況下的跟蹤精度不足問(wèn)題，但是個(gè)體位置變化隨機(jī)性大，搜索過(guò)程中功率震蕩較大，在1 s附近才趨于穩(wěn)定，收斂慢，且參數(shù)設(shè)置較為復(fù)雜，在變化的環(huán)境下可能會(huì)導(dǎo)致MPPT 失效；文獻(xiàn)[4]提到的和聲搜索算法對(duì)多峰值曲線(xiàn)全局尋優(yōu)能力強(qiáng)，但收斂速度依賴(lài)于算法參數(shù)的設(shè)置。

本文提出一種基于模糊和聲搜索(fuzzy harmony search，F(xiàn)HS)算法的光伏MPPT 控制方法，通過(guò)模糊控制的思想優(yōu)化和聲搜索算法MPPT，根據(jù)多峰值光伏電池功率特性曲線(xiàn)變化規(guī)律作為專(zhuān)家經(jīng)驗(yàn)，設(shè)置相應(yīng)的模糊規(guī)則，在MPPT 中自適應(yīng)地平衡全局搜索與局部搜索，提高搜索速度和精度。該算法全局功率搜索能力強(qiáng)，對(duì)MPP 的跟蹤精度高，收斂速度快且更加穩(wěn)定。

1 光伏電池陣列

1.1 光伏電池?cái)?shù)學(xué)模型

光伏電池陣列的輸出特性易受光照強(qiáng)度、環(huán)境溫度及工作負(fù)載影響[4]。單個(gè)PV 的等效電路圖如圖1 所示。

圖1 光伏電池等效電路圖

其輸出特性方程如式(1)：

式中：Ipv、Ido、Ish分別為光生電流、輸出電流、二極管反向飽和電流；Rs、Rsh分別為等效串聯(lián)和并聯(lián)電阻；K、q、A分別為玻爾茲曼常數(shù)、電子電荷常數(shù)、光伏電池特性常數(shù)；T為工作溫度。

1.2 光伏電池陣列輸出特性分析

為突出本文所研究的局部陰影情況(PSC)對(duì)光伏陣列輸出的影響，在四塊光伏電池上給定非均勻的光照強(qiáng)度仿真，在給定的標(biāo)準(zhǔn)溫度25 ℃下進(jìn)行仿真。為了驗(yàn)證本文的MPPT 算法是否能在光照突變后依舊能快速跟蹤到全局最大功率點(diǎn)，故設(shè)置兩組不同的光照強(qiáng)度分別為：工況一PSC1：1 000、1 000、800 和500 W/m2；工況二PSC2：500、500、400 和250 W/m2。對(duì)應(yīng)的P-V 曲線(xiàn)如圖2 所示。

圖2 光伏電池P-V特性曲線(xiàn)

由圖2 可知，光伏電池陣列在非均勻光照工況下，P-V 特性曲線(xiàn)出現(xiàn)明顯的多峰值現(xiàn)象，常規(guī)的MPPT 算法在多峰值尋優(yōu)過(guò)程中容易誤跟蹤到局部極值功率點(diǎn)，導(dǎo)致MPPT 難以發(fā)揮最佳的效果，降低了光伏電池陣列整體的能量輸出，因此需要研究全局性多峰值快速尋優(yōu)的MPPT 算法。

2 基于FHS 算法的MPPT 算法

2.1 標(biāo)準(zhǔn)和聲搜索算法(HS)

和聲搜索算法[5]是一種啟發(fā)式全局搜索算法，在和聲彈奏中，樂(lè)師們憑借記憶，通過(guò)重復(fù)調(diào)試和聲音調(diào)，彈奏出最優(yōu)美的和聲，該算法由和聲記憶庫(kù)和相關(guān)和聲調(diào)整參數(shù)組成[6]，算法基本流程如下：

(1)由和聲記憶庫(kù)(HM)產(chǎn)生初始解HMS個(gè)，優(yōu)美和聲評(píng)估函數(shù)為f(x)，見(jiàn)式(2)：

(2)算法每輪迭代在和聲音調(diào)保留概率(PHMCR)外產(chǎn)生隨機(jī)新的和聲音調(diào)xnew，見(jiàn)式(3)：

(3)以音調(diào)調(diào)節(jié)概率(PPAR)按一定的音調(diào)調(diào)節(jié)步長(zhǎng)(BW)進(jìn)行擾動(dòng)去調(diào)整現(xiàn)有和聲音調(diào)，見(jiàn)式(4)：

式中：r1、r2∈[0,1]。

(4)通過(guò)和聲(功率大小)評(píng)估，對(duì)比和聲記憶中最差功率值，若新的功率值優(yōu)于最差功率值則進(jìn)行替換，更新到和聲記憶庫(kù)中進(jìn)行功率輸出。

2.2 引入模糊自適應(yīng)因子

在本文的MPPT 算法中，光伏輸出功率對(duì)應(yīng)為和聲尋優(yōu)變量，占空比對(duì)應(yīng)為音調(diào)步長(zhǎng)調(diào)節(jié)變量，將模糊控制思想與和聲搜索算法結(jié)合運(yùn)用于步驟S3 中，根據(jù)多峰值光伏電池陣列的P-V 功率曲線(xiàn)的特性來(lái)設(shè)置模糊規(guī)則，引入模糊自適應(yīng)因子k，在MPPT 中對(duì)功率和聲f(x)=Vpv·Ipv進(jìn)行評(píng)估對(duì)比，通過(guò)模糊自適應(yīng)因子k對(duì)音調(diào)調(diào)節(jié)步長(zhǎng)BW進(jìn)行動(dòng)態(tài)調(diào)整，算法前期較大的BW有利于獲得更好的全局搜索能力，快速搜索到最優(yōu)和聲附近，隨著迭代次數(shù)的增加，根據(jù)設(shè)置的模糊規(guī)則自適應(yīng)調(diào)整BW，最終獲得最優(yōu)和聲。

本文所提出的模糊控制器總共有兩個(gè)輸入量NI和E及一個(gè)輸出量k。第一個(gè)輸入量為迭代次數(shù)NI，由于非均勻光照下所進(jìn)行尋優(yōu)的光伏陣列P-V 曲線(xiàn)通常是前期平緩穩(wěn)定，中后期是呈多峰狀，考慮其曲線(xiàn)特性，NI的模糊論域設(shè)為[0,150]，子集為{L,LM,M,MH,H}，表示迭代次數(shù)NI的遞增；前期音調(diào)調(diào)節(jié)步長(zhǎng)BW較大有利于大范圍搜索，能提高算法的搜索速度和全局能力；隨著迭代次數(shù)NI增加，此時(shí)BW值應(yīng)自適應(yīng)減小，進(jìn)行小范圍精細(xì)搜索。

第二個(gè)輸入量為光伏輸出功率誤差E，由式(5)定義：

式中：Ei為第i次迭代的功率誤差；Pi為第i次迭代搜索到的最大功率值；Pi－1為第i－1 次迭代的較差功率值。二者的誤差輸入到模糊控制器中。為了實(shí)現(xiàn)在不同功率誤差下的快速跟蹤，E的模糊論域設(shè)為[0,0.5]，子集為{L,LM,M,MH,H}，表示功率誤差的遞增；當(dāng)E較大，表示此時(shí)功率點(diǎn)離最大功率點(diǎn)較遠(yuǎn)，輸出較大的k值調(diào)整BW，在大范圍內(nèi)進(jìn)行和聲搜索，保證其全局功率點(diǎn)搜索能力；當(dāng)E較小，表示此時(shí)功率點(diǎn)已接近最大功率點(diǎn)，應(yīng)使k值減小，降低BW值，縮小搜索范圍，讓其實(shí)現(xiàn)精準(zhǔn)的局部搜索，保證其搜索精度。

k定義為模糊自適應(yīng)因子，為模糊控制器的輸出量，用于和聲音調(diào)調(diào)節(jié)步長(zhǎng)BW的自適應(yīng)調(diào)整，達(dá)到平衡搜索速度與跟蹤精度，提高算法的全局搜索能力的目的，關(guān)系如下：

模糊控制器輸出k的模糊論域?yàn)閇0,1]，子集為{VL,L,ML,LM,M,HM,MH,H,VH}，表示BW對(duì)應(yīng)調(diào)節(jié)為很低、低、中低、低中、中、高中、中高、高、很高。

算法初期，當(dāng)離最大功率點(diǎn)遠(yuǎn)時(shí)，需要較大的BW值進(jìn)行全局搜索，應(yīng)使k值較大，隸屬度函數(shù)選擇較緩三角形，取值稍大，搜索范圍廣；隨著迭代的進(jìn)行，算法逐漸靠近全局最大功率點(diǎn)，此時(shí)需要快速減少BW值，減小搜索震蕩，應(yīng)使k值較小，隸屬度函數(shù)選用較陡三角形，取值精細(xì)，搜索精度高。

綜上所述，NI、E、k的隸屬度函數(shù)設(shè)置如圖3 所示。

圖3 NI、E、k的隸屬度函數(shù)

根據(jù)上述經(jīng)驗(yàn)，模糊控制遵循的原則如表1 所示。

表1 模糊規(guī)則表

改進(jìn)后結(jié)合模糊自適應(yīng)優(yōu)化參數(shù)的和聲音調(diào)更新公式為：

2.3 基于FHS 算法的MPPT 算法

根據(jù)上述內(nèi)容，為解決非均勻光照下的多峰值MPPT 尋優(yōu)問(wèn)題，同時(shí)克服MPPT 尋優(yōu)過(guò)程中的收斂速度與跟蹤精度無(wú)法兼顧的問(wèn)題，提出了基于FHS 算法的MPPT 算法，并將其運(yùn)用于光伏MPPT 中。

基于FHS 的MPPT 控制框圖如圖4 所示。首先，測(cè)量模型中光伏電池陣列的Vpv和Ipv作為輸入量，輸入到HS-MPPT算法模塊中，根據(jù)和聲功率評(píng)估函數(shù)f(x)=Vpv·Ipv的功率值對(duì)比進(jìn)行迭代尋優(yōu)，輸出當(dāng)前占空比D，用于調(diào)節(jié)電壓值Vpv，獲得迭代后和聲功率值；同時(shí)，用HS 算法模塊中的迭代次數(shù)NI與功率誤差E作為輸入量，從模糊控制器中輸出自適應(yīng)因子k來(lái)調(diào)節(jié)BW，對(duì)占空比D進(jìn)行自適應(yīng)調(diào)整，通過(guò)升壓電路獲得最優(yōu)電壓，進(jìn)行全局最大功率點(diǎn)跟蹤。

圖4 基于FHS的光伏MPPT結(jié)構(gòu)模型

基于FHS 算法的MPPT 算法具體步驟如下：

S1.參數(shù)初始化，包括HMS、PHMCR、PPAR、BW和NI。

S2.初始化和聲記憶庫(kù)HM，得到HMS個(gè)初始和聲音調(diào)；式(2)中，初始和聲音調(diào)xi表示初始占空比D，為了保證算法的遍歷性，使占空比D在[0,1]區(qū)間內(nèi)均勻分布，設(shè)置HMS=4個(gè)初始和聲音調(diào)，每個(gè)和聲音調(diào)的初始占空比為{x1,x2,x3,x4}=[0,0.3,0.6,0.9]。

S3.產(chǎn)生新解xnew。以PHMCR的概率在HM的占空比中進(jìn)行初步功率評(píng)估，以(1－PHMCR)的概率獲得HM外的隨機(jī)占空比，然后通過(guò)式(6)的模糊因子以PPAR的概率對(duì)和聲記憶庫(kù)的占空比進(jìn)行步長(zhǎng)為BW的調(diào)節(jié)，獲得新的和聲音調(diào)占空比。

S4.更新和聲記憶庫(kù)HM。評(píng)估函數(shù)為f(xi)=P，新產(chǎn)生和聲功率進(jìn)行評(píng)估，計(jì)算當(dāng)前占空比下功率大??；若當(dāng)前功率值優(yōu)于最差值則進(jìn)行替換，代表搜索到了更優(yōu)的功率點(diǎn)，并將其更新至HM；否則，直接進(jìn)入下一次和聲創(chuàng)作。

S5.檢查算法迭代次數(shù)有無(wú)達(dá)到NI，判斷是否重復(fù)步驟S3 和步驟S4。

圖5 為基于FHS 算法的MPPT 算法流程框圖。

圖5 算法流程框圖

3 仿真分析

3.1 仿真模型

為驗(yàn)證基于FHS 的MPPT 控制策略的高效性和穩(wěn)定性，在MATLAB/Simulink 建立了光伏MPPT 仿真模型。光伏電池模塊參數(shù)為Uoc=43.6 V，Isc=8.35 A，Um=35 V，Im=7.6 A。使用擾動(dòng)觀察法(P&O)、粒子群算法(PSO)、和聲搜索算法(HS)與模糊和聲搜索算法(FHS)分別在局部陰影和光照突變工況下進(jìn)行靜態(tài)與動(dòng)態(tài)的仿真效果對(duì)比。圖6 為基于FHS-MPPT的仿真實(shí)驗(yàn)?zāi)Ｐ汀?/p>

圖6 基于FHS-MPPT的仿真實(shí)驗(yàn)?zāi)Ｐ?/p>

3.2 靜態(tài)仿真分析

靜態(tài)仿真模擬的是光伏電池陣列在局部陰影的非均勻光照下的最大功率點(diǎn)跟蹤，對(duì)多個(gè)最大功率點(diǎn)進(jìn)行全局尋優(yōu)。為驗(yàn)證算法在多峰值P-V 曲線(xiàn)尋優(yōu)的有效性，在標(biāo)準(zhǔn)溫度25 ℃工況PSC1 下進(jìn)行仿真，由圖2 的工況PSC1 的P-V 特性可知，此時(shí)全局最大功率點(diǎn)為684 W，跟蹤精度為算法收斂后跟蹤的平均最大功率點(diǎn)與該工況下的全局最大功率點(diǎn)的比值。各MPPT 算法跟蹤效果如圖7 所示。

圖7 靜態(tài)仿真MPPT效果圖

由圖7 和表2 可知，P&O-MPPT 搜索到P-V 曲線(xiàn)的第一個(gè)波峰505 W 附近，明顯陷入局部最優(yōu)，跟蹤精度約為73.88%，損失功率高達(dá)179 W，雖然收斂用時(shí)最短，為0.14 s，但跟蹤精度最低；PSO-MPPT 算法前期震蕩較大，收斂速度較慢，用時(shí)1.12 s，且收斂后仍存在8 W 左右的小幅搜索震蕩，跟蹤精度約為98.53%；本文提出的FHS-MPPT 算法能做到迅速收斂的同時(shí)兼顧高的跟蹤精度，用時(shí)僅0.25 s 就收斂到功率點(diǎn)681 W處，跟蹤精度高達(dá)99.56%，與全局最大功率點(diǎn)僅相差3 W。相較于傳統(tǒng)算法P&O-MPPT，本文所提算法獲得的功率值有176 W 的大幅提升，能夠兼顧收斂速度和跟蹤精度，在實(shí)現(xiàn)快速跟蹤的同時(shí)也能保持高的跟蹤精度；相較于智能算法PSO-MPPT和HS-MPPT，本文所提算法收斂速度更快，跟蹤到的最大功率值分別提高了7 和2 W，跟蹤精度有所提高，且搜索震蕩小，加強(qiáng)了光伏系統(tǒng)的穩(wěn)定性，提高了太陽(yáng)能的利用率。

表2 靜態(tài)仿真性能對(duì)比表

3.3 動(dòng)態(tài)仿真分析

實(shí)際上，由于受鳥(niǎo)類(lèi)、落葉、移動(dòng)云層等因素的動(dòng)態(tài)影響，光伏電池陣列所處的環(huán)境是復(fù)雜多變的，故設(shè)置光照突變仿真組分析驗(yàn)證算法應(yīng)對(duì)光照突變的跟蹤能力，模擬的是在原有的局部陰影光照工況下，光照強(qiáng)度突然發(fā)生突變的狀況；當(dāng)2 s時(shí)，光照強(qiáng)度從靜態(tài)工況PSC1突變到PSC2的500、500、400和250 W/m2，其P-V 特性曲線(xiàn)對(duì)比見(jiàn)圖2，突變后工況PSC2 的全局最大功率點(diǎn)為375 W，跟蹤效果仿真對(duì)比波形如圖8。

圖8 動(dòng)態(tài)仿真MPPT效果圖

由圖8 和表3 可知，在光照突變情況下，P&O 算法跟蹤效果不佳；PSO 算法的震蕩幅度大，跟蹤精度僅有96.35%，光照突變后跟蹤精度降低為90.40%，且無(wú)法穩(wěn)定收斂；本文所使用FHS 算法能一直保持較高精度的穩(wěn)定跟蹤，光照突變前后的跟蹤精度都高達(dá)99.56%，大大提高了MPPT 的穩(wěn)定性，使算法在動(dòng)態(tài)光照下也能穩(wěn)定快速地跟蹤到最大功率點(diǎn)。

表3 動(dòng)態(tài)仿真性能對(duì)比表

4 結(jié)論

本文通過(guò)把模糊控制的思想結(jié)合到和聲搜索算法MPPT中，引入模糊自適應(yīng)因子k，提出了一種基于模糊和聲搜索(fuzzy harmony search，F(xiàn)HS)算法的光伏MPPT 控制方法，經(jīng)過(guò)靜態(tài)和動(dòng)態(tài)的仿真驗(yàn)證表明，該算法能夠兼顧和聲搜索算法全局性強(qiáng)的優(yōu)點(diǎn)，同時(shí)通過(guò)模糊自適應(yīng)提高了算法的收斂速度，在非均勻光照下能夠保持較高的跟蹤精度，收斂速度快，功率震蕩??；在動(dòng)態(tài)光照工況下也能快速、穩(wěn)定、準(zhǔn)確地搜索到全局最大功率點(diǎn)，為光伏MPPT 多峰值功率尋優(yōu)提供了穩(wěn)定有效的解決方案。