(華南理工大學(xué) 電力學(xué)院，廣州 510000)

0 引言

非侵入式負(fù)荷監(jiān)測(cè)分解[1-2](Non-Intrusive Load Monitoring，NILM)以最小的安裝維護(hù)成本，將總功率曲線分解為單獨(dú)電氣設(shè)備的功率曲線，實(shí)現(xiàn)了對(duì)用戶用電行為的捕捉，有利于制定電網(wǎng)運(yùn)行計(jì)劃，也能為消費(fèi)者找到最優(yōu)的用能方式。在基于暫態(tài)過程的NILM方法中，變點(diǎn)識(shí)別的準(zhǔn)確性和魯棒性[3]，對(duì)于正確捕獲設(shè)備啟/停時(shí)刻，截取設(shè)備暫態(tài)信息起到關(guān)鍵作用，從而為特征提取與負(fù)荷識(shí)別的準(zhǔn)確性提供充分的保證[4]。以往NILM的研究多關(guān)注于負(fù)荷印記的種類，以及識(shí)別負(fù)荷的方法，卻忽略了變點(diǎn)檢測(cè)的重要性。但隨著用電設(shè)備日益多樣化，特別是多狀態(tài)型、連續(xù)變化型負(fù)荷的增多，變點(diǎn)檢測(cè)方法的抗干擾能力和精確性，對(duì)NILM顯得更為重要。

變點(diǎn)，即一段數(shù)據(jù)中的異常點(diǎn)，在負(fù)荷分解中則定義為負(fù)荷的開關(guān)事件。變點(diǎn)識(shí)別最早由陳希孺教授提出[5]，經(jīng)過幾十年的發(fā)展，目前在信號(hào)處理的領(lǐng)域已經(jīng)有較為成熟的應(yīng)用。然而，在非侵入式負(fù)荷分解方面，變點(diǎn)識(shí)別的應(yīng)用尚不廣泛。文獻(xiàn)[6-7]對(duì)于變點(diǎn)的檢測(cè)僅僅依靠對(duì)于負(fù)荷總有功功率的變化量，然而這對(duì)于大功率電器與小功率電器混合運(yùn)行時(shí)的檢測(cè)效果比較差。文獻(xiàn)[7]提出了利用最優(yōu)化算法來實(shí)現(xiàn)對(duì)不同負(fù)荷組合的運(yùn)行狀態(tài)的識(shí)別，但只針對(duì)于負(fù)荷穩(wěn)定運(yùn)行時(shí)的工作狀態(tài)。文獻(xiàn)[8-9]提出了基于最大似然估計(jì)的思想來檢測(cè)變點(diǎn)，即認(rèn)為變點(diǎn)的前后服從兩個(gè)不同的分布，它們的統(tǒng)計(jì)量會(huì)有顯著的差異，從而達(dá)到識(shí)別變點(diǎn)的目的。然而最大似然估計(jì)的方法需要估計(jì)變點(diǎn)前后的分布參數(shù)，而目前對(duì)于變點(diǎn)前后的穩(wěn)態(tài)功率所服從的分布尚不能有一個(gè)準(zhǔn)確的定論。

最小二乘(Least Square, LS)法最早是由勒讓德(A.M.Legendre)于1805年提出[10]。它不須對(duì)模型中的隨機(jī)誤差的分布有特定的假設(shè), 且計(jì)算尚不很復(fù)雜,準(zhǔn)確度高。而非侵入式負(fù)荷分解中的變點(diǎn)識(shí)別的前后分布參數(shù)未知，但可視作負(fù)荷的穩(wěn)態(tài)運(yùn)行，可以構(gòu)建均值模型；另一方面，其計(jì)算量小，計(jì)算速度快，可以實(shí)現(xiàn)在線提取負(fù)荷暫態(tài)特征和負(fù)荷數(shù)據(jù)的傳輸。

然而隨著用電負(fù)荷多樣化，單純地利用最小二乘法計(jì)算統(tǒng)計(jì)量無法應(yīng)對(duì)大小功率負(fù)荷混合，連續(xù)變化型負(fù)荷等復(fù)雜的情況。而且在較高采樣頻率下，直接計(jì)算最小二乘法計(jì)算量非常大，這大大影響了負(fù)荷暫態(tài)特征提取的實(shí)時(shí)性。因此本文提出了基于模擬退火尋優(yōu)方式的改進(jìn)最小二乘法，通過最優(yōu)化算法，能快速尋優(yōu)，從而節(jié)省大量的計(jì)算成本。而且考慮到上述多電器同時(shí)運(yùn)行時(shí)的復(fù)雜情況，提出了自適應(yīng)閾值，從而確保較高的識(shí)別率。

