汪世平，鄒志楊，陳洪濤，蔡 慧

（1.國電南瑞科技股份有限公司，江蘇 南京 211106；2.中國計量大學(xué)機電工程學(xué)院，浙江 杭州 310018）

0 引言

智能電能表作為供電量與耗電量的計量記錄裝置，其運行狀況不但直接影響電網(wǎng)公司的發(fā)展和經(jīng)營效益，還影響用戶用電的公平性與安全性［1?3］。目前在發(fā)電、輸電、配電、用電等重要的環(huán)節(jié)中，智能電能表得到了大面積的應(yīng)用。智能電能表作為智能電網(wǎng)中的關(guān)鍵部分，在電網(wǎng)公司進行實時計量、控制負荷、響應(yīng)電力需求等方面起到重要作用［4］。

電網(wǎng)公司大力推廣大數(shù)據(jù)技術(shù)在智能電網(wǎng)的全面建設(shè)，各種新興技術(shù)在電力行業(yè)領(lǐng)域得到廣泛應(yīng)用，其中變電站數(shù)據(jù)、配電數(shù)據(jù)、用電數(shù)據(jù)、營銷數(shù)據(jù)等數(shù)據(jù)庫的建立也日趨完善，為大數(shù)據(jù)技術(shù)研發(fā)提供了良好的研究環(huán)境。

國家電網(wǎng)有限公司在2018 年1 月頒發(fā)了《國家電網(wǎng)公司關(guān)于營銷大數(shù)據(jù)應(yīng)用工作的指導(dǎo)意見》，要推動電力營銷服務(wù)領(lǐng)域的業(yè)務(wù)創(chuàng)新，大數(shù)據(jù)技術(shù)研究是其中的重要戰(zhàn)略步驟［5］。目前國家電網(wǎng)有限公司完成了相應(yīng)數(shù)據(jù)管理平臺的建設(shè)和投運，逐步發(fā)揮大數(shù)據(jù)技術(shù)在用電數(shù)據(jù)管理分析方面的作用［6］。但是，目前的數(shù)據(jù)分析技術(shù)還不能充分地利用智能電能表數(shù)據(jù)達到理想的分析效果，仍需要不斷地進行研究實踐。因此采用大數(shù)據(jù)技術(shù)對電能表數(shù)據(jù)進行分析是未來智能電網(wǎng)技術(shù)的發(fā)展趨勢，而目前日益完善的大數(shù)據(jù)平臺及相應(yīng)技術(shù)使其在發(fā)展上具有很大的潛力。目前電力公司的主要研究方向之一就是如何利用大數(shù)據(jù)技術(shù)對智能電能表的用電數(shù)據(jù)進行分析，從而更好地對智能電能表狀態(tài)進行評估。

文獻［7］針對現(xiàn)行計量裝置狀態(tài)評估指標的不確定性和各評估方法結(jié)果的不一致性，將電能計量裝置狀態(tài)評估的各類數(shù)據(jù)、評估指標體系及評估方法組件化，動態(tài)搭建電能計量裝置狀態(tài)評估系統(tǒng)，狀態(tài)指標描述較為全面，但未涉及電能計量裝置各組成部分狀態(tài)指標的細節(jié)差異；文獻［8］中所述的基于灰色關(guān)聯(lián)度的智能電能表供應(yīng)商多目標綜合評價方法，在保留多目標綜合評價模型優(yōu)點的同時解決了傳統(tǒng)智能電表供應(yīng)商綜合評價中指標煩瑣，主觀性強的問題。但該方法只是從整體的電能表批次進行分析，忽略了電能表個體的狀態(tài)評價。

結(jié)合變異系數(shù)賦值法和K?means 聚類算法提出了新的智能電能表狀態(tài)評價方法，通過地區(qū)影響指數(shù)、全事件、計量異常、電能表過載和時鐘電池欠壓6個指標進行分析，得到智能電能表的評價結(jié)果，并根據(jù)結(jié)果分析方法的合理性。

