傅質馨,溫順潔,朱俊澎,袁 越
(1.河海大學 能源與電氣學院,南京 211100;2.河海大學 可再生能源發電技術教育部工程研究中心,南京 210098)
負荷特性研究作為電力系統研究領域的一項重要內容,能有效地提高電網運行的經濟性及科學性[1]。電力負荷特性受到多種隨機因素的影響,具有一定的不確定性,氣象因素對電力負荷的影響尤為顯著[2]。從實際應用角度來看,溫度的變化會使用戶對電能的需求量大幅度變動,增大電網的峰谷差和最高用電負荷,影響電網的發電與調控,導致停電事故發生的可能性上升。因此對溫度敏感負荷進行辨識至關重要,可以擴展到負荷模型的具體應用研究上。
溫度因素對負荷特性的影響最為顯著。文獻[3]基于基準負荷比較法測算地區冬夏兩季的空調負荷,通過二次多項式擬合,量化溫度對冬夏兩季空調負荷的影響。文獻[4]利用線性和非線性回歸分析方法對陜西地區夏季降溫負荷與氣溫之間的關系進行定性定量分析研究,劃分溫度區間后得到夏季降溫負荷與氣溫變化之間的幅值增減函數關系。文獻[5]在將負荷進行標幺化處理的基礎上,建立基于靈敏度的基礎預測模型,計算由氣象變化引起的冬季取暖負荷。
對負荷資源構成辨識方法,文獻[6]基于差值思想,實現對母線分布式光伏的辨識。文獻[7]建立卷積神經網絡模型,結合支持向量機對變電站負荷進行分解,分解精度比常規單模型負荷分解方法提高10%,具有很大的優越性。文獻[8]采用獨立分量分析法將變電站負荷曲線分解為不同的負荷類型,如居民負荷或者商業負荷。以上研究一方面尚未考慮氣象因素對負荷辨識影響,往往導致辨識出現較大誤差,由此本文利用皮爾遜相關系數篩選溫度敏感負荷,同時采用多項式回歸模型量化溫度對負荷的影響;另一方面,以上研究多采用傳統的卷積神經網絡,存在訓練效率低、辨識精度不高等問題,由此本文采用殘差卷積神經網絡(residual convolution neural network,ResCNN)進行負荷辨識的方法。
本文首先利用基準負荷比較法,在區分工作日和非工作日的基礎上,對敏感負荷進行定性定量分析研究。其次,對于溫度敏感負荷,將每小時的負荷與溫度的多項式系數作為負荷特征加入到動態負荷特征庫中。最后,采用傳統的卷積神經網絡(convolution neural network,CNN)模型和本文所提出的ResCNN模型對負荷構成進行辨識對比實驗。
本文在基準負荷比較法[9]的基礎上,對溫度敏感負荷進行測算,具體過程如下:
(1)采用日平均負荷最小的日負荷曲線作為基準負荷;
(2)假定所選定的日平均負荷最小的負荷曲線為LW,d,h和LK,d,h,其中,W為工作日,K為非工作日,d為日,h為小時數;
(3)用工作日和非工作日最小負荷曲線的平均值代表工作日和非工作日的基礎負荷曲線,具體表達式如下

式中:n1為工作日天數;n2為非工作日天數;L基礎W,h為工作日基礎負荷曲線;L基礎K,h為非工作日基礎負荷曲線。
(4)工作日敏感負荷測算:工作日剩余負荷曲線L剩余W,d,h與基礎負荷曲線L基礎W,h相減,得到的差值即為工作日敏感負荷曲線L敏感W,d,h,公式如下

非工作日敏感負荷測算:非工作日剩余負荷曲線L剩余K,d,h與基礎負荷曲線L基礎K,h相減,得到的差值是非工作日敏感負荷曲線L敏感K,d,h,公式如下

本文所用數據為某地區2018年全年商業負荷數據及溫度數據,負荷數據采樣間隔1 h。基于該數據,分析氣溫與負荷的關聯程度,以此篩選得到溫度敏感負荷。
本文采用皮爾森相關系數用來描述變量間相關程度,取各企業每個時刻的敏感負荷與實時溫度進行相關性計算,其中與溫度強相關的敏感負荷是溫度敏感負荷,計算結果見表1。

表1 超市溫度敏感負荷與溫度相關系數Table 1 Temperature sensitive load and temperature correlation coefficient of supermarket
由表1可以看出,超市非工作日負荷的溫度相關性要大于工作日負荷的溫度相關性;在非工作日,負荷與溫度相關性較強時段主要集中在12:00—17:00,16:00時,溫度與負荷相關性最強;在工作日,負荷與溫度相關性較強時段主要集中在7:00—8:00,以及13:00—16:00,其中16:00左右,溫度與負荷相關性最強。旅館溫度敏感負荷與溫度相關系數如表2所示。

