# Remaining Useful Life Estimation Based on Asynchronous Multisource Monitoring Information Fusion

Authors: Yanyan Hu; Shuai Qi; Xiaoling Xue; Kaixiang Peng
Journal: Journal of Control Science and Engineering
 (2017)

---

## Abstract

An asynchronous RUL fusion estimation algorithm is presented for the hidden degradation process with multiple asynchronous monitoring sensors based on multisource information fusion. Firstly, a state-space type model is established by modeling the stochastic degradation as a Wiener process and transforming asynchronous indirectly observations in the fusion period to the fusion time. The statistical characteristics of involved noises and their correlations are analyzed. Secondly, the estimate of the hidden degradation state is obtained by applying Kalman filtering with correlated noises to the established state-space model, where the synchronized observations are fused. Also, the unknown model parameters are recursively identified based on the Expectation-Maximization (EM) algorithm with the Generic Algorithm (GA) adopted to solve the maximization problem. Finally, the probability distribution of RUL is obtained using the fused degradation state estimation and the updated identification result of the model parameters. Simulation results show that the proposed fusion method has better performance than the RUL estimation with single sensor.


---

## Body

## 1. Introduction

Modern engineering systems are becoming large in scale, huge in investment, and more and more sophisticated in structure with the development of science and technology. As a result, once an accident occurs in these systems, it would cause tremendous damage and enormous loss of property [1]. As time goes on, degradation of equipment in systems will inevitably occur due to complex operating environment and many other reasons, which leads to the reduction of system reliability and safety [2, 3]. Consequently, it is of great importance to predict the RUL of the system components in order to correctly evaluate the health state of the system and make appropriate maintenance plans to ensure the safe operation of the system and increase economic returns.The complexity of the system itself and the running environment makes it hard to establish a mechanism model for RUL estimation. In recent years, the RUL estimation based on monitoring data attracts a lot of research attentions due to its wide application range and the ability to quantify the uncertainties of the estimation results [4, 5]. Si et al. in [6] reviewed the statistical data driven approaches for RUL estimation. The existing methods were classified into two categories: direct condition monitoring data based approaches and indirect condition monitoring data based approaches, which can be further divided into stochastic filtering based methods [7–10], covariance based hazard model methods [11, 12], Wiener-process-based methods [13, 14], Gamma process based methods [15, 16], and Markovian-based methods [17] and others.Although there are fruitful researches on RUL estimation based on stochastic modeling with monitoring data, most of works are specific to single sensor. However, the complexity of the system itself and the running environment makes the RUL prediction have large uncertainties. At the same time, the information from single source tends to be quite limited. In response to the above issues, multiple sensors are usually utilized to monitor the condition of the system in order to reduce the uncertainties of the system and improve the accuracy of the RUL estimation [3]. Wei et al. in [18] proposed a multisensor information based Remaining Useful Life estimation method with anticipated performance for a kind of stochastic linear degradation with Wiener properties. The degradation process as well as the RUL distribution with a random threshold was identified recursively by adopting the distributed fusion filtering and a two-stage parameter estimation approach. Meanwhile, the sensor selection problem was also discussed by quantitatively analyzing the benefits of the fusion manipulation. Dong and He in [19] presented a hidden semi-Markov-model based methodology for multisensor equipment health diagnosis and prognosis. Sensor fusion was implemented by adjusting the weight assigned to each sensor based on the discriminant function analysis. A new health indicator was constructed in [20] by fusing the mutual information of multiple features extracted from the vibration signal to predict the RUL of Machinery. Li and He in [21] proposed a methodology to predict the RUL of both wheels and trucks by fusing data from three types of detectors, including wheel impact load detector, machine vision systems, and optical geometry detectors. Although multisource information is used in the above works to improve the performance of RUL estimation, observations from multiple sensors were assumed to be collected at the same time, that is, time consistent as shown in Figure 1(a). However, in practical engineering, observations from multiple sensors are usually not aligned in time due to different sampling periods, distinct initial sampling time, and many other reasons [22, 23]. This results in that multisensor observations are not synchronous but asynchronous as shown in Figure 1(b) and poses new challenges to the RUL estimation.Figure 1
Observations from multiple sensors.
(a)
Synchronous
(b)
AsynchronousConsequently, motivated by the above discussions, the asynchronous RUL fusion estimation problem is studied in this paper. We assume that the hidden degradation which is modeled as Wiener process with unknown model parameters is observed by an arbitrary number of asynchronous sensors, whose sampling frequencies and initial sampling times are all arbitrary. The asynchronous sensor observations are firstly synchronized to the fusion time and then fused using the Kalman filtering technology to get the fused estimate of the latent degradation state with the correlations between various noises introduced by the synchronization process analyzed. The unknown model parameters are also recursively identified using the synchronized observations based on theEM algorithm with the maximization problem solved byGA. Finally, the RUL is derived based on the fused estimates of the latent degradation state and unknown model parameters and the simulation example is provided to demonstrate the feasibility and effectiveness of the proposed algorithm.The organization of this paper is as follows. The fusion estimation problem of RUL for stochastic degradation process with multiple asynchronous sensors is formulated in Section2. The proposed asynchronous RUL fusion estimation algorithm is derived in Section 3. Section 4 gives simulation results and conclusions are drawn in Section 5.


