朱瑞金，龔雪嬌，張娟娟

(西藏農(nóng)牧學(xué)院電氣工程學(xué)院，西藏 林芝 860000)

0 引 言

光伏發(fā)電作為可再生清潔能源在電力系統(tǒng)中占比越來越大，由于光伏發(fā)電功率受到太陽輻射角、溫度、濕度等不確定因素的影響，造成輸出功率隨機性強，給光伏并網(wǎng)帶來了調(diào)度困難、安全性低等問題[1-2]。為保障電網(wǎng)安全穩(wěn)定運行，現(xiàn)主要采用拉閘限電等維護策略[3]，存在著耗能大，不利于光伏發(fā)電的繼續(xù)發(fā)展。近年來，國內(nèi)外研究者為實現(xiàn)常規(guī)能源與光伏發(fā)電之間的合理調(diào)度，提出了一系列對光伏發(fā)電功率提前預(yù)測算法[4-6]，以適應(yīng)其輸出功率對電網(wǎng)系統(tǒng)的影響。

光伏發(fā)電功率預(yù)測主要分為間接預(yù)測和直接預(yù)測兩種方式。間接預(yù)測方法是根據(jù)光伏發(fā)電輸入條件構(gòu)造信息特征來進行預(yù)測，如光照強度、太陽輻射強度、風(fēng)速、溫度、天氣狀態(tài)等外部條件。徐先峰等[7]提出利用雙向長短期記憶(BiLSTM)算法，充分挖掘氣象數(shù)據(jù)以實現(xiàn)光伏發(fā)電功率的短期預(yù)測，將光伏發(fā)電歷史功率曲線下的氣象條件與當(dāng)前氣象數(shù)據(jù)進行曲線擬合時，考慮到氣象條件的隨機性，存在可靠性和信號特性相關(guān)性低等問題。王志忠等[8]根據(jù)晴朗系數(shù)構(gòu)建基于相似日輸出功率的ARIMA 預(yù)測模型，以單純的晴朗系數(shù)作為輸入特性時，存在信號特征單一和樣本集要求高等問題。

直接預(yù)測方式是將歷史發(fā)電輸出功率作為研究對象建立功率預(yù)測模型。候松寶等[9]提出利用主成分分析和BP 神經(jīng)網(wǎng)絡(luò)進行短期光伏發(fā)電功率短期預(yù)測，由于受天氣狀況影響，當(dāng)外部條件變化時，功率預(yù)測精度和穩(wěn)定性將會下降。李攀等[10]利用稀疏貝葉斯學(xué)習(xí)算法，對輸出功率的影響因素進行權(quán)重分析，建立觀測器模型對輸出功率進行回歸分析。

基于以上分析，本文提出將間接預(yù)測法和直接預(yù)測法相結(jié)合的EEMD-MPE-LSSVM 光伏發(fā)電功率預(yù)測算法，使輸入條件參與輸出功率預(yù)測，提高了預(yù)測的準(zhǔn)確度。采用直接預(yù)測方法以光伏發(fā)電輸出功率為研究對象，EEMD-LSSVM 預(yù)測方法結(jié)合迭代誤差完成功率預(yù)測；考慮到間接預(yù)測方法中天氣因素的影響，將不同天氣類型進行信號復(fù)雜度量化分析，采用MPE 構(gòu)建輸入特征向量，并參與到LSSVM 曲線預(yù)測中。通過現(xiàn)場數(shù)據(jù)進行實例驗證，該算法提高了預(yù)測精度，具有一定適用性。

1 光伏發(fā)電系統(tǒng)輸出功率特性

光伏發(fā)電系統(tǒng)輸出功率具有非線性、不平穩(wěn)特性，影響因素除天氣類型以外，還考慮到氣象條件使得輸出功率隨機性增加。本文參照中國氣象局制定的GB/T 22164—2008 對33種天氣類型進行模糊分類，將相似度較高的天氣類型歸為一類，經(jīng)分類后的天氣主要有晴天、雨雪、多云和突變天氣4種類型。

在4種不同天氣類型下進行輸出功率統(tǒng)計，圖1為某光伏發(fā)電站從06:00-20:00，每15 min 采集一次光伏發(fā)電輸出功率的曲線，其縱坐標(biāo)為光伏發(fā)電輸出功率。

圖1 不同天氣類型下光伏發(fā)電輸出功率曲線

從圖1 中可以看出：陰天下受到太陽輻射最小，輸出功率也隨之減少；晴天下的發(fā)電功率穩(wěn)定；遇到突變天氣，發(fā)電功率呈現(xiàn)不穩(wěn)定現(xiàn)象；雨雪天氣下較晴天發(fā)電功率下降。在連續(xù)時間序列下，用排列熵對不同天氣類型進行量化，并作為SSVM 的輸入向量的一部分提高預(yù)測精度。

2 功率預(yù)測算法分析