1 改進(jìn)最小二乘法

1.1 均值變點(diǎn)模型

在用電設(shè)備進(jìn)行啟/?；虬l(fā)生狀態(tài)切換時(shí)刻前后，總用電負(fù)荷的電氣量會(huì)發(fā)生劇烈的變化(例如功率，電流)。這個(gè)過程定義為一個(gè)事件，且視為暫態(tài)過程。在這個(gè)過程前后一小段時(shí)間內(nèi)，可看做兩個(gè)穩(wěn)態(tài)階段，其電氣量不會(huì)發(fā)生劇烈變化，可認(rèn)為前后兩個(gè)穩(wěn)態(tài)分別服從兩個(gè)分布，且兩個(gè)分布的均值有較大變化，但方差變化較小，這定義為均值變點(diǎn)模型[11]。

設(shè)有一段時(shí)間總負(fù)荷數(shù)據(jù)X1,X2, …，Xn，n為樣本總量。假設(shè)其有q個(gè)變點(diǎn)，q<

Xi=mi+ei,i=1,2,...,n

(1)

μ1=…=μq(1)-1=M1,…,μq(q)=…=μn=Mq+1

其中：μ為每個(gè)X所服從分布的數(shù)學(xué)期望，而q個(gè)變點(diǎn)將樣本分為q+1個(gè)區(qū)間。M為每段區(qū)間內(nèi)的平均值。理想情況下，如果Mi≠M(fèi)i+1，則認(rèn)為是一個(gè)真正的變點(diǎn)。

ei為隨機(jī)誤差 , 假定為相互獨(dú)立，數(shù)學(xué)期望為0，方差σ2<∞。

最小二乘法的目標(biāo)函數(shù)為：確定q個(gè)變點(diǎn)的具體位置，使得總目標(biāo)函數(shù)T最小。

(2)

(定義Xq0=X1,Xq(q-1)=Xn)

1.2 模擬退火算法

接下來需要確定每個(gè)q的位置，從而確定T的最小值。由于樣本數(shù)據(jù)龐大，且有多個(gè)變點(diǎn)，無法直接通過求微分等符號(hào)數(shù)學(xué)手段求出，直接計(jì)算其數(shù)值解又需要花費(fèi)太多時(shí)間。因此采用啟發(fā)式搜索的模擬退火算法(Simulated Annealing, SA)。

模擬退火算法對(duì)初值不敏感，而且不容易陷入局部最優(yōu)的陷阱，收斂速度較快。它的核心是Metropolis準(zhǔn)則和退火過程兩部分。Metropolis采用概率來接受新狀態(tài)，而不是使用完全確定的規(guī)則。假設(shè)前一狀態(tài)為x(n)，系統(tǒng)受到一定擾動(dòng)，狀態(tài)變?yōu)閤(n+1)，系統(tǒng)能量由E(n)變?yōu)镋(n+1)，定義接受概率為P：

(3)

當(dāng)狀態(tài)轉(zhuǎn)移之后，如果能量減小，那么轉(zhuǎn)移被接受(P=1)。如果能量增大，就說明函數(shù)偏離全局或局部最優(yōu)更遠(yuǎn)，此時(shí)進(jìn)行概率操作，在[0,1]區(qū)間內(nèi)產(chǎn)生一個(gè)隨機(jī)數(shù)ξ，如果ξ

Metropolis保證了算法能搜索到全局最優(yōu)解，但直接使用Metropolis會(huì)導(dǎo)致尋優(yōu)的速度太慢，以致于失去了最優(yōu)化算法的意義，其關(guān)鍵在于退火參數(shù)的選擇。本文選擇指數(shù)式下降溫度的退火控制策略。使得算法不至于退火太快，導(dǎo)致算法在局部最優(yōu)即結(jié)束迭代。另外，為了減小不必要的搜索時(shí)間，提高算法效率，設(shè)定當(dāng)尋優(yōu)的結(jié)果趨于穩(wěn)定時(shí)終止迭代。

1.3 “變點(diǎn)有無”假設(shè)檢驗(yàn)

檢驗(yàn)原假設(shè)H0: 不存在變點(diǎn)；H1：存在變點(diǎn)。針對(duì)這個(gè)問題，引入樣本未除以自由度的方差：

(4)

