基于灰色關(guān)聯(lián)分析居民屋頂光伏陣列的狀態(tài)檢測

高 天

(安徽理工大學(xué) 電氣與信息工程學(xué)院, 安徽 淮南 232000)

光伏發(fā)電是我國重要的戰(zhàn)略性新興產(chǎn)業(yè)，隨著屋頂分布式光伏電站的急速發(fā)展，電站后期面臨的一些維護(hù)問題也逐漸涌現(xiàn)出來。早期建設(shè)的光伏電站基本上以大中型集中式光伏電站為主，這類大中型集中式光伏電站可以采用基于物理特性的診斷方法，比如：紅外熱成像物理檢測方法[1]，也可以采用基于能量損失的診斷方法[2]等。但是前者檢測設(shè)備昂貴、維護(hù)成本高，后者十分依賴仿真模型的準(zhǔn)確性[3]，并不適合居民小型光伏電站維護(hù)。隨著居民屋頂小型光伏發(fā)電的大力發(fā)展，研究一種高性價比的小型光伏系統(tǒng)故障狀態(tài)檢測方法很有必要。本文基于MATLAB建立4×4的小型光伏陣列模型，模擬實(shí)際屋頂小型光伏系統(tǒng)發(fā)電，實(shí)現(xiàn)光伏陣列各支路的狀態(tài)檢測。

1 光伏陣列仿真模型

太陽能光伏電池的發(fā)電原理被稱為光生伏特效應(yīng)，光生伏特效應(yīng)是一種使導(dǎo)體或半導(dǎo)體在光照作用下與金屬接觸的不同位置形成電位差的現(xiàn)象。第1步是從光子變化成電子即光能量變化成電能量，第2步則是形成電壓，有了電壓差便會流通電流，形成回路[4]。光伏電池的等效電路如圖1所示。

圖1 光伏電池等效電路圖

根據(jù)光伏電池等效電路圖，可以得到光伏電池輸出電流I為

I=Iph-Id-Ish

(1)

式中：Iph是光生電流，在穩(wěn)定的光照條件下，光伏電池可以看作是穩(wěn)定產(chǎn)生電流的恒流源；Id為半導(dǎo)體內(nèi)部P-N結(jié)的結(jié)電流；Ish是流經(jīng)等效并聯(lián)電阻的電流；Rsh等效并聯(lián)電阻代表半導(dǎo)體材料內(nèi)部固有缺陷或者表面不清潔等原因形成的并聯(lián)電阻，通常約為幾千歐姆[5]；I是光伏電池的輸出電流；Rs等效串聯(lián)電阻代表半導(dǎo)體內(nèi)部電阻和電極電阻構(gòu)成的串聯(lián)電阻，正常小于1 Ω；R為電路的負(fù)載；U代表輸出電壓。

進(jìn)一步可推出：

(2)

式中：I0為P-N結(jié)的反向飽和電流；q為電子電荷，1.6×10-19C；UR為太陽能電池的輸出電壓；A為P-N理想特性因子，理想情況取數(shù)值1；K為玻爾茲曼常數(shù)，1.380 6×10-23J/K；T為太陽電池工作絕對溫度值。由于Rsh數(shù)值很大，導(dǎo)致Ish電流很小，通常為了便于搭建模型，一般會忽略Ish[6]。最終公式簡化為

(3)

根據(jù)光伏電池輸出特性的方程，通過Matlab/Simulink軟件搭建出可調(diào)整輻照度(模型中字母S代表輻照度)和溫度(模型中字母T代表溫度)的光伏組件仿真模型，如圖2和圖3所示。

圖2 光伏組件仿真模型

圖3 光伏組件模型的三參數(shù)法底層封裝

Matlab模型搭建完成需要驗(yàn)證模型的正確性。本文取光伏組件的標(biāo)稱值與光伏組件模型的仿真結(jié)果進(jìn)行比較(兩種結(jié)果均在標(biāo)準(zhǔn)測試情況下取值，輻照度取1 000 W/m2，溫度取25 ℃)，如表1所示。

表1 組件標(biāo)稱值與仿真數(shù)據(jù)對比

由表1得出，在標(biāo)準(zhǔn)情況下，光伏組件仿真模型得到的數(shù)據(jù)與組件標(biāo)稱值之間的誤差小于2%，可以應(yīng)用于組件理論研究[7]。由于居民屋頂可以利用的面積較小，所以光伏系統(tǒng)組件的數(shù)量較少。本文通過Matlab/simulink軟件搭建4×4組件排布的光伏陣列，與實(shí)際屋頂小型光伏系統(tǒng)的大小較為吻合。光伏陣列仿真模型圖如圖4所示。

