一种基于图神经网络的遥测时序数据异常检测方法及系统与流程

1.本发明涉及遥测时序数据异常检测技术领域，特别是涉及一种基于图神经网络的遥测时序数据异常检测方法及系统。


背景技术：

2.卫星是一个庞大、复杂的高精度控制系统，并且众多分系统和子系统存在耦合与交错。近年来，随着航天技术的发展，对卫星运行的稳定性和可靠性要求越来越高。卫星在轨运行的空间环境十分恶劣，一旦卫星发生在轨故障，若不及时处理则可能导致系统服务、可用性变差，甚至无法使用，带来巨大的经济损失。及时有效的异常检测和故障定位能保障卫星安全可靠的运行，因此卫星故障诊断系统在卫星管理过程中发挥着愈来愈重要的作用。其中，卫星遥测时序数据流异常检测可以对卫星的故障进行早期预警，是其中最为关键的一部分，因此，针对卫星遥测时序数据异常检测具有重要意义。
3.国内外已经有了一些关于遥测时序数据异常检测的研究成果，根据异常检测的算法原理，异常检测算法大体可以分为五类：基于聚类的方法、基于分类的方法、基于邻域的方法、基于统计的方法和基于回归的方法。由于卫星数据具有数据量大、有噪声、变化缓慢和维度高等特点，而深度学习对于处理高维数据效果很好，深度学习方法开始大量用于异常检测领域。
4.但是这些深度学习方法仍然存在一些缺陷，目前异常检测领域中的深度学习算法模型大多数都忽视了变量之间和时间序列之间潜在的相互依赖性，即每时刻的卫星某一参数不仅仅取决于它的历史数据，还与其它参数存在依赖性，不考虑它们的依赖性对预测模型精度会造成严重的影响，例如向量自回归模型和高斯过程模型是通过假设变量存在线性相关性，但是当针对非线性模式数据效果就会变差。


技术实现要素：

