周 宇 劉海璇 胡衛(wèi)豐 胥 崢 劉志明 袁曉玲

（1.國網(wǎng)江蘇省電力有限公司 南京 210008）（2.中國電力科學研究院 南京 210003）（3.河海大學能源與電氣學院 南京 211100）

1 引言

影響用戶用電行為的因素眾多，需要從時間、空間、用戶類型等多個維度，分析自然、社會各種潛在因素與用電行為的關(guān)聯(lián)關(guān)系，并深度挖掘因素集與用戶用電行為的關(guān)聯(lián)強度［1］。影響用戶用電行為模式的各個因素關(guān)聯(lián)強度亦不相同，為此，需要利用關(guān)聯(lián)挖掘技術(shù)分析影響因素集合與相關(guān)性或因果結(jié)構(gòu)［2］。從而針對特定的用戶行為，用電行為集合之間的頻繁模式、關(guān)聯(lián)性、從眾多因素中發(fā)現(xiàn)提煉與其關(guān)聯(lián)的影響要素，建立精細化用電模式與影響因素的關(guān)聯(lián)模型，從而進一步理解用戶用電行為［3］。

目前，隨著泛在電力物聯(lián)網(wǎng)進程的提出和不斷推進，眾多信息系統(tǒng)和平臺在有效管理各類智能設(shè)備、促進智能電網(wǎng)建設(shè)、提升公司電網(wǎng)運營管理水平的同時［4］，電網(wǎng)企業(yè)存儲了海量數(shù)據(jù)，這些數(shù)據(jù)直觀上沒有表現(xiàn)出任何有價值的信息，但采用數(shù)據(jù)挖掘的方法，數(shù)據(jù)中隱藏的信息就會被提取出來［5～6］。文獻［7］考慮城市微氣象和節(jié)假日等影響因素，建立用電行為差異化模型，精細化挖掘用戶用電行為特征。文獻［8～9］利用數(shù)據(jù)挖掘和數(shù)理統(tǒng)計的相關(guān)技術(shù)，對一個地區(qū)不同行業(yè)的用電行為進行分析，提取出各個行業(yè)的用電特征和客戶的群體特性。文獻［10］以高斯濾波技術(shù)消除初始數(shù)據(jù)中的噪聲，結(jié)合k-means、SOM等聚類算法建立聚類分析模型，獲得隱藏在數(shù)據(jù)集中典型用電模式及其合適的聚類數(shù)目。文獻［11］和文獻［12］中的聚類方法對所有用戶都提取同一組特征量對用戶負荷進行聚類研究，無法充分掌握用戶的負荷特性。文獻［13］提出基于差異化特征量提取的分層聚類方法實現(xiàn)用戶用電行為的分類識別。

本文利用UCI（University of California，Irvine）數(shù)據(jù)庫［14］公布的法國克拉馬地區(qū)某用戶的用電數(shù)據(jù)進行聚類分析?；贙均值聚類算法計算時間短、速度快、聚類結(jié)果容易解釋的優(yōu)點［15～16］，通過合理的選擇聚類初始K值來提高聚類的準確度，精確挖掘出用戶的不同用電行為，為客戶用電精細化管理和提供優(yōu)質(zhì)用電服務(wù)等提供支持。

2 基于K均值聚類算法的居民用戶電力負荷特征分析

為實現(xiàn)對居民電力負荷特征的分析，首先要準確提取用戶的用電特征曲線。本文的居民用戶電力負荷特征分析框架如圖1所示。首先對原始的居民用戶用電數(shù)據(jù)進行預(yù)處理，這些預(yù)處理包括缺失數(shù)據(jù)處理、數(shù)據(jù)歸一化。然后對預(yù)處理完后的數(shù)據(jù)采用K均值算法提取的用戶負荷特征曲線，根據(jù)不同日類型下的負荷特征曲線，分析其用電行為特征。將聚類分析技術(shù)應(yīng)用到電力大數(shù)據(jù)中，有助于電力公司進行需求側(cè)管理、用戶細分、用電優(yōu)化建議等［17～19］。

圖1 電力負荷特征分析框架

2.1 K均值算法概述

K均值聚類算法是劃分聚類分割的方法，其工作原理是：首先隨機從給定的數(shù)據(jù)集中選中K個點，每個點代表每個簇的初始聚類中心，然后計算剩余各個樣本到聚類中心的歐式距離，將其歸為離它最近的那個簇，接著重新計算每一個簇的平均值，整個過程不斷重復(fù)，直到平方誤差準則函數(shù)最?。?0］。平方誤差準則定義為