1 基于變異系數(shù)和K?means的狀態(tài)評價方法

1.1 設(shè)計思路

首先根據(jù)前期統(tǒng)計分析研究得到與智能電能表狀態(tài)評價影響較大的6 項指標，即電能表可靠度、地區(qū)因素、全事件、計量異常、電能表過載和時鐘電池欠壓［9］，然后通過兩種方法對這些指標進行分析，一種是通過變異系數(shù)法加權(quán)處理，通過不同指標對于電能表運行狀態(tài)的影響情況，賦予各個指標相應(yīng)的權(quán)重，最后以評分的形式反映每個電能表的運行狀態(tài)；另一種是將這些指標作為電能表的特征，將電能表視為空間中的點，這些指標即為點的坐標，采用聚類方法將這些點分類，從而得到不同的電能表評價狀態(tài)。將這兩種方法的評價結(jié)果相結(jié)合，最后得到最終的電能表運行狀態(tài)評價結(jié)果。流程如圖1所示。

圖1 狀態(tài)評價方法設(shè)計流程

1.2 指標數(shù)據(jù)

1）電能表可靠度。

電能表可靠度計算公式為

式中：Mr為電能表可靠度指標；f(i)為當前批次電能表中第i個月故障表數(shù)量；N為當前批次電能表的總數(shù)量；t為當前月份。

2）地區(qū)因素。

地區(qū)因素指標計算公式為

式中：Mt為地區(qū)因素指標；Hx為第x個市安裝的電能表數(shù)量；H為安裝電能表總數(shù)量；Jx為第x個市的故障電能表數(shù)量；J為故障電能表總數(shù)量。

3）全事件。

全事件指標計算公式為

式中：Mq為全事件指標；A2為電能表檢定結(jié)果為故障；Ci(i=1，2，…，11)為電能表停走發(fā)生、電能表飛走發(fā)生、電能表超差發(fā)生、電能表反接、電能表時間超差發(fā)生、電能表停電次數(shù)變化發(fā)生、電能表斷相次數(shù)發(fā)生、磁場異常發(fā)生、電能表互感器倍率更改、開蓋記錄和對時這11 種事件；P(A2|Ci)表示在Ci條件下A2發(fā)生的概率。

4）計量異常。

計量異常指標計算公式為

式中：Ma為計量異常指標；Bj(j=1，2，…，6)為電能表示值不平、電能表飛走、電能表倒走、電能表停走、反向電量異常和時鐘異常；yk(k=1，2，…，5)為電壓越限、電壓失壓、電流過流、電壓斷相和潮流反向事件與異常關(guān)聯(lián)度；P(A2|Bj)表示在Bj條件下A2發(fā)生的概率。

5）電能表過載。

電能表過載指標計算公式為

式中：Ml為電能表過載指標；WN為電能表正常額定工作下24 h 計量的電量；KW為6 個月內(nèi)日用電量超過其電能表標準計量電量的天數(shù)比例；WO為6 個月內(nèi)日用電量超過電能表標準計量電量部分的平均日用電量。

6）時鐘電池欠壓。

時鐘電池欠壓指標計算公式為

式中：Mc為時鐘電池欠壓指標；z為6 個月內(nèi)發(fā)生的時鐘欠壓次數(shù)。

1.3 變異系數(shù)法原理

變異系數(shù)賦值法是直接利用各項指標所包含的信息，通過計算得到指標的權(quán)重，是一種客觀賦權(quán)的方法。根據(jù)指標數(shù)據(jù)變化對于電能表評價結(jié)果的影響，進一步分析該指標對于結(jié)果評價的重要程度。體現(xiàn)在數(shù)值上就是當該指標數(shù)據(jù)的變異程度越大，對該指標的賦值也越大［10?11］。

本文選擇的變異系數(shù)為標準差，其主要的計算步驟如下：

1）假設(shè)有m個待評估對象，一共有n個評估指標，指標的評估矩陣X可表示為

