基于計(jì)算智能的竊電檢測(cè)模型研究

林華城,葉泳泰,陳錦迅,賴佛強(qiáng),陸建巧

(廣東電網(wǎng)有限責(zé)任公司惠州供電局,廣東 惠州 516000)

0 引言

隨著網(wǎng)絡(luò)、大數(shù)據(jù)、物聯(lián)網(wǎng)、通信技術(shù)[1-2]的不斷發(fā)展,電力系統(tǒng)每時(shí)每刻都在不停地接收新信息。同時(shí),電力公司可通過(guò)高級(jí)量測(cè)體系(advanced metering infrastructure,AMI)[3]監(jiān)控細(xì)粒度時(shí)間間隔內(nèi)的電力消耗,從而更有效地監(jiān)控電力系統(tǒng)中的異常情況。

一般情況下,電力系統(tǒng)中存在2種類(lèi)型的損耗[4],分別是技術(shù)損耗和非技術(shù)損耗。技術(shù)損耗由線路、變壓器和其他設(shè)備中的電阻元件加熱引起。非技術(shù)損耗主要由竊電、電表故障或計(jì)費(fèi)錯(cuò)誤引起。為此,可通過(guò)采集到的大數(shù)據(jù)分析用戶的用電行為,建立異常用電行為檢測(cè)模型。這樣不僅可以減少竊電行為,還可為電能表故障檢測(cè)提供新的解決思路。

目前,AMI計(jì)量數(shù)據(jù)在竊電檢測(cè)中的應(yīng)用方案主要有兩種,分別是基于數(shù)據(jù)統(tǒng)計(jì)[5-6]的方法和基于機(jī)器學(xué)習(xí)[7-8]的方法。基于數(shù)據(jù)統(tǒng)計(jì)的方法實(shí)施相對(duì)簡(jiǎn)單,但只能判斷低壓站發(fā)生了竊電,無(wú)法準(zhǔn)確定位非法用戶。如果要確定可疑用戶,則需要手動(dòng)逐個(gè)檢查區(qū)域內(nèi)的所有用戶。這樣做效率較低,且對(duì)檢測(cè)過(guò)程中檢測(cè)人員的素質(zhì)提出了更高要求。文獻(xiàn)[9]提出了1種基于差分進(jìn)化支持向量機(jī)(support vector machine,SVM)的識(shí)別方法。文獻(xiàn)[10]提出了1種基于貝葉斯優(yōu)化和改進(jìn)XGBoost模型的竊電檢測(cè)方法。然而,機(jī)器學(xué)習(xí)方法也存在一定缺陷,如訓(xùn)練超參數(shù)確定困難、算法易陷入局部最優(yōu)等。隨著深度學(xué)習(xí)技術(shù)的發(fā)展,部分學(xué)者將深度學(xué)習(xí)方法引入竊電檢測(cè),從而獲得更高的準(zhǔn)確度。然而,深度學(xué)習(xí)方法訓(xùn)練復(fù)雜,且非常依賴樣本數(shù)據(jù)。同時(shí),如果無(wú)法有效處理AMI不平衡樣本集問(wèn)題,將會(huì)嚴(yán)重影響訓(xùn)練效果。

為改善上述問(wèn)題,本文提出了1種將無(wú)監(jiān)督學(xué)習(xí)和有監(jiān)督學(xué)習(xí)相結(jié)合的竊電檢測(cè)模型。本文首先基于Wasserstein距離、相似性約束和真實(shí)約束,訓(xùn)練生成對(duì)抗網(wǎng)絡(luò)(generative adversarial network,GAN),以生成符合竊電特征的高精度測(cè)量數(shù)據(jù);然后,為了提高相似性度量的有效性,在相似性度量過(guò)程中綜合考慮了數(shù)值特征和形態(tài)特征,并采用動(dòng)態(tài)時(shí)間扭曲(dynamic time warping,DTW)來(lái)度量形態(tài)特征的相似性;最后,使用SVM-K近鄰(K-nearest neighbor,KNN)進(jìn)行竊電檢測(cè)。

