張秀釗，王志敏，錢(qián)紋，胡凱

（云南電網(wǎng)有限責(zé)任公司電網(wǎng)規(guī)劃建設(shè)研究中心，昆明 650011）

0 前言

2015 年在《中共中央國(guó)務(wù)院關(guān)于進(jìn)一步深化電力體制改革的若干意見(jiàn)》中[1]，明確了新一輪電改，將破除電網(wǎng)企業(yè)“獨(dú)買(mǎi)獨(dú)賣(mài)”的模式，電力改革的新模式為電動(dòng)汽車(chē)的發(fā)展帶來(lái)了新的挑戰(zhàn)和機(jī)遇，由于電動(dòng)汽車(chē)用戶(hù)需求和行為的不確定性與相互差異，未來(lái)大規(guī)模電動(dòng)汽車(chē)充電負(fù)荷具有時(shí)間和空間上的隨機(jī)性、間歇性和波動(dòng)性等不確定特點(diǎn)，將給電網(wǎng)的安全運(yùn)行和優(yōu)化調(diào)度帶來(lái)困難。為了分析電動(dòng)汽車(chē)充電負(fù)荷帶來(lái)的影響，以及為電動(dòng)汽車(chē)廣泛接入電網(wǎng)的調(diào)控策略制訂打下基礎(chǔ)，提供理論支持，要求建立有效的電動(dòng)汽車(chē)充電負(fù)荷預(yù)測(cè)模型[2]。

電動(dòng)汽車(chē)負(fù)荷預(yù)測(cè)主要分為兩類(lèi)，一類(lèi)是采用數(shù)學(xué)模型模擬電動(dòng)汽車(chē)充電行為，從而得出電動(dòng)汽車(chē)負(fù)荷預(yù)測(cè)值的方法[3-8]，此類(lèi)方法在綜合考慮充電負(fù)荷的時(shí)空特性時(shí)數(shù)學(xué)模型太過(guò)復(fù)雜，難以保證預(yù)測(cè)精度。另一類(lèi)是基于歷史數(shù)據(jù)采用統(tǒng)計(jì)學(xué)習(xí)模型進(jìn)行預(yù)測(cè)的方法，用模型學(xué)習(xí)歷史數(shù)據(jù)的潛在規(guī)律，從而達(dá)到較好的預(yù)測(cè)效果。

電動(dòng)汽車(chē)充電負(fù)荷預(yù)測(cè)的傳統(tǒng)負(fù)荷預(yù)測(cè)方法有回歸分析法、相似日法等；現(xiàn)代預(yù)測(cè)方法有基于神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)法、基于小波分析的預(yù)測(cè)法以及支持向量機(jī)（Support Vector Machine，SVM）預(yù)測(cè)法等。文獻(xiàn)[9]提出一種基于大數(shù)據(jù)的電動(dòng)汽車(chē)充電負(fù)荷預(yù)測(cè)方法；文獻(xiàn)[10]提到一種基于時(shí)間序列距離測(cè)量方法的短期充電樁負(fù)荷預(yù)測(cè)模型，該方法根據(jù)負(fù)荷序列的和時(shí)間序列的距離測(cè)量來(lái)進(jìn)行充電負(fù)荷的預(yù)測(cè)；文獻(xiàn)[11-12]采用支持向量機(jī)對(duì)電動(dòng)汽車(chē)進(jìn)行短期負(fù)荷預(yù)測(cè)，并將最后結(jié)果與蒙特卡洛方法進(jìn)行比較，該方法結(jié)果明顯優(yōu)于蒙特卡洛方法；文獻(xiàn)[13]提出一種快速充電站的規(guī)劃建設(shè)方法，首先采用共享最近鄰聚類(lèi)法(SNN) 來(lái)確定充電站位置以及充電站覆蓋范圍，然后采用排隊(duì)論預(yù)測(cè)充電裝容量進(jìn)行定容；文獻(xiàn)[14]提 出Modified Pattern Sequence Forecasting (MPSF) 和 Time Weighted Dot Product Nearest Neighbor (TWDP NN)兩個(gè)算法在智能手機(jī)應(yīng)用程序上進(jìn)行電動(dòng)汽車(chē)充電負(fù)荷預(yù)測(cè)，其中MPSF算法用于預(yù)測(cè)，TWDPNN 用于使得預(yù)測(cè)算法加速，該方法能很好的在手機(jī)上對(duì)電動(dòng)汽車(chē)負(fù)荷進(jìn)行有效的預(yù)測(cè)；文獻(xiàn)[15]用時(shí)間序列模型對(duì)電動(dòng)汽車(chē)負(fù)荷進(jìn)行了預(yù)測(cè)，并采用了伯利克分校的實(shí)際充電負(fù)荷數(shù)據(jù)驗(yàn)證了模型的有效性；文獻(xiàn)[16]提出了一種基于數(shù)據(jù)新鮮度和交叉熵的組合預(yù)測(cè)模型，并驗(yàn)證了該方法比單一模型的有效性。

