孫景釕，胡長(zhǎng)洪，項(xiàng)燁鋆，趙 碚，劉津源，陳夢(mèng)翔，蔡昌春

（1.國(guó)網(wǎng)浙江省電力有限公司溫州供電公司，浙江 溫州 325000；2.江蘇省輸配電裝備技術(shù)重點(diǎn)實(shí)驗(yàn)室（河海大學(xué)），江蘇 常州 213022）

0 引言

電力系統(tǒng)短期負(fù)荷預(yù)測(cè)是電力系統(tǒng)運(yùn)行與調(diào)度的基礎(chǔ)，預(yù)測(cè)精度直接影響電力系統(tǒng)運(yùn)行的安全性和經(jīng)濟(jì)性。影響負(fù)荷波動(dòng)的因素眾多，其中氣象因素尤為重要［1-2］，氣象因素直接影響用戶(hù)用電習(xí)慣從而導(dǎo)致負(fù)荷波動(dòng)。因此，從負(fù)荷波動(dòng)的機(jī)理出發(fā)，分析氣象因素對(duì)負(fù)荷變化的影響對(duì)于提高短期負(fù)荷預(yù)測(cè)精度具有重要意義。

近年來(lái)，國(guó)內(nèi)外學(xué)者針對(duì)負(fù)荷預(yù)測(cè)及影響因素開(kāi)展了大量的研究工作。文獻(xiàn)［3］針對(duì)通過(guò)可視化降維方法解決負(fù)荷預(yù)測(cè)中高維氣象數(shù)據(jù)的擁擠和結(jié)構(gòu)易變問(wèn)題，從而提高負(fù)荷預(yù)測(cè)精度。文獻(xiàn)［4-5］利用費(fèi)歇信息進(jìn)行氣象因素建模，挖掘負(fù)荷-氣象因素關(guān)聯(lián)性使得負(fù)荷預(yù)測(cè)中氣象因素的處理更加客觀。文獻(xiàn)［6-8］通過(guò)分析母線(xiàn)負(fù)荷與氣象因素相關(guān)特性，創(chuàng)建了基于數(shù)值天氣預(yù)報(bào)和負(fù)荷分類(lèi)預(yù)測(cè)的母線(xiàn)負(fù)荷預(yù)測(cè)模型，并提出了基于氣象預(yù)報(bào)的母線(xiàn)負(fù)荷預(yù)測(cè)方法，但是對(duì)于氣象因素沒(méi)有給出明確的劃分和選擇方法。

為了更加精準(zhǔn)描述氣象因素對(duì)負(fù)荷預(yù)測(cè)的影響，聚類(lèi)分析法常常被用來(lái)提取數(shù)據(jù)樣本特征以分析氣象因素的影響［9-12］。文獻(xiàn)［13］針對(duì)多類(lèi)樣本數(shù)據(jù)提出多核模糊聚類(lèi)算法，選取子核函數(shù)及其參數(shù)用于構(gòu)造多核函數(shù)，從而增大不同類(lèi)別樣本間的差別。文獻(xiàn)［14］提出了一種電力短期負(fù)荷場(chǎng)景中改進(jìn)多核模糊C 均值聚類(lèi)算法，提升預(yù)測(cè)聚類(lèi)能力，但是在核函數(shù)各參數(shù)選擇偏于經(jīng)驗(yàn)。

