劉春蕾 龐鵬飛 石紋赫 孔令號(hào) 黃 洵 戚 軍
(1.國網(wǎng)河北省電力有限公司保定供電分公司 保定 071000;2.浙江工業(yè)大學(xué)信息工程學(xué)院 杭州 310000)
近年來,負(fù)荷監(jiān)測(cè)作為智能用電技術(shù)體系的重要組成部分[1],實(shí)現(xiàn)了電力公司與居民用戶之間的雙向互動(dòng)機(jī)制,是目前各國研究的熱點(diǎn)[2]。負(fù)荷監(jiān)測(cè)技術(shù)通過對(duì)綜合用電數(shù)據(jù)進(jìn)行深度挖掘,不僅能夠?yàn)殡娋W(wǎng)側(cè)提供負(fù)荷構(gòu)成信息,為電網(wǎng)的優(yōu)化調(diào)度提供助力,還能夠反饋于用戶,解析用戶電能消費(fèi)結(jié)構(gòu)[3],有效提高電能的利用率和峰谷電價(jià)引導(dǎo)下的削峰填谷積極性。
目前負(fù)荷監(jiān)測(cè)可大致分為侵入式監(jiān)測(cè)(Intrusive load monitoring, ILM)和非侵入式監(jiān)測(cè)(Non-intrusive load monitoring, NILM)。ILM 需要在設(shè)備上安裝監(jiān)測(cè)傳感器,雖然監(jiān)測(cè)精度高且無需復(fù)雜的識(shí)別算法,但其硬件安裝成本高,且在更加注重個(gè)人隱私的信息化時(shí)代,用戶接受程度也較低[4-6]。與之相對(duì),由HART 于20 世紀(jì)80 年代首先提出的NILM 技術(shù)只需要在配電進(jìn)線處安裝監(jiān)測(cè)設(shè)備,對(duì)監(jiān)測(cè)數(shù)據(jù)進(jìn)行負(fù)荷分解,即可實(shí)現(xiàn)對(duì)各設(shè)備的監(jiān)測(cè)。該方法硬件成本低,用戶接受程度高,在工業(yè)設(shè)備監(jiān)測(cè)以及常用居民負(fù)荷監(jiān)測(cè)等領(lǐng)域得到了廣泛的應(yīng)用[7-11]。文獻(xiàn)[8-9]基于低頻電力數(shù)據(jù),采用穩(wěn)態(tài)功率信息構(gòu)建負(fù)荷識(shí)別特征,特征提取速度快,但是實(shí)際工程應(yīng)用中容易出現(xiàn)特征重疊問題導(dǎo)致識(shí)別準(zhǔn)確率降低。為提高復(fù)雜場(chǎng)景下的識(shí)別正確率,可采用包含更多信息量的負(fù)荷特征。文獻(xiàn)[10]基于電流時(shí)域波形和頻域諧波信息,采用主成分分析法降維構(gòu)建負(fù)荷特征。文獻(xiàn)[11]則采用差異分析方法構(gòu)建易于區(qū)分的融合特征作為負(fù)荷特征。這兩種方法都能提高負(fù)荷識(shí)別準(zhǔn)確率,但與此同時(shí)負(fù)荷識(shí)別的速度也明顯降低,在實(shí)際工程應(yīng)用中影響負(fù)荷識(shí)別的實(shí)時(shí)性。
本文基于對(duì)各類負(fù)荷特征的對(duì)比分析,綜合考慮負(fù)荷特征提取代價(jià)及識(shí)別效果構(gòu)建雙層特征組,采用基于支持向量機(jī)(Support vector machine,SVM)多分類算法的負(fù)荷動(dòng)態(tài)識(shí)別方法按需識(shí)別負(fù)荷類型。最后通過Matlab 仿真及實(shí)際工程應(yīng)用對(duì)所提NILM 動(dòng)態(tài)識(shí)別方法的性能進(jìn)行測(cè)試驗(yàn)證。
負(fù)荷特征作為負(fù)荷識(shí)別的依據(jù),直接影響著NILM 識(shí)別算法的識(shí)別準(zhǔn)確度及速度。表1[12-16]總結(jié)了實(shí)際工程中較為常用的負(fù)荷特征的提取代價(jià):采樣頻率是線路電壓、電流測(cè)量裝置的采樣頻率;監(jiān)測(cè)窗口寬度是提取負(fù)荷特征所需的測(cè)量數(shù)據(jù)長(zhǎng)度;計(jì)算復(fù)雜度表示負(fù)荷特征提取算法的復(fù)雜程度;負(fù)荷特征存儲(chǔ)量是記錄負(fù)荷特征所需要的硬件空間。
表1 負(fù)荷特征綜合分析表
從表1 可以看出,穩(wěn)態(tài)特征中除電流諧波特征外,穩(wěn)態(tài)功率特征(P、Q)以及電壓-電流軌跡特征對(duì)采樣頻率要求較低,監(jiān)測(cè)窗口寬度普遍較小,電壓-電流軌跡及電流諧波特征的計(jì)算復(fù)雜度和負(fù)荷特征存儲(chǔ)量明顯高于穩(wěn)態(tài)功率特征;暫態(tài)特征的采樣頻率普遍較高,監(jiān)測(cè)寬度窗口寬度取決于具體的負(fù)荷類型,一般不超過100 s,除S 變換特征外,其余暫態(tài)電流特征計(jì)算復(fù)雜度都較低且負(fù)荷特征存儲(chǔ)量較??;時(shí)間特征的采樣頻率和計(jì)算復(fù)雜度較低且負(fù)荷特征存儲(chǔ)量較小,但所需監(jiān)測(cè)窗口寬度往往較大。在實(shí)際工程應(yīng)用中,時(shí)間特征在負(fù)荷停止運(yùn)行前往往無法識(shí)別負(fù)荷類型,實(shí)時(shí)性較差;穩(wěn)態(tài)功率特征容易出現(xiàn)特征重疊問題,通常與暫態(tài)特征結(jié)合使用。這種結(jié)合方法要求NILM 識(shí)別方法采用更多的特征或通過復(fù)雜的特征融合方法構(gòu)造新特征,從而導(dǎo)致特征提取代價(jià)較高以及實(shí)時(shí)性較差。
綜合考慮工程應(yīng)用中對(duì)識(shí)別準(zhǔn)確率及速度的要求,本文所提負(fù)荷動(dòng)態(tài)識(shí)別方法包含兩層特征組,如表2 所示。負(fù)荷識(shí)別時(shí),由簡(jiǎn)入繁,動(dòng)態(tài)選用負(fù)荷特征,首先選用第一層特征組,針對(duì)部分無法用第一層特征組區(qū)分的相似負(fù)荷,進(jìn)一步采用第二層特征組進(jìn)行區(qū)分。每一層特征組中的多類負(fù)荷特征類型,根據(jù)負(fù)荷特征提取代價(jià)進(jìn)一步區(qū)分優(yōu)先級(jí),根據(jù)實(shí)際應(yīng)用場(chǎng)景選取其中的一類或多類負(fù)荷特征。
表2 分層特征組
基于上述負(fù)荷特征分層思想,本文所設(shè)計(jì)的負(fù)荷動(dòng)態(tài)識(shí)別方法由分類模型訓(xùn)練和在線負(fù)荷識(shí)別兩部分構(gòu)成,詳細(xì)流程如圖1 所示。其中,分類模型訓(xùn)練基于已標(biāo)定負(fù)荷類型的訓(xùn)練/測(cè)試數(shù)據(jù)集,提取各層特征進(jìn)行模型訓(xùn)練和檢驗(yàn),形成分層的負(fù)荷識(shí)別模型供在線負(fù)荷識(shí)別使用;在線負(fù)荷識(shí)別部分實(shí)時(shí)采集電壓電流數(shù)據(jù)進(jìn)行投切事件監(jiān)測(cè),然后按需提取特征輸入至對(duì)應(yīng)層的負(fù)荷識(shí)別模型進(jìn)行判斷,最后輸出負(fù)荷類型識(shí)別結(jié)果。針對(duì)流程圖中的關(guān)鍵步驟,下文將詳細(xì)展開介紹。
步驟1:獲取負(fù)荷投切事件數(shù)據(jù)。針對(duì)于待識(shí)別的負(fù)荷群,試運(yùn)行一段時(shí)間并記錄共m次投切事件ETp(p=1,2,3,…,m)的暫態(tài)電壓vp(n)、暫態(tài)電流ip(n),然后根據(jù)實(shí)際負(fù)荷投切情況標(biāo)定ETp的負(fù)荷類型并記為yp。
步驟2:設(shè)置負(fù)荷特征層及準(zhǔn)確率閾值δ。首先根據(jù)實(shí)際工程需求設(shè)置識(shí)別準(zhǔn)確率閾值δ,然后根據(jù)表2 所述特征選取優(yōu)先級(jí)并結(jié)合實(shí)際工程應(yīng)用場(chǎng)景,確定第一層特征組及第二層特征組所采用的特征。
步驟3:提取第一層負(fù)荷特征。根據(jù)步驟2 所設(shè)置的負(fù)荷特征層,結(jié)合表1 所述提取方法,提取所有負(fù)荷投切事件的第一層特征數(shù)據(jù)集FL=[FL1…FLp…FLm]T,其中FLp為第p次投切事件的第一層特征向量。
步驟4:訓(xùn)練第一層識(shí)別模型??紤]到訓(xùn)練數(shù)據(jù)集的樣本數(shù)量較少且負(fù)荷類型往往不止兩類,采用小樣本環(huán)境下識(shí)別準(zhǔn)確率較高的SVM 多分類OAA(One-against-ALL)方法。識(shí)別模型的構(gòu)建步驟如下[17-19]所示。
(1) 針對(duì)c種負(fù)荷類型,OAA 共需要構(gòu)建c個(gè)SVM 模型。首先訓(xùn)練第一類負(fù)荷的SVM 模型,即初始化λ=1。
(2) 基于訓(xùn)練數(shù)據(jù)集中的負(fù)荷類型標(biāo)記yi以及第一層負(fù)荷特征集FL構(gòu)建負(fù)荷特征樣本訓(xùn)練集S以及測(cè)試集S',其中若ETp為第λ類負(fù)荷,則將目標(biāo)值yp設(shè)置為λ,若否則令yp設(shè)置為-1。