式中，k為聚類個數(shù)，ti為第i類中的樣本的個數(shù)，mi是第i類中樣本的均值。

K均值聚類算法是典型的基于距離的聚類算法，采用點與點之間的距離作為相似性評價指標，即認為兩個對象的距離越近，相似度就越大。距離算法一般采用歐氏距離，如式（2）。

式中，xi是樣本X的第i個變量值；yi是樣本Y的第i個變量值。這種算法認為簇是由距離靠近的對象組成，因此把得到緊湊且獨立的簇作為最終目標。

2.2 數(shù)據(jù)預(yù)處理

2.2.1 數(shù)據(jù)清洗

進行聚類分析的數(shù)據(jù)來自生產(chǎn)、生活、商業(yè)中的實際數(shù)據(jù)，在現(xiàn)實世界中，由于各種內(nèi)部或外部影響，導(dǎo)致數(shù)據(jù)會出現(xiàn)缺漏或者異常數(shù)據(jù)的存在。為了提高最終結(jié)果的可信性和可解釋性，在進行聚類之前對這些壞數(shù)據(jù)進行剔除和修復(fù)。

針對數(shù)據(jù)集中出現(xiàn)的缺失值，利用向前向后移動平均法對缺失數(shù)據(jù)補齊。設(shè)缺失數(shù)據(jù)Li是某一天負荷曲線的第i個數(shù)據(jù)點，修正后的數(shù)據(jù)Li"為

式中，Li-h和Li+g分別代表Li向前h個數(shù)據(jù)和向后g個數(shù)據(jù)；h1和g1一般可取5～10。

2.2.2 數(shù)據(jù)標準化

為了更明顯地體現(xiàn)用戶用電行為的動態(tài)變化，對用戶的原始負荷數(shù)據(jù)進行標準化處理：對原始數(shù)據(jù)進行歸一化處理，使最終結(jié)果落入［0，1］區(qū)間，歸一化公式如式（4）：

式中xi和x*分別表示第i個采樣時刻實際的用電負荷和標準化之后的用電負荷，xmax和xmin分別表示每一天樣本數(shù)據(jù)中的負荷最大值和最小值，n為每一天樣本中的負荷采樣點數(shù)。

2.3 初始K值的確定

雖然K均值算法具有簡單、效果好的優(yōu)點，但算法中的K值是事先給定的，而K值的大小往往難以估計，不同的K值對聚類結(jié)果的準確性也會產(chǎn)生影響。目前常用的確定K值的方法有輪廓圖法、誤差平方和法等［21］。本文采用輪廓圖法作為初始K值的確定方法。

輪廓圖是利用輪廓系數(shù)做出的圖形。輪廓系數(shù)是聚類效果好壞的一種評價方式，最早由Peter J.Rousseeuw在1986年提出，它結(jié)合內(nèi)聚度和分離度兩種因素，用來在相同原始數(shù)據(jù)的基礎(chǔ)上評價不同算法、或者算法不同運行方式對聚類結(jié)果所產(chǎn)生的影響［22］。輪廓圖上第i點的輪廓值定義為

式中，a是第i點與同一個簇中其他點的平均距離。b是向量，其元素表示第i點與其他不同簇中各點的平均距離。S(i)的取值范圍是［-1，1］，S(i)的值接近于1說明點i更傾向于當前的類，S(i)接近0表示點i傾向于在兩個類之間，S(i)接近-1表示點i傾向于其他某一個類。

平均輪廓值就是將所有點的輪廓系數(shù)求和取平均，平均輪廓值越大，表示聚類效果越好。

3 算例分析

本文利用UCI數(shù)據(jù)庫公布的2007年某智能小區(qū)的一戶居民的用電數(shù)據(jù)集進行分析。數(shù)據(jù)集包含該家庭每隔1分鐘一次電量消耗測量數(shù)據(jù)。

3.1 初始K值的確定

在進行初始聚類數(shù)K值確定之前，首先對預(yù)處理后的負荷數(shù)據(jù)進行每60分鐘求和一次，所得的數(shù)據(jù)即為日24點數(shù)據(jù)，處理完后的數(shù)據(jù)采用輪廓圖法確定其最佳聚類類別數(shù)。圖2分別表示聚類類別為2、3、4時的輪廓圖。表1為不同K值對應(yīng)的平均輪廓值，根據(jù)2.3節(jié)輪廓圖以及平均輪廓值的定義，確定選擇此數(shù)據(jù)集類別數(shù)為2。

