CN104569576B - 智能电网非法用电行为分布式检测系统 - Google Patents

Abstract

本发明公开一种智能电网非法用电行为分布式检测系统，包括智能电表数据采集模块、智能电表数据传输模块、智能电表数据存储模块、智能电表数据安全验证模块和智能电表核心计算模块；以及检测方法。本发明实现了对智能电网的全局覆盖，依托于智能电表完成非法用电检测，在系统部署完成后，相较原有集中式检测系统，大大扩展了系统的覆盖度,也更能应对突发情况，计算和检测的精度更高。

Description

智能电网非法用电行为分布式检测系统

技术领域

本发明所涉及的是智能电网数据分析领域中的非法用电行为检测系统，具体是一种高效率快速响应智能电网非法用电行为分布式检测系统及其检测方法。

背景技术

智能电网是针对下一代电网改造提出的规划设想，其核心意义在于引入智能设备，更新原有机械装置，使得电网全局信息交互更为方便直接，调控配置也智能化。直观表现有，引入智能电表，实现电力企业和电网用户间的双向信息交互。

非法用电行为是指用户不按照合同使用电能，一般包含电力欺诈，偷电漏电等，电力行业中，因这种行为产生的耗费被称作非技术性损耗。这一行为对于电网产生了巨大危害。据统计，全世界各国每年生产的电能中，高达20％的份额为非技术性损耗。非法用电行为主要有两种攻击手段：一种为物理攻击，即通过对电表、变压器、传输线路等进行物理性篡改(短接电表电路、破坏变压器内加压煤油等)，使得自身电能消费不被记录；另一种为网络攻击，即针对电网内的智能设备(智能电表、智能变压器等)进行网络式入侵，篡改内部程序，更改计数器设置，同样来逃避电能消耗的计费。

传统电网中，非法用电行为的检测通常会有如下几个阶段：1.一线电力巡查人员或者终端用户报备可疑用户、可疑片区及可疑行为；2.电力企业派出技术人员采集可疑用户或片区的电力数据，为了提高检测的准确度，通常会采样2-3个月的数据；3.数据分析人员针对采样数据进行周密分析，以此判断是否非法用电。该方案对于人力资源耗费巨大，且时间周期长，精度不高，对于电网的覆盖也十分有限。

智能电网时代，人们借助智能设备，使得电力数据的采集更为快捷，但相比传统电网，智能设备的引入，使得非法用电行为的攻击手段更为多样，检测也更为困难。非法用电行为检测作为领域内热点问题，其解决方法实施存在诸多难点，主要表现为：

1)全网覆盖：电网作为关乎国计民生的骨干网络，节点众多，网络复杂，全局而言，异常行为的产生几乎是每时每刻，对于全网的监控是一大难点。

2)快速响应：智能电网产生的数据量庞大，非法用电行为的分布面又十分广泛，如何开发算法，使得检测机制能够快速响应是又一难点。

3)应对复杂多变的攻击行为：电网，尤其是智能电网引入的智能设备，使得针对电网的非法攻击行为在原有物理方式之外，增加了网络攻击的可能，网络攻击呈现复杂多变的趋势，非法用电者呈现的行为模式也难以界定。如何应对不断“创新”的非法用电行为是另一难点。

通过对现有技术的检索发现，现有的检测方法大多为集中式的机器学习分类方法。如美国Toledo大学电子工程及计算机科学系的Soma Shekara Sreenadh ReddyDepuru、Lingfeng Wang等人提出的基于支持向量机的非法用电行为检测系统，首先对现有电网内数据进行采样，获取非法用电用户和诚信用户的电力消费数据，对数据样本进行标定并分离特征值后，输入支持向量机，对数据完成建模后，生成数据分类处理程序，之后才将系统投入运行，将电网内各用户一段时间内的电力消费数据交由支持向量机生成的分类处理程序检测判断，如果程序将其分类为非法用电行为，则标定出该用户，否则则略过该数据。此方法有相当的局限性：1.对于稠密网络，机器学习检测方法运行时间长，效率低下；2.检测精度并不高，按现有应用实例来看，准确度在60％左右，最后仍然需要数据分析人员来重新验证。3.训练样本依赖度高，且针对其特征值不能及时更新，无法做到“实时”的行为检测。其他关于非法用电行为检测的尝试，采用的方法有统计学方法、决策树、神经网络算法及遗传算法等，这些方法都或多或少具有上述局限性。

发明内容

本发明针对现有系统存在的上述不足，提出一种基于线性方程简化求解的分布式高效率快速反应的智能电网非法用电行为检测系统，借助智能电网中的智能设备的可计算能力，分布式地就地检测用户的非法用电行为，避免了数据集中式处理带来的延迟耗费，并能稳定地运行处理稠密电网拓扑，同时天然地覆盖全局电网，“实时”地分析检测智能电网中的非法用电行为。

