Formal Method for Behavior Verification and Data Validation of Station Interlocking System

随着中速磁悬浮列车的广泛商业化应用^[1]，作为其关键核心技术的电磁悬浮控制技术已经成熟稳定^[2]. 商业化运营过程中随着规模和复杂性的迅速增加，对中速磁悬浮列车悬浮系统的安全性和可靠性的需求不断增长. 在过去几十年中，对磁浮列车悬浮系统的性能监测尤其是异常检测的研究引起了极大的兴趣^[3-4].

磁浮列车悬浮系统的异常指在列车实际运行中，悬浮系统出现工作状态与期望状态不匹配，但没有导致故障或失效的情况. 悬浮系统的异常具有出现时刻未知，持续时间短，通过控制器调节可以恢复稳定的特点. 对悬浮系统异常的精确检测可以提前获取系统的运行情况，有利于提前安排预防性维护计划，降低工程事故率^[5].

针对磁悬浮列车的异常检测研究通常可以分为针对轨道的检测和针对悬浮系统的检测. 针对磁浮列车轨道的异常检测问题，Zhang等^[6]提出基于小波与分形相结合的高速磁悬浮列车长轨道磁异常检测方法，以发现高速磁悬浮列车长定子铁芯的局部短路故障；Deng等^[7]对高温超导磁悬浮列车运行时永磁导轨的不规则性进行测量和表征，并设计了一种不规则性检测设备；罗茹丹^[8]提出了“2+1”的检测方法实现对长定子行波主漏磁场异常情况的动态检测；林国斌等^[9]提出一种基于车轨状态监测的悬浮冗余控制系统，用于在控制列车稳定悬浮的同时对车轨状态进行检测，并通过对悬浮间隙和电磁铁振动情况的分类和学习确定状态类型，具有实时性好、容错性高等优点；杨杰等^[10]提出一种基于悬挂式磁悬浮轨道交通系统的轨道维护设备，在对轨道状态进行取样和检测后对轨道异常实现诊断. 文献[11]收集了振动条件下高温超导钉扎系统的水平和垂直加速度数据.

针对磁浮列车悬浮系统的异常检测问题，朱跃欧等^[12]提出一种磁浮列车悬浮控制器的异常预警方法与检测系统，通过分析悬浮间隙值所在区间，实现了车辆异常状况的在线检测与预警；王平等^[13]提出了一种基于改进典型相关分析的中低速悬浮系统异常检测方法；通过运营线数据验证了该方法能获得更高的检测率；王平等^[14]提出了一种基于超球体高斯分布的悬浮系统异常检测方法,该方法对系统异常检测的结果较为精确. 以新型悬挂式永磁磁浮列车研究对象，马政^[15]设计并实现了一套基于物联网的永磁磁浮列车远程监控系统,实现了对列车悬浮状态的异常监控；Wang等^[16]引入了支持向量机方法，利用悬浮传感器采集的实时信号来检测磁悬浮列车悬浮系统的异常状态；Ma等^[17]将磁悬浮列车悬浮控制系统收集的数据用于训练机器学习模型，并根据数据确定磁悬浮列车的实际运行状态，以达到实时监控的目的.

然而，上述方法存在2个问题：1） 大部分方法是基于单变量的异常检测，没有充分利用系统的运行数据. 磁浮列车悬浮系统的异常检测大多都是只基于悬浮间隙，而实际运行中电流、加速度数据同样包含系统的异常信息. 因此，如何应用列车悬浮系统运行时产生的多变量数据进行异常检测是一个关键问题. 2） 需要有异常的先验信息，而实际系统的异常先验信息是不完全可知的，这无疑会降低异常检测的准确率；此外，面向实际工程系统的异常检测模块或系统在设计时需考虑异常漏检率和异常误检率的折中，这需要从数理的角度进行分析.

基于上述分析，本文将基于参数化残差的异常检测方法引入悬浮系统的异常检测问题，以最低异常漏检率为设计目标，提出了磁悬浮列车悬浮系统基于最低漏检率的残差生成器设计方法及异常检测方法. 结合磁悬浮列车悬浮系统的运行数据，验证了该残差生成器能够在满足期望误报率的情况下，实现异常漏检率的最小化.

