Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Hybrid Framework Model Based on Wavelet Neural Network with Improved Fruit Fly Optimization Algorithm for Traffic Flow Prediction
Previous Article in Journal
From Galactic Bars to the Hubble Tension: Weighing Up the Astrophysical Evidence for Milgromian Gravity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Dynamic GLR-Based Fault Detection Method for Non-Gaussain Dynamic Processes

1
College of System Engineering, National University of Defense Technology, Changsha 410073, China
2
School of Electrical and Information Engineering, Tianjin University, Nankai District, Tianjin 300072, China
3
School of Automation, Central South University, Changsha 410083, China
*
Authors to whom correspondence should be addressed.
Symmetry 2022, 14(7), 1332; https://doi.org/10.3390/sym14071332
Submission received: 21 May 2022 / Revised: 20 June 2022 / Accepted: 22 June 2022 / Published: 28 June 2022
(This article belongs to the Section Engineering and Materials)

Abstract

:
Non-Gaussian dynamic processes are ubiquitous due to the presence of non-Gaussian distributed variables. Therefore, fault detection of non-Gaussian dynamic processes plays a vital role to maintain the safe operation of systems and symmetry of data distribution. In this paper, a dynamic generalized likelihood ratio (DGLR)-based fault detection method is proposed for non-Gaussian dynamic processes. Different from the conventional principal component analysis (PCA)-based, dynamic PCA-based, and PCA-based GLR fault detection methods, the novelty of the proposed method is that the GLR is extended to non-Gaussian dynamic processes, and the randomized algorithm is integrated for threshold setting to attenuate the influence of non-Gaussian. The application scope of these methods is also discussed. The proposed method is compared with four existing fault detection methods on a numerical simulation and the continuous stirred-tank reactor (CSTR) process. The achieved results show that the proposed method is able to significantly improve the detection performance in terms of fault detection rate and prompt response to faults.

1. Introduction

Fault detection is becoming increasingly important to maintain high quality products, operation safety of processes, and symmetry of data distribution. In recent years, considerable attention has been paid to research on solving fault detection problems. Model-based and data-driven methods are two common types [1,2,3,4,5,6]. In the model-based approaches, accurate physical or mathematical models are needed. On the other hand, in data-driven methods, only the availability of historical process data is required. Data-driven fault detection techniques have been widely used due to the simple application form and fewer requirements on development. The commonly used data-driven methods include generalized likelihood ratio (GLR)-based [7,8], multivariate analysis (MVA)-based, such as principal component analysis (PCA) [9,10,11,12], partial least square (PLS) [13,14,15], canonical correlation analysis (CCA) [16,17,18].
PCA is one of the most MVA techniques used for FD, which considers a single dataset. The successful application of such a method can be found in a wide range of applications, for example, dynamical, time varying, non-Gaussian and nonlinear processes [19]. Several variants of the standard PCA have been developed. These variants include dynamic PCA which is used to find dynamical linear relationships between the process variables [20], moving window PCA which handles time varying features [21], kernel density based-PCA which is used in non-Gaussian fault detection [22]. For GLR-based fault detection methods, some variants have also been developed in the fields of time varying [23], non-Gaussian [24], etc. However, the dynamic characteristic in process is rarely considered in the existing methods, and the better detection performance could be achieved by considering the auto-correlation characteristic in dynamic processes.
Furthermore, the successful application of the GLR-based method requires that the process data follows a Gaussian distribution. In practice, processes with non-Gaussian features put forward more challenges for fault detection. To this end, some variants of the existing fault detection methods have been developed. Commonly, there are two types of methods to deal with the fault detection problem for non-Gaussian processes. The first ones are either to use some methods, which are free of distribution, such as independent component analysis (ICA) [19] and support vector machine-based methods [25,26], or to extract high-order statistics and then use the standard method based on the obtained statistics [27]. For example, in [26], the ICA method is firstly used to get the independent components and then the support vector data description method is applied to generate a suitable threshold. The authors in [27] used various statistics to quantify process characteristics, such as non-Gaussian, furthermore monitoring these statistics instead of process variables themselves to perform fault detection. The second type of methods are first to estimate a probability distribution of the monitored variables or statistics, and then determine an appropriate threshold based on the resulting distribution. These methods can be referred to as distribution estimation-based methods [28]. A great number of methods have been developed for the estimation of distribution, including Gaussian Mixture Models (GMM) [29,30], kernel-based approaches [31,32] and sequential quantile estimation [33]. Motivated by the success of the second ones, in this paper, we use the same strategy. Although these existing methods are successful in this application domain, their performance in fault detection is commonly constrained by the determination of kernel structure and method-specified parameters, for example, the bandwidth parameter for a Gaussian kernel [31]. Therefore, due to the ability to iterative updating threshold, a randomized algorithm-based threshold learning method is used to enhance the dynamic GLR (DGLR)-based method to deal with the fault detection task in non-Gaussian dynamic processes.
Motivated by the above analysis, a DGLR-based fault detection method combined with the threshold learning method is proposed for non-Gaussian dynamic processes. The contribution of this work is four-fold: (1) to develop a DGLR-based detection statistic for non-Gaussian dynamic processes; (2) to improve DGLR-based fault detection performance by iteratively learning the suitable threshold by a randomized algorithm; (3) to compare the DGLR-based fault detection methods with the GLR-based, PCA-based, DPCA-based, and PCA-based GLR ones [34]. Based on our best knowledge, there are few works to compare these methods with the purpose to clarify the application scope of these method and to guide the practitioners to select a suitable fault detection method; (4) to assess the DGLR-based fault detection performance by comparing it with the GLR-based, PCA-based, DPCA-based, and PCA-based GLR methods using a numerical simulation and the continuous stirred-tank reactor (CSTR) process.
Notation 1.
The notation used in this paper is standard. R n denotes the n-dimensional Euclidean space consisting of n × 1 vectors with real components, R n × m is the set of all n × m real matrices, and diag ( , , ) is a square diagonal matrix. A ( : , i ) represents the i-th column of A . I n is an n × n identity matrix. x N μ x , Σ x denotes that x is a normally distributed random vector with mean μ x and covariance Σ x . E ( · ) denotes the expectation operator. χ 2 m stands for the chi-square distribution with m degrees of freedom. Let pr ( χ 2 > χ α 2 m ) = α be the probability that χ 2 > χ α 2 m equals α (significance level).