本发明是通过以下技术方案实现的：

一种智能电网非法用电行为分布式检测系统，其特点在于，包括智能电表数据采集模块、智能电表数据传输模块、智能电表数据存储模块、智能电表数据安全验证模块和智能电表核心计算模块；

所述的智能电表数据采集模块负责收集电表附属用户的电能消费信息；所述的智能电表数据传输模块通过智能电表间的无线传输连接，构建基础的数据接收和发送单元，用以传递或者接收相关的计算结果，并借助数据安全验证模块检测数据来源及数据本身的可靠性；所述的智能电表数据存储模块负责将数据按指定形式、结构存储在电表内存中，并向外提供数据存储调用接口；所述的智能电表数据安全验证模块，用于对数据完整性进行检测以及访问权限进行控制，维护数据传递过程中的一致性及数据传输参与双方的身份可靠性；所述的智能电表核心计算模块，计算本地附属用户的非法用电检测参数及可推导的其他用户的智能电表数据，并调用数据传输模块将计算结果传输出去。

所述的智能电表数据采集模块包括：数据采样模块和数据存储调用模块；所述的数据采样模块连接智能电表数据接口，按照设定的采样周期，获取周期内的智能电表数据和；所述的数据存储调用模块将采样数据以键值对数组的形式组织，将之存储在智能电表内存中，键值对是指时间周期与相应电力数据的一一对应。

所述的智能电表数据传输模块包括：数据接收模块和数据发送模块，所述的数据接收模块负责和数据发送方建立数据传输通路，接收其发送的数据，并调用智能电表数据存储模块，以对象式文件的形式完成数据存储；所述的数据发送模块负责和数据接收方建立数据传输通路，调用智能电表数据存储模块获取指定数据内容，向其发送数据。

所述的智能电表数据存储模块包括：数据读取模块和数据写入模块，所述的数据读取模块负责接收数据文件名，查阅文件名和数据内容存储地址对照表，从智能电表内存中的实际物理地址处，读取相对应的数据内容；所述的数据写入模块负责接收数据文件名、数据类型、数据组织形式等参数信息和数据内容本身，将实际数据内容按照指定的数据类型和组织形式写入智能电表内存。

所述的数据安全验证模块包括：数据完整性计算模块和数据源身份验证模块，所述的数据完整性计算模块获取数据中的校验和数据信息，计算当前数据内容对应的校验和值，与数据自带的校验和值进行比对，如果相同，则说明数据的完整性并未受到破坏，如果不同，则丢弃该数据；所述的数据源身份验证模块获取数据中发送方信息，经由接收方私钥解密后，鉴定数据源是否可靠，如果解密后获取的数字序列在接收方验证库可以查到，则数据源可靠，否则丢弃数据。

所述的智能电表核心计算模块包括：LU矩阵分解模块和回溯数值替换模块，所述的LU矩阵分解模块通过智能电表数据传输模块获取本地计算所需的基础数据后，调用智能电表数据安全验证模块验证安全后，计算矩阵的相关列数据，完成计算后，将计算结果通过智能电表数据传输模块发送出去；所述的回溯数值替换模块通过智能电表数据传输模块获取本地计算所需的基础数据后，调用智能电表数据安全验证模块验证安全后，计算矩阵的相关行数据及本身参数，完成计算后，将计算结果通过智能电表数据传输模块发送出去。

一种智能电网非法用电行为分布式检测方法，其特点在于，该方法包括如下步骤：

①部署过程：在指定电力小区的各住户处及变电站处分别安装具备上述的智能电网非法用电行为分布式检测系统的智能电表，并对各智能电表进行编号SM₀至SM_m，其中，SM₀为变电站处智能电表，SM₁至SM_m为第一住户处智能电表至第m住户处智能电表；

②数据采样：从t₁至t_n时刻，各住户处智能电表由各自的智能电表数据采样模块分别采样对应住户的电能消耗值，变电站处智能电表号SM₀由智能电表数据采样模块采样小区总电能消耗值，变电站处智能电表号SM₀将t_j时刻的小区总电能消耗值经智能电表数据传输模块发送给对应智能电表SM_j，j为1到n之间的任意正整数；

③矩阵分解过程：变电站处智能电表SM₀的智能电表核心计算模块计算下三角矩阵L协作计算向量y的第一个元素y1，并通过智能电表数据传输模块传输给第一住户处智能电表SM₁，第一住户处智能电表SM₁经由数据传输模块接收下三角矩阵L协作计算向量y的第一个元素y1，并经智能电表数据安全验证模块验证数据完整性及其来源后，计算下三角矩阵L与上三角矩阵U的第1列，并将该计算结果存储在第一住户处智能电表SM₁的内存中，同时传输给第二住户处智能电表SM₂；依次进行，第j个住户处的智能电表SM_j计算下三角矩阵L与上三角矩阵U的第j列，j为1到m之间的任意正整数，直至第m住户处智能电表SM_m，完成下三角矩阵L与上三角矩阵U第1至m列的计算，计算结果分别存储在各住户处智能电表的本地内存中；