圖4 光伏陣列仿真模型圖

光伏陣列共有4條支路，每一條支路串聯(lián)4個組件，4個支路再并聯(lián)組成1個陣列。根據(jù)搭建的光伏陣列模型，可以得到陣列輸出的I-U和P-U特性曲線。設(shè)置溫度為25℃，輻照度為1 000 W/m2，輸出的特性曲線如圖5所示。

(a)I-U曲線 (b)P-U曲線

由圖5可以看出，仿真光伏陣列的輸出功率較為理想，仿真效果較好，可以應(yīng)用于理論研究中。

2 光伏陣列狀態(tài)檢測

2.1 光伏陣列輸出特性

光伏陣列的輸出特性與光照強(qiáng)度和環(huán)境溫度有關(guān)，設(shè)輻照度1 000 W/m2不變，環(huán)境溫度變化時輸出的I-U和P-U特性曲線如圖6所示。

(a)I-U曲線 (b)P-U曲線

由圖6可知，環(huán)境溫度為-10 ℃時，輸出電流最低，輸出電壓最高，最大功率最大。隨著環(huán)境溫度逐漸升高，輸出電流在逐漸增加，輸出電壓在逐漸降低，同時最大功率點(diǎn)隨著環(huán)境溫度的升高呈現(xiàn)下降趨勢[8]。

環(huán)境溫度保持25 ℃不變，光照強(qiáng)度變化時輸出的I-U和P-U特性曲線如圖7所示。

(a)I-U曲線 (b)P-U曲線

由圖7可知，隨著光照強(qiáng)度的增加，輸出電壓和輸出電流逐漸增加，呈現(xiàn)正向線性規(guī)律變化趨勢。同時隨著光照強(qiáng)度的增加，最大功率以及最大功率點(diǎn)的電壓和電流均在增加[9]。

2.2 狀態(tài)檢測

2.2.1 灰色關(guān)聯(lián)理論

灰色關(guān)聯(lián)分析是對一個系統(tǒng)發(fā)展變化態(tài)勢定量描述和比較的方法，基本思想是通過確定參考數(shù)據(jù)列和若干個比較數(shù)據(jù)列的幾何形狀相似程度來判斷其聯(lián)系是否緊密，它反映了曲線間的關(guān)聯(lián)程度[10]。

由于光伏陣列故障時輸出的特性曲線與理想情況下特性曲線存在較為明顯的相關(guān)性和差異性?？紤]理想情況光伏陣列任一支路輸出的P-U特性曲線作為母序列，故障情況下各支路輸出的P-U特性曲線作為子序列，通過分析其間的相關(guān)程度，根據(jù)相關(guān)性的大小判定出故障線路。

(1)確定母序列與子序列。

確定反映系統(tǒng)行為特征的參考數(shù)列和影響系統(tǒng)行為的比較數(shù)列。反映系統(tǒng)行為特征的數(shù)據(jù)序列，稱為參考數(shù)列，也就是母序列。影響系統(tǒng)行為的因素組成的數(shù)據(jù)序列，稱比較數(shù)列，也就是子序列。

參考序列：

Y=Y(k)|k=1,2...n

(4)

比較序列：

Xi=Xi(k)|k=1,2...n,i=1,2...m

(5)

式中：k對應(yīng)時間段；i對應(yīng)比較數(shù)列中的一行。

(2)變量的無量綱化。

由于系統(tǒng)中各因素列中的數(shù)據(jù)可能因量綱不同，不便于比較或者在比較時難以得到正確的結(jié)論。因此在進(jìn)行灰色關(guān)聯(lián)度分析時，一般都要進(jìn)行數(shù)據(jù)的無量綱化處理。

初值化處理：

(6)

均值化處理：

(7)

(3)計(jì)算關(guān)聯(lián)系數(shù)

關(guān)聯(lián)系數(shù)是比較數(shù)列與參考數(shù)列在各個時刻(即曲線中的各點(diǎn))的關(guān)聯(lián)程度值。

關(guān)聯(lián)系數(shù)ξi(k)：

(8)

記Δi(k)=|y(k)-xi(k)|，公式可簡化為：

(9)

ρ∈(0，∞)，稱為分辨系數(shù)。ρ越小，分辨力越大，一般ρ的取值區(qū)間為(0,1)，具體取值可視情況而定。當(dāng)ρ≤0.546 3時，分辨力最好，通常取ρ=0.5[11]。

(4)計(jì)算關(guān)聯(lián)度