## 2. Problem Formulation

Wiener process is widely used to model the stochastic degradation of a system attributed to its good mathematical characteristics such as the infinite separability. In general, a linear Wiener-process-based degradation model can be represented as(1)Xt=ϕ+ηt+σBt,where X(t) is the random variable representing the degradation at time t, ϕ is the initial degradation state, η and σ are, respectively, the drift and diffusion coefficient, and B(t) is the standard Brownian motion reflecting the stochastic dynamics of the degradation process.We assume the degradation process (1) is observed by a number of M asynchronous sensors. Denote tk-1 and tk, respectively, as the previous and oncoming fusion time instants. Since these asynchronous sensors may have different sampling periods and the sampling intervals may even be nonuniform, more than one measurement could probably be obtained by a given sensor during the fusion interval (tk-1,tk]. Denote tkm as the time when the latest measurement ykm, which is the nearest to tk, is observed by sensor m, where ykm is the indirect observation of the degradation process described by(2)ykm=hmxtkm+vtkm,where xtkm is the realization of Xt at time tkm and vtkm~N0,rm is the observation noise.Consequently, the objective of this paper is to predict the distribution of RUL of degradation process (1) based on the cumulative measurements Yk={yk′i,k′=1,…,k,i=1,…,M}, where the RUL of a system at time tk is defined as(3)lk=inf⁡lk:Xtk+lk≥ω∣X0<ω,where ω is the predefined threshold level.


## 3. The Proposed RUL Fusion Estimation Algorithm

### 3.1. Fusion Estimation of the Degradation State