人工智能、深度學(xué)習(xí)方法能夠充分挖掘負(fù)荷變化的本質(zhì)特征，其在負(fù)荷預(yù)測(cè)領(lǐng)域的應(yīng)用越來(lái)越廣泛。文獻(xiàn)［15-16］提出了一種GRU-NN（門(mén)控循環(huán)單元神經(jīng)網(wǎng)絡(luò)）的短期負(fù)荷預(yù)測(cè)方法，通過(guò)建模學(xué)習(xí)負(fù)荷數(shù)據(jù)內(nèi)部動(dòng)態(tài)變化規(guī)律，融合天氣影響因素分析負(fù)荷變化的內(nèi)在聯(lián)系。文獻(xiàn)［17-18］提出了在LSTM（長(zhǎng)短期記憶）神經(jīng)網(wǎng)絡(luò)中采用Attention 機(jī)制進(jìn)行輸入特征分析，提高歷史訓(xùn)練數(shù)據(jù)特征的可辨識(shí)性。文獻(xiàn)［19-20］利用改進(jìn)深度系數(shù)自編碼器和極限學(xué)習(xí)機(jī)結(jié)合的方法進(jìn)行短期負(fù)荷預(yù)測(cè)，實(shí)現(xiàn)氣象因素的稀疏處理和低維壓縮編碼。然而，傳統(tǒng)LSTM 神經(jīng)網(wǎng)絡(luò)在負(fù)荷預(yù)測(cè)中往往存在誤差累計(jì)現(xiàn)象，從而影響預(yù)測(cè)精度。

本文針對(duì)復(fù)雜氣象因素下負(fù)荷預(yù)測(cè)困難的問(wèn)題，分別從負(fù)荷數(shù)據(jù)處理和預(yù)測(cè)方法兩方面開(kāi)展工作。提出了基于多核C 均值模糊聚類(lèi)的負(fù)荷影響因素聚類(lèi)壓縮，利用核函數(shù)提取影響負(fù)荷波動(dòng)的氣象因素動(dòng)態(tài)特征；在傳統(tǒng)LSTM結(jié)構(gòu)的基礎(chǔ)上，引入反饋環(huán)節(jié)，利用正反向計(jì)算結(jié)合消除LSTM神經(jīng)網(wǎng)絡(luò)訓(xùn)練過(guò)程的累計(jì)誤差，提出基于深度學(xué)習(xí)的多層雙向LSTM 神經(jīng)網(wǎng)絡(luò)負(fù)荷預(yù)測(cè)模型，由此提升負(fù)荷預(yù)測(cè)精度。最后，以歷史負(fù)荷的聚類(lèi)數(shù)據(jù)為訓(xùn)練樣本驗(yàn)證本文方法的有效性和合理性，實(shí)驗(yàn)結(jié)果表明本文方法能夠充分考慮歷史負(fù)荷數(shù)據(jù)中的氣象因素特征，提高負(fù)荷預(yù)測(cè)的準(zhǔn)確性。

1 基于多核模糊C 均值聚類(lèi)的負(fù)荷氣象影響分析

1.1 負(fù)荷氣象因素關(guān)系

影響負(fù)荷波動(dòng)的氣象因素較多，實(shí)際負(fù)荷數(shù)據(jù)和氣象因素間存在強(qiáng)耦合、非線(xiàn)性關(guān)系。分析負(fù)荷與氣象影響因素的相關(guān)性是剖析負(fù)荷特征的主要方式。斯皮爾曼相關(guān)系數(shù)法是分析氣象因素與負(fù)荷相關(guān)性的常用方法，可以更直觀的看出每一類(lèi)影響因素對(duì)負(fù)荷的影響程度，能夠精準(zhǔn)的刻畫(huà)變量間的關(guān)聯(lián)程度。斯皮爾曼相關(guān)系數(shù)法是通過(guò)內(nèi)部重復(fù)值來(lái)比較內(nèi)部變量間的單調(diào)性，用+1或-1來(lái)表示，斯皮爾曼相關(guān)系數(shù)公式如下所示：

1.2 多核C均值模糊聚類(lèi)

核函數(shù)法是數(shù)據(jù)特征提取的常用方法，可以彌補(bǔ)單一核方法提取負(fù)荷數(shù)據(jù)特性存在特征描述單一的不足。利用多核函數(shù)將氣象因素的低維非線(xiàn)性關(guān)系轉(zhuǎn)變?yōu)楦呔S線(xiàn)性問(wèn)題從而提高聚類(lèi)對(duì)多維空間數(shù)據(jù)的聚合能力。在聚類(lèi)過(guò)程中采用基于模糊準(zhǔn)則自動(dòng)調(diào)節(jié)不同核函數(shù)的權(quán)重系數(shù)，提升組合核函數(shù)的特征學(xué)習(xí)能力和聚類(lèi)后系統(tǒng)描述的泛化能力。