1.   残差生成器与异常检测

1.1   参数化残差生成器

系统出现异常时的状态空间模型^[18]为

式中：${{\boldsymbol{x}}}(k)$为时刻$k$的系统状态；${{\boldsymbol{u}}}(k)$为时刻$k$的系统输入；${{\boldsymbol{y}}}(k)$为时刻$k$的系统输出；A、B、C、D为一定维数的系统参数矩阵；${{{{\boldsymbol{B}}}}_{\rm{f}}}$和${{{{\boldsymbol{D}}}}_{\rm{f}}}$为异常参数矩阵；${\boldsymbol{f}}(k)$为异常向量；${{\boldsymbol{w}}}(k)$为系统的过程噪声；${{\boldsymbol{v}}}(k)$为系统的输出噪声.

则状态检测量${{\boldsymbol{z}}}(k)$为

式中：$s$为正整数，表示截断数据的长度；${{{\boldsymbol{u}}}_s}(k)$、${\boldsymbol{y}}_s(k)$为长度为$s + 1$的输入、输出数据矩阵；${{{{\boldsymbol{\varPsi}} }}_{{s}}}$为基于数据驱动的系统信息矩阵^[19]；${{{{\boldsymbol{\varPsi}} }}_{{s}}}^ \bot$为${{{{\boldsymbol{\varPsi }}}}_{{s}}}$的正交补矩阵.

当系统出现异常时的状态检测量为

式中：${{{{\boldsymbol{H}}}}_{{\rm{w}},s}}$为过程噪声的Hankel矩阵；${{{{\boldsymbol{W}}}}_{{\rm{k}},s}}$为数据驱动的过程噪声矩阵；${{{{\boldsymbol{V}}}}_{{\rm{k}},s}}$为数据驱动的过程噪声矩阵；${{{H}}_{{\rm{f}},s}}$为异常向量${\boldsymbol{f}}(k)$的Hankel矩阵；${{{{\boldsymbol{F}}}}_s}$为基于数据驱动的异常矩阵.

从式（3）可以看出，当系统出现异常时，${{\boldsymbol{z}}}(k)$不仅包含噪声信息，同时也包含了异常信息^[20].

在基于模型的异常检测或故障检测领域，残差是通过比较实际系统的输出${{\boldsymbol{y}}}(k)$与估计输出${\hat {\boldsymbol{y}}}(k)$得到的，即定义${{\boldsymbol{r}}}(k) = {{\boldsymbol{y}}}(k) - {\hat {\boldsymbol{y}}}(k)$. 残差数据中包含了系统扰动和系统异常的全部信息，残差产生器设计的主要目标是设法让异常检测系统对异常敏感，同时对其他信号具有鲁棒性.

结合Youla参数化和基于输入输出数据的状态检测量，残差生成器参数化形式为

式中：${{\boldsymbol{g}}}$为参数向量，是残差生成器和异常检测的核心^[21].

出于实际工程运营需要，通常在系统运行前就需要完成${{\boldsymbol{g}}}$与${{{{\boldsymbol{\varPsi }}}}_{{{{s}}}}}^ \bot$的设计，并且在一段时间内都不能改变其参数.

1.2   基于残差的异常检测

由于残差数据中包含了系统扰动和系统异常的全部信息，因此，在系统运行时，可以通过对残差进行分析、处理，以及基于某种准则或阈值，来进一步掌握系统的运行状态信息，判断系统的运行状态.

当系统运行时，残差生成器实时获得残差数据${{\boldsymbol{r}}}(k)$，并将残差数据送入残差评估单元，计算残差评估函数${\boldsymbol{J}}({\boldsymbol{r}})=J_k ({\boldsymbol{r}})$. 在设置评估函数阈值${J_{{\rm{th}}}}$后，根据式（5）判断系统是否发生异常

2.   基于最低漏检率的残差生成器设计

2.1   基于均值和方差的置信集