過(guò)去的電動(dòng)汽車(chē)充電負(fù)荷預(yù)測(cè)統(tǒng)計(jì)學(xué)習(xí)方法更多都是只考慮時(shí)間維度的預(yù)測(cè)方法，沒(méi)有考慮電動(dòng)汽車(chē)充電負(fù)荷所包含的復(fù)雜空間性。因此，綜合考慮負(fù)荷的時(shí)間及空間雙重動(dòng)態(tài)變化，才能更好的進(jìn)行時(shí)空動(dòng)態(tài)預(yù)測(cè)。深度學(xué)習(xí)（Deep Learning, DL）作為機(jī)器學(xué)習(xí)領(lǐng)域一個(gè)重要的研究熱點(diǎn)，已經(jīng)在圖像分析、語(yǔ)音識(shí)別、自然語(yǔ)言處理、視頻分類(lèi)等領(lǐng)域取得了令人矚目的成功。在預(yù)測(cè)上也取得了一定成功，文獻(xiàn)[17] 提出了一種深度時(shí)空殘差網(wǎng)絡(luò)用于預(yù)測(cè)北京市的空間人群流動(dòng)情況取得了很好的效果。

綜上所述，提出一種深度多步時(shí)空神經(jīng)網(wǎng)絡(luò)DMSTN (Deep Multi-Step Spatio-Temporal Neural Network)，該網(wǎng)絡(luò)可以很好的學(xué)習(xí)到充電樁負(fù)荷的時(shí)空動(dòng)態(tài)規(guī)律，從而從空間和時(shí)間上整體的對(duì)負(fù)荷進(jìn)行預(yù)測(cè)。

1 電動(dòng)汽車(chē)時(shí)空動(dòng)態(tài)負(fù)荷矩陣

電動(dòng)汽車(chē)負(fù)荷具有時(shí)間和空間上的隨機(jī)性，為了更好的預(yù)測(cè)這種時(shí)空動(dòng)態(tài)性，需要將充電樁上的負(fù)荷進(jìn)行時(shí)空維度的刻畫(huà)。根據(jù)充電樁位置將充電樁充電負(fù)荷用二維矩陣表示，并整理成時(shí)長(zhǎng)T 的時(shí)空序列D{D1,D2,...DT},D ∈RT×X×Y，Dt為時(shí)間t 的電動(dòng)汽車(chē)空間負(fù)荷矩陣見(jiàn)式(1)，其中是坐標(biāo)為(x,y)點(diǎn)的負(fù)荷量。

圖1 充電樁分布圖

根據(jù)10 個(gè)充電樁的經(jīng)緯度分布（見(jiàn)圖1）建立一個(gè)負(fù)荷矩陣，構(gòu)建二維負(fù)荷矩陣的方法為：

1）構(gòu)建坐標(biāo)軸，確定所有充電樁的坐標(biāo)；

2）計(jì)算出每個(gè)充電樁負(fù)荷覆蓋的范圍，每個(gè)充電樁的覆蓋范圍是自己坐標(biāo)為中心的一個(gè)正方（正方形邊長(zhǎng)由充電樁覆蓋范圍確定）；

3）將所有充電樁負(fù)荷覆蓋范圍內(nèi)填上該時(shí)刻充電樁的負(fù)荷量并累加，得到該時(shí)刻的負(fù)荷矩陣。

2 建立時(shí)空動(dòng)態(tài)負(fù)荷預(yù)測(cè)模型