1 基于計(jì)算智能的竊電檢測(cè)模型

基于計(jì)算智能的竊電檢測(cè)模型結(jié)構(gòu)如圖1所示。

圖1 基于計(jì)算智能的竊電檢測(cè)模型結(jié)構(gòu)圖

典型的低壓配電網(wǎng)由多個(gè)相互連接的單元組成。電力通過(guò)高壓線從發(fā)電廠輸送到變電站,并從變電站輸送到工業(yè)、商業(yè)和住宅區(qū)。在此過(guò)程中,通過(guò)AMI可測(cè)量、收集、存儲(chǔ)、分析和使用客戶數(shù)據(jù)。收集的數(shù)據(jù)包括各種大、中、小型典型變壓器用戶以及380 V、220 V低壓居民用戶的數(shù)據(jù)。收集的信息包括數(shù)據(jù)項(xiàng),如電能數(shù)據(jù)、事件記錄和其他數(shù)據(jù)。通過(guò)對(duì)收集的數(shù)據(jù)及信息進(jìn)行分析,可以獲取用戶的電力消費(fèi)信息和消費(fèi)行為信息。在低壓變電站中,能量損失是基于能量平衡失配產(chǎn)生的。假設(shè)所有電表讀數(shù)正常,而低壓站被視為1個(gè)節(jié)點(diǎn),則根據(jù)基爾霍夫定律,子電表讀數(shù)與網(wǎng)損之和等于總電表讀數(shù)。因此,當(dāng)網(wǎng)損過(guò)大時(shí),可認(rèn)為竊電的可能性非常高。

基于計(jì)算智能的竊電檢測(cè)模型執(zhí)行過(guò)程包含3個(gè)關(guān)鍵過(guò)程。

①確定可疑站點(diǎn)。竊電檢測(cè)模型基于改進(jìn)的模糊C均值(improved fuzzy C-means,IFCM)聚類(lèi)算法對(duì)用戶的歷史數(shù)據(jù)進(jìn)行聚類(lèi),從而獲取用戶的用電特性曲線。

②基于相似性約束和真實(shí)性約束,使用GAN生成符合竊電特征的高精度測(cè)量數(shù)據(jù)。

③竊電用戶定位。綜合考慮待測(cè)曲線的數(shù)字特征、形態(tài)特征以及特征曲線,將平衡數(shù)據(jù)集代入改進(jìn)的SVM-KNN模型進(jìn)行訓(xùn)練,從而準(zhǔn)確識(shí)別非法竊電用戶。

2 可疑站點(diǎn)檢測(cè)和客戶特征分析

2.1 站點(diǎn)檢測(cè)方法

本文假定AMI系統(tǒng)獲得的竊電數(shù)據(jù)涉及時(shí)間n內(nèi)的m個(gè)用戶,其數(shù)據(jù)形式可由矩陣描述。同時(shí),本文令同一用戶在不同時(shí)期的數(shù)據(jù)為xj。不同用戶在時(shí)間n的數(shù)據(jù)向量X可描述如下。

(1)

式中:xnm為第n個(gè)用戶在第m個(gè)測(cè)量期間由智能儀表測(cè)量的值。

為消除數(shù)據(jù)多尺度影響,本文基于最大-最小歸一化函數(shù)對(duì)數(shù)據(jù)進(jìn)行處理,并將數(shù)據(jù)映射為[0,1]。歸一化函數(shù)如式(2)所示。

(2)

式中:x為實(shí)際測(cè)量數(shù)據(jù);xmax為樣本數(shù)據(jù)的最大值;xmin為樣本數(shù)據(jù)的最小值;x*為歸一化后的用電量數(shù)據(jù)。

2.2 特征選取

用戶的用電行為多種多樣。對(duì)于不同的用戶,通常選擇不同的特征集來(lái)分析用戶的用電行為具有不同的分析結(jié)果。然而,特征空間中包含的冗余信息導(dǎo)致分析結(jié)果效果不佳。因此,有必要去除重疊和冗余信息,從而提高分析性能。