From the degradation process (1), we know that the transition of the degradation state satisfies(4)xtj=xti+ηtj-ti+wtj,ti,where w(tj,ti)=σ(B(tj)-B(ti))~N(0,σ2(tj-ti)). By substituting (4) into (2) with tj and ti in (4) replaced with tk and tkm, respectively, it follows that(5)ykm=hmxtk-hmηtk-tkm-hmwtk,tkm+vtkm=hmxtk-hmηtk-tkm+ξkm,where ξkm=v(tkm)-hmw(tk,tkm). Then, by defining yk=yk1,…,ykMT, ξk=ξk1,…,ξkMT, h=h1,…,hMT, and vk=[v(tk1),…,v(tkM)]T, the asynchronous indirect observations yk1,yk2,…,ykM collected from M asynchronous sensors in tk-1,tk can be regarded as an augmented measurement at tk with(6)yk=hxtk-ηTkh+ξk,where(7)Tk=tk-tk1⋱tk-tkM,ξk=vk+h1wtk,tk1⋮hMwtk,tkMand correspondingly(8)cov⁡ξk=cov⁡vk+cov⁡wkwith(9)cov⁡vk=rj1⋱rjM,cov⁡wk=h1h1σ2tk-tk1h1h2σ2tk-tk2⋯h1hMσ2tk-tkMh2h1σ2tk-tk2h2h2σ2tk-tk2⋯h2hMσ2tk-tkM⋮⋮⋮⋮hMh1σ2tk-tkMhMh2σ2tk-tkM⋯hMhMσ2tk-tkM.Obviously, we see that cov⁡ξk is a function of parameters {σ2,rm,hm}, and we denote it as cov⁡ξk=Σkσ2,rm,hm for simplicity.From (4), we can directly have(10)xtk=xtk-1+ηtk-tk-1+wtk,tk-1.It follows that(11)covwtk-tk-1,ξkm=covwtk-tk-1,-hmwtk-tkm=-hmσ2tk-tkm,(12)cov⁡wtk-tk-1,ξk=-σ2h1tk-tk1h2tk-tk2⋮hMtk-tkMT=-σ2TkhTThen, by applying the Kalman filtering with correlated noises to the degradation process (10) and the augmented observation equation (6) [24, 25], the optimal estimate of the degradation state x at time tk can be obtained in a recursive form as follows.Firstly, the predicted estimate and its error covariance are given by(13)x^k∣k-1=Extk∣Yk-1=x^k-1∣k-1+ηtk-tk-1,Pk∣k-1=var⁡xtk∣Yk-1=Pk-1+σ2tk-tk-1.Secondly, the updated estimate and its error covariance are got by fusing the asynchronous indirect observationsyk=yk1,…,ykMT in (tk-1,tk] as(14)x^k=x^k∣k-1+Pkh-σ2TkhPk∣k-1TMkσ2,hm,rm-1yk-hx^k∣k-1+ηTkh,Pk=11/Pk∣k-1+h-σ2Tkh/Pk∣k-1TMkσ2,hm,rm-1h-σ2Tkh/Pk∣k-1T,where x^k=Extk∣Yk, Pk=var⁡xtk∣Yk, and(15)Mkσ2,hm,rm=cov⁡ξk-cov⁡wtk-tk-1,ξkTcov⁡wtk-tk-1,ξkPk∣k-1.Remark 1.
Note from (8), (9), and (12) that we know that M is also a function of unknown parameters σ2,hm,rm.


### 3.2. Unknown Parameter Fusion Identification

