基于观测模块的贝叶斯网络轨道客车牵引系统故障诊断方法及系统技术方案

本发明专利技术涉及故障分析技术领域，公开一种基于观测模块的贝叶斯网络轨道客车牵引系统故障诊断方法及系统，以利用多个观测点的信号实现精准故障诊断。方法包括：选取观测点，对各观测点单独从相对应的初始数据集中提取时频域特征并降维得到特征变量；将各观测点的状态变量、降维后的特征变量和故障模式变量作为贝叶斯网络模型的输入，根据各特征变量和故障模式变量与观测点状态变量之间的关联矩阵将所有输入划分成与观测点数量相等个数的观测模块；学习各观测模块的最优图，并将各最优图合并成一个模块化贝叶斯网络模型；实时采集各观测点的传感器数据，利用推理引擎更新模块化贝叶斯网络模型各节点的后验概率，根据各节点的后验概率进行故障判断。

【技术实现步骤摘要】
基于观测模块的贝叶斯网络轨道客车牵引系统故障诊断方法及系统
本专利技术涉及牵引传动系统故障分析
，尤其涉及一种基于观测模块的贝叶斯网络轨道客车牵引系统故障诊断方法及系统。
技术介绍
牵引传动系统作为轨道客车信息控制系统的关键系统之一，长时间地在高温高压等恶劣环境下运行，承受轨道不平顺的激励与自身引起的各种振动、碰撞和冲击，是高速列车的高发故障源。因此，对轨道客车牵引系统进行实时故障诊断能够有效监控列车运行状态、预测列车故障，可有效降低列车故障的发生率，具有重要的理论研究价值和实际应用价值。模块化贝叶斯网络的思想早已有研究者提出，但都是基于物理拓扑的模块。2014年申请的专利CN201410308885.1公开了一种专利技术名称为“模块化复杂装备贝叶斯网络故障预测方法”，该方法是根据复杂装备的物理结构进行模块划分，针对故障模式进行建模，然后将各个模块的故障模式用一个高级节点封装起来，根据共享节点将高级节点连接在一起。文献“孟倩倩.基于模块化贝叶斯网络的工业报警根源分析[D].北京化工大学,2018.”公开了一种基于大规模系统子模块划分的报警溯源分析方法。该方法仅根据流程工业生产系统的生产单元进行模块划分，模块化贝叶斯网络的网络节点只针对系统单元。目前对于牵引传动系统的故障诊断大多基于单个观测点的信息，然而牵引系统包含不同设备和大量元部件，故障模式复杂多变、故障种类多、故障原因多，不同位置的故障类型可能对某一单个观测点具有相同的表现形式，这会导致故障诊断发生误判。因此，现急需提供一种基于多观测点且不降低诊断效率的故障诊断方法及系统。
技术实现思路
本专利技术目的在于公开一种基于观测模块的贝叶斯网络轨道客车牵引系统故障诊断方法及系统，以在不降低诊断效率的条件下利用多个观测点的信号实现牵引传动系统的故障诊断。为实现上述目的，本专利技术公开一种基于观测模块的贝叶斯网络轨道客车牵引系统故障诊断方法，包括：步骤S1、选取观测点，在各所述观测点部署传感器采集各观测点包含正常状态数据和故障状态的初始数据集；对各观测点单独从相对应的初始数据集中提取时频域特征建立特征变量库D1，并对其进行降维得到特征变量库D2；步骤S2、将各观测点的状态变量、降维后的特征变量和故障模式变量作为贝叶斯网络模型的输入，建立输入的各特征变量和故障模式变量与观测点状态变量之间的关联矩阵，根据所述关联矩阵将所有输入划分成与观测点数量相等数量的观测模块；步骤S3、利用结构学习算法学习每个观测模块的最优图，并将各个观测点的最优图合并成一个完整的网络；网络各节点之间的信度条件概率参数采用最大似然算法进行学习，得到模块化贝叶斯网络模型；步骤S4、实时采集各观测点的传感器数据，利用推理引擎更新模块化贝叶斯网络模型各节点的后验概率，根据各节点的后验概率判断相应的节点是否发生故障。优选地，当所述特征变量库D1增加新的特征变量时，若新增变量映射到所述特征变量库D2，则依次更新所述关联矩阵、观测模块、最优图及贝叶斯网络模型；同理，当所述故障模式变量增加新的变量时，照样依次更新所述关联矩阵、观测模块、最优图及贝叶斯网络模型。优选地，所述步骤S1具体包括：S11：根据牵引系统的传感器分布，确定观测点的个数N；采集牵引系统历史正常状态数据及故障状态数据组成初始数据集，包括由恒定采样率为fs的电流传感器采集的电流瞬时值In,i＝1,2,...,N；其中，将三相电流传感器作为一个观测点，该观测点包括a、b和c相定子电流瞬时值In_a、In_b和In_c；S12：利用特征提取方法提取观测点In的第m个特征变量Fnm,n＝1,2,...,N,m＝1,2,...,M，建立特征变量库D1；S13：利用降维算法将特征变量库D1里的特征映射到低维的特征变量库D2，计算过程如下：对第n个观测点的特征Fnm进行去均值处理：式中，μ为均值计算函数，Fnm是D1中观测点In的第m个特征；计算主成分：式中，A是特征变量从D1映射到D2的变换矩阵，是D1中观测点In第m个特征对应的主成分；计算并提取D1中累计贡献率大于等于σ的前J个主成分作为D2的特征变量，计算公式为：式中，λnm是的特征值，σ为维度选取的阈值，J为降维之后的维数，Gn为In前J个主成分的累计贡献率，是低维特征变量库D2中In的第j个特征变量。优选地，在所述步骤S2的关联矩阵Cnk中，若第k个降维后的特征变量或故障模式变量对第n个观测点状态变量有影响，则对应的关联矩阵元素值为1，若无影响，则为0；k的取值范围包括所有和Mz；其中，将Xtrain划分成N个观测模块的分布情况为：优选地，所述步骤S3具体包括以下步骤：S31、利用贝叶斯信息标准评分函数对各个观测模块的所有图进行评分，然后计算出评分最高的最优图G*n；S32、优化G*n，将搜索空间限制在通过添加一条边，删除一条边或反转现有边来获得的相邻图上，利用贪心搜索算法确定是否接受相邻图；若是，则更新G*n为相邻图；若否，则继续搜索，直至搜索完所有相邻图为止；S33、获得最终的G*n后，将所有观测模块的G*n基于关联矩阵Cnk组合成完整的贝叶斯网络结构，基于观测模块得到的G*n包括以下网络节点：与所有故障模式Mz一一对应的故障节点；与所有特征一一对应的特征节点；与所有观测点一一对应的观测点状态节点；S34、利用最大似然估计方法计算贝叶斯网络结构中各节点间的条件概率参数，所述条件概率参数包括概率密度函数、对数似然函数和条件概率参数。优选地，所述步骤S4具体包括：S41、实时采集各观测点的传感器数据，计算相应的特征变量，得到输入证据Xtest；S42、已知贝叶斯网络节点为Ai＝{A1,A2,...,AL}，将得到的输入证据Xtest输入到贝叶斯网络模型中，利用信念传播推理引擎进行信度推理，计算传播到第i个节点Ai的信息，并根据该信息计算节点的置信度；完成贝叶斯网络中各节点的置信度更新；S43、根据S42得到的置信度，即各个节点的后验概率，设置检测阈值δ判断该节点是否发生故障；当故障节点的后验概率大于或等于δ，则判断该节点发生故障；当故障节点的后验概率小于δ，则判断该节点未发生故障。作为一个总的技术构思，本专利技术还提供了一种基于观测模块的贝叶斯网络轨道客车牵引系统故障诊断系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述任一方法的步骤。本专利技术具有以下有益效果：本专利技术考虑了轨道客车牵引系统多个观测点的信息，通过模块化的思想保障了故障诊断的效率，模块化的空间多变量信息融合的方法，能够快速、准确地诊断出牵引传动系统各个位置的故障，且基于观测模块的各节点(包括观测状态节点、特征变量节点和故障模式节点)具有良好的可扩展性和可变性，具有本文档来自技高网...