计算下三角矩阵L与上三角矩阵U的第j列的具体步骤如下：

1)下三角矩阵L第j列计算：

对于下三角矩阵L第j列第i个数据l_i,j，i为1到m间的任意正整数，

当i<j,下三角矩阵L第j列第i个数据l_i,j等于0；

当j＝1,下三角矩阵L第j列第i个数据l_i,j等于当前智能电表在t_i时刻采样存储的电能消耗值p_ti,j除以当前智能电表在t₁时刻采样存储的电能消耗值p_t1,j；

当2<＝j<＝m且i>j，首先计算当前智能电表在t_i时刻采样存储的电能消耗值p_ti,j与下三角矩阵L的前j列子矩阵和上三角矩阵U的前j行子矩阵的乘积和间的差值diff，接着计算该差值diff除以上三角矩阵U的第j列第j个数据值所得的商，即为下三角矩阵L第j列第i个数据l_i,j的计算结果。

参照上述过程，逐个计算下三角矩阵L第j列的各个数据值；

2)上三角矩阵U第j列计算：

对于上三角矩阵U第j列第i个数据u_i,j，i为1到m间的任意正整数，

当i>j,计算上三角矩阵U第j列第i个数据u_i,j等于0；

当i＝1,计算上三角矩阵U第j列第i个数据u_i,j等于当前智能电表在t_i时刻采样存储的电能消耗值p_ti,j；

当2<＝i<＝m且i<j，计算上三角矩阵U第j列第i个数据u_i,j等于当前智能电表在t_i时刻采样存储的电能消耗值p_ti,j减去下三角矩阵L的前i列子矩阵和上三角矩阵U的前i行子矩阵的乘积和间所得的差值；

参照上述过程，逐个计算上三角矩阵U第j列的各个数据值。

④回溯数值替换过程：第m住户处智能电表SM_m的智能电表核心计算模块从数据存储模块中获取存储于智能电表内存的矩阵分解过程计算结果L和U的1至m列数据信息，计算上三角矩阵U第m列第m-1个元素值u_m-1,m以及对应住户的诚信参数k_m，并将这两个值经由数据传输模块传递给变电站处智能电表SM₀和下一个智能电表，即第m-1住户处智能电表SM_m-1，第m-1住户处智能电表SM_m-1通过安全验证模块验证数据完整性及其来源，并从数据存储模块中获取存储于智能电表内存的矩阵分解过程计算结果L和U的1至m-1列数据信息，开始计算上三角矩阵U第m-1列第m-2个元素值u_m-2,m-1以及对应住户的诚信参数k_m-1，其将该计算结果经由传输模块传输给变电站处智能电表SM₀和下一个智能电表，即第m-2住户处智能电表SM_m-2；依次进行至第一住户处智能电表SM₁，完成全部住户的诚信参数的计算；

计算上三角矩阵U第m列第m-1个元素值u_m-1,m以及对应住户的诚信参数k_m的具体步骤如下：

1)上三角矩阵U第m列第m-1个元素值u_m-1,m等于当前智能电表在t_m-1时刻采样存储的电能消耗值p_tm-1,m减去下三角矩阵L的前m-1列子矩阵和上三角矩阵U的前m-1行子矩阵的乘积和所得的差值；

2)对应住户的诚信参数k_m的计算：

首先计算t_m时刻变电站处智能电表采样存储的小区电能耗费值P_tm减去下三角矩阵L的前m列子矩阵和上三角矩阵U的前m行子矩阵的乘积和所得的差值D，接着计算差值D减去上三角矩阵U第m+1列到第n列的子矩阵与第m+1个住户到第n个住户的诚信参数向量间的乘积和所得的差值E，最后计算差值E除以上三角矩阵U的第m列第m个数据值所得的商，该商值即为对应住户的诚信参数k_m。

⑤统计计算过程：变电站处智能电表SM₀获取全部住户的诚信参数后，依次检查其值是否大于1：

当k_j＝1，则用户j为诚信用户，没有窃电；

当k_j>1，则用户j实际消费量大于记录量，即存在偷电行为；

当0<k_j<1，则用户j实际消费量小于记录量，即电表故障。

与现有技术相比，本发明的有益效果是：1)本系统实现了对智能电网的全局覆盖，依托于智能电表完成非法用电检测，在系统部署完成后，即自然地实现了电网的全局覆盖，相较原有集中式检测系统，大大扩展了系统的覆盖度；2)本系统可进行实时监控，对电网用户电力数据进行实时采样监控，相较原有系统的漫长采样周期而言，效率更高，也更能应对突发情况；3)本系统的检测精度更高，线性方程组的求解过程中，将所覆盖电网内的全部用户各时间段的数据情况都计算在内，相较原有系统只计算单一用户的若干参数，计算和检测的精度更高。

