基于混合方法的并網(wǎng)光伏發(fā)電量預(yù)測(cè)模型的研究

嘉興職業(yè)技術(shù)學(xué)院，浙江嘉興 314016

一、引言

當(dāng)今世界，伴隨著煤炭、石油等化石能源面臨枯竭狀態(tài)，能源的短缺以及環(huán)境的惡化成為了制約全球經(jīng)濟(jì)發(fā)展的瓶頸。太陽(yáng)能作為自然界中最重要也是最基本的能源，以其分布廣泛、清潔、可再生等特點(diǎn)得到了非常廣泛的應(yīng)用[1]。其中，光伏發(fā)電技術(shù)成為了全球各國(guó)利用太陽(yáng)能的主要方式之一。但是，光伏發(fā)電具有波動(dòng)性、隨機(jī)性和不可控性等缺點(diǎn)，它將對(duì)電網(wǎng)造成不可逆的沖擊[2]。因此，對(duì)光電功率進(jìn)行預(yù)測(cè)能提高電力系統(tǒng)運(yùn)行的安全性和穩(wěn)定性，減輕其對(duì)電網(wǎng)的不利影響，為光伏電站的運(yùn)營(yíng)和維護(hù)提供了可靠的技術(shù)參考。

目前，歐洲各國(guó)對(duì)光伏發(fā)電預(yù)測(cè)的研究非常重視，并取得了重要的成果。以德國(guó)為例，其算法已經(jīng)實(shí)現(xiàn)了中尺度的預(yù)測(cè)，即可對(duì)未來(lái)三天的光伏發(fā)電量進(jìn)行計(jì)算，并達(dá)到了較小的誤差[3]。與此相比，國(guó)內(nèi)的研究集中在短期預(yù)測(cè)上，所用方法傾向于單一的統(tǒng)計(jì)算法，大多尚處于理論研發(fā)階段，還沒(méi)有實(shí)用的預(yù)測(cè)系統(tǒng)[4]。

本文以浙江晶科能源有限公司與德國(guó)Magdeburg大學(xué)、Fraunhofer光伏研究中心合作開發(fā)的10.21m2并網(wǎng)光伏系統(tǒng)為研究對(duì)象，以光伏監(jiān)控系統(tǒng)數(shù)據(jù)庫(kù)的歷史發(fā)電數(shù)據(jù)和氣象數(shù)據(jù)為參考，基于物理和統(tǒng)計(jì)相結(jié)合的預(yù)測(cè)方法，設(shè)計(jì)完成一套可落地應(yīng)用到行業(yè)的光伏預(yù)測(cè)系統(tǒng)。通過(guò)混合方式獲得的光伏系統(tǒng)發(fā)電量預(yù)測(cè)模型，既考慮了系統(tǒng)本身的結(jié)構(gòu)特性和所遵循的物理規(guī)律，又結(jié)合溫度、光照輻射強(qiáng)度等環(huán)境因素進(jìn)行仿真，是一種具有理論和實(shí)際價(jià)值的方法。

二、氣象因素對(duì)光伏發(fā)電的影響

光伏發(fā)電隨時(shí)間和環(huán)境的變化具有較明顯的周期性、波動(dòng)性和間斷性[5]。影響光伏系統(tǒng)發(fā)電量的因素有很多。其中，氣象因素起到了最重要的影響。為了獲得相對(duì)精確的預(yù)測(cè)數(shù)值，首先需要掌握光伏發(fā)電和氣候因素之間的關(guān)系，并對(duì)此作出定量的分析[6]。對(duì)于此既定光伏組件來(lái)說(shuō)，其歷史數(shù)據(jù)自身包含了光伏組件的信息，如轉(zhuǎn)換效率、安裝角度等。而對(duì)于異變的氣象條件，哪些是影響光伏發(fā)電的主要因素，哪些是次要因素，需要通過(guò)一定量的數(shù)據(jù)進(jìn)行統(tǒng)計(jì)、篩選和分析[7]。

1、太陽(yáng)輻射強(qiáng)度對(duì)光伏發(fā)電的影響

太陽(yáng)輻射強(qiáng)度指的是單位時(shí)間內(nèi)垂直投射在地球單位面積上的太陽(yáng)輻射能量[8]，它直接決定了光伏電池的發(fā)電特性。圖1選取在25oC（24.8oC ～25.2oC）恒溫條件下的測(cè)量數(shù)值，獲取不同輻射強(qiáng)度下，光伏系統(tǒng)的發(fā)電量曲線。