2. Background and Problem Formulation

2.1. The Basics of GLR-Based Fault Detection Technique

Consider the following fault detection problem using a GLR-based technique. Given a general model
y = y * + f R m
where y * N ( 0 , Σ ) represents the statistical features of the process and m is the dimension of the variable. Since f = 0 denotes the fault-free case, our task consists of detecting a fault f 0 with N number of available measurements of y 1 , , y N .
The fault detection task can be solved by testing the following hypotheses based on the available data y [7]
H 0 , null hypothesis : f = 0 , fault - free , H 1 , alternative hypothesis : f 0 , faulty .
The probability density functions (pdf) of y * and y are respectively given as
P 0 ( y ) = 1 ( 2 π ) m det Σ exp 0.5 y T Σ 1 y ,
P 1 ( y ) = 1 ( 2 π ) m det Σ exp 0.5 ( y T E ( y ) ) Σ 1 ( y E ( y ) ) ,
The log likelihood ratio is defined as
s ( y ) = 2 ln P 1 ( y ) P 0 ( y ) = y T Σ 1 y ( y E ( y ) ) T Σ 1 ( y E ( y ) )
To increase the confidence of the decision-making procedure, generally more samples are required. Using N samples of data y , Equation (4) is extended as
S 1 N = k = 1 N 2 ln P 1 ( y k ) P 0 ( y k ) = k = 1 N y k T Σ 1 y k = 1 N ( y k E ( y ) ) T Σ 1 ( y k E ( y ) ) = 2 N y ¯ T Σ 1 E ( y ) N E ( y ) T Σ 1 E ( y ) , y ¯ = 1 N k = 1 N y k = N ( y ¯ T Σ 1 y ¯ ( y ¯ E ( y ) ) T Σ 1 ( y ¯ E ( y ) ) )
Evidently, the maximum of S 1 N is
N y ¯ T Σ 1 y ¯
when y ¯ = E ( y ) is achieved. Since E ( y ) is generally unknown, it can be replaced by its maximum likelihood estimate
E ( y ) = 1 N k = 1 N y k
In practice Σ is also unknown, which needs to be identified from the data. It is straightforward that
lim N Σ ^ = lim N 1 N 1 k = 1 N ( y k E ( y ) ) ( y k E ( y ) ) T = Σ
which gives an asymptotically unbiased estimate of the covariance matrix. Thus if N is sufficiently large, the unknown parameter Σ could be approximated by its estimate
Σ ^ = 1 N 1 k = 1 N ( y k E ( y ) ) ( y k E ( y ) ) T

2.2. Randomized Algorithm-Based Threshold Setting

The randomized method has been widely used to the analysis and design of a robust control system [35,36]. Recently, this method has been used for threshold setting for the non-Gaussian process because it is independent of the probability distribution [28]. This successful application can be attributed to the iterative update of the threshold by means of the estimation of false alarm rate (FAR). The basic idea behind is that the required threshold should guarantee a desired false alarm rate, which is predefined. Let J t h , p F A R and p ^ F A R be the threshold, allow FAR and the estimated FAR, respectively A lowest threshold which is satisfied with the given false alarm rate is obtained by the following Algorithm 1 from [28].
Algorithm 1: Threshold-setting
Given allowed FAR ( 0 , 1 ) and δ ( 0 , 1 ) , let ϵ > 0 be some constant satisfying
FAR ϵ > 0
 and Δ > 0 be the iteration tolerance.
S1: Set J t h = J 0 ( > 0 ) ;
S2: Choose integer N according to the one-sided Chernoff inequality with ϵ ( 0 , 1 ) , δ ( 0 , 1 )
N 1 2 ϵ 2 log 1 δ
S3: Estimate FAR using the method given in [37];
S4: If p ^ F A R p F A R ϵ then return J t h and exit;
S5: Else J t h = J t h + Δ go to Step 3.
As given in Theorem 1 in [28], for a sufficiently small Δ , the estimated threshold satisfies J t h , m i n J t h J t h , m i n + Δ with J t h , m i n the lowest threshold. Since Δ is sufficiently small, the estimated threshold approaches the lower threshold, i.e., J t h J t h , m i n .

3. The Proposed Method