As we know, in real applications, the only information we have is the indirect observationsYk={yk′i,k′=1,…,k,i=1,…,M} from the M asynchronous sensors, in that, to obtain the fused estimate x^k of degradation state according to (13) and (14), we need to estimate unknown parameters Θ={η,σ2,hm,rm} at first. In this paper, theEM algorithm is implemented to identify Θ based on the online asynchronous sensor observations.EM algorithm is a widely used iterative algorithm for maximum likelihood parameter estimation with unobserved latent variables. Each iteration ofEM algorithm involves two steps: the expectation step (E-step) and the maximization step (M-step). In theE-step, a function for the expectation of the log-likelihood is established based on the current estimate of the unknown parameter. In the M-step, the parameter estimate is updated through maximizing the expected function found in the E-step, and then the updated parameter estimate is utilized to determine the distribution of the latent variables in theE-step of next iteration. By this way,EM algorithm iterates between the two steps until convergence [14, 18].The joint log-likelihood function at timetk is defined as(16)LkΘ=ln⁡pyk,xk∣Θ=ln⁡pxk∣Θ+ln⁡pyk∣xk,Θ,where {yk,xk} is the complete data set. Based on the Gauss white noise assumptions, we have yk~Nhxk-ηTkh,Σkσ2,rm,hm and xk~Nxk-1+ηtk-tk-2,σ2tk-tk-1; further it follows that(17)LkΘ=ln⁡12πσ2tk-tk-1exp⁡-xk-xk-1-ηtk-tk-122σ2tk-tk-1+ln⁡2π-M/2Σkσ2,rm,hm-1/2exp⁡-12yk-hxk+ηTkhTΣkσ2,rm,hm-1yk-hxk+ηTkhT=-12ln⁡2πσ2tk-tk-1-xk-xk-1-ηtk-tk-122σ2tk-tk-1-M2ln⁡2π-12ln⁡Σkσ2,rm,hm-12yk-hxk+ηTkhTΣkσ2,rm,hm-1yk-hxk+ηTkh.As we said, theEM algorithm consists of two steps:E-step andM-step in each iteration. (i)
E-step(18)ϑΘ∣Θ^kl=EXtk∣Yk,Θ^klln⁡pyk,xk∣Θ∝-12ln⁡σ2-12σ2tk-tk-1xk2^l+xk-1∣k2^l-2xkxk-1∣k^l-2ηtk-tk-1x^kl-x^k-1∣kl+η2tk-tk-12-12ln⁡Σkσ2,rm,hm-12ykTΣkσ2,rm,hm-1yk-2ykTΣkσ2,rm,hm-1hx^kl+xk2^lhTΣkσ2,rm,hm-1h+ηhTTkΣkσ2,rm,hm-12yk-2hx^kl+ηTkh,where x^kl=Extk∣Yk,Θ^kl, xk2^(l)=Ex2tk∣Yk,Θ^kl=x^kl2+Pk(l), x^k-1∣k(l)=Extk-1∣Yk,Θ^kl, xk-1∣k2^(l)=Ex2tk-1∣Yk,Θ^kl=x^k-1∣kl2+Pk-1∣k(l), xkxk-1∣k^(l)=Extkxtk-1∣Yk,Θ^kl=Pk,k-1∣k(l)+x^k-1∣k(l)x^k(l), and(19)x^k-1∣kl=x^k-1l+Pk-1lPk∣k-1lx^kl-x^k∣k-1l,Pk-1∣kl=Pk-1l-Pk-1lPk∣k-1l2Pk∣k-1l-Pkl,Pk,k-1∣kl=PklPk-1lPk∣k-1l,where x^k∣k-1(l), Pk∣k-1(l), x^k(l), and Pkl can be calculated from (13) and (14) with Θ=η,σ2,hm,rm replaced with Θkl={ηk(l),σ2k(l),hmk(l),rmk(l)}, respectively [24].(ii)
M-step: at time tk, the estimate of Θ for the l+1thEM iteration is calculated by(20)Θ^kl+1=arg⁡maxΘ⁡ϑΘ∣Θ^kl.In order to get the value of Θ^kl+1,GA [26, 27] is adopted to solve the maximization problem of (20).GA is a heuristic algorithm inspired by the process of natural selection, which is widely used for solving optimization problems. The evolution ofGA is an iterative process. It starts with an initial population, usually generated randomly, and the fitness of each individual in the population is evaluated, where the fitness function is the objective function of the optimization problem being solved. The new generation, which is used for the next iteration, is then formed through bioinspired operators, such as selection, mutation, and crossover, based on the fitness. Generally, the algorithm terminates when either a maximum number of generations has been generated or a satisfactory fitness level has been reached for the population.The iteration of theEM algorithm is terminated if the convergence condition Θ^(l+1)-Θ^(l)≤ς is satisfied, where ς is a preset threshold.


### 3.3. RUL Fusion Estimation

From (3), we can see that the RUL is defined as the first passage time of stochastic degradation process Xt. As we know, the first passage time of Wiener process is inverse Gauss distribution, and thus the probability distribution function flk∣Yk of RUL can be obtained by(21)flk∣Yk=ω-x^kσ2+Pkη2πPk+σ2tk3exp⁡-ω-x^k-ηlk22Pk+σ2lk,where Elk∣Yk=ω-x^k/η and var⁡lk∣Yk=ω-x^kσ2+Pkη/η3.


## 3.1. Fusion Estimation of the Degradation State

