The Scientific World Journal

Volume 2014 (2014), Article ID 104052, 14 pages

http://dx.doi.org/10.1155/2014/104052

Research Article

Anomaly Monitoring Method for Key Components of Satellite

Jian Peng,¹ Linjun Fan,¹ Weidong Xiao,¹ and Jun Tang²

¹College of Information Systems and Management, National University of Defense and Technology, Changsha 410073, China

²Department of Telecommunications and Systems Engineering, University of Barcelona, 08202 Barcelona, Spain

Received 29 October 2013; Accepted 8 December 2013; Published 22 January 2014

Academic Editors: J. M. Corchado Rodríguez and M. Riera-Guasp

Copyright © 2014 Jian Peng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper presented a fault diagnosis method for key components of satellite, called Anomaly Monitoring Method (AMM), which is made up of state estimation based on Multivariate State Estimation Techniques (MSET) and anomaly detection based on Sequential Probability Ratio Test (SPRT). On the basis of analysis failure of lithium-ion batteries (LIBs), we divided the failure of LIBs into internal failure, external failure, and thermal runaway and selected electrolyte resistance () and the charge transfer resistance () as the key parameters of state estimation. Then, through the actual in-orbit telemetry data of the key parameters of LIBs, we obtained the actual residual value () and healthy residual value () of LIBs based on the state estimation of MSET, and then, through the residual values ( and ) of LIBs, we detected the anomaly states based on the anomaly detection of SPRT. Lastly, we conducted an example of AMM for LIBs, and, according to the results of AMM, we validated the feasibility and effectiveness of AMM by comparing it with the results of threshold detective method (TDM).

1. Introduction

As primary means of human space exploration, aerospace technology is playing an increasingly important role in military, politics, economy, technology, and many other fields. Meanwhile, as highly requirement of satellite, the system of satellite becomes more and more complex. In such complex and huge investment system, there is an urgent need to improve the reliability and security of satellite system. Otherwise, a failure could lead to catastrophic consequences and enormous economic losses [1]. According to the statistics, the number of spacecrafts including satellite and space station had been launched is 764 between 1990 and 2001, but there are 121 spacecrafts that have been failed accounting for 15.8% of the total number of spacecrafts [2].

Fault diagnosis is an important problem that has been researched within diverse research areas and application domains. Many fault diagnosis techniques have been specifically developed for certain application domains, while others are more generic [3]. According to the view of German fault diagnosis authority Professor P. M. Frank, fault diagnosis methods can be divided into three kinds: based on knowledge method; based on signal process method; based on analytical model method.

As rapidly developing artificial intelligence and computer technology, it provides a new theoretical foundation to fault diagnosis. As a result, the knowledge-based fault diagnosis method (KFDM) shows strong vitality and superiority. Moreover, because of its intelligence, the KFDM is becoming the focus of fault diagnosis and has made ​ ​a lot of achievements. In accordance with the intelligent method used, the KFDM can be divided into based-on-neural-network fault diagnosis method (NNFDM), based-on-expert-system fault diagnosis method (ESFDM), based-on-consistency-checking fault diagnosis method (CCFDM), and based-on-qualitative-model fault diagnosis method (QMFDM). Although the NNFDM [4–7] requires a lot of historical data to train network prior, it needs to address the problems of data collection. In addition, the scale of neural network (number of layers and hidden nodes) is always determined by experiment or experience, which does not have mature theoretical guidance. The serious drawback of ESFDM [8, 9] is that it only can detect predictable failure, and the bottleneck of its diagnosis ability is the source of expert knowledge. Moreover, there is a logical reasoning “combinatorial explosion” problem. The drawbacks of CCFDM [10–12] are when dealing with uncertain situations, its adaptive ability is very poor and bring poor false alarm rate and missed alarm rate; the preset logic decision is not necessarily the optimal approach; the production rules can only express predictable failure. The drawbacks of QMFDM [13–19] are: the accuracy of diagnosis entirely depends on the model; the potential failure diagnosis capability is weak; the models currently used are discrete model, so first of all the successive system model should be discretized, which would affect the accuracy of diagnosis. In a word, the KFDM is not suitable for key components of satellite.

According to the input and output signals of system which directly can be measured, the based on signal processing fault diagnosis method (SPFDM) [21–25] finds out the relationship between these signals and the source of failures. The SPFDM extracts the amplitude, phase, spectral, and any other characteristic values by using ​ ​the related functions, higher-order statistics, spectrum, autoregressive moving average, wavelet technology, and so on and detects the anomalies with these characteristic values. The advantages of SPFDM are that it is easy to be achieved and has good effectiveness, but the disadvantage of SPFDM is that the ability to detect the potential early anomaly is not enough. So the SPFDM is also not suitable for key components of satellite.