由圖1可以得知，并網(wǎng)光伏系統(tǒng)的輸出功率與太陽(yáng)輻射強(qiáng)度的變化趨勢(shì)較為一致。當(dāng)太陽(yáng)輻射強(qiáng)度增強(qiáng)時(shí)，光伏電池的輸出功率升高；當(dāng)太陽(yáng)輻射強(qiáng)度減弱時(shí)，光伏電池的輸出功率也逐漸減小。因此，太陽(yáng)輻射強(qiáng)度是發(fā)電量預(yù)測(cè)模型的一個(gè)重要輸入變量。

2、溫度對(duì)光伏發(fā)電的影響

光伏電池的溫度特性是指環(huán)境溫度和電池自身溫度對(duì)光伏電池發(fā)電性能的影響[9]。實(shí)際工程中，光伏電池的模塊溫度是影響其發(fā)電量的直接因素。圖2選取太陽(yáng)輻照強(qiáng)度為1000W/m2(990W/m2～1100W/m2)恒定值時(shí)的歷史測(cè)量數(shù)據(jù)，獲取不同模塊溫度條件下，光伏系統(tǒng)的發(fā)電量曲線。

由圖2可知，在相同輻射強(qiáng)度環(huán)境下，隨著光伏模塊溫度的升高反而會(huì)引起光伏電池發(fā)電量的逐漸下降。這是由于溫度影響了光伏電池的開路電壓和短路電流。因此，溫度是發(fā)電量預(yù)測(cè)模型的另一個(gè)重要參數(shù)。

三、基于混合方法的預(yù)測(cè)模型

光伏系統(tǒng)出力特性模型是實(shí)現(xiàn)發(fā)電功率準(zhǔn)確預(yù)測(cè)的關(guān)鍵。目前，國(guó)內(nèi)外在光伏發(fā)電預(yù)測(cè)方面主要有物理建模和統(tǒng)計(jì)建模兩種方法[10]。

物理方法著眼于從太陽(yáng)能到電能轉(zhuǎn)化過(guò)程中涉及的各組成裝置的內(nèi)部結(jié)構(gòu)及其物理特性，建立相應(yīng)的數(shù)學(xué)模型。其仿真精度高，但模型復(fù)雜，參數(shù)難解。

統(tǒng)計(jì)方法則是把光伏系統(tǒng)當(dāng)做是一個(gè)“黑箱”，基于歷史氣象和發(fā)電數(shù)據(jù)擬合出功率預(yù)測(cè)模型[11]。此方法相對(duì)簡(jiǎn)便，但卻誤差較大。

本文采用的是物理和統(tǒng)計(jì)相結(jié)合的混合方法。首先根據(jù)系統(tǒng)各組成部分的物理特性建立數(shù)學(xué)參數(shù)模型，再利用電站的歷史數(shù)據(jù)擬合出該模型的參數(shù)，從而對(duì)光伏發(fā)電量進(jìn)行預(yù)測(cè)。

此實(shí)驗(yàn)平臺(tái)的并網(wǎng)光伏系統(tǒng)主要由光伏組件陣列和逆變器構(gòu)成。光伏組件陣列在太陽(yáng)光照的情況下，通過(guò)光生伏特效應(yīng)將照射到電池表面的太陽(yáng)能轉(zhuǎn)化成直流電，通過(guò)逆變器轉(zhuǎn)換成交流電輸出[12]，供交流負(fù)載使用。

1、光伏組件建模

光伏電池陣列是光伏系統(tǒng)的最基本單元和核心部件，其輸出功率具有非線性特征。由前文可知，太陽(yáng)輻射強(qiáng)度和電池溫度是影響其輸出功率的最主要因素[13]。氣象環(huán)境的變化導(dǎo)致了光伏組件無(wú)法持續(xù)工作在最大輸出功率點(diǎn)處，從而在一定程度上降低了光伏系統(tǒng)的能量轉(zhuǎn)換效率。因此，通常在光伏系統(tǒng)中加入MPPT跟蹤，保證光伏系統(tǒng)的輸出功率維持在最大值[14]。

根據(jù)德國(guó)Fraunhofer研究所Sauer教授的研究，光伏組件的直流輸出功率可由最大輸出功率點(diǎn)處工作電流IMPP和工作電壓UMPP之間的經(jīng)驗(yàn)公式表示，并與電池溫度呈一定的線性比例，如下所示：