2.1 集成經(jīng)驗?zāi)B(tài)分解(EEMD)

EEMD 分解方式是在經(jīng)驗?zāi)B(tài)分解(EMD)基礎(chǔ)上提出的，EMD 算法對信號序列進行分解時，每個模態(tài)函數(shù)(instrinsic mode function,IFM)必須滿足兩個條件[11]：

1)極值點和過零點的個數(shù)誤差不超過1。

2)上下包絡(luò)線的均值接近0。

為充分表征原始信號序列的頻率特性，當(dāng)剩余分量rn(t)為單調(diào)信號或沒有足夠極值點時結(jié)束迭代。以x(t)為原始信號序列，EMD 分解如下式所示：

考慮到當(dāng)信號序列極值點不足時，EMD 分解存在模態(tài)混疊的問題，即一種模態(tài)分量可能存在不同頻率特性。為此提出EEMD 分解方式，在分解之前加入n組均值為零的高斯白噪聲hi(t)，如下式所示：

再對每組信號進行EMD 分解：

式中：Cij(t)——第i組模態(tài)分解；

ri(t)——第i組剩余分量；

j——分解層數(shù)。

最后根據(jù)高斯白噪聲均值為零的特性，對i組模態(tài)分量進行求均值處理，即：

EEMD 分解中Cj(t)代替式(1)中IFMi(t)模態(tài)分量，成為新的模態(tài)分量。

2.2 多尺度排列熵(MPE)

為表征不同信號序列的復(fù)雜度，引入MPE 算法計算熵值，從而實現(xiàn)復(fù)雜度的量化目標(biāo)。在熵值計算之前，需要對原始信號序列X={x(i)|i=1,2,···,n}進行相空間重構(gòu)，假設(shè)時延參數(shù)等于τ，嵌入維度等于m，則重構(gòu)矩陣為：

其中j=1,2,···,K；K=n–(m–1)τ

對重構(gòu)矩陣中每一行進行降序或升序排列，計算每種排序方式在m！種排序中的概率分布P1，P2,···,PK，按照下式計算其多尺度排列熵，并進行歸一化處理：

由熵值定義可知：熵值越大，表示信號序列越復(fù)雜。

2.3 最小二乘支持向量機(LSSVM)

支持向量機(SVM)算法引用最小二乘損失函數(shù)，將最優(yōu)問題轉(zhuǎn)換為線性方程，從而實現(xiàn)樣本數(shù)據(jù)的預(yù)測[12]。具體算法步驟如下所示：

1)構(gòu)造目標(biāo)函數(shù)：對于線性不可分問題，加入核函數(shù)φ (xi)和懲罰因子c。

式中：y——故障類別標(biāo)簽，即yi∈[–1,1]；

ξ——誤差分量。

2)利用拉格朗日函數(shù)，對目標(biāo)函數(shù)進行求解，當(dāng)拉格朗日乘子等于αi時，優(yōu)化問題表述為：

3)利用KKT 條件對式(9)進行求解，最后求得回歸函數(shù)為：

LSSVM 具有小樣本、自適應(yīng)預(yù)測的優(yōu)勢，采用徑向基函數(shù)進行線性化處理，減少了輸入?yún)?shù)，提高了運算速度。

3 光伏發(fā)電功率預(yù)測流程

本文采用EEMD-MPE-LSSVM 模型進行光伏發(fā)電功率預(yù)測。在實際應(yīng)用中受到氣象條件的影響，首先用排列熵表征對天氣類型進行量化；然后利用EEMD 分解方式提取細節(jié)分量，并對分量回歸預(yù)測后進行重構(gòu)；最后計算初始訓(xùn)練樣本迭代誤差序列，從而實現(xiàn)發(fā)電功率的整體預(yù)測。具體算法流程如圖2 所示。

圖2 EEMD-MPE-LSSVM 光伏發(fā)電功率預(yù)測流程圖

1）考慮到同一天有可能存在不同天氣狀況，選擇4 個天氣類型相似的臨近日發(fā)電功率計算排列熵，并作為LSSVM 輸入特征向量。

2）當(dāng)構(gòu)造初始迭代誤差序列時，在歷史數(shù)據(jù)中選擇q個數(shù)據(jù)作為訓(xùn)練樣本集，以第a個數(shù)據(jù)點為起點，利用滾動時間窗口預(yù)測方法得到q–a個功率預(yù)測值。在與真實值比較后得到預(yù)測的一級誤差，并重復(fù)上述步驟，最終獲取x(e1)～x(ep)為1 級～p級誤差序列。在實際應(yīng)用中參考文獻[13]，p一般不超過3。

3）EEMD 分解結(jié)合排列熵特征對每一模態(tài)函數(shù)IFM 進行LSSVM 預(yù)測，將重構(gòu)后的預(yù)測結(jié)果加入迭代誤差序列得到最終預(yù)測值。