附图说明

图1系统架构

图2小区电力规划图

图3LU矩阵分解过程

图4回溯替代过程

图5系统运行流程图

具体实施方式

下面对本发明的实施过程作举例详细说明，本实施例是在附图2所示小区中运行并以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程。

如图1所示，本实施例包括：(1)智能电表数据采集模块、(2)智能电表数据传输模块、(3)智能电表数据存储模块(4)智能电表数据安全验证模块、(5)智能电表核心计算模块。其中：智能电表数据采集模块负责收集电表附属用户的电能消费信息，其分布式地部署于各智能电表中，收集记录用户一定时间内的电能消费和；智能电表数据传输模块通过智能电表间的无线传输连接，构建基础的数据接收和发送单元，用以传递或者接收相关的计算结果，并借助数据安全验证模块检测数据来源及数据本身的可靠性；智能电表数据存储模块负责将数据按指定形式、结构存储在电表内存中，其向外提供数据存储调用接口；智能电表数据安全验证模块集成数据完整性检测及访问权限控制功能，维护数据传递过程中的一致性及数据传输参与双方的身份可靠性；智能电表核心计算模块，依照核心算法，计算本地附属用户的非法用电检测参数及可推导的其他用户的智能电表数据，并调用数据传输模块将计算结果传输出去。

所述的智能电表数据采集模块包括：数据采样模块、数据存储调用模块。其中：数据采样模块连接智能电表数据接口，按照设定的采样周期，获取周期内的智能电表数据和；数据存储调用模块将采样数据以键值对数组的形式组织，调用智能电表数据存储模块，将之存储在智能电表内存中，键值对是指时间周期与相应电力数据的一一对应。

所述的智能电表数据传输模块包括：数据接收模块、数据发送模块。数据传输通路均基于无线通信信道建立，其中：数据接收模块负责和数据发送方(任意其他智能电表)建立数据传输通路，接收其发送的数据，并调用智能电表数据存储模块，以对象式文件的形式完成数据存储；数据发送模块负责和数据接收方(本智能电表指定的传输目的智能电表)建立数据传输通路，调用智能电表数据存储模块获取指定数据内容，向其发送数据。

所述的智能电表数据存储模块包括：数据读取模块、数据写入模块。其中：数据读取模块负责接收数据文件名，查阅文件名和数据内容存储地址对照表，从智能电表内存中的实际物理地址处，读取相对应的数据内容；数据写入模块负责接收数据文件名、数据类型、数据组织形式等参数信息和数据内容本身，将实际数据内容按照指定的数据类型和组织形式写入智能电表内存

所述的数据安全验证模块包括：数据完整性计算模块、数据源身份验证模块。其中：数据完整性计算模块获取数据中的校验和数据信息，根据消息摘要算法第五版(MD5)计算当前数据内容对应的校验和值，与数据自带的校验和值进行比对，如果相同，则说明数据的完整性并未受到破坏，如果不同，则丢弃该数据；数据源身份验证模块获取数据中发送方信息，主要为一串经由当前数据接收方公钥加密的数字序列，经由接收方私钥解密后，鉴定数据源是否可靠，如果解密后获取的数字序列在接收方验证库可以查到，则数据源可靠，否则丢弃数据。

所述的智能电表核心计算模块包括：LU矩阵分解模块、回溯数值替换模块。二者运行基础为本系统的核心算法，稍后会进行介绍。其中：LU矩阵分解模块通过智能电表数据传输模块获取本地计算所需的基础数据后，计算矩阵的相关列数据，完成计算后，将计算结果通过智能电表数据传输模块发送出去；回溯数值替换模块通过智能电表数据传输模块获取本地计算所需的基础数据后，计算矩阵的相关行数据及本身参数，完成计算后，将计算结果通过智能电表数据传输模块发送出去。

本发明的核心算法是基于如下的数学模型：

如图2所示是某个小区的电力规划图。其中，小区变电站、每家每户都安装了智能电表。这些智能电表通过无线网络或有线网络传递实时的用电信息。假设，小区包含N个用户，我们定义信息采集周期为SP，即每个采样周期内，N+1(包含总电表)个智能电表都会统计各自在前一周期内的用电量，分别标记为P_t,j(用户j在时刻t的电力消费值)和TP_t(变电站电表在时刻t记录的用电量)。K_j是每个用户诚信度，则K_j*P_t,j代表用户j从(t-SP)～t时刻期间的实际耗电量。(下面的描述将以Smart Meter的缩写SM代称智能电表)

所有用户上报的用电量与变电站电表的计量值之间，存在以下等式，

其中，我们的目标是获知每个用户的K值(诚信参数)。