本文選取的基本用戶用電行為特征包括統(tǒng)計(jì)特征、時(shí)間序列特征和關(guān)系特征。其中:統(tǒng)計(jì)特征包括日用電量數(shù)據(jù)、年用電量數(shù)據(jù)、季節(jié)用電量數(shù)據(jù)、日最大和最小負(fù)荷、平均負(fù)荷率等;時(shí)間序列特征包括高峰小時(shí)耗電率、谷功率系數(shù)等;關(guān)系特征包括房屋面積和家庭成員人數(shù)等。為去除冗余特征,本文提出了IFCM特征選擇算法。該算法具體過(guò)程如下。

①選擇特征數(shù)為1,確定特征個(gè)數(shù)為1時(shí)的最高聚類(lèi)評(píng)價(jià)標(biāo)準(zhǔn)。

②選擇特征數(shù)為2,根據(jù)所選特征選擇新特征,確定具有最高聚類(lèi)評(píng)估標(biāo)準(zhǔn)的特征。

③選擇特征數(shù)為i,根據(jù)所選特征選擇(i+1)個(gè)特征。其中,(i+1)特征為具有最高的聚類(lèi)評(píng)價(jià)標(biāo)準(zhǔn)的特征。

④重復(fù)上述步驟,直到選擇(n+1)個(gè)特征,且具有最高聚類(lèi)評(píng)估標(biāo)準(zhǔn)的特征為所選特征。

3 基于GAN的數(shù)據(jù)增強(qiáng)

本文使用GAN[11]生成竊電數(shù)據(jù)。由于測(cè)量數(shù)據(jù)是一維時(shí)間序列,因此本文設(shè)計(jì)了基于一維卷積層的GAN結(jié)構(gòu)。同時(shí),本文基于Wasserstein距離、相似性約束和真實(shí)約束,生成符合竊電特征的高精度測(cè)量數(shù)據(jù)。所生成的樣本與現(xiàn)有樣本相結(jié)合,可以獲得大量樣本。基于GAN的數(shù)據(jù)增強(qiáng)結(jié)構(gòu)總體框架如圖2所示。

模型首先選擇現(xiàn)有的少量竊電數(shù)據(jù)作為訓(xùn)練集。本文令竊電數(shù)據(jù)為pd(X),數(shù)據(jù)中存在的1組隨機(jī)變量z滿足高斯分布pz(z)。GAN可從已知分布中采樣,生成滿足原始數(shù)據(jù)分布的新數(shù)據(jù)。訓(xùn)練過(guò)程中,生成器G負(fù)責(zé)學(xué)習(xí)樣本分布的規(guī)律并生成新樣本。G由神經(jīng)網(wǎng)絡(luò)組成。其輸入為先驗(yàn)分布pz(對(duì)應(yīng)z),輸出為G(z)。需注意,生成數(shù)據(jù)的目標(biāo)是生成盡可能真實(shí)的數(shù)據(jù),即生成數(shù)據(jù)的分布規(guī)律pg(z)與樣本數(shù)據(jù)pd(X)擬合。如生成器損失函數(shù)為Ez～pz{-D[G(z)]},則目標(biāo)函數(shù)fD定義如下。

fD=minEz～pz{-D[G(z)]}

(3)

式中:D為鑒別器,負(fù)責(zé)確定輸入數(shù)據(jù)是否真實(shí)。

D是1個(gè)神經(jīng)網(wǎng)絡(luò),但其輸入是實(shí)際數(shù)據(jù)或生成器生成的數(shù)據(jù)。鑒別器的主要任務(wù)是區(qū)分2種數(shù)據(jù),因此其輸出是1個(gè)介于0和1之間的標(biāo)量,即屬于實(shí)際數(shù)據(jù)或生成數(shù)據(jù)的概率。D的損失函數(shù)可以定義為Ex～pd[-D(x)]+Ez～pz{-D[G(z)]}。D的目標(biāo)函數(shù)fE定義為:

fE=maxEx-pd[-D(x)]+Ez～pz{-D[G(z)]}

(4)