The based on analytical model fault diagnosis method (AMFDM) [26–29] is the earliest and most thorough, most proven methods, which can be divided into state estimation method, equivalent space method, and parameter estimation method. However, due to the fact that the particularity and complexity of space environment are complex and other uncertain factors are numerous, it is difficult to establish an accurate analytical model; as a result the equivalent space method and parameter estimation method are unsuitable for the fault diagnosis problems of satellite. Therefore, there is only the state estimation method that can solve the fault diagnosis problems of satellite. Multivariate State Estimation Techniques (MSET) [30] is a nonlinear, nonparametric modeling method that was originally developed by Argonne National Laboratory (ANL) for high-sensitivity proactive fault monitoring applications in commercial nuclear power applications. Some papers [31, 32] had used the MSET to detect the anomaly of internal combustion engine, electronic products, and so on. Nevertheless, the MSET is not applied in the fault diagnosis of satellite so far. In view of the fact that the data of key components of satellite is too few, it is difficult to obtain large amount of historical data and sufficient knowledge on abnormal components and rules. Consequently, the traditional based-on-analysis model methods are not suitable for fault diagnosis of satellite.

Satellite power system is the heart of entire satellite system, whose mission is to provide electrical power to the satellite, as well as to store, distribute, and control the electrical power [33]. Because of its superior weight ratio of the energy, the volume ratio of the energy, and cycle life performance, lithium-ion batteries (LIBs) have big advantages in comparison to the traditional energy storage batteries, such as nickel-cadmium battery and nickel-hydrogen battery.

To sum up, grounded on the fault diagnosis problems of LIBs in satellite, we propose a fault diagnosis method for key components of satellite, which is called Anomaly Monitoring Method (AMM), based on study of Multivariate State Estimation Techniques (MSET) and Sequential Probability Ratio Test (SPRT). The AMM is made up of state estimation based on MSET and anomaly detection based on SPRT, which is developed by improving MSET and creatively combining MSET and SPRT.

The main contributions of this paper are as follows.

(i)We propose a fault diagnosis method for key components of satellite, which is called Anomaly Monitoring Method (AMM), and it is made up of state estimation based on MSET and anomaly detection based on SPRT.(ii)On the basis of the study of MSET, we improve the method of MSET, considering the fact that it is difficult to obtain large amount of historical data and sufficient knowledge on abnormal components and rules.(iii)According to practical engineering applications, we detailed introduce the key technologies of state estimation based on MSET (i.e., data normalization process, selection of training data , selection and optimization of nonlinear operator, and structure of memory matrix ).(iv)We conduct quantitative measurement experiment to validate the feasibility and effectiveness of AMM by comparing it with the result of threshold detective method (TDM).

The rest of this paper is organized as follows. The analysis failure of LIBs is presented in Section 2. Section 3 discusses the state estimation based on MSET. Section 4 describes the anomaly detection based on SPRT. Section 5 provides the result of AMM and compares it with the result of TDM. Section 6 presents the conclusion and recommendations for future work.

2. Analysis Failure of LIBs

Compared with the traditional batteries, such as nickel-cadmium battery and nickel-hydrogen battery, LIBs are different from the general chemistry power supply, whose charging and discharging operations are done by insertion and extraction of their positive and negative. When the LIBs are charging, the positive electrode releases the lithium-ion in the electrolyte and this process is called extraction. And then the negative electrode intakes the lithium-ion from the electrolyte and this process is called insertion. Similarly, when the LIBs are discharging, the procedure occurs contrary to the above process. Therefore, due to the iterative insertion and extraction of lithium-ion, the process of charge and discharge of LIBs seems like a swaying rocking chair, that is why the LIBs are always called “rocking chair battery.”

On the basis of analysis failure of LIBs, we know that the failure of LIBs can be divided into internal failure, external failure, and thermal runaway. The internal failure is caused by three main aspects [34]: degradation of the electrode performance; loss of the electrolyte; aging of the diaphragm. These previous aspects have resulted in the increase of internal resistance which directly leads to the LIBs failure. But the external failure of LIBs is the result of many factors. Not only the effect of the structure, but also the method of usage. And the involving factors are the use procedure, the environment (temperature-sensitive), the charge and discharge mechanism, the lithium-ion electric stress, and so on. As the thermal reaction between the electrolyte and the electrode material, which is called “thermal runaway,” always leads to a devastating accident, the safety floor of LIBs is determined by their initial temperature [34]. According to the analysis of the internal failure, external failure, and thermal runaway of LIBs, the key performance parameters of LIBs are internal resistance, discharge current, operating voltage, and temperature. However, LIBs as an energy carrier can be equivalent to a constant voltage source and a resistor in the ideal electrical calculations. The study shows that the status and performance of LIBs are closely related to the impedance change [34].

Therefore, the study of impedance of LIBs has become a hot issue in manufacturing, testing, and surveillance of LIBs. Through the working principle of LIBs research, the LIBs’ impedance equivalent circuit model is shown in Figure 1, where represents the inductance, represents the electrolyte resistance, represents the charge transfer resistance, represents the capacitance of the electrode polarization which is the electric double layer capacitor, and represents the Warburg impedance. In general, and can be represented by . So the internal resistance of the LIBs can be represented by the electrolyte resistance () and the charge transfer resistance (). In practical engineering applications, and can be measured by electrochemical impedance spectroscopy method (EIS).

