首页/文章/ 详情

Information Fusion|重庆大学攻克图冗余传播难题,提出RUL预测新框架

1小时前浏览1

    本期聚焦机器剩余使用寿命预测新进展。剩余使用寿命预测长期受图冗余信息不可控传播与细粒度子图结构忽视两大技术瓶颈制约,全局视角表征不足、子图特征融合效率低更是领域普遍短板。重庆大学秦毅团队近日提出基于信息瓶颈的层次子图神经网络(hierarchical subgraph neural network with information bottleneck, HSNN-IB)框架,为多源传感器信息融合驱动的RUL预测提供了全新技术路径。

    该框架采用混合度量图划分算法,融合节点连接强度与确定性自动识别中心节点,构建具有代表性的子图;通过层次门控子图特征更新策略,先执行局部更新再选择性跨子图整合信息,有效提升节点表示质量;配套多子图融合信息瓶颈模块,压缩并整合子图融合表示,抑制冗余信息并保留任务相关特征。实验证实,该方法在公开基准数据集上优于当前最优方法,并在实际风电机组监测数据上验证了其实际适用性。

    论文基本信息

    论文题目:

    Hierarchical subgraph neural network with information bottleneck for machine remaining useful life prediction  

    论文期刊:Information Fusion

    论文日期:2026年

    论文链接:

    https://doi.org/10.1016/j.inffus.2025.103953

    作者:Linjie Zheng (a), Junyu Qi (b), Junyu Qi (a), Yi Qin*(a)

    机构:

    a. State Key Laboratory of Mechanical Transmission for Advanced Equipment, Chongqing University, Chongqing 400044, China; 

    b. Reutlingen Research Institute (RRI), Reutlingen 72762, Germany.

    通讯作者简介:秦毅,重庆大学教授,于重庆大学机械电子工程系获工学博士学位。长期致力于智能故障诊断与装备健康管理、动力学建模与数字孪生、智能无人系统以及智能制造等领域的创新研究。秦毅教授在国内外权威学术期刊及高水平会议上发表论文230余篇,研究成果涵盖《IEEE Transactions on Industrial Informatics》《Mechanical Systems and Signal Processing》《Reliability Engineering and System Safety》等知名期刊,学术影响力广泛。他长期担任国家重点研发计划、长江学者奖励计划、国家自然科学基金等多项国家级与省部级科研项目评审专家,深度参与国际学术交流与期刊建设。在学术组织方面,秦毅教授是中国机械工程学会、中国振动工程学会的高级会员,并担任IEEE可靠性学会重庆分会副主席,在智能故障诊断与动态测试技术领域积累了深厚的学术声望。(来源:https://faculty.cqu.edu.cn/YiQin)

    目录

    1 摘要

    2 引言

    3 所提方法

    3.1 图构建

    3.2 基于混合度量的图划分算法

    3.3 层次门控子图特征更新策略

    3.4 多子图融合信息瓶颈模块

    4 实验

    4.1 评估指标

    4.2 案例研究一:CMAPSS数据集

    4.3 案例研究二:SCADA数据集

    4.4 模型分析

    5 结论

    1 摘要

    图神经网络在多源传感器信息融合驱动的剩余使用寿命(RUL)预测中日益受到关注。然而,现有方法大多依赖整合整个图的结构特征与节点表示以获得全局视角,而忽略了更细粒度的子图结构。这往往导致图中冗余信息的不可控传播,限制了预测性能。

    为解决该问题,本文提出了一种基于信息瓶颈的层次子图神经网络(HSNN-IB)用于RUL预测。HSNN-IB采用所提出的基于混合度量的图划分算法,该算法融合节点连接的强度与确定性来自动识别中心节点,从而构建具有代表性的子图。随后,提出了一种层次门控子图特征更新策略,该策略先执行局部更新,再选择性地跨子图整合信息,从而有效提升节点表示质量。在此基础上,开发了一种多子图融合信息瓶颈模块,用于压缩并整合子图融合表示,有效抑制冗余信息并保留任务相关特征。最后,在公开基准数据集上的对比实验证明了所提方法优于当前最优方法,同时在实际风电机组监测数据上的实验验证了其实际适用性。

    关键词:剩余使用寿命预测;图神经网络;图划分;信息瓶颈;多传感器信息融合

    2 引言

    机械设备在长期运行过程中会逐渐出现性能退化。故障发生后的被动维护策略不仅会导致非计划停机与生产损失,还可能引发连锁安全事故,对人员生命与资产安全构成严重威胁。在此背景下,剩余使用寿命(RUL)预测技术已成为机械设备高效维护的关键支柱。

    RUL预测方法可大致分为两类:基于模型的方法与数据驱动的方法。基于模型的方法依赖物理建模,通常需要对设备运行机理有深入理解。因此,其在应用于复杂系统时面临重大挑战。相比之下,数据驱动方法利用大量监测数据,通过统计分析、机器学习或深度学习技术构建预测模型,且不需要精确的物理先验知识,因而展现出更强的适应性。特别是深度学习方法——如卷积神经网络(CNN)、长短期记忆(LSTM)网络和门控循环单元(GRU)网络——因其在特征提取与时序建模方面的优势,在RUL预测中受到广泛关注。

    然而,许多现有的基于深度学习的RUL预测方法主要关注挖掘传感器数据中的时序依赖关系,而相对忽略了不同传感器之间的空间相关性。实际上,传感器之间的空间依赖关系在表征设备局部结构状态与多源信息交互方面起着至关重要的作用。忽视这一特性可能导致模型无法充分捕捉关键退化模式,从而影响预测精度与鲁棒性。因此,挖掘空间依赖信息对于构建更全面、更精细的RUL预测模型具有重要意义。

    近年来,越来越多的研究开始采用图神经网络(GNN)来建模传感器之间的空间依赖关系。在GNN中,多传感器系统被表示为图结构,其中节点对应单个传感器,边描述它们之间的空间或功能关系。因此,GNN能够有效捕捉传感器之间的结构化交互信息。与传统方法相比,GNN利用消息传递机制聚合邻域节点信息,实现多源异构传感器数据的深度特征提取与融合。该方法显著提升了设备退化建模的精度,为RUL预测提供了新的技术途径。

    尽管GNN已成功应用于建模多个传感器之间的空间依赖关系,但现有的基于GNN的RUL预测方法仍主要依赖图级信息传播进行特征更新。然而,随着传播范围扩大,目标节点不断接收来自远处结构区域的特征信号,其中不可避免地包含大量任务无关甚至干扰性的噪声特征。这一现象最终导致节点表示逐渐模糊或语义偏移。此类局限性不仅降低了模型的判别能力,还严重损害了其在复杂下游任务中的泛化性能。

    为解决该问题,本文提出了一种基于信息瓶颈的层次子图神经网络(hierarchical subgraph neural network with information bottleneck, HSNN-IB)用于RUL预测。与依赖全图级信息传播策略的传统图神经网络不同,HSNN-IB设计了一种基于混合度量的图划分算法(hybrid metric-based graph partitioning algorithm, HMGPA),用于将原始图划分为多个子图,该算法能够有效定位信息传播范围。具体而言,本文提出了两个子图中心节点选择准则,即中心节点应具有强连接强度与高连接确定性。在此基础上,构建混合度量对节点中心性进行排序并选择中心节点,从而将整个图划分为若干子图,每个子图聚焦于与中心节点直接相关的局部结构。为实现高效特征学习,提出了一种层次门控子图特征更新(hierarchical gated subgraph feature update, HGSFU)策略,该策略包含两个特征更新阶段。在子图内更新阶段,节点特征仅基于原始图中的局部连接进行更新,使模型能够首先学习结构可靠的局部语义表示。在子图间信息交互阶段,动态控制不同子图之间的信息流动强度,确保仅传播任务相关表示,同时抑制来自远处结构的噪声扩散。通过这种层次化更新设计,模型不仅增强了节点表示对局部结构的敏感性,还隔离了无关区域的影响,并构建了层次化特征。此外,设计了一种多子图融合信息瓶颈模块(multi-subgraph fusion information bottleneck module, MFIBM),用于压缩和优化来自多个子图的融合表示,迫使网络仅保留对下游任务最关键的特征。本文的主要贡献总结如下:

    • 提出了一种新颖的基于混合度量的图划分算法(HMGPA),该算法联合考虑节点连接强度与连接确定性,构建混合度量以量化和排序节点中心性。基于该排序,算法智能选择最优子图中心节点,实现对图的精确划分。

    • 提出了一种层次门控子图特征更新(HGSFU)策略,以增强局部结构学习并有效调控全局信息流。

    • 首次设计了一种多子图融合信息瓶颈模块(MFIBM),用于压缩和整合来自多个子图的表示,显著提升模型的RUL预测性能。

    • 开发了一种基于信息瓶颈的层次子图神经网络(HSNN-IB)的新RUL预测方法。其优越性与实际适用性通过公开基准数据集和实际风电机组的多源监测数据得到验证。

    3 所提方法

    所提出的基于HSNN-IB的RUL预测方法的整体流程如图1所示。其主要部分,包括图构建、基于混合度量的图划分算法以及层次门控子图特征更新策略,将在以下小节中详细阐述。

         

    图1 所提出的用于RUL预测的HSNN-IB整体框架

    3.1 图构建

    针对多源传感器数据,不同传感器通道通常具有不同的物理含义,且它们之间存在复杂且潜在的依赖关系。图结构为捕捉多维数据中的交互提供了一种有效的建模途径,其中节点表示实体(如传感器通道),边描述实体之间的关系或相似性。在本研究中,原始多维传感器数据被转换为图表示,以建模多个传感器源之间的潜在关系。图可表示为           ,其中            表示节点集 合,           表示边权重集 合。           表示原始图中节点的总数。令            表示经过最小-最大归一化操作后的多源传感器数据,节点            与            之间的初始边权重            计算如下:

             

    其中            表示            与            之间的余弦相似度。鉴于归一化数据的余弦相似度落在            范围内,选取该区间的中点值            作为阈值。该值能够保留关键边,同时有效滤除弱相关和噪声边。该值的有效性已通过后续敏感性分析实验得到验证。采用 LSTM 网络提取            的初始特征表示           :

             

    其中            表示 LSTM 网络的参数。

    3.2 基于混合度量的图划分算法

    为了减少信息冗余并促进高效的层次化特征学习,在更新每一层特征之前,将原始图            划分为            个具有代表性的子图           。该划分过程有助于约束信息传播范围,从而提升局部特征表示的质量。这些子图的节点互斥,且每个子图中的最大节点数设为           。由于算法工作流程的独特设计,子图中的实际节点数可能小于           ,从而有效实现冗余节点的自动排除。因此,所有子图中的节点总数可能小于原始图的节点总数。           和            均为预定义参数。具体过程如下。

    令              表示原始图              中的一个节点。则其邻居节点集 合定义为:

               

                的计算公式如下:

               

               

               

    其中              表示节点              在第              层的特征表示,             表示节点              与              之间的相似度,             表示批归一化,             表示              卷积层,             表示 Sigmoid 激活函数。

    以上公式共同构成相似度计算网络。本文采用可学习的相似度计算网络,因为与固定相似度度量相比,它更适合处理在多个网络层中持续更新的节点特征。可学习的相似度函数能够自适应地建模节点之间的潜在动态关系,从而产生更准确的边权重。因此,下游预测任务的性能可以得到有效提升。

    在图划分阶段,为确保所选子图具有代表性,需要为每个子图选择中心节点。中心节点的选择准则定义如下:

    1) 节点应具有高连接强度,即它应尽可能与更多邻居节点直接相连,以确保其代表性。

    2) 节点与其邻居之间的连接应表现出高确定性,即连接权重的分布应具有低信息熵,反映信息传输的结构稳定性与可靠性。为量化第一个准则(即节点的连接强度),节点的归一化平均连接强度定义为:

               

    可以很容易地注意到,             反映了图结构中节点与其邻域之间的整体交互强度,即它可作为节点作为中心节点的连通覆盖能力的度量。为量化第二个准则,本文利用信息熵来刻画节点邻接中边权重分布的确定性。具体而言,对于给定节点             ,其邻接边权重分布的信息熵写为:

               

    其中              表示第              层节点              与              之间的归一化边权重,反映该边在节点              所有连接中的相对重要性。信息熵越小意味着边权重分布越集中,表明该节点的结构确定性越高。通过整合这两个指标,本文构建中心性评分函数如下:

               

    其中              是一个调节参数(默认设为0.5),用于平衡连接强度与连接确定性之间的权重。基于所有节点的中心性得分,按降序排列如下:

                 

    假设所选节点集 合记为               ,初始化时将其设为空集。第                个子图的构造过程如下,从未使用的节点中,选择当前排名最高的节点作为中心节点:

                 

    然后,获取该节点的所有邻居节点如下:

                 

    接下来,基于与中心节点相连的前                个最大边权重,从邻居节点中选择至多                个节点,并将其指定为该子图内中心节点的邻居节点:

                 

    其中                表示从剩余邻居节点中选出的前                个节点的集 合,按边权重                排序。由于某些节点对的边权重为零,这可能导致某些节点的剩余邻居不足(即其数量小于               )。随后,第                个子图的节点集构造如下:

                 

    随后,已使用节点集 合按下式更新:

                 

    最后,子图的边集按如下方式构成:

                 

    3.3 层次门控子图特征更新策略

                     

    图2展示了所提出的HGSFU策略的工作流程,该策略在两个连续阶段中运行:第一阶段完成子图内节点特征的更新,在此基础上,第二阶段实现跨子图的信息传播。

           

    图2 层次门控子图特征更新策略的工作流程

    3.3.1 子图内节点特征更新

    完成图划分后,每个节点的局部特征在子图内进行更新。局部邻域信息通常携带最直接且高置信度的结构信息。通过优先进行子图内的局部聚合,模型能够构建稳定且具有判别性的局部表示,从而为后续的子图间信息交换奠定基础。令              表示节点              在第              个子图第              层的特征表示。采用标准消息传递机制在每个子图内实现第              层的特征更新,其公式如下:

               

    其中,             表示子图内节点特征更新网络,该网络由两个全连接层组成,每个全连接层后接一个 Tanh 激活函数。

    3.3.2 门控跨子图信息传播机制

    为实现全局信息流动,有必要进一步融合与协调跨子图的信息。为此,本文提出了一种门控跨子图信息传播机制,该机制允许每个子图内的节点吸收来自其他子图中与其相关的特征,同时抑制无关干扰。首先,对于每一对跨子图节点              与             ,跨子图边权重              由下式计算:

                 

    其中,               表示子图间相似度计算网络。接下来,对来自子图                的所有节点进行加权求和,以获得中间聚合表示,如下所示:

                 

    为增强跨子图信息传播的可控性与适应性,引入两个门控因子                和               ,其定义如下:

                 

                 

    其中                表示特征拼接,              、              、               和                为可训练网络参数。               用于控制信息融合,而                用于抑制特征冗余。最终,节点                在第                层的特征表示                由下式更新:

                 

    完成层次子图特征更新后,所有子图被重新组合为完整图,并重新计算整个图的边权重,为下一层的图划分做准备。

    3.4. 多子图融合信息瓶颈模块

    在完成子图内局部特征更新与跨子图门控信息传播后,设计了一种MFIBM,以进一步整合全局信息并抑制冗余特征对预测精度的干扰。该模块旨在压缩并提取多个子图中与RUL预测任务高度相关的信息,从而获得面向下游应用的低冗余联合表示。

    给定子图               ,其子图表示                由下式构成:

                   

    在完成子图构建、子图内节点特征更新以及跨子图信息传播之后,所得表示                  包含了局部邻域中与退化模式高度相关的结构信息,以及不同子图之间的相关特征。这表明,在输入MFIBM之前,                 已经编码了从原始数据中提取的与RUL预测密切相关的核心信息。由于MFIBM的优化过程只能压缩和精炼已有的输入信息,而无法生成额外的任务相关内容,因此模型的预测能力取决于输入中是否存在与RUL相关的信息。若输入缺乏足够的RUL相关信息,模型将难以具备良好的预测性能。事实上,在多个数据集上获得的结果进一步验证了                  确实包含足以服务于RUL预测的信息。因此,                 可被视为一种“充分信息源表示”。

    为从多个子图表示中获得压缩且任务相关的联合表示,基于信息瓶颈原理,MFIBM的优化目标形式化如下:

                   

    其中超参数                  作为信息压缩与任务相关性之间的权衡系数。

    在RUL预测任务中,模型的核心目标是从传感器数据中提取能够反映设备退化趋势与健康状态的关键特征,即“任务相关特征”。若特征表示                  与目标变量                 (RUL)高度相关,则对                  的观测应能显著降低                  的不确定性。互信息                  是量化这种不确定性降低的有价值指标,形式上刻画了“已知                  所提供的关于                  的信息增益”。因此,最大化                  本质上鼓励模型增强与RUL预测最相关特征的表示。

    此外,从                  的定义来看,其等价于最大化条件对数似然                 ,即:

                   

    在RUL预测场景中,假设给定表示                  时,RUL标签                  服从均值为                 、方差为常数                  的高斯分布,则优化目标可表示为:

                   

    忽略与模型参数无关的常数项,上述表达式的最大化等价于最小化预测值与真实值之间的期望平方差(即均方误差(MSE)):

                   

    因此,在常见且合理的高斯观测假设下,公式中的预测项在优化目标上与RUL任务中常用的MSE损失一致。然而,若仅优化预测项,模型可能将输入中的所有信息编码到                  中,这会导致过拟合和泛化能力减弱。公式中的压缩项                  约束了                  的容量,鼓励模型优先选择能显著增强                 (即降低任务损失)的信息子集。这样,任务相关特征得以保留,而冗余信息被丢弃,从而提升泛化能力。在RUL预测背景下,这意味着模型必须主动识别并移除冗余输入信息,仅保留最关键和有用的特征。

    因此,MFIBM的优化目标在理论上与RUL预测任务是兼容的。由于                  的计算具有挑战性,基于已知形式的分布                  计算一个上界以实现可处理的优化,其推导如下:

                   

    该不等式由Kullback-Leibler(KL)散度的非负性引入。在实践中,采用标准高斯分布                  作为变分分布。对于                 ,假设其服从高斯分布,其参数由神经网络估计。计算过程如下:

                   
                 

                     

                     

    其中                   ,                   是一个单层全连接层。利用重参数化技巧,联合表示嵌入                    通过下式计算:

                     

    随后,预测的RUL标签由下式获得:

                     

    最后,所提出模型的优化目标表示为:

                     

    由于                    和                    均服从高斯分布,上述方程中的KL散度具有闭式解:

                     

    表1展示了所提出模型的伪代码,涵盖了数据输入、子图构建和特征更新等主要过程。

    表1 所提出模型的伪代码

                 

    4 实验

    4.1 评估指标

    为定量评估所提方法的性能,本文采用均方根误差(RMSE)和评分函数(SF)作为评估指标,其公式如下:

           

           

    RMSE衡量预测值与实际值之间的偏差,因此适用于评估整体预测精度。SF指标对提前预 测和滞后预测施加不同程度的惩罚,从而反映模型在不同预测阶段的性能。RMSE或SF越低,表示预测精度越高。两种指标的结合使用能够对模型预测性能提供更全面的评估。

    4.2 案例研究一:CMAPSS数据集

    4.2.1 CMAPSS数据集的介绍与预处理

    CMAPSS数据集由NASA开发,是PHM领域广泛使用的基准数据集。它模拟了涡扇发动机在不同飞行条件下的退化过程。该数据集包含四个子数据集(FD001–FD004),每个子数据集代表不同的运行条件和故障模式,详见表2。每个子集包括训练集和测试集。训练集包含来自多台发动机的完整全寿命运行数据,而测试集仅包含部分运行数据。每台发动机由21个传感器监测。由于并非所有传感器都携带 meaningful 的退化信息,本研究选择编号为2、3、4、7、8、9、11、12、13、14、15、17、20和21的传感器用于RUL预测。在每个子数据集中,所有传感器读数均归一化至        范围。在本研究中,监测数据采用长度为60、步长为1的滑动时间窗口分割为多个样本。

    表2 CMAPSS数据集的细节

             

    4.2.2. 对比方法与实验设置

    为有效验证所提方法的性能与优越性,将其实验结果与九种主流方法进行对比分析。这九种方法包括基于CNN的方法(MS-DCNN、AGCNN和STRUL)、基于LSTM的方法(DA-LSTM)、基于Transformer的方法(TATFA),以及基于GNN的方法(MGCAL-UQ、ASTHGNN、GCN-PR和KEST-GCN)。所提方法在PyTorch中实现,采用Adam优化器更新模型参数。学习率设为0.001,训练轮数为200,批次大小设为512。在四个子数据集中,              和              的值分别设为3和5。表3展示了所提模型在CMAPSS数据集上的网络架构与参数设置。

    表3 所提出的 HSNN-IB 在 CMAPSS 数据集上的网络架构

             

    以FD001子数据集为例,图3展示了使用所提方法的图划分结果。超参数              在四个子数据集上分别设为[2e-5, 1e-5, 1e-5, 1e-5]。

             
    图3 FD001 子数据集上的图划分结果          

    4.2.3 性能对比与分析

    多种方法在CMAPSS数据集上的实验结果总结于表4。如表所示,所提方法在所有子数据集上取得了七项最优指标和一项次优指标,整体性能显著优于其他基线方法。基于CNN的方法(如MS-DCNN)在复杂工况下表现不佳,因其无法有效捕捉信号的时序特征。例如,在FD004子数据集上,MS-DCNN的RMSE高达22.22,与其他方法存在显著性能差距。尽管DA-LSTM能够建模信号中的时序依赖关系,但缺乏捕捉多传感器间空间依赖的能力,限制了其相对于所提方法的性能。

    表4 各种方法在CMAPSS数据集上的实验结果

             

    与主流基于GNN的方法(如GCN-PR和KEST-GCN)相比,所提方法引入了子图分割策略以实现层次化特征更新。此外,它结合了信息瓶颈机制以有效融合子图特征,从而在抑制冗余信息的同时增强关键特征的表示。因此,所提方法在所有子数据集上均展现出优越的预测性能。图4展示了四个子数据集上预测RUL值与实际RUL值的对比。预测RUL与实际RUL的高度一致性进一步验证了所提方法在复杂工况下的有效性与可行性。

             

    图4 所有测试发动机上实际RUL与预测RUL的对比

    4.2.4 消融实验

    为评估各组件对所提模型整体性能的影响,本文通过设计三种变体对HSNN-IB进行消融研究:模型1:不含图划分策略的HSNN-IB。模型2:不含HGSFU策略和MFIBM的HSNN-IB。取而代之,采用常规信息聚合方法进行跨子图消息传递,并使用拼接生成多子图融合表示。模型3:不含MFIBM的HSNN-IB,同样采用拼接方法进行多子图表示融合。

    消融结果如图5所示,清晰表明HSNN-IB优于所有三种变体。模型2与模型3的对比表明,所提出的HGSFU策略能够实现更有效的子图级特征学习。模型3与HSNN-IB模型之间的性能差异表明,MFIBM有效压缩了冗余信息并提升了多子图融合表示的质量。模型1表现最弱,凸显了图划分策略的重要性。这些发现共同证实,所提模型的每个组件均对提升其整体性能有实质性贡献。

         

    图5 消融实验结果

    4.3 案例研究二:SCADA数据集

    4.3.1 SCADA数据集的介绍与预处理

    为验证所提方法在实际工程场景中的适用性,本文使用位于中国西北某风场的多台2 MW风电机组的真实运行数据评估其性能。风电机组的主要结构如图6所示。

         

    图6 实验风电机组的主要部件

    数据由SCADA系统以1分钟间隔采集,覆盖2018年1月1日至2019年12月29日期间。共安装26个传感器,部分传感器的功能与参数列于表5。

    表5 部分传感器的描述

         

    图7展示了特定风电机组在选定时间段内的运行数据。本研究选择编号为7、15和18的风电机组作为实验对象。这些机组在数据采集期间发生故障,并在停机维护后重新投入运行。为构建样本数据集,提取每台机组每天中午12:00记录的运行数据,并采用滑动窗口方法生成实验样本。

         

    图7 部分传感器监测数据

    4.3.2 实验设置

    选取7号和18号风电机组的运行数据作为训练集,15号风电机组的数据作为测试集。时间滑动窗口长度设为60,最大RUL上限设为55。这三台风电机组维护后采集的监测数据也将用于训练和测试,由这些数据划分的样本标签设为最大RUL。所提方法使用PyTorch机器学习框架实现。采用Adam优化器更新模型参数,学习率为0.001,批次大小为512,总训练轮数为600。模型超参数          、          和          分别设为4、7和5e-5。

    4.3.3 性能对比与分析

    为全面评估所提方法的有效性,与多个基准模型进行对比实验,包括BiLSTM、SA-LSTM、BiGRU、GCN和GGCN。如表6所示,所提方法在RMSE和SF两方面均显著优于基线方法,表明其在RUL预测中具有更高的精度。

    表6 SCADA数据集上各种方法的实验结果

         

    图8展示了某次实验中HSNN-IB分割出的四个子图。可以看出,节点1、3和15在此图划分过程中被丢弃,这证明了HSNN-IB自适应消除冗余节点的能力。

         

    图8 SCADA数据集上的图划分结果

    此外,图9展示了HSNN-IB与另外三种对比方法的RUL预测结果。可以看出,HSNN-IB预测的RUL与实际RUL更为接近,表明拟合效果更好。这些实验结果共同验证了HSNN-IB在RUL预测中的优越性能,并凸显了其在实际工业应用中实现更准确、更稳定RUL预测的潜力。

         

    图9 SCADA数据集上的RUL预测结果

    4.4 模型分析

    4.4.1 超参数          的敏感性分析

    在节点中心性评分函数中,超参数          用于平衡连接强度与连接确定性两个组成部分,在节点中心性评估中起关键作用。为验证          设置的合理性及模型对其的鲁棒性,在四个FD001–FD004子数据集和SCADA数据集上开展了参数敏感性分析实验,系统研究了不同          下的模型性能。结果如图10所示。

         

    图10 超参数            对模型性能的影响

    总体而言,          的变化对模型性能有一定影响。具体来说,两种评估指标在大多数情况下呈适度波动;然而,当          接近0.1或0.9时,模型在所有数据集上的预测性能均有一定程度下降。相比之下,当参数          在0.5左右时,模型表现相对更好。更具体地,当          取值在0.5–0.7范围内时,模型在FD002子数据集上取得优越的预测结果;当          在0.4–0.6范围内时,其在FD003子数据集上也表现出良好性能。

    这一趋势在四个CMAPSS子数据集和SCADA数据集上均一致观察到,表明当          在0.5左右时,模型表现出更好的稳定性与泛化能力。基于上述发现,本文默认将          设为0.5,以在“连接强度”与“连接确定性”之间实现合理平衡。此外,为提升特定任务中的性能,可在此基线值附近进行微调。

    4.4.2 超参数          的敏感性分析

    参数          作为MFIBM模型中信息压缩项的权重系数,其主要作用是约束潜变量从输入特征中保留的信息量,从而在特征压缩与信息保留之间实现权衡。为研究该超参数对模型整体性能的影响,在CMAPSS数据集和SCADA数据集上进行了系统的敏感性实验。具体实验结果如图11所示。

           

    图11 超参数 λ 对模型性能的影响

     从整体趋势来看,随着          值的不断增大,模型在多个数据集上的预测精度呈现波动的下降趋势。这一现象可能是因为当信息压缩项的权重设置过高时,一些对预测任务至关重要的有效特征也被过度压缩,导致关键判别信息丢失,从而引起模型性能下降。值得注意的是,在不同数据集上均可观察到类似规律,表明          对模型性能的影响具有一定程度的普遍性。基于实验结果与分析,为防止有用信息被过度压缩而导致性能损失,建议将          设置在相对较小的值范围内,以避免因过度压缩而造成性能下降。

    4.4.3 超参数          和          的敏感性分析

    超参数          控制图结构中划分的局部子图数量。若          值过小,模型性能可能因局部结构表达能力不足而受限;若          值过大,则可能引入冗余划分并增加计算开销。超参数          定义每个子图允许的最大节点数,从而控制子图规模。考虑到CMAPSS与SCADA数据集在传感器规模上的差异,针对这些数据集采用不同的          组合进行敏感性实验。对于CMAPSS数据集,使用          、          、          和          ;对于SCADA数据集,使用          、          、          和          。

    实验结果如图12所示,表明模型性能在不同          组合下存在一定波动。具体而言,CMAPSS数据集在          时获得最高预测性能,而SCADA数据集在          时达到最优。基于上述结果,我们建议根据具体数据集在合理范围内选择适中的子图数量和节点规模,以平衡模型的表达能力与计算效率。综上,          和          的适当设置有利于预测性能的提升。

         

    图12 超参数 K 和 M 对模型性能的影响

    4.4.4 初始边阈值的敏感性分析

    为考察公式中初始边阈值对图构建及模型性能的影响,在CMAPSS和SCADA数据集上进行了敏感性实验。该阈值直接决定初始图中保留的边数量及空间依赖关系的强度。阈值过低会引入大量弱相关或噪声连接,而阈值过高则可能导致图过于稀疏,损害有效的节点交互与特征传播。因此,有必要评估不同阈值下的模型性能。实验结果如图13所示。

         

    图13 始边阈值对模型性能的影响

    在FD001和FD003子集中,大多数节点相似度大于0.8。因此,当阈值设置为低于0.8时,所得图结构差异很小,模型预测性能基本保持不变。据此,重点关注这两个子集中阈值超过0.8的情况。如图13所示,当阈值为0.5或0.8时,模型表现几乎相同;然而,当阈值进一步增至0.85或0.9时,图变得过于稀疏,导致模型性能下降。

    对于其他数据集,阈值对模型性能具有类似影响。当阈值较小(如0.3)时,大量弱相关或噪声边被保留,导致预测性能不佳。当阈值增至0.7以上时,一些关键边被移除,造成性能下降,这在SCADA数据集上尤为明显。基于上述分析,公式中的初始边阈值最终设为0.5。该值有效保留了重要的空间依赖关系,同时避免了噪声边的干扰,从而确保了稳定且优越的预测性能。

    4.4.5 不同子图表示计算策略的对比分析

    为评估不同子图表示计算方法在本任务中的有效性,本文设计了若干实验,系统比较三种方法在所有CMAPSS子数据集和SCADA数据集上的性能,包括平均池化(Model4)、拼接操作(Model5)以及基于多头自注意力的加权求和(Model6)。需要指出的是,由于所提出的HMGPA方法生成的子图中节点数量不一致,Model5在构建子图表示时需要采用零填充以满足后续线性层的输入要求。实验结果如图14所示。

         

    图14 不同子图表示计算方法的比较

    从结果可以看出,Model4在所有五次实验中均取得最优或接近最优的综合性能,展现出强大的泛化能力与稳定性。相比之下,另外两种策略在不同方面表现出一定局限性:Model5仅在FD002子数据集上表现略好,但在大多数情况下性能较差,这可能是由于填充操作引入了噪声并增加了学习难度。Model6在FD001上取得了最低的RMSE,但在SCADA数据集上性能显著下降。这可能是由于其复杂权重学习机制存在稳定性不足和泛化能力较弱的问题。

    此外,表7从参数量、内存占用和计算时间方面对这三种方法进行了比较。可以看出,拼接操作显著增加了参数量,而加权求和操作则导致计算时间增加。相比之下,平均池化在不显著增加参数量或计算时间的情况下取得了优异性能。基于上述比较结果,本文最终选择平均池化作为子图表示的计算方法。

    表7 不同子图表示计算方法的计算复杂度与效率

         

    6 结论

    本文提出了一种新颖的HSNN-IB用于RUL预测。所提出的模型采用HMGPA自动提取具有代表性的子图结构,有效捕捉关键的局部拓扑信息。通过引入HGSFU策略,模型显著增强了对局部结构的敏感性,同时有效抑制了冗余信息在子图间的传播。此外,利用MFIBM压缩并优化多个子图的融合表示,从而仅保留对下游预测最任务相关的信息。在CMAPSS和SCADA数据集上进行的大量实验证明了所提方法的优越性与工程适用性。尽管HSNN-IB能够分割出具有代表性的子图,但仍无法实现子图结构的优化。因此,我们未来的研究方向之一是探索图结构优化方法。

    编辑:Jin

    校核:李正平、陈凯歌、Leo、曹希铭、赵学功、白亮、任超、海洋、Tina、陈宇航、Kira、肖鑫鑫、赵诚、张优
    该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除


    来源:故障诊断与python学习
    ACTMechanicalSystem振动复合材料旋转机械电子ANSA海洋声学理论电机材料数字孪生控制装配
    著作权归作者所有,欢迎分享,未经许可,不得转载
    首次发布时间:2026-11-06
    最近编辑:1小时前
    故障诊断与python学习
    硕士 签名征集中
    获赞 87粉丝 154文章 345课程 0
    点赞
    收藏
    作者推荐

    面向安全关键工程系统的可信智能体 AI 框架:可审计的无故障标签机械故障诊断(下)

    接上篇(故障诊断大模型 | 面向安全关键工程系统的可信智能体 AI 框架:可审计的无故障标签机械故障诊断(上)),论文提出了一种由感知、诊断与验证三个智能体协作的可信诊断框架。系统仅利用正常工况数据和轴承几何参数,通过调用物理分析工具,在无需故障类别标签的条件下开展诊断,并记录谐波、频率偏差与边带等关键证据。论文将可审计推理、诊断纠错与组件失效后的规则回退纳入统一流程,为智能体走向可靠的工业应用提供了值得关注的思路。下篇主要介绍实验和讨论,更多实验细节可查看论文附录链接。论文基本信息论文题目:A trustworthy agentic AI framework for auditable label-free machinery fault diagnosis in safety-critical engineering systems论文日期:2026(online)发表期刊:Advanced Engineering Informatics 77 (2027), 105205论文链接:https://doi.org/10.1016/j.aei.2026.105205作者:Yuntong Chen,Binhao Liu,Yingqi Li,Ziang Wang,Liping Ma,Jianyu Liu,Xitian Tian,Lijiang Huang(∗)作者单位:School of Mechanical Engineering, Northwestern Polytechnical University, Xi’an, Shaanxi 710072, China.目录摘要符号与缩略语1 引言2 相关工作2.1 轴承故障诊断的深度学习方法2.2 基于大语言模型的轴承故障诊断2.3 工程信息学中的多智能体大语言模型系统2.4 研究空白总结3 方法3.1 问题定义与框架概述3.1.1 问题定义3.1.2 框架架构3.2 基于自编码器的正常基线学习3.3 感知智能体:大语言模型自适应异常评估3.3.1 动机3.3.2 多指标融合3.3.3 两级决策架构3.4 诊断智能体:基于物理的工具链推理3.4.1 设计理念3.4.2 三层工具架构3.4.3 推理过程3.4.4 跨域泛化3.5 验证智能体:证据一致性与反思纠正3.5.1 动机3.5.2 三步证据一致性协议3.5.3 反思纠正循环3.6 多智能体编排流程3.6.1 流程执行3.6.2 设计属性3.6.3 计算成本3.7 框架的可信性属性4 实验与结果4.1 实验设置4.2 域内评估4.3 消融研究4.4 感知智能体评估4.5 向未见目标数据集的无故障标签迁移4.6 诊断推理分析4.6.1 系统性审计记录评估4.6.2 案例一:标准正确诊断4.6.3 案例二:边带消歧4.6.4 包络谱证据可视化4.6.5 错误分析4.6.6 频率匹配容差分析4.6.7 轴承参数扰动分析4.7 与基于大语言模型的诊断方法比较5 结论4. 实验与结果4.1. 实验设置为评估基于物理的多智能体诊断推理框架,我们在四个滚动轴承数据集上实验:CWRU [22] 为域内评估源域,Paderborn [47]、JNU、MFPT 为跨数据集无故障标签迁移的三个目标域。这些数据集在轴承几何、采集设备、故障生成机理及运行条件上显著不同,为域内准确率与跨域泛化提供严格测试环境。检测率对应可靠性工程中的检测概率(POD)。详细配置见表 2。实验环境:NVIDIA RTX 4090 GPU(24 GB)、AMD Ryzen R7 9800X3D CPU(5.6 GHz)、64 GB 内存。框架采用 PyTorch 2.1 实现。CWRU 数据集来自电机驱动机械系统,使用电火花加工(EDM)在内外圈引入三种尺寸的单点缺陷:0.007、0.014 和 0.021 英寸。使用驱动端加速度计的 12000 Hz 数据,覆盖 0–3 HP 四种电机负载,对应转速 1797、1772、1750、1730 RPM。每个 LLM 的测试量为 7 类 × 4 负载 × 5 个重复片段,共 140 个样本。Paderborn 提供更具挑战的场景,既有 EDM 和钻孔制造的人工缺陷,也有真实运行条件下加速寿命试验形成的故障。信号从 64000 Hz 重采样到 12000 Hz。四个健康轴承 K001–K004 仅用于自编码器训练;测试包含 23 个轴承、五类状态:健康 2 个、人工外圈故障 3 个、真实磨损外圈故障 7 个、人工内圈故障 3 个、真实磨损内圈故障 8 个。关键的是,训练或提示构建都不使用这两个数据集的任何带标签故障数据,自编码器仅学习正常信号。江南大学 JNU 数据集包含 600、800 和 1000 RPM 的变速工况,这是 CWRU 与 Paderborn 不具备的独特挑战。正常、外圈及滚珠故障使用 N205 轴承, ;内圈故障使用 NU205, 。振动采样率 50000 Hz,重采样至 12000 Hz。每种转速取一个正常文件,共三个文件训练自编码器;在正常、内圈、外圈三类的九个文件上测试,每文件三个重复片段,共 27 个样本。MFPT 数据集由机械故障预防技术学会提供,使用以 1500 RPM 运行的 NICE 轴承。正常信号采样率为 97656 Hz,故障信号为 48828 Hz,均重采样到 12000 Hz。三个正常文件用于训练;20 个测试文件包括正常 3 个、内圈 7 个、外圈 10 个,每文件三个重复片段,共 60 个样本。MFPT 在轴承诊断文献中广泛使用,故障类别定义明确。为考察模型通用性,我们评估三个先进 LLM 作为推理骨干:阿里巴巴 Qwen-Plus、OpenAI GPT-4.1 和 Anthropic Claude Sonnet 4.6。监督基线在 CWRU 全部七类故障标签上训练,包括 1D-CNN、ResNet-1D、WDCNN [10](使用宽首层卷积核抑制噪声)、1D-Transformer(patch 嵌入与四层编码器),以及具有梯度反转的域对抗神经网络 DANN [48]。前四个模型使用全部四种负载,按标准监督学习训练 50 个 epoch,采用 Adam,学习率 0.001、权重衰减 、余弦退火调度。DANN 额外接收各目标域无标签振动信号,进行 100 个 epoch 域对抗训练; 按标准公式 调度,域损失权重为 0.5。为减轻不同试验台幅值尺度造成的域偏移,所有五个基线在训练和推理、源域和目标域均统一使用逐样本 z 分数归一化。五个模型在 CWRU 上都达到近乎完美的准确率(至少 99%),然后不重训,直接用于目标域。该设计确保基线失败反映真实域偏移坍缩,而非训练不足。表 2 实验配置 4.2. 域内评估首先使用三个 LLM 在 CWRU 上评估完整系统。各模型处理完全相同的输入,包括自编码器输出、工具实现及提示模板,唯一区别是推理骨干。结果见表 3。三种 LLM 在五次独立运行中的平均总体准确率均超过 92%,表明框架在不同推理骨干间具有良好通用性。Qwen-Plus 最高,为 ;GPT-4.1 为 ;Claude Sonnet 4.6 为 。模型间仅相差 2.7 个百分点,运行间标准差不超过 1.2 个百分点,共同表明性能主要由物理工具链、感知机制和验证协议等架构决定,而非特定 LLM 或随机推理路径。每个值为在相同 140 样本上独立运行五次的均值 ± 标准差,温度为 0。较小的总体标准差(0.6–1.2 个百分点)说明诊断输出稳定。类别方差集中在 OR014,三个 LLM 的标准差均为 3.5 个百分点,与 4.6.5 节 BPFO/BPFI 谐波混叠分析一致:推理路径恰在谐波数相同的边界样本上略有变化。表 3 CWRU 上多 LLM 诊断准确率(%, ) 三模型有一致规律。IR014 和 OR021 达到完全或近乎完全准确,IR014 均为 ,说明中大尺寸故障产生强且无歧义的谱特征,不依赖后端即可可靠识别。OR014 最困难,三者平均准确率分别为 85.0%、75.0%、75.0%,标准差却完全相同(3.5 个百分点),说明困难来自确定性的物理原因,即该故障尺寸下谐波重叠,而非随机 LLM 行为。按第 4.6 节分析,0.014 英寸故障下 BPFO 高阶与 BPFI 低阶谐波重叠,是 CWRU 6205 轴承的基本几何性质,而非框架局限。物理推理的根本优势是特征频率随轴转速线性缩放。我们用四个 CWRU 负载的全部六个两两组合验证(表 4)。对每一对,系统先在第一转速诊断,根据转速比预测特征频移,再验证第二转速包络谱中是否出现预测频率。六对均通过,观测频移比与理论转速比误差不超过 0.3%,确认框架依据物理原理而非模式记忆。跨模型结果也间接支持提示鲁棒性。三个后端的预训练语料、指令遵循及推理策略差异显著,却处理相同提示模板;仅 2.7 个百分点的范围,说明能力主要由结构化物理工具链而非提示措辞细节决定。提示指导工具选择和解释,FFT、包络解调及谐波匹配等数值结果不依赖提示措辞。表 4 跨工况验证(六对) 4.3. 消融研究使用 Qwen-Plus 在 CWRU( )上进行系统消融,沿两个正交方向组织。第一方向比较完整系统与“纯确定性”流程,隔离 LLM 编排的总体贡献。后者采用完全相同的物理工具,但将所有 LLM 决策替换为人工规则:按最大 z 分数选频带,按谐波数最大值分类,并内置 BSF/FTF 审慎规则及规则式边带消歧。第二方向移除各子组件,包括感知智能体、诊断智能体三个子模块及验证智能体两个子模块,定位贡献来源。另设移除全部工具、只保留 LLM 的合理性检查。结果见表 5 和图 4。表 5 消融研究 图 4 CWRU 消融结果(Qwen-Plus, ,一次代表性运行)。(a)各配置总体准确率;(b)七类故障的分类别准确率;(c)相较完整系统的性能下降( )。缩写:TL,全部工具;BA,频带归因;PD,纯确定性;PA,感知智能体;FM,频率匹配;KB,知识库;VA,验证智能体;SD,边带检测。消融揭示三层发现。第一层:LLM 编排比相同工具的规则编排提高 12.2 个百分点。 完整系统 95.0% 减去纯确定性 82.8%,对应表 5 移除 LLM 后的 。纯确定性明显弱于完整系统,但远高于任何无工具配置,确立 LLM 的真实贡献:作为能够选取分析参数、化解歧义的上下文感知编排器,而不是特征提取器。第二层:提升集中在物理歧义样本。 对明显且占主导的故障特征类别(OR007、OR021、IR014),确定性流程已接近完美,达到每类 19–20/20,LLM 仅增加 0–2 个样本;OR014 则由 7/20 恢复到 17/20,提高 50 个百分点,贡献了整体差距的大部分。按 4.6.5 和图 8(c),OR014 恰好存在 BPFO 三倍频与 BPFI 二倍频在 0.7% 内重叠,导致人工最大值规则无法解决的谐波数近似并列。LLM 结合边带、感知频带指导及验证输出改善该类,这正是上下文推理区别于规则逻辑的多证据综合。第三层:各子模块均非冗余。 诊断子模块中频带归因影响最大,移除后下降 17.9 个百分点。没有感知频带指导,诊断智能体只能盲选参数,常选择不理想的包络谱带通范围。内圈类别影响最明显:IR007 为 12/20、IR014 为 15/20、IR021 为 14/20,因为冲击故障特征需要精确选择高频频带才能有效包络解调。验证智能体贡献 10.0 个百分点,其中边带检测单独贡献 7.1 个百分点。OR014 再次受影响最大,无验证时从 17/20 降为 3/20,确认边带消歧是 LLM 编排解决谐波混叠的主要机制。仅 LLM 的合理性检查为 12.9%:通过常识正确识别 18/20 个正常样本,但全部故障类别为 0%,证明仅靠文本描述不能完成振动诊断。此项排除“LLM 仅凭足够先验即可无计算支持分类”的简单假设;决定性消融仍是纯确定性对照,而非此检查。4.4. 感知智能体评估为验证 LLM 自适应异常感知(贡献 C2),在 Paderborn 比较三种感知配置。该跨域挑战要求在一种轴承(SKF 6203,K001–K004)上训练的自编码器 MSE 阈值,泛化到运行特征不同的测试轴承。配置 A:固定 MSE。 当 MSE 比值大于 1.0,即重构误差超过正常基线的 ,触发异常,代表标准自编码器异常检测。配置 B:固定 MSE + 峭度。 在 A 基础上加入规则:MSE 比值大于 1.0,或滤波峭度大于 20,即触发异常。峭度在 2000–10000 Hz 带通滤波信号上计算,以提取冲击故障内容并抑制低频旋转成分。配置 C:LLM 自适应。 完整系统由感知 LLM 综合 MSE 比值、频带 z 分数、峭度、峰值因子及跨域上下文,进行自适应异常评估,即本文方法。总体与分类别结果见表 6。表 6 感知智能体配置比较(Paderborn,23 个轴承)。上表为总体指标,下表为分类别结果;Acc 为分类准确率,Det 为检测率,单位均为 %。 结果呈现感知复杂度提升及其对下游诊断的影响。仅依赖 MSE 的 A 在跨域下严重退化:由于自编码器训练于不同轴承类型,大多数故障轴承的 MSE 比值低于 1.0,通常为 0.4–0.8 倍,使阈值不可靠。如图 5,Paderborn 健康与故障轴承的比值分布大量重叠,大多数故障位于阈值以下,与源域可获得的清晰分离形成对照。因此 A 检测率仅 61.9%,健康准确率为 0%;两健康轴承均误判,其中一个触发 MSE 阈值,而系统缺少推翻误报的上下文推理能力。B 引入反映冲击故障的物理指标峭度,将检测率提高至 95.2%。Paderborn 故障轴承滤波峭度范围为 7.27–301.87,多数超过 20。但激进检测牺牲特异性:B 分类准确率仅 43.5%,因为无法区分故障类型。峭度规则标记异常,却不指导分析频带,诊断智能体因而常把外圈判为内圈,真实磨损外圈准确率仅 14.3%。C 取得最佳总体平衡:分类准确率 63.8%、检测率 90.5%、误报率 33.3%。跨域感知使 LLM 认识 MSE 不可靠,主要依赖峭度与频带 z 分数,将 MSE 作为次要参考。C 可正确识别 66.7% 健康轴承,而 A 为 0%,同时保持高故障检出。B 与 C 检测率相差 4.7 个百分点,反映 C 更保守、更智能的阈值判断:并非高峭度一律报异常,而是综合其他指标评估是否真正意味着故障。 图 5 Paderborn 各轴承类别的自编码器 MSE 比值分布。多数故障轴承比值低于异常阈值(1.0 处虚线),说明固定 MSE 阈值为何在跨域时失效。4.5. 向未见目标数据集的无故障标签迁移最严格评估是跨数据集无故障标签迁移。对 Paderborn、JNU、MFPT 各目标域,自编码器仅用该域正常信号重新训练;物理工具链、提示模板及所有智能体逻辑与 CWRU 完全相同,只更新轴承几何与转速。这一设置没有可直接比较的基线:表 1 的现有 LLM 方法在训练或微调中都需要故障标签,因此不适用。加入在 CWRU 全标签上训练、源域准确率近乎完美的五个监督基线,是为了展示域偏移坍缩的严重性与普遍性。DANN 还额外接收目标域无标签信号进行对抗训练,信息严格多于其他基线。结果见表 7 与图 6。表 7 三个目标域上的跨数据集无故障标签评估 注:全部基线逐样本 z 分数归一化,以减轻幅值尺度域偏移。DANN 额外接收无标签目标域信号,训练 100 个 epoch,余弦学习率调度、权重衰减 、域损失权重 0.5,信息严格多于其他基线。“全部 → X”表示坍缩为单一恒定类别、无判别能力;“部分”表示输出分散但无判别力;× 表示诊断坍缩。Paderborn 的 43.5%–47.8% 恰对应多数类比例(外圈 43.5%、内圈 47.8%),等同于始终预测多数类。JNU 的 25.0% 同样对应正常样本比例;DANN 的 8.3% 低于四个均衡类别的 25% 随机猜测水平,表明负迁移。全部监督模型在 CWRU 七类全标签上训练,测试按文件以 20 片段多数投票评估。 图 6 三个目标数据集 Paderborn、JNU、MFPT 的跨数据集迁移评估。(a)故障检测率;(b)分类准确率。结果有三个发现。第一,无任何故障标注样本,框架在三个域上分别取得 90.5%、96.3%、100.0% 检测率,验证正确指定几何与转速即可迁移诊断能力。第二,五个监督基线,包括获得目标域无标签数据的 DANN,均无法泛化,表现为恒定类别坍缩,或分散但无判别力的输出。Paderborn 的表面准确率恰等于多数类比例,表明域内模式记忆无法跨越设备边界。第三,虽然分类准确率随频谱条件变化,MFPT 为 96.7%,JNU 为 63.0%——其低转速将故障谐波压缩至重叠频段——但检测率始终超过 90%,说明即使下游类型判别受限,多指标感知仍可靠识别异常。比较的公平性需要说明。DANN 用正常及故障混合的无标签目标域信号进行对抗特征对齐,信息严格多于只用正常信号的本文框架。尽管如此,DANN 在三域仍坍缩,准确率为 43.5%、8.3%、25.0%。这说明跨域失败并非目标数据不足,而是诊断知识编码于绑定源域分布的权重所导致的根本局限。本文用运动学方程、谐波匹配等物理公式编码知识,这些规则在设计上不随域偏移变化。结果还揭示超越准确率指标的结构性认识:轴承状态评估本质是具有非对称运行后果的两级决策。第一级故障检测回答“是否需要关注”,与安全直接相关;漏检可能造成非计划停机、连锁损坏或安全事故。第二级分类回答“哪个部件损坏”,用于维护规划;即使分类错误,也会触发检查,在拆解后发现真实状态,带来的是增量成本而非灾难性后果。本文在无故障标签条件下,以 90.5%–100.0% 检测率应对安全关键层级。Paderborn 与 JNU 分类准确率较低,与黑箱错误具有本质差别:每个诊断携带匹配谐波、理论偏差、边带结果等证据链,使工程师可独立评估歧义。分类不确定性由无声失败转化为可审计建议,适于保留人类专业判断的安全关键系统。与监督基线及 LLM 方法比较(4.3、4.7)还产生更广泛的架构启示:监督模型的知识编码于源域权重,分布不成立时知识不可访问;本文将知识编码于运动学方程、包络解调与谐波匹配等物理关系,避免这一问题。与监督 LLM 的参考比较显示,域内准确率与完全监督方法差距在 1.7 个百分点以内,大域差距下跨域检测更稳健。实际部署只需要正常数据、几何与转速,无需故障样本、重训或微调。4.6. 诊断推理分析4.6.1. 系统性审计记录评估为超越个案评估审计记录质量,我们分析全部 140 个 CWRU 测试样本的完整输出日志(Qwen-Plus,一次代表性运行),见表 8。表 8 审计记录完整性统计(CWRU,Qwen-Plus, ) 结果证实,框架在完整测试集系统性生成可审计证据链。边带分析主要在内圈类触发:IR007 为 3/20、IR014 为 8/20、IR021 为 6/20,与内圈缺陷产生轴频幅值调制的机理一致。12 个中置信度诊断对应谐波证据模糊样本,为人工复核提供明确不确定性信号。唯一不完整样本为超时,说明即使失败也可追溯,而非无声发生。为深入展示推理,我们给出两个案例、物理证据包络谱可视化以及系统性错误模式分析。4.6.2. 案例 1:标准正确诊断考察 1797 RPM 的 OR007,即直接明确的外圈故障。感知智能体在 3500–4000 Hz 频带发现极端异常,z 分数大于 160 万,触发诊断。诊断智能体依次计算特征频率:BPFO 107.4 Hz、BPFI 162.2 Hz、轴频 29.95 Hz;对异常高频段计算包络谱;提取 15 个峰;执行确定性匹配。结果无歧义:BPFO 匹配五阶谐波,依次为 107.5、215.5、323.0、430.5、538.0 Hz,平均匹配误差 0.28%,单个误差 0.1%–0.36%,见图 8(a)。BPFI 和 BSF 各仅一阶匹配,FTF 为零。知识库确认至少五阶谐波表示第 3 阶段(进展期)外圈故障。整个诊断仅需六次工具调用,50.6 s 完成。此例体现物理推理范式:LLM 不直接处理原始振动,而是编排由几何频率计算、信号处理到确定性匹配的工具链,综合输出为连贯诊断。平均 0.28% 误差体现以计算工具而非模式识别支撑推理时的高精度。4.6.3. 案例 2:边带消歧考察 1797 RPM 的 IR007(图 7 右),存在谐波歧义。包络谱和匹配后,BPFI、BPFO、BSF 各有两阶谐波,单纯计数无法打破三方并列。LLM 没有默认选最常见故障,而是主动调用边带检测,识别 132.0 Hz(BPFI 减轴频)和 191.5 Hz(BPFI 加轴频)两处清晰边带,间隔恰与 对应。物理依据是内圈缺陷随轴旋转,每转一周通过承载区一次,使 BPFI 信号以轴频幅值调制,在 产生边带。原文将其描述为内圈故障特有、外圈及滚动体缺陷不具备的特征。框架识别并应用该原理,将模糊诊断提升为高置信内圈故障,展示超越简单频率匹配的推理。该案例直接对应表 5 移除边带检测下降 7.1 个百分点:若无此机制,涉及谐波歧义的 7.1% 样本将依赖默认启发式或随机机会,而非物理推理。 图 7 1797 RPM 下两个代表性样本的诊断推理链:标准外圈故障 OR007(左);需边带消歧的模糊内圈故障 IR007(右)4.6.4. 包络谱证据可视化图 8 展示四个代表案例的包络谱。(a)(b)对应上述案例:OR007 有五阶 BPFO 谐波,匹配误差小于 0.4%,是明确外圈特征;IR007 有三阶 BPFI 谐波,并伴随 BPFI ± 轴频边带,确认旋转内圈缺陷的幅值调制。(c)揭示主要错误来源:OR014 的 BPFO 三倍频 322.1 Hz 与 BPFI 二倍频 324.4 Hz 仅相差 0.7%,引发 4.6.5 节的系统性误分类。(d)健康轴承没有占主导的故障谐波,作为基线参考。 图 8 1797 RPM 下四个代表样本的包络谱。(a)OR007:清晰识别五阶 BPFO 谐波,误差小于 0.4%;(b)IR007:三阶 BPFI 谐波及轴频边带(BPFI ± ),确认幅值调制;(c)OR014:BPFO 三倍频(322.1 Hz)与 BPFI 二倍频(324.4 Hz)相差 0.7% 的混叠区域,是内外圈误分类根因;(d)正常:噪声底之上无占主导的故障谐波。4.6.5. 错误分析CWRU 实验的七个误分类样本(错误率 5.0%)揭示两种失效模式,见表 9。表 9 错误分类(Qwen-Plus,CWRU,140 样本中 7 个误分类) 主要类型为 BPFO/BPFI 谐波混淆,占七个错误中的五个,且全部影响 OR014。0.014 英寸故障下,BPFO 三阶约 ,与 BPFI 二阶约 相差 0.7%。采用 2% 容差的匹配工具,会合理地将约 322–323 Hz 的观测峰同时匹配至 BPFO 三阶和 BPFI 二阶,产生额外 BPFI 计数,偏向内圈诊断。这是 CWRU 6205 几何频率比引起的基本物理限制,而非框架缺陷。可行缓解包括收紧容差(代价是漏掉真实匹配),或按幅值加权谐波评分,留待未来研究。单个误报案例说明统计上的不可避免性:每个样本提取 15 个峰,对四类故障频率、各五阶谐波应用 2% 窗口,噪声峰偶然对齐模板的概率不可忽略。正常类误报率为 ,与此参数空间预期偶然匹配率一致,说明增加幅值阈值或峰显著性等判别特征可能进一步降低误报。4.6.6. 频率匹配容差分析第 3.4.2 节的 容差反映三种实际因素。第一,不同负载下转速波动;CWRU 四工况为 1730–1797 RPM,范围 3.7%,实际特征频率即使在同一数据集中也随转速变化。第二,轴承几何公差与内部游隙变化,使名义参数理论频率偏离实测谱频率。第三,有限 FFT 分辨率产生离散化下限:12 kHz、16384 点时, 约为 BPFO 107.4 Hz 的 0.7%。因此 是在容纳综合不确定性与限制偶然匹配间的工程折中。从 4.6.5 错误分析可评估其敏感性。收紧至 会减少偶然落入窗口的噪声峰,可能降低正常类误报;但不能消除主要 OR014 错误,因为两组谐波仅相差 0.7%,仍落在两侧 窗口内。放宽至 会增加各类噪声偶然匹配,抬高误报及类间混淆。当前 取得 95.1% 总体准确率,主要错误来自 CWRU 6205 特定几何巧合,而非系统性框架限制。第 4.6.7 节参数扰动提供进一步定量证据:几何参数扰动 至 时,理论频率变化相当于在真实参数不确定性下测试容差窗口。表 10 确认,在符合工业参数精度的扰动下保持鲁棒。针对不同轴承几何细致优化容差,是未来部署定制方向。4.6.7. 轴承参数扰动分析物理推理依赖几何参数和转速计算理论频率。实际中,制造公差、长期磨损或文档不完整会使名义值偏离真实值。为评估鲁棒性,我们在 CWRU 上系统扰动。每个扰动幅度为 1%、3%、5%,对滚珠直径 、节圆直径 、滚动体数 (四舍五入至最近整数)以及转速 独立施加随机乘性扰动: 各幅度使用五个独立随机种子重复测试,每次用 Qwen-Plus 诊断全部 140 个样本。结果见表 10。表 10 参数扰动分析(CWRU,Qwen-Plus,每次 ,每级五次) 第一,分类准确率随扰动平稳下降: 仅下降 4.1 个百分点, 下降 7.8 个百分点,没有突发坍缩。 代表厂家数据表精度,此时与名义基线在统计上无法区分(均 95.1%),确认 匹配容差充分吸收小不确定性。较大扰动下主要是转速影响,因它成比例平移所有特征频率;转速偏差超过匹配容差时,高阶谐波落出窗口,匹配计数减少、分类歧义增加。第二,更重要的是,即使 ,检测率仍超过 97%。这来自感知—诊断分离:感知依赖直接由信号计算的 MSE 比值、频带 z 分数与峭度,完全独立于几何参数。扰动只影响诊断频率匹配,削弱分类但不削弱上游检测。结果验证两级设计:即便参数不确定性远超典型工业条件,框架仍可靠识别故障存在,并保持合理分类准确率。4.7. 与基于 LLM 的诊断方法比较为将本文置于新兴 LLM 轴承诊断研究中,我们与 Tao 等 [16]、BearLLM [18]、DiagLLM [17] 作参考比较。由于这些方法需故障标签训练/微调、无法直接用于无标签条件,因此采用原论文报告结果,而非受控复现。表 11 从三个维度汇总比较。表 11 与基于 LLM 的轴承故障诊断方法的参考比较 标签效率。 本文在 CWRU 无故障标签达到 95.1%,仅比 Tao 等使用全标签和 LoRA 的 96.85% 低 1.7 个百分点,且是更难的七类而非四类设置。BearLLM 达 100%,但需九个数据集故障标签;其 CWRU 留一跨数据集准确率为 90.3%,即使使用其他八个有标签数据集仍比本文低 4.8 个百分点。这种较小准确率差距对应部署需求的质变:设备投用后即可部署,不必等待故障数据。跨域鲁棒性。 在最接近工业部署的有限数据情形,Tao 等即使有三个源域完整标签和 10% 目标域故障数据,准确率仍降至 64.6%–86.9%。本文监督更少,在相同目标域故障检测率达到 90.5%–100%,说明物理推理在大域差距下比学习模式表示更平稳地退化。需要指出,本比较来自不同论文,数据划分、类别定义和评估协议不同,应视为展示标签效率与诊断能力的范式级权衡,而非受控的一对一比较。在相同实验条件下复现,仍是未来基准研究的重要方向。5. 结论本文提出可信多智能体 LLM 框架,应对安全关键机械状态监测部署智能体 AI 的保障要求。通过将 LLM 从模式分类器重新定位为确定性物理工具链的推理编排器,实现完全无故障标签诊断,并把四种可信属性直接落实于架构:LLM 推理—信号计算分离、可审计推理溯源、安全单调反思纠正及组件失效平稳降级。主要结论如下。无故障标签的物理推理在 CWRU 取得 准确率。其中确定性工具链贡献 82.8%,LLM 编排额外贡献 12.2 个百分点,这是 4.3 节单次运行消融结果。在三个未见目标域的无标签迁移中,平均故障检测率 95.6%;五个监督基线均未能泛化,输出退化为恒定类别或多数类先验。带 LLM 自适应感知的三智能体架构,将跨域分类准确率从 21.7% 提高至 63.8%,误报率从 100% 降至 33.3%,在分布偏移下有效平衡灵敏度与特异性。多层证据一致性验证贡献 10.0 个百分点总体准确率,边带消歧单独贡献 7.1 个百分点,有效缓解 BPFO/BPFI 谐波重叠引起的系统性内外圈误分类。与三种 LLM 方法的参考比较显示,本文无需故障标签,域内准确率与完全监督替代方法差距在 1.7 个百分点以内,并在大域差距下提供更稳健跨域故障检测。此外,三 LLM 后端五次运行均值为 92.4%–95.1%,差距仅 2.7 个百分点,标准差 0.6–1.2 个百分点,确认能力主要由物理工具链架构而非某一 LLM 决定。需要指出若干局限。每样本 30–90 s 的端到端时间主要来自 API 延迟,因此框架定位为定期检查和离线维护支持,而非连续实时监测。框架尚未验证多故障同时存在的复合工况,例如外圈点蚀与保持架磨损同时发生。此时不同谱特征可能重叠,谐波逻辑会报告多个类型及各自计数。多证据综合或许能标记并发故障,但边带消歧原为内外圈二分设计,面对多种故障谱叠加的行为未经测试。通过幅值加权谐波评分和多故障假设检验扩展至复合诊断,是优先方向。防回退约束虽防止真实故障回退为正常,却引入互补风险:若初始诊断因偶然谱峰把健康轴承判为故障,约束会阻止系统自行纠正误报。这反映漏检代价远高于误报的安全关键非对称性。可考虑置信度门控放宽:若原诊断低置信且证据不完整,例如只匹配单个低 SNR 谐波,就放宽约束;也可设计人工覆盖路径,让工程师依据实体检查证据撤销已锁定诊断。与 LLM 方法的参考比较采用不同论文、不同划分、类别及协议,应解释为范式级权衡,而非受控对照。LLM 随机性虽有界,运行间不超过 1.2 个百分点,工业部署仍可能受益于多次诊断投票。最后,物理推理依赖准确几何参数和轴转速,长期磨损可能使其偏离名义值。4.6.7 节仅确认 下平稳退化,尚未覆盖严重磨损或错误文档引起的更大偏差。未来工作聚焦六个方向:以齿轮啮合频率公式替换运动学方程,将物理工具链扩展到多部件旋转机械。融合振动、声发射、电机电流等多传感器多模态信息,丰富异常指标与谐波消歧证据。引入 RAG,积累历史诊断案例以持续学习。构建物理与数据驱动混合架构:无标签框架作为初筛,故障数据可用后由微调模型细化分类。迁移至 7B–13B 轻量级本地 LLM,实现边缘部署;框架仅传输文本提示而非原始信号,自然适应工业“数据留在本地”的严格隐私要求。与制造执行系统(MES)和数字孪生平台集成,对智能产线多设备进行系统级状态监测。编辑:Kira校核:李正平、陈凯歌、Leo、曹希铭、赵学功、白亮、任超、海洋、陈宇航、Tina、王金、赵诚、肖鑫鑫、张优该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除来源:故障诊断与python学习

    未登录
    还没有评论
    课程
    培训
    服务
    行家
    VIP会员 学习计划 福利任务
    下载APP
    联系我们
    帮助与反馈