From the degradation process (1), we know that the transition of the degradation state satisfies(4)xtj=xti+ηtj-ti+wtj,ti,where w(tj,ti)=σ(B(tj)-B(ti))~N(0,σ2(tj-ti)). By substituting (4) into (2) with tj and ti in (4) replaced with tk and tkm, respectively, it follows that(5)ykm=hmxtk-hmηtk-tkm-hmwtk,tkm+vtkm=hmxtk-hmηtk-tkm+ξkm,where ξkm=v(tkm)-hmw(tk,tkm). Then, by defining yk=yk1,…,ykMT, ξk=ξk1,…,ξkMT, h=h1,…,hMT, and vk=[v(tk1),…,v(tkM)]T, the asynchronous indirect observations yk1,yk2,…,ykM collected from M asynchronous sensors in tk-1,tk can be regarded as an augmented measurement at tk with(6)yk=hxtk-ηTkh+ξk,where(7)Tk=tk-tk1⋱tk-tkM,ξk=vk+h1wtk,tk1⋮hMwtk,tkMand correspondingly(8)cov⁡ξk=cov⁡vk+cov⁡wkwith(9)cov⁡vk=rj1⋱rjM,cov⁡wk=h1h1σ2tk-tk1h1h2σ2tk-tk2⋯h1hMσ2tk-tkMh2h1σ2tk-tk2h2h2σ2tk-tk2⋯h2hMσ2tk-tkM⋮⋮⋮⋮hMh1σ2tk-tkMhMh2σ2tk-tkM⋯hMhMσ2tk-tkM.Obviously, we see that cov⁡ξk is a function of parameters {σ2,rm,hm}, and we denote it as cov⁡ξk=Σkσ2,rm,hm for simplicity.From (4), we can directly have(10)xtk=xtk-1+ηtk-tk-1+wtk,tk-1.It follows that(11)covwtk-tk-1,ξkm=covwtk-tk-1,-hmwtk-tkm=-hmσ2tk-tkm,(12)cov⁡wtk-tk-1,ξk=-σ2h1tk-tk1h2tk-tk2⋮hMtk-tkMT=-σ2TkhTThen, by applying the Kalman filtering with correlated noises to the degradation process (10) and the augmented observation equation (6) [24, 25], the optimal estimate of the degradation state x at time tk can be obtained in a recursive form as follows.Firstly, the predicted estimate and its error covariance are given by(13)x^k∣k-1=Extk∣Yk-1=x^k-1∣k-1+ηtk-tk-1,Pk∣k-1=var⁡xtk∣Yk-1=Pk-1+σ2tk-tk-1.Secondly, the updated estimate and its error covariance are got by fusing the asynchronous indirect observationsyk=yk1,…,ykMT in (tk-1,tk] as(14)x^k=x^k∣k-1+Pkh-σ2TkhPk∣k-1TMkσ2,hm,rm-1yk-hx^k∣k-1+ηTkh,Pk=11/Pk∣k-1+h-σ2Tkh/Pk∣k-1TMkσ2,hm,rm-1h-σ2Tkh/Pk∣k-1T,where x^k=Extk∣Yk, Pk=var⁡xtk∣Yk, and(15)Mkσ2,hm,rm=cov⁡ξk-cov⁡wtk-tk-1,ξkTcov⁡wtk-tk-1,ξkPk∣k-1.Remark 1.
Note from (8), (9), and (12) that we know that M is also a function of unknown parameters σ2,hm,rm.


## 3.2. Unknown Parameter Fusion Identification

