一种部分可观测信息下维修和备件供应联合优化方法

摘要

本发明公开了一种部分可观测信息下维修和备件供应联合优化方法，首先，基于系统的实时多维状态监测数据，建立向量自回归模型，并计算整体数据集的残差；其次，基于残差数据建立能反映系统退化过程的隐马尔科夫模型，并用期望最大化算法估计隐马尔科夫模型中未知的状态参数和观测参数；然后，基于贝叶斯定理实时更新系统在每个采样时刻的剩余使用寿命的概率密度函数；最后，将备件的交付时间考虑为随机变量而非传统的常量，以平均费用率最小化为目标，动态更新最优的部件更换时间和备件订购时间。本发明可充分利用系统的多维状态监测信息和预测信息，指导维修工程师根据备件的实际需要灵活的调整相应的备件交付时间的分布，确定最优的维修决策。

说明书

技术领域

本发明属于机械系统的维修工程领域，具体涉及一种部分可观测信息下维修和备件供应联合优化方法。

背景技术

不可修机械系统的维修和备件优化问题是实施维修活动的关键环节。在现有的维修模型中，广泛使用基于历史失效统计数据的可靠度分布用于更换模型和备件库存模型，解决诸如在满足经济性指标的情况下部件的最优更换时间和备件订购时间。此外，还需要优化备件的订购数量以及库存量，在满足备件可用度最大化的同时使得储存费用最小化。

传统的方法是使用基于部件的历史失效时间来确定寿命的可靠度函数，其反应了部件总体的统计特性，而未考虑各部件之间潜在的物理失效过程。在统计寿命分布的基础上，进行部件更换和备件订购的序贯优化决策。然而，该方法并未考虑部件个体的差异性，也无法描述单个个体的寿命分布函数。对于新研设备及大型或昂贵的关键设备，大多缺乏历史失效数据，因此也难以获得寿命的分布函数。近年来，先进的状态监测技术使得个体的剩余寿命预测成为可能，为基于预测信息的维修决策奠定了基础。

传统的维修和备件决策模型常常存在一些不够合理的假设，如假设部件的寿命分布函数已知、备件总是可用的或不考虑备件的交付时间等。虽有文献研究了基于寿命预测信息的部件更换和备件订购序贯优化决策，但其不是全局最优解，在工程实践中难以实施。

发明内容

发明目的：本发明提供一种部分可观测信息下维修和备件供应联合优化方法，可充分利用系统的状态监测信息制定维系计划，指导维修工程师根据实际需要灵活的调整备件交付时间的分布。

技术方案：本发明所述的一种部分可观测信息下维修和备件供应联合优化方法，包括以下步骤：

(1)基于被监测的机械系统的实时多维状态监测数据，选取健康数据部分建立向量自回归模型，并计算整体数据集的残差，使得预处理的数据满足正态性和独立性；

(2)基于残差数据建立反映系统退化过程的隐马尔科夫模型，并用期望最大化算法估计隐马尔科夫模型中未知的状态参数

(3)基于贝叶斯定理实时更新系统在每个采样时刻t

(4)将备件的交付时间考虑为随机变量而非传统的常量，以平均费用率最小化为目标，动态更新最优的部件更换时间和备件订购时间。

进一步地，所述步骤(1)实现过程如下：

对于传感器采集的多维状态监测数据，其表示为

其中，ε

利用估计的VAR模型参数

其中，

进一步地，所述步骤(2)实现过程如下：

建立能反应被监测的机械系统退化的3状态，包括状态0、状态1和状态2 的隐马尔科夫模型，其中状态0代表健康状态，状态1代表不健康状态，状态2 代表失效状态；使用期望最大化算法估计隐马尔科夫模型中的状态参数

Π

进一步地，所述步骤(3)实现过程如下：

由贝叶斯定理，在每个采样点，后验概率Π

在采样时刻t

概率密度函数为：

在采样时刻t

其中，E(CC)为一个寿命周期内的期望总费用，E(CL)为一个寿命周期的期望时间长度，使得一个周期的平均费用率最小化的值即为对应的最优更换时间

进一步地，所述步骤(4)的实现过程如下：

在采样时刻t

五种情况下的期望备件短缺时间分别用ES1、ES2、ES3、ES4和ES5表示，各条件下期望备件短缺时间的计算方法分别如下：

ES3＝0 (14)

ES5＝0 (16)

五种情况下的期望备件短缺时间分别用EH1、EH2、EH3、EH4和EH5表示，各条件下期望备件持有时间的计算方法分别如下：