假設(shè)聚類(lèi)模型樣本集為{x1，x2，…，xN}，xk∈RN，利用非線(xiàn)性映射函數(shù)Φ將低維空間的樣本映射到高維特征空間中{Φ(x1)，Φ(x2)，…，Φ(xN)}。高維特征空間中基本核函數(shù)可以用低維空間的點(diǎn)積表示：

在滿(mǎn)足Mercer 條件下將多個(gè)核函數(shù)構(gòu)造凸線(xiàn)性組合以提高組合核函數(shù)的泛化能力，線(xiàn)性組合關(guān)系表示為：

式中：Kk為基本核函數(shù)；β為不同核函數(shù)的權(quán)重系數(shù)；M為核函數(shù)的個(gè)數(shù)。

利用核函數(shù)將負(fù)荷數(shù)據(jù)樣本中的氣象因素從高維空間向低維空間映射，并提取各種影響因素的特征空間，實(shí)現(xiàn)每類(lèi)樣本特征的劃分和分類(lèi)。文中多核C均值模糊聚類(lèi)的目標(biāo)函數(shù)為：

式中：d為第i個(gè)樣本聚類(lèi)數(shù)據(jù)到第j類(lèi)聚類(lèi)中心的Euclid距離，由式（8）計(jì)算：

聚類(lèi)中心Φ(Vi)在特征空間H中表示為：

由上述可知，多核模糊C 均值聚類(lèi)模型聚類(lèi)過(guò)程中需要明確核函數(shù)參數(shù)值σ、模糊指數(shù)m、收斂閾值ε和聚類(lèi)簇?cái)?shù)c等四個(gè)參數(shù)。通過(guò)初始化中心矩陣及隸屬度矩陣使其滿(mǎn)足歸一化條件，利用迭代計(jì)算獲得最終權(quán)重和隸屬度，迭代過(guò)程的終止條件為隸屬度矩陣||U(t)-U(t-1)||＜ε。

2 基于多層雙向LSTM 神經(jīng)網(wǎng)絡(luò)的負(fù)荷預(yù)測(cè)

2.1 雙向LSTM神經(jīng)網(wǎng)絡(luò)

傳統(tǒng)單向LSTM 神經(jīng)網(wǎng)絡(luò)模型在短期負(fù)荷預(yù)測(cè)中由于存儲(chǔ)歷史與當(dāng)前的信息，導(dǎo)致在訓(xùn)練過(guò)程中容易出現(xiàn)累積誤差等問(wèn)題。如圖1所示，雙向LSTM 神經(jīng)網(wǎng)絡(luò)模型輸出由兩個(gè)信息傳遞相反的LSTM 循環(huán)層構(gòu)成，前向?qū)影凑諘r(shí)間順序傳遞信息，后向?qū)影磿r(shí)間順序逆向傳遞信息。在神經(jīng)網(wǎng)絡(luò)訓(xùn)練過(guò)程中，分前向訓(xùn)練和后向訓(xùn)練兩個(gè)LSTM 模塊，前向訓(xùn)練輸入序列是數(shù)據(jù)樣本，后向訓(xùn)練的輸入序列的反向樣本，網(wǎng)絡(luò)的輸出由前向輸出和后向輸出共同決定。

圖1 雙向LSTM神經(jīng)網(wǎng)絡(luò)記憶單元結(jié)構(gòu)

前向?qū)訌? 時(shí)刻到t時(shí)刻正向計(jì)算得到并保存每個(gè)時(shí)刻向前隱含層的輸出；后向?qū)幽鏁r(shí)間序列反向計(jì)算得到并保存每個(gè)時(shí)刻向后隱含層的輸出。最后融合前向?qū)雍秃笙驅(qū)訉?duì)應(yīng)時(shí)刻的輸出結(jié)果，雙向LSTM神經(jīng)網(wǎng)絡(luò)的表達(dá)式為：