關(guān)聯(lián)系數(shù)是比較數(shù)列與參考數(shù)列在各個時刻的關(guān)聯(lián)程度值，它的數(shù)不止一個。信息過于分散不便于進(jìn)行整體性比較，有必要將各個時刻的關(guān)聯(lián)系數(shù)集中為一個值，即求其平均值，作為比較數(shù)列與參考數(shù)列間關(guān)聯(lián)程度的數(shù)量。

關(guān)聯(lián)度ri公式如下：

(10)

(5)關(guān)聯(lián)度排序

根據(jù)計(jì)算得到的比較數(shù)列關(guān)聯(lián)度，按照大小排序，關(guān)聯(lián)度高的數(shù)列，說明其與參考數(shù)列更相似。假如計(jì)算得到4個比較數(shù)列的關(guān)聯(lián)度，其關(guān)聯(lián)度大小排序?yàn)閞1>r2>r3>r4，則參考數(shù)列就與比較數(shù)列r1的相似度最高，與比較數(shù)列r4的相似度最低。

根據(jù)各支路關(guān)聯(lián)度，設(shè)置關(guān)聯(lián)度對應(yīng)的狀態(tài)評判標(biāo)準(zhǔn)，關(guān)聯(lián)度0.9～1(包括0.9)為健康狀態(tài)、0.8～0.9(包括0.8)為良好狀態(tài)、0.7～0.8(包括0.7)為合格狀態(tài)、0.7以下為故障狀態(tài)[12]。

2.2.2 故障狀態(tài)檢測

在光伏陣列實(shí)際運(yùn)行時，會存在影響輸出特性的干擾因素。比如，陰影遮擋、支路短路、支路開路、組件老化和灰塵遮擋等因素，都會導(dǎo)致光伏陣列的輸出異常[13]。這些干擾因素不僅會影響光伏陣列的整體輸出效率，甚至還會產(chǎn)生安全隱患，嚴(yán)重時可能導(dǎo)致火災(zāi)等事故的發(fā)生。

圖8 光伏陣列布局圖

(1)斷路故障

在Matlab/Simulink仿真中可以通過設(shè)置支路串聯(lián)電阻阻值為無窮，模擬支路斷路故障。設(shè)置支路1斷路和支路1與支路3斷路2種斷路故障，與無故障時的輸出特性曲線進(jìn)行對比，如圖9所示。

(a)I-U曲線 (b)P-U曲線

由圖9可知，隨著光伏陣列支路斷路情況的不斷惡化，陣列的輸出電流呈現(xiàn)線性遞減，同時最大功率點(diǎn)逐漸降低,說明斷路故障對整體的光伏陣列會造成較大的惡劣影響，嚴(yán)重影響發(fā)電效率。

對支路1和支路3斷路時進(jìn)行灰色關(guān)聯(lián)分析，檢測能否找出故障支路。仿真設(shè)置支路1和支路3斷路，其他支路正常運(yùn)行。由Matlab/Simulink仿真得到4條支路的曲線數(shù)據(jù)，數(shù)據(jù)處理結(jié)果如表2所示。

表2 斷路故障分析

由表2可知，支路1和支路3為仿真設(shè)置的斷路故障，關(guān)聯(lián)度分析得出，支路1和支路3的關(guān)聯(lián)度最低，與母序列誤差達(dá)到約30%，誤差較大，根據(jù)關(guān)聯(lián)度評判標(biāo)準(zhǔn)可以判定為故障支路，與預(yù)期結(jié)論一致。

(2)短路故障

仿真時通過將支路1組件并聯(lián)阻值為0的電阻，模擬組件短路的情況。設(shè)置4種情況，分別是組件1短路、組件1和組件2短路、組件1～組件3短路、組件1～組件4短路，如圖10所示。

(a)I-U曲線 (b)P-U曲線

由圖10可知，當(dāng)組件1短路時，由于陣列電流失配，導(dǎo)致I-U和P-U曲線出現(xiàn)雙階梯雙峰狀，P-U曲線存在2個最大功率點(diǎn)，同時輸出功率較正常情況下降；當(dāng)組件1和組件2短路時，I-U和P-U曲線的拐點(diǎn)出現(xiàn)的更早，陣列失配更加嚴(yán)重，其輸出功率衰減更快；當(dāng)組件1、組件2和組件3短路時，曲線拐點(diǎn)出現(xiàn)的趨勢更早，輸出功率下降更嚴(yán)重；當(dāng)支路1的4個組件都短路時，這時并沒有拐點(diǎn)出現(xiàn)，因?yàn)楣夥嚵兄?整個支路都被短路，反而陣列不會失配，但是光伏陣列短路電流下降，同時輸出功率大幅衰減。

