首页/文章/ 详情

开源论文推荐 | 一种基于独立样本的个体泛化框架,旨在构建更合理的故障诊断基准

8月前浏览1288

    域偏移是故障诊断中不可避免的现象,本期给大家推荐一篇沈卫明教授团队的论文:一种基于独立样本的个体泛化框架,旨在构建更合理的故障诊断基准。

    本文探讨了传统数据驱动故障诊断评估基准中一个被忽视的问题,即训练集和测试集中通常来源于同一个设备个体。个体差异(例如制造误差,装配差异)可能导致潜在的被忽视的域偏移问题。为了建立更合理的故障诊断基准,本研究探索了一种基于独立样本的个体泛化框架。


    论文基本信息

    论文题目:An individual generalization framework based on independent samples towards a more reasonable fault diagnosis benchmark

    论文期刊:Computers in Industry

    论文日期:2025

    论文链接

    https://doi.org/10.48550/arXiv.2506.05854

    作者:Yiming He , Weiming Shen 

    机构:State Key Laboratory of Digital Manufacturing Equipment and Technology, Huazhong University of Science and Technology, Wuhan 430074, China

    作者简介:

    沈卫明,加拿大工程院院士,华中科技大学教授,IEEEFellow。现兼任 IET CIM 主编,IEEE TS MC、TASE 等杂志副主编。主要从事协同设计与制造、智能制造、生产调度、物联网、大数据、智能体理论与应用等方面的研究。在国际杂志及会议上发表论文 450 多篇,其中SCI检索 200 多篇,谷歌论文被引 28000余次(其中单篇论文被引最高 900 余次),谷歌学术 H 指数 83。获IEEE加拿大唐纳工业领袖奖和加拿大联邦大楼最佳能效奖等奖项(来源:华中科技大学教师主页)  

    目录

    摘要

    1 引言

    2 相关工作

    3 所提个体泛化框架

    4 实验结果与讨论

        4.1 一维卷积分块

        4.2 多尺度分块交互

        4.3 粗粒度时序定位器

    5 实验与讨论

        5.1 案例1:帕德博恩大学轴承数据集

        5.2 案例2:实际生产线整机级电机数据集

        5.3 分析与讨论

    6 总结

    摘要

    域偏移是工业信号故障诊断中不可避免的现象。本文探讨了传统数据驱动故障诊断评估基准中一个被忽视的问题,即训练集和测试来源于同一个设备个体。为了建立更合理的故障诊断基准,研究了一种使用独立测试个体的个体泛化框架。此外,还应用了一种改进的轻量级Transformer,以增强动态全局特征提取和无关信息过滤。在帕德博恩大学轴承数据集和从实际生产线收集的机器级电机数据集上进行了综合实验。结果表明,传统基准无法有效评估对故障无关特征的筛选能力和对新个体的泛化能力。所提出的轻量级Transformer实现了最高的泛化性能,具有巨大的应用潜力  

    关键词评估基准故障诊断独立样本Transformer个体泛化

    1  引言

    智能故障诊断(IFD)是提升现代制造系统可靠性、稳定性与安全性的关键维护环节。随着深度学习(DL)技术的发展与落地,数据驱动式故障诊断已成为该领域的主流研究方向(Peng 等,2024;Xiao 等,2023)。  

    目前已有多个实验台被搭建,用于采集旋转机械的传感器信号,以此验证各类基于深度学习的故障诊断方法。行业内公认的评估基准是采用滑动窗口对不同时段的信号进行连续采样,进而构建训练集与测试集(Wang 等,2024a)。当前几乎所有故障诊断方法均沿用这一评估基准,并基于公开数据集开发了大量算法模型(Yu 等,2023;Li 等,2025)。  

    但鲜有研究关注该评估基准的合理性。由于深度学习技术最初广泛应用于图像和自然语言处理领域,其数据集构建基准便被直接迁移至故障诊断领域,未经过针对性的充分考量与深入探索。考虑到振动信号的周期性,同一设备不同时段的信号应被视为相同或相似样本,这类样本的干扰信息相近,数据分布偏移极小。在源设备个体上完成的模型验证,无法评估模型对实际工业场景中复杂干扰因素的抵抗能力,例如制造不一致性、装配差异以及数据采集环境中的噪声等(He 和 Shen,2023)。S mith 与 Randall 针对凯斯西储大学(CWRU)轴承数据集展开研究,得出了一项重要结论(S mith 和 Randall,2015):他们在传感器信号中观测到部分与轴承故障无关的脉冲信息,这类信息疑似与试验台的机械松动,以及每次轴承测试时的反复拆装操作相关。面对此类特殊的信号特征,现有评估基准可能引发两大问题  

    • 1)故障诊断精度并不一定可靠:现有评估基准的一大潜在缺陷是验证集与测试集均来源于同一设备个体,从而忽略了上述干扰信息。多数基于设备个体标签的分类结果(例如 10 个轴承对应 10 个不同标签),仅能说明深度学习模型可区分不同的设备个体,无法保证模型实现了真正的故障分类,尤其是在干扰特征较为明显的场景下。    

    • 2新设备个体泛化性能差:对于新的设备个体(即便属于同一类别),域偏移现象不可避免,在更复杂的实际工业环境以及存在装配误差的多部件系统中,这一现象会尤为突出。事实上,针对新个体的故障诊断可被视为一类全新的域泛化任务,需要对特征进行校正或抑制无关特征,在当前故障标签稀缺的行业现状下,这一需求更为迫切。    

    然而这些问题在后续研究中并未得到足够重视。一项重要的基准研究(Hendriks 等,2022)探讨了现有评估框架的潜在缺陷,并进一步验证了上述观点。另一项近期的基准研究也再次指出,相较于在源设备上开展诊断,实现不同个体设备间的泛化诊断具备更高的实际意义与应用价值(Zhao 等,2024)。因此,亟需更合理地构建数据集,确保训练集与测试集的设备个体无重叠,使其更贴近真实工业场景。此外,考虑到已有性能优异的物理模型方法可实现单一轴承的故障诊断,建议将研究目标转向更复杂的系统 —— 多部件系统中不同设备个体间的域差异可能会表现得更为显著(He 等,2023a)。  

    本文进一步探究了不同设备个体间被忽视的域差异,并提出了一套面向数据驱动智能故障诊断的个体泛化框架。同时,设计了一款改进型轻量级 Transformer 模型作为该框架的可行解决方案。具体而言,本文通过特征可视化技术对比了传统基准与独立样本评估基准的差异;在 Transformer 模型中集成了一系列即插即用的改进模块,包括一维卷积分块(1-DCP)、多尺度分块交互(MPI)以及粗粒度时序定位器(CTL),这些模块的设计初衷是增强模型的动态全局特征提取能力与无关信息过滤能力。研究首先基于具有代表性的公开实验室数据集 —— 帕德博恩大学轴承数据集(Less meier 等,2016),完成了对该泛化框架的测试;随后将所提模型部署于实际工业电机生产线,实现了多部件系统的复合故障诊断。实验结果表明,在新的评估基准下,所提方法的性能显著优于其他对比的先进深度学习模型。本文的主要贡献可总结如下  

    • 1)本文揭示一种被忽视的潜在域偏移现象:不同设备个体固有的制造与装配差异,加之复杂的动态数据采集环境,即便对于同类型设备,也会导致不可避免的域差异。    

    • 2提出个体泛化框架,将其作为更合理的故障诊断基准,以实现新设备个体上的可靠故障分类;并基于特征可视化技术,进一步论证与剖析了传统故障诊断基准的局限性。    

    • 3)设计一系列即插即用的模块嵌入先进的深度学习模型,以增强动态全局特征提取能力,降低个体差异的影响,提升模型泛化性能;通过帕德博恩大学轴承数据集与实际生产线上的工业电机整机数据集,验证了该模型和基准的有效性。    

    论文剩余部分的结构安排如下:

    • 第 2 节简要综述智能故障诊断的主流方法;

    • 第 4 节揭示设备间的域偏移现象并介绍所提泛化框架;

    • 第 5 节阐述所提模型的详细结构;实验步骤与相关分析在第 5.3 节呈现;

    • 第 6 节对全文进行总结。

    2 相关工作

    卷积神经网络(CNN)已然成为智能故障诊断(IFD)领域的主流模型架构,诸多 CNN 变体已在各类实验室数据集上完成验证。Zhang 等人提出的宽第一层卷积核 CNN(WDCNN)(Zhang 等,2017)是经典且易于实现的模型,被广泛用于对比实验。Li 等人提出了一种新颖的基于注意力机制的深度元迁移学习方法(ADMTL),可利用目标域的少量标注样本完成诊断任务(Li 等,2023)。Lu 等人设计了嵌入先验知识的卷积自编码器(PKECA),适用于跨转速、时变转速及小样本诊断场景(Lu 等,2025)。Ye 等人研发了多传感器残差卷积融合网络(MRCFN),通过多传感器特征信息的互补与校准实现高精度诊断(Ye 等,2025)。Wang 等人基于经典 ResNet 网络,设计了多源域特征适配网络(MDFAN),用于时变工况下的轴承故障诊断(Wang 等,2022)。这些研究通过借鉴和改进图像与自然语言处理领域的先进结构及算法,大幅提升了深度学习模型的分类性能与计算效率。  

    近年来,一种名为 Transformer 的新型注意力机制架构革新了深度学习领域。部分成功应用案例(Yoon 等,2023)已证实其强大的全局信息提取与多任务泛化能力,例如广受欢迎的 ChatGPT(Radford 等,2018)。学者们开始探索 Transformer 在故障诊断领域的性能表现(Ding 等,2022;Jin 等,2022)。Lian 等人提出了融合声振信号的残差注意力引导视觉 Transformer(RAViT),通过捕获并利用序列的长程依赖关系提升诊断效率(Lian 等,2025)。Wang 等人提出了基于 Transformer 架构的轻量级模型 C-ECAFormer,在强噪声干扰的小样本故障诊断中表现优异(Wang 等,2023)。Zhou 等人设计了轻量级 Transformer,以及带类间排斥判别的同质广义对比学习策略(HGCL-ICRD),用于变工况下的轴承故障诊断(Zhou 等,2025)。  

    Chen 等人提出轻量级卷积 Transformer(LCT),用于工业机器人的复合故障诊断(Chen 等,2023)。Weng 等人设计了多传感器融合 Transformer(MsFT),并在两个实验室数据集上完成验证(Weng 等,2023)。Wang 等人提出了特征增强型端到端故障诊断模型 PSparseFormer,在强噪声干扰下的特征提取与增强任务中效果显著(Wang 等,2024b)。但如前文分析,基于传统评估基准的实验所获得的分类精度,无法保证模型在未见过的新设备上的故障分类可靠性。在多数工业场景中,实际需求是对新设备个体进行诊断,而非在同一设备上反复训练与测试,全局域偏移与动态测试噪声干扰不可避免。此外,这一需求也限制了部分基于目标域信息引导训练的迁移学习与域适配技术的应用。  

    当前研究尚未充分挖掘 Transformer 在故障诊断领域的潜力,尤其是在被忽视的域泛化挑战方面。相较于 CNN 架构对空间变化敏感、长程特征感知能力薄弱的缺陷(Arkin 等,2023),Transformer 具备卓越的全局信息交互能力与注意力增强的特征过滤机制,非常适合在不同个体的域泛化任务中动态感知全局特征、抵抗多类干扰因素。

    3 所提个体泛化框架

    本节进一步分析和探讨个体差异对故障诊断任务的影响,并借助可视化技术对相关现象进行验证。  

    如引言所述,不同设备个体之间存在当前研究中被忽视的、由额外信息干扰引发的域差异。即便是同类型的目标诊断设备,依然存在域泛化需求 —— 故障诊断(即便在相同工况下)不能简单视为源域特征提取任务,因为诊断的对象始终是新的设备个体,而非在注入故障的实验设备上反复开展测试。  

    众所周知,旋转机械的信号具有明显的周期性。在同一设备上,不同时刻通过滑动窗口采集到的关键信息差异极小。图 1 展示了帕德博恩大学轴承数据集中 K005 和 K006 两种正常轴承的信号可视化结果:K005-1 和 K005-2 分别代表 K005 轴承的第一次和第二次测试数据,在传统评估基准中二者会分别作为训练集和测试集;K006 时域波形中的红色矩形标注的是 K005 轴承信号的幅值范围。可以看到,K005 和 K006 的时域信号存在显著差异,但 K005-1 和 K005-2 的波形却高度相似;在频域层面,即便是幅值低于 0.005 的噪声分量,二者也呈现出相似性。尽管 K005 和 K006 均为正常轴承,但二者在时域和频域均存在明显差异。

     

    图1. 正常轴承的时域和频域信号

     

    此外,本文对 K005、KI18 轴承不同时段的样本,以及 K006、KI21 轴承的样本进行了 T-SNE 可视化,结果如图 2 所示。能够发现,在传统任务设置下,训练样本和测试样本的域分布几乎没有差异;而在本文所提任务设置下,同类型的两个轴承样本则出现了分布差异(KI006 样本分布更分散,KI21 样本分布更集中)。

     

    图2. 不同任务设置下部分轴承样本的T-SNE可视化

     

    综上,传统诊断基准划分的测试集不能简单当作新样本,来评估故障诊断模型的真实泛化能力。当个体差异不可避免且愈发显著时,传统基准的缺陷会愈发突出,进而阻碍数据驱动故障诊断方法的推广与落地。图 3 对比了传统故障诊断框架与本文所提的独立样本故障诊断框架,其中     代表第    台设备的第    个样本。所提框架的显著改进在于,训练集和测试集之间不存在共享的设备个体信息。这种面向故障诊断的独立样本评估基准具有重要价值,且值得大力推广,后续实验将对此做进一步验证。

     

    图 3. 传统故障诊断框架与所提独立样本故障诊断框架的对比

     

    4 实验结果与讨论  

    本节提出了一系列先进的即插即用模块,并将其集成到 Transformer 基模型中,形成泛化诊断的解决方案(如图 4 (b) 所示)。该模型主要包含三个模块:一维卷积分块(1-DCP)、多尺度分块交互(MPI)和粗粒度时序定位器(CTL)。

     

    图 4. 所提模型与视觉 Transformer(ViT)的对比

     

    4.1 一维卷积分块

    最初用于自然语言处理的 Transformer 模型(Han 等,2021)可处理一维数据,看似适用于单通道信号。但实际中,信号输入长度通常可达数千个数据点,直接处理会产生极高的计算资源开销,因此不建议直接应用。这一问题在 Transformer 拓展至图像分类领域时同样存在(Dosovitskiy 等,2020)。研究人员采用分块分割方法解决该难题,实现了对二维数据的处理,但这类构造方式通常只是简单分割,仅能获取固定数量的空间位置信息。视觉 Transformer(ViT)的经典分块构造模式如图 4 (a) 所示。 本文设计的一维卷积分块(1-DCP),旨在为振动信号处理提供更灵活的标记器。具体而言,针对多通道信号,1-DCP 从时序维度同步处理所有通道,以此构建二维卷积分块(也可称为单向卷积分块)。该设计充分发挥了卷积的优势 —— 相较于线性运算,卷积具备更强的空间关联特征提取能力,且计算资源消耗更低;同时,它还结合了卷积的局部特征感知优势与 Transformer 的全局特征感知优势。此外,本文还提出了重叠卷积分块模式,能够更灵活地控制分块的类型与数量,从而捕获更丰富的信息。该设计的核心动机是兼顾输入信号的周期性与空间位置可变性。1-DCP 的数学表达式如下:

       

    式中,    代表第    个分块的第    维数据;卷积核权重    与信号点    进行卷积运算。卷积分块的尺寸为[     ,    ],其中    为分块长度,    为信号通道数。在每个分块区域内,会执行    次全域卷积以构建卷积分块     ,最终得到输出     。分块数量N可通过下述公式计算:

       

       为输入信号的长度。

    4.2 多尺度分块交互

    多尺度信息是振动信号的典型特征。 多尺度分块交互(MPI) 模块被设计为即插即用组件,用于进一步增强卷积分块之间的交互能力。MPI 模块由三个并行分支构成(如图 5 所示),卷积分块可基于空洞卷积实现不同尺度下的交互,其数学表达式为:

       

    分块交互器的尺寸为[     ,     ,    ],其中    为参与交互的分块序列间隔。需注意,为保证特征维度的一致性,需采用 “same” 填充方式,且     需与    保持一致。对所有分支的多尺度分块进行求和并计算均值:

       

    式中,     为并行尺度分支的数量。MPI 模块的最终融合输出为     。此外,模型还引入了内外位置编码(IEPC,He 和 Shen,2023),该编码可同时对标记间和标记内的位置信息进行编码:

       

    式中,    为随机初始化的可学习矩阵(考虑到通过滑动窗口构建的信号样本具有动态性);     为与     维度相同的矩阵,代表内外位置编码。

     

    图 5. 多尺度分块交互(MPI)模块的结构

    4.3 粗粒度时序定位器  

    为进一步强化关键有效特征并实现模型轻量化设计,本文将  粗粒度时序定位器(CTL) 设计为特征剪枝模块。其结构包含层归一化(LN)、残差连接和多头自注意力模块(MSA,Vaswani 等,2017),具体结构如图 6 所示。最终输出会通过 Tanh 函数映射至特定范围:

       

    式中,    。随后,计算每个分块的特征值之和作为权重:

    提取权重排名前    的分块索引    ,并获取对应的原始分块:

       

    剪枝后的分块会被输入至 Transformer 的骨干网络以进行进一步特征提取,其中多层感知机(MLP)的神经元数量分别设置为D和2D。模型采用带 10% 丢弃率(dropout)的一维全局平均池化(GAP)作为决策特征层,替代了传统 Transformer 框架中的展平(flatten)层与 MLP 层。 对于采用独热标签的单一故障诊断任务,模型使用基于 softmax 激活函数和类别交叉熵损失的分类器;对于多标签复合故障诊断任务,则采用基于 sigmoid 激活函数和二元交叉熵(BCE)损失函数的分类器(Ruby 和 Yendapalli,2020)。该分类器会对每个类别进行独立的二分类评估,将复合故障诊断任务转化为多个二分类任务。

     

    图 6. 粗粒度时序定位器(CTL)模块的结构

     

    5 实验与讨论

    本节在两个数据集的两种框架下对多个深度学习模型进行评估,并从多维度展开分析。实验代码基于 TensorFlow 2.6 框架编写,运行于搭载 NVIDIA GeForce RTX 3070 Ti 显卡的计算机。  

    选取的对比方法均为改进型卷积神经网络(CNN)模型,这类模型均提供完整源码,且通过特殊结构增强了泛化能力,便于研究人员快速评估与优化;同时,与部分大型深度模型相比,所选模型规模相近,可归为轻量化模型。在实际工业场景中,新目标设备个体的任何信息(无论是否标注)都难以在训练阶段获取,因此部分基于迁移学习和域适配的先进算法并不适用 —— 这类算法需借助目标域数据来引导迁移或域对齐过程。所有模型均采用 Adam 优化器训练 20 个 epoch,学习率设置为 0.01。本文选用的对比深度学习模型如下:  

    • 1)WDCNN(Zhang 等,2017):经典且受广泛关注的 CNN 模型,采用宽第一层卷积核,适用于滚动轴承故障诊断。    

    • 2TICNN(Zhang 等,2018):基于训练干扰的先进抗噪 CNN 模型,针对滚动轴承故障诊断任务设计。    

    • 3SRDCNN(Zhuang 等,2019):基于空洞卷积的先进 CNN 模型,用于滚动轴承故障诊断。    

    • 4MSRFCNN(He 等,2023b):结合空洞卷积与多阶段残差的先进 CNN 模型,适配整机级电机故障诊断。    

    • 5MS MCNN(He 等,2023a):搭载多尺度融合模块的先进 CNN 模型,可用于工业机器人与滚动轴承故障诊断。    

    • 6基准模型(Baseline):配备一维卷积分块(1-DCP)标记器的 Transformer 模型,核心结构参数参考 He 和 Shen(2023)的研究,具体如表 1 所示。    

    • 7B+CTL:在基准模型中嵌入粗粒度时序定位器(CTL)模块,用于消融实验。    

    • 8B+MPI:在基准模型中嵌入多尺度分块交互(MPI)模块,用于消融实验。    

    • 9B+MPI+CTL:在基准模型中同时嵌入 MPI 与 CTL 模块,用于消融实验。

    表 1 模型结构超参数    

     
     

    5.案例1帕德博恩大学轴承数据集

    帕德博恩大学轴承数据集(Less meier 等,2016)是具有代表性的实验室数据集,可用于泛化诊断验证。实验采用多个 6203 型滚动轴承(包含人工模拟损伤和实际损伤,如图 7 所示),该数据集在相关基准研究中被高度推荐(Hendriks 等,2022)。选取多组实际损伤轴承的振动信号,验证独立样本评估基准的有效性;具体而言,将 K006、KA22 和 KI21 轴承的信号作为目标域测试集(从未出现在训练集中),其余轴承作为训练集(如表 2 所示)。为每个轴承随机采样 1000 个样本,每个样本包含 5120 个振动数据点。首先探究一维卷积分块(1-DCP)中步长s和重叠率对模型性能的影响,通过不同结构的 Transformer 基模型获取未见过设备诊断精度(UDA);根据实验结果(如表 3 所示),最终确定步长s取值为 256。

     

    图 7. 帕德博恩大学轴承数据集的实验环境

    表 2 帕德博恩大学轴承数据集的轴承详情    

     

    表 3 步长与重叠率对泛化性能的影响

     

    开展对比实验,分别获取传统基准下的  源设备诊断精度(SDA) 和独立样本评估基准下的 UDA,10 次实验的平均结果如表 4 图 8 所示。可以看到,在传统基准下,所有测试的深度学习模型均能轻松实现接近 100% 的源设备诊断精度;但在新轴承个体测试中,模型无法维持预期的高精度,部分 CNN 模型的 UDA 出现大幅下降 ——WDCNN、TICNN 和 MSRFCNN 的 UDA 分别仅为 69.73%、73.64% 和 65.32%;通过先进结构提升感受野的 SRDCNN 和 MS MCNN,UDA 分别达到 85.44% 和 77.73%。得益于 Transformer 特有的全局注意力机制,嵌入本文所提模块的系列 Transformer 基模型展现出更优泛化性能;消融实验进一步表明,两个设计模块分别实现了 2.45% 和 2.08% 的 UDA 提升,同时让模型训练更稳定,最小标准差仅为 0.968。值得注意的是,10 次实验中 CNN 模型的训练过程极不稳定,导致 UDA 波动较大(例如 TICNN 和 SRDCNN 的部分实验结果 UDA 甚至超过所提框架),这说明 CNN 的特征提取具有较强随机性且缺乏过滤机制 —— 由于 CNN 的感受野局限于局部特征,且信号输入的空间位置具有动态性,其无法保证提取的每个局部特征都具备通用性和故障关联性。

    表 4 轴承故障诊断对比实验结果

     
     

    图 8. 所有模型的轴承诊断精度箱线图

     

    此外通过 T-SNE 可视化和混淆矩阵,从类别层面探究该现象(如图 9 图 10 所示)。图 9 中 “I-F” 代表内圈故障样本,“O-F” 代表外圈故障样本;可以看到,部分 CNN 模型存在明显的个体聚类现象(尤其是 WDCNN 和 TICNN),直观反映出其提取的特征包含与个体相关、与类别无关的干扰信息 —— 该现象越显著,模型越难提取类别相关特征,跨设备泛化能力也就越差。故障个体间的域偏移问题更为严重,部分模型对外圈故障 L14 的误报率较高。而所提模型将 15 个轴承个体划分为 3 个区域,类内距离最小,同类别轴承可聚集形成重叠的公共簇。图 10 的混淆矩阵进一步阐明各类别的分类结果:多数网络对未见过的正常样本识别精度优异,但对新故障类别(尤其是内圈故障)易与外圈故障混淆,这说明部分对比 CNN 模型的决策特征并非精准的故障关联特征。所提模型能够从有限的源域信息中提取更准确的故障关联信息,充分证明其具备更高的应用潜力。

     

    图 9. 所有模型决策层的 T-SNE 可视化。I-F 代表内圈故障样本;O-F 代表外圈故障样本    

     

    图 10. 所有模型的混淆矩阵

     

    5.案例 2:实际生产线整机级电机数据集  

    本节提供一个实际生产线的应用案例。整机级电机是更复杂的多部件旋转系统,通常存在复合故障和更复杂的工业噪声。选取某厂商维修检测现场 18 台 1000rpm 的电机开展测试,通过 4 个通道采集融合振动信号,采样频率为 25600Hz;将 1-14 号电机作为训练集,15-18 号电机作为未见过的目标域测试集(如表 5 所示)。该数据采集环境无精密实验平台,导致域偏移和个体差异问题更为严重(如图 11 所示),即使是原始数据也能清晰观测到个体差异现象,因此独立样本诊断任务的重要性进一步凸显。针对复合故障标签,采用 F1 分数评估未见过设备 F1 值(UDF1),对比实验结果如表 6 所示。

    表 5 电机数据集详情    

     
     

    图 11. 实验平台与个体差异

    表 6 电机故障诊断对比实验结果    

     
     

    与案例 1 的轴承实验结果类似,电机故障诊断的源设备 F1 值(SDF1)与未见过设备 F1 值(UDF1)差异显著:WDCNN 的 UDF1 达 82.53%,超过 TICNN 的 66.80%;SRDCNN、MSRFCNN 和 MS MCNN 的 UDF1 分别为 84.53%、76.97% 和 86.27%。需要注意的是,Transformer 基架构的整体泛化性能优于 CNN 架构;对于所提模型,两个设计模块分别带来 1.85% 和 2.31% 的泛化性能提升,进一步验证了该特殊设计的优越性与通用性。


     

    5.3 分析与讨论

    基于上述两个案例,从以下方面对智能故障诊断技术展开进一步讨论与总结:  

    • 1)个体差异的普遍性:轴承和整机级电机中均普遍存在个体差异,且会对数据驱动模型产生显著影响。如图 2 图 9 所示,即使是同类型设备,不同个体间的数据仍存在分布差异,这是传统基准此前忽视的潜在缺陷。强烈建议未来在数据集构建和算法开发中纳入个体泛化考量:构建新数据集时,应采集多台设备的数据而非单一设备;验证算法性能时,训练集与测试集不应来自同一设备个体。    

    • 2传统基准的局限性:传统基准无法有效验证模型对采集过程和测试环境引入的故障无关特征的屏蔽能力(如图 12 所示,所有 SDA 和 SDF1 均接近 100%)。在同一设备不同时段诊断中取得的高性能,无法保证模型的工程可靠性。部分基于传统诊断基准的算法模型可能仍具备有效性,但建议研究人员重新验证其独立样本泛化能力和实际工业场景适用性。    

    • 3复杂系统诊断的挑战性:整机级电机数据集的诊断任务比实验室轴承数据集更具挑战性。尽管电机数据集规模大于轴承数据集,但其性能指标仍低于实验室轴承 —— 这可能是由于现场信号受更多干扰因素影响,且需处理更复杂的复合故障。    

    • 4个体差异的缓解思路:增加同类设备数量可弱化个体差异,或通过设计多试验台训练数据主动掩盖个体差异;但在故障诊断领域,由于有效样本个体稀缺且标注成本高昂,此类增加数据样本的方案目前不具备实用性。    

    • 5Transformer 模型的潜力:Transformer 基模型是极具前景的网络架构,其优异的全局特征感知能力和注意力机制,非常适合处理个体差异引发的全局特征偏移和无关干扰。

    图 12. 源设备 F1 值(SDF1)与未见过设备 F1 值(UDF1)的对比

    6 总结

    本文探究了一种此前被忽视的、由动态数据采集过程、复杂测试环境以及固有装配差异引发的域偏移现象。即便对于状态类型相同的设备,不同个体间的数据分布仍存在差异。该现象在复杂整机级设备的故障诊断任务中表现得更为显著,使得此类诊断任务相比实验室轴承故障诊断具备更高的挑战性。仅在源设备个体上验证得到的高性能结果,无法保障模型的工程可靠性。  

    为此,本文提出了一种基于独立样本的个体泛化框架,并设计了一款改进型轻量化 Transformer 作为可行的解决方案,该模型融合了卷积神经网络(CNN)与 Transformer 架构的优势。研究还设计并嵌入了即插即用型模块,以实现多尺度信息交互与特征剪枝。  

    在公开轴承数据集与实际生产线工业电机数据集上开展的全面实验表明,独立样本评估基准具备重要的实践价值,同时也验证了所提模型的优越性。  

    在未来的研究工作中,我们将针对多台设备个体(而非单一设备)采集综合性数据集,且所有算法验证均会采用个体泛化框架。我们将进一步探究数据集规模对个体差异的影响,并致力于解决复杂整机级系统的复合故障诊断难题,尤其是未知复合故障组合的诊断任务。此外,个体差异的作用机制与影响因素也将是未来的研究方向,相关研究成果将有助于研发更可靠、更先进的特征提取模型


     

    开源代码获取:

    点击文章左下角阅读原文即可获取。

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

    来源:故障诊断与python学习
    ACT振动UDF旋转机械通用海洋理论电机机器人多尺度数字孪生控制DAP
    著作权归作者所有,欢迎分享,未经许可,不得转载
    首次发布时间:2025-12-27
    最近编辑:8月前
    故障诊断与python学习
    硕士 签名征集中
    获赞 87粉丝 152文章 327课程 0
    点赞
    收藏
    作者推荐

    寿命预测 SCI一区论文学习 |基于非参数退化模型和不完全数据的锂离子电池剩余使用寿命预测

    本期分享发表于一区Top期刊RESS的论文,基于非参数退化模型和不完全数据的锂离子电池剩余使用寿命预测。针对实际应用中监测数据存在的数据稀疏或碎片化问题,文章提出了一种非参数化退化建模方法,通过函数主成分分析(FPCA)技术,自动构建由退化数据驱动的状态依赖型退化模型,将RUL预测问题转化为迭代优化问题。通过丰田-MIT-斯坦福电池实验数据验证表明提出的方法能提供更准确可靠的预测结果。论文基本信息论文题目:Remaining useful life prediction of lithium-ion battery with nonparametric degradation modeling and incomplete data论文期刊: Reliability Engineering and System Safety论文时间:2025年作者: Naipeng Li (a), Mingyang Wang (a), Yaguo Lei (a), Bin Yang (a), Xiang Li (a), Xiaosheng Si (b)机构:a: Key Laboratory of Education Ministry for Modern Design and Rotor-Bearing System, Xi’an Jiaotong University, Xi’an, Shaanxi 710049, China;b: Zhijian Laboratory, Rocket Force University of Engineering, Xi’an 710025, China.论文下载:https://doi.org/10.1016/j.ress.2024.110721摘要准确可靠地预测电池剩余使用寿命(Remaining Useful Life, RUL)对保障电池安全运行至关重要。现有大多数RUL预测方法都依赖完整密集的监测数据,但在实际应用中,监测数据往往存在不完整性——即数据稀疏或碎片化,特别是在复杂严苛的工况下,这给电池RUL预测带来了巨大挑战。我们前期研究提出了一种非参数化退化建模方法,专门用于处理碎片化数据的RUL预测问题。其核心思路是通过函数主成分分析(Functional Principal Component Analysis, FPCA)技术,自动构建由退化数据驱动的状态依赖型退化模型。该方法将RUL预测问题转化为迭代优化问题。本文进一步完善了这一非参数化退化建模方法,并将其应用于锂离子电池在数据不完整情况下的RUL预测。主要改进包括:在梯度下降优化过程中引入新的伪Huber损失函数以提升RUL预测精度;提出结合贝叶斯更新的迭代优化算法来估计函数主成分(Functional Principal Component, FPC)分数分布,从而提供RUL的不确定性估计而非仅给出点估计值。通过丰田-MIT-斯坦福电池实验数据对所提方法进行验证,结果表明该方法相比其他预测方法能提供更准确可靠的预测结果。关键词: 锂离子电池(LIBs),剩余使用寿命预测, 不完整数据, 非参数退化建模目录1 引言 2 不完整数据下锂电池的RUL预测 2.1 数据预处理 2.2 离线训练 2.3 在线预测 3 实验验证 3.1 实验数据介绍 3.2 离线训练流程演示 3.3 在线预测过程演示 3.4 RUL预测结果的比较 4 结论 注:本文只选中原论文部分进行分享,若想进一步拜读,请下载原论文。 小编能力有限,如有翻译不恰之处,请多多指正~ 1 引言锂离子电池(Lithium-ion Batteries, LIBs)凭借其长寿命、低自放电、高能量密度和高输出电压等优势,在电动汽车(Electric Vehicles, EVs)、便携式电子设备和储能系统等领域得到广泛应用[1,2]。然而由于多种物理化学机制的影响,锂离子电池的性能会随着充放电循环次数的增加而逐渐下降[3]。若无法实时监测电池的劣化状态,可能导致设备性能下降,甚至引发爆炸、火灾等灾难性事故[4]。为保障LIBs的安全高效运行并制定及时维护计划,建立可靠且精准的电池性能衰减模型以预测其剩余使用寿命(Remaining Useful Life, RUL)至关重要[5]。 近年来,针对锂离子电池(LIBs)的循环寿命预测方法已发展出多种,主要分为两大类:基于模型的方法和数据驱动方法[6,7]。其中,基于模型的方法主要通过构建电化学模型[8,9]或等效电路模型[10,11]来预测电池的劣化行为。这类模型具有较强的可解释性,能够准确描述电池内部的电化学反应过程。但其存在明显局限:首先需要预先掌握电池工作原理,其次由于实际应用中存在非线性电化学过程、多重劣化机制及复杂运行条件,要构建全面精准的模型极具挑战性[12]。数据驱动方法通过机器学习或统计分析技术,无需预先掌握作用机制即可从历史数据中构建降解模型。支持向量机(Support Vector Machine, SVM)[13,14]和深度神经网络(Deep Neural Networ, DNN)[15–17]等机器学习技术具有高度灵活性,在RUL预测中表现良好。例如,张等人[17]提出了一种深度学习方法,用于预测电动飞机所用锂离子电池的荷电状态(State of Charge, SOC)。该方法将AI技术整合到电动飞机电池的预测与健康管理中,标志着电动航空智能化发展迈出重要一步。但大多数机器学习方法仅能提供RUL的点估计值,无法描述预测结果的不确定性。为此,学界已开发出多种统计方法,例如维纳过程[1,18]、伽马过程[19]、逆高斯过程[20]等。这些方法属于参数化建模技术,需要选择合适的函数形式来描述退化行为[21]。然而,退化过程本身具有动态性和复杂性,这给精确匹配函数的选择带来了重大挑战[22]。为避免手动构建模型,非参数建模技术因其灵活性优势备受关注。这类技术能直接从历史数据中自动构建构建模型。主要非参数劣化建模方法包括高斯过程(Gaussian Process, GP)回归、样条函数回归和函数主成分分析(Functional Principal Component Analysis, FPCA)等。高斯过程回归将劣化过程描述为遵循联合高斯分布的随机变量集 合[23,24]。样条函数回归通过组合常用正弦基函数和B样条基函数等基函数来构建劣化曲线[25]。 FPCA 是主成分分析处理函数数据的改进版本[26-29],其核心思想是将退化轨迹分解为均值函数、可分解为正交函数主成分(Functional Principal Components, FPCs)及其对应 FPC 分数的偏差函数,以及噪声项的叠加[30]。目前常用的非参数化电池性能衰减建模方法需要对反映电池性能的健康指标进行完整密集的监测。但对于LIB而言,其直接健康指标——放电容量——必须通过完整的充放电循环才能获得,但在实际应用中往往难以实现[3]。因此,每个充放电循环都难以获取精确的放电容量数据,导致数据存在缺失[31]。这类缺失数据可分为两类:稀疏数据和碎片化数据。如图1(a) 所示,稀疏数据表现为在离散时间点间歇性观测到的少量数据点。由于LIB的设计寿命较长,其性能衰减过程通常较为缓慢,因此在漫长使用寿命期间持续监测需要耗费大量时间。为此,实际应用中常采用稀疏采样策略。此外,考虑到网络带宽和存储容量等限制因素,制造商通常会以稀疏格式将电池数据传输至工业云系统[32]。碎片化数据的显著特征在于,整个LIB生命周期中仅记录随机观测片段,如图1(b) 所示。这种现象主要源于传感器故障或数据传输中断。例如,电动汽车因成本限制采用带宽受限的控制器局域网(Controller Area Network, CAN)进行数据交换。高频数据传输可能导致网络拥塞和数据中断[33]。 图1 两种不同类型不完整数据的示例:(a) 稀疏数据,(b) 片段数据。在利用不完整数据构建LIB退化模型时,主要障碍之一是无法从历史数据中确定退化过程的函数形式,这会显著影响参数化建模方法的预测准确性[34]。此外,常用的非参数化方法通过将健康状态或其代理指标作为时间函数来建立退化模型。但由于电池整个生命周期内数据都是随机采样,因此无法从不完整数据记录中确定退化时间尺度。这使得时间依赖型退化建模方法在此场景下难以应用。李等人[29]提出了一种基于两步轴变换的依赖状态退化建模方法来解决这些问题。该方法将剩余寿命预测问题转化为迭代优化问题,通过多维度数据变换实现精准建模。RUL预测值通过匹配观测数据与模型估计值进行优化。该方法采用平方误差和(Sum of Squared Errors, SSE)作为损失函数,但未考虑异常值对预测精度的影响。此外,该方法仅提供点估计值,无法量化预测结果的不确定性。针对上述两个问题,本文进一步改进了状态相关性劣化建模方法,并将其改进版本应用于LIBs在稀疏数据和碎片数据条件下的RUL预测。与现有大多数数据驱动方法不同,本文采用两步轴变换技术将退化数据映射到状态相关坐标系,有效解决了当前方法在处理不完整数据时难以准确捕捉劣化时间的难题。随后,通过基于 FPCA 的非参数建模技术构建了退化模型。相较于传统基于模型的方法,本方法无需依赖先验知识,既实现了模型构建的灵活性,又避免了预设函数形式的束缚,从而在实际应用中展现出更强的适应性和扩展性。模型构建完成后,通过结合贝叶斯更新的迭代优化算法,可实现在线预测电池剩余寿命(RUL),最终达成电池健康状态的动态监测与精准寿命预测。具体而言,本文基于文献[29]的核心改进点可归纳为以下内容:1)基于伪Huber损失函数,我们为迭代优化过程构建了一个新的目标函数。相较于SSE损失函数,这个函数对异常值具有更强的鲁棒性,能提供更精准的RUL预测结果。此外,与SSE相比,它能有效避免因SSE中较大误差导致的剧烈梯度波动。这种特性不仅提升了梯度优化过程的稳定性,还能防止模型过度受异常值影响,使其特别适用于数据噪声和异常值普遍存在的实际应用场景。2)我们开发了一种结合贝叶斯更新的迭代优化算法,用于估算 FPC 评分的分布。与现有RUL预测方法不同,该方法将预测问题转化为迭代优化问题,使预测过程更加直观。同时,该算法充分考虑了 FPC 评分的随机性特征。基于估算的 FPC 评分分布,我们采用自助法近似计算RUL的概率密度函数(Probability Density Function,PDF),这种量化预测不确定性的方法在维护决策中具有显著优势,而文献[29]仅能提供RUL的点估计值。3)本研究系统探讨了两种不完整数据类型——稀疏数据与片段数据,并对其进行了系统性分析。该方法充分考虑了不同类型数据不完整性的独特特征,即使在处理高度稀疏或片段数据时,仍能保持较高的预测精度。相比之下,我们先前的研究[29]仅专注于片段数据处理,未涉及稀疏性带来的挑战,而本文提出的方法在应对各类不完整数据场景时更具全面性。通过实验验证,该方法在这两种场景下的有效性得到了充分证明。本文其余部分安排如下:第2节介绍针对数据不完整锂离子电池的剩余使用寿命预测方法;第3节采用丰田-MIT-斯坦福数据集验证该方法;第4节总结全文。2 不完整数据下的RUL预测在本节中,我们提出了一种基于 FPCA 的预测方法,用于处理数据不完整的LIBs。图2概述了所提出方法的流程图,该方法由三个过程组成:数据预处理、离线训练和在线预测。首先,采用两步轴变换技术对监测数据进行预处理,并将RUL重新表述为状态依赖函数。在离线训练过程中,基于历史数据使用 FPCA 技术构建降解模型。然后,在在线预测过程中,开发了一种迭代优化算法来估计 FPC 分数的后验分布。最后,通过自助法预测RUL的概率密度函数。每一步的具体细节将在下文中介绍。 图2 基于 FPCA 的预测方法流程图2.1 数据预处理假设我们已经收集了 个故障单元的放电容量观察数据。对于第 个单元,记录了放电容量序列 及其相关的时间序列 ,其中 表示样本总数。这些故障单元的放电性能已低于预设的故障判定标准。在多数情况下,该标准设定为放电容量降至初始值80%以下时触发的性能退化阈值[3]。通常,故障单元的寿命终止时间戳(记为 )可通过历史维护日志或事件记录追溯。1)时间变量转换在数据不完整的情况下,通常无法获取退化过程的初始时间,这给时间依赖性退化建模方法带来了相当大的障碍。为解决这一挑战,我们采用了一种时间变量转换方法。鉴于每个失效单元的 已知,可计算每个观测值对应的RUL,即 。由此,时间被转换成了剩余寿命RUL。2) 轴旋转在实际应用中,当LIB的放电能力低于故障阈值时,就需要进行更换[3];即,在故障阈值以下将不再获得任何进一步的观测数据。这类数据被称为截断退化数据[35]。然而, FPCA 技术无法直接应用于此类数据,因为 FPCA 需要所有单元共享相同的x轴刻度以避免观测数据的选择偏差。这一问题已在[36]中得到了彻底解决。为了使截断退化数据适应 FPCA 的要求,通过轴旋转将 轴和 轴互换,考虑到每个单元的健康状态都处于相同的状态刻度范围内 ,其中 是初始容量值, 是故障阈值。因此,预处理后的数据更适合使用 FPCA 进行分析,并将RUL重新定义为健康状态的函数。鉴于预测的主要目标是根据健康阶段的退化水平来预测RUL,将RUL表示为状态水平的函数对于预测任务来说更为直接。2.2 离线训练1)基于 FPCA 的LIBs退化模型按照上述数据预处理策略,对于 的单元 ,其RUL 表达式如下: 其中 是所有单元退化数据的潜在平均趋势, 表示单元 与平均趋势的偏差,是均值为零,协方差 的随机变量。 是一个均值为零、方差为 的噪声项,假设其与 独立。值得注意的是,上述状态依赖模型仅适用于退化曲线呈现单调行为的情况。然而,在存在噪声数据时,这一假设可能不成立。为了确保单调性,可以在数据预处理前使用单调基函数对噪声数据进行平滑处理。根据Karhunen-Lo‘eve定理[26], 可展开为: 其中 对于 是相应的特征函数, 是有序的非负特征值。鉴于特征函数正交且构成域 内的函数基, 可以如下投影到函数基: 其中 对于 为 FPC 分数,具有以下性质 。在实际应用中,只有前几个FPC定义了偏差变化的最重要部分。因此,我们在第 个分量处截断等式(3)。等式(1)进一步重新表述如下: 2) 模型各成分的估计均值函数 根据可用的不完整数据使用局部线性平滑器[26]进行估计。 通过最小化以下项来进行估计: 其中 是一个带宽为 的核函数,通常选择高斯核函数。然后, 的估计值 。令 。已知 和 相互独立,因此协方差函数具有以下关系: 其中 在 时为1,否则为0。因此,只有当 时,选择 进行协方差面平滑处理,因为对角元素包含了测量噪声的方差。 的估计通过最小化以下项来定义: 其中 是带宽为 的对偶核函数。因此, 的估计为 。方差函数 也可以通过使用局部线性平滑器来估计: 其中 是 的对角线曲线。 的估计值为 。方差 可近似为:当时,计算公式为: ;否则为0。通过求解以下方程可估算特征函数及相应特征值: 在实际应用中,我们建议选用前 个函数主成分,其累计贡献率需超过98%。2.3 在线预测1) FPC分数更新为了模拟在役机组的退化过程,准确估计 FPC 分数至关重要。理论上,FPC 分数可以通过 求解。然而,不完整数据不足以准确估计 FPC 分数。为解决这一问题,开发了一种结合贝叶斯更新的迭代优化算法来更新 FPC 分数。我们首先忽略 FPC 分数的随机性,将其视为一系列 确定的值。根据等式(4)可知,RUL 服从正态分布 ,其中 , 是观测时间戳。 , ,其中 。 令 , FPC 分数的估计为: 假设我们有一个在役单元,其放电容量观测值不完全,记为 和时间戳 。由于在役单元的放电容量仅部分退化,因此未记录到失效时间 。因此,真实的剩余使用寿命 RUL 无法获得,这实际上是我们预测任务的主要关注点。显然, 对于 总是成立的。因此,每个观察点的RUL与最新时间戳的RUL直接对应: 其中 表示最新时间戳的RUL,该变量需进行预测。令 ,并将其代入等式(11)。 FPC 分数的估计值可作为依赖于 的函数确定。因此, FPC 分数更新问题转化为一个优化问题,即寻找最优的 以实现 与估计的 之间的最接近匹配。这里, 表示实际观察到的退化轨迹,而 RUL 预测的任务集中在确定其在坐标系中的垂直位置。 是通过计算 FPC 分数获得的估计轨迹。迭代优化过程专注于在坐标系内不断调整估计曲线相对于观察曲线的垂直位置,以最小化两者之间的误差。这种方法能够精确定位 ,从而有助于准确估计 RUL。然而,真实数据中常存在噪声和异常值,而SSE损失函数对这些异常值一视同仁,导致平方运算使误差贡献被夸大,从而可能扭曲预测精度。相比之下,伪Huber损失函数[37]提供了一种更稳健的处理异常值的方法,并在多个重要方面优于SSE:<1> 从物理特性来看,伪Huber损失函数在保持小误差拟合精度的同时,对大误差的敏感度较低。这种对异常值的鲁棒性在降解数据分析中至关重要,因为降解数据中常存在噪声或极端值,而伪Huber损失函数能更准确地捕捉潜在的退化趋势。<2> 数学上,伪Huber损失函数定义为: ,其中 表示误差, 是控制平滑度的超参数。对于接近零的误差,伪Huber损失函数的表现类似于均方误差(MSE)损失;而对于较大的误差,则接近均绝对误差(MAE)损失。这种平衡使其非常适合具有显著噪声的复杂数据集,该损失函数既能精确拟合微小偏差,又能在面对显著偏差或异常值时保持稳健性。在此,我们基于伪Huber损失函数构建目标函数。该函数对异常值具有鲁棒性,并且在接近零时平滑,这有助于减少优化过程中的振荡,特别是在梯度下降算法中。设 。目标函数的表述如下: 其中 是控制目标值周围损失函数平滑度的超参数。RUL的先验估计可通过梯度下降算法确定。迭代方案如下所示: 其中, , 。 迭代优化过程的初始值 被设置为零。当 或迭代步数时,迭代停止。这里, 表示收敛精度, 表示最大迭代步数。目标函数通过迭代最小化得到最优的 ,该值可用于等式(11)中推导 FPC 分数的估计。为了量化预测结果的不确定性,需要考虑 FPC 分数的随机性。假设 FPC 分数服从正态分布。 FPC 分数的先验分布 为 ,其中 是一个 的向量,元素为零, 。给定第 个单位的观测值,可以使用以下贝叶斯更新理论来更新先验分布: 是 关于 的条件分布,其服从正态分布 。 的后验分布为正态分布 ,其中 通过上述迭代优化步骤输入估计的 可获得后验分布。2)RUL估计的PDF在实际应用中,RUL 的点估计无法为锂离子电池的预测性维护提供可靠且科学的指导。相比之下,RUL 的概率密度函数在预测任务中起着关键作用,因为它可以量化预测中的不确定性。本文中,我们提供了一种基于 FPC 评分分布的自助法 [34] 来估计 RUL 的概率密度函数。具体步骤如下:首先,从上述提供的 FPC 分数分布中抽取样本 , 。然后,为每个 FPC 分数样本生成一个模拟的退化轨迹,具体如下: 最后,获得了一系列采样的RUL值 。由于RUL的概率密度函数遵循正态分布,可以使用基于这些采样值的最大似然估计算法计算该分布的参数,即均值和方差。因此,可以得到RUL的概率密度函数。此外,RUL的 置信区间(Confidence Interval, CI)可以近似为 ,其中L 和 ,使用经验的 和 分位数。3 实验验证3.1 实验数据介绍本节中,我们使用了一个基准实验数据集来验证所提出的方法。该数据集由丰田、麻省理工学院和斯坦福大学提供,包含标称容量为1.1Ah的磷酸铁锂(LFP)电池的退化数据[38]。数据集根据不同的快速充电协议分为三个批次。具体来说,我们选择了标记为“2018-04-12”的批次来验证我们提出的方法。由于实验过程中处理不当,部分电池提前终止。因此,我们排除了异常电池数据,保留了能够准确反映实际情况的数据。总共使用了25个电池数据。 在这25个电池单元中,随机选取20个单元构建退化模型。其余单元作为在役单元,用于评估我们方法的预测准确性。为了评估我们方法在不同不完整数据场景下的性能,我们模拟了包括稀疏数据场景和片段数据场景在内的不同不完整数据组。1)在稀疏场景中,我们以四种缺失率(75%、80%、85%、90%)对每个单元进行稀疏采样观察。例如,90%表示从每个单元中随机删除90%的观测值,剩余10%的观测值可用于模型训练。2)在片段数据场景中,以不同缺失率(10%、20%、30%、40%)删除连续数据段。对于片段数据,剩余部分由一至两个连续数据段组成。稀疏数据与碎片化数据的代表性样本分别如图3(a) 和(b)所示。对于在役设备,放电容量的退化起始值为区间[1.04 Ah,1.06 Ah]内的随机值。 图3 训练单元的电池退化数据:(a) 缺失率为90%的稀疏数据,(b) 缺失率为40%的片段数据3.2 离线训练流程演示为了构建状态依赖的退化模型,原始数据通过第2.1节中描述的两步轴变换进行预处理。图4显示了预处理后的数据,其中 轴表示健康状态, 轴表示RUL。可以使用预处理后的数据推导出均值函数、协方差函数 和方差函数 。然后,可以求解 的特征函数和特征值。我们以缺失率为40%的片段数据为例。协方差函数 及其对应的方差函数 的估计结果如图5(a) 所示。通过观察可以发现,方差曲线的显著部分超出了协方差曲面的对角线。这表明测量噪声导致的方差与协方差存在本质区别——通过剔除估计值中的对角线元素即可区分。图5(b) 展示了前两个特征函数(即函数主成分,FPCs)及其对应的占比。FPC1呈现渐进式衰减趋势,可解释95.77%的变分,成为方差的主要来源。根据预设标准,我们选取前两个FPCs,其合计解释了超过98%的变分。 图4 预处理退化数据:(a) 缺失率为90%的稀疏数据,(b) 缺失率为40%的片段数据 图5 FPCA 的估计值:(a)协方差和方差函数,(b)前两个 FPC 函数3.3 在线预测过程演示模型构建完成后,根据在役单元的观察数据获得RUL。使用第2.3节中开发的迭代优化算法估计 FPC 分数的后验分布。基于首次预测时间(First Prediction Times, FPT)[1,39]的 FPC 分数更新分布预测RUL概率密度函数。为了便于统计分析预测结果,当每个单元的实际RUL达到200个周期时,开始进行RUL预测。RUL及其随退化过程的预测将持续到 ,间隔为20个周期。图6展示了第一个测试单元在不同不完整度数据情景下的RUL概率密度函数及其期望值。具体而言,我们以丢失率为90%的稀疏数据情况和丢失率为40%的片段数据情况为例。从结果可以看出,随着有更多的退化数据可用,预测的RUL的期望值逐渐接近实际的RUL。同时,RUL的分布随着操作时间的推移逐渐收敛,这表明结果的不确定性变小,精度变高。 图6 测试单元1在不完整数据场景下的RUL预测:(a) 缺失率为90%的稀疏数据,(b) 缺失率为40%的片段数据3.4 RUL预测结果的比较为了展示所提出方法的优越性,本节选择了四个比较方法。方法1也是一种依赖状态的非参数方法,通过两步轴变换实现[29]。但它忽略了 FPC 分数的随机性,只能提供RUL的点估计值。此外,方法1的目标函数基于SSE损失函数构建,容易受到异常值的影响,有时会导致优化算法收敛速度变慢。方法2是一种基于 FPCA 的非参数方法,同样将退化过程建模为依赖状态的模型[36]。但其时间尺度仍然是退化时间,没有考虑在数据不完整的情况下降解过程初始时间的损失。方法3是一种最先进的参数化方法[40],考虑了降解过程的单元间变异性。通过计算预测RUL的绝对误差(Absolute Error, AE)来评估预测准确性。我们提出的方法、方法1、方法2和方法3的结果分别展示在图7和图8中。 图7 采用四种不同方法在稀疏数据缺失率为90%的情况下对RUL的预测误差:(a)所提出方法,(b)方法1,(c)方法2,(d)方法3 图8 采用四种不同方法对缺失率为40%的片段数据案例中RUL的预测误差进行评估:(a) 提出的方法,(b) 方法1,(c) 方法2,(d) 方法3可以看出,我们提出的方法在四个方法中具有最小的预测误差中位数和方差。随着RUL的减少,我们方法的预测结果变得更加准确和稳定。方法1忽略了 FPC 分数的随机性,只能提供RUL的点估计值。方法2的一个缺点是它仍然使用原始时间尺度来描述退化过程,由于缺乏退化过程的初始时间,观察时间并不能准确代表退化时间,导致退化建模过程中的偏差。方法3由于函数形式估计不准确和缺乏退化时间,表现最差。为了进一步分析缺失率对RUL预测准确性的影响,我们使用了两个指标,即MAE和均方根误差(Root Mean Squared Error, RMSE)[41]。图9展示了在不同缺失率下五个测试单元获得的MAE和 RMSE 。比较分析表明,我们提出的方法在各种不完整数据场景中始终优于其他方法,展现出更高的预测准确性和稳定性。 图9 四种方法在不同场景下的预测结果:(a)稀疏情况下的MAE值,(b)片段数据情况下的MAE值,(c)稀疏情况下的 RMSE 值,(b)片段数据情况下的 RMSE 值4 结论为了处理不完整的数据,本文提出了一种基于FPCA的预后方法来预测LIB的RUL。为了解决不完整数据场景下退化过程初始时间缺失的问题,开发了一种数据预处理技术,将RUL重新表述为状态相关函数。基于FPCA构建了一个非参数退化模型,可以在应用过程中消除对微调参数或经验预定义模型的需要。此外,本文开发了一种结合贝叶斯更新的迭代优化算法来更新FPC分数的分布。最后,开发了一种引导方法来量化预测结果的不确定性,这可以为预测性维护提供比点估计更可靠的信息。使用电池退化数据集评估所提出方法的预测结果。综合分析表明,相较于其他典型非参数建模方法和前沿参数化方法,本文提出的方法在精度方面具有显著优势。该方法在不同缺失率下保持高度稳定,能够捕捉跨单元的性能退化规律,特别适用于各类电池的RUL预测。其核心优势在于强大的适应性,可灵活处理多种退化行为模式。更重要的是,该方法能量化RUL的不确定性,有效提升预测性维护水平,具备良好的实际应用价值。尽管该方法效果显著,但其主要局限在于仍需通过完整的充放电循环来获取电池健康指标,这限制了其在实时应用中的实用性——因为完整充放电数据并非总能获取。未来研究可重点开发基于部分充放电数据的电池健康评估技术,从而提升该方法在动态现实场景中的适用性。 编辑:Tina校核:李正平、陈凯歌、曹希铭、赵学功、白亮、任超、海洋、陈宇航、陈莹洁、赵栓栓、王金、赵诚、肖鑫鑫、张优该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除来源:故障诊断与python学习

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