Figure 1: LIBs impedance equivalent circuit model.

Therefore, according to analysis failure of LIBs, we select and as the key parameters of state estimation based on MSET.

3. State Estimation Based on MSET

MSET is an advanced pattern recognition technology, which achieves the state estimation by measuring the similarity among various monitoring parameters in the normal operating range. The basic process of state estimation based on MSET is shown in Figure 2.

Figure 2: State estimation process based on MSET [20].

3.1. Steps of State Estimation

The basic steps of state estimation based on MSET can be roughly described as follows.

Step 1 (get observation matrix of monitoring parameters). On the basis of actual in-orbit telemetry data, the observation matrix of monitoring parameters is selected from it, which is denoted by . The observation matrix is ​ ​a matrix of , where represents the number of monitoring parameters and represents the number of the states. The actual in-orbit telemetry data matrix defined by MSET is shown in Figure 3.

Figure 3: Actual in-orbit telemetry data matrix in MSET [20].

Step 2 (select training data). The Training data is a matrix that consists of many healthy historic states, which is denoted by the letter . The training data is a matrix of , where also represents the number of monitoring parameters and represents the number of the healthy states. The selection requirements of training data are the states of training data must include the dynamic change states of monitoring parameters; the training data cannot contain any unhealthy states.

Step 3 (create memory matrix). The memory matrix is selected from training data in accordance with appropriate selecting rules, which is denoted by the letter . The memory matrix is a matrix of , where represents the number of monitoring parameters and represents the number of the states contained in the memory matrix .

The steps of create memory matrix are as follows:

(1)artificially determine according to the number of monitoring parameters;(2)select the extreme states of all monitoring parameters from training data . If the number of extreme states is equal to , then the memory matrix is completely created; otherwise, go to ;(3)if the number of extremely states is less than (the number of extremely states cannot be more than ), then continue to follow the next selection rules: firstly, calculate the Euclidean norm of the remaining training data and then get the vector ; secondly, arrange the vector by descending or ascending, and then get the vector ; lastly, according to equidistant sampling method, select the remaining states of memory matrix from the vector .

Step 4 (create remaining training data). After the memory matrix has been selected, those nonselected states in the training data form a new matrix called remaining training data which is defined as the letter . The remaining training data is a matrix of , where represents the number of monitoring parameters and is the number of the states contained in the remaining training data . So the remaining training data , the memory matrix , and the remaining training data have the following equation relationship:

Step 5 (calculate estimation matrix for the observation matrix ). The estimate of observation matrix , defined as , is the expected value calculated from the healthy data. This estimate has the same data format as the observation, which is an -element vector that is thought to be the weighted (linear) combination of states in memory matrix . The calculation equation is

where is a weight vector, which decides a similarity measure between each state in memory matrix and the estimation matrix . The formula for this vector is derived from the least square method by minimizing the error vector:

when the is constrained by minimization (set ), the least square error estimation solution of weight vector is:

So through (2), (3), and (4), the estimation matrix for the observation matrix is

Step 6 (calculate the estimation matrix of the remaining training data ). Similarly, the estimation matrix for the remaining training data defined as , which is the estimates of all the remaining data . The calculation equation is

where is also a weight vector, which decides a similarity measure between each state in memory matrix and the estimation matrix . The same as above, the estimation matrix for the remaining training data is

Step 7 (calculate residual value). MSET then calculates the residual value between the estimates () and remaining training data . Because all these remaining training data are healthy, the residual value presents the features of healthy states of the product and is called healthy residual, which is defined as . MSET also calculates the residual value between the estimates () and the observation matrix ; this residual value shows the actual states of the product and is called actual residual, which is defined as . The definition of residual value is the value of the estimated value minus the actual value. Therefore, the calculation equations of the actual residual and healthy residual are

In summary, through the introduction above, the basic steps of the state estimates based on MSET are clear. But in practical engineering applications, the state estimation based on MSET is more complex than the above-mentioned steps, which involves, in the data normalization process, the selection of training data , the selection and optimization of nonlinear operator, and the structure of memory matrix . So next there is a detailed introduction of the state estimation based on MSET in practical engineering applications.

3.2. Data Normalization Process

In practical engineering application, the selected parameters of product might not be the same order of magnitude. And due to the fact that the weight vector represents a similarity measure between the observation matrix and the memory matrix , the value of weight vector is completely determined by high order of magnitude parameter data, such as selecting , , and capacity as the key parameters, shown in Table 1. Obviously, the internal resistance data and the capacity data are clearly not the same order of magnitude. So at this time the value of weight vector is completely determined by the capacity data, and the internal resistance data have lost the meaning of existence. Moreover, according to Section 3.4, if the observation matrix is not normalized, the nonlinear operator is not available to be optimized. So, in order to solve these problems, the data of each monitoring parameter firstly should be normalized, just as shown in Table 2.