传统的异常检测方案通常采用基于假设检验的方法进行残差评估，即假设残差的分布知识是准确已知的，或者可以从历史数据中估计出来^[22]. 然而，对于磁浮列车悬浮系统而言，使用历史过程数据获得高度可信的概率分布估计比较繁琐，经验估计与真实概率分布的偏差可能会导致异常检测结果不可靠^[23].

尽管难以获得悬浮系统的高度可信概率分布，但在无异常时的${{\boldsymbol{z}}}(k)$可以视为属于同一个集合. 考虑到悬浮系统运行数据的统计特性，可以依据数据的均值和协方差矩阵对其进行建模.

设无异常情况下，${{{\boldsymbol{z}}}}(k)$满足

式中：${S_{\rm{h}}}$为无异常数据的置信集；${\zeta _{{\rm{z}}}{\rm{}}}$为所有${\boldsymbol{z}}(k)$的分布集合；${\bar {\boldsymbol{z}}_{\rm{h}}}$和${{\boldsymbol{\varSigma}} _{\rm{h}}}$为无异常时状态检测变量的均值和方差，${\bar {\boldsymbol{z}}_{\rm{h}}}$和${{\boldsymbol{\varSigma}} _{\rm{h}}}$是通过对健康运行数据进行计算得来的，通常可以选择刚投入运行时无异常的运行数据，也可以是某一次无异常运行时的数据；E为${{{\boldsymbol{z}}}}(k)$的均值；V为${{{\boldsymbol{z}}}}(k)$的方差.

在悬浮系统运行初期，由于缺乏异常数据，难以分析出异常的数据分布特性. 因此，通常会依据正常运行时的数据特点预先设计一个异常检测系统，以实现对大部分潜在异常的检测. 假设系统运行时会出现$M$个典型异常，且发生第$i$个异常时的异常信号设为${{\boldsymbol{f}}_i}(k)，i = 1,2, \cdots ,M$. 为不失一般性，假设每个异常发生的概率相同. 采用均值和协方差矩阵对异常进行建模. 设异常情况时的${{\boldsymbol{z}}}(k)$满足

式中：${S _{{{\boldsymbol{f}}_i}}}$为异常数据的置信集；${\bar {\boldsymbol{z}}_{{{\boldsymbol{f}}_i}}}$为${{{\boldsymbol{z}}}_{{{\boldsymbol{f}}_i}}}(k)$的均值，${{{\boldsymbol{z}}}_{{{\boldsymbol{f}}_i}}}(k)$为第i个异常时的状态检测变量；${S _{{{\boldsymbol{f}}_i}}}$为${{{\boldsymbol{z}}}_{{{\boldsymbol{f}}_i}}}(k)$的方差.

当有可用的异常数据时，可以直接使用该数据建立置信集. 若无异常数据或异常数据少时，则可以自定义M个异常情况下的均值和协方差，进而构建置信集，当有新可用的异常数据时，记录数据并离线更新置信集.

2.2   基于最低漏检率的异常检测

针对第 $i$个异常，基于异常检测的残差信号为

式中：${{{\boldsymbol{g}}}_i}$为第$i$个异常时的参数向量.

因此，第$i$个异常时的残差评估函数${\boldsymbol{J}}({\boldsymbol{r}})$与阈值$J_{{\rm{th}}}$为

式中：${{{\boldsymbol{r}}}_{i,h}}$为第$i$个无异常场景的残差.

定义异常误检率${P_{\rm{F}}} = \Pr \left\{ {J({\boldsymbol{r}}) > 1\left| {{{{\boldsymbol{f}}}_i}(k) = {\boldsymbol{0}}} \right.} \right\}$，定义异常漏检率为${P_{\rm{M}}} = \Pr \left\{ {J({\boldsymbol{r}}) \leqslant 1\left| {{{{\boldsymbol{f}}}_i}(k) \ne {\boldsymbol{0}}} \right.} \right\}$. 理想情况下，异常检测系统应该对所有潜在的异常检测且不存在异常的误报，即${P_{\rm{F}}} = {P_{\rm{M}}} = 0$，在实际中却难以实现. 当设计的${P_{\rm{M}}}$越小时，往往会导致${P_{\rm{F}}}$较高，反之亦然. 因此，异常检测系统通常会采取${P_{\rm{F}}}$和${P_{\rm{M}}}$的折中.