式中：xab為第b(b=1，2…，m)個評估對 象的第a(a=1，2，…，n)個評估指標的特征數(shù)據(jù)。

2）計算第a個指標的平均值為

3）計算第a個指標的標準差σa為

4）計算第a個指標的變異系數(shù)va為

5）將求得的變異系數(shù)進行歸一化處理，計算第a個指標的客觀權(quán)重βa為

1.4 K?means聚類原理

K?means 算法的核心思路是先隨機選取數(shù)據(jù)集中的數(shù)據(jù)作為初始聚類中心，并計算其他數(shù)據(jù)到這幾個數(shù)據(jù)點的距離，將這幾個數(shù)據(jù)劃分到距離它們最近的聚類中心，在遍歷完所有數(shù)據(jù)后，以各個類數(shù)據(jù)的平均值作為新的聚類中心，并再次重復(fù)上述操作，直到滿足某個閾值或達到預(yù)先設(shè)定的迭代次數(shù)后停止［12?15］。

K?means 算法的具體步驟主要分為如下幾個步驟：

1）根據(jù)已知的數(shù)據(jù)集，以k個數(shù)據(jù)作為初始的聚類中心C，其中這k個數(shù)據(jù)使任意選取的；

2）計算得到除聚類中心以外的數(shù)據(jù)樣本與聚類中心的歐式距離；

3）以歐式距離作為判斷依據(jù)，將數(shù)據(jù)樣本劃分到與其距離最近的聚類中心所屬的簇中；

4）計算得到每一個簇中數(shù)據(jù)樣本的均值，并以此作為各個簇新的聚類中心，計算這個數(shù)據(jù)集的誤差平方和；

5）判斷整個數(shù)據(jù)集的總誤差平方和是否保持不變或在一個很小的范圍波動，若是，則結(jié)束聚類，輸出最后聚類的結(jié)果；否則，返回步驟2）進行按步驟順序進行循環(huán)，直到滿足要求或者達到設(shè)定的迭代次數(shù)。

在實際計算中，將全事件、計量異常、電能表過載和時鐘電池欠壓以式（3）—式（6）轉(zhuǎn)化為1～100 的指標數(shù)據(jù)，實現(xiàn)數(shù)量級的統(tǒng)一，然后進行聚類算法計算。歐式距離公式為

式中：z為數(shù)據(jù)樣本；Ep為第p個聚類中心；Q為數(shù)據(jù)樣本的維度；zq、Epq為z和Ep的第q個特征。

2 案例分析

2.1 變異系數(shù)法分析

選取浙江省內(nèi)全部拆回檢定的電能表在評價時間之前6 個月的運行數(shù)據(jù)，將這些數(shù)據(jù)經(jīng)過相應(yīng)的分析處理得到智能電能表評價的指標數(shù)據(jù)，其中包括智能電能表可靠度、地區(qū)因素、全事件、計量異常、電能表過載和時鐘電池欠壓。可靠度指標和地區(qū)因素是根據(jù)全省電能表從安裝到評價時間為止的數(shù)據(jù)分析得到的，而全事件、計量異常、電能表過載和時鐘電池欠壓是以批次電能表在全年每個月的平均指標數(shù)據(jù)，采用變異系數(shù)賦值法得到每個指標的權(quán)重。并根據(jù)評價時間之前6 個月的數(shù)據(jù)得到每個智能電能表當前狀態(tài)的指標數(shù)據(jù)，最終結(jié)合可靠度指標和地區(qū)因素計算得到每個智能電能表當前的運行狀態(tài)。

規(guī)定智能電能表狀態(tài)評價表示為

式中：R為每個電能表的評價結(jié)果；w1—w4為其相應(yīng)權(quán)重，由變異系數(shù)賦值法得到。