Table 1: Actual observation matrix .

Time

…

0.078
0.088
0.078
0.085
0.078
0.084
0.078
0.084
0.078
0.083
…

0.061
0.065
0.060
0.062
0.061
0.062
0.060
0.062
0.059
0.062
…

Capacity
2.035
2.025
2.013
2.013
2.000
2.014
2.013
1.968
1.968
1.957
…

Table 2: Observation matrix after data normalization process.

Time

…

−0.809
1.738
−0.939
0.897
−0.844
0.664
−0.927
0.655
−0.956
0.522
…

−0.245
2.229
−0.586
0.540
−0.456
0.328
−0.905
0.300
−1.432
0.228
…

Capacity
1.291
0.909
0.466
0.464
−0.013
0.487
0.045
−1.202
−1.225
−1.635
…

3.3. Selection and Update of Training Data

According to the introduction of training data in Section 3.1, the training data is a matrix that consists of many healthy historic states. However, in practical engineering applications, LIBs of satellite do not have historical healthy data. Therefore, the traditional MSET be cannot suitable for the fault diagnosis of satellite in practical engineering applications. In order to resolve this problem, we improve the MSET on the basis of the traditional MSET. So we select the training data from the actual in-orbit telemetry data and consider the initial data of the actual in-orbit telemetry data as the training data , shown in Figure 4. But in order to ensure that the states contained in training data are healthy states, the selected portion of states should be small. And then due to the fact that the trading data selected originally does not contain all dynamic change states of monitoring parameters, the trading data selected originally should be updated. There are two methods for training data update.

(1)Nonquantitative Update Method. With the continuation of anomaly detection, get the healthy states detected by AMM and add them to the training data selected originally. As shown in Figure 5, the (101–150) states are the healthy states detected by AMM in the first time and the (151–200) states are the healthy states detected by AMM in the second time, and get all the healthy states (101–200) detected by AMM and add them to the training data .(2)Quantitative Update Method. With the continuation of anomaly detection, only get the healthy states in the floating window every time detected by AMM and add them to the training data selected originally. As shown in Figure 6, the (151–200) states are the healthy states detected by AMM in the second time, and only get the healthy states (151–200) detected by AMM in the second time and add them to the training data selected originally.

Figure 4: Selection of training data from the actual in-orbit telemetry data.

Figure 5: Update of training data with nonquantitative update method.

Figure 6: Update of training data with quantitative update method.

Comparing the two update methods, the more states contained in the nonquantitative update method, the better state estimation effect of the nonquantitative update method might be. However, in practical engineering applications, the actual in-orbit telemetry data of LIBs is extremely large, so the calculated amount of the nonquantitative update method would be vast. Actually by using the quantitative update method, not only the training data reaches the purpose that contains dynamic change states of monitoring parameters but also the amount of calculation will be greatly reduced. Therefore, considering the calculated amount in practical engineering applications, the quantitative update method is better than the nonquantitative update method.

3.4. Selection and Optimization of Nonlinear Operator

Through (4), we know that if the weight vector is existent, the must be reversible. However, a necessary but not sufficient condition of the reversible is that the number of columns of the memory matrix should be less than the number of lines of the memory matrix ; in other words, the number of states contained in memory matrix should be less than the number of monitoring parameters. However, in practical engineering applications, this condition is very difficult to be fulfilled. In order to provide adequate statistical information of monitoring parameters, the memory matrix must include a large number of states. So the number of states contained in memory matrix cannot be less than the number of monitoring parameters. Actually the cannot be reversible, that is because of the correlation of different states contained in the memory matrix . Therefore, in order to solve this problem in practical engineering applications, MSET method introduction of the nonlinear operator is defined as , which is a typical operator that calculates normalized similarity between different data vectors.

The most common nonlinear operators are [35]

(1)the Euclidean norm (),(2)city block distance (),(3)linear correlation coefficient (),(4)the Gaussian operator, and so on.

In this paper, we use the most commonly used nonlinear operator—the Gaussian operator. The Gaussian operator equations can be described as follows:

where is the filter coefficient or also known as bandwidth.

Then, through the selection of nonlinear operator, (4), (5), and (7) can be described as follows:

However, if the states contained in the memory matrix have strong correlation, the state estimation will be very badly. This is because the cardinal number of the will become tremendous, and then the estimated value will also become tremendous. So, under these circumstances, the state estimation based on MSET is equivalent of a magnifying glass. If the observation matrix includes noise, then it is difficult to obtain a stable state estimation, and this restricted memory matrix is like a noise amplifier, just as shown in Figure 7.

Figure 7: Memory matrix with strong correlation.

If the states contained in memory matrix do not have strong correlation, the state estimation based on MSET will be good, and the estimated value will not amplify the noise, just as shown in Figure 8.