因此,整個(gè)對(duì)抗過(guò)程的目標(biāo)函數(shù)定義為:

Ez～pz{log{1-D[G(z)]}}

(5)

本文選取最小Wasserstein距離為目標(biāo)訓(xùn)練GAN,從而有效提高GAN訓(xùn)練的穩(wěn)定性。Wasserstein距離定義如下。

(6)

式中:∏(pd,pg)為聯(lián)合分布γ的集合,pd和pg分別為鑒別器和生成器中數(shù)據(jù)的邊際分布;W(pd,pg)為γ(x,y)期望的下確界,即pg到pd的擬合需要x到y(tǒng)的距離。

由于很難直接計(jì)算任意分布之間的Wasserstein距離,因此本文采用對(duì)偶形式:

(7)

經(jīng)過(guò)訓(xùn)練后,生成器可以生成大量樣本,且滿足分布要求。為了確保生成的測(cè)量數(shù)據(jù)的真實(shí)性,必須同時(shí)滿足真實(shí)性和相似性約束。真實(shí)性約束用于確保生成的數(shù)據(jù)能夠接近真實(shí)情況。真實(shí)性損失Lr為:

Lr=W[G(z;θ(G));θ(D)]

(8)

式中:G(z;θ)為生成器的生成數(shù)據(jù);W(～;θ(D))為生成數(shù)據(jù)與實(shí)際樣本之間的Wasserstein距離。

同理,生成的數(shù)據(jù)應(yīng)盡可能與實(shí)際數(shù)據(jù)相似。

因此,相似性損失Ls為:

(9)

(10)

(11)

4 竊電用戶定位

本文提出的竊電用戶定位包括2個(gè)關(guān)鍵環(huán)節(jié),分別為基于相似性度量確定可疑用戶,以及基于改進(jìn)SVM初步檢測(cè)可疑用戶并輸出非法用戶。竊電用戶定位執(zhí)行流程如圖3所示。

圖3 竊電用戶定位執(zhí)行流程圖

圖3中,D1和D2為設(shè)置的閾值,且D1D2時(shí),用戶大概率為可疑用戶,應(yīng)手動(dòng)驗(yàn)證。當(dāng)DW

4.1 相似性度量

時(shí)間序列的相似性包括值和形態(tài)2個(gè)方面。為了獲得特征曲線的值和形態(tài),本文分別使用歐幾里德距離、DTW度量值和形態(tài)特征相似性。為了簡(jiǎn)單而準(zhǔn)確地描述曲線的形態(tài)特征,例如不同時(shí)期的上升、下降和穩(wěn)定性,本文使用直線斜率表示該時(shí)期的形態(tài)特征。因此,長(zhǎng)度為n的時(shí)間序列可縮減為形態(tài)序列(n-1)。

本文令2個(gè)獨(dú)立的時(shí)間序列分別為X=(x1,x2,…,xn-1)和Y=(y1,y2,…,ym-1)。因此,X和Y的值相似性可由距離矩陣度量。其中,矩陣中的每個(gè)元素由歐幾里德距離表示。

(12)

式中:i=1,2,…,n;j=1,2,…,m。

同理,X和Y的形態(tài)相似性可由DTW度量。然而,DTW不是隨機(jī)選擇,其需要滿足邊界條件、連續(xù)性和單調(diào)性約束。在滿足這3個(gè)約束后,可得到許多路徑。本文需要選擇1條使最終獲得的總距離最小化的路徑。為此,本文基于動(dòng)態(tài)規(guī)劃方法構(gòu)造累積距離γ。累積距離γ(i,j)表示當(dāng)前網(wǎng)格點(diǎn)的距離D(i,j)與可到達(dá)該點(diǎn)的最小相鄰元素的累積距離之和。

γ(i,j)=D(q,c)+min[γ(i-1,j-1)

γ(i-1,j)γ(i,j-1)]

(13)

4.2 改進(jìn)的SVM-KNN模型

SVM-KNN模型執(zhí)行流程如圖4所示。

圖4 SVM-KNN模型執(zhí)行流程圖