圖2 不同K值下的輪廓圖

表1 不同K值對應(yīng)的平均輪廓值

3.2 聚類結(jié)果的分析

本文以數(shù)據(jù)集中2007年一月份的數(shù)據(jù)為例，隨機選取初始聚類中心，得到兩類居民日負荷特征曲線如圖3所示。

圖3 居民一月日負荷曲線及特征曲線

根據(jù)圖3的聚類結(jié)果，該用戶的一月的日負荷特征曲線主要有兩類。這兩類特征曲線變化趨勢有明顯的差異。

結(jié)合圖4中的空調(diào)熱水器負荷和廚房負荷可以看出：第一類負荷特征曲線在11：00～14：00出現(xiàn)一個負荷高峰，與空調(diào)熱水器負荷高峰出現(xiàn)的時間一致。第二類負荷特征曲線有兩個用電負荷高峰，分別出現(xiàn)在上午8：00～10：00和晚上的18：00～22：00，在這兩個用電高峰期間，空調(diào)熱水器負荷基本一直處于大功率運行的狀態(tài)中，而廚房負荷運行時間分散且功率較小，說明該用戶用電追求舒適性，空調(diào)熱水器類等對用戶用電舒適性有較強相關(guān)性負荷在該家庭中占有很大的比重。

圖4 負荷曲線的熱力圖

由于空調(diào)熱水器等負荷對用戶的舒適性影響很大，且廚房負荷的使用時間較為固定，考慮到洗衣房負荷對用戶舒適性影響較小且用電時間較為分散，為了緩解用電高峰時段電力系統(tǒng)的調(diào)峰壓力，將聚類一和聚類二中用電高峰時段的洗衣房負荷轉(zhuǎn)移到用電低谷期。仿真分析采用文獻［23］國內(nèi)實施的用電峰谷時段劃分，轉(zhuǎn)移后得到的結(jié)果如圖5所示。

圖5 負荷轉(zhuǎn)移前后聚類結(jié)果

由圖5可以看出，轉(zhuǎn)移后的負荷特征曲線的用電峰值比負荷轉(zhuǎn)移前都有所下降。第一類特征曲線在11：00～13：00出現(xiàn)的負荷高峰轉(zhuǎn)移到了夜間1：00～4：00。而第二類特征曲線19：00出現(xiàn)的負荷轉(zhuǎn)移到了凌晨負荷比較低的時段。

4 用戶用電行為特征分析

為了進一步分析用戶的用電行為，本文提取節(jié)假日和不同季節(jié)典型負荷特性曲線，對用戶的用電行為特征進一步分析。

4.1 節(jié)假日用電行為

法國的主要節(jié)假日有元旦和勞動節(jié)共6天法定節(jié)假日以及萬圣和圣誕等四個宗教節(jié)日。文中將這21天的節(jié)假日負荷數(shù)據(jù)利用K均值聚類方法提取出節(jié)假日的特征曲線，當K取2時，平均輪廓值為0.4847，大于K取其它值時的平均輪廓值，最終得到聚類結(jié)果如圖6。

圖6 節(jié)假日用電負荷曲線及特征曲線

根據(jù)圖6的聚類結(jié)果可以看出，節(jié)假日有兩類特征曲線，第一類用電負荷特征曲線有明顯的用電高峰區(qū)和用電低谷區(qū)，用電高峰出現(xiàn)在晚上的19：00～21：00，凌晨1：00～6：00為用電低谷區(qū)；第二類用電負荷曲線一天之內(nèi)的波動比第一類小，用電高峰出現(xiàn)在中午12：00和凌晨1：00。

為了對節(jié)假日用戶用電行為作進一步分析，利用熱力圖將這21天的空調(diào)熱水器負荷、廚房負荷和洗衣房負荷曲線表示出來，得到結(jié)果如圖7所示。

圖7 負荷曲線的熱力圖