式中：st為t時(shí)刻的前向隱藏層狀態(tài)；為t時(shí)刻的反向隱藏層狀態(tài)；ot為t時(shí)刻輸出層的值；g和f為激活函數(shù)；xt為輸入向量；V、U、W為隱藏層到輸出層、輸入層到隱藏層、隱藏層之間的權(quán)重矩陣；V′、U′、W′為對(duì)應(yīng)的反向權(quán)重矩陣。

2.2 多層堆疊雙向LSTM神經(jīng)網(wǎng)絡(luò)

在傳統(tǒng)雙向LSTM 神經(jīng)網(wǎng)絡(luò)的基礎(chǔ)上構(gòu)建深度化多層機(jī)制和反向反饋機(jī)制，強(qiáng)化負(fù)荷數(shù)據(jù)本質(zhì)特征，避免學(xué)習(xí)過(guò)程中訓(xùn)練記憶的模糊化，減少訓(xùn)練誤差。本文提出的多層堆疊雙向LSTM 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)如圖2所示，相對(duì)單層網(wǎng)絡(luò)多層堆疊模型能夠?qū)υ缙陬A(yù)測(cè)數(shù)據(jù)加深記憶，形成正、反記憶網(wǎng)絡(luò)雙向訓(xùn)練，從負(fù)荷數(shù)據(jù)的不同時(shí)序提取數(shù)據(jù)特征。多層堆疊形態(tài)的LSTM 神經(jīng)網(wǎng)絡(luò)各層之間的關(guān)聯(lián)如下所示：

圖2 多層雙向LSTM神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

2.3 基于數(shù)據(jù)聚類(lèi)的負(fù)荷預(yù)測(cè)流程

通過(guò)對(duì)歷史負(fù)荷數(shù)據(jù)進(jìn)行基于氣象因素的聚類(lèi)分析，構(gòu)建多層雙向LSTM 神經(jīng)網(wǎng)絡(luò)的短期負(fù)荷預(yù)測(cè)模型，具體流程如圖3所示，負(fù)荷預(yù)測(cè)過(guò)程分兩個(gè)階段，分別為負(fù)荷訓(xùn)練樣本數(shù)據(jù)的聚類(lèi)分析和負(fù)荷預(yù)測(cè)。

圖3 負(fù)荷預(yù)測(cè)流程

第一階段：負(fù)荷訓(xùn)練數(shù)據(jù)的多核均值聚類(lèi)。構(gòu)建基于高斯核和線(xiàn)性核組合的多核模糊C 均值聚類(lèi)方法，實(shí)現(xiàn)多特征氣象因素樣本數(shù)據(jù)集的低維非線(xiàn)性空間映射至高維線(xiàn)性空間，通過(guò)指標(biāo)評(píng)價(jià)構(gòu)建負(fù)荷樣本訓(xùn)練數(shù)據(jù)。

第二階段：多層雙向LSTM 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)?；谟?xùn)練數(shù)據(jù)構(gòu)建短期負(fù)荷預(yù)測(cè)模型進(jìn)行負(fù)荷預(yù)測(cè)，在神經(jīng)網(wǎng)絡(luò)訓(xùn)練過(guò)程中利用預(yù)測(cè)誤差評(píng)價(jià)指標(biāo)優(yōu)化神經(jīng)網(wǎng)絡(luò)參數(shù)，評(píng)價(jià)指標(biāo)主要有平均絕對(duì)百分比誤差MAPE、平均絕對(duì)誤差MAE和均方根誤差RMSE。

式中：yi為實(shí)際測(cè)量負(fù)荷數(shù)據(jù)；預(yù)測(cè)模型輸出負(fù)荷數(shù)據(jù)；M為負(fù)荷樣本數(shù)量。