In this section, a DGLR-based fault detection method is proposed for non-Gaussian dynamic processes. Firstly, a DGLR-based test statistic is built for detection purposes. Then, it is well known in probability theory that the complete information of a Gaussian distribution can be described by the mean value and the covariance. Therefore, if the measured variable y follows a Gaussian distribution, the threshold setting can be achieved using the standard distribution table. However, in a non-Gaussian case, the distribution of the test statistics inevitably deviates from the standard distribution, e.g., χ 2 distribution. In this case, the threshold, which is set based on Gaussian assumption, will decrease the detection performance, e.g., lower detection rate or higher false alarms. An alternative solution to this problem is to set an appropriate threshold. Therefore, in Section 3.1, the DGLR method is integrate with the RA-based threshold-setting algorithm for the purpose of fault detection of non-Gaussian processes.

3.1. DGLR with RA-Based Threshold Setting Algorithm

In the DGLR-based fault detection method, the solution of fault detection problem consists of two procedures:
  • Off-line training. Using the stacking data to identify the unknown parameters, i.e., the mean value E ( y ) and the covariance matrix Σ ;
  • On-line implementation. Detecting faults with on-line data.
In the first procedure, with N recorded data available, the data can be augmented and stacked in the following manner:
y = y k T y k 1 T y k N T y k + 1 T y k T y k N + 1 T y k + p T y k + p 1 T y k + p N T
where y k T denotes the recorded data at time k, and p + 1 represents the number of samples. The mean value E ( y ) and the covariance matrix Σ can be estimated according to Equations (7) and (8), respectively. For on-line implementation, first, we collect on-line measurement data y k + i = [ y k + i , y k + i + 1 , , y k + i + p ] T , i = 1 , , n , and then calculate
y ¯ n = 1 n i = 1 n y k + i , Δ y ¯ = y ¯ n E ( y )
For fault detection purposes, a test statistic in a maximum likelihood ratio-like form can be used as
J D G L R = n ( Δ y ¯ T Σ ^ 1 Δ y ¯ )
By summarizing the previous analysis, an extension of DGLR using RA-based threshold setting is proposed to detect underlying faults subject to non-Gaussian processes. The step-by-step procedure of the DGLR with RA-based threshold setting algorithm is illustrated in Algorithm 2.
Algorithm 2: DGLR-based fault detection method with RA-based threshold setting.
Off-line training
S1: Computation of
E ( y ) = 1 N i = 1 N y i , Σ ^ = 1 N 1 i = 1 N ( y i E ( y ) ) ( y i E ( y ) ) T
S2: Determine the corresponding thresholds J t h , n g using Algorithm 1 with a given significance level α , in which the statistic, J D G L R is estimated according to Equation (12);
On-line implementation
S3: Collect real-time measurement y k + i , i = 1 , , n and calculate
y ¯ n = 1 n i = 1 n y k + i , Δ y ¯ = y ¯ n E ( y )
S4: Build test statistic
J D G L R = n ( Δ y ¯ T Σ ^ 1 Δ y ¯ )
S5: Check the decision logic:
J D G L R > J t h , n g faulty J D G L R J t h , n g fault - free .
Remark 1.
The extension of PCA-based, DPCA-based, and PCA-based GLR fault detection methods with RA technique can be achieved according to the procedures in Algorithm 2. Hence, they are not presented in this paper.

3.2. Comparison among GLR-Based, DGLR-Based, PCA-Based, DPCA-Based, and PCA-Based GLR Fault Detection Methods

Considering the applications of PCA-based, DPCA-based, and PCA-based GLR fault detection methods, in this subsection, we briefly discuss the relationship among the three fault detection methods with the GLR-based one to distinguish their scope.
In practice, the direct application of the DGLR-based method may be unavailable due to numerical reasons, e.g., the invertibility of the estimated covariance matrix. This fact leads to the application of PCA-based methods, in which the SVD (singular value decomposition) of the estimated covariance matrix is the core. The principle of DPCA is similar to that of PCA, the following introduces the basic PCA technique [5]:
Σ ^ = P Λ P T = P p c P r e s Λ p c 0 0 Λ r e s P p c T P r e s T
where P p c = p 1 , , p γ R m × γ and P r e s = p γ + 1 , , p m R m × ( m γ ) consist of the loading vectors, known as the principal components and residual components, respectively; γ represents the number of principal components, Λ p c = diag λ 1 , , λ γ and Λ r e s = diag λ γ + 1 , , λ m contain the corresponding eigenvalues, satisfying γ 1 λ γ > > λ γ + 1 λ m .
For fault detection, the Hotelling’s test statistic can be calculated with a single, on-line measurement
J P C A = ( y E ( y ) ) T P p c Λ p c 1 P p c T ( y E ( y ) )
From the calculation formula of two test statistics, it is clear that the data normalization plays a central role. It not only provides us with an estimation of the covariance matrix in off-line training, but also delivers the required residual signal for fault detection.
If the principal component γ < m , then the matrix P p c R m × γ is rank deficient, i.e., rank ( P p c ) = γ < m . From the fault detection viewpoint, the matrix P p c is not ‘all pass’ for faults, that is, there exists f 0 such that P p c T f = 0 . This situation is caused by the artificial design of rank deficient matrix P p c . If the principal component γ equals m, then the J P C A statistic reduces to Hotelling’s T 2 statistic, that is, when γ = m
J P C A = ( y E ( y ) ) T P Λ 1 P T ( y E ( y ) ) = ( y E ( y ) ) T Σ ^ 1 ( y E ( y ) )
As introduced in [5], the test statistic J D G L R in the form (12) is also called Hotelling’s T 2 test statistic. It is evident that the test statistic in form (14) is equivalent with the one in (12) under the single on-line measurement case.
It is worthwhile noting that, recently, a PCA-based GLR fault detection method was proposed in [34], in which the PCA technique is only used for establishing the mathematical process model and the GLR test is used to evaluate the residual signal. Table 1 presents a comparison between them to clarify their relationship.
Actually, the numerical problem of Σ ^ 1 is rare due to the presence of process noise and measurement error. Hence, in this paper, we assume that Σ ^ 1 is available.
Remark 2.
The other use of PCA is for dimensionality reduction purpose. However, for fault detection, dimension reduction is not always crucial [32].