灰色關(guān)聯(lián)分析支路1存在一個組件短路時的故障情況。各支路分別存在4個組件，各支路串聯(lián)電阻與組件參數(shù)如表3所示。處理仿真數(shù)據(jù)后得到結(jié)果，如表4所示。由表4可以看出，支路1的關(guān)聯(lián)度最低，對比其他支路明顯下降，由于該支路只有一個組件短路，所以下降幅度較小。其余支路由于串聯(lián)電阻阻值與母序列串聯(lián)電阻不一致，其關(guān)聯(lián)度有所變化，隨著支路串聯(lián)電阻的逐漸增大，關(guān)聯(lián)度呈現(xiàn)逐漸減小的趨勢。根據(jù)評判標(biāo)準(zhǔn)，可判定支路1狀態(tài)為合格，結(jié)論與預(yù)期一致。

表3 各支路串聯(lián)電阻與組件狀態(tài)

表4 短路故障分析

(3)陰影遮擋故障

光伏組件實(shí)際運(yùn)行時，可能存在組件灰塵過多或者附近較高建筑造成的陰影遮擋等問題，都會導(dǎo)致組件發(fā)電效率降低，所以需要進(jìn)行光照強(qiáng)度不同時的仿真模擬分析。環(huán)境溫度保持25 ℃不變，設(shè)置4種不同情況下的輻照度模擬陰影遮擋進(jìn)行輸出特性分析，光伏陣列I-U和P-U特性曲線如圖11所示。

(a)I-U曲線 (b)P-U曲線

由圖11(a)可知，第1種情況，支路1輻照度取400 W/m2，其余支路輻照度取1 000 W/m2不變，對比正常情況，其輸出電流降低；第2種情況，支路1和支路2輻照度取400 W/m2，其余支路輻照度不變，此時輸出電流對比正常情況的電流下降幅度更大；第3種情況下，支路1輻照度取400 W/m2，支路2輻照度取600 W/m2，支路3輻照度取800 W/m2，支路4輻照度取1 000 W/m2，此時光伏陣列輸出電流與第2種情況相似，因?yàn)閮煞N情況光伏陣列的輻照度總和是一致的，所以得到的I-U曲線較為相似，但是后一種情況的電流下降較快，分析原因?yàn)樵撝匪〉妮椪斩容^低以及4條支路的輻照度均不一樣，導(dǎo)致電流振蕩較大，下降較快。隨著輻照度的降低，光伏陣列的輸出電流呈現(xiàn)逐漸下降的趨勢。由圖11(b)可知，隨著輻照度升高，輸出功率與輻照度呈現(xiàn)一種正向線性變化，光伏陣列輸出功率逐漸遞增，最大功率點(diǎn)也在逐漸提高。在實(shí)際應(yīng)用中，應(yīng)該盡量避免光伏陣列附近建筑物造成的陰影遮擋，尤其是居民屋頂小型光伏系統(tǒng)，由于居民房屋附近環(huán)境較為復(fù)雜更要避免陰影遮擋，提高屋面利用率，提升發(fā)電量。

設(shè)計(jì)2組陰影遮擋情況進(jìn)行灰色關(guān)聯(lián)分析。第1組支路1輻照度取400 W/m2，其余支路取1 000 W/m2；第2組支路1取400 W/m2、支路2取600 W/m2、支路3取800 W/m2和支路4取1 000 W/m2。

(a)第1組陰影遮擋故障分析

光伏陣列設(shè)置的各支路串聯(lián)電阻與輻照度數(shù)值如表5所示。根據(jù)設(shè)置數(shù)值做仿真分析，各支路數(shù)據(jù)處理結(jié)果如表6所示。由表6可知，較其他支路關(guān)聯(lián)度，支路1由于輻照度較小，關(guān)聯(lián)度最低，誤差約10%，故健康狀態(tài)評為良好。

表5 各支路串聯(lián)電阻與輻照度的數(shù)值

表6 陰影遮擋故障分析

(b)第2組陰影遮擋故障分析

光伏陣列設(shè)置的各支路串聯(lián)電阻與輻照度數(shù)值如表7所示。各支路數(shù)據(jù)處理結(jié)果如表8所示。由表8可知，支路1和支路2由于輻照度為400 W/m2，關(guān)聯(lián)度較低，誤差較大。支路1和支路3串聯(lián)電阻相同，輻照度不同，對比得出輻照度較大的支路3關(guān)聯(lián)度更高，誤差更小；支路2和支路4輻照度相同，串聯(lián)電阻不同，依然是輻照度較大的支路4關(guān)聯(lián)度更高，說明輻照度對組件輸出功率影響較大，同時輻照度的變化，可以通過關(guān)聯(lián)度分析判定得出結(jié)論，支路1和支路2的健康狀態(tài)均為良好。