由圖7節(jié)假日期間空調(diào)熱水器負荷和廚房負荷可以看出，第一類用電曲線代表的節(jié)假日各類負荷的使用頻率明顯高于第二類用電曲線代表的節(jié)假日。進一步，通過對比兩類節(jié)假日廚房用電曲線，在圣誕節(jié)（屬于第一類節(jié)假日用電特征曲線）這種重大的宗教節(jié)日期間，廚房負荷遠高于第二類節(jié)假日的廚房負荷，說明用戶在這些節(jié)假日期間，用戶更可能在家里與親人朋友團聚；第二類節(jié)假日特征曲線有多天廚房負荷為0，說明用戶在元旦（屬于第二類節(jié)假日用電特征曲線）這種比較長的假期外出游玩或外出用餐。

由于聚類二中洗衣房負荷都處于閑置的狀態(tài)，只針對聚類一中的廚房負荷進行轉(zhuǎn)移，得到的結(jié)果如圖8。

圖8 負荷轉(zhuǎn)移前后聚類結(jié)果

根據(jù)圖8，用電高峰期的洗衣房負荷轉(zhuǎn)移到凌晨時段，使得節(jié)假日第一類用電負荷曲線在白天時段的用電波動變得更加平緩。

4.2 季節(jié)用電行為規(guī)律

為了分析該地區(qū)用電負荷隨氣候狀況變化的規(guī)律，本文結(jié)合該地區(qū)的氣候特征，分析用電行為與氣候之間的關(guān)系。文中采用傳統(tǒng)意義上的典型季節(jié)劃分：春季指每年的3～5月，夏季指每年的6～8月，秋季指每年的9～11月，冬季為12月和次年的1～2月。

圖9 各季節(jié)用電負荷特征曲線

當K=2時，四個季節(jié)的平均輪廓值均為最大，選取聚類結(jié)果中包含天數(shù)較多的那一類作為該季節(jié)的典型用電負荷曲線，如圖9所示。從圖中可以看出，無論哪一個季節(jié)，該用戶用電負荷的日變化均具有明顯的波動特征。凌晨2：00～6：00之間用電負荷較小，日最低用電負荷值出現(xiàn)在凌晨4點左右。從早晨6：00開始至上午8：00負荷逐漸增大，在8：00達到了日變化的最大值。下午1：00出現(xiàn)用電負荷的一個低值區(qū)，下午6：00之后，出現(xiàn)了一天中第二次用電負荷的快速增長，在夜間10：00左右達到了一天的最大值。由于該用戶位于法國克拉馬地區(qū)，該地區(qū)冬季溫度較低，夏季涼爽，氣溫年差較小，所以該用戶春、夏、秋三個季節(jié)的特征曲線變化趨勢基本一致，用戶第二次負荷增長出現(xiàn)的時間以及最大值出現(xiàn)時間基本一致。由于冬季氣溫較低以及天黑時間提前，所以冬季的第二次用電負荷增長出現(xiàn)的時間提前，最高值出現(xiàn)在20：00左右。

5 結(jié)語

1）本文利用K均值算法對居民用戶用電負荷數(shù)據(jù)進行聚類分析，為了尋求最優(yōu)初始聚類數(shù)K值，采用輪廓圖和平均輪廓值確定最優(yōu)的初始聚類數(shù)，該方法可以對分類的合理性進行判斷，解決了K均值算法對初始K值敏感的問題；

2）通過K均值算法提取用戶的日負荷特征曲線、節(jié)假日特征曲線以及各季節(jié)特征曲線，結(jié)合空調(diào)熱水器、廚房負荷和洗衣房負荷熱力圖，得出居民在不同情境下的用電消費習慣，對居民用戶的用電行為進行分析。

3）通過將對用戶舒適性影響不大的洗衣房負荷進行轉(zhuǎn)移，可以有效改善用戶的用電負荷曲線。

用戶用電行為特征分析結(jié)果有助于電力公司根據(jù)用戶的用電行為設(shè)計相應(yīng)的需求響應(yīng)激勵機制，提高居民用戶對需求響應(yīng)的參與度，使精細化和實時化的需求響應(yīng)工作更易實現(xiàn)。

隨著泛在電力物聯(lián)網(wǎng)技術(shù)的發(fā)展以及負荷監(jiān)測技術(shù)的不斷進步，更細粒度的電能消耗數(shù)據(jù)將會更容易獲得，基于這些數(shù)據(jù)的數(shù)據(jù)挖掘?qū)τ秒娯摵深A(yù)測、用戶的用電個性化管理以及用能診斷提供依據(jù)。