隨后將樣本隨機(jī)分成兩部分，前一段有j-1個(gè)數(shù)據(jù)，后者有n-j+1個(gè)數(shù)據(jù)。然后計(jì)算這兩段數(shù)據(jù)的方差并求其和：

(5)

閾值C由概率論中的中心極限定理[12]得到：

C=σ2(2loglogn+logloglogn-logπ+xa)

(6)

一般情況下分布σ2未知，因此未知，所以采用一個(gè)相合且漸進(jìn)無偏的估計(jì)來代替[13]：

(7)

由此計(jì)算得到自適應(yīng)的閾值C進(jìn)行檢測(cè),可以使假設(shè)檢驗(yàn)有α的檢驗(yàn)水平。

2 變點(diǎn)識(shí)別算法

2.1 變點(diǎn)識(shí)別算法流程

本文提出的變點(diǎn)識(shí)別方法以負(fù)荷有功功率為識(shí)別對(duì)象。在變點(diǎn)識(shí)別中負(fù)荷的有功功率是最為明顯的特征，不像無功等其他特征，電阻性的負(fù)荷開/關(guān)則不會(huì)產(chǎn)生變化，從而干擾了算法的識(shí)別。具體算法流程如圖1所示。

圖1 變點(diǎn)識(shí)別算法流程

算法的核心在于兩個(gè)循環(huán)，第一個(gè)循環(huán)檢測(cè)每個(gè)窗口變點(diǎn)數(shù)目，并給出粗略的變點(diǎn)時(shí)刻。第2個(gè)循環(huán)針對(duì)于假設(shè)：“只有一個(gè)變點(diǎn)”的區(qū)間進(jìn)行檢驗(yàn)，以得到精確的變點(diǎn)時(shí)刻。

2.2 LS-SA算法

為了利用最小二乘法來估計(jì)窗口內(nèi)的變點(diǎn)數(shù)目，首先先確定窗口長(zhǎng)度winlen和窗口內(nèi)至多有q個(gè)變點(diǎn)。q應(yīng)取得相當(dāng)小(減少算法計(jì)算量)，相對(duì)于窗口長(zhǎng)度而言，但q又需滿足大于變點(diǎn)個(gè)數(shù)。在本文中，選取q=5。

接下來，對(duì)于q從1到5的情形，分別由公式(2)計(jì)算T1min-T5min，考慮到窗口內(nèi)沒有變點(diǎn)的情況，定義沒有變點(diǎn)的目標(biāo)函數(shù)T0min為窗口數(shù)據(jù)中未除以自由度的方差。

為了求解對(duì)應(yīng)于T1min-T5min的變點(diǎn)位置M1-M5，本文采用最小二乘-模擬退火算法(LS-SA)，即采用模擬退火方式來快速尋找LS的最優(yōu)解，其流程如表1所示。

表1 LS-SA算法

例如在電動(dòng)機(jī)型負(fù)荷啟動(dòng)曲線中，會(huì)產(chǎn)生較大的過電流，隨后又緩慢下降到正常電流水平，導(dǎo)致誤識(shí)別出負(fù)變點(diǎn)。

在搜索到各個(gè)q值對(duì)應(yīng)的Mq和Tqmin之后，接下來尋找正確變點(diǎn)數(shù)目。

以在SA算法得到的T1min-T5min以及T0min，顯然T0min≥T1min≥…≥T5min。由公式(2)可以看出，在變點(diǎn)數(shù)qk，則增加q并不會(huì)明顯增大Tqmin，因?yàn)橐呀?jīng)劃分好各穩(wěn)態(tài)區(qū)間，再對(duì)區(qū)間劃分顯然方差并不會(huì)有明顯變化。

這樣問題轉(zhuǎn)化為尋找Tqmin序列變化的拐點(diǎn)。這里采用一個(gè)比值K=Tqmin/Tq-1min，目標(biāo)是尋找到剛好不滿足K>K0的最小的變點(diǎn)數(shù)目q。即小于q時(shí)滿足該不等式，否則不滿足。K0=1.1，為比較閾值。

2.3 LS-HT算法

在根據(jù)上一個(gè)循環(huán)算法A能得到變點(diǎn)粗略的位置，對(duì)于復(fù)雜負(fù)荷開關(guān)事件往往會(huì)誤識(shí)別出“變點(diǎn)”。在連續(xù)變化型負(fù)荷運(yùn)行中，由于其功率連續(xù)變化，在算法A求解最優(yōu)變點(diǎn)數(shù)目時(shí)，也往往會(huì)將其誤識(shí)別。因此本文采用最小二乘-假設(shè)檢驗(yàn)法(Least Square-Hypothetical Test, LS-HT)對(duì)識(shí)別的變點(diǎn)進(jìn)行檢查。