EH1＝0 (17)

EH2＝0 (18)

EH4＝0 (20)

一个寿命周期内的期望备件短缺时间ES为：

ES＝ES1+ES2+ES3+ES4+ES5 (22)

一个寿命周期内的期望备件持有时间EH为：

EH＝EH1+EH2+EH3+EH4+EH5 (23)

那么，一个寿命周期内的期望总费用E(CC)为：

一个寿命周期内的期望时间长度E(CL)为：

有益效果：与现有技术相比，本发明的有益效果：代替传统的基于失效时间分布的优化决策模型，本发明可充分利用系统的状态监测信息制定维系计划，指导维修工程师根据实际需要灵活的调整备件交付时间的分布，进而制定符合实际工程应用的最优维修决策。

附图说明

图1为本发明的流程图。

具体实施方式

下面结合附图对本发明做进一步详细说明。

如图1所示，本发明提供一种部分可观测信息下维修和备件供应联合优化方法，具体包括以下步骤：

步骤1：基于被监测的机械系统的实时多维状态监测数据，选取健康数据部分建立向量自回归模型，并计算整体数据集的残差，使得预处理的数据满足正态性和独立性。

对于传感器采集的多维状态监测数据，其表示为

其中ε

利用估计的VAR模型参数

其中，

步骤2：基于残差数据建立反映系统退化过程的隐马尔科夫模型，并用期望最大化算法估计隐马尔科夫模型中未知的状态参数

基于整体残差数据，假设系统的退化过程符合3状态(状态空间S＝{0，1，2}) 的隐马尔科夫模型(HMM)：健康状态(状态0)、警告状态或不健康状态(状态1)和失效状态(状态2)。状态0和状态1是不可观测的状态，即状态是隐藏的。只有状态2可直接被观测。假设系统总是起始于健康状态，即P(X

其中，λ

用ξ＝inf{t≥0：X

其中，

假设收集到的N组状态监测失效数据，用F

E-step：计算伪似然函数：

其中

M-step：选取Λ

每步更新的参数Λ

最大化每步更新的期望值，可得到每步更新的状态参数和观测参数的估计值，即各参数导数为0的点。每步更新的

每步更新的

其中：

在完成状态参数

Π

由贝叶斯定理，在每个采样点，后验概率Π

其中初始值Π

因此，在采样时刻t

概率密度函数为：

步骤3：将备件的交付时间考虑为随机变量而非传统的常量，以平均费用率最小化为目标，动态更新最优的部件更换时间和备件订购时间。

在采样时刻t

1)在备件订购时间点之前部件发生失效；

2)在备件已经发出订购和备件还未送达的时间点之间部件发生失效；

3)在备件已送达和最优更换时间点之间部件发生失效。

4)在备件已经发出订购要求，但在最优更换时间点之前备件并未送达。

5)备件已送达，在最优更换时间之后部件发生失效。

五种情况下的期望备件短缺时间分别用ES1、ES2、ES3、ES4和ES5表示，各条件下期望备件短缺时间的计算方法分别如下：

ES3＝0 (18)

ES5＝0 (20)

其中，l为交付时间；f(l)为交付时间的概率密度函数；

五种情况下的期望备件短缺时间分别用EH1、EH2、EH3、EH4和EH5表示，各条件下期望备件持有时间的计算方法分别如下：

EH1＝0 (21)

EH2＝0 (22)

EH4＝0 (24)

由上述分析可得，一个寿命周期内的期望备件短缺时间ES为：

ES＝ES1+ES2+ES3+ES4+ES5 (26)

一个寿命周期内的期望备件持有时间EH为：

EH＝EH1+EH2+EH3+EH4+EH5 (27)

那么，一个寿命周期内的期望总费用E(CC)为：

一个寿命周期内的期望时间长度E(CL)为：

综合以上分析，基于预测信息的更换和备件订购联合优化函数为：

上式中，使得一个周期的平均费用率最小化的值即为最优的更换时间

若优化得出的备件订购时间大于下次采样时间，那么不采取备件订购策略，继续进行下一次采样。若备件订购时间在下次采样之前，则根据优化的更换时间和备件订购时间采取维修措施。

虽然本发明已以较佳实施例揭露如上，然其并非用以限定本发明。本发明所属技术领域中具有通常知识者，在不脱离本发明的精神和范围内，当可作各种的更动与润饰。因此，本发明的保护范围当视权利要求书所界定者为准。