Figure 8: Memory matrix without strong correlation.

In summary, if the observation matrix includes noise, the states contained in memory matrix could not have strong correlation when using the nonlinear operator; otherwise, the state estimation based on MSET will be like a noise amplifier. However, in practical engineering applications, the observation matrix inescapable includes some noise and the states contained in memory matrix cannot be completely without strong correlation. So, in order to resolve this problem, we introduced the ridge regularization when inverse the . And then after that, (10) can be described as follows:

where is the ridge regularization parameter () and is the identity matrix.

The optimization problem of is a one-dimensional optimization problem, which can be solved by using nonlinear methods, such as conjugate gradient descent and so on. When the data of each monitoring parameter follow normal distributions with mean 0 and variance 1, the optimization bandwidth is and the optimization ridge regularization is [36].

3.5. Structure of Memory Matrix

On the basis of the introduction of creating memory matrix in Section 3.1, the selection of memory matrix can be roughly divided into two steps:

(1)determine the number of states contained in the memory matrix , ;(2)select states from the training data with appropriate rules.

The most important steps of creating memory matrix is to determine the number of states contained in the memory matrix which is closely related to the estimation performance. Generally, estimates produced by using a lesser number of states contained in the memory matrix tend to be less accurate and less reliable than those created by using more states. Under some circumstances, using too many states can capture signal noise and lead to undesirable effects (the “overtraining” syndrome), just as shown in Figure 7. In addition, consideration of the spend time of state estimation, the spend time of state estimation based on MSET is closely related to the number of states contained in the memory matrix , that is, the more states contained in the memory matrix , the longer time spent by state estimation [31].

On the other hand, the selection of memory matrix is directly related to the selection of remaining training data , according to (1). So if the number of states contained in the memory matrix is too big, it will inevitably lead to the fact that the number of states contained in the remaining training data is too small, which is obviously not conducive to the SPRT detection.

To sum up, there is a basic principle for the determination of the number of states contained in the memory matrix : the states contained in the memory matrix should be as few as possible when the states cover all the dynamic range of monitoring parameters.

4. Anomaly Detection Based on SPRT

Through the state estimation based on MSET, the actual residual value and healthy residual value have been obtained, and then the anomaly detection can be going based on these residual values. SPRT is a common technique for binary hypothesis testing, which has been widely used to conduct binary hypothesis testing in many applications. Wald’s SPRT [37] was developed based upon two assumptions: samples are independent and identically distributed; samples follow a priori known distribution function. As the SPRT method can detect not only the anomaly states outside the threshold but also the anomaly state in the threshold, so in this paper we detect the anomaly states based on SPRT.

However, the actual residual value and healthy residual value obtained by MSET are both matrices. So in view of this situation this paper puts forward a preprocessing method for the residual values before anomaly detection based on SPRT: reduce the dimension of the actual residual value and healthy residual value by introducing a weight vector , where is the weight ratio of the monitoring parameter on product performance, just as follows:

where and represent the residual value of the actual residual value and the healthy residual value separately when the monitoring parameter is on the state and and represent the residual value of the actual residual value () and the healthy residual value () after reducing the dimension separately on the state .

The actual residual value can be tested by mean and variance based on SPRT when the residuals ( and ) obey normal distribution. Then set the original hypothesis to be , and the alternative hypothesis is any one of :

where is the variance factor, and are the mean and variance of healthy residual value, and are the mean of anomaly states, and are the alternative hypothesis for the mean test, and and are the alternative hypothesis for the variance test.

In the classic Wald SPRT method, a sequential probability ratio is constructed as

where is the likelihood function for observing the sample sequence if the hypothesis is true and is the likelihood function for observing the sample sequence if the hypothesis is true. Deciding which hypothesis to accept is determined by comparing the ratio with some upper and lower bound values, which being related to and , by the following rules:

(1)if , accept ;(2)if , accept ;(3)if cannot make a decision, increase the sample size to and then repeat the comparison,

where is called prescribed first kind of error and is called prescribed second kind of error, and in practical engineering applications, and are also called the prescribed false alarm rate and prescribed miss alarm rate, respectively.

Next, the mean test of actual residual based on SPRT can be detailed described and as follows.

According to (13), we can get that the original hypothesis is , and the alternative hypothesis is . And then the likelihood function of actual residual can be described as follows:

So according to (14), the likelihood ratio is

And then the logarithmic function of the likelihood ratio is

Lastly, let substitute in the decision expression of SPRT (), namely,

The meaning of the in the upper equations does not represent the total length of the test sample but represents the length of the test sample when the decision is done.

When has been calculated, the SPRT Index of residual sample can be calculated, and the calculation steps are as follows.

Step 8. Set the logarithmic of the likelihood ratio of each state to be equal to .

Step 9. If , then set ; if , then set ; if , then set , where is the lower limit and is the upper limit.