對(duì)其劃分時(shí)間段，使得每段窗口數(shù)據(jù)中至多含有一個(gè)變點(diǎn)，且有足夠多的樣本數(shù)目。接下來進(jìn)行假設(shè)檢驗(yàn)，流程如表2所示。

表2 LS-ST

進(jìn)行假設(shè)檢驗(yàn)不僅可以對(duì)上一個(gè)循環(huán)的最優(yōu)化算法的結(jié)果進(jìn)行檢驗(yàn)，以防止多識(shí)別變點(diǎn)的情況；而且最優(yōu)化算法往往只能給出個(gè)在精確解附近的相似解，但后續(xù)的假設(shè)檢驗(yàn)針對(duì)于只有一個(gè)變點(diǎn)的區(qū)間，再次運(yùn)用最小二乘法可以準(zhǔn)確的確定變點(diǎn)的時(shí)刻。

3 算例仿真分析

3.1 非侵入式負(fù)荷測(cè)量系統(tǒng)

本文采用的監(jiān)測(cè)設(shè)備為廣州貫行電能技術(shù)有限公司的高級(jí)全信息采集單元CEIU-S-01，如圖2所示。CEIU為一款具備雙向通信的集成多功能電能計(jì)量?jī)x、電能質(zhì)量監(jiān)測(cè)儀、周波表、故障錄波等裝置全部功能的智能化終端裝置。將其通過網(wǎng)線連接至局域網(wǎng)的電腦，即可實(shí)時(shí)獲得測(cè)量的數(shù)據(jù)。

圖2 CEIU量測(cè)單元電氣圖

CEIU擁有7個(gè)通道接入測(cè)量數(shù)據(jù)，包括3相電壓，零序電壓，三相電流，輸入電源為市電220 V，本文采用其中一相電壓、電流作為實(shí)驗(yàn)測(cè)試。其電壓變比kU=5928256*515.0,電流變比kI=5928256*10.32。采樣頻率為8 kHz，即每個(gè)周波160個(gè)點(diǎn)。其采用基于UDP通信協(xié)議Modbus隧道的通信方式將電壓、電流波形數(shù)據(jù)向外發(fā)送。

本文通過圖3所示非侵入式負(fù)荷分解測(cè)量裝置獲得了實(shí)驗(yàn)所需的數(shù)據(jù)。測(cè)量裝置右邊通過插頭接入電源，電源通過端子排分流，一路供給CEIU用電，另一路通過電流互感器，供給負(fù)荷用電。電流互感器二次側(cè)出線和連接負(fù)荷的插座引出線，作為測(cè)量的電流和電壓信息接入CEIU，CEIU最后通過網(wǎng)線將測(cè)量到的數(shù)據(jù)發(fā)送送入電腦，以供后續(xù)變點(diǎn)識(shí)別分析。

圖3 非侵入式負(fù)荷測(cè)量系統(tǒng)

本文選取幾種常用電器(空調(diào)，微波爐，筆記本電腦，燈管，熱水壺等)進(jìn)行30次測(cè)量，采樣周期T=0.02 s,每次每種家電運(yùn)行時(shí)間4 min得到各電器運(yùn)行曲線如圖4所示。

圖4 實(shí)驗(yàn)家電運(yùn)行曲線

由有功功率運(yùn)行曲線可以看出，負(fù)變點(diǎn)較好識(shí)別，負(fù)荷關(guān)閉時(shí)，不像啟動(dòng)過程，總有功功率總是瞬間減小，可以視作無暫態(tài)過程。然而在負(fù)荷啟動(dòng)時(shí)，則情況較為復(fù)雜，根據(jù)其負(fù)荷啟動(dòng)時(shí)暫態(tài)功率的情況分為三類：Ⅰ，Ⅱ，Ⅲ型負(fù)載。各家電負(fù)荷特征如表3所示。

表3 實(shí)驗(yàn)家電特征庫(kù)

1)Ⅰ型負(fù)載：這種負(fù)載以電阻性負(fù)載為主，其變點(diǎn)識(shí)別比較簡(jiǎn)單，跟一般的用電負(fù)荷關(guān)閉時(shí)類似，其功率變化很快，可視作無暫態(tài)過程。

