一种基于MK突变检验的用电异常分析方法与流程

一种基于mk突变检验的用电异常分析方法
技术领域
1.本发明属于用电异常检测技术领域，具体涉及一种基于mk突变检验的用电异常分析方法。


背景技术：

2.随着经济发展，用电需求日益增长，窃电时有发生。窃电不仅损害供电企业的经济效益，更对电网安全带来严重影响。然而，窃电方式、窃电技术的多样化，给实地查处带来了巨大的挑战。
3.各相关单位对用电异常分析进行了大量研究，开发了多种不同的用电异常分析方法。然而，这些方法大都作为理论参考，难以实际应用。本发明旨在提供一种不同于现有的用电异常分析方法，可操作性好，检测准确。


技术实现要素：

4.本发明针对现有用电异常分析方法可操作性较差或者数据处理量大或者准确度较低的不足，提供一种基于mk突变检验的用电异常分析方法，可操作性好，检测准确。
5.为实现上述目的，本发明采用如下技术方案：一种基于mk突变检验的用电异常分析方法，所述基于mk突变检验的用电异常分析方法包括：
6.步骤s1、获取待分析用户的数据，所需数据包括用电数据和台区线损数据；
7.步骤s2、异常数据处理，筛选用户；
8.步骤s3、分类输出单项表和三相表用电异常可疑用户；
9.步骤s4、采用mk突变检验法得到可疑用户用电的突变日期；
10.步骤s5、采用mk突变检验法得到台区线损的突变日期；
11.步骤s6、根据可疑用户用电的突变日期和台区线损的突变日期，输出用电异常用户名单。
12.本发明的基于mk突变检验的用电异常分析方法，结合台区线损及用户用电量变化，运用mk突变检验法，可以输出窃电用户名单，可操作性好，结果较为准确；所需数据种类较少、数据处理要求较低；先筛选用户再以一定标准输出单相表和三相表用电异常可疑用户，需要处理的数据相对较少，可以提高处理速度。其中，部分步骤的先后顺序可以根据需要相应调整。
13.作为改进，步骤s1中，选择非光伏用户作为待分析用户，获取非光伏用户前若干天的负荷数据、前一年用电量以及台区线损数据。光伏用户火线电流是负的，光伏用户的数据表现和窃电用户的数据表现是一样的，会影响准确率，所以剔除掉光伏用户。
14.作为改进，步骤s2中，对零/火线电流长期为负值的情况，判断为由单相表零线火线进出反接导致，将零/火线电流转变为电流的绝对值，对这部分用户进行后续分析。
15.作为改进，步骤s2中，查找同一台区下零线电流与火线电流不均衡次数占比超过μ的单相表表计，对此类表计进行同台区下同一时刻点不同表计对应负荷数据比较，若存在
两个计量表，使得：(表计一的零线电流减去表计二的火线电流)的绝对值《k且(表计二的零线电流减去表计一的火线电流)的绝对值《k，满足该条件的时刻点占比值超过λ，则认为是由两个计量表零火线相互反接所致，将此类用户作为正常用户，从待分析用户清单中剔除，其中，μ、k和λ为相关阈值，可通过理论和/或实际数据得到。
16.作为改进，步骤s2中，提取单相表中上午a时至下午b时负荷数据，若a相电流与零线电流均为负值，且该表计下午c时至凌晨d时的a相电流与零线电流均非负，则标记此类用户，并从待分析用户清单中剔除。a与d相邻，b与c相邻。
17.作为改进，步骤s3中，对非光伏的单相表用户，筛选出前若干天零线电流减去火线电流＞k的记录，标记为电异常可疑用户，根据异常数占比情况，输出单相表用电异常可疑用户名单，其中，k为相关阈值。
18.作为改进，步骤s3中，三相表用电异常可疑用户输出过程如下：
19.步骤s31、计算用户用电量的平均值avg_use，设定阈值k1、k2，比较用户用电量与k1*avg_use、k2*avg_use的关系，输出用户用电量向量；
20.步骤s32、计算台区线损电量平均值为avg_loss，设置阈值k3、k4，比较台区线损电量与k3*avg_loss、k4*avg_loss的关系，输出台区线损向量；
21.步骤s33、计算用户用电量向量与台区线损向量在多维空间的夹角cos值；
22.步骤s34、通过空间向量夹角cos值，得到用户用电量变化与线损变化相反的记录数，当记录数大于相关阈值时，认为具有用电异常可疑，输出三相表用电异常可疑用户名单。
23.作为改进，步骤s4中，过程如下，包括：
24.步骤s41、根据可疑用户用电量时间序列数据构建秩序列，sk表示第i个时刻点大于j时刻点数值个数累计值，其中：
25.步骤s42、计算用电量顺序时间序列的秩序列ufk，，其中：uf1＝0,e(sk)为sk的平均值，var(sk)为sk的方差；
26.当x1，x2
…
xn相互独立且具有相同连续分布时，xn相互独立且具有相同连续分布时，
27.步骤s43、计算用电量逆序时间序列的秩序列ubk，ubk＝-ufk，k＝n,n-1,
…
,1，其中，ub1＝0；
28.步骤s44、确定显著性水平α，绘制ubk、ufk曲线以及临界值线，ubk，ufk两条曲线的交点即为突变开始日期，若在临界值线范围内，则交叉点对应的突变点显著，结合ufk曲线变化情况，输出用户用电量下降的突变日期名单。
29.作为改进，步骤s5中，采用权利要求8中的方法，输出台区线损升高的突变日期名单。
30.作为改进，步骤s6中，对比用户用电量mk突变点与台区线损mk突变点，判断台区线损突变点时刻点与该台区下用户用电量突变点时刻点是否接近，并根据台区线损电量在该时刻点前后趋于上升趋势，用户用电量趋于下降趋势,确认该用户的窃电可疑，最终输出用
电异常用户名单。
31.本发明的基于mk突变检验的用电异常分析方法的有益效果是：结合台区线损及用户用电量变化，运用mk突变检验法，可以输出窃电用户名单，可操作性好，结果较为准确；所需数据种类较少、数据处理要求较低；先筛选用户再以一定标准输出单相表和三相表用电异常可疑用户，需要处理的数据相对较少，可以提高处理速度。
附图说明
32.图1是本发明实施例一中应用本发明实施例一的基于mk突变检验的用电异常分析方法得到的某用电异常用户所在台区线损电量的ubk、ufk曲线(横坐标为时间，单位天，纵坐标为ubk、ufk值)。
33.图2是应用本发明实施例一的基于mk突变检验的用电异常分析方法得到的某用电异常用户用电量的ubk、ufk曲线(横坐标为时间，单位天，纵坐标为ubk、ufk值)。
34.图3是对该用电异常用户处理后所在台区线损电量的ubk、ufk曲线(横坐标为时间，单位天，纵坐标为ubk、ufk值)。
35.图4是对该用电异常用户处理后用电量的ubk、ufk曲线(横坐标为时间，单位天，纵坐标为ubk、ufk值)。
具体实施方式
36.下面结合本发明创造实施例的附图，对本发明创造实施例的技术方案进行解释和说明，但下述实施例仅为本发明创造的优选实施例，并非全部。基于实施方式中的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的其他实施例，都属于本发明创造的保护范围。
37.参见图1和图2，本发明的一种基于mk突变检验的用电异常分析方法，所述基于mk突变检验的用电异常分析方法包括：
38.步骤s1、获取待分析用户的数据，所需数据包括用电数据和台区线损数据；
39.步骤s2、异常数据处理，筛选用户；
40.步骤s3、分类输出单项表和三相表用电异常可疑用户；
41.步骤s4、采用mk突变检验法得到可疑用户用电的突变日期；
42.步骤s5、采用mk突变检验法得到台区线损的突变日期；
43.步骤s6、根据可疑用户用电的突变日期和台区线损的突变日期，输出用电异常用户名单。
44.本发明的基于mk突变检验的用电异常分析方法，结合台区线损及用户用电量变化，运用mk突变检验法，可以输出窃电用户名单，可操作性好，结果较为准确；所需数据种类较少、数据处理要求较低；先筛选用户再以一定标准输出单相表和三相表用电异常可疑用户，需要处理的数据相对较少，可以提高处理速度。
45.实施例一
46.参见图1和图2，本发明实施例一的一种基于mk突变检验的用电异常分析方法，所述基于mk突变检验的用电异常分析方法包括：
47.步骤s1、获取待分析用户的数据，所需数据包括用电数据和台区线损数据；
48.步骤s2、异常数据处理，筛选用户；
49.步骤s3、分类输出单项表和三相表用电异常可疑用户；
50.步骤s4、采用mk突变检验法得到可疑用户用电的突变日期；
51.步骤s5、采用mk突变检验法得到台区线损的突变日期；
52.步骤s6、根据可疑用户用电的突变日期和台区线损的突变日期，输出用电异常用户名单。
53.本实施例中，步骤s1中，选择非光伏用户作为待分析用户，获取非光伏用户前若干天的负荷数据、前一年用电量以及台区线损数据。
54.本实施例中，步骤s2中，对零/火线电流长期为负值的情况，判断为由单相表零线火线进出反接导致，将零/火线电流转变为电流的绝对值。
55.本实施例中，步骤s2中，查找同一台区下零线电流与火线电流不均衡次数占比超过μ的单相表表计，对此类表计进行同台区下同一时刻点不同表计对应负荷数据比较，若存在两个计量表，使得：(表计一的零线电流减去表计二的火线电流)的绝对值《k 且(表计二的零线电流减去表计一的火线电流)的绝对值《k，满足该条件的时刻点占比值超过λ，则认为是由两个计量表零火线相互反接所致，将此类用户作为正常用户，从待分析用户清单中剔除；其中，μ、k和λ为相关阈值。
56.本实施例中，步骤s2中，提取单相表中上午a时至下午b时负荷数据，若a相电流与零线电流均为负值，且该表计下午c时至凌晨d时的a相电流与零线电流均非负，则标记此类用户，并从待分析用户清单中剔除。
57.本实施例中，步骤s3中，对非光伏的单相表用户，筛选出前若干天零线电流减去火线电流＞k的记录，标记为电异常可疑用户，根据异常数占比情况，输出单相表用电异常可疑用户名单，其中，k为相关阈值。
58.本实施例中，步骤s3中，三相表用电异常可疑用户输出过程如下：
59.步骤s31、计算用户用电量的平均值avg_use，设定阈值k1、k2，比较用户用电量与k1*avg_use、k2*avg_use的关系，输出用户用电量向量；
60.步骤s32、计算台区线损电量平均值为avg_loss，设置阈值k3、k4，比较台区线损电量与k3*avg_loss、k4*avg_loss的关系，输出台区线损向量；
61.步骤s33、计算用户用电量向量与台区线损向量在多维空间的夹角cos值；
62.步骤s34、通过空间向量夹角cos值，得到用户用电量变化与线损变化相反的记录数，输出三相表用电异常可疑用户名单。
63.本实施例中，步骤s4中，过程如下，包括：
64.步骤s41、根据可疑用户用电量时间序列数据构建秩序列，sk表示第i个时刻点大于j时刻点数值个数累计值，其中：
65.步骤s42、计算用电量顺序时间序列的秩序列ufk，，其中：uf1＝0,e(sk)为sk的平均值，var(sk)为sk的方差；
66.当x1，x2
…
xn相互独立且具有相同连续分布时，
67.步骤s43、计算用电量逆序时间序列的秩序列ubk，ubk＝-ufk，k＝n,n-1,
…
,1，其中，ub1＝0；
68.步骤s44、确定显著性水平α，绘制ubk、ufk曲线以及临界值线，ubk，ufk两条曲线的交点即为突变开始日期，若在临界值线范围内，则交叉点对应的突变点显著，结合ufk曲线变化情况，输出用户用电量下降的突变日期名单。
69.本实施例中，步骤s5中，采用权利要求8中的方法，输出台区线损升高的突变日期名单。
70.本实施例中，步骤s6中，对比用户用电量mk突变点与台区线损mk突变点，判断台区线损突变点时刻点与该台区下用户用电量突变点时刻点是否接近，并根据台区线损电量在该时刻点前后趋于上升趋势，用户用电量趋于下降趋势,确认该用户的窃电可疑，最终输出用电异常用户名单。
71.参见图1，某用户所在台区线损电量的ubk、ufk曲线以及临界值线。，通过ubk,ufk曲线交点位置确定突变点及显著性。结合ufk变化趋势，输出台区线损升高的突变日期名单。
72.参见图2某用户用电量的ubk、ufk曲线以及临界值线。ubk、ufk两条曲线的交点即为突变开始日期，若在临界线范围内，则交叉点对应的突变点显著。结合ufk曲线变化情况，输出用户用电量下降的突变日期名单。
73.从图1和图2可以看出，台区线损突变点时刻点与该台区下用户用电量突变点时刻点相近，并且台区线损电量在该时刻点前后趋于上升趋势，用户用电量趋于下降趋势。因此，确认该用户具有可疑，将该用户列入异常用户名单，进行后续处理。
74.参见图3和图4，对该用电异常用户进行处理后，台区线损突变点时刻点与该台区下用户用电量突变点时刻点相近，并且台区线损电量在该时刻点前后趋于下降趋势，用户用电量趋于上升趋势。因此，确认该用户已无异常。
75.本发明实施例一的基于mk突变检验的用电异常分析方法的有益效果是：结合台区线损及用户用电量变化，运用mk突变检验法，可以输出窃电用户名单，可操作性好，结果较为准确；所需数据种类较少、数据处理要求较低；先筛选用户再以一定标准输出单相表和三相表用电异常可疑用户，需要处理的数据相对较少，可以提高处理速度。
76.以上所述，仅为本发明创造的具体实施方式，但本发明创造的保护范围并不局限于此，熟悉该本领域的技术人员应该明白本发明创造包括但不限于上面具体实施方式中描述的内容。任何不偏离本发明创造的功能和结构原理的修改都将包括在权利要求书的范围中。