3 仿真驗(yàn)證

為了驗(yàn)證本文方法的合理性和準(zhǔn)確性，本文采用配電網(wǎng)變電站實(shí)際運(yùn)行數(shù)據(jù)和相關(guān)氣象數(shù)據(jù)進(jìn)行仿真驗(yàn)證。負(fù)荷訓(xùn)練樣本數(shù)據(jù)的采樣精度為5 min和15 min兩個(gè)時(shí)間尺度，負(fù)荷訓(xùn)練數(shù)據(jù)為一個(gè)月的數(shù)據(jù)量，預(yù)測(cè)未來(lái)一天和星期進(jìn)行分析。在進(jìn)行氣象因素聚類(lèi)時(shí)，聚類(lèi)核函數(shù)參數(shù)σ=0.005、隸屬度函數(shù)模糊指數(shù)m=1.08、收斂閾值ε=0.000 1，聚類(lèi)迭代次數(shù)為300 次、聚類(lèi)簇為c=3。本文采用python 開(kāi)發(fā)LSTM 深度神經(jīng)網(wǎng)絡(luò)，電腦配置為I9-9900k-3.6GHz，16G內(nèi)存。

3.1 氣象因素聚類(lèi)分析

利用多核模糊C 均值聚類(lèi)模型分析負(fù)荷訓(xùn)練數(shù)據(jù)中的氣象影響因素，利用多核模糊C 均值聚類(lèi)模型對(duì)原始數(shù)據(jù)集樣本進(jìn)行聚類(lèi)處理。樣本數(shù)據(jù)中影響因素為平均溫度等8種氣象數(shù)據(jù)。通過(guò)多核模糊C 均值算法聚類(lèi)后的樣本數(shù)據(jù)分布結(jié)果如圖4所示，從圖4中可以看出，聚類(lèi)模型將多維氣象因素聚類(lèi)為兩類(lèi)，兩類(lèi)子簇之間有邊界清晰，其聚類(lèi)結(jié)果的誤差不超過(guò)4%?；谏鲜?中氣象數(shù)據(jù)的負(fù)荷相關(guān)性分析如表1所示，聚類(lèi)結(jié)果的各項(xiàng)指標(biāo)如表2 所示，多核模糊C 均值聚類(lèi)較FCM（模糊C 均值）、KFCM（單核模糊C 均值）和DBSCAN（密度聚類(lèi)算法）在NMI（標(biāo)準(zhǔn)化互信息）、ACC（準(zhǔn)確度）、RI（蘭德系數(shù)）和ARI（調(diào)整蘭德系數(shù)）等4 個(gè)指標(biāo)都具有更加合理的聚類(lèi)效果。

圖4 數(shù)據(jù)集1三維分布

表1 不同影響因素的斯皮爾曼相關(guān)系數(shù)

表2 不同聚類(lèi)算法比較

3.2 神經(jīng)網(wǎng)絡(luò)層數(shù)選擇

利用多層堆疊神經(jīng)網(wǎng)絡(luò)進(jìn)行負(fù)荷預(yù)測(cè)，此時(shí)不同層數(shù)的神經(jīng)網(wǎng)絡(luò)在不同程度上提取數(shù)據(jù)特征，神經(jīng)網(wǎng)絡(luò)層數(shù)和最終預(yù)測(cè)精度存在一定的關(guān)系。表3給出了神經(jīng)網(wǎng)絡(luò)模型層數(shù)與預(yù)測(cè)誤差之間的關(guān)系，本文神經(jīng)網(wǎng)絡(luò)的層數(shù)為3層。圖5給出了不同預(yù)測(cè)方法的結(jié)果對(duì)比，由圖5 可知多層雙向LSTM 模型的預(yù)測(cè)結(jié)果相較于傳統(tǒng)LSTM 模型、雙向LSTM神經(jīng)網(wǎng)絡(luò)以及BP網(wǎng)絡(luò)，預(yù)測(cè)結(jié)果更加接近負(fù)荷真實(shí)值。