5.本发明的目的是提供一种基于图神经网络的遥测时序数据异常检测方法及系统，以提高遥测时序数据异常检测的准确性。
6.为实现上述目的，本发明提供了如下方案：
7.一种基于图神经网络的遥测时序数据异常检测方法，包括：
8.对历史遥测时序数据进行预处理，得到预处理后的历史数据；
9.通过所述历史数据训练图神经网络遥测时序数据预测模型，得到训练好的预测模型；
10.利用小波方差计算待测遥测时序数据的周期；所述待测遥测时序数据包括多个时刻的待测遥测数据；
11.利用所述训练好的预测模型对所述待测遥测时序数据进行预测，得到预测遥测时序数据；所述预测遥测时序数据包括多个时刻的所述待测遥测数据的预测数据；
12.根据所述待测遥测时序数据的周期以及多个时刻的所述待测遥测数据的预测数
据确定预测周期数据；
13.计算当前时刻的待测遥测数据和所述当前时刻的待测遥测数据的预测数据之间的马氏距离，记为第一马氏距离；
14.计算之前时刻的待测遥测数据和所述之前时刻的待测遥测数据的预测周期数据之间的马氏距离，记为第二马氏距离；所述之前时刻的待测遥测数据包括上一时刻的待测遥测数据；
15.计算所述第二马氏距离的平均值和方差；
16.将所述平均值和所述方差进行放大设置阈值；
17.根据所述第一马氏距离和所述阈值判断所述当前时刻的待测遥测数据是否异常。
18.可选地，所述对历史遥测时序数据进行预处理，得到预处理后的历史数据，具体包括：
19.对所述历史遥测时序数据进行清洗、归一化处理、小波去噪以及灰色关联度分析，得到预处理后的历史数据。
20.可选地，所述图神经网络遥测时序数据预测模型包括图形学习层、图形卷积层以及时间卷积层。
21.可选地，所述待测遥测时序数据的的小波方差计算公式如下：
[0022][0023]
其中，var(a)为待测遥测时序数据的小波方差，w
f
(a,b)为小波变换系数，a为伸缩尺度因子，b为平移因子，d为积分。
[0024]
可选地，所述第一马氏距离的计算公式如下：
[0025][0026]
其中，re
score
(i)为第i个时刻的第一马氏距离，x(i)为第i个时刻的数据记录，为经过训练好的预测模型的第i个时刻的数据记录，为数据集的协方差矩阵的逆矩阵，x为待测遥测数据集。
[0027]
可选地，所述阈值的计算公式如下：
[0028]
threshold(t)＝k1μ
y
+k2σ
y
[0029]
其中，threshold(t)为t时刻的阈值，μ
y
为平均值，σ
y
为方差，y为第二马氏距离构成的数据集，k1为平均值的放大倍数，k2为方差的放大倍数。
[0030]
一种基于图神经网络的遥测时序数据异常检测系统，包括：
[0031]
预处理模块，用于对历史遥测时序数据进行预处理，得到预处理后的历史数据；
[0032]
训练模块，用于通过所述历史数据训练图神经网络遥测时序数据预测模型，得到训练好的预测模型；
[0033]
第一计算模块，用于利用小波方差计算待测遥测时序数据的周期；所述待测遥测时序数据包括多个时刻的待测遥测数据；
[0034]
预测模块，用于利用所述训练好的预测模型对所述待测遥测时序数据进行预测，得到预测遥测时序数据；所述预测遥测时序数据包括多个时刻的所述待测遥测数据的预测数据；
[0035]
预测周期数据确定模块，用于根据所述待测遥测时序数据的周期以及多个时刻的所述待测遥测数据的预测数据确定预测周期数据；
[0036]
第二计算模块，用于计算当前时刻的待测遥测数据和所述当前时刻的待测遥测数据的预测数据之间的马氏距离，记为第一马氏距离；
[0037]
第三计算模块，用于计算之前时刻的待测遥测数据和所述之前时刻的待测遥测数据的预测周期数据之间的马氏距离，记为第二马氏距离；所述之前时刻的待测遥测数据包括上一时刻的待测遥测数据；
[0038]
第四计算模块，用于计算所述第二马氏距离的平均值和方差；
[0039]
阈值设置模块，用于将所述平均值和所述方差进行放大设置阈值；
[0040]
判断模块，用于根据所述第一马氏距离和所述阈值判断所述当前时刻的待测遥测数据是否异常。
[0041]
可选地，所述图神经网络遥测时序数据预测模型包括图形学习层、图形卷积层以及时间卷积层。
[0042]
可选地，所述待测遥测时序数据的的小波方差计算公式如下：
[0043][0044]
其中，var(a)为待测遥测时序数据的小波方差，w
f
(a,b)为小波变换系数，a为伸缩尺度因子，b为平移因子，d为积分。
[0045]
可选地，所述第一马氏距离的计算公式如下：
[0046][0047]
其中，re
score
(i)为第i个时刻的第一马氏距离，x(i)为第i个时刻的数据记录，为经过训练好的预测模型的第i个时刻的数据记录，为数据集的协方差矩阵的逆矩阵，x为待测遥测数据集。
[0048]
根据本发明提供的具体实施例，本发明公开了以下技术效果：
[0049]
本发明通过图神经网络遥测时序数据预测模型学习多维时序数据之间的依赖关系，并通过周期窗口的动态阈值法对待测遥测时序数据进行异常检测，提高了遥测时序数据异常检测的准确性。
附图说明
[0050]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0051]
图1为本发明实施例提供的基于图神经网络的遥测时序数据异常检测方法的流程图；
[0052]
图2为本发明实施例提供的基于图神经网络的遥测时序数据异常检测方法原理图一；
[0053]
图3为本发明实施例提供的基于图神经网络的遥测时序数据异常检测方法原理图
二；
[0054]
图4为本发明实施例提供的小波去噪流程示意图；
[0055]
图5为本发明实施例提供的图神经网络遥测时序数据预测模型训练过程图。
具体实施方式
[0056]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0057]
本发明的目的是提供一种基于图神经网络的遥测时序数据异常检测方法及系统，以提高遥测时序数据异常检测的准确性。
[0058]
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0059]
如图1
‑
3所示，基于图神经网络的遥测时序数据异常检测方法包括：
[0060]
步骤101：对历史遥测时序数据进行预处理，得到预处理后的历史数据。具体包括：
[0061]
对所述历史遥测时序数据进行清洗、归一化处理、小波去噪以及灰色关联度分析，得到预处理后的历史数据。其中，小波去噪采用全局统一阈值σ＝mad/0.6745，mad是首层小波分解系数绝对值的中间值，0.6745是高斯噪声标准方差的调整系数，n是信号的尺度。图4为本发明实施例提供的小波去噪流程示意图。
[0062]
步骤102：通过所述历史数据训练图神经网络遥测时序数据预测模型，得到训练好的预测模型。其中，所述图神经网络遥测时序数据预测模型包括图形学习层、图形卷积层以及时间卷积层。如图5所示，训练过程如下：
[0063]
(1)将历史数据输入图神经网络遥测时序数据预测模型，图形学习层采用抽样的方法计算每批次节点之间的关系，学习邻接矩阵，节点之间关系计算公式如下所示：
[0064]
n1＝tanh(βe1θ1)
[0065]
n2＝tanh(βe2θ2)
[0066][0067]
其中，e1和e2表示节点，θ1和θ2表示学习过程中的模型参数，relu表示线性整流函数，n1和n2表示过度变量，b表示学习得到的邻接矩阵，β表示超参数。从计算公式中可以看出节点之间的关系是单向的。
[0068]
(2)图形卷积层提取图中节点之间的空间特征，将节点信息和其邻居信息相融合。
[0069]
(3)时间卷积层采用1d卷积滤波器提取时间特征，获取遥测时序数据的序列模式。
[0070]
(4)保存训练过程误差最小的模型。
[0071]
步骤103：利用小波方差计算待测遥测时序数据的周期；所述待测遥测时序数据包括多个时刻的待测遥测数据。小波方差图的最高点为待测遥测时序数据的周期。小波方差的计算公式如下：
[0072][0073]
其中，var(a)为待测遥测时序数据的小波方差，a为伸缩尺度因子，b为平移因子，d
为积分，w
f
(a,b)为小波变换系数，(a,b)为小波变换系数，为ψ(t)的共轭函数，ψ(t)为基小波函数，r为实数集，t为时刻，f(t)为信号函数。
[0074]
步骤104：利用所述训练好的预测模型对所述待测遥测时序数据进行预测，得到预测遥测时序数据；所述预测遥测时序数据包括多个时刻的所述待测遥测数据的预测数据。
[0075]
步骤105：根据所述待测遥测时序数据的周期以及多个时刻的所述待测遥测数据的预测数据确定预测周期数据。
[0076]
步骤106：计算当前时刻的待测遥测数据和所述当前时刻的待测遥测数据的预测数据之间的马氏距离，记为第一马氏距离。其中，所述第一马氏距离的计算公式如下：
[0077][0078]
其中，re
score
(i)为第i个时刻的第一马氏距离，x(i)为第i个时刻的数据记录，为经过训练好的预测模型的第i个时刻的数据记录，为数据集的协方差矩阵的逆矩阵，x为待测遥测数据集。
[0079]
步骤107：计算之前时刻的待测遥测数据和所述之前时刻的待测遥测数据的预测周期数据之间的马氏距离，记为第二马氏距离；所述之前时刻的待测遥测数据包括上一时刻的待测遥测数据。
[0080]
步骤108：计算所述第二马氏距离的平均值和方差；
[0081]
步骤109：将所述平均值和所述方差进行放大设置阈值；所述阈值的计算公式如下：
[0082]
threshold(t)＝k1μ
y
+k2σ
y
[0083]
其中，threshold(t)为t时刻的阈值，μ
y
为平均值，σ
y
为方差，y为第二马氏距离构成的数据集，k1为平均值的放大倍数，k2为方差的放大倍数。
[0084]
步骤110：根据所述第一马氏距离和所述阈值判断所述当前时刻的待测遥测数据是否异常。
[0085]
本实施例还提供了一种基于图神经网络的遥测时序数据异常检测系统，系统包括：
[0086]
预处理模块，用于对历史遥测时序数据进行预处理，得到预处理后的历史数据。
[0087]
训练模块，用于通过所述历史数据训练图神经网络遥测时序数据预测模型，得到训练好的预测模型。
[0088]
第一计算模块，用于利用小波方差计算待测遥测时序数据的周期；所述待测遥测时序数据包括多个时刻的待测遥测数据。
[0089]
预测模块，用于利用所述训练好的预测模型对所述待测遥测时序数据进行预测，得到预测遥测时序数据；所述预测遥测时序数据包括多个时刻的所述待测遥测数据的预测数据。
[0090]
预测周期数据确定模块，用于根据所述待测遥测时序数据的周期以及多个时刻的所述待测遥测数据的预测数据确定预测周期数据。
[0091]
第二计算模块，用于计算当前时刻的待测遥测数据和所述当前时刻的待测遥测数据的预测数据之间的马氏距离，记为第一马氏距离。
[0092]
第三计算模块，用于计算之前时刻的待测遥测数据和所述之前时刻的待测遥测数据的预测周期数据之间的马氏距离，记为第二马氏距离；所述之前时刻的待测遥测数据包括上一时刻的待测遥测数据。
[0093]
第四计算模块，用于计算所述第二马氏距离的平均值和方差。
[0094]
阈值设置模块，用于将所述平均值和所述方差进行放大设置阈值。
[0095]
判断模块，用于根据所述第一马氏距离和所述阈值判断所述当前时刻的待测遥测数据是否异常。
[0096]
进一步地，所述图神经网络遥测时序数据预测模型包括图形学习层、图形卷积层以及时间卷积层。
[0097]
进一步地，所述待测遥测时序数据的的小波方差计算公式如下：
[0098][0099]
其中，var(a)为待测遥测时序数据的小波方差，w
f
(a,b)为小波变换系数，a为伸缩尺度因子，b为平移因子，d为积分。
[0100]
进一步地，所述第一马氏距离的计算公式如下：
[0101][0102]
其中，re
score
(i)为第i个时刻的第一马氏距离，x(i)为第i个时刻的数据记录，为经过训练好的预测模型的第i个时刻的数据记录，为数据集的协方差矩阵的逆矩阵，x为待测遥测数据集。
[0103]
本实施例还可以通过以下代码实现：
[0104]
输入:
[0105]
训练集o；
[0106]
测试集g；
[0107]
节点集v；
[0108]
批次大小b；
[0109]
步长s；
[0110]
学习率γ；
[0111]
输出:
[0112]
测试集异常检测集合se；
[0113]
[0114]
[0115][0116]
本发明提供的基于图神经网络的无领域知识、动态阈值遥测时序数据异常检测方法及系统，是无监督的遥测时序数据异常检测方法及系统。传统的图神经网络需要明确定义图的结构，不能学习多维时间序列之间的依赖关系。本发明中的图形学习层，自适应将多维遥测时序数据学习为图邻接矩阵，用一种新颖的方法学习多维遥测时间序列数据之间的依赖关系，即变量之间和时间序列之间的潜在依赖性，将多维遥测时序数据转化为图。采用马氏距离计算样本异常得分，遥测时序数据采集过程中不确定干扰比较多，异常得分往往峰值比较多，并且卫星围绕地球旋转，卫星数据往往具有周期性，因此本文采用基于周期窗口的动态阈值法对卫星时序数据流异常检测，避免了遥测时序数据采集过程中的不确定干扰影响。
[0117]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0118]
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。