Step 10. If , then set , or set .

Step 11. Then the , where corresponds to the abnormal state and indicated by the diamond in the figure, where corresponds to the healthy state and indicated by the circle in the figure, where corresponds to the state that is unable to be determined and indicated by the dot in the figure. The steps of calculation also can be orderly expressed in Figure 9.

Figure 9: Calculation steps of SPRT Index.

5. AMM Applications for LIBs

The AMM consists of state estimation based on MSET and anomaly detection based on SPRT. In Section 3, we have detailed introduction of the steps of LIBs’ state estimation based on MSET. In Section 4, according to the actual residual value and healthy residual value of LIBs that are got by state estimation, we have detailed illustration of the mean test of actual residual based on SPRT. Next, we conduct an example of AMM for LIBs, and, according to the result of AMM, we validate the feasibility and effectiveness of AMM by comparing it with the result of TDM.

5.1. State Estimation for LIBs

According to the steps of state estimation based on MSET, first of all is to select parameters and get observation matrix of monitoring parameters . In Section 2, according to analysis failure of LIBs, we select and of LIBs as the key parameters of state estimation based on MSET. And then, on the basis of actual in-orbit telemetry data of LIBs, we get observation matrix of and . Meantime, select training data from actual in-orbit telemetry data of LIBs according to the introduction of Section 3.3, and the training data and observation matrix of and are shown in Figure 10, where dot represents the training data of triangle represents the observation matrix of , five-pointed star represents the training data of , and plus sign represents the observation matrix of .

Figure 10: Training data and observation matrix of and .

Secondly, according to the introduction of data normalization process in Section 3.2, the training data and observation matrix of and should be normalized with mean 0 and variance 1 of normal distribution before state estimation, so the data after normalization process are shown in Figure 11, where label represents the value of resistance value after normalization process.

Figure 11: Training data and observation matrix of and after normalization process.

Thirdly, along with the steps of state estimation based on MSET, we create memory matrix and remaining training data of and from training data in Figure 11, just as shown in Figure 12, where the left figure is the memory matrix of and and the right figure is the remaining training data of and .

Figure 12: Construction of memory matrix and the remaining training data of and .

And then through the construction of memory matrix and the remaining training data of and in the upper figure and the observation matrix in Figure 11, the weight vector , the estimation matrix of observation matrix , , and the estimation matrix of remaining training data , can be calculated. Just as shown in Figure 13, the left figure is the comparison of observed value and estimated value of and the right figure is the comparison of observed value and estimated value of .

Figure 13: Comparison of the observed value and estimated value of and .

Lastly, the actual residual value and healthy residual value of and can be calculated by MSET, just as shown in Figure 14, where left figure is the actual residual value of and and right figure is the healthy residual value of and .

Figure 14: Actual residual value and healthy residual value of and .

5.2. Anomaly Detection for LIBs

According to actual residual value and healthy residual value of LIBs obtained by MSET, the anomaly detection for LIBs can be carried out based on SPRT. But according to the introduction of Section 4, firstly, the residual matrix should be preprocessed before anomaly detection based on SPRT. On the basis of (12), the residual matrix is transformed into one-dimensional vector by introducing the weight vector , shown in Figure 15, where the left figure is the actual residual value after dimensionality reduction and the right figure is the healthy residual value after dimensionality reduction.

Figure 15: Actual residual value and healthy residual value of and after dimensionality reduction.

So after that, the mean test of residual value will begin. Secondly, we formulate the following assumptions:

Thirdly, we select and and then get the upper limit of detection and the lower limit of detection .

Lastly, in accordance with (14)–(18) and calculation steps of SPRT Index, the final SPRT Index of AMM is obtained, and the result of AMM is shown in Figure 16, where diamond represents the abnormal states (such as state 6, state 33, and state 74), dot represents the unable determined states and circle represents the healthy states.

Figure 16: Result of AMM.

However, based on the data in Figure 11, the result of TDM is shown in Figure 17, where the abnormal states are state 6, state 33, state 82, and so on. Comparing these two results, we demonstrate that the AMM can detect not only the abnormal states outside the threshold (such as state 6 and state 33) but also the abnormal states in the threshold (such as state 74); furthermore, the AMM detects the abnormal state ahead of eight time states alarm than the TDM.

Figure 17: Result of TDM.

6. Conclusions

This research proposes a fault diagnosis method for key components of satellite, called AMM, based on the study of MSET and SPRT, which is developed by improving MSET and creatively combining MSET and SPRT. Due to the importance of satellite power system in the entire satellite system and the superiority of LIBs, we select the LIBs as the subject investigated. Then, after analysis failure of LIBs, we select and as the key parameters of AMM. The AMM consists of state estimation based on MSET and anomaly detection based on SPRT. On the basis of fault diagnosis problems of LIBs in practical engineering applications, we improve the MSET and expound the key technologies of MSET at great length, that is, data normalization process, selection of training data , selection and optimization of nonlinear operator, and structure of memory matrix . And then according to the actual residual value and healthy residual value of LIBs that are got by state estimation, we have detailed illustration of the mean test of actual residual based on SPRT. Lastly, we conduct an example of AMM for LIBs, and, according to the result of AMM, we validate the feasibility and effectiveness of AMM by comparing it with the result of TDM.