式中，PDC—光伏電池直流輸出功率；

IMPP—最大功率點(diǎn)工作電流；

UMPP—最大功率點(diǎn)工作電壓；

αP—與溫度相關(guān)的功率因素，值為-0.0036K-1；

TZ—光伏電池模塊溫度；

TSTC—標(biāo)準(zhǔn)溫度，值為25oC。

其中，工作電流IMPP跟光照強(qiáng)度Gm呈二次項(xiàng)的比例，如下式所示：

式中，α—工作電流公式參數(shù)；

Gm—太陽(yáng)輻照強(qiáng)度。

公式中的參數(shù)τ在光照強(qiáng)度較大時(shí)可忽略。在本設(shè)計(jì)中所取數(shù)據(jù)的光照強(qiáng)度均大于100W/m2，因此參數(shù)τ忽略不計(jì)。

工作電壓UMPP跟光照強(qiáng)度Gm間的關(guān)系如下式所示：

式中，a，b，c，d—工作電壓公式參數(shù)，

另外，光伏電池模塊溫度TZ跟光照強(qiáng)度Gm間遵循“kGT模型”，如下式4所示：

式中，kGT，J—電池模塊溫度公式參數(shù)；

TA—環(huán)境溫度。

2、逆變器建模

逆變器是并網(wǎng)光伏系統(tǒng)的關(guān)鍵設(shè)備之一。它是光伏組件陣列與電網(wǎng)之間的接口，通過(guò)開關(guān)元件規(guī)律的進(jìn)行開與關(guān)的操作，實(shí)現(xiàn)將光伏陣列產(chǎn)生的直接電轉(zhuǎn)換為交流電輸出[15]，供負(fù)載或并網(wǎng)使用。

根據(jù)德國(guó)Schmidt和Sauer教授的研究，光伏逆變器的數(shù)學(xué)模型跟其物理特性密切相關(guān)[16]，如下所示：

其中，ηW—光伏逆變器轉(zhuǎn)換效率；

PAC—光伏電站交流輸出功率；

PV—逆變器功率損耗。

所有電子器件都不是理想的，都存在電阻，當(dāng)電流經(jīng)過(guò)，就會(huì)造成一定的功率損耗[17]，如下式所示：

其中，Pself—與功率無(wú)關(guān)的自身?yè)p耗；

vloss—與光伏組件直流功率線性相關(guān)的系數(shù)；

rloss—直流功率二次項(xiàng)的系數(shù)。

光伏并網(wǎng)逆變器的損耗主要發(fā)生在IGBT功率器件、濾波電感電容、接觸器等。

四、預(yù)測(cè)模型的訓(xùn)練

在科學(xué)研究中，很多時(shí)候物理量間的相互關(guān)系很難由物理理論推導(dǎo)出來(lái)。但因?yàn)樾枰?，科研人員必須知道這些數(shù)據(jù)之間的函數(shù)關(guān)系，所以就需要利用曲線擬合的方法，利用Matlab軟件得出這些物理量之間的近似函數(shù)表達(dá)式[18]。Matlab做曲線擬合的方式有兩種：內(nèi)建函數(shù)和曲線擬合工具箱。這個(gè)工具箱包括了用Matlab建立的圖形界面和文件函數(shù)。曲線擬合工具箱可以對(duì)數(shù)學(xué)模型進(jìn)行參數(shù)擬合，也可以通過(guò)其他插值方法進(jìn)行非參數(shù)擬合[19]。

本文使用德國(guó)Magdeburg大學(xué)光伏實(shí)驗(yàn)平臺(tái)在2010年采集的氣象數(shù)據(jù)和發(fā)電量數(shù)據(jù)，通過(guò)Matlab軟件的曲線擬合工具箱對(duì)光伏發(fā)電量預(yù)測(cè)模型進(jìn)行參數(shù)擬合。

1、光伏組件擬合曲線

