首页 专利交易 科技果 科技人才 科技服务 商标交易 会员权益 IP管家助手 需求市场 关于龙图腾
 /  免费注册
到顶部 到底部
清空 搜索

一种评价智能电能表运行状态的方法和系统 

申请/专利权人:中国电力科学研究院有限公司;国家电网有限公司;国网福建省电力有限公司电力科学研究院

申请日:2019-01-16

公开(公告)日:2024-06-28

公开(公告)号:CN109858785B

主分类号:G06Q10/0637

分类号:G06Q10/0637;G06Q50/06;G06F18/2321

优先权:

专利状态码:有效-授权

法律状态:2024.06.28#授权;2020.10.30#实质审查的生效;2019.06.07#公开

摘要:本发明提供一种评价智能电能表运行状态的方法和系统,所述方法和系统首先基于历史数据构建若干个贝叶斯网,然后通过计算所述贝叶斯网之间的相似度,并中选取若干个精度高、多样性强的贝叶斯网个体,最后采用凸函数证据理论方法对选取的贝叶斯网进行集成,从而得到智能电能表每种运行状态的概率值。所述方法和系统准确地实现了智能电能表的状态评价,为电网的安全运行、状态监测、数据抄读、电价下发等工作提供了极大的便利,有效地节省了定期回收检测电能表的巨大成本,为更换智能电能表提供了有力的技术保障。

主权项:1.一种评价智能电能表运行状态的方法,其特征在于,所述方法包括:步骤1、根据设定的评价智能电能表运行状态的n个指标生成指标集合X={X1…Xi…Xn},采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合F;步骤2、采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ,包括:步骤2-1、构建T×T的相似度矩阵S,其中Si,j表示贝叶斯网Bi和Bj的相似度,当贝叶斯网Bi和Bj各自根据其采样集中每组采样数据确定的运行状态评价结果为yk和y′k时,计算所述贝叶斯网Bi和Bj的相似度Si,j的公式为: 式中,m表示贝叶斯网Bi和Bj的采样集采样的样本的组数;步骤2-2、在集合F中随机选择K个点作为初始质心其中1≤l≤L,1≤p≤K,L为设置的最大迭代次数;步骤2-3、对于集合F中每个贝叶斯网Bi,根据相似度矩阵S,考察Bi与质心的相似度,将Bi分配到与其相似度最高的质心对应的簇中,其中,初始的簇为空集;步骤2-4、对于所述K个簇中的每一个簇,重新计算质心,以使所述质心所处的贝叶斯网与簇内其余贝叶斯网的相似度和最高,并令l=l+1以生成新的质心步骤2-5、当连续两次迭代中,所述K个簇均不发生变化或达到最大迭代次数L时,转至步骤2-6,否则转至步骤2-3继续执行;步骤2-6、从所述K个簇中的每一个簇里选择精度最高的贝叶斯网,以得到新的贝叶斯网集合其中所述精度最高的贝叶斯网是在所述簇中,当质心位于该贝叶斯网时,与簇内其余贝叶斯网的相似度和最高的贝叶斯网;步骤3、采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表每种运行状态的概率值,包括:步骤3-1、对于集合φ中的贝叶斯网φ1,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λ1sl的公式为: 式中,1≤l≤d,1≤i≤n+1,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤3-2、对于集合φ中的贝叶斯网φj,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λjsl的公式为: 式中,2≤j≤K,j的初始值是2,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤3-3、用下式将λ1sj和λjsl进行融合,并将结果赋给λ1sl: 其中Δ1和Δ2是预先设置的常数;步骤3-4、当j<K时,令j=j+1,并转至步骤3-2,当j=K时,转至步骤3-5;步骤3-5、最终得到的λ1sl即为智能电能表每种可能的运行状态的概率值。