為了預(yù)測(cè)未來(lái)K 個(gè)時(shí)間點(diǎn)的充電樁時(shí)空動(dòng)態(tài)負(fù)荷矩陣，需要建立一個(gè)根據(jù)過(guò)去S 個(gè)時(shí)間點(diǎn)觀(guān)測(cè)值，來(lái)預(yù)測(cè)未來(lái)K 個(gè)時(shí)間點(diǎn)的因果系統(tǒng)見(jiàn)式(2)，其中P(..) 代表一個(gè)因果系統(tǒng)，Dt代表第t 個(gè)時(shí)刻的預(yù)測(cè)負(fù)荷矩陣。

實(shí)驗(yàn)將STN (SPATIO-TEMPORAL NET WORK)[18]的結(jié)構(gòu)修改為多步時(shí)空動(dòng)態(tài)神經(jīng)網(wǎng)絡(luò)后用于構(gòu)建該因果系統(tǒng)。

2.1 深度多步時(shí)空動(dòng)態(tài)神經(jīng)網(wǎng)（DMSTN）

2.1.1 ConvLSTM層

LSTM 通常用于解決一維時(shí)間序列的預(yù)測(cè)問(wèn)題，無(wú)法考慮空間上的相關(guān)性。在此基礎(chǔ)之上文獻(xiàn)[18]提出了可以用于考慮時(shí)空兩個(gè)維度的網(wǎng)絡(luò)結(jié)構(gòu)ConvLSTM，ConvLSTM 結(jié)構(gòu)能夠?qū)W習(xí)到長(zhǎng)期的二維數(shù)據(jù)規(guī)律，這非常適合用于電動(dòng)汽車(chē)動(dòng)態(tài)負(fù)荷矩陣的預(yù)測(cè)。ConvLSTM 和LSTM 一樣，將上一層的輸出和新的輸入作為下一層的輸入。不同點(diǎn)在于加上了卷積操作之后，為不僅能夠得到時(shí)序關(guān)系，還能夠像卷積層一樣提取空間特征。這樣就得到時(shí)空兩個(gè)維度的特征，并且將狀態(tài)與狀態(tài)之間的切換也換成了卷積計(jì)算。針對(duì)有一個(gè)二維充電時(shí)空動(dòng)態(tài)負(fù)荷矩陣的時(shí)間序列的輸入XT+{X1,X2,...,XT}，其中Xt表示t 時(shí)刻的負(fù)荷矩陣。ConvLSTM 的具體過(guò)程見(jiàn)式(3)。

其中Wxf、bf表示遺忘門(mén)權(quán)重和偏置；Wxi、bi表示輸入門(mén)的權(quán)重和偏置；Wxc、bc表示更新值的權(quán)重和偏置；Wxo、bo表示更新值的權(quán)重和偏置；σ(g) 代表sigmiod 激活函數(shù)、tanh(g) 代表雙曲正切激活函數(shù)，“*”表示卷積計(jì)算，“⊙”代表Hadamard 乘法。

2.1.2 三維卷積層(3D-ConvNet)

三維卷積既可以捕獲空間信息也可以捕獲時(shí)間維度的信息。三維卷積核是將二維卷積核的擴(kuò)充為三維，將數(shù)據(jù)根據(jù)時(shí)間維度堆疊為三維數(shù)據(jù)后進(jìn)行卷積操作。針對(duì)有一個(gè)二維充電時(shí)空動(dòng)態(tài)負(fù)荷矩陣的時(shí)間序列的輸入XN={X1,X2,...,XN}，其中XN表示n 時(shí)刻的負(fù)荷矩陣，公式(4)給出了3D-ConvNet 層的與輸出。

有M 個(gè)卷積核的3D-ConvNet 層輸出由H1,...HM構(gòu)成，act( ●)表示激活函數(shù)。三維卷積可以很好的學(xué)習(xí)時(shí)空特性，尤其是短時(shí)的二維數(shù)據(jù)規(guī)律，這可以很好的學(xué)習(xí)到由電動(dòng)汽車(chē)用戶(hù)隨機(jī)充電行為引起的波動(dòng)。

2.1.3 融合層

在實(shí)際的電動(dòng)汽車(chē)時(shí)空動(dòng)態(tài)負(fù)荷中，既有長(zhǎng)期的充電負(fù)荷規(guī)律，也有短期由于天氣、事件以及個(gè)人充電行為等隨機(jī)性引起的短期充電負(fù)荷規(guī)律。為了同時(shí)學(xué)習(xí)到這兩種規(guī)律，需要將三維卷積層的輸出結(jié)果和ConvLSTM 層輸出結(jié)果結(jié)合起來(lái)，因此建立一種融合層，其表達(dá)式見(jiàn)(5)。