4. Simulated Examples

4.1. Fault Detection in Synthetic Data

In this example, a non-Gaussian dynamic process is first simulated by using Matlab. Then, the performance of the GLR-based, DGLR-based, PCA-based, DPCA-based, and PCA-based GLR fault detection methods is assessed through its application to detect faults in synthetic data.

4.1.1. Data Generation

Except for the fact that the noise sources are non-Gaussian, the model for data generation is the same as the one used in [38], which is given as:
y k = T q 1 u k + K q 1 d k + n k + f y , k ,
where
T q 1 = 0.05 q 3 1 0.95 q 1 0 0.7 q 3 1 0.3 q 1 0 0.02966 q 3 1 1.627 q 1 + 0.706 q 2 0.0627 q 6 1 0.937 q 1 0 0 0 0.235 q 5 1 0.765 q 1 0.5 q 2 1 q 1 + 0.25 q 2 0 0.5 q 5 0.4875 q 6 1 1.395 q 1 + 0.455 q 2 0 0 0.2 q 6 1 0.8 q 1
K q 1 = 1 0.1875 q 1 1 0.9875 q 1 1 0.1875 q 1 1 0.9875 q 1 1 0.1875 q 1 1 0.9875 q 1 1 0.1875 q 1 1 0.9875 q 1 T
y k and u k are the output and control input vectors, n k denotes the independent non-Gausian white noise, f y , k is the fault introduced in k-th sensor. The four controllers are given as
u k i = G i q 1 y k i , i = 1 , 2 , 3 , 4
where
G 1 q 1 = 3.2235 + 3.07 q 1 1 q 1 , G 2 q 1 = 0.6641 + 0.625 q 1 1 q 1 , G 3 q 1 = 0.6991 + 0.518 q 1 1 q 1 , G 4 q 1 = 0.444 + 3.70 q 1 1 q 1 .
The above model is used to simulate 1000 fault-free data samples. These data are used to estimate the required parameters for the three fault detection methods. The number of principal components is determined as four by using of the cumulative percent variance method. p is set to be five, which is determined as given in [20]. To validate the detection performance, three faults listed in Table 2 are introduced in this process.
After identifying the required parameters, one significant problem remaining in the training phase is to set the threshold. For threshold-setting, the sample number N = 2.65 × 10 4 is set using Equation (9) for ϵ = 0.01 , δ = 0.01 . The remaining steps are followed by Algorithm 2. Given a significance level of 0.01, Table 3 lists the thresholds used in this work.
In order to demonstrate the advantage of the proposed RA-based threshold setting, the J D G L R test statistic is used as an example. The result of this statistic is shown in Figure 1, where the red line represents the Gaussian assumption-based threshold and the green line represents the RA-based threshold. Figure 1 shows that the J D G L R test statistic is always below the threshold value, which means a zero false alarm rate. Unfortunately, zero FAR will lead to a lower fault detection rate (FDR). From Figure 1, we can see that the RA-based threshold makes the J D G L R statistic approach to a FAR of 0.01, which satisfies the given significance level. Due to the limited space, we only use this statistic as an example. Next, only the RA-based threshold is used for comparison purpose.

4.1.2. Comparing the Five Methods Using Faulty Data

The testing data set, which is simulated using the same model given in Equation (15), consists of 1000 data samples, which are completely independent from the training data. In this case, the possibility of only a single fault is considered. To assess the abilities of the various fault detection methods, Fault 2 was introduced at the sample of 400. To compare the performances of the GLR-based, DGLR-based, PCA-based, DPCA-based, and PCA-based GLR methods, the receiver operating characteristic (ROC) curves of the five methods are shown in Figure 2, which shows the FDR for different values of the fault alarm rate (FAR). The ROC curves provide a measure to compare the detection accuracy of all test statistics in three methods as well as their sensitivities to variations in the detection thresholds. Figure 2 shows that there is a trade-off between a high FDR and a low FAR. It can be seen that the PCA-based GLR test provides a higher FDR than the conventional PCA-based method. This fact is consistent with the conclusion in [34]. The J G L R statistic of GLR-based method has a similar detection performance as the PCA-based GLR test. Evidently, the DGLR-based method has the highest FDR than the other four test statistics. This clearly shows the advantages of the DGLR-based method over the other methods. Furthermore, the FDR performance of the DGLR-based method can be further improved by choosing the number of p. It should also be kept in mind that a large value of p can lead to a high detection delay. Therefore, there is also a trade-off between FDR and detection delay when determining an appropriate p.

4.2. Fault Detection in CSTR Process

4.2.1. Data Generation