SVM是1種典型的有監(jiān)督學(xué)習(xí)方法,主要用于分析給定數(shù)據(jù)并識(shí)別輸入值相對(duì)于輸出的模式或趨勢(shì)。傳統(tǒng)SVM為1個(gè)層次化的分類(lèi)模型,其最大的問(wèn)題是錯(cuò)誤積累。這影響了分類(lèi)的準(zhǔn)確性。如果使用有偏二叉樹(shù)進(jìn)行分類(lèi),則需要構(gòu)造一個(gè)誤差積累小、分類(lèi)精度高的決策樹(shù)。為了減少誤差累積的影響,本文采用投影向量方法來(lái)度量類(lèi)間的分離度,并在此基礎(chǔ)上構(gòu)造了1個(gè)有偏二項(xiàng)決策樹(shù)。同時(shí),由于數(shù)據(jù)遠(yuǎn)離超平面,SVM能夠準(zhǔn)確地進(jìn)行分類(lèi)。但當(dāng)距離接近超平面時(shí),分類(lèi)效果較低,在超平面附近容易發(fā)生誤分類(lèi)。為提高竊電用戶分類(lèi)精度,本文提出了1種改進(jìn)的SVM-KNN模型。SVM-KNN模型利用界面附近樣本提供的信息來(lái)提高分類(lèi)精度。

對(duì)識(shí)別樣本進(jìn)行分類(lèi)時(shí),需計(jì)算樣本與分類(lèi)超平面之間的距離。如果距離大于給定的閾值,則直接應(yīng)用SVM分類(lèi);否則,應(yīng)用KNN分類(lèi)。在KNN分類(lèi)中,每一類(lèi)的支持向量用來(lái)計(jì)算識(shí)別樣本與每個(gè)SVM之間的距離。

5 仿真與分析

5.1 數(shù)據(jù)集

仿真試驗(yàn)數(shù)據(jù)使用某電力公司提供的2019年至2020年約5 000個(gè)家庭用戶和企業(yè)用戶的半小時(shí)用電報(bào)告。數(shù)據(jù)經(jīng)用戶同意,主要從用戶家或辦公室中安裝的智能電表中獲取。每個(gè)用戶數(shù)據(jù),至少都有1個(gè)包含350天的半小時(shí)用電報(bào)告信息。為了不失一般性,本文假設(shè)所有樣本都屬于誠(chéng)實(shí)用戶。根據(jù)竊電的實(shí)際情況,本文建立了6種類(lèi)別的竊電情況。第一類(lèi)竊電是所有樣本乘以相同的隨機(jī)選擇系數(shù)。第二類(lèi)竊電是1種典型的“開(kāi)-關(guān)”攻擊,即在某些時(shí)間間隔內(nèi),耗電量報(bào)告為零。第三類(lèi)竊電是將耗電量乘以隨時(shí)間變化的隨機(jī)系數(shù)。第四類(lèi)竊電是第二類(lèi)和第三類(lèi)的結(jié)合。第五類(lèi)竊電是在高峰時(shí)段乘以相同的隨機(jī)選擇系數(shù)。第六類(lèi)竊電是1種隨機(jī)周期的“開(kāi)-關(guān)”攻擊,但持續(xù)時(shí)間短且不連續(xù),因此減少了總用電量。與第二類(lèi)竊電相比,由于時(shí)間段的隨機(jī)性,第六類(lèi)竊電檢測(cè)顯得更加困難。

5.2 站點(diǎn)檢測(cè)性能對(duì)比分析

5.2.1 特征選取

本節(jié)對(duì)用戶用電行為的常用特征選取策略進(jìn)行對(duì)比。表1所示為特征數(shù)量和準(zhǔn)確率變化統(tǒng)計(jì)結(jié)果。由表1可知:隨著特征數(shù)量的增加,聚類(lèi)的準(zhǔn)確率增加;當(dāng)特征數(shù)量超過(guò)4時(shí),聚類(lèi)的準(zhǔn)確率降低。因此,研究最終確定代表電力消費(fèi)行為的特征數(shù)量為4。特征指標(biāo)為負(fù)荷率、谷系數(shù)、高峰小時(shí)用電率和正常時(shí)段的用電量百分比。