定义，

1)如果K_j＝1，则用户j为诚信用户，没有窃电；

2)如果K_j>1，则用户j实际消费量大于记录量，即存在偷电行为；

3)如果0<K_j<1，则用户j实际消费量小于记录量，即电表故障。

求解n个变量(即k_1～n)，需要n个线性方程组。即假设每户的用电量(即p_t,i)是随时间变化的，则通过不同时段多次采集，即可获得一系列的方程组，如下所示：

从而，非法用电的行为检测变成了线性方程组的求解过程，其可简化为矩阵相乘：

Pk＝∑P′ (2)

针对上述的线性方程组，分布式算法如下：

将上述方程组求解简化为矩阵分解，将上述等式归约为矩阵P，将该矩阵P分解为下三角矩阵L与上三角矩阵U的乘积，意即：

∑P′＝LU (3)

对于上三角矩阵U，其中元素可以按照如下过程计算：

u_i，j＝0， i＞j

其中，为矩阵P中第j列中的第i个元素，同理，下三角矩阵L中的元素可以按照如下过程计算：

l_i，j＝0， i＜j

L和U能够互相协作式地计算求解，在其求解完成后，有：

Ly＝∑P′ (6)

Uk＝y (7)

对于第j-1个智能电表SM_j-1来说，其计算y_j由以下过程推出：

最终，每个SM_j计算其对应诚信参数K_j,使用回溯替代法：k_n＝y_n/u_n,n,以及

对于j<n的情况。

故而算法整体需要两个过程：分布式LU矩阵分解和回溯替代过程。伪代码实现如附图3和4所示。

SM₀计算y₁＝P_t1，之后传递该值给SM₁，SM₀不需要计算L U中的任意元素，SM1接收y₁后，其计算U中的第1列和L中的第1列以及y₂，之后传递y₁，y₂，L中的第1列给SM₂.对于SM_j，它获取y₁到y_j的所有y值，和L中的第1列至第j-1列，一次计算U中的第j列，L中的第j列，和y_j+1,之后SM_j将y₁到y_j+1和L的1至j列的值传递给SM_j+1。最终SM_n计算完毕L和U中的所有元素值，以及所有的y值，进而发通知给总智能电表，让其开启回溯替换过程。对于回溯替换过程，从SM_n开始，参照公式8，开始计算，u_n-1,n以及k_n，并将这两个值传递给SM_n-1，SM_n-1依照上数值计算k_n-1以及u_n-2,n-1，依次类推，直到所有的参数值都被计算出来，每个电表计算出各自参数后，将值传递给总电表计算，总电表可根据诚信参数值来鉴别所有的非法用电用户。

系统运行流程

如附图5所示，本发明分布式非法用电检测系统运行包含5个步骤，分别为：系统部署过程、数据采样过程、矩阵分解过程、回溯数值替换过程、统计计算过程。下面就这5个步骤的详细情况做说明。

所述的系统部署过程，包含软件模块安装过程、软件运行初始化过程。其中：软件模块安装过程是将上述采样、传输、存储、安全、计算等智能电表处理模块安装在住户智能电表和变电站处智能电表(此处称为电力数据收集器)中；软件运行初始化过程是将系统各模块运行的参数初始化，具体为电力数据收集器中计算模块初始化自身序号为SM₀，并根据住户总数n来标定其他各智能电表序号为SM₁至SM_n，并将该序号传递给各对应智能电表。完成系统部署过程后，系统开始运行数据采样过程。

所述的数据采样过程，是系统开始运行后的第一项任务，根据初始化时设定的住户总数n，各住户智能电表记录t₀至t_n n个时刻的住户30分钟内电能消耗值(单位千瓦时KWH)，电力数据收集器同样记录t₀至t_n n个时刻的小区电力总消费值。当完成所有n个时刻的数据采样后，电力数据收集器通过数据传输模块将t_j时刻的小区电力总消费值P_tj发送给智能电表SM_j。传输完毕后，系统开始运行矩阵分解过程。

所述的矩阵分解过程，是系统核心计算的第一部分，其从电力数据收集器SM₀开始，其智能电表核心计算模块根据上述核心算法计算对应参数y₁为t₁时刻的电力总消费值，将该数据通过智能电表数据传输模块传输给智能电表SM₁；从此处开始，核心计算模块运行稍有不同，对于智能电表SM_j，其核心计算模块从上一个智能电表SM_j-1处获取计算相关数据，并经由数据安全验证模块验证数据源为SM_j-1且数据本身是完整的，根据核心算法，其开始计算矩阵L的第j列数据和矩阵U的第j列数据，完成计算后，其将计算结果经由数据存储模块存储在智能电表本地内存中，并将该计算结果L的第1至j列数据和U的第1至j列数据经由数据传输模块发送给下一个智能电表SM_j+1；如此分布式运行，直到最后一个智能电表SM_n完成其计算。之后，系统开始运行回溯数值替代过程。