根據(jù)前文可知，光伏組件的直流輸出功率與光照輻射強(qiáng)度和組件溫度相關(guān)。首先，只考慮光照強(qiáng)度這一唯一的影響因素。在Matlab的圖形界面里輸入光照強(qiáng)度值和光伏組件的轉(zhuǎn)換效率，按下“Fitting”按鈕，點(diǎn)擊“new fi t”窗口填寫擬合曲線的類型，然后點(diǎn)擊“Apply”，完成函數(shù)的擬合[20]。擬合得到的與太陽(yáng)輻射強(qiáng)度相關(guān)的光伏組件轉(zhuǎn)換效率如下式所示，其擬合曲線為圖3所示：

式中，ηPV(Gm) —基于光照強(qiáng)度影響下的光伏組件轉(zhuǎn)換效率。

除此之外，光伏電池組件的直流輸出功率還與組件溫度息息相關(guān)。將太陽(yáng)輻射強(qiáng)度作為橫坐標(biāo)數(shù)據(jù)，環(huán)境溫度與光伏組件溫度間的差值作為縱坐標(biāo)數(shù)據(jù)，利用擬合工具箱獲得這兩者之間的關(guān)系擬合曲線，如式（8）所示，擬合曲線如圖4所示：

2、逆變器擬合曲線

根據(jù)逆變器的轉(zhuǎn)換效率模型，將逆變器輸入端的直流功率作為Matlab擬合圖形界面里的橫坐標(biāo)數(shù)據(jù)，逆變器的轉(zhuǎn)換效率作為縱坐標(biāo)數(shù)據(jù)，使用Curve Fitting Toolbox獲得擬合曲線的參數(shù)。其擬合公式如式（9），擬合曲線如圖5所示：

綜上所述，該并網(wǎng)光伏系統(tǒng)的直流和交流功率預(yù)測(cè)擬合曲線可表示為式（10）和式（11）：

五、精確度分析

為了評(píng)價(jià)模型的精確性，將會(huì)用到兩個(gè)概念：誤差的平均值和標(biāo)準(zhǔn)差。計(jì)算公式如下：

系統(tǒng)的平均誤差能很好的描述誤差大小，而系統(tǒng)的標(biāo)準(zhǔn)差能反應(yīng)一個(gè)數(shù)據(jù)集的離散程度，平均誤差相同時(shí)，標(biāo)準(zhǔn)差未必相同。當(dāng)誤差平均值和誤差的標(biāo)準(zhǔn)差同時(shí)都相對(duì)較小時(shí)，可以判定此模型有較好的性能。

通過(guò)并網(wǎng)光伏系統(tǒng)發(fā)電量預(yù)測(cè)的擬合公式求得系統(tǒng)的交流輸出功率，將此計(jì)算的預(yù)測(cè)值與實(shí)際測(cè)量值進(jìn)行比較，就可以判斷這個(gè)并網(wǎng)光伏系統(tǒng)模型的精確度。其相對(duì)誤差平均值和標(biāo)準(zhǔn)差分別為0.0627和0.1293。其誤差的正態(tài)分布曲線如圖6所示。

由圖6可知，運(yùn)用混合方法獲得的并網(wǎng)光伏系統(tǒng)發(fā)電量預(yù)測(cè)模型具有較高的精確度。它在考慮光伏組件自身物理特性的同時(shí)，也考慮了外界環(huán)境對(duì)光伏組件性能的影響。

六、結(jié)論

總而言之，通過(guò)對(duì)并網(wǎng)光伏系統(tǒng)發(fā)電量預(yù)測(cè)的建模，可以事先預(yù)測(cè)光伏系統(tǒng)的產(chǎn)能。通過(guò)計(jì)算值與實(shí)際測(cè)量值的對(duì)比分析，可以對(duì)光伏系統(tǒng)性能進(jìn)行判定。當(dāng)計(jì)算值和測(cè)量值之間的誤差愈來(lái)愈大時(shí)，說(shuō)明此并網(wǎng)光伏系統(tǒng)的性能愈來(lái)愈差。原因有很多，其中光伏系統(tǒng)組件的老化是最可能也是最重要的影響因素。因此，對(duì)并網(wǎng)光伏系統(tǒng)進(jìn)行建模能夠?qū)ζ溥M(jìn)行長(zhǎng)期的后期監(jiān)測(cè)和維護(hù)。此外，對(duì)光伏系統(tǒng)產(chǎn)能進(jìn)行準(zhǔn)確的預(yù)測(cè)，有利于電網(wǎng)部門及時(shí)調(diào)整計(jì)劃，有利于提高太陽(yáng)能的利用效率，有利于獲得良好的經(jīng)濟(jì)、社會(huì)效益。