全文数据:一种评价智能电能表运行状态的方法和系统技术领域本发明涉及电能计量领域,并且更具体地,涉及一种评价智能电能表运行状态的方法和系统。背景技术智能电能表正逐步进入每个家庭,实现用电信息采集的完全覆盖,其在智能电网的建设中的重要作用正在被凸显出来。智能电能表寿命有限,1、2级有功电能表检定周期一般小于8年,这说明每只电能表使用8年将被收回检测。若仍能继续使用,也将经历“两拆两装”的繁琐流程,再度使用不超过2年时间,故当前回收的智能电能表均采取报废处理的方式,造成很大的资源浪费。若能使用计算机技术对智能电能表的状态进行预测,及时获取智能电能表的状态信息,则可以最大化的对智能电能表进行使用。发明内容为了解决现有技术中对智能电能表的状态无法进行准确评估,致使大量仍能使用的电能表被报废的技术问题,本发明提供一种评价智能电能表运行状态的方法,所述方法包括:根据设定的评价智能电能表运行状态的n个指标生成指标集合X={X1…Xi…Xn},采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合F;采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ;采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表每种运行状态的概率值。进一步地,所述指标集合X选取误差限值X1、安装日期X2、生产厂家X3、到货建档时间X4、运行时长X5、拆除状态X6、是否故障X7、批次号X8、批次数量X9、批次换货数量X10、批次换货率X11、批次故障数量X12、批次故障率X13、实验室误差X14作为评价智能电能表运行状态的指标。进一步地,所述方法在根据设定的评价智能电能表运行状态的n个指标生成指标集合X={X1…Xi…Xn},采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合F之前还包括当所述评价智能电表运行状态的属性指标的取值是连续值时,将所述属性指标的取值离散化,其中:步骤1、对于具有连续取值的指标x,当其在所述历史数据中有k个不重复的连续属性值xi时,计算所述指标x的每个不重复的连续属性值xi出现的概率Pxi,其中1≤i≤k,k的初始值为k0;步骤2、基于所述概率Pxi计算指标x的信息熵步骤3、将所述指标x的不重复的连续属性值划分至k个区间,每个区间对应所述指标x的一个不重复的值,并记录划分点;步骤4、任意选择相邻的两个区间合并,计算合并后的k-1个区间中每个区间出现的概率,并根据所述每个区间的概率计算指标x的信息熵Hk-1,使合并前后的熵之差Hk-Hk-1最小,若出现两组或两组以上相邻区间同时满足条件,则随机选择一组;步骤5、根据合并后的区间的信息熵Hk-1,计算Sj=k0-1×Hk-1-k-2×Hk0,其中,j的初始值为1;步骤6、当Sj≤Sj-1,此时每个区间对应指标x的一个离散值,否则令k=k-1,j=j+1,转步骤4,其中,S0的值为0。进一步地,所述采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合包括:步骤1、对智能电能表运行状态历史数据进行第t次随机采样,共采集m组,得到包含m组样本的采样集Dt,其中,t的初始值是1,1≤t≤T;步骤2、用采样集Dt训练第t个贝叶斯网Bt,所述贝叶斯网Bt包含n+1个变量,其中X1到Xn表示所述评价智能电能表运行状态的n个指标,Xn+1表示电能表状态,有d种可能取值,其中,所述用采样集Dt训练第t个贝叶斯网Bt包括:步骤2-1、将贝叶斯网Bt中包含的n+1个变量作为n+1个节点,令第i个结点Xi的父结点集πi的初始值为空,令式中,ri表示Xi可能取值个数,qi表示Xi的父结点的可能取值个数,Nijk是采样集Dt中Xi取其第k个值,且πi取其第j个值,即时的训练数据个数,Nij是采样集Dt中πi取其第j个值的训练数据个数,i的初始值是1,1≤i≤n+1;步骤2-2、根据公式Z={X1…Xi…Xn+1}-Xi-πi计算集合Z,并使fi,πi∪{Z}达到最大的结点;步骤2-3、令Pnew=fi,πi∪{Z};步骤2-4、当Pnew>Pold,使Pold=Pnew,令πi=πi∪{Z},并转至步骤2-2继续执行,当Pnew<Pold时,πi即为结点Xi的父结点,令i=i+1,并转到步骤2-1,当i=n+1且πn+1即为结点Xn+1的父结点时,转至步骤3;步骤3、根据公式F=F∪{Bt}确定集合F,并在t<T时,令t=t+1,转到步骤1,其中,集合F的初始值为空集。进一步地,所述采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ包括:步骤1、构建T×T的相似度矩阵S,其中Si,j表示贝叶斯网Bi和Bj的相似度,当贝叶斯网Bi和Bj各自根据其采样集中每组采样数据确定的运行状态评价结果为yk和y′k时,计算所述贝叶斯网Bi和Bj的相似度Si,j的公式为:式中,m表示贝叶斯网Bi和Bj的采样集采样的样本的组数;步骤2:在集合F中随机选择K个点作为初始质心其中1≤l≤L,1≤p≤K,L为设置的最大迭代次数;步骤3:对于集合F中每个贝叶斯网Bi,根据相似度矩阵S,考察Bi与质心的相似度,将Bi分配到与其相似度最高的质心对应的簇中,其中,初始的簇为空集;步骤4:对于所述K个簇中的每一个簇,重新计算质心,以使所述质心所处的贝叶斯网与簇内其余贝叶斯网的相似度和最高,并令l=l+1以生成新的质心步骤5:当连续两次迭代中,所述K个簇均不发生变化或达到最大迭代次数L时,转至步骤6,否则转至步骤3继续执行;步骤6:从所述K个簇中的每一个簇里选择精度最高的贝叶斯网,以得到新的贝叶斯网集合其中所述精度最高的贝叶斯网是在所述簇中,当质心位于该贝叶斯网时,与簇内其余贝叶斯网的相似度和最高的贝叶斯网。进一步地,所述采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表各种运行状态的概率值包括:步骤1、对于集合φ中的贝叶斯网φ1,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λ1sl的公式为:式中,1≤l≤d,1≤i≤n+1,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤2、对于集合φ中的贝叶斯网φj,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λjsl的公式为:式中,2≤j≤K,j的初始值是2,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤3、用下式将λ1sj和λjsl进行融合,并将结果赋给λ1sl:其中Δ1和Δ2是预先设置的常数;步骤4、当j<K时,令j=j+1,并转至步骤2,当j=K时,转至步骤5;步骤5、最终得到的λ1sl即为智能电能表每种可能的运行状态的概率值。根据本发明的另一方面,本发明提供一种评价智能电能表运行状态的系统,所述系统包括:第一集合单元,其用于根据设定的评价智能电能表运行状态的n个指标生成指标集合X={X1…Xi…Xn},采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合F;第二集合单元,其用于采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ;状态确定单元,其用于采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表每种运行状态的概率值。进一步地,所述第一集合单元的指标集合X选取误差限值X1、安装日期X2、生产厂家X3、到货建档时间X4、运行时长X5、拆除状态X6、是否故障X7、批次号X8、批次数量X9、批次换货数量X10、批次换货率X11、批次故障数量X12、批次故障率X13、实验室误差X14作为评价智能电能表运行状态的指标。进一步地,所述系统还包括离散化单元,其用于当所述评价智能电表运行状态的属性指标的取值是连续值时,将所述属性指标的取值离散化,其中:步骤1、对于具有连续取值的指标x,当其在所述历史数据中有k个不重复的连续属性值xi时,计算所述指标x的每个不重复的连续属性值xi出现的概率Pxi,其中1≤i≤k,k的初始值为k0;步骤2、基于所述概率Pxi计算指标x的信息熵步骤3、将所述指标x的不重复的连续属性值划分至k个区间,每个区间对应所述指标x的一个不重复的值,并记录划分点;步骤4、任意选择相邻的两个区间合并,计算合并后的k-1个区间中每个区间出现的概率,并根据所述每个区间的概率计算指标x的信息熵Hk-1,使合并前后的熵之差Hk-Hk-1最小,若出现两组或两组以上相邻区间同时满足条件,则随机选择一组;步骤5、根据合并后的区间的信息熵Hk-1,计算Sj=k0-1×Hk-1-k-2×Hk0,其中,j的初始值为1;步骤6、当Sj≤Sj-1,此时每个区间对应指标x的一个离散值,否则令k=k-1,j=j+1,转步骤4,其中,S0的值为0。进一步地,所述第一集合单元采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合包括:步骤1、对智能电能表运行状态历史数据进行第t次随机采样,共采集m组,得到包含m组样本的采样集Dt,其中,t的初始值是1,1≤t≤T;步骤2、用采样集Dt训练第t个贝叶斯网Bt,所述贝叶斯网Bt包含n+1个变量,其中X1到Xn表示所述评价智能电能表运行状态的n个指标,Xn+1表示电能表状态,有d种可能取值,其中,所述用采样集Dt训练第t个贝叶斯网Bt包括:步骤2-1、将贝叶斯网Bt中包含的n+1个变量作为n+1个节点,令第i个结点Xi的父结点集πi的初始值为空,令式中,ri表示Xi可能取值个数,qi表示Xi的父结点的可能取值个数,Nijk是采样集Dt中Xi取其第k个值,且πi取其第j个值,即时的训练数据个数,Nij是采样集Dt中πi取其第j个值的训练数据个数,i的初始值是1,1≤i≤n+1;步骤2-2、根据公式Z={X1…Xi…Xn+1}-Xi-πi计算集合Z,并使fi,πi∪{Z}达到最大的结点;步骤2-3、令Pnew=fi,πi∪{Z};步骤2-4、当Pnew>Pold,使Pold=Pnew,令πi=πi∪{Z},并转至步骤2-2继续执行,当Pnew<Pold时,πi即为结点Xi的父结点,令i=i+1,并转到步骤2-1,当i=n+1且πn+1即为结点Xn+1的父结点时,转至步骤3;步骤3、根据公式F=F∪{Bt}确定集合F,并在t<T时,令t=t+1,转到步骤1,其中,集合F的初始值为空集。进一步地,所述第二集合单元采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ包括:步骤1、构建T×T的相似度矩阵S,其中Si,j表示贝叶斯网Bi和Bj的相似度,当贝叶斯网Bi和Bj各自根据其采样集中每组采样数据确定的运行状态评价结果为yk和y′k时,计算所述贝叶斯网Bi和Bj的相似度Si,j的公式为:式中,m表示贝叶斯网Bi和Bj的采样集采样的样本的组数;步骤2:在集合F中随机选择K个点作为初始质心其中1≤l≤L,1≤p≤K,L为设置的最大迭代次数;步骤3:对于集合F中每个贝叶斯网Bi,根据相似度矩阵S,考察Bi与质心的相似度,将Bi分配到与其相似度最高的质心对应的簇中,其中,初始的簇为空集;步骤4:对于所述K个簇中的每一个簇,重新计算质心,以使所述质心所处的贝叶斯网与簇内其余贝叶斯网的相似度和最高,并令l=l+1以生成新的质心步骤5:当连续两次迭代中,所述K个簇均不发生变化或达到最大迭代次数L时,转至步骤6,否则转至步骤3继续执行;步骤6:从所述K个簇中的每一个簇里选择精度最高的贝叶斯网,以得到新的贝叶斯网集合其中所述精度最高的贝叶斯网是在所述簇中,当质心位于该贝叶斯网时,与簇内其余贝叶斯网的相似度和最高的贝叶斯网。进一步地,所述状态确定单元采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表各种运行状态的概率值包括:步骤1、对于集合φ中的贝叶斯网φ1,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λ1sl的公式为:式中,1≤l≤d,1≤i≤n+1,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤2、对于集合φ中的贝叶斯网φj,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λjsl的公式为:式中,2≤j≤K,j的初始值是2;步骤3、用下式将λ1sj和λjsl进行融合,并将结果赋给λ1sl:其中Δ1和Δ2是预先设置的常数;步骤4、当j<K时,令j=j+1,并转至步骤2,当j=K时,转至步骤5;步骤5、最终得到的λ1sl即为智能电能表每种可能的运行状态的概率值。本发明技术方案提供的评价智能电能表运行状态的方法和系统首先基于历史数据构建若干个贝叶斯网,然后通过计算所述贝叶斯网之间的相似度,并中选取若干个精度高、多样性强的贝叶斯网个体,最后采用凸函数证据理论方法对选取的贝叶斯网进行集成,从而得到智能电能表每种运行状态的概率值。所述方法和系统准确地实现了智能电能表的状态评价,为电网的安全运行、状态监测、数据抄读、电价下发等工作提供了极大的便利,有效地节省了定期回收检测电能表的巨大成本,为更换智能电能表提供了有力的技术保障。附图说明通过参考下面的附图,可以更为完整地理解本发明的示例性实施方式:图1为根据本发明优选实施方式的评价智能电能表运行状态的方法的流程图;图2为根据本发明优选实施方式的评价智能电能表运行状态的系统的结构示意图。具体实施方式现在参考附图介绍本发明的示例性实施方式,然而,本发明可以用许多不同的形式来实施,并且不局限于此处描述的实施例,提供这些实施例是为了详尽地且完全地公开本发明,并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中,相同的单元元件使用相同的附图标记。除非另有说明,此处使用的术语包括科技术语对所属技术领域的技术人员具有通常的理解含义。另外,可以理解的是,以通常使用的词典限定的术语,应当被理解为与其相关领域的语境具有一致的含义,而不应该被理解为理想化的或过于正式的意义。图1为根据本发明优选实施方式的评价智能电能表运行状态的方法的流程图。如图1所示,本优选实施方式所述的评价智能电能表运行状态的方法100从步骤101开始。在步骤101,对预先设置的评价智能电能表运行状态的n个属性指标中具有连续值的属性指标进行离散化。在步骤102,根据设定的评价智能电能表运行状态的n个指标生成指标集合X={X1…Xi…Xn},采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合F。在步骤103,采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ。在步骤104,采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表每种运行状态的概率值。在本优选实施方式中,所述指标集合X选取误差限值X1、安装日期X2、生产厂家X3、到货建档时间X4、运行时长X5、拆除状态X6、是否故障X7、批次号X8、批次数量X9、批次换货数量X10、批次换货率X11、批次故障数量X12、批次故障率X13、实验室误差X14作为评价智能电能表运行状态的指标。优选地,所述对预先设置的评价智能电能表运行状态的n个属性指标中具有连续值的属性指标进行离散化包括:步骤1、对于具有连续取值的指标x,当其在所述历史数据中有k个不重复的连续属性值xi时,计算所述指标x的每个不重复的连续属性值xi出现的概率Pxi,其中1≤i≤k,k的初始值为k0;步骤2、基于所述概率Pxi计算指标x的信息熵步骤3、将所述指标x的不重复的连续属性值划分至k个区间,每个区间对应所述指标x的一个不重复的值,并记录划分点;步骤4、任意选择相邻的两个区间合并,计算合并后的k-1个区间中每个区间出现的概率,并根据所述每个区间的概率计算指标x的信息熵Hk-1,使合并前后的熵之差Hk-Hk-1最小,若出现两组或两组以上相邻区间同时满足条件,则随机选择一组;步骤5、根据合并后的区间的信息熵Hk-1,计算Sj=k0-1×Hk-1-k-2×Hk0,其中,j的初始值为1;步骤6、当Sj≤Sj-1,此时每个区间对应指标x的一个离散值,否则令k=k-1,j=j+1,转步骤4,其中,S0的值为0。优选地,所述采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合包括:步骤1、对智能电能表运行状态历史数据进行第t次随机采样,共采集m组,得到包含m组样本的采样集Dt,其中,t的初始值是1,1≤t≤T;步骤2、用采样集Dt训练第t个贝叶斯网Bt,所述贝叶斯网Bt包含n+1个变量,其中X1到Xn表示所述评价智能电能表运行状态的n个指标,Xn+1表示电能表状态,有d种可能取值,其中,所述用采样集Dt训练第t个贝叶斯网Bt包括:步骤2-1、将贝叶斯网Bt中包含的n+1个变量作为n+1个节点,令第i个结点Xi的父结点集πi的初始值为空,令式中,ri表示Xi可能取值个数,qi表示Xi的父结点的可能取值个数,Nijk是采样集Dt中Xi取其第k个值,且πi取其第j个值,即时的训练数据个数,Nij是采样集Dt中πi取其第j个值的训练数据个数,i的初始值是1,1≤i≤n+1;步骤2-2、根据公式Z={X1…Xi…Xn+1}-Xi-πi计算集合Z,并使fi,πi∪{Z}达到最大的结点;步骤2-3、令Pnew=fi,πi∪{Z};步骤2-4、当Pnew>Pold,使Pold=Pnew,令πi=πi∪{Z},并转至步骤2-2继续执行,当Pnew<Pold时,πi即为结点Xi的父结点,令i=i+1,并转到步骤2-1,当i=n+1且πn+1即为结点Xn+1的父结点时,转至步骤3;步骤3、根据公式F=F∪{Bt}确定集合F,并在t<T时,令t=t+1,转到步骤1,其中,集合F的初始值为空集。优选地,所述采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ包括:步骤1、构建T×T的相似度矩阵S,其中Si,j表示贝叶斯网Bi和Bj的相似度,当贝叶斯网Bi和Bj各自根据其采样集中每组采样数据确定的运行状态评价结果为yk和y′k时,计算所述贝叶斯网Bi和Bj的相似度Si,j的公式为:式中,m表示贝叶斯网Bi和Bj的采样集采样的样本的组数;步骤2:在集合F中随机选择K个点作为初始质心其中1≤l≤L,1≤p≤K,L为设置的最大迭代次数;步骤3:对于集合F中每个贝叶斯网Bi,根据相似度矩阵S,考察Bi与质心的相似度,将Bi分配到与其相似度最高的质心对应的簇中,其中,初始的簇为空集;步骤4:对于所述K个簇中的每一个簇,重新计算质心,以使所述质心所处的贝叶斯网与簇内其余贝叶斯网的相似度和最高,并令l=l+1以生成新的质心步骤5:当连续两次迭代中,所述K个簇均不发生变化或达到最大迭代次数L时,转至步骤6,否则转至步骤3继续执行;步骤6:从所述K个簇中的每一个簇里选择精度最高的贝叶斯网,以得到新的贝叶斯网集合其中所述精度最高的贝叶斯网是在所述簇中,当质心位于该贝叶斯网时,与簇内其余贝叶斯网的相似度和最高的贝叶斯网。优选地,所述采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表各种运行状态的概率值包括:步骤1、对于集合φ中的贝叶斯网φ1,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λ1sl的公式为:式中,1≤l≤d,1≤i≤n+1,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤2、对于集合φ中的贝叶斯网φj,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λjsl的公式为:式中,2≤j≤K,j的初始值是2,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤3、用下式将λ1sj和λjsl进行融合,并将结果赋给λ1sl:其中Δ1和Δ2是预先设置的常数;步骤4、当j<K时,令j=j+1,并转至步骤2,当j=K时,转至步骤5;步骤5、最终得到的λ1sl即为智能电能表每种可能的运行状态的概率值。图2为根据本发明优选实施方式的评价智能电能表运行状态的系统的结构示意图。如图2所示,本优选实施方式所述的评价智能电能表运行状态的系统200包括:离散化单元201,其用于对预先设置的评价智能电能表运行状态的n个属性指标中具有连续值的属性指标进行离散化。第一集合单元202,其用于根据设定的评价智能电能表运行状态的n个指标生成指标集合X={X1…Xi…Xn},采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合F。第二集合单元203,其用于采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ;状态确定单元204,其用于采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表每种运行状态的概率值。优选地,所述第一集合单元202的指标集合X选取误差限值X1、安装日期X2、生产厂家X3、到货建档时间X4、运行时长X5、拆除状态X6、是否故障X7、批次号X8、批次数量X9、批次换货数量X10、批次换货率X11、批次故障数量X12、批次故障率X13、实验室误差X14作为评价智能电能表运行状态的指标。优选地,所述离散化单元201对预先设置的评价智能电能表运行状态的n个属性指标中具有连续值的属性指标进行离散化包括:步骤1、对于具有连续取值的指标x,当其在所述历史数据中有k个不重复的连续属性值xi时,计算所述指标x的每个不重复的连续属性值xi出现的概率Pxi,其中1≤i≤k,k的初始值为k0;步骤2、基于所述概率Pxi计算指标x的信息熵步骤3、将所述指标x的不重复的连续属性值划分至k个区间,每个区间对应所述指标x的一个不重复的值,并记录划分点;步骤4、任意选择相邻的两个区间合并,计算合并后的k-1个区间中每个区间出现的概率,并根据所述每个区间的概率计算指标x的信息熵Hk-1,使合并前后的熵之差Hk-Hk-1最小,若出现两组或两组以上相邻区间同时满足条件,则随机选择一组;步骤5、根据合并后的区间的信息熵Hk-1,计算Sj=k0-1×Hk-1-k-2×Hk0,其中,j的初始值为1;步骤6、当Sj≤Sj-1,此时每个区间对应指标x的一个离散值,否则令k=k-1,j=j+1,转步骤4,其中,S0的值为0。优选地,所述第一集合单元202采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合包括:步骤1、对智能电能表运行状态历史数据进行第t次随机采样,共采集m组,得到包含m组样本的采样集Dt,其中,t的初始值是1,1≤t≤T;步骤2、用采样集Dt训练第t个贝叶斯网Bt,所述贝叶斯网Bt包含n+1个变量,其中X1到Xn表示所述评价智能电能表运行状态的n个指标,Xn+1表示电能表状态,有d种可能取值,其中,所述用采样集Dt训练第t个贝叶斯网Bt包括:步骤2-1、将贝叶斯网Bt中包含的n+1个变量作为n+1个节点,令第i个结点Xi的父结点集πi的初始值为空,令式中,ri表示Xi可能取值个数,qi表示Xi的父结点的可能取值个数,Nijk是采样集Dt中Xi取其第k个值,且πi取其第j个值,即时的训练数据个数,Nij是采样集Dt中πi取其第j个值的训练数据个数,i的初始值是1,1≤i≤n+1;步骤2-2、根据公式Z={X1…Xi…Xn+1}-Xi-πi计算集合Z,并使fi,πi∪{Z}达到最大的结点;步骤2-3、令Pnew=fi,πi∪{Z};步骤2-4、当Pnew>Pold,使Pold=Pnew,令πi=πi∪{Z},并转至步骤2-2继续执行,当Pnew<Pold时,πi即为结点Xi的父结点,令i=i+1,并转到步骤2-1,当i=n+1且πn+1即为结点Xn+1的父结点时,转至步骤3;步骤3、根据公式F=F∪{Bt}确定集合F,并在t<T时,令t=t+1,转到步骤1,其中,集合F的初始值为空集。优选地,所述第二集合单元203采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ包括:步骤1、构建T×T的相似度矩阵S,其中Si,j表示贝叶斯网Bi和Bj的相似度,当贝叶斯网Bi和Bj各自根据其采样集中每组采样数据确定的运行状态评价结果为yk和y′k时,计算所述贝叶斯网Bi和Bj的相似度Si,j的公式为:式中,m表示贝叶斯网Bi和Bj的采样集采样的样本的组数;步骤2:在集合F中随机选择K个点作为初始质心其中1≤l≤L,1≤p≤K,L为设置的最大迭代次数;步骤3:对于集合F中每个贝叶斯网Bi,根据相似度矩阵S,考察Bi与质心的相似度,将Bi分配到与其相似度最高的质心对应的簇中,其中,初始的簇为空集;步骤4:对于所述K个簇中的每一个簇,重新计算质心,以使所述质心所处的贝叶斯网与簇内其余贝叶斯网的相似度和最高,并令l=l+1以生成新的质心步骤5:当连续两次迭代中,所述K个簇均不发生变化或达到最大迭代次数L时,转至步骤6,否则转至步骤3继续执行;步骤6:从所述K个簇中的每一个簇里选择精度最高的贝叶斯网,以得到新的贝叶斯网集合其中所述精度最高的贝叶斯网是在所述簇中,当质心位于该贝叶斯网时,与簇内其余贝叶斯网的相似度和最高的贝叶斯网。优选地,所述状态确定单元204采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表各种运行状态的概率值包括:步骤1、对于集合φ中的贝叶斯网φ1,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λ1sl的公式为:式中,1≤l≤d,1≤i≤n+1,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤2、对于集合φ中的贝叶斯网φj,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λjsl的公式为:式中,2≤j≤K,j的初始值是2;步骤3、用下式将λ1sj和λjsl进行融合,并将结果赋给λ1sl:其中Δ1和Δ2是预先设置的常数;步骤4、当j<K时,令j=j+1,并转至步骤2,当j=K时,转至步骤5;步骤5、最终得到的λ1sl即为智能电能表每种可能的运行状态的概率值。已经通过参考少量实施方式描述了本发明。然而,本领域技术人员所公知的,正如附带的专利权利要求所限定的,除了本发明以上公开的其他的实施例等同地落在本发明的范围内。通常地,在权利要求中使用的所有术语都根据他们在技术领域的通常含义被解释,除非在其中被另外明确地定义。所有的参考“一个所述该[装置、组件等]”都被开放地解释为所述装置、组件等中的至少一个实例,除非另外明确地说明。这里公开的任何方法的步骤都没必要以公开的准确的顺序运行,除非明确地说明。