表3 多層堆疊雙向LSTM層數(shù)與預(yù)測(cè)誤差

圖5 不同神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)的對(duì)比

3.3 基于聚類(lèi)分析的負(fù)荷預(yù)測(cè)

通過(guò)聚類(lèi)分析將影響因素根據(jù)其對(duì)負(fù)荷的影響進(jìn)行聚類(lèi)劃分，進(jìn)一步的對(duì)比不同算法的負(fù)荷預(yù)測(cè)結(jié)果。經(jīng)過(guò)聚類(lèi)分析后，將訓(xùn)練數(shù)據(jù)的影響因素分為a和b兩類(lèi)，分別進(jìn)行基于MBLSTM（多層雙向長(zhǎng)短期記憶）神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)，結(jié)果如圖6 所示，圖7 給出了不同LSTM 神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)結(jié)果對(duì)比，預(yù)測(cè)誤差的定量分析如表4所示。由表4可知在考慮氣象因素并對(duì)其進(jìn)行聚類(lèi)分析后，誤差從1.44%降低至0.25%。通過(guò)氣象因素聚類(lèi)，能夠有效甄別影響因素相似的負(fù)荷數(shù)據(jù)，由此提高訓(xùn)練數(shù)據(jù)的相似性，但是并沒(méi)有增加計(jì)算負(fù)擔(dān)，考慮氣象因素比不考慮氣象因素的計(jì)算時(shí)間增加了0.96 s。

表4 負(fù)荷預(yù)測(cè)誤差對(duì)比

圖6 a、b類(lèi)數(shù)據(jù)預(yù)測(cè)結(jié)果

圖7 考慮氣象因素聚類(lèi)前后負(fù)荷預(yù)測(cè)

為進(jìn)一步驗(yàn)證本文方法的合理性，利用冬季負(fù)荷數(shù)據(jù)進(jìn)行仿真分析，此時(shí)負(fù)荷預(yù)測(cè)的結(jié)果如圖8所示。利用聚類(lèi)劃分處理訓(xùn)練樣本數(shù)據(jù)并以此進(jìn)行預(yù)測(cè)大大降低了氣象因素對(duì)負(fù)荷訓(xùn)練、預(yù)測(cè)過(guò)程的影響。從表5可以看出，未考慮氣象因素的短期負(fù)荷預(yù)測(cè)模型的MAPE 誤差值在4.15%，由于冬季負(fù)荷的敏感性較低，整體上冬季負(fù)荷預(yù)測(cè)的誤差較夏季大，考慮氣象因素分類(lèi)整體預(yù)測(cè)誤差為0.65%。

圖8 考慮聚類(lèi)氣象因素負(fù)荷預(yù)測(cè)結(jié)果

表5 負(fù)荷預(yù)測(cè)誤差對(duì)比

4 結(jié)語(yǔ)

本文提出了基于多層模糊C 均值聚類(lèi)和MBLSTM神經(jīng)網(wǎng)絡(luò)的短期負(fù)荷預(yù)測(cè)方法。通過(guò)多層模糊C 均值聚類(lèi)將負(fù)荷氣象影響因素進(jìn)行低維空間映射，實(shí)現(xiàn)影響因素的精確分類(lèi)和聚合，明確負(fù)荷氣象因素的耦合關(guān)系；進(jìn)一步的，提出多層堆疊形式的雙向LSTM 神經(jīng)網(wǎng)絡(luò)，通過(guò)構(gòu)建雙向計(jì)算機(jī)制和多層堆疊形式，提高負(fù)荷數(shù)據(jù)間的前后關(guān)聯(lián)，挖掘負(fù)荷序列數(shù)據(jù)關(guān)聯(lián)信息，從而提高負(fù)荷預(yù)測(cè)精度。仿真結(jié)果表明，本文所提方法能夠提高負(fù)荷影響因素的區(qū)分度，有效提高負(fù)荷預(yù)測(cè)精度。