基于GRU-VAE的无监督航迹异常检测方法.pdf

第 卷 第 期 年 月指挥控制与仿真 文章编号：（）基于 的无监督航迹异常检测方法李 磊，张 静，欧阳齐铖，周明康（中国人民解放军战略支援部队信息工程大学，河南 郑州；中国通信建设第四工程局有限公司，河南 郑州）摘 要：针对海量无行为模式标签航迹数据的目标行为异常检测问题，提出一种基于门控循环单元的变分自编码器模型（，）无监督航迹异常检测方法。该方法通过检测航迹异常发现目标的行为异常，分为模型训练阶段和异常检测阶段两步实施：在模型训练阶段，构建了以 和 为主要结构元素的基于门控循环单元的变分自编码器模型，利用无异常信息标签的历史航迹数据对 模型进行训练，根据训练集的航迹点重构损失分布情况，采用正态分布法或百分位数法划定置信区间为航迹点重构损失门限；在异常检测阶段，该模型对实时航迹数据集进行检测，将重构损失超出航迹点重构损失门限的航迹点视为异常航迹点，当航迹序列中的异常航迹点占比超出占比阈值时，判定为异常航迹序列，结合数据异常情况向一线人员发送目标的异常行为信息。数据实验结果表明，模型最高 分数达 ，查全率达。本方法对异常航迹的检测具有高灵敏度和低漏警率，可满足战场态势认知需求。关键词：数据挖掘；航迹数据；异常检测；无监督学习中图分类号：文献标志码：，（，；，）：，：，（），：；收稿日期：修回日期：作者简介：李磊（），男，硕士研究生，助理工程师，研究方向为智能信息处理。张 静（），女，博士，副教授。当前，外军舰机对我周边军事侦察活动日渐频繁，国家海上安全形势不容乐观，对海上高价值目标的技术侦察指挥控制是决定战争胜负的关键。如图 所示，在基于 环的战场侦察指挥控制链中，战场态势认知至关重要。战场态势认知是支撑指挥决策的重要依据。只有基于对战场当前“态”的全面理解，才能实现对战场未图 侦察指挥控制决策链 李 磊，等：基于 的无监督航迹异常检测方法第 卷来“势”的准确预测与推演，从而提升战场态势理解效率。航迹异常检测技术是海上态势情报获取的关键技术之一。充分挖掘多种渠道获取的航迹数据，掌握目标时空特征和运动模式，通过检测异常航迹，提前发现目标异常行为，可为高效侦察指挥决策提供有力支撑。通过实地调研发现，当前，在基于航迹特征的舰船异常行为检测业务上存在以下两点困难：）海量航迹数据大多无法准确标注其行为模式信息，导致正常航迹数据和异常航迹数据样本量极不均衡，监督学习法在此数据条件下极易出现过拟合现象，模型性能受到严重制约；）战场态势瞬息万变，伴随各类舰船目标行为模式的不断变化，新的异常航迹类型不断出现，过去标注的异常航迹难以表征新出现的异常航迹类型。因此，针对海量无行为模式标签航迹数据的目标行为异常检测问题，本文方法的思路是对历史航迹数据集进行深入挖掘，建立舰船目标正常运动规律模型，通过比较实时航迹数据和模型输出的偏差发现异常航迹。本文把利用历史航迹数据构建常态模型的无监督异常检测方法作为主要研究方向。在该研究方向上，基于无监督网络模型的异常检测方法值得关注。文献以自编码器（，）作为水位异常检测框架核心，通过学习正常数据的特征分布，将输出的重构损失作为异常分数并设定阈值，进而实现对水位数据的异常检测。变分自编码器（，）以时间序列中潜在变量分布参数的重建概率作为异常检测的度量，不受数据结构限制，比自编码器在数据重构上更有优势。文献采用变分自编码器模型对脑电数据进行异常检测，降低癫痫发作带来的安全风险。然而，以上两种模型针对航迹数据的时序建模能力还有待提高。长短时记忆网络（，）和门控循环单元（，）作为循环神经网络的重要变体，在时间序列建模问题上具有优势。在其他领域中，以上两种自编码器模型常配合使用以增强其时序建模能力。文献提出一种基于长短时记忆网络模型 的无监督异常检测模型，在此基础上，文献将 和 相结合，进一步提高了模型在时序数据异常检测上的表现。文献将 作为 编码器的输入层和解码器的输出层，进一步 提 升 了 模 型 对 时 序 数 据 的 重 构 能 力。而的参数量比 更少，保持了 的优异性能，同时其结构更加简单，过拟合风险更低，文献将其应用于飞机振动数据的异常检测，并取得了良好效果。本文构建了以 和 为主要结构元素的无监督神经网络模型，提出了一种基于 模型的无监督航迹异常检测方法。在模型训练阶段，使用历史航迹数据集对 模型进行训练，对数据集中所有航迹序列进行重构，将重构航迹序列和原始航迹序列的平均绝对距离作为模型输出的重构损失，基于输出重构损失的分布类型确定其在任意置信度下的置信区间，并将其作为航迹点重构损失门限；在异常检测阶段，该模型对实时航迹数据进行检测，将损失超出航迹点重构损失门限的航迹点视为异常航迹点，将异常航迹点占比超出阈值的航迹序列判定为异常航迹序列，再结合数据特征和重构损失异常情况向一线人员发送目标的异常行为信息。数据实验结果表明，模型平均 分数高达 ，查全率达。对比现有相关研究成果，本方法在查准率、查全率和模型 分数上比实验性能稍次的基于长短时记忆网络的变分自编码器模型（，）方法分别高出 、和 ，对异常航迹具有高灵敏度和低漏警率，符合战场态势认知需求。问题描述针对海量无异常标签数据的航迹异常检测问题，本节在相关定义的基础上，提出对应研究思路，并对异常航迹序列进行分类。基本定义定义：航迹点向量及特征，分别对应的是航迹点在时间戳 的经度（）、纬度（）、速度（）和航向（）特征。其中，是指按照时间顺序排列的航迹点序数，经度和纬度坐标是航迹点向量的位置属性特征，速度和航向是航迹点向量的运动属性特征。定义：航迹序列 （，）是指由 个航迹点向量根据时间戳排列得到的序列。由定义，航迹点序列 可以表示为式（）：，（）定义：航迹数据集 ，表示由 条航迹序列组成的航迹数据集合。第 期指挥控制与仿真 定义：异常航迹点向量若实时航迹序列 经模型检测输出的重构损失序列 中存在重构损失值 ，则重构损失值 对应的航迹点向量 为异常航迹点向量。定义：异常航迹序列设异常航迹点向量占所有航迹点向量总数之比的阈值为，若由 个航迹点向量组成的实时航迹序列中共包含 个异常航迹点向量，当 时，则该实时航迹序列 （，）为异常航迹序列。研究思路针对海量无异常标签数据的航迹异常检测问题，本节提出基于航迹数据重构模型的无监督航迹异常检测方法系统模型，如图 所示。图 研究思路 结合 节基本定义，设由给定时间和空间范围内的历史航迹序列组成的历史航迹数据集为，实时航迹数据集为，集合 中的 是与历史航迹数据集行为模式不同的异常航迹序列，所有集合中的航迹数据均无异常行为模式标签。输入历史航迹数据集 训练航迹数据重构模型后，能够输出对应的重构损失，并据此确定航迹点重构损失门限为，。将实时航迹数据集 中任意一条包含 个航迹点向量的实时航迹序列 （，）输入航迹数据重构模型，能够输出包含 个重构损失值的重构 损 失 序 列 为（，）。研究思路分两步：）航迹数据重构模型通过训练学习集合 的舰船历史航迹运动模式，根据历史航迹数据重构损失分布，确定航迹点重构损失门限，；）将集合 中的实时航迹序列逐条输入模型，经航迹异常检测方法发现异常航迹序列。异常航迹序列的分类异常航迹序列可按照航迹点向量特征的分类划分为两类：）运动异常，目标实际行为体现为舰船航速超出历史范围或短时间内变化过快，以及在侦察区域出现的短时间内航向变化过快；）位置异常，目标实际行为体现为偏离历史航迹或出现在非法位置。本方法可检测包含其中一类异常或同时存在这两类异常的航迹序列。如图 所示。图 异常航迹序列的分类 方法描述基于 模型的无监督航迹异常检测方法总体流程如图 所示，主要分两步：）模型训练阶段，使用预处理后得到的历史航迹数据集训练 航迹数据异常检测模型，采用平均绝对误差计算重构损失，基于输出重构损失的分布类型确定重构损失的置信区间，进而确定航迹点重构损失门限；）异常检测阶段，运用训练好的模型对预处理后的实时航迹数据集进行检测，对异常航迹序列进行深入研判，确定其异常类型。航迹数据集预处理为减少数据不同量纲对深度神经网络训练的影响，使模型有效拟合训练数据，准确学习正常航迹的运动模式和规律，需要对历史航迹数据集进行预处理。历史航迹数据集的预处理流程如图 所示，分为以下 步。：根据舰船类型标签，在航迹数据集中筛选出单个目标的完整航迹序列以及各时间戳下的航迹点 李 磊，等：基于 的无监督航迹异常检测方法第 卷图 方法总体思路 图 预处理流程 向量；：针对参与训练的航迹数据集，对同一目标类型的航迹序列在相同时空范围内的重复航迹予以剔除，避免深度神经网络模型训练的计算冗余或将其误识别为静止点；：对航迹序列中航迹点向量特征值缺失的特征进行均值插值，避免数据缺失问题影响模型对常态运动模式的准确学习，均值计算式如式（）所示：（）其中，为非缺失航迹点特征值 的总个数。：对所有的航迹点向量特征值进行 均值标准化处理，即将航迹数据集转换为均值为、方差为 的数据集，避免在模型训练和重构损失的计算过程中由于量纲不同、自身数值相差较大所引起的误差。先计算所有同一特征值的标准差，如式（）所示：（）（）其中，为该航迹序列中航迹点向量总个数，之后运用式（）计算标准化处理后的特征值：（）式中，指同一航迹数据集或航迹序列中航迹点向量的某特征原始值，为标准化处理后得到的结果。对需要进行检测且可能存在异常航迹的实时航迹数据集进行除 之外的预处理步骤。航迹异常检测模型本节对 和 的理论基础进行描述，并在此基础上搭建本文的 模型。）变分自编码器 变分自编码器是一种无监督学习的深度神经网络模型，它基于反向传播算法和梯度下降最优化方法，利用输入数据作为监督指导神经网络训练学习实现重构输出。如图 所示，变分自编码器由 个编码器和 个解码器组成，其基本思想是经过编码器函数（）的每一个点向量 可用各航迹点向量特征值上的均值 和方差 表示输出的概率分布，然后，通过变分自编码器的变分推断过程，用一个近似分布（）代替概率分布，并通过贝叶斯准则得到解码器输入（），最后，经解码器（）输出重构的航迹点向量。图 模型结构 在变分推断过程中，采用（，）散度函数来度量两个分布之间的相似性，模型的目标函数如式（）所示：（，；）（）（）（）（）第 期指挥控制与仿真 其中，是期望的自编码重建输出特征向量；是正则化项，使得（）向标准正态分布对标。为了解决神经网络反向梯度求解问题，中间采用了重参数化的操作，样本 从分布（，）中获取变量，通过图 所示的期望 和方差 计算得到，见式（）。（）联立式（）和式（），最终的目标函数如式（）所示。（，；）（）（）（）（）门控循环单元门控循环单元的主要功能是在数据重构过程中捕捉航迹数据的时序特性，提升数据重构的准确性。作为循环神经网络的一种变形单元，门控循环单元解决了长期记忆和反向传播中的梯度消失和爆炸问题，在达到相同模型性能的前提下，其计算量小，训练效率更高。图 为 运算示意图，它主要的功能由更新门和重置门实现，这些门结构决定了数据信息的输入、存储、遗忘和输出。图 运算示意图 图 中，表示在时刻状态，为时间步的数据输入，和 为非线性激活函数。式（）（）展示了 的运算过程。和 表示更新门和重置门结构的权重和偏置，其余运算符号如图中所示。式（）中更新门 控制从上一时刻状态 中保留的向量特征信息以及从候选状态中接收的向量特征信息，以确保数据的不断更新。（，）（）式（）中的重置门 决定了候选状态 对上一时刻 的依赖程度。（，）（）（，）（）（）（）由式（）和式（）可见，使用更新门控制输入和遗忘的平衡。当 时，当前状态 和上一状态为非线性关系；当 时，和 为线性关系。）模型的搭建 模型是在 框架下，将编码器输入层和解码器输出层神经元替换为 门控循环单元，提升了模型对航迹数据时序特征的学习能力，克服了传统循环神经网络在特征提取过程中的梯度消失和爆炸问题，实现了对数据特征的正确编码和解码。模型通过学习航迹数据的特定分布，输出得到精准的重构数据。模型的基本结构如图 所示，数据重构过程分为以下三步。首先，航迹序列 （，）直接输入模型，编码器的 单元捕捉航迹序列点向量的时序特征，经 线性激活层分别输出特征值对应的均值 和方差；然后，在 框架下进行重参数化操作，对点向量的隐变量分布 进行估计；最后，单元译码器根据隐变量 完成对数据的时序性重构，输出重构航迹序列，实现端到端的数据精确重构。图 模型基本结构图 模型训练阶段模型训练阶段流程如图 所示。首先，将历史航迹数据集 输入模型进行 训 练；然 后，模 型 输 出 训 练 集 重 构 损 失 集，设定重构损失的置信度为 ；之后，对进行正态性检验，结合设置的置信度，若满足正态分布，则采用正态分布法计算航迹点重构损失门限，否则，采用百分位数法计算航迹点重构损失门限；最后，保存模型及其对应的航迹点重构损失门限，。）重构损失计算方法航迹序列重构损失的计算采用平均绝对误差计算方法，如式（）所示。李 磊，等：基于 的无监督航迹异常检测方法第 卷图 方法工作流程 （）其中，和 是重构航迹序列和原始航迹序列，和 分别是第 个时间戳上的重构航迹点向量和原始航迹点向量，是航迹序列中的航迹点向量总个数，是航迹点向量中的特征个数。经计算得到重构损失序列。）航迹点重构损失门限的确定确定航迹点重构损失门限的目的是通过模型训练掌握历史航迹数据隐含的目标行为模式，为异常检测阶段提供常态基准参考。下面根据历史航迹数据输入模型得到的重构损失分布，采用正态分布法和百分位数法，对需要检测的实时航迹数据输出的重构损失值进行置信区间的估计，将置信区间作为航迹点重构损失门限。图 为航迹 点重构损失 门限示意 图。假设服从标准正态分布，将重构损失置信度（）下的置信区间作为航迹点重构损失门限，表示为参与检测的实时航迹序列中的每个航迹点向量对应重构损失值有（）的概率会落在区间，中。若训练集重构损失序列 中重构损失值的分布能够通过正态分布检验，则采用正态分布法进行置信区间的估计。图 航迹点重构损失门限示意图 设重构损失 是一个服从正态分布的随机变量，即 （，），其中，为总体均值，为总体标准差。模型是深度神经网络模型，需要大样本量数据进行训练，相应输出重构损失序列中的样本数量必然大于。故由大数定律和中心极限定理可知，数据样本量足够大时，样本均值、方差与总体的均值和方差相等，即样本均值 ，样本标准差 ，之后对该随机变量 作线性运算，将其转换为标准正态分布：（，）（）应用标准正态分布下的 分位点的定义可得 （）根据式（），可以求得航迹点重构损失在任意置信度 下的置信区间为 ，（）其中，为区间的显著性水平，界值可通过查标准正态分布表获得，该置信区为航迹点重构损失门限。若重构损失分布未能通过严格的正态性检验，则采用百 分 位 数 法 确 定。首 先，将 重 构 损 失 序 列中的所有重构损失值按升序排列；然后，合并相同的数值并计算重构损失值 在重构损失序列中出现次数占全部重构损失值的比例，进而得到每个重构损失值的出现概率；最后，在已知置信度为 的条件下，确定航迹点重构损失出现可能性最大的区域，对称分配 个百分位数 和，从而确定航迹点重构损失门限为，。异常检测阶段将对保存的模型及其对应的航迹点重构损失门限进行调用。异常检测阶段异常检测阶段主要是对实时航迹序列进行检测，图 所示为具体的异常检测阶段流程。第 期指挥控制与仿真 图 异常检测阶段流程图 首先，将含有 条航迹序列的实时航迹数据集，输 入 训 练 好 的 模型，模型逐个遍历检测数据集 中所有航迹序列；然后，模型将输出 对应的重构损失序列（，），若序列中不存在重构损失值超出航迹点重构损失门限，即该航迹序列中不包含异常航迹点向量，则判定为正常航迹序列，继续遍历检测下一个航迹序列，反之，若重构损失序列存在异常航迹点向量，则记录异常航迹点向量个数，再计算异常航迹点向量和航迹中所有航迹点向量之比，若该比值小于异常点占比阈值，则判定为正常航迹序列，否则，判定为异常航迹序列；之后，对异常航迹序列及其对应的重构损失序列分别进行动态可视化，通过比较各航迹序列中各特征和重构损失随时间的变化情况，对航迹的异常类型进行判断；最后，完成对数据集中所有实时航迹序列的检测。整个工作流程中，训练数据集无须人工标注，模型输入无标签数据进行无监督训练，允许训练数据中存在少量的异常航迹序列，适用于正常样本和异常样本数量悬殊的情况。实验与分析首先，对本方法的模型训练阶段进行展示，采用无标签的历史航迹数据集作为训练集对 模型进行训练，记录并展示模型在训练过程中目标函数值的变化情况，查看重构损失分布情况，并确定航迹点重构损失门限；然后，将重构损失置信度设为 ，异常航迹点占比阈值 设为，分别采用不同的重构损失计算方法计算对应的重构损失序列并分别划定阈值上下限，对测试集进行异常检测实验，探究重构损失计算方法对异常检测模型性能的影响；接着，复现现有的基于神经网络模型数据重构的异常检测方法，使用本实验数据集进行异常检测实验，与本方法模型性能进行对比；之后，在上一步实验的基础上，探究 编码器隐藏神经元个数对 模型性能的影响；最后，探究重构损失置信度的变化对模型性能的影响。实验条件）设备条件实验设备为 台联想 工作站，内存为 ，主频为 ，两张显存各为 的 显卡。在 位操作系统上，使用 编程语言，在 深度学习框架中搭建 模型。）数据集选取本实验采用从 年 月 日至当年 月 日间的 数据作为实验对象，将地理空间范围限定在经度（，）和纬度（，）的美国墨西哥湾海域上，将超出 数据正常属性范围的数据进行剔除，再从中筛选包含客船、油船和货船三种船型的共计 条航迹数据作为实验数据集。按照时间先后顺序截选，本实验将 年 月 日至 月 日的航迹数据作为训练集，月 日至 日的航迹数据作为测试集，具体划分如表 所示。表 实验数据集 舰船类型训练集测试集正常航迹异常航迹客船 油船货船 总计 为了评价无监督异常检测方法的性能，必须对测试集数据进行标注，划分正常和异常的航迹序列，以检验方法的异常检测能力。李 磊，等：基于 的无监督航迹异常检测方法第 卷其中，训练集包含 条航迹序列数据，无监督模型的训练无须对数据进行标注，经预处理后输入神经网络模型进行无监督训练。测试数据集共包含 条航迹序列。经标注，测试集包含正常航迹数据 条，异常航迹数据 条，按照舰船类型划分为三组测试集，模型性能指标参数将在三组测试集中取平均值。测试集的标注采用以下两种方式：）采用无监督的隔离森林算法对数据进行划分，标注正常和异常标签；）人工修改航迹序列的位置属性或运动属性，修改后的航迹数据均标注为异常数据。经标注后的异常航迹序列数据包含 节定义的所有航迹异常类型。）神经网络模型参数设置为确保实验不受其他因素的干扰，对 模型参数进行统一设置。神经网络模型的相关参数如表 所示。其中，梯度优化函数采用 优化算法，学习率设置为 ，模型的目标函数按照式（）设置。表 模型参数设置 参数设置值数据输入维度训练集：验证集（比值）批尺寸 激活函数优化函数学习率 模型目标函数式（）由于需要探究 隐含层神经元个数对模型性能的影响，故在前面的实验中，将 隐含层神经元个数设置为固定值。评估指标采用查准率、查全率 和 分数对异常检测模型性能进行评估，如式（）（）所示。（）（）（）其中，表示被模型预测为异常的异常样本，表示被模型预测为正常的正常样本，表示被模型预测为异常的正常样本，表示被模型预测为正常的异常样本。采用虚警率 和漏警率 对异常检测方法性能进行评估，如式（）、（）所示。（）（）实验过程）模型训练阶段本节对本方法的模型训练阶段进行展示，并示例说明重构损失的阈值上下限的划定方法。先在 深度学习框架下搭建 网络模型，后将训练集输入模型进行训练。图 展示了模型训练过程中，训练集和验证集目标函数值随训练轮数的变化情况，可见模型在 轮内稳定收敛，对训练数据集进行了充分拟合。图 模型损失值变化 图 利用高斯核密度估计曲线和直方图对三种船型训练集的重构损失序列的分布情况进行展示。下面以计算客船训练集的航迹点重构损失门限过程为例。对客船航迹数据集的重构损失序列进行对数变换后，模型能够通过正态性检验，且采用平均绝对误差计算方法得到的重构损失序列分布均可拟合正态分布，下面采用正态分布法确定航迹点重构损失门限。本文中置信度设定为 ，以采用绝对平均误差计算方法计算重构损失序列为例，输出的重构损失序列中共有 个航迹点向量，重构损失均值为 ，标准差为 ，置信度为 时，查标准正态分布表有 ，根据本方法可算得重构损失的置信区间为 ，阈值上限为 ，阈值下限为 ，表示参与检测的航迹序列中，每个航迹点向量所对应的重构损失值有 的概率会落在区间 ，中。）重构损失计算方法对模型性能的影响为选择更合适的重构损失计算方法，实验探究了重构损失计算方法对模型异常检测性能的影响。除平第 期指挥控制与仿真 图 三种船型训练集的重构损失序列直方图和核密度曲线展示 均绝对误差计算方法外，分别使用平均绝对百分比误差和均方误差计算重构损失，分别如式（）和式（）所示。（）（）（）其中，和 是重构航迹序列和原始航迹序列，和 分别是第 个时间戳上的重构和原始航迹点向量，是航迹序列中的航迹点个数，是航迹点向量中的特征个数。实验将本方法工作流程中的异常航迹点占比阈值 设为，即只要存在异常航迹点向量则将该航迹序列判定为异常航迹。在航迹点重构损失置信度为 的条件下，分别记录模型在三组测试数据集上的平均查准率、平均查全率和平均 分数，如表 所示。表 模型性能（单位：）（）计算方法查准率召回率 得分平均绝对百分比误差 均方误差 平均绝对误差 表 记录了在不同重构损失计算方法下的模型平均性能，采用平均绝对误差计算重构损失的模型平均查准率、平均查全率和平均 分数分别比采用均方误差计算方法的模型高 、和 。可见，基于该实验数据集，本方法采用平均绝对距离计算重构损失的模型性能最好。）类似方法性能比较为了验证本方法在性能上的优势，本节对现有的基于数据重构模型的异常检测方法进行了复现，表 记录了性能对比情况。表 基于数据重构模型的异常检测方法性能对比 方法 模型平均查准率 平均查全率 平均 分数 本文方法 如表 所示，文献，分别提出了基于自编码器网络、深度长短时记忆网络以及将两种网络相结合的方法，但在本实验数据集上的实验性能相对较差，模型的 分数均未超过。文献、文献和本文方法采用了变分自编码器后，模型性能提升明显，分数均在以上。本方法将 门控循环单元作为变分自编码器的输入层和输出层，使得本文提出的 模型在查准率、查全率和模型 分数上分别比文献提出的 模型高 、和 ，可见引入 模块可有效提高模型的异常检测性能。综上，本方法采用的 模型在解决航迹异常检测问题上具有性能优势。）编码器隐藏神经元个数对模型性能的影响对比为了对比 编码器隐藏神经元个数对模型性能的影响，基于上一步实验得到的性能最佳模型，本节采用平均绝对误差计算重构损失，并在重构损失置信 李 磊，等：基于 的无监督航迹异常检测方法第 卷度为 的条件下划定航迹点重构损失门限，在本实验数据集中进行异常检测实验，探究 编码器隐藏神经元个数变化对模型性能的影响。实验得到的模型平均 分数对比情况如图 所示。图 平均 分数随 隐含层神经元个数变化情况 由图 可见，当 编码器隐藏神经元个数由 开始逐渐增加时，模型平均 分数随之提高，可见适当增加 编码器隐藏神经元个数对模型性能具有提升作用。当数目为 时，模型性能最好，分数高达 。隐含层神经元个数继续增多，使得模型参与训练的参数增多，无法很好地拟合数据集，导致模型平均 分数逐渐降低。当个数为 时模型平均 分数为 ，比隐含层神经元个数为 时的模型得分还要低 。）重构损失置信度对模型性能的影响根据本方法提出的异常检测方法流程，理论上重构损失置信度越高，其门限则越大，对应的查准率随置信度增大而增大，查全率对应减小。为了进一步探究重构损失置信度的设置对模型性能的影响，本节实验对不同重构损失置信度下的模型性能进行对比，进一步总结本方法的实际应用要点。本节实验将异常航迹点占比阈值 设为，分别在置信度为、和 的条件下设定航迹点重构损失门限，对三组测试数据集进行异常检测实验，记录模型的平均查准率、平均查全率、平均 分数、平均漏警率和平均虚警率的变化情况，实验结果如图 和 所示。根据以上实验结果，可得出以下结论：模型的 分数随着置信度的增大而增大，在重构损失置信度为 时 分数最高，为 。结果表明，设置的重构损失置信度越大，模型的异常检测性能越好。根据本文提出的异常检测方法，理论上设置的图 模型查准率、查全率和 分数随置信度的变化情况 ，图 模型虚警率和漏警率随置信度变化情况 航迹点重构损失置信度越高，其重构损失门限则越大，对应的查准率随置信度增大而增大，查全率随置信度增大而减小，实验结果验证了这一点。模型平均漏警率和重构损失置信度呈正相关，所有参与实验的置信度对应的平均漏警率均在 以下，可见本模型对异常航迹的检测较为灵敏。模型平均虚警率和重构损失置信度呈负相关，然而最大和最小虚警率之间相差 ，可见重构损失置信度的设置对模型虚警率影响较大。综上，针对一线侦察应用场景，虽然重构损失置信度过低会使方法的虚警率较高，耗费人员精力对检测出的异常航迹进行甄别，但本方法较低的漏警率可使在一线侦察业务流程中不错过任何一个可疑的运动目标，及时从异常航迹中发现战场态势变化，符合一线单位侦察需求。）异常航迹类型研判过程本节选取两条在前文异常检测实验中被判定为异常的航迹作为实验素材，展示本方法对异常航迹的研判过程。依据前文关于航迹异常检测阶段的理论描第 期指挥控制与仿真 述，重构损失置信度设置为 ，将异常航迹点占比阈值 设为，即存在异常航迹点就将其判定为异常航迹序列，对异常航迹进行深入研判。按照航迹的异常类型，分别对不同类型的异常航迹研判过程进行展示。运动异常航迹的检测在前文实验中，号为 的货船航迹被检测为异常航迹，异常类型的研判过程分为以下三步。首先，对该航迹序列的重构损失和向量特征随时间的变化情况进行可视化，结果如图 所示。红色方框标注了超出航迹点重构损失门限的部分，对应可以查看超出部分所处的时间段。同时展示的航迹序列所有向量特征随时间的变化情况无法发现导致重构损失值超限的对应特征。图 重构损失和向量特征随时间变化情况 然后，对该船航迹序列和历史航迹数据集中具有相近航道的航迹序列在地图上可视化展示，结果如图 所示。红色航迹线是该船在地图上呈现的实际航迹，绿色区域是历史航迹数据集中具有相近航道的航迹覆盖区域。由图 可知，该货船航迹基本在绿色区域中活动，与历史航迹基本吻合，无明显位置异常迹象。最后，操纵动态可视化图，逐个展示航迹向量特征随时间的变化情况。图 为航向变化率随时间的变化情况，可发现航向变化率曲线的异常突起时间段和重构损失超限时间段相一致。综上，可以判定该货船航迹序列为运动异常航迹序列，异常原因为航向在短时间内的快速变化。位置异常航迹的检测在前文实验中，号为 的货船航迹被检测为异常航迹，异常类型的研判过程同上节。图 运动异常航迹的可视化 首先，通过图 查看了该航迹序列的重构损失和向量特征随时间的变化情况，可以看到，在 年 月 日晚 点 分至次日的 点 分，重构损失超出了阈值上限。然后，查看航迹序列在地图上的可视化情况，如图 李 磊，等：基于 的无监督航迹异常检测方法第 卷图 航向变化率随时间的变化情况 图 重构损失和向量特征随时间变化情况 所示，并和历史航迹数据集中具有相近航道的航迹进行对比，可以看到，有部分航迹点位置偏离历史航迹所经区域，初步可判定为是位置异常航迹。最后，通过图 查看各航迹向量特征值随时间的变化情况，经度的变化和重构损失的超限时间范围相吻合。综上研判，该航迹序列判定为位置异常航迹，异常表现为偏离历史航迹所在区域。图 位置异常航迹的可视化 第 期指挥控制与仿真 图 经度随时间的变化情况 结束语针对海量无异常信息标签航迹数据的目标行为异常检测问题，提出一种基于 模型的无监督航迹异常检测方法。该方法引入 门控循环单元提升 变分自编码器模型的时序建模能力，利用无异常信息标签的历史航迹数据对模型进行训练，采用正态分布法或百分位数法将指定置信度下的置信区间作为航迹点重构损失门限，结合实时航迹序列的重构损失值超限和异常航迹点占比情况检测出异常航迹。通过实验证明本方法可实现对无行为模式标签的航迹数据进行异常检测，与现有的基于数据重构模型的无监督航迹异常检测方法相比，具有性能优势，对异常航迹的检测具有高灵敏度和低漏警率，基本满足战场态势认知需求。参考文献：胡晓峰，贺筱媛，饶德虎，等 基于复杂网络的体系作战指挥与协同机理分析方法研究 指挥与控制学报，（）：，（）：崔令飞，郭永红，邵鹏志，等 运用时空大数据增强战场态 势 感 知 能 力 国 防 科 技，（）：，（）：段玉先，刘昌云，魏文凤 战场态势感知关键技术研究 火力与指挥控制，（）：，（）：，（）：，（）：，：，：杨帆，何正伟，何帆 基于 神经网络的船舶异常行为检测方法 武汉理工大学学报（交通科学与工程版），（）：，（），（）：，：，李 磊，等：基于 的无监督航迹异常检测方法第 卷 ，（）：，：彭井桐，祝永新，汪辉，等 基于 的 神经网络飞行数据异常检测 微电子学与计算机，（）：，（）：，：方积乾，陆盈 现代医学统计学 北京：人民卫生出版社，：，何晓群 现代统计分析方法与应用 北京：中国人民大学出版社，：，（）：，（）：，（），：（责任编辑：张培培）