权利要求:1.一种评价智能电能表运行状态的方法,其特征在于,所述方法包括:根据设定的评价智能电能表运行状态的n个指标生成指标集合X={X1…Xi…Xn},采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合F;采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ;采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表每种运行状态的概率值。2.根据权利要求1所述的方法,其特征在于,所述指标集合X选取误差限值X1、安装日期X2、生产厂家X3、到货建档时间X4、运行时长X5、拆除状态X6、是否故障X7、批次号X8、批次数量X9、批次换货数量X10、批次换货率X11、批次故障数量X12、批次故障率X13、实验室误差X14作为评价智能电能表运行状态的指标。3.根据权利要求2所述的方法,其特征在于,所述方法在根据设定的评价智能电能表运行状态的n个指标生成指标集合X={X1…Xi…Xn},采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合F之前还包括当所述评价智能电表运行状态的属性指标的取值是连续值时,将所述属性指标的取值离散化,其中:步骤1、对于具有连续取值的指标x,当其在所述历史数据中有k个不重复的连续属性值xi时,计算所述指标x的每个不重复的连续属性值xi出现的概率Pxi,其中1≤i≤k,k的初始值为k0;步骤2、基于所述概率Pxi计算指标x的信息熵步骤3、将所述指标x的不重复的连续属性值划分至k个区间,每个区间对应所述指标x的一个不重复的值,并记录划分点;步骤4、任意选择相邻的两个区间合并,计算合并后的k-1个区间中每个区间出现的概率,并根据所述每个区间的概率计算指标x的信息熵Hk-1,使合并前后的熵之差Hk-Hk-1最小,若出现两组或两组以上相邻区间同时满足条件,则随机选择一组;步骤5、根据合并后的区间的信息熵Hk-1,计算Sj=k0-1×Hk-1-k-2×Hk0,其中,j的初始值为1;步骤6、当Sj≤Sj-1,此时每个区间对应指标x的一个离散值,否则令k=k-1,j=j+1,转步骤4,其中,S0的值为0。4.根据权利要求3所述的方法,其特征在于,所述采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合包括:步骤1、对智能电能表运行状态历史数据进行第t次随机采样,共采集m组,得到包含m组样本的采样集Dt,其中,t的初始值是1,1≤t≤T;步骤2、用采样集Dt训练第t个贝叶斯网Bt,所述贝叶斯网Bt包含n+1个变量,其中X1到Xn表示所述评价智能电能表运行状态的n个指标,Xn+1表示电能表状态,有d种可能取值,其中,所述用采样集Dt训练第t个贝叶斯网Bt包括:步骤2-1、将贝叶斯网Bt中包含的n+1个变量作为n+1个节点,令第i个结点Xi的父结点集πi的初始值为空,令式中,ri表示Xi可能取值个数,qi表示Xi的父结点的可能取值个数,Nijk是采样集Dt中Xi取其第k个值,且πi取其第j个值,即时的训练数据个数,Nij是采样集Dt中πi取其第j个值的训练数据个数,i的初始值是1,1≤i≤n+1;步骤2-2、根据公式Z={X1…Xi…Xn+1}-Xi-πi计算集合Z,并使fi,πi∪{Z}达到最大的结点;步骤2-3、令Pnew=fi,πi∪{Z};步骤2-4、当Pnew>Pold,使Pold=Pnew,令πi=πi∪{Z},并转至步骤2-2继续执行,当Pnew<Pold时,πi即为结点Xi的父结点,令i=i+1,并转到步骤2-1,当i=n+1且πn+1即为结点Xn+1的父结点时,转至步骤3;步骤3、根据公式F=F∪{Bt}确定集合F,并在t<T时,令t=t+1,转到步骤1,其中,集合F的初始值为空集。5.根据权利要求4所述的方法,其特征在于,所述采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ包括:步骤1、构建T×T的相似度矩阵S,其中Si,j表示贝叶斯网Bi和Bj的相似度,当贝叶斯网Bi和Bj各自根据其采样集中每组采样数据确定的运行状态评价结果为yk和y′k时,计算所述贝叶斯网Bi和Bj的相似度Si,j的公式为:式中,m表示贝叶斯网Bi和Bj的采样集采样的样本的组数;步骤2、在集合F中随机选择K个点作为初始质心其中1≤l≤L,1≤p≤K,L为设置的最大迭代次数;步骤3、对于集合F中每个贝叶斯网Bi,根据相似度矩阵S,考察Bi与质心的相似度,将Bi分配到与其相似度最高的质心对应的簇中,其中,初始的簇为空集;步骤4、对于所述K个簇中的每一个簇,重新计算质心,以使所述质心所处的贝叶斯网与簇内其余贝叶斯网的相似度和最高,并令l=l+1以生成新的质心步骤5、当连续两次迭代中,所述K个簇均不发生变化或达到最大迭代次数L时,转至步骤6,否则转至步骤3继续执行;步骤6、从所述K个簇中的每一个簇里选择精度最高的贝叶斯网,以得到新的贝叶斯网集合其中所述精度最高的贝叶斯网是在所述簇中,当质心位于该贝叶斯网时,与簇内其余贝叶斯网的相似度和最高的贝叶斯网。6.根据权利要求5所述的方法,其特征在于,所述采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表各种运行状态的概率值包括:步骤1、对于集合φ中的贝叶斯网φ1,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λ1sl的公式为:式中,1≤l≤d,1≤i≤n+1,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤2、对于集合φ中的贝叶斯网φj,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λjsl的公式为:式中,2≤j≤K,j的初始值是2,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤3、用下式将λ1sj和λjsl进行融合,并将结果赋给λ1sl:其中Δ1和Δ2是预先设置的常数;步骤4、当j<K时,令j=j+1,并转至步骤2,当j=K时,转至步骤5;步骤5、最终得到的λ1sl即为智能电能表每种可能的运行状态的概率值。7.一种评价智能电能表运行状态的系统,其特征在于,所述系统包括:第一集合单元,其用于根据设定的评价智能电能表运行状态的n个指标生成指标集合X={X1…Xi…Xn},采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合F;第二集合单元,其用于采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ;状态确定单元,其用于采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表每种运行状态的概率值。8.根据权利要求7所述的系统,其特征在于,所述第一集合单元的指标集合X选取误差限值X1、安装日期X2、生产厂家X3、到货建档时间X4、运行时长X5、拆除状态X6、是否故障X7、批次号X8、批次数量X9、批次换货数量X10、批次换货率X11、批次故障数量X12、批次故障率X13、实验室误差X14作为评价智能电能表运行状态的指标。9.根据权利要求8所述的系统,其特征在于,所述系统还包括离散化单元,其用于当所述评价智能电表运行状态的属性指标的取值是连续值时,将所述属性指标的取值离散化,其中:步骤1、对于具有连续取值的指标x,当其在所述历史数据中有k个不重复的连续属性值xi时,计算所述指标x的每个不重复的连续属性值xi出现的概率Pxi,其中1≤i≤k,k的初始值为k0;步骤2、基于所述概率Pxi计算指标x的信息熵步骤3、将所述指标x的不重复的连续属性值划分至k个区间,每个区间对应所述指标x的一个不重复的值,并记录划分点;步骤4、任意选择相邻的两个区间合并,计算合并后的k-1个区间中每个区间出现的概率,并根据所述每个区间的概率计算指标x的信息熵Hk-1,使合并前后的熵之差Hk-Hk-1最小,若出现两组或两组以上相邻区间同时满足条件,则随机选择一组;步骤5、根据合并后的区间的信息熵Hk-1,计算Sj=k0-1×Hk-1-k-2×Hk0,其中,j的初始值为1;步骤6、当Sj≤Sj-1,此时每个区间对应指标x的一个离散值,否则令k=k-1,j=j+1,转步骤4,其中,S0的值为0。10.根据权利要求9所述的系统,其特征在于,所述第一集合单元采用bagging算法从智能电能表运行状态历史数据中抽样所述指标集合X的数据集并训练生成T个贝叶斯网的集合包括:步骤1、对智能电能表运行状态历史数据进行第t次随机采样,共采集m组,得到包含m组样本的采样集Dt,其中,t的初始值是1,1≤t≤T;步骤2、用采样集Dt训练第t个贝叶斯网Bt,所述贝叶斯网Bt包含n+1个变量,其中X1到Xn表示所述评价智能电能表运行状态的n个指标,Xn+1表示电能表状态,有d种可能取值,其中,所述用采样集Dt训练第t个贝叶斯网Bt包括:步骤2-1、将贝叶斯网Bt中包含的n+1个变量作为n+1个节点,令第i个结点Xi的父结点集πi的初始值为空,令式中,ri表示Xi可能取值个数,qi表示Xi的父结点的可能取值个数,Nijk是采样集Dt中Xi取其第k个值,且πi取其第j个值,即时的训练数据个数,Nij是采样集Dt中πi取其第j个值的训练数据个数,i的初始值是1,1≤i≤n+1;步骤2-2、根据公式Z={X1…Xi…Xn+1}-Xi-πi计算集合Z,并使fi,πi∪{Z}达到最大的结点;步骤2-3、令Pnew=fi,πi∪{Z};步骤2-4、当Pnew>Pold,使Pold=Pnew,令πi=πi∪{Z},并转至步骤2-2继续执行,当Pnew<Pold时,πi即为结点Xi的父结点,令i=i+1,并转到步骤2-1,当i=n+1且πn+1即为结点Xn+1的父结点时,转至步骤3;步骤3、根据公式F=F∪{Bt}确定集合F,并在t<T时,令t=t+1,转到步骤1,其中,集合F的初始值为空集。11.根据权利要求10所述的系统,其特征在于,所述第二集合单元采用聚类算法从集合F中的T个贝叶斯网中选择K个精度最高的贝叶斯网组成新的集合φ包括:步骤1、构建T×T的相似度矩阵S,其中Si,j表示贝叶斯网Bi和Bj的相似度,当贝叶斯网Bi和Bj各自根据其采样集中每组采样数据确定的运行状态评价结果为yk和y′k时,计算所述贝叶斯网Bi和Bj的相似度Si,j的公式为:式中,m表示贝叶斯网Bi和Bj的采样集采样的样本的组数;步骤2、在集合F中随机选择K个点作为初始质心其中1≤l≤L,1≤p≤K,L为设置的最大迭代次数;步骤3、对于集合F中每个贝叶斯网Bi,根据相似度矩阵S,考察Bi与质心的相似度,将Bi分配到与其相似度最高的质心对应的簇中,其中,初始的簇为空集;步骤4、对于所述K个簇中的每一个簇,重新计算质心,以使所述质心所处的贝叶斯网与簇内其余贝叶斯网的相似度和最高,并令l=l+1以生成新的质心步骤5、当连续两次迭代中,所述K个簇均不发生变化或达到最大迭代次数L时,转至步骤6,否则转至步骤3继续执行;步骤6、从所述K个簇中的每一个簇里选择精度最高的贝叶斯网,以得到新的贝叶斯网集合其中所述精度最高的贝叶斯网是在所述簇中,当质心位于该贝叶斯网时,与簇内其余贝叶斯网的相似度和最高的贝叶斯网。12.根据权利要求11所述的系统,其特征在于,所述状态确定单元采用凸函数证据理论方法,对集合φ中的K个贝叶斯网进行集成以确定智能电能表各种运行状态的概率值包括:步骤1、对于集合φ中的贝叶斯网φ1,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λ1sl的公式为:式中,1≤l≤d,1≤i≤n+1,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点,Xi∈ChildrenY表示电能表的属性指标Xi是运行状态Y的子结点,πY是运行状态Y的父结点;步骤2、对于集合φ中的贝叶斯网φj,当x1,...,xn表示实际观测到的评价智能电能表运行状态的n个指标的数值,Yl表示智能电能表的运行状态d种可能取值中的第l种,则计算智能电能表的运行状态Yl的概率λjsl的公式为:式中,2≤j≤K,j的初始值是2;步骤3、用下式将λ1sj和λjsl进行融合,并将结果赋给λ1sl:其中Δ1和Δ2是预先设置的常数;步骤4、当j<K时,令j=j+1,并转至步骤2,当j=K时,转至步骤5;步骤5、最终得到的λ1sl即为智能电能表每种可能的运行状态的概率值。

百度查询: 中国电力科学研究院有限公司;国家电网有限公司;国网福建省电力有限公司电力科学研究院 一种评价智能电能表运行状态的方法和系统

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。