4 實驗部分

4.1 不同天氣類型的排列熵計算

本文選用某裝機容量為50 MW 光伏發(fā)電站實驗數(shù)據(jù)來源，隨意抽取4種不同天氣類型，并通過臨近日發(fā)電功率計算排列熵如表1 所示。

表1 不同天氣類型下的排列熵

4.2 初始訓(xùn)練樣本功率序列EEMD 分解

光伏發(fā)電輸出功率受到多個不確定因素的影響，以發(fā)電系統(tǒng)前20 天歷史數(shù)據(jù)為訓(xùn)練集，畫出相應(yīng)發(fā)電功率曲線圖。由于輸出功率呈現(xiàn)非線性、不平穩(wěn)特征，為充分表征曲線特性，采用EEMD 分解方式將原始曲線分解成具有不同頻率特性的模態(tài)函數(shù)分量IFM1～IFM4，圖3 為20 天內(nèi)光伏發(fā)電功率曲線及EEMD 分解結(jié)果。

由圖3 可知，IFM4 符合沒有足夠極值點的EEMD分解終止條件。每個模態(tài)函數(shù)分量在保持功率曲線特性同時，降低了LSSVM 計算復(fù)雜度，提高了計算精度。

圖3 訓(xùn)練集樣本發(fā)電功率EEMD 分解結(jié)果

4.3 EEMD-MPE-LSSVM 功率預(yù)測結(jié)果

光伏發(fā)電功率受到太陽輻射角、氣象條件、風(fēng)速、濕度、溫度等多個條件的影響，因此為預(yù)測特定時刻t的發(fā)電功率，以模態(tài)分量結(jié)合不同天氣類型排列熵作為預(yù)測特征集進行樣本訓(xùn)練。分別將預(yù)測之后的IFM1、IFM2、IFM3、IFM4 進行重構(gòu)，并記為y1。然后將初始迭代誤差x(e1)、x(e2)、x(e3)加入初始預(yù)測值y1。最后，最終t時刻的光伏發(fā)電功率預(yù)測值為y=y1+e1+e2+e3。

為驗證該預(yù)測模型的有效性，本文選擇春、夏、秋、冬4 個季節(jié)的4種不同天氣類型共計16種發(fā)電功率進行預(yù)測。以春季晴天為例，分別在EEMDLSSVM 和EEMD-MPE-LSSVM 兩種算法下進行光伏發(fā)電預(yù)測，經(jīng)樣本訓(xùn)練后的預(yù)測值與真實值分布如圖4 所示。

通過對圖4 中(a)和(b)的對比，EEMD-MPELSSVM 功率預(yù)測模型在相似日的基礎(chǔ)上，量化了氣象條件，并作為有效特征加入預(yù)測，提高了預(yù)測模型的精度。本文統(tǒng)計了對光伏發(fā)電功率50 天內(nèi)的預(yù)測值和有效值，分別在EEMD-LSSVM 和EEMDMPE-LSSVM 兩種不同的算法下評估其平均相對誤差(MRE)和均方根誤差(RMES)，結(jié)果如表2 所示。從表中可知在用MPE 量化天氣因素后，其預(yù)測結(jié)果的MRE 和RMES 分別為1.56%、3.14%，均低于EEMD-LSSVM 預(yù)測誤差。

表2 50 天內(nèi)預(yù)測誤差

圖4 春季晴天光伏發(fā)電功率預(yù)測結(jié)果

考慮到天氣狀況、太陽輻射角、溫度、濕度、氣象條件等不確定因素的影響，利用多尺度排列熵量化天氣類型，并作為輸入特征向量細化了LSSVM預(yù)測，證明了該算法的有效性。

5 結(jié)束語

本文針對光伏發(fā)電功率曲線特性，提出了結(jié)合天氣因素的EEMD-MPE-LSSVM 光伏發(fā)電預(yù)測算法，通過實驗論證得出以下結(jié)論：

1）光伏發(fā)電功率曲線數(shù)據(jù)進行EEMD 分解，通過提取出曲線細節(jié)分量，充分表征了功率曲線的局部特征，使得LSSVM 預(yù)測算法更加精準(zhǔn)。

2）MPE 用于衡量天氣因素的復(fù)雜度，同時參與到光伏發(fā)電功率預(yù)測中，進一步提高預(yù)測精度，并提供有效的細節(jié)分量。

3）光伏發(fā)電功率曲線具有隨機性、非線性、不平穩(wěn)等特性，通過迭代誤差修正發(fā)電功率預(yù)測值。

本文提出的光伏發(fā)電預(yù)測模型兼顧不同天氣類型的影響，經(jīng)小樣本訓(xùn)練后，能夠?qū)夥l(fā)電功率實現(xiàn)準(zhǔn)確預(yù)測。為實現(xiàn)光伏并網(wǎng)的合理調(diào)度和安全運行，提供了可行性參考方案。