其 中hC、hL分 別 是3D-ConvNet 層 和ConvLSTM 層的輸出，表示網(wǎng)絡(luò)的權(quán)重和偏置。

2.2 時(shí)空動(dòng)態(tài)神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

為了消除滾動(dòng)預(yù)測(cè)帶來(lái)的誤差， 將STN(SPATIO-TEMPORAL NETWORK) 最后的輸出部分改為3D-ConvNet，使得該網(wǎng)絡(luò)能同時(shí)輸出多個(gè)時(shí)間步數(shù)的時(shí)空動(dòng)態(tài)負(fù)荷矩陣，因?yàn)榫W(wǎng)絡(luò)是同時(shí)預(yù)測(cè)的多步矩陣序列是短時(shí)間序列，所以3D-ConvNet 層作為輸出層可以更好的捕捉此類(lèi)規(guī)律。假設(shè)輸入為Ds={Dt-s,Dt-s+1,...Dt}，模型的用M 表示，該模型的所有參數(shù)用Θ 表示，則模型的表示見(jiàn)式(6)。

圖2 多步時(shí)空神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

模型的整個(gè)網(wǎng)絡(luò)結(jié)構(gòu)見(jiàn)圖2，采用將ConvLSTM 層和3D-ConvNet 層分別輸出后用融合層進(jìn)行融合的結(jié)構(gòu)單元，能夠全面的學(xué)習(xí)到長(zhǎng)期規(guī)律和短期規(guī)律，多個(gè)單元的堆疊使網(wǎng)絡(luò)的學(xué)習(xí)能力更強(qiáng)，3D-ConvNet 層作為輸出實(shí)現(xiàn)多步預(yù)測(cè)。

3 時(shí)空動(dòng)態(tài)負(fù)荷預(yù)測(cè)

3.1 算法流程

為了得到準(zhǔn)確的負(fù)荷矩陣預(yù)測(cè)結(jié)果，需求解出該模型的最優(yōu)參數(shù)，將該預(yù)測(cè)模型的目標(biāo)函數(shù)設(shè)為(7),然后采用adam 自適應(yīng)隨機(jī)梯度下降算法求解該模型的最優(yōu)參數(shù)。

圖3 算法流程圖

模型的整個(gè)流程見(jiàn)圖3，首先根據(jù)充電樁經(jīng)緯度將建立充電樁平面分布圖；然后根據(jù)歷史負(fù)荷數(shù)據(jù)，按時(shí)間順序在分布圖上做出時(shí)空動(dòng)態(tài)負(fù)荷矩陣；再將每個(gè)矩陣數(shù)據(jù)進(jìn)行歸一化；接著將數(shù)據(jù)輸入進(jìn)模型進(jìn)行訓(xùn)練，根據(jù)訓(xùn)練結(jié)果調(diào)整模型的超參數(shù)，直到模型訓(xùn)練達(dá)到滿(mǎn)意的效果；最后將圖片數(shù)據(jù)反歸一化并輸出，得到最終預(yù)測(cè)結(jié)果。

3.2 評(píng)價(jià)指標(biāo)

1）RMSE(Root Mean Square Error)：

2）MAPE(Mean Absolute Percentage Error)：

3）Coefficient of Determination：

4 算例

采用某區(qū)域內(nèi)實(shí)際充電樁充電負(fù)荷數(shù)據(jù)進(jìn)行仿真實(shí)驗(yàn)，實(shí)驗(yàn)平臺(tái)在深度學(xué)習(xí)框架Keras上進(jìn)行，計(jì)算機(jī)條件是CPU:酷睿i7-7700、內(nèi)存：16G、GPU：1080Ti 11G。

圖4 充電負(fù)荷熱量圖

圖4給出部分時(shí)刻充電負(fù)荷矩陣的可視化示例，顏色越深代表此處的負(fù)荷越大。并根據(jù)式(11) 將矩陣歸一化，然后用DMSTN 模型進(jìn)行訓(xùn)練和預(yù)測(cè)，將過(guò)去500 個(gè)時(shí)刻的負(fù)荷矩陣當(dāng)作訓(xùn)練集，提前6 小時(shí)預(yù)測(cè)未來(lái)24 個(gè)小時(shí)的負(fù)荷矩陣，時(shí)間間隔為一小時(shí)，且模型采用過(guò)去48 小時(shí)的負(fù)荷矩陣當(dāng)作輸入，未來(lái)6 小時(shí)的負(fù)荷矩陣當(dāng)作輸出，最終實(shí)現(xiàn)同時(shí)預(yù)測(cè)未來(lái)6個(gè)時(shí)間步的多步預(yù)測(cè)。