表2 旅館溫度敏感負荷與溫度相關系數Table 2 Temperature sensitive load and temperature correlation coefficient of hotel
由表2可以看出,旅館非工作日負荷的溫度相關性要明顯強于工作日負荷的溫度相關性;在非工作日,負荷與溫度相關性較強時段主要集中在11:00—17:00,19:00時,溫度與負荷相關性最強;在工作日期間,負荷與溫度相關性較強時段主要集中在11:00—17:00,19:00左右,溫度與負荷相關性最強。
由皮爾遜相關系數法篩選得到的與溫度相關性強的負荷時段作為溫度敏感負荷。
回歸分析不僅可以反映出自變量對因變量的影響程度,還可利用得到的回歸模型對數據趨勢進行預測。
本文采用回歸方式處理溫度敏感負荷特征構建動態負荷特征庫,體現了溫度敏感負荷與溫度之間的關系,可基于不同溫度下得到相對應該溫度下的溫度敏感負荷特征庫。
為了進一步分析溫度敏感負荷與實時溫度變化的規律,量化溫度因素的影響程度,這里根據詳細的歷史數據,提出線性擬合公式如下

式中:L敏感為敏感負荷;T為對應溫度;a為自變量T的系數;b為常數項。
對溫度敏感負荷進行分時段擬合,從表1、表2可以看出非工作日超市溫度敏感負荷高峰時段主要集中在12:00—17:00,工作日超市溫度敏感負荷高峰時段集中在13:00—16:00;工作日旅館溫度敏感負荷高峰時段主要集中在11:00—17:00,非工作日旅館溫度敏感負荷高峰時段主要集中在11:00—17:00。對以上幾個時段分別進行線性擬合,擬合關系如下:
①非工作日超市溫度敏感負荷與實時溫度,呈現正相關的關系。線性方程為y=3.637x-53.3,其中,擬合評價指標R2=0.770 8。②工作日超市溫度敏感負荷與實時溫度,呈現正相關的關系。線性方程為y=3.031x-17.65,其中,擬合評價指標為R2=0.488 4。③非工作日旅館溫度敏感負荷與實時溫度,呈現正相關的關系。線性關系為y=4.402x-64.59,其中,擬合評價指標為R2=0.8133。④工作日旅館溫度敏感負荷與實時溫度,呈現正相關的關系。線性關系為y=4.01x-52.84,其中,擬合評價指標為R2=0.674 7。
對溫度敏感負荷與實時溫度進行二次曲線擬合,擬合公式為

式中:L敏感為敏感負荷;T為對應溫度。
對以上幾個時段分別進行多項式擬合,擬合關系如下:
(1)非工作日超市溫度敏感負荷與實時溫度,呈現正相關的關系。溫度區間為12~38℃,二次多項式方程為y=0.102 8x2-1.586x+8.486,其中,擬合評價指標R2=0.802。
(2)工作日超市溫度敏感負荷與實時溫度,呈現正相關的關系。溫度區間為10~40℃,二次多項式方程為y=0.094 06x2-166x+36.23,其中,擬合評價指標R2=0.516。
(3)非工作日旅館溫度敏感負荷與實時溫度,呈現正相關的關系。溫度區間為15~37℃,二次多項式方程為y=0.066 74x2+0.903 4x-21.33,其中,擬合評價指標R2=0.819 5。
(4)工作日旅館溫度敏感負荷與實時溫度,呈現正相關的關系。溫度區間為15~35℃,二次多項式方程為y=0.097 1x2-0.947 7x+6.692,其中,擬合評價指標R2=0.692 7。
利用兩種擬合方法所得到的擬合評價指標對比如表3所示。

表3 兩種擬合方式擬合度對比Table 3 Comparison of two fitting methods
通過對該地區溫度敏感負荷與實時溫度的挖掘分析可得以下結論:商業區溫度敏感負荷多指空調負荷,且商業負荷的溫度敏感負荷有共同的特點,負荷高峰時段多集中在人們出行時段,因此休息日的相關系數要高于工作日的相關系數。其次,無論工作日還是休息日,對于與實時溫度的擬合,非線性擬合度要明顯高于線性擬合度,因此多項式回歸更適用于對溫度敏感性負荷與溫度之間關系的描述。
本文提出溫度敏感負荷辨識問題的數學模型可表示如下

式中:Ltotal(1×T)為由多類溫度敏感負荷所構成的總負荷;T為采樣點數;Lsub(n×T)為當前負荷特征庫中所有溫度敏感負荷種類;n為負荷種類數目;A(1×n)為判斷典型負荷有無的向量,由0或1組成,為神經網絡待求量。
為了求解典型負荷的標簽量,將上式進行變形,得到典型負荷標簽量的求解公式如下