實驗案例分析的是2018 年6 月份評價的11 個批次電能表的數(shù)據(jù)，每個批次包括拆回的電能表和仍在運行的電能表。其中指標權(quán)重以批次電能表在2017 年全年每個月的平均指標數(shù)據(jù)計算得到的，由于智能電能表拆回需要一定的流程和時間，因此以評價結(jié)果后3 個月內(nèi)拆回表的情況對評價結(jié)果的好壞進行評估。其評價結(jié)果如表1 和圖2、圖3 所示。

表1 電能表狀態(tài)評價表

表1 展示了各個分數(shù)段的具體數(shù)量情況，其中故障數(shù)量是根據(jù)當前評價時間點之后三個月內(nèi)拆回檢測得到的故障表數(shù)量，正常數(shù)量是根據(jù)當前評價時間點之后三個月內(nèi)拆回檢測得到的正常表數(shù)量。圖2 和圖3 是分別對表1 中的電能表狀態(tài)評價所有表數(shù)量分布和故障表數(shù)量分布可視化。

圖2 電能表狀態(tài)評價分布

圖3 故障電能表數(shù)量分布

從表1和圖2中可以看出，根據(jù)該方法得到的電能表評分在80 分以下的電能表數(shù)量占電能表總數(shù)量的22.09%，評分在80 分以上的電能表占電能表總數(shù)量的77.91%，同時在表1 和圖3 中可以看出評分80分以下電能表的故障數(shù)量占故障電能表總數(shù)量的73.25%。這說明該方法符合實際情況和客觀規(guī)律。

2.2 K?means算法分析

該評價方法是將智能電能表視為空間中的點，以電能表可靠度、地區(qū)因素、全事件、計量異常、電能表過載和時鐘電池欠壓這些電能表狀態(tài)評價的指標視為該空間中點的坐標，將這些點的坐標作為K?means算法的輸入，根據(jù)點和點之間的距離可以將距離相近的點聚集在同一區(qū)域，從而實現(xiàn)分類效果。

選用上文的智能電能表批次的數(shù)據(jù)進行分析，參考變異系數(shù)法的分析結(jié)果，基于數(shù)據(jù)特征和便于結(jié)果比對分析的原則下，本方法確定K?means 聚類算法的類別為5 類。聚類的結(jié)果如表2 和圖4—圖6所示。

表2 電能表狀態(tài)評價表

圖4 電能表K?means聚類結(jié)果三維分布

圖4 為K?means 聚類結(jié)果分布圖，x，y坐標是智能電能表的指標數(shù)值，z坐標是智能電能表在變異系數(shù)法得到的狀態(tài)分值，5 種顏色分別代表不同的類別。表2 展示了各個分數(shù)段的具體數(shù)量情況，其中故障數(shù)量是根據(jù)當前評價時間點之后3 個月內(nèi)拆回檢測得到的故障表數(shù)量，正常數(shù)量是根據(jù)當前評價時間點之后3 個月內(nèi)拆回檢測得到的正常表數(shù)量。圖5和圖6 是分別對表2 中的電能表狀態(tài)分類中所有表數(shù)量分布和故障表數(shù)量分布可視化。

圖5 所有電能表K?means聚類分布

圖6 故障表的K?means聚類分布

從表2 和圖5 中可以看出，Ⅰ類電表的數(shù)據(jù)最多，因為分析的是包括拆回檢定的和仍在運行的智能電能表，所以正常電能表占大多數(shù)。顯然Ⅰ類應(yīng)該是正常的電能表類別，Ⅱ～Ⅴ類作為非正常電能表類別，此處的分類是為了和變異系數(shù)法做對比。根據(jù)K?means 聚類算法得到的結(jié)果中有77.25%的電能表均被劃分為正常電能表，其中非正常電能表的類別中包含了71.81%的實際故障電能表，使用K?means算法的電能表狀態(tài)評價中分析得到的非正常電能表中包含了大部分實際故障電能表，符合實際情況，說明該評價結(jié)果具有一定的科學(xué)性和合理性。

2.3 變異系數(shù)法和K?means算法結(jié)合分析