本文研究的悬浮系统的最低异常漏检率可被表述为${P_{\rm{F}}}$固定时的最小化${P_{\rm{M}}}$，即考虑当${P_{\rm{F}}}$固定且可接受时，如何使得${P_{\rm{M}}}$最小化^[24]. 根据式（8）可知，${{{\boldsymbol{g}}}_i}$是${{\boldsymbol{r}}_i}(k)$生成的关键，${{{\boldsymbol{g}}}_i}$的取值会影响异常检测的${P_{\rm{F}}}$和${P_{\rm{M}}}$. 因此，基于最低漏检率的残差生成器的设计最终转化为求得参数向量${{{\boldsymbol{g}}}_i}({{{\boldsymbol{g}}}_i} \ne {\boldsymbol{0}})$，使得

式中：${\alpha _i} \in (0,1)$，${\beta _i} \in (0,1)$，分别为${P_{\rm{M}}}$和${P_{\rm{F}}}$的上界， ${\beta _i}$在设计时应预先给出^[25].

进一步，异常检测式（10）可以被重构为

式中：${k_{\beta_i}} = \sqrt {1/{\beta _i}}$；${\tilde{\boldsymbol{z}}}_{{\boldsymbol{f}}_i} = {\bar {\boldsymbol{z}}_{{\boldsymbol{f}}_i}} - {\bar {\boldsymbol{z}}_{\rm{h}}}$.

当获得最优${{{\boldsymbol{g}}}_i}$后，只需将${{{\boldsymbol{g}}}_i}$乘以一个非零常数，即可满足${k_{\beta_i}}\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}{\varSigma _{\rm{h}}}{{{\boldsymbol{g}}}_i}} \leqslant 1$. 因此，先分析$a({{{\boldsymbol{g}}}_i}) > 0$时$a({{{\boldsymbol{g}}}_i})$的最大值.

令

式中：${{\boldsymbol{\varSigma }} _{{\boldsymbol{f}}_i}} \ne 0$；$b({{{\boldsymbol{g}}}_i}) > c({{{\boldsymbol{g}}}_i}) > 0$.

则$d({{{\boldsymbol{g}}}_i}) = c({{{\boldsymbol{g}}}_i})/b({{{\boldsymbol{g}}}_i}) < 1$，因此，可以得到式（13）的关系.

对$e({{{\boldsymbol{g}}}_i})$求关于$a({{{\boldsymbol{g}}}_i})$的导数，得

因此，$e({{\boldsymbol{g}}_i})$关于$a({{{\boldsymbol{g}}}_i})$是单调增加的，所以式（11）可以变为求${{{\boldsymbol{g}}}_i} = \arg \mathop {\max }\limits_{{{\boldsymbol{g}}_i} \ne 0} e({{{\boldsymbol{g}}}_i})$的问题^[26]，即

通过上述分析，基于最低漏检率的残差生成器的设计问题最终转换为求解式（15）的广义特征值和特征向量^[27]，具体求解步骤如算法1所示.

1） 算法 1：最优参数向量的求解

步骤 1　对$\sqrt {{{{\boldsymbol{\varSigma }}} _{{{\boldsymbol{f}}_i}}}}$做 SVD 分解，$\sqrt {{{{\boldsymbol{\varSigma }}} _{{\boldsymbol{f}}_i}}} = {{{{\boldsymbol{U}}}}_i}\left[ { {{{{{\boldsymbol{S}}}}_i}}\;\; 0 } \right] {{{\boldsymbol{V}}}}_i^{\rm{T}}$，其中，$\lambda _{{\rm{max}}}$为最大奇异值；${\boldsymbol{U}}_i、{\boldsymbol{S}}_i、{\boldsymbol{V}}_i$均为SVD分解后得出的矩阵；${\boldsymbol{v}}_i$为${\boldsymbol{V}}_i$中的元素向量.

步骤 2　通过求解方程组

得出${{{\boldsymbol{v}}}_i}$.