式中:pinv為求解偽逆矩陣的運算,若pinv(Lsub(n×T))存在唯一解,則輸入Ltotal(1×T),數學計算便可得到A(1×n)。
但僅僅采用傳統數學方法,往往會使pinv(Lsub(n×T))多解或無解。因此,需要嘗試采用有強大非線性映射能力的方法來求解該問題。
本文提出利用非線性擬合方式計算高峰負荷時段每小時的擬合系數,構建動態的溫度敏感負荷特征庫。負荷特征庫如表4所示。

表4 負荷特征庫Table 4 Load signature database
該負荷特征庫是一個隨溫度動態變化的數據庫,本文選取夏季典型日溫度,將典型日24點溫度數據代入到特征庫中,得到每個特征下24點負荷數據。典型溫度敏感負荷特征庫如表5所示。

表5 典型溫度敏感負荷特征庫Table 5 Fypical temperature sensitive load signature database
4.3.1 改進的卷積神經網絡模型
傳統的卷積神經網絡達到一定深度后,增加層數并不能帶來進一步地分類性能地提高,反而會導致網絡收斂變得很慢,分類準確率也會有所下降。利用ResCNN,可通過殘差結構使用多個有參層來學習輸入輸出之間的殘差表示,而非直接學習輸入輸出之間的映射。
本文所采取的改進卷積神經網絡結構采用結合殘差結構的CNN網絡來提取溫度敏感特征。殘差卷積神經網絡結構如圖1所示。

圖1 殘差卷積神經網絡結構Fig.1 Residual convolution neural network structure
(1)數據輸入層
訓練集:80 000條母線負荷序列Ltotal(1×T)以及所含溫度敏感負荷標簽向量A(1×n)。
測試集:20 000條母線負荷序列Ltotal(1×T)。
(2)卷積層
卷積層計算式[10]為

(3)批量歸一化層
在網絡中加入批量歸一化層不僅極大提升了訓練速度,收斂過程大大加快,還能增加分類效果。
(4)激活函數層
卷積層的激活函數選取ReLU相比傳統tanh函數,能有效克服梯度消失的問題以及加快了訓練的速度。
(5)殘差結構模塊
殘差塊[11]可表示為

式中:h(x l)=W l′x。Wl′為1×1卷積操作;F(x l,W l)為殘差部分,由3個卷積操作構成。
(6)扁平化層
將拼接輸出扁平化,傳入全連接層。
(7)全連接層
全連接層激活函數為sigmoid,輸出維度為負荷特征庫中負荷種類數目。
4.3.2 仿真分析
所用數據為某地區2018年全年商業負荷數據和溫度數據,通過以上方法得到的溫度敏感特征庫將作為神經網絡的輸入,輸入的母線負荷序列由該特征庫中的負荷所構成。
網絡訓練過程的目標是以最小化損失函數的方式學習每層網絡的參數。使用二分類交叉熵(binary cross entropy,BCE)作為損失函數來計算輸出層中的誤差。在反向傳播期間,使用Adam優化器更新模型參數,這是處理大型數據集的常用選擇。
殘差卷積神經網絡的相關參數學習率、反向傳播損失函數、批數據數量、優化算法、評價指標分別設 置 為0.000 1、binary_crossentropy、64、Adam、Accuracy。
卷積層卷積核個數為24,步長為1,填充方式為same,激活函數為ReLU。對于分類問題全連接層的激活函數為Sigmoid,而對于回歸問題全連接層的激活函數為Softmax。對于本文的分類問題,激活函數應選擇Sigmoid。
基于上述模型參數的設置,本文分工作日和非工作日分別對傳統卷積神經網絡和改進的卷積神經網絡進行辨識仿真分析,采用準確率來作為評價指標。特征庫數據采用4種不同縮放方式,仿真結果如表6和表7所示。

表6 非工作日辨識準確率對比結果Table 6 Comparison results of non-working days recognition accuracy%

表7 工作日辨識準確率對比結果Table 7 Comparison results of working days recognition accuracy%
由以上仿真分析結果可得,本文所提出的ResCNN模型能很好地根據當前總負荷識別出所含溫度敏感負荷的種類。4種不同的數據縮放方式中,ResCNN模型均比傳統的CNN測試結果平均辨識準確率提高10%。
本文對溫度敏感負荷辨識方法展開研究。首先,本文基于基準負荷比較法對敏感負荷進行測算,利用皮爾遜相關系數篩選得到與溫度強相關的溫度敏感負荷。其次,針對溫度敏感負荷,比較兩種回歸分析方法,得出非線性回歸分析更適于表征溫度敏感負荷特征。最后,本文提出利用溫度敏感負荷與溫度間的非線性擬合系數,構建動態的溫度敏感負荷特征庫,提出用ResCNN結構實現對溫度敏感負荷的辨識,辨識結果表明,較傳統的CNN結構,無論對于何種數據縮放方式,辨識準確率均提高10%。D