As we know, in real applications, the only information we have is the indirect observationsYk={yk′i,k′=1,…,k,i=1,…,M} from the M asynchronous sensors, in that, to obtain the fused estimate x^k of degradation state according to (13) and (14), we need to estimate unknown parameters Θ={η,σ2,hm,rm} at first. In this paper, theEM algorithm is implemented to identify Θ based on the online asynchronous sensor observations.EM algorithm is a widely used iterative algorithm for maximum likelihood parameter estimation with unobserved latent variables. Each iteration ofEM algorithm involves two steps: the expectation step (E-step) and the maximization step (M-step). In theE-step, a function for the expectation of the log-likelihood is established based on the current estimate of the unknown parameter. In the M-step, the parameter estimate is updated through maximizing the expected function found in the E-step, and then the updated parameter estimate is utilized to determine the distribution of the latent variables in theE-step of next iteration. By this way,EM algorithm iterates between the two steps until convergence [14, 18].The joint log-likelihood function at timetk is defined as(16)LkΘ=ln⁡pyk,xk∣Θ=ln⁡pxk∣Θ+ln⁡pyk∣xk,Θ,where {yk,xk} is the complete data set. Based on the Gauss white noise assumptions, we have yk~Nhxk-ηTkh,Σkσ2,rm,hm and xk~Nxk-1+ηtk-tk-2,σ2tk-tk-1; further it follows that(17)LkΘ=ln⁡12πσ2tk-tk-1exp⁡-xk-xk-1-ηtk-tk-122σ2tk-tk-1+ln⁡2π-M/2Σkσ2,rm,hm-1/2exp⁡-12yk-hxk+ηTkhTΣkσ2,rm,hm-1yk-hxk+ηTkhT=-12ln⁡2πσ2tk-tk-1-xk-xk-1-ηtk-tk-122σ2tk-tk-1-M2ln⁡2π-12ln⁡Σkσ2,rm,hm-12yk-hxk+ηTkhTΣkσ2,rm,hm-1yk-hxk+ηTkh.As we said, theEM algorithm consists of two steps:E-step andM-step in each iteration. (i)
E-step(18)ϑΘ∣Θ^kl=EXtk∣Yk,Θ^klln⁡pyk,xk∣Θ∝-12ln⁡σ2-12σ2tk-tk-1xk2^l+xk-1∣k2^l-2xkxk-1∣k^l-2ηtk-tk-1x^kl-x^k-1∣kl+η2tk-tk-12-12ln⁡Σkσ2,rm,hm-12ykTΣkσ2,rm,hm-1yk-2ykTΣkσ2,rm,hm-1hx^kl+xk2^lhTΣkσ2,rm,hm-1h+ηhTTkΣkσ2,rm,hm-12yk-2hx^kl+ηTkh,where x^kl=Extk∣Yk,Θ^kl, xk2^(l)=Ex2tk∣Yk,Θ^kl=x^kl2+Pk(l), x^k-1∣k(l)=Extk-1∣Yk,Θ^kl, xk-1∣k2^(l)=Ex2tk-1∣Yk,Θ^kl=x^k-1∣kl2+Pk-1∣k(l), xkxk-1∣k^(l)=Extkxtk-1∣Yk,Θ^kl=Pk,k-1∣k(l)+x^k-1∣k(l)x^k(l), and(19)x^k-1∣kl=x^k-1l+Pk-1lPk∣k-1lx^kl-x^k∣k-1l,Pk-1∣kl=Pk-1l-Pk-1lPk∣k-1l2Pk∣k-1l-Pkl,Pk,k-1∣kl=PklPk-1lPk∣k-1l,where x^k∣k-1(l), Pk∣k-1(l), x^k(l), and Pkl can be calculated from (13) and (14) with Θ=η,σ2,hm,rm replaced with Θkl={ηk(l),σ2k(l),hmk(l),rmk(l)}, respectively [24].(ii)
M-step: at time tk, the estimate of Θ for the l+1thEM iteration is calculated by(20)Θ^kl+1=arg⁡maxΘ⁡ϑΘ∣Θ^kl.In order to get the value of Θ^kl+1,GA [26, 27] is adopted to solve the maximization problem of (20).GA is a heuristic algorithm inspired by the process of natural selection, which is widely used for solving optimization problems. The evolution ofGA is an iterative process. It starts with an initial population, usually generated randomly, and the fitness of each individual in the population is evaluated, where the fitness function is the objective function of the optimization problem being solved. The new generation, which is used for the next iteration, is then formed through bioinspired operators, such as selection, mutation, and crossover, based on the fitness. Generally, the algorithm terminates when either a maximum number of generations has been generated or a satisfactory fitness level has been reached for the population.The iteration of theEM algorithm is terminated if the convergence condition Θ^(l+1)-Θ^(l)≤ς is satisfied, where ς is a preset threshold.