步骤 3　依据 ${{{\boldsymbol{g}}}_i} = {{{{\boldsymbol{U}}}}_i}{{{\boldsymbol{S}}}}_i^{ - 1}{{{\boldsymbol{v}}}_i}$计算${{{\boldsymbol{g}}}_i}$.

步骤 4　依据 ${{\hat {\boldsymbol{g}}}_i} = \dfrac{{{{\boldsymbol{g}}}_i}}{{{k_{{\beta _i}}}\sqrt {{{\boldsymbol{g}}}_i^{\rm{T}}{{{{\boldsymbol{\varSigma}}}} _{\rm{h}}}{{{\boldsymbol{g}}}_i}} }}$计算出满足式（10）的最优参数向量${{\hat {\boldsymbol{g}}}_i}$.

在得出最优参数向量后，第$i$个异常的参数化残差生成器可以通过式（16）构建.

令${{\hat {\boldsymbol{r}}}_{i,{\rm{h}}}} = {{\hat {\boldsymbol{g}}}_i}^{\rm{T}}{\bar {\boldsymbol{z}}_{\rm{h}}}$，异常检测器的残差评估函数和异常检测阈值为

完整的异常检测见算法 2

2） 算法 2：磁悬浮列车悬浮系统异常检测算法

① 离线辨识

步骤 1　采集无异常情况下的输入输出数据，基于子空间辨识法离线辨识出矩阵${{{{\boldsymbol{\varPsi}} }}_{\boldsymbol{s}}}^ \bot$.

步骤 2　利用无异常的输入输出数据，计算得到${\bar {\boldsymbol{z}}_{\rm{h}}}$和${{\boldsymbol{\varSigma}} _{\rm{h}}}$，构建健康状况下的置信集${S_{\rm{h}}}$.

步骤 3　建立M个异常的置信集${S_{{{\boldsymbol{f}}_i}}}$. 当有M个异常时的输入输出数据，计算每个异常的${\bar {\boldsymbol{z}}_{{{\boldsymbol{f}}_i}}}$和${{{{\boldsymbol{\varSigma}} }}_{{\boldsymbol{f}}_i}}$，建立${S _{{{\boldsymbol{f}}_i}}}$. 如果无异常数据，预设每个异常的${\bar {\boldsymbol{z}}_{{{\boldsymbol{f}}_i}}}$和${{{{\boldsymbol{\varSigma}}}} _{{{\boldsymbol{f}}}_i}}$，依此建立${S _{{{\boldsymbol{f}}_i}}}$.

步骤 4　确定每个异常下的${\beta _i}$，根据算法1计算出${{\hat {\boldsymbol{g}}}_i}$和${{\hat {\boldsymbol{r}}}_{i,h}}$.

② 在线异常检测

步骤 5　在线收集时间间隔为$\left[ {k-s,k} \right]$的输入输出数据，并构造数据矩阵$\left[ {{{{\boldsymbol{u}}}_s}(k)}\quad{{{{\boldsymbol{y}}}_s}(k)} \right]^{\rm{T}}$.

步骤 6　依据式（16）计算${{{\boldsymbol{r}}}_i}\left( k \right)$，依据式（17）执行异常检测并输出异常检测的结果.

步骤 7　令$k = k + 1$，重复步骤5和步骤6，直到中止时刻.

3.   异常检测实验验证

本节结合长沙磁浮快线某辆磁浮列车的悬浮系统的健康数据和异常数据，验证所提出的基于最低漏检率的残差生成器在悬浮系统异常检测中的有效性.

3.1   悬浮系统的3类异常

以长沙磁浮快线的运营磁浮列车为例，如图1所示，作为典型的常导电磁型磁浮列车^[28]，通常选用间隙作为其悬浮系统异常检测的检测变量，通过比较间隙值和异常阈值的大小来判断是否出现异常. 磁浮列车悬浮系统的工程异常阈值设置为标准间隙值 ±4 mm. 该线路的标准悬浮工作间隙为9 mm，工程异常阈值分别设为5 mm和13 mm.

图  1  长沙磁浮快线