所述的回溯数值替代过程，是系统核心计算的第二部分，其从最后一个智能电表SM_n开始，其核心计算模块参照核心算法中公式8和9，从数据存储模块中获取存储于智能电表内存的矩阵分解过程计算结果，开始计算u_n-1,n以及对应住户的诚信参数k_n，并将这两个值经由数据传输模块传递给下一个智能电表SM_n-1和电力数据收集器SM₀；依次类推，对于对于智能电表SM_j，其核心计算模块从上一个智能电表SM_j+1处获取计算相关数据，并经由数据安全验证模块验证数据源为SM_j+1且数据本身是完整的，根据核心算法中公式8和9，从数据存储模块中获取存储于智能电表内存的矩阵分解过程计算结果，计算u_j-1,j以及对应住户的诚信参数k_j，将这两个值经由数据传输模块传递给下一个智能电表SM_j-1和电力数据收集器SM₀；如此分布式运行，直到第一个智能电表SM₁完成其计算。之后，系统开始运行统计计算过程。

所述的统计计算过程，是系统运行的最后一个过程，其运行在电力数据收集器SM₀中，其将所有的大于1的诚信参数值和对应智能电表的序列号输出，至此完成一轮非法用电行为检测，之后电力数据收集器向其他住户的智能电表发出重置信号，系统重新进入数据采样过程。

本实施例的具体操作过程包含系统部署过程、数据采样过程、矩阵分解过程、回溯数值替换过程、统计计算过程，如下：

1)系统部署过程：在指定电力小区(图2所示)中部署本检测系统，各户智能电表中安装检测系统实例，并按照拓扑，指定各实例的序列号，住户A-D的智能电表序列号为SM1至SM4，SM0标注变电站处的电力数据采集器。系统开始运行采样过程。

2)数据采样过程：从t1时刻(此处为上午8点)开始，SM1至SM4经由智能电表数据采样模块采样对应住户的30分钟电能消耗值(单位千瓦时)，SM0经由智能电表数据采样模块采样小区总电能消耗值，采样至t4时刻(此处为上午10点)结束，SM0将tj时刻的小区总电能消耗值经由智能电表数据传输模块发送给对应智能电表SMj。系统开始运行矩阵分解过程。

3)矩阵分解过程：自SM₀开始，运行智能电表核心计算模块，参照核心算法计算y1，完成计算后，其将y1经由智能电表数据传输模块传输给SM1，SM1经由数据传输模块接收y1，经由安全验证模块验证数据完整性和来源为SM0后，根据公式(4)计算矩阵L和U的第1列，其将该计算结果经由传输模块传输给SM2，并经由存储模块将结果存储在SM1内存中；依次进行，SM1至SM2，SM2至SM3，SM3至SM4，完成L和U全部列数据的计算。系统开始运行回溯数值替代过程。

SM₃智能电表计算完毕自身诚信参数k₀-k₃后，将其发送给变压站处的智能电表SM₄，SM₄判断各参数大小，如果其大于1，则标定该用户为非法用电用户，如果其小于1，则标定该用户为受害用户。SM₄将上述标定信息经由警报界面信息生成模块转化后提供给一线运维人员查看。

4)回溯数值替换过程：自SM4开始，运行智能电表核心计算模块，参照核心算法中公式8和9，从数据存储模块中获取存储于智能电表内存的矩阵分解过程计算结果L和U的1至4列数据信息，开始计算u_3,4以及对应住户的诚信参数k₄，完成计算后，其将u_3,4经由智能电表数据传输模块传输给SM3和SM0，SM3经由数据传输模块接收u_3,4，经由安全验证模块验证数据完整性和来源为SM4后，参照核心算法中公式8和9，从数据存储模块中获取存储于智能电表内存的矩阵分解过程计算结果L和U的1至3列数据信息，开始计算u_2,3以及对应住户的诚信参数k₃，其将该计算结果经由传输模块传输给SM2和SM0，；依次进行，SM2至SM1，完成k₁至k₄全部诚信参数的计算。系统开始运行统计计算过程。

5)统计计算过程：SM0获取k₁至k₄全部诚信参数后，检查其值是否大于1，将k值大于1的智能电表用户标定为非法用电用户，至此完成本轮非法用电检测。系统持续运行，进行数据采样过程，并不停歇。

通过该实施例，我们提出的一种高效率快速响应的分布式非法用电行为检测系统，以其覆盖全局电网的特性，快速检测智能电网中非法用电行为，分布式的形式使其计算效率大大提高，避免了集中式机器学习方法处理的若干缺陷。

最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (8)