表7 各支路串聯(lián)電阻與輻照度的數(shù)值

表8 陰影遮擋故障分析

(4)組件老化故障

光照強(qiáng)度和環(huán)境溫度保持不變，光伏組件老化可以通過調(diào)整支路串聯(lián)電阻的大小進(jìn)行模擬。光照強(qiáng)度取1 000 W/m2，環(huán)境溫度取25 ℃，支路1串聯(lián)電阻分別設(shè)為4 Ω、7 Ω、10 Ω和13 Ω時進(jìn)行光伏陣列組件老化仿真分析。輸出特性曲線如圖12所示。

(a)I-U曲線 (b)P-U曲線

由圖12可知，隨著組件逐漸老化，開路電壓和短路電流并沒有變化，但I(xiàn)-U輸出特性曲線的切線斜率絕對值逐漸變大，說明輸出電流衰減的速率加快。由P-U曲線圖可以看出，隨著組件老化程度的加深，最大輸出功率衰減越來越嚴(yán)重，同時最大功率點(diǎn)的電流和電壓值也在逐漸變小，說明光伏組件的老化對光伏陣列的輸出功率有較大的影響，在實(shí)際運(yùn)行中，一旦檢測到光伏組件老化程度加深，應(yīng)該立即更換組件，減少發(fā)電量損失[14]。

組件老化是較難辨別的故障，因?yàn)榻M件老化是一個漫長持續(xù)的過程，實(shí)際測試中可能短期無法判定組件的老化狀態(tài)。為測試灰色關(guān)聯(lián)分析故障線路的效果，設(shè)置兩組支路組件老化情況進(jìn)行分析對比，輻照度和環(huán)境溫度此時取定值保持不變。

(a)第1組

各支路串聯(lián)電阻數(shù)值如表9所示。各支路數(shù)據(jù)處理結(jié)果如表10所示。由表10可知，支路1串聯(lián)電阻值最大，關(guān)聯(lián)度最??；支路2串聯(lián)電阻值最小，關(guān)聯(lián)度最大。隨著支路2到支路4的串聯(lián)電阻逐漸增大，關(guān)聯(lián)度逐漸較小，與串聯(lián)電阻的變化一致，說明灰色關(guān)聯(lián)分析可以辨別串聯(lián)電阻的微小變化，支路1和支路3的健康狀態(tài)評估為良好。

表9 各支路串聯(lián)電阻數(shù)值

表10 組件老化故障分析

(b)第2組

其余支路串聯(lián)電阻不變，支路2串聯(lián)電阻數(shù)值變化為10 Ω。各支路串聯(lián)電阻數(shù)值如表11所示。各支路數(shù)據(jù)處理結(jié)果如表12所示。由表12可知，隨著支路2串聯(lián)電阻大幅增加，其關(guān)聯(lián)度大幅降低，誤差增加至約22%，與上表得出結(jié)論一致，支路2健康狀態(tài)評定為合格，說明關(guān)聯(lián)度分析可以評估組件的老化程度，為光伏陣列的狀態(tài)評估提供一定的依據(jù)。

表11 各支路串聯(lián)電阻的數(shù)值

表12 組件老化故障分析

3 結(jié) 語

基于Simulink/Matlab軟件搭建4×4組件排布的光伏陣列模型，分析光伏陣列I-U和P-U輸出特性曲線，仿真模擬光伏陣列各支路出現(xiàn)故障的情景，測試灰色關(guān)聯(lián)分析故障效果。光伏陣列支路1和支路3斷路時，通過關(guān)聯(lián)度分析可以找出斷路支路，評估出支路的健康狀態(tài)；支路1中存在1個組件短路時，可以判斷出光伏陣列的短路支路；各支路的輻照度變化時，通過關(guān)聯(lián)度的變化可以辨別出陰影遮擋對光伏陣列的影響程度，評估出健康狀態(tài)；組件老化故障分析時，設(shè)置的串聯(lián)電阻阻值較大，易于分析關(guān)聯(lián)度的變化趨勢，還需要考慮當(dāng)電阻阻值變化更微小時，能否分辨出組件的老化程度。通過灰色關(guān)聯(lián)分析，可以較為便捷地評估出光伏陣列各支路的健康狀態(tài)，為居民屋頂小型光伏系統(tǒng)的健康狀態(tài)檢測提供一定的依據(jù)。