由于單一的評價方法可能會導(dǎo)致實驗結(jié)果出現(xiàn)偏差，在分析了算法特性和電能表的特征數(shù)據(jù)分布特點后，輸出變異系數(shù)法和K?means 算法均認為正常表或故障表的結(jié)果，以此構(gòu)建新的智能電能表狀態(tài)評價方法。統(tǒng)計變異系數(shù)法和K?means 算法以及兩者結(jié)合的方法對于同一批電能表狀態(tài)評價的輸出結(jié)果，得到各個算法對于正常電能表和非正常電能表（即故障電能表）的判斷情況。

表3—表5 的數(shù)據(jù)對象為同一批電能表，即上文提及的2018 年6 月份評價的11 個批次的電能表，這其中實際正常表和故障表數(shù)量是以評價時間點后3個月實際拆回檢定結(jié)果統(tǒng)計的。以表3 為例，分別從正常表和故障表的角度分析變異系數(shù)法預(yù)測結(jié)果，該批電能表的故障表數(shù)量為2 026個，檢驗得到其中實際故障表2 026 個，實際正常表0 個；預(yù)測正常表數(shù)量為977 個，檢驗得到其中實際故障表759 個，實際正常表218個。

表3 變異系數(shù)法判斷統(tǒng)計表

表4 K?means算法判斷統(tǒng)計表

表5 K?means+變異系數(shù)法判斷統(tǒng)計表

從上文看出，K?means和變異系數(shù)法結(jié)合的方法在保持故障表預(yù)測高準確率的情況下，大幅提高了對于正常表的預(yù)測準確率，說明這兩種算法結(jié)合后的性能比其中任意一種單一的算法的性能更好。

對智能電能表進行狀態(tài)評價的最終目的是通過對智能電能表的運行狀態(tài)進行評估，找出具有較高可能性發(fā)生故障的電能表進行排查，通過對于變異系數(shù)評價方法、基于K?means 算法的評價方法以及變異系數(shù)法和K?means 結(jié)合的電能表評價方法的結(jié)果進行定量比較，進一步說明綜合電能表評價方法的合理性。

采用二分類模型中的精確率、召回率和F1?score衡量評價模型的好壞。1）精確率。

式中：Pr為精確率；TP為預(yù)測為故障的電能表中實際為故障的電能表數(shù)量；FN為預(yù)測為故障的電能表實際為正常電能表的數(shù)量。

2）召回率。

式中：Re為召回率；FE為預(yù)測為正常電能表實際為故障電能表的數(shù)量。

（3）F1?score。

式中：Fs為F1?score，為精確率和召回率的調(diào)和平均數(shù)，精確率反映了模型對負樣本的區(qū)分程度，召回率反映了模型對正樣本的識別力，F(xiàn)1?score是兩者的綜合，F(xiàn)1?score值越大，說明分類模型越穩(wěn)健。對于變異系數(shù)評價法、基于K?means 算法的評價方法以及變異系數(shù)法和K?means 方法結(jié)合的評價方法的二分類模型參數(shù)進行比較，如表6所示。

表6 各算法的結(jié)果比較

由表6 可以看出，K?means 和變異系數(shù)法結(jié)合的評價方法在召回率和F1?score上都有所提高，這是由于這兩種方法結(jié)合后縮小了評價的范圍，但同時進一步提高了評價的準確率和召回率，由于供電公司比較看重評價為故障表的電能表實際是否為故障表的比例，因此該評價方法的結(jié)果更符合供電公司要求。

3 結(jié)語

介紹變異系數(shù)算法和K?means 算法的原理，并基于這兩種算法對電能表的狀態(tài)進行評價，在比較變異系數(shù)法和K?means 算法的評價結(jié)果后，結(jié)合電能表評價參數(shù)的特性，融合了變異系數(shù)法和K?means算法構(gòu)建了新的狀態(tài)評價方法，通過數(shù)據(jù)分析證明其科學(xué)行和可行性，為智能電表的狀態(tài)評價提供新的思路。該方法在實際實驗中得到了供電公司的認可。