To validate the proposed method, a CSTR process is used with several typical faulty scenarios. CSTR is widely used in chemical processes and the Matlab simulation model used in this study is similar to that used in [39]. The schematic of this CSTR is shown in Figure 3, where the reactor temperature T is controlled by manipulating the coolant flow rate Q c . As can be seen, it consists of three inputs ( C i , T i , and T c i ) and four outputs (C, T, T C , and Q c ). The corresponding model is given as follows:
d C d t = Q V C i C a k C + ν 1
d T d t = Q V T i T a Δ H r k C ρ C p b U A ρ C p V T T c + ν 2
d T c d t = Q c V c T c i T c + b U A ρ c C p c V c T T c + ν 3
where v i represents process noise, k is an Arrhenius-type rate constant, k = k 0 exp E R T . The parameters in the above model are listed in Table 4. Five typical faults are used to validate the proposed fault detection method, which have been descripted and listed in Table 5. b 0 , a 0 , T c , 0 , C 0 , and Q c , 0 in table are nominal values. Fault 1 and Fault 2 represent catalyst decay and heat transfer fouling, and Faults 3–5 are the sensor faults on each of the three measured variables.
In existing literature, there are several ways to check the non-Gaussianity. Among them, the probability plot of the variables is the most commonly used. We plot the quantiles of the A phase of converter current as an example. Figure 4 shows the distribution of the variable, where ‘+’ denotes the sample and the dotted line shows the locus of zero-mean samples which are normally distributed. It can be seen that the samples do not match the zero-mean normal distribution, that is, the samples from the example system obey the non-Gaussian distribution.

4.2.2. Comparing the Five Methods Using Faulty Data

The CSTR data consist of two blocks: the training and test data blocks. The normal operating data will be referred to as the training data. Then, parameters E ( y ) and Σ ^ used in the proposed method are estimated from the training data. After obtaining the necessary parameters, one remaining issue in the training phase is to determine an appropriate threshold. For threshold-setting, the sample number N = 2.65 × 10 4 is chosen by means of Equation (9) for ϵ = 0.01 , δ = 0.01 . Hence, the thresholds used in this work are determined as J t h , G L R = 16.20, J t h , D G L R = 56.60, J t h , P C A = 15.40, J t h , D P C A = 50.20 and J P G = 66.01 with a significance level of 0.01. The number of principal components and time lag p are set to be six and five, respectively. The total samples of each operation run in the test data block are 500 and the various faults are introduced only at sample 200. This means that for each of the faults, the process is fault-free for the first 200 samples before the system becomes abnormal at the introduction of the fault.
To demonstrate the effectiveness of the proposed method, the monitoring performance will be discussed. The detection sensitivity of a fault detection method is commonly quantified by calculating the FDR, which will be used for discussing the detection sensitivity of the proposed methods. The response of the fault detection method is commonly represented by the detection delay (DD), which is the time period it takes to detect a fault after occurrence of the fault. As the desired FAR is given for threshold-setting, the two indicators, FDR and DD, are used to assess the detection performance of fault detection methods. From the result discussion of Figure 2, the detection performance of the proposed method is compared with the GLR-based, PCA-based, DPCA-based, and PCA-GLR-based methods by using all faults described above. The superiority of the test statistic J D G L R over the other test statistics considered in this paper is shown in Table 6 with respect to FDR. Due to the non-Gaussian characteristic, the test statistic J D G L R with RA achieves the better performance with the higher FDR value than all faults compared. Detection delays of all test statistics are presented in Table 7. The unit of DD is the same as the sample interval. As shown in Table 7, the J D G L R with RA test statistics is able to detect most of these faults earlier than the other ones. This point alos demonstrates the advantage of the test statistics J D G L R .
For demonstrating the advantages of DGLR with the RA approach, the detection results of both methods for Faults 3 and 5 are shown in Figure 5 and Figure 6, respectively. The solid line represents the test statistic, the red line is the threshold based on a Gaussian assumption, and the blue line means the threshold based on the RA approach. In both figures, from top to bottom, there are J G L R test statistics of the GLR-based method, J D G L R of the DGLR-based method, J P C A of the PCA-based method, J D P C A of the DPCA-based method, and J P G of the PCA-based GLR method, respectively. FDRs of the GLR-based, DGLR-based, PCA-based, DPCA-based, and PCA-based GLR methods are 6.7%, 36.3%, 5.4%, 14.87%, and 21.54%, respectively. DDs of the GLR-based, DGLR-based, PCA-based, DPCA-based, and PCA-based GLR methods are 1, 2, 1, 2, and 2 min, respectively. All figures clearly show that all test statistics are able to detect the faults. In addition, the test statistics with the RA approach result in a higher FDR and smaller detection delay than the standard test statistics without RA. It should be noted that the performances of J D G L R test statistic with respect to FDR and DD are better than its companions in both faulty cases.

5. Conclusions