Figure  1.  Changsha Maglev Express

下载: 全尺寸图片 幻灯片

磁浮列车悬浮系统经常会遭受到3类异常：间隙突变异常、砸轨异常和加速度传感器异常. 间隙突变异常的现象是指在运行中间隙有明显的波动，但没有超过工程异常阈值，此时加速度、电流也对应有一定的波动，出现两个尖峰，但电压几乎没有波动. 该异常通常是由于磁浮列车经过轨道接缝引起的，如图2所示.

图  2  间隙突变异常时间隙、电流、加速度和电压

Figure  2.  Gap, current, acceleration, and voltage at gap mutation anomaly

下载: 全尺寸图片 幻灯片

砸轨异常的现象是间隙有剧烈波动，在一段时间内间隙值超过工程异常阈值并且需要经历5 s左右的恢复期，此时加速度、电流、电压均在一段时间内出现剧烈波动，如图3所示.

图  3  砸轨异常时的间隙、电流、加速度和电压

Figure  3.  Gap, current, acceleration, and voltage at rail smashing anomaly

下载: 全尺寸图片 幻灯片

加速度计异常的现象为间隙没有超过工程异常阈值，也没有较为明显的波动，但加速度的波动十分明显且变为负值，电压也出现两个极小值尖峰. 由于实际加速度计存在着1g的偏置，当加速度计输出信号为负时，表明此刻的加速度已经低于$- 10\;{\rm{m/{s^2}}}$，如图4所示.

图  4  加速度计异常时的间隙、电流、加速度和电压

Figure  4.  Gap, current, acceleration, and voltage at accelerometer anomaly

下载: 全尺寸图片 幻灯片

3.2   异常检测验证与分析

以磁浮列车悬浮系统的第13个悬浮点的运行数据作为检测对象，如图5所示.

图  5  第13个悬浮点的运行数据

Figure  5.  Operational data of the 13th suspension point

下载: 全尺寸图片 幻灯片

列车正线行驶时存在着站内悬浮和站间行驶两个运行工况，考虑到上述3类异常均出现在站间行驶运行工况，因此，以速度为划分基准，分段提取站间行驶工况的运行数据. 本节选择一个无异常站间行驶数据段、3个包含间隙异常的数据段、2个包含砸轨异常的运行数据段和包含3个加速度传感器异常的数据段作为检测数据.

选择悬浮系统的电流数据作为输入数据，间隙数据和加速度数据作为输出数据，如图6所示. 依据正常运行时的输入输出数据建立${S_{\rm{h}}}$并离线辨识出${{{{\boldsymbol{\varPsi}} }}_s}^ \bot$.

图  6  辨识所用的输入输出数据

Figure  6.  Input and output data used for identification

下载: 全尺寸图片 幻灯片

依据3种异常的输入输出数据建立异常置信集${S _{{{\boldsymbol{f}}_i}}}$. 出于工程运营需求，设期望的故障误报率${P_{\rm{F}}} = 5{\text{%}}$，执行算法1可得

基于最低漏检率的残差生成器为

依据式（17）算出9个检测数据段的$J({{\boldsymbol{r}}_i})$并送入检测单元. 图7为无异常数据段的检测结果. 从图中可以看出，当系统无异常时，$J({{\boldsymbol{r}}_i})$小于0.4，均在异常检测阈值${J_{{\rm{th}}}}$之下，不存在异常误检.

图  7  正常运行数据段的检测结果

Figure  7.  Detection results of normal operational data segment

下载: 全尺寸图片 幻灯片

含间隙突变异常的3个数据段的检测结果如图8所示. 在每个数据段，大部分数据都低于异常检测阈值${J_{{\rm{th}}}}$，但是在出现异常的时刻，存在着$J({{{\boldsymbol{r}}}_i})$大于${J_{{\rm{th}}}}$的情况，个数分别为8、6和5，说明该方案能够实现对间隙异常的检测. 此外，间隙突变异常的$J({{{\boldsymbol{r}}}_i})$的峰值略高于检测阈值${J_{{\rm{th}}}}$，一般在2以下，且不存在第二个大于${J_{{\rm{th}}}}$的次尖峰. 表1为间隙突变异常的检测位置与异常实际发生位置的对比，从表中可以看出，间隙突变异常1的检测时刻落后实际异常发生时刻 0.2 s，其他两个数据段可以及时检测出异常.

图  8  间隙突变异常的检测结果

Figure  8.  Detection results of gap mutation anomalies

下载: 全尺寸图片 幻灯片

表  1  间隙突变异常的检测位置与异常实际发生位置

Table  1.  Detection position and actual occurrence position of gap mutation anomalies

异常	检测位置/s	实际位置/s
间隙异常 1	5.9	5.7
间隙异常 2	568.7	568.5
间隙异常 3	21.1	21.1

下载: 导出CSV 

| 显示表格

含砸轨异常的两个数据段的检测结果如图9所示. 从图中可以看出，该方法能够实现对砸轨异常的检测. 图9（a）中有20个样本的$J({{{\boldsymbol{r}}}_i})$>${J_{{\rm{th}}}}$，图9（b）中有20个样本的$J({{{\boldsymbol{r}}}_i})$>${J_{{\rm{th}}}}$，此外，砸轨异常位置的$J({{{\boldsymbol{r}}}_i})$存在着两个高于${J_{{\rm{th}}}}$的峰值. 表2为砸轨异常的检测位置与异常实际发生位置对比. 从表2可以看出，砸轨异常的检测时刻落后于实际异常的发生时刻，落后约0.2 s. 结合图3可以发现，造成这一现象的主要原因是砸轨异常发生时刻数据变化不明显，间隙和加速度都有一个缓慢减小的阶段，没有大范围变化.

图  9  砸轨异常的检测结果

Figure  9.  Detection results of rail smashing anomalies

下载: 全尺寸图片 幻灯片

表  2  砸轨异常的检测位置与异常实际发生位置

Table  2.  Detection position and actual occurrence position of rail smashing anomalies

异常	检测位置/s	实际位置/s
砸轨异常 1	7.1	7.0
砸轨异常 2	20.2	20.0

下载: 导出CSV 

| 显示表格

含加速度计异常的3个数据段的检测结果如图10所示. 在每个数据段，大部分数据都低于异常检测阈值${J_{{\rm{th}}}}$，但是在出现加速度计异常的时刻，存在着$J({{{\boldsymbol{r}}}_i})$>${J_{{\rm{th}}}}$的情况，个数分别为5、7和4，且$J({{{\boldsymbol{r}}}_i})$的值较大，说明该方案能够实现对加速度计异常的检测. 表3为加速度计异常的检测位置与异常实际发生位置对比. 从表3中可以看出，一旦出现加速度计异常，本文所提出的方法可以立即检测出该异常，及时性较好.

图  10  加速度计异常的检测结果

Figure  10.  Detection results of accelerometer anomalies

下载: 全尺寸图片 幻灯片

表  3  加速度计异常的检测位置与异常实际发生位置

Table  3.  Detection position and actual occurrence position of accelerometer anomalies

异常	检测位置/s	实际位置/s
加速度计异常 1	52.2	52.2
加速度计异常 2	52.1	52.1
加速度计异常 3	40.5	40.5

下载: 导出CSV 

| 显示表格

经过上述分析，基于残差生成器的异常检测方法可以用于磁悬浮列车悬浮系统的异常检测. 基于最低漏检率的异常检测方法能够实现磁悬浮列车悬浮系统3类典型异常的完全检测，不存在漏检的情况且不存在误检.

4.   结　论

本文针对电磁悬浮系统的异常检测问题开展了研究，以基于参数化残差的异常检测为核心，依据数据的均值和方差建立置信集，设计了基于最低异常漏检率的残差生成器，依据生成的残差数据在线进行异常检测. 最后，以长沙磁浮快线的运行数据为例，对悬浮系统常见的3类异常进行检测，实验结果表明，在给定异常误报率为5%时，本文所提出的方法能够实现对3类典型异常的全部检测，没有异常的漏检和误检. 检测时间对比结果进一步说明，本文提出的方法检测实时性较好.