將預(yù)測(cè)結(jié)果根據(jù)式(12) 進(jìn)行反歸一化，然后將圖片輸出。

為了找出模型的最優(yōu)單元堆疊個(gè)數(shù)和體現(xiàn)模型的優(yōu)越性，將分別用堆疊單元8-12 層的5種情況與STN 預(yù)測(cè)模型對(duì)比。表1 給出了兩個(gè)模型在不同單元層數(shù)下的5 次實(shí)驗(yàn)24 個(gè)預(yù)測(cè)時(shí)刻的平均結(jié)果評(píng)價(jià)，包含了MAPE、MSE、R2三種評(píng)價(jià)指標(biāo)。從中可以看出兩個(gè)模型都是在單元層數(shù)為12 時(shí)模型能夠取得最好的效果，以上三個(gè)指標(biāo)的都是DMSTN 效果更好，可以看出DMSTN 模型具有更高的預(yù)測(cè)準(zhǔn)確率。

選擇了最優(yōu)單元層數(shù)12 以后，最優(yōu)預(yù)測(cè)結(jié)果和真實(shí)值進(jìn)行對(duì)比。圖中給出了24 個(gè)時(shí)刻的預(yù)測(cè)負(fù)荷矩陣，可以看出預(yù)測(cè)結(jié)果與真實(shí)結(jié)果比較一致，說(shuō)明了該算法的有效性。

表1 DMSTN和STN的比較

圖5 誤差對(duì)比

然后將預(yù)測(cè)誤差MAE 與真實(shí)值的對(duì)比見(jiàn)圖6，黃色代表DMSTN 誤差，紅色代表STN 誤差。圖中誤差是提前6 個(gè)小時(shí)預(yù)測(cè)誤差，所以誤差周期大概為6?？梢钥闯鯠MSTN 的每6 個(gè)預(yù)測(cè)點(diǎn)誤差大小不會(huì)隨著預(yù)測(cè)時(shí)間變長(zhǎng)而慢慢變高很多，而STN 的誤差是一直在遞增的，除了前幾個(gè)點(diǎn)的誤差和DMSTN 相差不大，后面的誤差越來(lái)越大，這充分說(shuō)明了DMSTN 模型的多步預(yù)測(cè)，可以有效的降低滾動(dòng)預(yù)測(cè)帶來(lái)的誤差。

為了驗(yàn)證DMSTN 算法的優(yōu)越性，將其他5種算法做上述相同實(shí)驗(yàn)，得出表2 中6 種算法提前6 個(gè)小時(shí)預(yù)測(cè)未來(lái)24 小時(shí)的RMSE 數(shù)據(jù)，可以看出DMSTN 算法對(duì)比其他幾種算法精度高出很多，由此驗(yàn)證該算法的有效性。

表2 不同算法間的比較

5 結(jié)束語(yǔ)

本文提出一種DMSTN 模型，用于解決電動(dòng)汽車(chē)短期時(shí)空動(dòng)態(tài)負(fù)荷問(wèn)題，取得了一定的效果，對(duì)比STN 模型有更高的準(zhǔn)確率。在12層單元時(shí)DMSTN 的MAE 和RMSE 誤差平均值比11 層單元的STN 低18.32% 和33.72%，DMSTN 的R2得分比STN 高出8.53%?？梢钥闯鯠MSTN 模型的預(yù)測(cè)精確度具有明顯提升。對(duì)比其他幾種算法進(jìn)度也有很大提高。

該方法不僅可以消除滾動(dòng)預(yù)測(cè)帶來(lái)的誤差提高預(yù)測(cè)精度，而且可以預(yù)測(cè)電動(dòng)汽車(chē)充電負(fù)荷的時(shí)空動(dòng)態(tài)，能給電網(wǎng)帶來(lái)時(shí)空二維的負(fù)荷信息，更大程度的幫助電網(wǎng)解決充電汽車(chē)大量入網(wǎng)帶來(lái)的問(wèn)題。