技术领域
[0001] 本公开属于机械运行
稳定性在线监测领域，尤其涉及一种基于图论的旋转机械振动信号监测方法及系统。
背景技术
[0002] 本部分的陈述仅仅是提供了与本公开相关的背景技术信息，不必然构成在先技术。
[0003] 机械系统运行稳定性在线监测是确保机械系统正常运行，
预防故障和实现机械系统自动化的关键技术。时频分析方法因其能同时获取监测信号的时域信号和频域信息而得到广泛应用，其中现有的方法大多通过提取时
频谱特征(如谱峰值，谱峭度等)和训练分类器来实现机械系统的运行稳定性监测。
[0004] 而在实际应用中，提取何种特征往往依赖于先验知识，同时需要事先收集大量训练样本来保证监测效果。
发明内容
[0005] 根据本公开的一个或多个
实施例的一个方面，提供一种基于图论的旋转机械振动信号监测方法，其无需先验知识、特征提取以及训练分类器，则可进行在线无监督的旋转机械系统运行状态监测。
[0006] 本公开的一种基于图论的旋转机械振动信号监测方法，包括：
[0007] 采集旋转机械系统运行过程中的振动信号；
[0010] 构建每个窗口内功率谱相对应的图结构；
[0011] 将当前
时间窗口所对应的图结构与自定义中值图之间的权重距离，作为异常度指标；其中，权重距离为两个图结构所有相对应的加权边的权重差与相应加权边的权重最大值的比值的累加和；
[0012] 对异常度指标进行假设检验，判断机械系统运行过程是否出现异常。
[0013] 在一个或多个实施例中，构建每个窗口内功率谱相对应的图结构的具体步骤，包括：
[0014] 将每个
频率采样点作为图结构
节点，所有频率采样点两两之间的连线作为图结构的加权边；
[0015] 计算加权边两端节点之间的功率谱幅值差作为加权边的权重。
[0016] 在一个或多个实施例中，计算机械正常运行时所对应的正常图结构集合中，每一个图结构元素与剩余所有图结构元素的权重距离并求和，并且将集合中与剩余所有图结构的权重距离之和最小的图结构元素定义为中值图。
[0017] 在一个或多个实施例中，依据3σ原则对异常度指标进行假设检验。
[0018] 在一个或多个实施例中，若机械系统运行发生异常则进行报警。
[0019] 在一个或多个实施例中，若机械正常稳定运行，则更新正常图结构集合并进行下一次窗口内振动信号的监测。
[0020] 根据本公开的一个或多个实施例的另一个方面，提供一种基于图论的旋转机械振动信号监测系统，其无需先验知识、特征提取以及训练分类器，则可进行在线无监督的旋转机械系统运行状态监测。
[0021] 本公开的一种基于图论的旋转机械振动信号监测系统，包括：
[0022]
加速度
传感器，其被配置为采集旋转机械系统运行过程中的振动信号并传送至
控制器；
[0023] 控制器，所述控制器包括
存储器和处理器，所述存储器存储有
计算机程序，所述程序被处理器执行时能够实现以下步骤：
[0024] 接收振动信号；
[0025] 加窗处理振动信号；
[0026] 计算每个窗口内振动信号的功率谱；
[0027] 构建每个窗口内功率谱相对应的图结构；
[0028] 将当前时间窗口所对应的图结构与自定义中值图之间的权重距离，作为异常度指标；其中，权重距离为两个图结构所有相对应的加权边的权重差与相应加权边的权重最大值的比值的累加和；
[0029] 对异常度指标进行假设检验，判断机械系统运行过程是否出现异常。
[0030] 在一个或多个实施例中，所述处理器，还被配置为：
[0031] 将每个频率采样点作为图结构节点，所有频率采样点两两之间的连线作为图结构的加权边；
[0032] 计算加权边两端节点之间的功率谱幅值差作为加权边的权重。
[0033] 在一个或多个实施例中，所述处理器，还被配置为：
[0034] 计算机械正常运行时所对应的正常图结构集合中，每一个图结构元素与剩余所有图结构元素的权重距离并求和，并且将集合中与剩余所有图结构的权重距离之和最小的图结构元素定义为中值图。
[0035] 在一个或多个实施例中，所述处理器，还被配置为：
[0036] 依据3σ原则对异常度指标进行假设检验。
[0037] 在一个或多个实施例中，所述处理器，还被配置为：
[0038] 若机械系统运行发生异常，则输出报警信号。
[0039] 在一个或多个实施例中，所述处理器，还被配置为：
[0040] 若机械正常稳定运行，则更新正常图结构集合并进行下一次窗口内振动信号的监测。
[0041] 本公开的有益效果是：
[0042] (1)本公开通过当前时间窗口所对应的图结构与中值图之间的权重距离，作为异常度指标；对异常度指标进行假设检验，判断机械系统运行过程是否出现异常，通过对所监测振动信号进行在线分析，精确检测和识别机械系统状态发生异常的变化点位。
[0043] (2)本公开无需先验知识、无需进行特征提取、无需训练分类器，实现了在线无监督的旋转机械系统运行状态监测的目的。
附图说明
[0044] 构成本公开的一部分的
说明书附图用来提供对本公开的进一步理解，本公开的示意性实施例及其说明用于解释本公开，并不构成对本公开的不当限定。
[0045] 图1是本公开的一种基于图论的旋转机械振动信号监测方法
流程图。
[0046] 图2(a)是振动信号的功率谱。
[0047] 图2(b)是振动信号对应的图结构。
[0048] 图2(c)是振动信号对应的图结构的邻接矩阵表示。
[0049] 图3是原始振动信号。
[0050] 图4是本公开的一种基于图论的旋转机械振动信号监测结果。
具体实施方式
[0051] 应该指出，以下详细说明都是例示性的，旨在对本公开提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常理解的相同含义。
[0052] 需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本公开的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
[0053] 图1是本公开的一种基于图论的旋转机械振动信号监测方法流程图。
[0054] 如图1所示，本实施例的一种基于图论的旋转机械振动信号监测方法，包括：
[0055] (1)采集旋转机械系统运行过程中的振动信号。
[0056] 具体地，可通过加速度传感器来采集旋转机械系统运行过程中的振动信号。
[0057] 其中，旋转机械系统运行过程中的振动信号为正常振动信号或异常振动信号；
[0058] 正常振动信号代表旋转机械系统运行正常；
[0059] 异常振动信号代表旋转机械系统运行异常。
[0060] (2)加窗处理所采集的振动信号。
[0061] 针对所采集的振动信号进行实时加窗处理，一般选用矩形窗或汉宁窗。
[0062] (3)计算每个窗口内振动信号的功率谱。
[0063] 针对窗口内振动信号进行基于周期图谱的功率谱估计。
[0064] 如图2(a)所示，将所提取的功率谱记为Pm，其中m为时间。
[0065] 值得注意的是Pm包含n＝T*fs/2个频率成分，其中fs为
采样频率，T为窗函数长度。
[0066] (4)构建每个窗口内功率谱相对应的图结构。
[0067] 在具体实施中，构建每个窗口内功率谱相对应的图结构的具体步骤，包括：
[0068] (4.1)将每个频率采样点作为图结构节点，将所有频率采样点(共n个)两两之间进行连线，作为图结构的加权边；
[0069] (4.2)计算加权边两端节点之间的功率谱幅值差作为加权边的权重di,j，其中i,j为节点编号，如图2(b)所示。
[0070] 将上述构造的图模型数学表达为邻接矩阵。具体地，将权重di,j放在矩阵中第i行，第j列，进而将一个图结构表达为一个邻接矩阵，如图2(c)所示。
[0071] (5)将当前时间窗口所对应的图结构与自定义中值图之间的权重距离，作为异常度指标；其中，权重距离为两个图结构所有相对应的加权边的权重差与相应加权边的权重最大值的比值的累加和。
[0072] 其中，计算机械正常运行时所对应的正常图结构集合中，每一个图结构元素与剩余所有图结构元素的权重距离并求和，并且将集合中与剩余所有图结构的权重距离之和最小的图结构元素定义为中值图。
[0073] 具体地，将当前监测时刻前的z个正常图结构看作正常图结构集合N＝[G1,G2,...,Gz],并计算该集合的中值图 (例如：z取5)，用其表示机械系统正常运行状态，中值图计算公式为，
[0074]
[0075] 其中：M(G,Gk)为计算图G与图Gk间权重距离。
[0076] 权重距离计算公式如下：
[0077]
[0078] 其中，
[0079]
[0080] 上式中di,j为图G中节点i与节点j之间加权边的权重，d'i,j为图Gk中节点i与节点j之间加权边的权重。
[0081] (6)将当前时间窗口所对应的图结构与中值图之间的权重距离，作为异常度指标。
[0082] 将当前时刻所得图结构与中值图 之间计算权重距离。
[0083] (7)对异常度指标进行假设检验，判断机械系统运行过程是否出现异常。
[0084] 具体地，根据3σ准则对所得加权边距离D进行如下假设检验，
[0085] H0:无异常:D∈S,
[0086] HA:出现异常:
[0087] 其中S＝[μ-3σ,μ+3σ]为置信区间，μ为前h个所测距离的均值，σ为前h个所测距离的方差，h值推荐取30。
[0088] 原始振动信号，如图3所示；采用本实施例的基于图论的旋转机械振动信号监测的监测结果，如图4所示。若当前所测距离D不属于置信区间S则认为系统出现异常并进行报警，若当前所测距离D属于置信区间S则认为系统运行正常并更新正常图结构集合。
[0089] 本实施例通过当前时间窗口所对应的图结构与中值图之间的权重距离，作为异常度指标；对异常度指标进行假设检验，判断机械系统运行过程是否出现异常，通过对所监测振动信号进行在线分析，精确检测和识别机械系统状态发生异常的变化点位。
[0090] 本实施例无需先验知识、无需进行特征提取、无需训练分类器，实现了在线无监督的旋转机械系统运行状态监测的目的。
[0091] 本公开还提供了一种基于图论的旋转机械振动信号监测系统，包括：
[0092] (1)加速度传感器，其被配置为采集旋转机械系统运行过程中的振动信号并传送至控制器。
[0093] (2)控制器，所述控制器包括存储器和处理器，所述存储器存储有计算机程序，所述程序被处理器执行时能够实现以下步骤：
[0094] (a)接收振动信号。
[0095] 其中，旋转机械系统运行过程中的振动信号为正常振动信号或异常振动信号；
[0096] 正常振动信号代表旋转机械系统运行正常；
[0097] 异常振动信号代表旋转机械系统运行异常。
[0098] (b)加窗处理振动信号。
[0099] 一般选用矩形窗或汉宁窗。
[0100] (c)计算每个窗口内振动信号的功率谱。
[0101] 针对窗口内振动信号进行基于周期图谱的功率谱估计。
[0102] 如图2(a)所示，将所提取的功率谱记为Pm，其中m为时间。
[0103] 值得注意的是Pm包含n＝T*fs/2个频率成分，其中fs为采样频率，T为窗函数长度。
[0104] (d)构建每个窗口内功率谱相对应的图结构。
[0105] 在具体实施中，构建每个窗口内功率谱相对应的图结构的具体步骤，包括：
[0106] (d1)将每个频率采样点作为图结构节点，将所有频率采样点(共n个)两两之间进行连线，作为图结构的加权边；
[0107] (d2)计算加权边两端节点之间的功率谱幅值差作为加权边的权重di,j，其中i,j为节点编号，如图2(b)所示。
[0108] 将上述构造的图模型数学表达为邻接矩阵。具体地，将权重di,j放在矩阵中第i行，第j列，进而将一个图结构表达为一个邻接矩阵，如图2(c)所示。
[0109] (e)将当前时间窗口所对应的图结构与自定义中值图之间的权重距离，作为异常度指标；其中，权重距离为两个图结构所有相对应的加权边的权重差与相应加权边的权重最大值的比值的累加和。
[0110] 其中，计算机械正常运行时所对应的正常图结构集合中，每一个图结构元素与剩余所有图结构元素的权重距离并求和，并且将集合中与剩余所有图结构的权重距离之和最小的图结构元素定义为中值图。
[0111] 具体地，将当前监测时刻前的z个正常图结构看作正常图结构集合N＝[G1,G2,...,Gz],并计算该集合的中值图 (例如：z取5)，用其表示机械系统正常运行状态，中值图计算公式为，
[0112]
[0113] 其中：M(G,Gk)为计算图G与图Gk间权重距离。
[0114] 权重距离计算公式如下：
[0115]
[0116] 其中，
[0117]
[0118] 上式中di,j为图G中节点i与节点j之间加权边的权重，d'i,j为图Gk中节点i与节点j之间加权边的权重。
[0119] (f)将当前时间窗口所对应的图结构与中值图之间的权重距离，作为异常度指标。
[0120] 将当前时刻所得图结构与中值图 之间计算权重距离。
[0121] (g)对异常度指标进行假设检验，判断机械系统运行过程是否出现异常。
[0122] 具体地，根据3σ准则对所得加权边距离D进行如下假设检验，
[0123] H0:无异常:D∈S,
[0124] HA:出现异常:
[0125] 其中S＝[μ-3σ,μ+3σ]为置信区间，μ为前h个所测距离的均值，σ为前h个所测距离的方差，h值推荐取30。
[0126] 原始振动信号，如图3所示；采用本实施例的基于图论的旋转机械振动信号监测的监测结果，如图4所示。
[0127] 在一个或多个实施例中，所述处理器，还被配置为：
[0128] 若机械系统运行发生异常，则输出报警信号。
[0129] 在一个或多个实施例中，所述处理器，还被配置为：
[0130] 若机械正常稳定运行，则更新正常图结构集合并进行下一次窗口内振动信号的监测。
[0131] 本实施例通过当前时间窗口所对应的图结构与中值图之间的权重距离，作为异常度指标；对异常度指标进行假设检验，判断机械系统运行过程是否出现异常，通过对所监测振动信号进行在线分析，精确检测和识别机械系统状态发生异常的变化点位。
[0132] 本实施例无需先验知识、无需进行特征提取、无需训练分类器，实现了在线无监督的旋转机械系统运行状态监测的目的。
[0133] 本领域内的技术人员应明白，本公开的实施例可提供为方法、系统、或计算机程序产品。因此，本公开可采用
硬件实施例、
软件实施例、或结合软件和硬件方面的实施例的形式。而且，本公开可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
[0134] 本公开是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方
框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程
数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中
指定的功能的装置。
[0135] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0136] 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0137] 本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random AccessMemory，RAM)等。
[0138] 上述虽然结合附图对本公开的具体实施方式进行了描述，但并非对本公开保护范围的限制，所属领域技术人员应该明白，在本公开的技术方案的
基础上，本领域技术人员不需要付出创造性劳动即可做出的各种
修改或
变形仍在本公开的保护范围以内。