1.一种智能电网非法用电行为分布式检测系统，包括智能电表数据采集模块(1)、智能电表数据传输模块(2)、智能电表数据存储模块(3)、智能电表数据安全验证模块(4)和智能电表核心计算模块(5)；所述的智能电表数据采集模块负责收集电表附属用户的电能消费信息；所述的智能电表数据传输模块通过智能电表间的无线传输连接，构建基础的数据接收和发送单元，用以传递或者接收相关的计算结果，并借助数据安全验证模块检测数据来源及数据本身的可靠性；所述的智能电表数据存储模块负责将数据按指定形式、结构存储在电表内存中，并向外提供数据存储调用接口；所述的智能电表数据安全验证模块，用于对数据完整性进行检测以及访问权限进行控制，维护数据传递过程中的一致性及数据传输参与双方的身份可靠性；所述的智能电表核心计算模块，计算本地附属用户的非法用电检测参数及可推导的其他用户的智能电表数据，并调用数据传输模块将计算结果传输出去；其特征在于，所述的智能电表核心计算模块(5)包括：LU矩阵分解模块和回溯数值替换模块，所述的LU矩阵分解模块通过智能电表数据传输模块获取本地计算所需的基础数据后，调用智能电表数据安全验证模块验证安全后，计算矩阵的相关列数据，完成计算后，将计算结果通过智能电表数据传输模块发送出去；所述的回溯数值替换模块通过智能电表数据传输模块获取本地计算所需的基础数据后，调用智能电表数据安全验证模块验证安全后，计算矩阵的相关行数据及本身参数，完成计算后，将计算结果通过智能电表数据传输模块发送出去。

2.根据权利要求1所述的智能电网非法用电行为分布式检测系统，其特征在于，所述的智能电表数据采集模块(1)包括：数据采样模块和数据存储调用模块；所述的数据采样模块连接智能电表数据接口，按照设定的采样周期，获取周期内的智能电表数据和；所述的数据存储调用模块将采样数据以键值对数组的形式组织，将之存储在智能电表内存中，键值对是指时间周期与相应电力数据的一一对应。

3.根据权利要求1所述的智能电网非法用电行为分布式检测系统，其特征在于，所述的智能电表数据传输模块(2)包括：数据接收模块和数据发送模块，所述的数据接收模块负责和数据发送方建立数据传输通路，接收其发送的数据，并调用智能电表数据存储模块，以对象式文件的形式完成数据存储；所述的数据发送模块负责和数据接收方建立数据传输通路，调用智能电表数据存储模块获取指定数据内容，向其发送数据。

4.根据权利要求1所述的智能电网非法用电行为分布式检测系统，其特征在于，所述的智能电表数据存储模块(3)包括：数据读取模块和数据写入模块，所述的数据读取模块负责接收数据文件名，查阅文件名和数据内容存储地址对照表，从智能电表内存中的实际物理地址处，读取相对应的数据内容；所述的数据写入模块负责接收数据文件名、数据类型、数据组织形式参数信息和数据内容本身，将实际数据内容按照指定的数据类型和组织形式写入智能电表内存。

5.根据权利要求1所述的智能电网非法用电行为分布式检测系统，其特征在于，所述的数据安全验证模块(4)包括：数据完整性计算模块和数据源身份验证模块，所述的数据完整性计算模块获取数据中的校验和数据信息，计算当前数据内容对应的校验和值，与数据自带的校验和值进行比对，如果相同，则说明数据的完整性并未受到破坏，如果不同，则丢弃该数据；所述的数据源身份验证模块获取数据中发送方信息，经由接收方私钥解密后，鉴定数据源是否可靠，如果解密后获取的数字序列在接收方验证库可以查到，则数据源可靠，否则丢弃数据。

6.一种智能电网非法用电行为分布式检测方法，其特征在于，该方法包括步骤：

①部署过程：在指定电力小区的各住户处及变电站处分别安装具备权利要求1所述的智能电网非法用电行为分布式检测系统的智能电表，并对各智能电表进行编号SM₀至SM_m，其中，SM₀为变电站处智能电表，SM₁至SM_m为第一住户处智能电表至第m住户处智能电表；

②数据采样：从t₁至t_n时刻，各住户处智能电表由各自的智能电表数据采样模块分别采样对应住户的电能消耗值，变电站处智能电表号SM₀由智能电表数据采样模块采样小区总电能消耗值，变电站处智能电表号SM₀将t_j时刻的小区总电能消耗值经智能电表数据传输模块发送给对应智能电表SM_j，j为1到n之间的任意正整数；