In this work, a DGLR-based method has been proposed for non-Gaussian dynamic processes, which combines the standard GLR method and the RA approach to iteratively learn the suitable threshold, which releases the assumption of Gaussian distributed variables. Furthermore, the DGLR-based approach has been compared with GLR-based, PCA-based, DPCA-based, and PCA-based GLR fault detection approaches to clarify the application scope of these methods. The major difference between them lies in the inverse of the estimated covariance matrix. In addition, the detection performance of the DGLR-based method has been compared with the aforementioned methods using a numerical simulation of a non-Gaussian process and the CSTR process. The comparison results show that the DGLR-based approach is better than the other methods. For instance, the average FDRs of the GLR-based, DGLR-based, PCA-based, DPCA-based, and PCA-based GLR approaches are 49.84%, 72.52%, 28.16%, 47.82%, 62.69% in CSTR process. The average DDs of the GLR-based, DGLR-based, PCA-based, DPCA-based, and PCA-based GLR methods are 13, 9, 12.8, 20.4, and 20.4 min in CSTR process. The J D G L R test statistics of the DGLR-based method shows the best detection performance compared with all other test statistics considered in this paper. Because the detection performance could be affected by kinds of factors, the robustness of the proposed method will be validated by using real data collected in industrial processes in the future.

Author Contributions