【技术保护点】
1.一种基于观测模块的贝叶斯网络轨道客车牵引系统故障诊断方法，其特征在于，包括：/n步骤S1、选取观测点，在各所述观测点部署传感器采集各观测点包含正常状态数据和故障状态的初始数据集；对各观测点单独从相对应的初始数据集中提取时频域特征建立特征变量库D1，并对其进行降维得到特征变量库D2；/n步骤S2、将各观测点的状态变量、降维后的特征变量和故障模式变量作为贝叶斯网络模型的输入，建立输入的各特征变量和故障模式变量与观测点状态变量之间的关联矩阵，根据所述关联矩阵将所有输入划分成与观测点数量相等数量的观测模块；/n步骤S3、利用结构学习算法学习每个观测模块的最优图，并将各个观测点的最优图合并成一个完整的网络；网络各节点之间的信度条件概率参数采用最大似然算法进行学习，得到模块化贝叶斯网络模型；/n步骤S4、实时采集各观测点的传感器数据，利用推理引擎更新模块化贝叶斯网络模型各节点的后验概率，根据各节点的后验概率判断相应的节点是否发生故障。/n

【技术特征摘要】
1.一种基于观测模块的贝叶斯网络轨道客车牵引系统故障诊断方法，其特征在于，包括：
步骤S1、选取观测点，在各所述观测点部署传感器采集各观测点包含正常状态数据和故障状态的初始数据集；对各观测点单独从相对应的初始数据集中提取时频域特征建立特征变量库D1，并对其进行降维得到特征变量库D2；
步骤S2、将各观测点的状态变量、降维后的特征变量和故障模式变量作为贝叶斯网络模型的输入，建立输入的各特征变量和故障模式变量与观测点状态变量之间的关联矩阵，根据所述关联矩阵将所有输入划分成与观测点数量相等数量的观测模块；
步骤S3、利用结构学习算法学习每个观测模块的最优图，并将各个观测点的最优图合并成一个完整的网络；网络各节点之间的信度条件概率参数采用最大似然算法进行学习，得到模块化贝叶斯网络模型；
步骤S4、实时采集各观测点的传感器数据，利用推理引擎更新模块化贝叶斯网络模型各节点的后验概率，根据各节点的后验概率判断相应的节点是否发生故障。


2.根据权利要求1所述的方法，其特征在于，当所述特征变量库D1增加新的特征变量时，若新增变量映射到所述特征变量库D2，则依次更新所述关联矩阵、观测模块、最优图及贝叶斯网络模型；或者
当所述故障模式变量增加新的变量时，则依次更新所述关联矩阵、观测模块、最优图及贝叶斯网络模型。


3.根据权利要求1或2所述的方法，其特征在于，所述步骤S1具体包括：
S11：根据牵引系统的传感器分布，确定观测点的个数N；采集牵引系统历史正常状态数据及故障状态数据组成初始数据集，包括由恒定采样率为fs的电流传感器采集的电流瞬时值In,i＝1,2,...,N；其中，将三相电流传感器作为一个观测点，该观测点包括a、b和c相定子电流瞬时值In_a、In_b和In_c；
S12：利用特征提取方法提取观测点In的第m个特征变量Fnm,n＝1,2,...,N,m＝1,2,...,M，建立特征变量库D1；
S13：利用降维算法将特征变量库D1里的特征映射到低维的特征变量库D2，计算过程如下：
对第n个观测点的特征Fnm进行去均值处理：



式中，μ为均值计算函数，Fnm是D1中观测点In的第m个特征；
计算主成分：



式中，A是特征变量从D1映射到D2的变换矩阵，是D1中观测点In第m个特征对应的主成分；
计算并提取D1中累计贡献率大于等于σ的前J个主成分作为D2的特征变量，计算公式为：





...

【专利技术属性】
技术研发人员：陈志文，陈文英，彭涛，樊欣宇，阳春华，伍晓赞，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南;43