一种考虑状态依赖时滞的剩余寿命预测方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN112231925A(43)申请公布日2021.01.15(21)申请号202011159533.6(22)申请日2020.10.27(71)申请人山东科技大学地址266590山东省青岛市黄岛区前湾港路579号(72)发明人周东华席霄鹏纪洪泉钟麦英高明王建东(74)专利代理机构青岛智地领创专利代理有限公司37252代理人肖峰(51)Int.Cl.G06F30/20(2020.01)G06F111/08(2020.01)G06F111/10(2020.01)G06F119/04(2020.01)权利要求书2页说明书5页附图2页(54)发明名称一种考虑状态依赖时滞的剩余寿命预测方法(57)摘要本发明公开了一种考虑状态依赖时滞的剩余寿命预测方法，属于预测与健康管理领域。本发明包括以下步骤：输入一组退化数据，初始化退化模型的状态依赖时滞结构和参数；基于隐式欧拉格式下的离散化模型描述，对原始数据进行重采样；利用极大似然算法实现未知参数的估计；结合维纳过程变换理论降低首达时间分析的复杂度，并通过一步外推法模拟未来的退化轨迹，进而求得剩余寿命分布的近似解析解；最后，输出各个监测时刻下的剩余寿命分布。本方法适用于处理具有时变滞后效应的分形退化过程，主要应用在大型高炉炉壁、冷却壁侵蚀分析与维护。CN112231925ACN112231925A权利要求书1/2页1.一种考虑状态依赖时滞的剩余寿命预测方法，其特征在于预测工业设备剩余寿命，包括如下步骤：步骤1：输入一组长度为N的目标设备的退化数据{x1,x2,…,xN}，将初值置为零，并根据如下状态依赖时滞模型结构初始化最大时滞T：其中，t是监测时刻，η是漂移系数，ξ是非线性系数，σ是扩散系数，BH(t)是赫斯特指数为H的标准分形布朗运动，Ts(x(t))是时滞函数，且有x(t)∈[0,∞)；步骤2：在隐式欧拉格式下对状态依赖时滞模型进行离散化处理，离散化表示形式如下：其中，其中，k是离散监测时刻的序号，a是加权因子，ρ是重采样间隔，n是重采样时刻的序号，τ是原始采样间隔，y1:NΔ是重采样样本，a(·)是非线性漂移函数，εn和是高斯噪声，εn～N(0,ρ)，h(·)是分形扩散项在弱收敛条件下的尺度函数，I{·}是示性函数，H是赫斯特指数，Γ(·)是Gamma函数，s，r，m，c，J，χ和是用于求解h(·)的内部变量；步骤3：调用Nelder-Mead单纯形法搜索模型未知参数集Θ＝{η,ξ,σ,H}的极大似然解：2CN112231925A权利要求书2/2页步骤4：在tk时刻下，生成M条模拟退化轨迹基于维纳过程变换理论，剩余寿命的概率密度函数表示为：其中，lk是剩余寿命，其中，为失效阈值；步骤5：最终输出各个监测时刻下目标设备的剩余寿命估计结果值。2.根据权利要求1所述的一种考虑状态依赖时滞的剩余寿命预测方法，其特征在于，所述工业设备是面向机理复杂且存在状态依赖时滞与特定分形特征的工业设备。3.根据权利要求1所述的一种考虑状态依赖时滞的剩余寿命预测方法，其特征在于，所述工业设备包括但不限于高炉炉壁、冷却壁。3CN112231925A说明书1/5页一种考虑状态依赖时滞的剩余寿命预测方法技术领域[0001]本发明属于预测与健康管理领域，具体涉及一种考虑状态依赖时滞的剩余寿命预测方法。背景技术[0002]剩余寿命通常定义为退化过程距离首次达到失效阈值的剩余时间，可通过估计其数学期望或概率分布来优化现行的维护策略。对于炼铁、炼油等复杂工业过程，基于监测数据建立合理的退化模型则是预测系统剩余寿命的重要基础。[0003]近年来，分形布朗运动驱动的退化建模方法为解决非马尔可夫过程的剩余寿命预测问题提供了新的思路，得到了诸多专家学者的关注(Zhang,2017；Wang,2020；Song,2020)。在布朗运动的基础上，分形布朗运动及其扩展形式引入了长短期记忆效应，构建了历史退化状态与未来演变趋势的统计关联，相比传统的随机游走过程更适用于刻画大型高炉炉壁、冷却壁温度等性能变量的有偏扩散特征。[0004]然而，现有方法在进行退化分析时大多忽略了潜在的滞后效应。不同于全寿命周期中全局的长短期记忆效应，时滞则更侧重于反映近邻区间内退化状态的局部依赖关系。注意到仅有的一些时滞系统寿命预测研究均是采用人工智能的方法来弱化退化过程的马尔可夫特性，并未建立更为直观的时滞模型(Zhang,2015；Liu,2016；Rai,2017)。特别地，上述方法一般假定时滞定常，无法应对不稳定工况下时滞时变且与当前状态相关的情形。此问题可归结为一类状态依赖时滞问题，即需要根据实时数据动态调整退化速率变化规律的滞后程度。发明内容[0005]针对现有技术中存在的上述问