Conceptualization, X.P. and L.G.; methodology, Z.C.; software, L.G.; validation, L.G. and Z.C.; formal analysis, X.P.; investigation, X.P.; writing—original draft preparation, Y.J.; writing—review and editing, X.P.; visualization, Y.J.; supervision, X.P. and Z.C.; funding acquisition, Z.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China (#62173349, #U20A20186), in part by the science and technology innovation Program of Hunan Province in China (#2021RC4054).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationship that could have appeared to influence the work reported in this paper.

References

  1. Maleki, M.R.; Amiri, A.; Castagliola, P. Measurement errors in statistical process monitoring: A literature review. Comput. Ind. Eng. 2017, 103, 316–329. [Google Scholar] [CrossRef]
  2. Zhao, C.; Chen, J.; Jing, H. Condition-Driven Data Analytics and Monitoring for Wide-Range Nonstationary and Transient Continuous Processes. IEEE Trans. Autom. Sci. Eng. 2021, 18, 1563–1574. [Google Scholar] [CrossRef]
  3. He, Z.; Chen, Z.; Zhou, H.; Wang, D.; Xing, Y.; Wang, J. A visualization approach for unknown fault diagnosis. Chemom. Intell. Lab. Syst. 2018, 172, 80–89. [Google Scholar] [CrossRef]
  4. Huang, K.; Zhang, L.; Yang, C.; Gui, W.; Hu, S. Unified Stationary and Nonstationary Data Representation for Process Monitoring in IIoT. IEEE Trans. Instrum. Meas. 2022, 71, 1–12. [Google Scholar] [CrossRef]
  5. Ding, S.X. Data-Driven Design of Fault Diagnosis and Fault-Tolerant Control Systems; Springer: London, UK, 2014. [Google Scholar]
  6. Luo, H.; Yang, X.; Krueger, M.; Ding, S.; Peng, K. A plug-and- play monitoring and control architecture for disturbance compensation in rolling mills. IEEE/ASME Trans. Mechatron. 2018, 23, 200–210. [Google Scholar] [CrossRef]
  7. Basseville, M.; Nikiforov, I. Detection of Abrupt Changes: Theory and Application; Prentice-Hall: New York, NY, USA, 1993. [Google Scholar]
  8. Gao, Z.W.; Cecati, C.; Ding, S.X. A Survey of Fault Diagnosis and Fault-Tolerant Techniques, Part I: Fault Diagnosis with Model-Based and Signal-Based Approaches. IEEE Trans. Ind. Electron. 2015, 62, 3757–3767. [Google Scholar] [CrossRef] [Green Version]
  9. Chaouch, H.; Charfeddine, S.; Ben Aoun, S.; Jerbi, H.; Leiva, V. Multiscale monitoring using machine learning methods: New methodology and an industrial application to a photovoltaic system. Mathematics 2022, 10, 890. [Google Scholar] [CrossRef]
  10. Chen, H.; Jiang, B.; Lu, N.; Mao, Z. Deep PCA Based Real-Time Incipient Fault Detection and Diagnosis Methodology for Electrical Drive in High-Speed Trains. IEEE Trans. Veh. Technol. 2018, 67, 4819–4830. [Google Scholar] [CrossRef]
  11. Cao, Y.; Yuan, X.; Wang, Y.; Gui, W. Hierarchical hybrid distributed PCA for plant-wide monitoring of chemical processes. Control. Eng. Pract. 2021, 111, 104784. [Google Scholar] [CrossRef]
  12. Liu, Q.; Kong, D.; Qin, S.J.; Xu, Q. Map-Reduce Decentralized PCA for Big Data Monitoring and Diagnosis of Faults in High-Speed Train Bearings. IFAC-PapersOnLine 2018, 51, 144–149. [Google Scholar] [CrossRef]
  13. Peng, K.X.; Zhang, K.; Li, G.; Zhou, D.H. Contribution rate plot for nonlinear quality-related fault diagnosis with application to the hot strip mill process. Control Eng. Pract. 2013, 21, 360–369. [Google Scholar] [CrossRef]
  14. Zhang, K.; Hao, H.Y.; Chen, Z.W.; Ding, S.X.; Peng, K.X. A Comparison and evaluation of key performance indicator-based multivariate statistics process monitoring approaches. J. Process. Control 2015, 33, 112–126. [Google Scholar] [CrossRef]
  15. Yin, S.; Wang, G.; Gao, H. Data-driven process monitoring based on modified orthogonal projections to latent Structures. IEEE Trans. Control. Syst. Technol. 2016, 24, 1480–1487. [Google Scholar] [CrossRef]
  16. Chen, Z.W.; Ding, S.X.; Zhang, K.; Li, Z.B.; Hu, Z.K. Canonical correlation analysis-based fault detection methods with application to alumina evaporation process. Control Eng. Pract. 2016, 46, 51–58. [Google Scholar] [CrossRef]
  17. Chen, Z.; Yang, C.; Peng, T.; Dan, H.; Li, C.; Gui, W. A Cumulative Canonical Correlation Analysis-Based Sensor Precision Degradation Detection Method. IEEE Trans. Ind. Electron. 2019, 66, 6321–6330. [Google Scholar] [CrossRef]
  18. Chen, H.; Li, L.; Shang, C.; Huang, B. Fault Detection for Nonlinear Dynamic Systems with Consideration of Modeling Errors: A Data-Driven Approach. IEEE Trans. Cybern. 2022, 1–11. [Google Scholar] [CrossRef]
  19. Ge, Z.Q.; Song, Z.H.; Gao, F.R. Review of recent research on data-based process monitoring. Ind. Eng. Chem. Res. 2013, 52, 3543–3562. [Google Scholar] [CrossRef]
  20. Ku, W.F.; Storer, R.H.; Georgakis, C. Disturbance detection and isolation by dynamic principal component analysis. Chemom. Intell. Lab. Syst. 1995, 30, 179–196. [Google Scholar] [CrossRef]
  21. Wang, X.; Kruger, U.; Irwin, G.W. Process Monitoring Approach Using Fast Moving Window PCA. Ind. Eng. Chem. Res. 2005, 44, 5691–5702. [Google Scholar] [CrossRef]
  22. Chen, Q.; Wynne, R.; Goulding, P.; Sandoz, D. The application of principal component analysis and kernel density estimation to enhance process monitoring. Control Eng. Pract. 2000, 8, 531–543. [Google Scholar] [CrossRef]
  23. Zhang, Q.; Basseville, M. Statistical detection and isolation of additive faults in linear time-varying systems. Automatica 2014, 50, 2527–2538. [Google Scholar] [CrossRef]
  24. Kini, K.R.; Madakyaru, M. Monitoring multivariate process using improved Independent component analysis-generalized likelihood ratio strategy. IFAC Pap. 2020, 53, 392–397. [Google Scholar] [CrossRef]
  25. Tang, M.; Yang, C.; Gui, W. Fault detection based on cost-sensitive support vector machine for alumina evaporation process. Control Eng. China 2011, 18, 645–649. [Google Scholar]
  26. Jiang, Q.; Yan, X.; Lv, Z.; Guo, M. Independent component analysis-based non-Gaussian process monitoring with preselecting optimal components and support vector data description. Int. J. Prod. Res. 2014, 52, 3273–3286. [Google Scholar] [CrossRef]
  27. He, Q.P.; Wang, J. Statistics pattern analysis: A new process monitoring framework and its application to semiconductor batch processes. AIChE J. 2011, 57, 107–121. [Google Scholar] [CrossRef]
  28. Chen, Z.W.; Ding, S.X.; Peng, T.; Yang, C.H.; Gui, W.H. Fault Detection for Non-Gaussian Processes Using Generalized Canonical Correlation Analysis and Randomized Algorithms. IEEE Trans. Ind. Electron. 2018, 65, 1559–1567. [Google Scholar] [CrossRef]
  29. Bishop, C.M. Pattern Recognition and Machine Learning (Information Science and Statistics); Springer: Secaucus, NJ, USA, 2006. [Google Scholar]
  30. Chen, T.; Morris, J.; Martin, E. Probability density estimation via an infinite Gaussian mixture model: Application to statistical process monitoring. J. R. Stat. Soc. Ser. C (Appl. Stat.) 2006, 55, 699–715. [Google Scholar] [CrossRef] [Green Version]
  31. Odiowei, P.; Cao, Y. Nonlinear dynamic process monitoring using canonical variate analysis and kernel density estimations. IEEE Trans. Ind. Inform. 2010, 6, 36–45. [Google Scholar] [CrossRef] [Green Version]
  32. Gonzalez, R.; Huang, B.; Lau, E. Process monitoring using kernel density estimation and Bayesian networking with an industrial case study. ISA Trans. 2015, 58, 330–347. [Google Scholar] [CrossRef]
  33. Tschumitschew, K.; Klawonn, F. Incremental quantile estimation. Evol. Syst. 2010, 1, 253–264. [Google Scholar] [CrossRef]
  34. Harrou, F.; Nounou, M.N.; Nounou, H.N.; Madakyaru, M. Statistical fault detection using PCA-based GLR hypothesis testing. J. Loss Prev. Process. Ind. 2013, 26, 129–139. [Google Scholar] [CrossRef]
  35. Chamanbaz, M.; Dabbene, F.; Tempo, R.; Venkataramanan, V.; Wang, Q.G. Sequential randomized algorithms for convex optimization in the presence of uncertainty. IEEE Trans. Autom. Control 2016, 61, 2565–2571. [Google Scholar] [CrossRef] [Green Version]
  36. Faradonbeh, M.K.S.; Tewari, A.; Michailidis, G. Randomized algorithms for data-drivn stabilization of stochastic linear systems. In Proceedings of the 2019 IEEE Data Science Workshop (DSW), Minneapolis, MN, USA, 2–5 June 2019; pp. 170–174. [Google Scholar]
  37. Chen, Z.W. Data-Driven Fault Detection for Industrial Processes: Canonical Correlation Analysis and Projection based Methods. Ph.D. Thesis, University of Duisburg-Essen, Duisburg, Germany, 2016. [Google Scholar]
  38. McNabb, C.A.; Qin, S.J. Fault diagnosis in the feedback-invariant subspace of closed-loop systems. Ind. Eng. Chem. Res. 2005, 44, 2359–2368. [Google Scholar] [CrossRef]
  39. Pilario, K.E.S.; Cao, Y. Canonical variate dissimilarity analysis for process incipient fault detection. IEEE Trans. Ind. Inform. 2018, 14, 5308–5315. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Detection results of DGLR in fault-free case.
Figure 1. Detection results of DGLR in fault-free case.
Symmetry 14 01332 g001
Figure 2. ROC curves of the five fault detection methods.
Figure 2. ROC curves of the five fault detection methods.
Symmetry 14 01332 g002
Figure 3. The schematic diagram of CSTR.
Figure 3. The schematic diagram of CSTR.
Symmetry 14 01332 g003
Figure 4. Probability distribution of the variable.
Figure 4. Probability distribution of the variable.
Symmetry 14 01332 g004
Figure 5. Detection results of Fault 3.
Figure 5. Detection results of Fault 3.
Symmetry 14 01332 g005
Figure 6. Detection results of Fault 5.
Figure 6. Detection results of Fault 5.
Symmetry 14 01332 g006
Table 1. Comparison between (D)GLR-based, (D)PCA-based and PCA-based GLR fault detection methods.
Table 1. Comparison between (D)GLR-based, (D)PCA-based and PCA-based GLR fault detection methods.
MethodApplication ConditionsNumber of Test StatisticsParameters
GLR(DGLR)-based Σ ^ 1 is available J G L R and J D G L R m, n, E ( y ) and Σ ^
PCA(DPCA)-based Σ ^ 1 is unavailable J P C A and J D P C A m, γ , P p c , P r e s   E ( y ) and Σ ^
PCA-based GLRBoth cases, in which PCA model is available J P G m, γ , P r e s , E ( y ) and Σ ^
Table 2. Faults introduced in process.
Table 2. Faults introduced in process.
Fault IDsDescriptionValue of δ
1 y 1 , i = y 1 , 0 + δ 0.2
2 y 1 , i = y 1 , 0 + δ t 0.005
3 y 1 , i = y 1 , 0 + δ N (0, 0.04)
Table 3. Thresholds of all test statistics.
Table 3. Thresholds of all test statistics.
MethodTest StatisticsRA-Based ThresholdGaussian-Based Threshold
GLR J G L R J t h , G L R = 13.3011 20.0901
DGLR J D G L R J t h , D G L R = 22.6015 32.0001
PCA J P C A J t h , P C A = 9.0010 13.2775
DPCA J D P C A J t h , D P C A = 4.0105 4.1991
PCA-based GLR J P G J t h , P G = 10.4016 13.2775
Table 4. Constant parameters in CSTR model.
Table 4. Constant parameters in CSTR model.
ParameterDescription100.0Unit
QInlet flow rate150.0L/min
VTank volume10.0L
V c Jacket volume0.7L
Δ H r Heat of reaction 2.0 × 10 5 cal/mol
UAHeat transfer coefficient 7.0 × 10 5 cal/min/K
k 0 Pre-exponential factor to k 7.2 × 10 10 min−1
E R Activation energy 1.0 × 10 4 K
ρ , ρ c Fluid density1000g/L
C p , C p c Fluid heat capacity1.0cal/g/K
Table 5. Brief description of typical CSTR faults.
Table 5. Brief description of typical CSTR faults.
Fault IDsDescription of FaultsValue of δ
1 b = b 0 e x p ( δ t ) 0.001
2 a = a 0 + δ 1.4
3 T c = T c , 0 + δ 0.7
4 C = C 0 + δ 0.005
5 Q c = Q c , 0 + δ 5
Table 6. Detection results with respect to FDR.
Table 6. Detection results with respect to FDR.
Fault IDs J GLR J DGLR J PCA J DPCA J PG
179.74%85.38%16.39%61.16%86.49%
299.35%96.75%99.35%96.76%99.35%
354.34%95.77%8.03%15.85%85.53%
49.03%48.37%11.57%50.48%20.57%
56.75%36.36%5.46%14.88%21.54%
Table 7. Detection results with respect to DD.
Table 7. Detection results with respect to DD.
Fault IDs J GLR J DGLR J PCA J DPCA J PG
14234164747
222222
323274747
41841844
512122
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pan, X.; Gao, L.; Jiao, Y.; Chen, Z. A Dynamic GLR-Based Fault Detection Method for Non-Gaussain Dynamic Processes. Symmetry 2022, 14, 1332. https://doi.org/10.3390/sym14071332

AMA Style

Pan X, Gao L, Jiao Y, Chen Z. A Dynamic GLR-Based Fault Detection Method for Non-Gaussain Dynamic Processes. Symmetry. 2022; 14(7):1332. https://doi.org/10.3390/sym14071332

Chicago/Turabian Style

Pan, Xiaogang, Long Gao, Yuanyuan Jiao, and Zhiwen Chen. 2022. "A Dynamic GLR-Based Fault Detection Method for Non-Gaussain Dynamic Processes" Symmetry 14, no. 7: 1332. https://doi.org/10.3390/sym14071332

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop