一种基于时空卷积神经网络的智能电表故障预测方法与流程

1.本发明主要涉及智能电表故障预测的技术领域，具体为一种基于时空卷积神经网络的智能电表故障预测方法。


背景技术：

2.目前，国家电网公司等系统内部主要依靠现场运维人员，通过周期性的巡检对现场的智能电表故障情况进行统计，该方法耗时耗力，并且不能及时发现智能电能表故障的情况。
3.部分学者展开了基于机器学习的方法对智能电能表故障进行预测，主要包括：
4.(1)基于无监督的聚类算法，如k-means聚类，pca等，该类算法的优点在于无需进行人为的标签添加，根据原始数据的空间分布情况就能进行分类，将故障数据与正常数据进行区分。但是该类方法无法准确对故障类型进行划分，而且故障识别精度很难进一步提升。
5.(2)基于关联规则的算法，如aprior算法，fp-growth算法等，该算法将智能电能表的事件记录以及相关运行数据的关联性进行分析，量化表示故障信息与诊断结果的关联关系，并按照置信度对模型进行评价。该类方法的优点在于没有复杂的理论推倒，易于实现；但是在运算过程中会产生大量的候选项集，对硬件要求较高，而且所用评价标准没有考虑到各个属性的重要程度，适应面较窄。
6.(3)基于人工智能的神经网络算法，如bp神经网络等，该算法将智能电能表的故障数据进行训练，获得良好的分类模型，对实时数据进行准确的故障预测。
7.但是现有故障预测方法存在一下不足之处：
8.1.定期巡检的方法需要耗费大量的人力物力，而且周期较长，并不能及时发现现场的故障情况。
9.2.现有基于人工智能的方法需要对数据进行人为的特征选择，与冗余数据的清除，否则会对分类模型的性能造成严重的影响，并且分类模型的鲁棒性和泛化能力较差。
10.3.现有基于人工智能的方法忽略了数据的时间特性与空间特性，分类模型并没有良好的捕捉到特征变量间的相互影响关系，和特征变量本身的时间特性，这会降低模型的分类能力。


技术实现要素：

11.为了改善上述技术方案的问题，本发明提供一种基于时空卷积神经网络的智能电表故障预测方法。
12.本发明采用如下的技术方案：一种基于时空卷积神经网络的智能电表故障预测方法，包括以下步骤：
13.(1)从智能电能表历史数据故障库中获得智能电能表的故障类型和故障相关数据；
14.(2)对智能电能表的历史数据进行清洗，包括异常值剔除和归一化处理；
15.(3)将相关变量的时间序列组成矩阵形式，使输入数据既包含时间信息又包含空间信息，并且基于cnn的强大特征捕捉能力，训练智能电能表的故障识别模型；
16.(4)采用准确率(accuracy)，精确率(precision)，召回率(recall)和综合评价指标f1作为智能电能表故障识别模型的评价标准；
17.(5)将实时数据通过数据清洗、构建成输入矩阵后，代入已训练完成的时空卷积神经网络智能电能表故障预测模型，对故障进行预测；
18.(6)当预测出现故障时，通过现场运维人员对故障的准确性进行判断，并且进行解决；同时通过运维人员确定的故障标签，作为新的样本数据进入历史故障库中进行重新训练，对智能电能表的故障预测模型参数进行实时更新，保证模型精度的可靠性和实时性。
19.进一步的，所述步骤1中：所述智能电能表故障相关数据包括总电量数，峰段电量数，平段电量数，谷段电量数，电流，电压，电表时间，终端时间；
20.所述智能电能表故障类别包括：正常、采集异常、电能表示值不平，电能表飞走，电能表反向有功示值异常，电能表倒走，电能表时间异常，电能表潜动异常。
21.进一步的，所述步骤2中：采集得到的数据难免会存在异常值，异常值会导致模型出现偏差，精度降低。采用hempel准则对异常值进行去除，hempel准则应用对异常值更加敏感的中值和中值绝对偏差代替常用的均值和均值偏差。根据该方法，对于每个特征变量x＝[x1,x2,
…
,xn]，满足下式的值即可被定义为异常值，
[0022][0023]
式中，xi表示特征变量中的第i个变量，x
0.5
表示变量中值。由于异常数据相对较少，本发明采用直接删除的方法对异常值进行处理，然后对每个特征变量依据下式进行归一化处理；
[0024][0025]
进一步的，所述步骤3中，还包括一下步骤：
[0026]
2)基于时间和空间特性构建输入矩阵
[0027]
定义融入时间和空间特性的输入矩阵为x，如下所示：
[0028][0029]
式中，s表示特征变量的个数，t表示选取的时间长度。
[0030]
2)将输入矩阵与cnn相结合并进行训练
[0031]
采用了2个卷积层，2个池化层，1个全连接层的模型结构进行建模，具体训练过程分为正向转播，反向转播。
[0032]
进一步的，，所述正向传播：
[0033]
11)卷积层c1[0034]
cnn网络结构的输入为x，则卷积层c1的传播过程如下式所示：
[0035][0036][0037]
式中，x和分别表示输入特征图(feature map)与输出特征图，q1表示输出的特征映射个数；表示卷积运算；表示卷积核；表示偏置；表示relu(rectified linear unit)激活函数使用非线性激活函数能够更好地捕捉多个输入变量之间的非线性关系。
[0038]
12)池化层p1[0039]
对于池化层p1的每一个特征映射组将其划分很多区域则其最大池化(取一个区域内所有神经元的最大值)可以表示为：
[0040][0041]
式中，为区域内每个神经元的激活值。
[0042]
13)卷积层c2[0043]
卷积层c2的传播过程如下式所示：
[0044][0045][0046]
式中，表示relu激活函数。
[0047]
14)池化层p2[0048]
池化层p2的传播过程如下式所示：
[0049][0050]
15)全连接层fc
[0051]
全连接层的作用则是对提取的特征进行非线性组合以得到输出，特征图在全连接层中会失去空间拓扑结构，被平铺(flatten)为向量。如下所示：
[0052][0053]
则全连接层输入为x
fc-1
，全连接层的正向传播过程可以表示为：
[0054]zfc
＝w
fc
·
x
fc-1
b
fc
[0055][0056]
式中，w
fc
表示全连接层权重，b
fc
表示全连接层偏置，f
fc
(
·
)表示sigmoid激活函
数，f
fc
(x)＝1/(1 exp(-x))。
[0057]
16)损失函数
[0058]
损失函数表示如下式所示：
[0059][0060]
式中，n表示训练样本个数，表示模型预测结果，y表示实际结果。
[0061]
进一步的，所述反向传播：
[0062]
cnn的反向传播就是对参数w
fc
，b
fc
，的更新。
[0063]
21)当前层为全连接层fc
[0064]
误差项的传播如下所示：
[0065][0066]
式中，表示点积运算，f
′
fc
(x)＝f
fc
(x) (1-f
fc
(x))。
[0067]
全连接层参数w
fc
梯度计算：
[0068][0069]
全连接层参数b
fc
梯度计算：
[0070][0071]
误差项传递到池化层p2，可以表示为
[0072][0073]
22)当前层为池化层p2[0074]
误差项由池化层p2传递到卷积层c2，如下所示：
[0075][0076]
式中，up(
·
)表示上采样运算，且
[0077][0078]
则卷积层c2层的参数梯度计算表示为：
[0079][0080][0081]
23)当前层为卷积层c2[0082]
误差项由卷积层c2传递到池化层p1，具体表示为：
[0083][0084]
式中，rot180(
·
)表示旋转180度。
[0085]
24)当前层为池化层p1[0086]
误差项由池化层p1传递到卷积层c1，具体表示为：
[0087][0088]
式中，
[0089][0090]
则卷积层c1层的参数梯度计算表示为：
[0091][0092][0093]
25)参数更新
[0094]
设置学习率α＝0.01，根据下式对模型参数w
fc
，b
fc
，，进行更新。
[0095][0096][0097][0098][0099][0100][0101]
式中，n
batch
表示批量样本个数。
[0102]
当训练误差满足时，停止模型训练，并保存模型参数。
[0103]
进一步的，所述步骤4中，所述当有1个模型评价指标小于95％时，则通过调整cnn超参数重新进行训练，直至达到模型评价要求。
[0104]
分类评价指标的定义：
[0105]
tp(true positive)：将正类预测为正类数
[0106]
tn(true negative)：将负类预测为负类数
[0107]
fp(false positive)：将负类预测为正类数
[0108]
fn(false negative)：将正类预测为负类数
[0109]
准确率(accuracy)的计算公式如下所示：
[0110][0111]
精确率(precision)的计算公式如下所示：
[0112][0113]
召回率(recall)的计算公式如下所示：
[0114][0115]
综合评价指标(f1)的计算公式如下所示：
[0116][0117]
在本发明的第二方面，提供了一种电子设备，包括：一个或多个处理器；存储装置，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-7中任一所述的方法。
[0118]
在本发明的第三方面，提供了一种计算机可读介质，其上存储有计算机程序，其中，所述计算机程序被处理器执行时实现如权利要求1-7中任一所述的方法。
[0119]
与现有技术相比，本发明的有益效果为：
[0120]
本发明解决了现有智能电能表故障预测算法准确率低，需要人为参与度高的问题。发明了一种将时间特性与空间特性融入输入数据的方法，并与cnn相结合，大大提高了智能电能表的故障预测精度，并且无需人为参与，降低了时间成本，是一种新型低复杂度，高精度的智能电能表故障预测方法。
[0121]
以下将结合附图与具体的实施例对本发明进行详细的解释说明。
附图说明
[0122]
图1为本发明的一些实施例中的基于时空卷积神经网络的智能电表故障预测方法的流程示意图；
[0123]
图2为本发明的时空卷积神经网络训练智能电能表故障识别模型示意图；
[0124]
图3为本发明的智能电能表故障预测多模型对比示意图。
具体实施方式
[0125]
为了便于理解本发明，下面将参照相关附图对本发明进行更加全面的描述，附图中给出了本发明的若干实施例，但是本发明可以通过不同的形式来实现，并不限于文本所描述的实施例，相反的，提供这些实施例是为了使对本发明公开的内容更加透彻全面。
[0126]
实施例一：
[0127]
请着重参照附图1-2，一种基于时空卷积神经网络的智能电表故障预测方法，包括以下步骤：
[0128]
(1)从智能电能表历史数据故障库中获得智能电能表的故障类型和故障相关数据；
[0129]
(2)对智能电能表的历史数据进行清洗，包括异常值剔除和归一化处理；
[0130]
(3)将相关变量的时间序列组成矩阵形式，使输入数据既包含时间信息又包含空间信息，并且基于cnn的强大特征捕捉能力，训练智能电能表的故障识别模型；
[0131]
(4)采用准确率(accuracy)，精确率(precision)，召回率(recall)和综合评价指标f1作为智能电能表故障识别模型的评价标准；
[0132]
(5)将实时数据通过数据清洗、构建成输入矩阵后，代入已训练完成的时空卷积神经网络智能电能表故障预测模型，对故障进行预测；
[0133]
(6)当预测出现故障时，通过现场运维人员对故障的准确性进行判断，并且进行解决；同时通过运维人员确定的故障标签，作为新的样本数据进入历史故障库中进行重新训练，对智能电能表的故障预测模型参数进行实时更新，保证模型精度的可靠性和实时性。
[0134]
请着重参照附图1-2，所述步骤1中：所述智能电能表故障相关数据包括总电量数，峰段电量数，平段电量数，谷段电量数，电流，电压，电表时间，终端时间；
[0135]
所述智能电能表故障类别包括：正常、采集异常、电能表示值不平，电能表飞走，电能表反向有功示值异常，电能表倒走，电能表时间异常，电能表潜动异常。
[0136]
请着重参照附图1-2，所述步骤2中：采集得到的数据难免会存在异常值，异常值会导致模型出现偏差，精度降低。采用hempel准则对异常值进行去除，hempel准则应用对异常值更加敏感的中值和中值绝对偏差代替常用的均值和均值偏差。根据该方法，对于每个特征变量x＝[x1,x2,
…
,xn]，满足下式的值即可被定义为异常值，
[0137][0138]
式中，xi表示特征变量中的第i个变量，x
0.5
表示变量中值。由于异常数据相对较少，本发明采用直接删除的方法对异常值进行处理，然后对每个特征变量依据下式进行归一化处理；
[0139][0140]
请着重参照附图1-2，其特征在于，所述步骤3中，还包括一下步骤：
[0141]
3)基于时间和空间特性构建输入矩阵
[0142]
定义融入时间和空间特性的输入矩阵为x，如下所示：
[0143][0144]
式中，s表示特征变量的个数，t表示选取的时间长度。
[0145]
2)将输入矩阵与cnn相结合并进行训练
[0146]
采用了2个卷积层，2个池化层，1个全连接层的模型结构进行建模，具体训练过程
分为正向转播，反向转播。
[0147]
请着重参照附图1-2，所述正向传播：
[0148]
11)卷积层c1[0149]
cnn网络结构的输入为x，则卷积层c1的传播过程如下式所示：
[0150][0151][0152]
式中，x和分别表示输入特征图(feature map)与输出特征图，q1表示输出的特征映射个数；表示卷积运算；表示卷积核；表示偏置；表示relu(rectified linear unit)激活函数使用非线性激活函数能够更好地捕捉多个输入变量之间的非线性关系。
[0153]
12)池化层p1[0154]
对于池化层p1的每一个特征映射组将其划分很多区域则其最大池化(取一个区域内所有神经元的最大值)可以表示为：
[0155][0156]
式中，为区域内每个神经元的激活值。
[0157]
13)卷积层c2[0158]
卷积层c2的传播过程如下式所示：
[0159][0160][0161]
式中，表示relu激活函数。
[0162]
14)池化层p2[0163]
池化层p2的传播过程如下式所示：
[0164][0165]
15)全连接层fc
[0166]
全连接层的作用则是对提取的特征进行非线性组合以得到输出，特征图在全连接层中会失去空间拓扑结构，被平铺(flatten)为向量。如下所示：
[0167][0168]
则全连接层输入为x
fc-1
，全连接层的正向传播过程可以表示为：
[0169]zfc
＝w
fc
·
x
fc-1
b
fc
[0189]
误差项由池化层p2传递到卷积层c2，如下所示：
[0190][0191]
式中，up(
·
)表示上采样运算，且
[0192][0193]
则卷积层c2层的参数梯度计算表示为：
[0194][0195][0196]
23)当前层为卷积层c2[0197]
误差项由卷积层c2传递到池化层p1，具体表示为：
[0198][0199]
式中，rot180(
·
)表示旋转180度。
[0200]
24)当前层为池化层p1[0201]
误差项由池化层p1传递到卷积层c1，具体表示为：
[0202]
[0203]
式中，
[0204][0205]
则卷积层c1层的参数梯度计算表示为：
[0206][0207][0208]
25)参数更新
[0209]
设置学习率α＝0.01，根据下式对模型参数w
fc
，b
fc
，，进行更新。
[0210][0211][0212][0213][0214][0215][0216]
式中，n
batch
表示批量样本个数。
[0217]
当训练误差满足时，停止模型训练，并保存模型参数。
[0218]
请着重参照附图1-2，所述步骤4中，所述当有1个模型评价指标小于95％时，则通过调整cnn超参数重新进行训练，直至达到模型评价要求。
[0219]
分类评价指标的定义：
[0220]
tp(true positive)：将正类预测为正类数
[0221]
tn(true negative)：将负类预测为负类数
[0222]
fp(false positive)：将负类预测为正类数
[0223]
fn(false negative)：将正类预测为负类数
[0224]
准确率(accuracy)的计算公式如下所示：
[0225][0226]
精确率(precision)的计算公式如下所示：
[0227][0228]
召回率(recall)的计算公式如下所示：
[0229][0230]
综合评价指标(f1)的计算公式如下所示：
[0231][0232]
本发明的一种基于时空卷积神经网络的智能电表故障预测方法的实施发式：
[0233]
1.样本数据选取
[0234]
选择智能电能表故障库中的故障数据，故障类别共8类，分别是：正常、采集异常、电能表示值不平，电能表飞走，电能表反向有功示值异常，电能表倒走，电能表时间异常和电能表潜动异常。
[0235]
智能电能表故障相关数据共8种：包括总电量数，峰段电量数，平段电量数，谷段电量数，电流，电压，电表时间，终端时间。
[0236]
2.输入矩阵构造
[0237]
数据的采样频率是秒级，分别选取每个故障数据连续运行30分钟的数据，采用1分钟的窗口对数据进行步长为1的滑动处理。
[0238]
每类故障可以得到1800个样本，输入为60
×
8的矩阵。
[0239]
将前20分钟连续的数据即每类故障1200个样本作为训练，后10分钟连续的数据即每类故障600个样本作为预测。
[0240]
3.时空卷积神经网络模型结构
[0241]
本发明所采用的时空卷积神经网络模型结构见下表所示，共采用2个卷积层，2个池化层，1个全连接层。
[0242]
本发明采用的核函数是3
×
3，这是因为3
×
3是最小的能够捕获像素八邻域信息的尺寸，而且多个3
×
3的卷积层比一个大尺寸卷积层使用了更多的非线性激活函数，判决性更强，同时使用的参数更少，模型复杂度更低。
[0243]
本发明采用3
×
3核函数堆叠的方法对模型结构进行试凑，当模型具有2层卷积层的时候，模型精度已经可以满足故障预测的需求，过高的层数在模型精度提升并不大，同时会造成大量的计算冗余，如表1所示。
[0244]
表1：时空卷积神经网络模型结构
[0245][0246]
4.智能电能表故障预测结果
[0247]
本发明将目前智能电能表故障预测常用的bp神经网络，随机森林，svm，朴素贝叶斯与本发明方法作对比，结果如附图3所示。
[0248]
从表2可以看出，本发明方法与其他4个模型相比，无论是在准确率，精确率，召回率还是f1，都是表现最好的模型，这也说明了本文算法的有效性，其次较好的模型是朴素贝叶斯，在准确率，召回率和f1比其他3个模型较高，性能较差的模型是bp神经网络。
[0249]
表2智能电能表故障预测多模型对比
[0250][0251]
表3给出了基于时空卷积神经网络的方法在每种故障的预测结果：
[0252]
表3本发明方法各故障的预测结果
[0253][0254]
[0255]
从上表3中可知，采集异常与电能表示值不平这2个故障预测准确率较低，其他故障的预测准确率相差不多，其中，电能表反向有功示值异常预测准确率最好。后续的研究方向，可以针对预测准确率低的故障开展详细研究，找到准确率低的具体原因，通过输入特征的修改和模型结构超参数的调整提高模型的整体性能。
[0256]
上述结合附图和实施例对本发明进行了示例性描述，显然本发明具体实现并不受上述方式的限制，只要采用了本发明的方法构思和技术方案进行的这种非实质改进，或未经改进将本发明的构思和技术方案直接应用于其他场合的，均在本发明的保护范围之内。