③矩阵分解过程：变电站处智能电表SM₀的智能电表核心计算模块计算下三角矩阵L协作计算向量y的第一个元素y1，并通过智能电表数据传输模块传输给第一住户处智能电表SM₁，第一住户处智能电表SM₁经由数据传输模块接收下三角矩阵L协作计算向量y的第一个元素y1，并经智能电表数据安全验证模块验证数据完整性及其来源后，计算下三角矩阵L与上三角矩阵U的第1列，并将该计算结果存储在第一住户处智能电表SM₁的内存中，同时传输给第二住户处智能电表SM₂；依次进行，第j个住户处的智能电表SM_j计算下三角矩阵L与上三角矩阵U的第j列，j为1到m之间的任意正整数，直至第m住户处智能电表SM_m，完成下三角矩阵L与上三角矩阵U第1至m列的计算，计算结果分别存储在各住户处智能电表的本地内存中；

④回溯数值替换过程：第m住户处智能电表SM_m的智能电表核心计算模块从数据存储模块中获取存储于智能电表内存的矩阵分解过程计算结果L和U的1至m列数据信息，计算上三角矩阵U第m列第m-1个元素值u_m-1,m以及对应住户的诚信参数k_m，并将这两个值经由数据传输模块传递给变电站处智能电表SM₀和下一个智能电表，即第m-1住户处智能电表SM_m-1，第m-1住户处智能电表SM_m-1通过安全验证模块验证数据完整性及其来源，并从数据存储模块中获取存储于智能电表内存的矩阵分解过程计算结果L和U的1至m-1列数据信息，开始计算上三角矩阵U第m-1列第m-2个元素值u_m-2,m-1以及对应住户的诚信参数k_m-1，其将该计算结果经由传输模块传输给变电站处智能电表SM₀和下一个智能电表，即第m-2住户处智能电表SM_m-2；依次进行至第一住户处智能电表SM₁，完成全部住户的诚信参数的计算；

⑤统计计算过程：变电站处智能电表SM₀获取全部住户的诚信参数后，依次检查其值是否大于1：

当k_j＝1，则用户j为诚信用户，没有窃电；

当k_j>1，则用户j实际消费量大于记录量，即存在偷电行为；

当0<k_j<1，则用户j实际消费量小于记录量，即电表故障。

7.根据权利要求6所述的智能电网非法用电行为分布式检测方法，其特征在于，所述的步骤③中计算下三角矩阵L与上三角矩阵U的第j列的具体步骤如下：

1)下三角矩阵L第j列计算：

对于下三角矩阵L第j列第i个数据l_i,j，i为1到m间的任意正整数，

当i<j,下三角矩阵L第j列第i个数据l_i,j等于0；

当j＝1,下三角矩阵L第j列第i个数据l_i,j等于当前智能电表在t_i时刻采样存储的电能消耗值p_ti,j除以当前智能电表在t₁时刻采样存储的电能消耗值p_t1,j；

当2<＝j<＝m且i>j，首先计算当前智能电表在t_i时刻采样存储的电能消耗值p_ti,j与下三角矩阵L的前j列子矩阵和上三角矩阵U的前j行子矩阵的乘积和间的差值diff，接着计算该差值diff除以上三角矩阵U的第j列第j个数据值所得的商，即为下三角矩阵L第j列第i个数据l_i,j的计算结果；

参照上述过程，逐个计算下三角矩阵L第j列的各个数据值；

2)上三角矩阵U第j列计算：

对于上三角矩阵U第j列第i个数据u_i,j，i为1到m间的任意正整数，

当i>j,计算上三角矩阵U第j列第i个数据u_i,j等于0；

当i＝1,计算上三角矩阵U第j列第i个数据u_i,j等于当前智能电表在t_i时刻采样存储的电能消耗值p_ti,j；

当2<＝i<＝m且i<j，计算上三角矩阵U第j列第i个数据u_i,j等于当前智能电表在t_i时刻采样存储的电能消耗值p_ti,j减去下三角矩阵L的前i列子矩阵和上三角矩阵U的前i行子矩阵的乘积和间所得的差值；

参照上述过程，逐个计算上三角矩阵U第j列的各个数据值。

8.根据权利要求6所述的智能电网非法用电行为分布式检测方法，其特征在于，所述的步骤④中计算上三角矩阵U第m列第m-1个元素值u_m-1,m以及对应住户的诚信参数k_m的具体步骤如下：

1)上三角矩阵U第m列第m-1个元素值u_m-1,m等于当前智能电表在t_m-1时刻采样存储的电能消耗值p_tm-1,m减去下三角矩阵L的前m-1列子矩阵和上三角矩阵U的前m-1行子矩阵的乘积和所得的差值；

2)对应住户的诚信参数k_m的计算：

首先计算t_m时刻变电站处智能电表采样存储的小区电能耗费值P_tm减去下三角矩阵L的前m列子矩阵和上三角矩阵U的前m行子矩阵的乘积和所得的差值D，接着计算差值D减去上三角矩阵U第m+1列到第n列的子矩阵与第m+1个住户到第n个住户的诚信参数向量间的乘积和所得的差值E，最后计算差值E除以上三角矩阵U的第m列第m个数据值所得的商，即为对应住户的诚信参数k_m。