## 3.3. RUL Fusion Estimation

From (3), we can see that the RUL is defined as the first passage time of stochastic degradation process Xt. As we know, the first passage time of Wiener process is inverse Gauss distribution, and thus the probability distribution function flk∣Yk of RUL can be obtained by(21)flk∣Yk=ω-x^kσ2+Pkη2πPk+σ2tk3exp⁡-ω-x^k-ηlk22Pk+σ2lk,where Elk∣Yk=ω-x^k/η and var⁡lk∣Yk=ω-x^kσ2+Pkη/η3.


## 4. Simulation Results

In this section, a simulation example is provided to illustrate the feasibility and effectiveness of the proposed asynchronous RUL fusion estimation algorithm.The stochastic degradation is formulated by (1) with ϕ=0, η=0.3, and σ2=0.01. We assume the degradation process is monitored by two asynchronous sensors with sampling periods T1=0.16 hours and T2=0.22 hours. The fusion period is Tf=0.2 hours. The observation coefficients of the two sensors are h1=1.03 and h2=0.76 with measurement noise variances r1=0.46 and r2=0.55, respectively. Obviously, the two sensors are asynchronous because of different sampling frequencies.The proposed asynchronous RUL fusion estimation algorithm is used to estimate the RUL of the hidden degradation process. Figure2 gives the actual degradation path and the estimated one. It can be seen from Figure 2 that the estimated degradation path has good ability to track the actual changes.Figure 2
The actual degradation path and the estimated.In the simulation, the failure thresholdw is set to 11.3, and the true failure time (i.e., the first passage time) of the degradation path given in Figure 2 is at 40 hours. To illustrate the effectiveness of the proposed fusion algorithm, the Probability Distribution Functions (PDFs) of RUL at different time instants are compared in Figure 3, where the PDFs drawn by solid lines are obtained by fusing observations from both of the asynchronous sensors, while the dashed and the dot dashed lines are results of single sensor 1 and sensor 2, respectively. In order to show this more clearly, the PDFs of RUL at a specific time (36 hours) are also given in Figure 4. It can be seen from Figures 3 and 4 that at each fusion time instant the fused PDF is more concentrated than that of single sensor. Meanwhile, the mean of the fused estimate is nearer to the true RUL (denoted with the asterisk). This means that the fused estimate has smaller estimation error as well as smaller error variation, in other words, better estimation performance, than the estimate obtained from single sensor observations. Meanwhile, it can be noticed in Figure 3 that the estimation error decreases with the increase of time.Figure 3
The PDFs of RUL at different time instants.Figure 4
The PDFs of RUL at 36 hours.In addition, it can be seen from Figure1 that “synchronous” can be regarded as a special case of “asynchronous.” Therefore, the proposed asynchronous fusion estimation algorithm is compared with the synchronous fusion estimation method in [18] under synchronous scenario in Figure 5, where T1=T2=Tf=0.2 hours. From Figure 5, we can see that the proposed algorithm outperforms the fusion algorithm in [18], since it has smaller mean error and more concentrated PDFs of RUL.Figure 5
Comparison of the proposed algorithm and the approach in [18] under synchronous scenario.


## 5. Conclusion

In this paper, an asynchronous RUL fusion estimation algorithm has been proposed for the latent degradation process with multiple asynchronous monitoring sensors. The asynchronous indirect observations are firstly synchronized to the fusion time and then fused using the Kalman filtering technology to get the estimate of the latent degradation state with the correlations between the involved noise considered and the unknown model parameters identified by theEM andGA. Finally, the fused estimate of the latent degradation state and the updated model parameters are used to get the probability distribution of the RUL. Simulation results demonstrate that fusion of multisource monitoring information from multiple asynchronous sensors can reduce the uncertainty in the systems and that improves the RUL estimation preformation.Due to the noise correlation introduced by the synchronization process,GA is adopted in this paper to update the estimate of unknown model parameters in theEM algorithm. Further work can focus on deriving an analytical solution by decollating the involved noises.



---

*Source: 4139563-2017-08-13.xml*