表1 特征數(shù)量和準(zhǔn)確率變化統(tǒng)計(jì)結(jié)果

5.2.2 站點(diǎn)檢測(cè)

圖5所示為IFCM和傳統(tǒng)模糊C均值(fuzzy C-means,FCM)聚類(lèi)目標(biāo)函數(shù)對(duì)比結(jié)果。

圖5 IFCM和傳統(tǒng)FCM聚類(lèi)目標(biāo)函數(shù)對(duì)比結(jié)果

由圖5可知,IFCM迭代次數(shù)較少,且每個(gè)點(diǎn)到聚類(lèi)中心的代數(shù)和較少。因此與FCM相比,IFCM可以有效提高算法的分類(lèi)效果和迭代時(shí)間。

5.3 數(shù)據(jù)生成對(duì)比分析

考慮到正常用戶和非正常用戶之間數(shù)據(jù)的不平衡,本節(jié)基于一維卷積層的GAN生成竊電數(shù)據(jù)。為了驗(yàn)證所提SVM-KNN模型的性能,本節(jié)比較了合成少數(shù)類(lèi)過(guò)采樣技術(shù)(synthetic minonity over-sampling technique,SMOTE)和分界線-SMOTE(borderline-SMOTE,B-SMOTE)數(shù)據(jù)增強(qiáng)算法的分類(lèi)性能。不同數(shù)據(jù)增強(qiáng)算法分類(lèi)性能對(duì)比結(jié)果如表2所示。

由表2可知,無(wú)論有無(wú)噪聲,SVM-KNN生成的樣本都能有效提高分類(lèi)器的分類(lèi)精度。

5.4 竊電檢測(cè)性能對(duì)比分析

本節(jié)將SVM-KNN與傳統(tǒng)SVM、決策樹(shù)-SVM(decision tree-SVM,DT-SVM)、KNN、卷積神經(jīng)網(wǎng)絡(luò)(convolutional neural network,CNN)等模型在不同噪聲數(shù)據(jù)中進(jìn)行綜合對(duì)比。表3所示為不同方法綜合性能對(duì)比結(jié)果。對(duì)比指標(biāo)包括識(shí)別準(zhǔn)確率、召回率和F分?jǐn)?shù)。

表3 不同方法綜合性能對(duì)比結(jié)果

由表3可知,各模型在無(wú)噪聲數(shù)據(jù)集中的指標(biāo)均優(yōu)于有噪聲數(shù)據(jù)集;隨著數(shù)據(jù)集中的噪聲數(shù)據(jù)增加,各指標(biāo)均有不同程度下降。綜合對(duì)比后可知,所提模型具有較高的魯棒性,在驗(yàn)證干擾數(shù)據(jù)集中表現(xiàn)較為優(yōu)異。對(duì)比結(jié)果進(jìn)一步驗(yàn)證了所提模型對(duì)電力行業(yè)竊電行為檢測(cè)具有較高的準(zhǔn)確性和穩(wěn)定的識(shí)別率。

6 結(jié)論

本文對(duì)電力行業(yè)竊電檢測(cè)進(jìn)行了研究與分析,建立了1種基于計(jì)算智能的竊電檢測(cè)模型。首先,本文基于IFCM對(duì)用戶的歷史數(shù)據(jù)進(jìn)行聚類(lèi),從而獲取用戶的用電特性曲線。其次,本文基于相似性約束和真實(shí)性約束,使用GAN生成符合竊電特征的高精度測(cè)量數(shù)據(jù)。最后,本文綜合考慮待測(cè)曲線的數(shù)字特征、形態(tài)特征以及特征曲線,將平衡數(shù)據(jù)集代入改進(jìn)的SVM-KNN模型進(jìn)行訓(xùn)練,從而準(zhǔn)確識(shí)別非法竊電用戶。該模型為電力部門(mén)分析用戶用電行為以及竊電檢測(cè)提供了借鑒。