The investigation in this paper is conducted from an engineering perspective. Extensive applications demonstrate that the model and algorithm are feasible, effective, and advanced. Studies on AMM can offer theoretical and technical guidance for fault diagnosis.

Recommendations for future work are threefold: An in-depth analysis failure of LIBs should be conducted and more influencing parameters should be considered; More studies on the key technologies of MSET and fault diagnosis methods should be carried out; to serve as guidance for prognostics and health management using our methods.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

X. P. Hu, Q. D. Han, and J. H. Li, Data Mining in Anomaly Detection, National University of Defence and Technology, 2009, (Chinese).

L. X. Lin, The Failure Analysis and Research of Spacecrafts' Control, Guide and Navigation System during 1990–2001, International Space, 2005, (Chinese).

V. Chandola, A. Banerjee, and V. Kumar, “Anomaly detection: a survey,” ACM Computing Surveys, vol. 41, no. 3, article 15, 2009.

Q. Chen and Z. Li, “The fault diagnosis based on immune neural network,” in Proceedings of the 3rd International Symposium on Information Science and Engineering (ISISE '10), pp. 454–457, Shanghai, China, December 2010.

D. Chen and P. Zhao, “Study of the fault diagnosis method based on RBF neural network,” in Proceedings of the 2nd International Conference on Artificial Intelligence, Management Science and Electronic Commerce (AIMSEC '11), pp. 4350–4353, Zhengzhou, China, August 2011.

W. J. Xu and J. Yao, “Integrated quality diagnosis algorithm method based on neural network and sensitivity analysis to input parameters,” Journal of Networks, vol. 8, no. 6, pp. 1307–1314, 2013.

J. Reyes, M. Vellasco, and R. Tanscheit, “Monitoring and diagnosis of multiple sensors by auto-associative neural networks,” Automatic Control, vol. 23, no. 2, pp. 121–137, 2012.

P. W. Tse, E. Y. Li, J. C. Chan, and J. T. Leung, “Automatic generator health assessment system that embedded with advanced fault diagnosis and expert system,” in Proceedings of the Prognostics and System Health Management Conference (PHM '10), Macau, China, January 2010.

Z. Zhou, M. Zhuang, X. Lu, L. Hu, and G. Xia, “Design of a real-time fault diagnosis expert system for the EAST cryoplant,” Fusion Engineering and Design, vol. 87, no. 12, pp. 2002–2006, 2012.

D. DeCoste, “Learning envelopes for fault detection and state summarization,” in Proceedings of the 2000 IEEE Aerospace Conference, pp. 337–344, Big Sky, Mont, USA, March 2000.

M. A. Swartwout, C. A. Kitts, and R. K. Batra, “Persistence-based production rules for on-board satellite automation,” 1999 IEEE Aerospace Applications Conference, pp. 273–281, 1999.

G. J. Mansingh, L. Rao, B. Osei, M. Kweku, and A. Mills, “Profiling internet banking users: a knowledge discovery in data mining process model based approach,” in Information Systems Frontiers, pp. 1–23, Springer, 2013.

A. Kumar, A. Srivastava, N. Goel, and A. Nayak, “Model-free approach and methodology for data anomaly detection for real time diagnostic solution,” in Proceedings of the 23rd Canadian Conference on Electrical and Computer Engineering, Calgary, Canada, May 2010.

S. Hiratsuka, H. H. Lu, and A. Fusaoka, “Model-based monitoring and diagnosis chip for embedded systems,” in Proceedings of the 12th International Conference Artificial Intelligence: Methodology, Systems, and Applications, Lecture Notes in Computer Science, pp. 255–264, Varna, Bulgaria, 2006.

C. M. Meyer, C. Fulton, B. Maul, A. Chicatelli, and H. Cannon, “Propulsion IVHM technology experiment overview,” in Proceedings of the 2003 IEEE Aerospace Conference, pp. 859–868, Big Sky, Mont, USA, 2003.

S. C. Hayden, A. J. Sweet, and S. E. Christa, “Livingstone model-based diagnosis of earth observing one,” in Proceedings of the AIAA 1st Intelligent Systems Technical Conference—Collection of Technical Papers, pp. 98–108, Chicago, Ill, USA, September 2004.

S. C. Hayden, A. J. Sweet, S. E. Christa, D. Tran, and S. Shulman, “Advanced diagnostic system on earth observing one,” in Proceedings of the AIAA Space 2004 Conference and Exposition—A Collection of Technical Papers, pp. 2209–2222, San Diego, Calif, USA, September 2004.

M. Schwabacher, J. Samuels, and L. Brownston, “The NASA integrated vehicle health management technology experiment for X-37,” in Proceedings of the Components and Systems Diagnostics, Progrnostics, and Health Management II, Proceedings of SPIE, pp. 49–60, The International Society for Optical Engineering, Orlando, Fla, USA, April 2002.

A. Kumar, A. Srivastava, N. Goel, and A. Nayak, “Model-free approach and methodology for data anomaly detection for real time diagnostic solution,” in Proceedings of the 23rd Canadian Conference on Electrical and Computer Engineering (CCECE '10), Calgary, Canada, May 2010.

S. Cheng and M. Pecht, “Multivariate State Estimation Technique for remaining useful life prediction of electronic products,” in Proceedings of the Papers from the AAAI Fall Symposium—Artificial Intelligence for Prognostics, pp. 26–32, Arlington, Va, USA, November 2007.

M. Chen, S. Zhang, Y. Xing, and Y. Zhang, “Combined attitude control method of small satellite in the case of reaction wheel failed,” Journal of Harbin Institute of Technology, vol. 39, no. 5, pp. 811–816, 2007.

Q. Wu and M. Saif, “Robust fault detection and diagnosis for a multiple satellite formation flying system using second order sliding mode and wavelet networks,” in Proceedings of the American Control Conference (ACC '07), pp. 426–431, New York, NY, USA, July 2007.

B. Banerjee, L. Kraemer, and W. Solano, “Particle filtering for diagnosis and prognosis of anomalies in rocket engine tests,” in Proceedings of the 2011 AIAA Infotech at Aerospace Conference and Exhibit, St. Louis, Mo, USA, 2011.

A. Borisov, G. Runger, and E. Tuv, “Target-aware anomaly detection and diagnosis,” in Proceedings of the 8th International Conference on Informatics in Control, Automation and Robotics (ICINCO '11), pp. 14–23, Noordwijkerhout, Netherlands, July 2011.

H. Ding, H. Sun, and K. Hou, “Abnormal ECG signal detection based on compressed sampling in Wearable ECG sensor,” in Proceedings of the International Conference on Wireless Communications and Signal Processing (WCSP '11), Nanjing, China, November 2011.

H. A. Talebi, R. V. Patel, and K. Khorasani, “Fault detection and isolation for uncertain nonlinear systems with application to a satellite reaction wheel actuator,” in Proceedings of the IEEE International Conference on Systems, Man, and Cybernetics (SMC '07), pp. 3140–3145, Montreal, Canada, October 2007.

K. Matsuoka, “Pitfalls in radar diagnosis of ice-sheet bed conditions: lessons from englacial attenuation models,” Geophysical Research Letters, vol. 38, no. 5, Article ID L05505, 2011.

L. Li, L. Ma, and K. Khorasani, “A dynamic recurrent neural network fault diagnosis and isolation architecture for satellite's Actuator/thruster failures,” in Proceedings of the 2nd International Symposium on Neural Networks (ISNN '05), pp. 574–583, Chongqing, China, June 2005.

N. Cheifetz, A. Samé, P. Aknin, and E. De Verdalle, “A pattern recognition approach for anomaly detection on buses brake system,” in Proceedings of the 14th IEEE International Intelligent Transportation Systems Conference (ITSC '11), pp. 266–271, Washington, DC, USA, October 2011.

K. C. Gross, R. M. Singer, S. W. Wegerich, and J. P. Herzog, “Application of a model-based fault detection system to nuclear plant signals,” International Conference on Intelligent Systems Applications to Power Systems, Seoul, Republic of Korea, 1997.

K. Vaidyanathan and K. C. Rross, “MSET performance optimization for detection of software aging,” in Proceedings of the International Symposium on Software Reliability Engineering (ISSRE '03), San Diego, Calif, USA, 2003.

E. Liu and D. Zhang, “Diagnosis of component failures in the space shuttle maine engines using Bayesian belief networks: a feasibility study,” in Proceedings of the 14th International Conference on Tools with Artificial Intelligence, pp. 181–188, Washington, DC, USA, November 2002.

H. L. Li and Y. Jia, “Failure mechanism analysis of lithium-ion battery,” in Proceedings of the 12th National Reliability Physics Symposium, Dujiangyan, China, 2007, (Chinese).

D. Kang and H. Shin, “Investigation on cell impedance for high-power lithium-ion batteries,” Journal of Solid State Electrochemistry, vol. 11, no. 10, pp. 1405–1410, 2007.

L. B. Christopher, E. U. Robert, and H. J. Wesley, “System modeling and instrument calibration verification with a nonlinear state estimation technique,” in Proceedings of the Maintenance and Reliability Conference, Knoxville, Tenn, USA, 1998.

H. J. Wesley and U. Alexander, Mset Performance Optimization through Regularization, Nuclear Engineering and Technology, 2005.

A. Wald, Sequential Analysis, Wiley, New York, NY, USA, 1947.