2)Ⅱ型負(fù)載：這類大部分涵蓋了電動(dòng)機(jī)類型的負(fù)荷，其啟動(dòng)過程有明顯的暫態(tài)過程，在啟動(dòng)瞬間會(huì)產(chǎn)生較大的過電流。以空調(diào)為例，空調(diào)壓縮機(jī)在啟動(dòng)過程中有較大的過電流，隨后立馬下降，然后其功率P再緩慢上升直至穩(wěn)定。

由空調(diào)啟動(dòng)曲線可以看出，在空調(diào)啟動(dòng)時(shí)會(huì)有很大的過電流(功率)產(chǎn)生，而且如果只是檢測(cè)功率變化ΔP，則會(huì)誤識(shí)別出很多個(gè)變點(diǎn)，從而將空調(diào)啟動(dòng)曲線誤認(rèn)為是多個(gè)負(fù)荷啟動(dòng)事件。然而通過尋找最小二乘法的最優(yōu)解并進(jìn)行檢驗(yàn)，可以排除不正確的變點(diǎn)。

3)Ⅲ型負(fù)載：Ⅲ型負(fù)載啟動(dòng)和運(yùn)行曲線比較復(fù)雜，其變化較為劇烈。本文以筆記本電腦作研究，其在運(yùn)行時(shí)功率會(huì)隨著運(yùn)行模式變化其ΔP變化也不相同，簡(jiǎn)單的以一段時(shí)間的平均功率變化難以識(shí)別這種情況。而基于改進(jìn)最小二乘法的變點(diǎn)識(shí)別方法，則關(guān)注于一段時(shí)間內(nèi)的總誤差和分段誤差的差的大小，從而避免了只考慮均值所帶來的偏差，就能正確識(shí)別其正確的變點(diǎn)時(shí)刻。

3.2 算例分析

為了評(píng)價(jià)本文使用事件檢測(cè)方法，采用F-score評(píng)估度量[13]。TP(True Positive)是實(shí)際發(fā)生事件也被正確識(shí)別的情況，F(xiàn)P(False Positive)是實(shí)際存在變點(diǎn)但未被識(shí)別，F(xiàn)N(False Negative)表示實(shí)際無負(fù)荷事件發(fā)生但被識(shí)別為變點(diǎn)。

接下來計(jì)算精度指標(biāo)[12]：PR代表的是提取出的變點(diǎn)數(shù)中正確信息條數(shù)所占的比例，RE值是所有樣本的信息條數(shù)中所占的比例。在評(píng)價(jià)中，P值和R值都是越高越好，但事實(shí)上這兩者在某些情況下是矛盾的，因此，就需要綜合考慮二者，F(xiàn)M值就是兩者的綜合指標(biāo)。定義：

PR=TP/(TP+FP)

(10)

RE=TP/(TP+FN)

(11)

FM=2×(PR×RE)/(PR+RE)

(12)

本文對(duì)實(shí)驗(yàn)家電各進(jìn)行30次開關(guān)測(cè)量，即TP+FP=60，選取采樣周期T=0.02 s，窗口長(zhǎng)度winlen=600，即12 s長(zhǎng)度的時(shí)間，進(jìn)行驗(yàn)證。并傳入數(shù)據(jù)對(duì)算法進(jìn)行驗(yàn)證，得到的結(jié)果如表4所示：

表4 某商業(yè)寫字樓負(fù)荷識(shí)別結(jié)果

從表4來看，I型負(fù)荷最為簡(jiǎn)單，其變點(diǎn)基本上能識(shí)別出來；而II型負(fù)荷由于暫態(tài)過電流的存在，使得FN≠0，即誤識(shí)別出一些“變點(diǎn)”；而在III型負(fù)荷中，雖然其負(fù)荷變化較為劇烈，但本文提出的算法識(shí)別率還是能達(dá)到70%以上。

4 結(jié)論

本文提出了一種基于均值變點(diǎn)模型的識(shí)別方法，通過滑動(dòng)窗口，利用最小二乘法計(jì)算相關(guān)目標(biāo)函數(shù)，以確定變點(diǎn)的個(gè)數(shù)。此算法簡(jiǎn)單、抗干擾能力強(qiáng)。通過實(shí)際數(shù)據(jù)實(shí)驗(yàn)表明它能夠自動(dòng)檢測(cè)、辨識(shí)負(fù)荷開關(guān)而產(chǎn)生電氣量變點(diǎn)，從而為NILM算法提供準(zhǔn)確的負(fù)荷事件識(shí)別，為之后負(fù)荷特征的提取和識(shí)別打下堅(jiān)實(shí)的基礎(chǔ)。