Multimode Process Monitoring Approach Based on Moving Window

Dec 15, 2017 - Due to the existence of multiple operating modes, traditional fault detection techniques are ill-suited for complex industrial processe...
4 downloads 7 Views 2MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Multi-mode Process Monitoring Approach Based on Moving Window Hidden Markov Model lin wang, ChunJie Yang, and Youxian Sun Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.7b03600 • Publication Date (Web): 15 Dec 2017 Downloaded from http://pubs.acs.org on December 22, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Industrial & Engineering Chemistry Research is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Multi-mode Process Monitoring Approach Based on Moving Window Hidden Markov Model Lin Wang, Chunjie Yang*, Youxian Sun State Key Laboratory of Industrial Control Technology, Institute of Industrial process Control, Zhejiang University, Hangzhou, Zhejiang 310027, P. R. China E-mail: [email protected] Abstract: Due to the existence of multiple operating modes, traditional fault detection techniques are ill-suited for complex industrial process. Although there are more and more literatures concerning this problem, only a few of them are based on hidden Markov model (HMM). But there is no exploration to concern the unknown mode in industrial process based on it. This artical proposes a novel process monitoring approach based on moving window HMM (MVHMM) for real-time multi-mode process monitoring with unknown mode. First, a hidden Markov model is built by training set. Instead of just considering the posterior probability of one single sample, moving window is introduced to utilize the independency of samples for improving the accuracy of on-line mode identification. Besides, an MVHMM-based threshold statistic is defined to identify the unknown mode. And various known modes which include stable modes and transitions are seperated based on Viterbi algorithm. Second, a new monitoring scheme is developed for fault detection of each mode. The effectiveness of the proposed approach is validated by Tennessee Eastman (TE) chemical process and a numerical simulation example. Key Words: Hidden Markov model, Multi-mode process, Transition process, Fault detection

1 INTRODUCTION Effective process monitoring and fault diagnosis are critical to improve plant safety, profitability and reliability.1, 2 Due to the wide utilization of sensors, a large amount of data have been collected, which promoted rapid development of data-based industrial monitoring techniques in the past several decades, such as partial least squares (PLS) and principal component analysis (PCA).3–7 Typical data-based multivariate statistical process monitoring (MSPM) method like PCA assumes that the correlation between the variables has the same process characteristics, while actual industrial production process often does not satisfy this assumption. Most of the productive processes have more than one operating mode, which means that there are multiple stable working points. And the correlation between variables of diverse stable working points is different. There are many reasons for multi-mode conditions: changes in raw materials, external environments, process loads, or equipment depreciation will alter the process conditions. Each of the above cases will cause multi-mode process. And the dominant variables and process characteristics are not always the same in each operating mode. The corresponding statistical properties of the process like mean, variance and correlation of variables are often distinct. However, the traditional MSPM methods have not well considered mode shifts. In recent years, a lot of scholars have done researches on multi-mode process. Srinvasan first estimated the Euclidean distance of the sample points to identify the mode of the sample, and then use dynamic PCA (DPCA) to approximate the sample clustering.8 Choi proposed an adaptive MSPM approach to update the mean and covariance of samples to achieve the purpose of updating the model.9 Liu developed an adaptive Takagi Sugeno fuzzy model in the subspace of PCA in order to fit the time-varying characteristics of This work is supported by the 863 Program (Grant 2012AA041709) and the National Natural Science Foundation of China (Grant No. 61290321)

multi-mode process.10 He and Wang used k-nearest neighbor (kNN) rule in multi-mode process.11 Yu combined finite Gaussian mixture model (GMM) with Bayesian inference mechanism, and successfully applied that method to continuous stirred tank heater (CSTH) process and TE process.12 Ge developed a just-in-time-learning (JITL) strategy for establishing the current monitoring model.13 The Bayesian method is widely used to distinguish different modes of samples.14 Ng chose a PCA with the least differences as the most suitable sub-model by calculating the differences between T 2 and SP E in each PCA model.15 Ge proposed a Bayesian regularization method based on probabilistic PCA (PPCA) for multi-mode process monitoring.16 Tan proposed a novel monitoring method based on the similarity of data feature, which constructs different models to obtain most of the characteristics of process variables.17 Zhu offered an adjoined multi-ICA-PCA model for multi-mode process monitoring.18 Xie combined the moving window and the Gaussian mixture model very well. The method can characterize the multi-mode process well on the one hand, and on the other hand, it can accurately identify the changes of the process and realize the online self-adaptation monitoring.19 Compared with traditional multi-mode process monitoring methods, HMM could efficiently deal with inherent uncertainty and dynamic characteristic in industrial process as a typical dynamic Bayesian network (DBN) model.20 While, dynamics and uncertainty are two main features in the process of practical production, which should be paid widespread attention.21 So HMM has been gradually used to handle the problems of multi-mode process monitoring in recent years. Yu proposed an HMM-based process monitoring approach for multi-mode and nonlinear process.22 Rashid proposed a new HMM based on ICA approach for fault detection.1 Ning proposed a HMM-SPA approach to monitor multi-mode process.23 Wang proposed a probability ratio strategy based on HMM to identify the transitions in multi-mode process.24, 25 To the best of our knowledge, there is no research on the

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 2 of 11

Page 3 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Based on the Baum-Welch algorithm, the parameters λ = {A, B, π} of HMM could be decided by given sequence of hidden states, the corresponding observation sequence and the inial state probability distribution.27, 28 It is important to note that all the modes of training samples are known in advance, which means that the unknown mode only appears in the testing data. At the same time, observation probability density distribution matrix bi (Yt ) = P(Yt |qt = Si ) can be gained. The element in the matrix means the probability that current observation belongs to each operating mode. For an observation, the probability which it belongs to its own mode is far greater than the probability that it belongs to other modes. YT ∗n = {Y1 , Y2 , · · · , YT }

(7)

where T means the sequence number of observations in training set and n represents the number of variables in training set. Once HMM has been trained, the variance σt2 of posterior probability of each observation in training set can be calculated as follows: L P

σt2

= var P(Yt |qt = Si ) =

i=1

Here observation sequence in the kth moving window for testing is as follows Ykw∗n = {Yk−w+1 , Yk−w+2 , · · · , Yk }

(11)

where w signifies the width of the moving window and k signifies the moving window index. It should be noted that when moving window is used to intercept the dynamic information of the process data, it is necessary to select the appropriate window length according to the dynamic characteristics of the specific process. For a fast dynamic response process, the length of the window should be chosen shorter, which is conducive to the rapid mode identification. For some relatively slow dynamic response processes, a longer window should be selected, which could contain more complete process dynamic information. This is beneficial to accurate mode identification. After identifying unknown mode and known mode, the observation probability density distribution matrix is utilized to separate current observations into suitable operating modes. The corresponding optimum state in the kth moving win dow Q∗w = qk−w+1 , qk−w+2 , · · · , qk can be obtained by maximizing the conditional probability P (Q|Y, λ) according to Viterbi algorithm.27

(P(Yt |qt = Si ) − Pt )2 L

. (8)

Q∗w = arg max P (qk−w+1 , qk−w+2 , · · · , qk | Qw

(12) Yk−w+1 , Yk−w+2 , · · · , Yk , λ)

Pt =

PL

i=1

P(Yt |qt = Si ) L

(9)

As long as the current observation belongs to one of the known modes, the variance of posterior probability is a relatively large value. On the contrary, if the observation does not belong to any known mode, the probability which it belongs to each known mode is all a small value. Hence, its variance is tiny. Based on the above analysis, the variance of posterior probability of each observation can be an indication to identify the unknown mode. Because all the modes in training set are known, we could get the predetermined threshold Pt ∗ by minimizing the variance of posterior probability of all training samples. If the variance of posterior probability is under Pt ∗ , it may be an indication that there is an unknown mode present in the process. The threshold can be decided as follows: Pt ∗ = α ∗ min(σt2 )

(10)

where α is a tolerant parameter of Pt . Actually, it should be a particularly small value to ensure the effectiveness of threshold Pt ∗ . Instead of just considering the posterior probability of one single sample in testing set, we introduce moving window to utilize the independency of samples. At each moment, the moving window contains the dynamic measurements of each variable. By accurately analyzing the dynamic characteristics of the measured data within the window, the dynamic characteristic information of the operation process is obtained. This could avoid the shortcomings of insufficient information which is contained in a single sample, and facilitate the accurate on-line mode recognition of the samples.

where Q∗w expresses the optimal hidden state sequence for k the corresponding observation sequence Yw∗n . Here HMM is used for estimating the sequence of operation condition. Each hidden state in the sequence represents a specific operation condition. Hence, any monitoring sample can be divided into an appropriate operation condition. With the identified operation type from HMM, the operating mode is confirmed. 3.2 On-line Fault Detection Using a New Monitoring Index Scheme When there are multiple operational modes in the production process, there is no sudden change from a steady operating mode to another stable operating mode.29 Just considering the steady mode is not comprehensive, it is also necessary to take into account the gradual transition between the stable operating modes. The process characteristics of the transition mode are different from those of the steady mode. The transition mode will show a dynamic gradient trend, which is not only reflected in the change of the process variable, but also in the changes of the correlations between process variable. In the view that the process characteristics of transitions and the stable operation modes are entirely different, a new monitoring scheme is developed for multi-mode process monitoring with transitions after mode recognition. In the scheme, two kinds of monitoring indexes are build separately for process monitoring. They can be automatically switched according to the mode of sample. The new monitoring index scheme combines two kinds of probability information: multivariate Gaussian probability density and likelihood probability.

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 4 of 11

Page 5 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 6 of 11

Page 7 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

Table 4: Process Fault for TEP Fault No. IDV(1) IDV(2) IDV(3) IDV(4) IDV(5) IDV(6) IDV(7) IDV(8) IDV(9) IDV(10) IDV(11) IDV(12) IDV(13) IDV(14) IDV(15) IDV(16-20) IDV(21)

Fault Description Step in A/C feed ratio, B composition constant (stream 4) Step in B composition, A/C ratio constant (stream 4) Step in D feed temperature (stream 2) Step in reactor cooling water inlet temperature Step in condenser cooling water inlet temperature A feed loss (step change in stream 1) C header pressure loss-reduced availability Random variation in A+B+C feed composition (stream 4) Random variation in D feed temperature (stream 2) Random variation in C feed temperature (stream 2) Random variation in reactor cooling water inlet temperature Random variation in condenser cooling water inlet temperature Slow drift in reaction kinetics Sticking reactor cooling water valve Sticking condenser cooling water valve Unknown The valve for steam 4 was fixed at the steady state position

Table 3: Six Different Operating Modes of TEP

500

mode 1 2 3 4 5 6

400

300

MGDLLP

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 11

200

G/H mass ratio 50/50 10/90 90/10 50/50 10/90 90/10

Production rate (steam 11) 7038 kg/h G and 7038 kg/h H 1408 kg/h G and 12669 kg/h H 10000 kg/h G and 1111 kg/h H Maximum Maximum Maximum

100

0

-100 0

20

40

60

80

100

120

140

160

180

200

sample

Fig. 13: Monitoring result of HMM in case 4 of numerical example Table 2: Monitored Variables of TEP Variable No. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

Variable Description A feed (stream 1) D feed (stream 2) E feed (stream 3) A and C feed (stream 4) Recycle flow (stream 8) Reactor feed rate (stream 6) Reactor temperature Purage rate (stream 9) Product separator temperature Product separator pressure Product separator underflow (stream 10) Stripper pressure Stripper temperature Compressor power Reactor cooling water outlet temperature Separator cooling water outlet temperature

In Fig. 15, 16, the monitoring results of fault 4 in mode 3 are displayed seperately by PCA and proposed approach. As can be seen from the Fig. 15, obviously it is completely un-

able to monitor the occurrence of fault in the process. But in Fig. 16, which represents the detection ability of the developed approach, the fault could be detected around the 270th sample. According to the result which is gained from TEP, the developed method could effectively identify the operating mode for multimode process with transitions even when fault happened. Then we can compare its monitoring performance with the PCA in TEP simulation cases, which validates the developed method is superior to PCA in detecting the faults under multi-mode process with transitions. Table 5: Result of mode recognition of TEP method MVHMM

5

errors 33

missclassification rate 0.825%

CONCLUSIONS

A new process monitoring approach based on MVHMM for multi-mode process with unknown mode is developed in this paper. Compared with the existing works, the proposed approach takes into account the identification of unknown operating mode in industrial process for the first time. And moving window is introduced to utilize the independency of samples for improving the accuracy of on-line mode identification. Then, a new monitoring scheme is suggested after mode recognition. For stable mode, monitoring index MDGLLP which combined the multivariate Gaussian probability density and likelihood probability is used for fault detection. And monitoring index MGD is used for fault detection of transition mode. Compared with the traditional

ACS Paragon Plus Environment

Page 9 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[5] Qin, S.J. Statistical process monitoring: basics and beyond. J. Chemometr. 2003, 17(8-9), 480-50. [6] Kano, M.; Hasebe S. J.,; Hashimoto, I.;Ohno,H . A new multivariate statistical process monitoring method using principal component analysis. Comput. Chem. Eng. 2001, 25(7-8), 1103–1113. [7] Wang, X.; Kruger, U.;Lennox, B. Recursive partial least squares algorithms for monitoring complex industrial processes. Control Eng. Pract. 2003, 49(6), 969–976. [8] Srinivasan, R.; Wang, C.; Ho, W. K.; Lim, K. W. Dynamic principal component analysis based methodology for clustering process states in agile chemical plants. Ind. Eng. Chem. Res. 2004, 43(9), 2123-2139. [9] Choi, S. W.; Martin, E. B.; Morris, A. J.; Lee, I. B. Adaptive multivariate statistical process control for monitoring timevarying processes. Ind. Eng. Chem. Res. 2006, 45(9), 31083118. [10] Liu, J. L. Modeling a large-scale nonlinear system using adaptive Takagi Sugeno fuzzy model on PCA subspace. Ind. Eng. Chem. Res. 2007, 46(3), 788-800. [11] He, Q. P. Fault detection using the k-Nearest neighbor rule for semiconductor manufacturing processes. IEEE T. Semiconduct. M. 2007, 20(4), 345-354. [12] Yu, J.; Qin, S. J. Multimode process monitoring with Bayesian inference-based finite Gaussian mixture models. AIChE J . 2008, 54(7), 1811-1829. [13] Ge, Z. Q.; Song, Z. H. Online monitoring of nonlinear multiple mode processes based on adaptive local model approach. Control Eng. Pract. 2010, 56(11), 2838-2849. [14] Liu, J. L. Fault detection and identification using modified Bayesian classification on PCA subspace. Ind. Eng. Chem. Res. 2009, 48(6), 3059-3077. [15] Ng, Y. S.; Srinivasan, R. An adjoined multi-model approach for monitoring batch and transient operations. Comput. Chem. Eng. 2009, 33(4), 887-902. [16] Ge, Z. Q.; Song, Z. H. Mixture Bayesian regularization method of PPCA for multimode process monitoring. AIChE J. 2008, 16(12), 1427-1437. [17] Tan S.,; Wang, F. L.; Peng, J.; Chang, Y. Q.; Wang, S. Multimode process monitoring based on mode identification. Ind. Eng. Chem. Res. 2012, 51(1)), 374-388. [18] Zhu, Z. B.; Z. Song, H.; Palazoglu, A. Process pattern construction and multi-mode monitoring. J. Process Control. 2012, 22(1), 247-262. [19] Xie, X.; Shi, H. B. Dynamic multimode process modeling and monitoring using adaptive Gaussian mixture models. Ind. Eng. Chem. Res. 2012, 51(15), 5497-5505. [20] Murphy, K.P. Dynamic Bayesian networks: representation, inference and learning, UC Berkeley. 2002, 13, 303-306. [21] Chen, X. R.; Ge, Z. Q. Switching LDS-based approach for process fault detection. Chemometr. Intell. Lab. 2015, 146, 169-178. [22] Yu, J. B. Hidden Markov models combining local and global information for nonlinear and multimodal process monitoring. J. Process Control. 2010, 20(3), 344-359. [23] Ning, C.; Chen, M. Y.; Zhou, D. H. Hidden Markov modelbased statistics pattern analysis for multimode process monitoring: an index-switching scheme. Ind. Eng. Chem. Res. 2014, 53(27), 11084-11095. [24] Wang, F.; Tan, S.; Shi,H. B. Hidden Markov model-based approach for multimode process monitoring. Chemometr. Intell. Lab. 2015, 148(15), 51-59. [25] Wang, F.; Tan, S.; Yang, Y. W.; Shi, H. B. Hidden Markov model-based fault detection approach for a multimode process. Ind. Eng. Chem. Res. 2016, 55(16), 4613-4621. [26] Yan, Z. B.; Huang, C. C.; Yao, Y. Semi-supervised mix-

ture discriminant monitoring for chemical batch process. Chemometr. Intell. Lab. 2014, 134, 10-22. [27] Rashid, L. R. A Tutorial on the hidden Markov models and selected applications in the speech recognition. Proc. IEEE 1989, 77, 257-286. [28] Baum, L. E.; Petrie, T. Statistical inference for probabilistic functions of finite state Markov chains. Ann. Math. Stat. 1966, 37(6), 1554-1563. [29] Srinivasan, R.; Doan, X. T. Online monitoring of multi-phase batch processes using phase-based multivariate statis- tical process control. Comput. Chem. Eng. 2008, 32(1-2), 230-243. [30] Ge, Z. Q.; Zhao, L. P.; Yao, Y.; Song, Z. H.; Gao, F. R. Utilizing transition information in online quality prediction of multiphase batch processes. J. Process Control. 2012, 22(3), 599-611. [31] Deng, X. G.; Tian, X. M. Nolinear process fault pattern recognition using statistics kernel PCA similarity factor. Neurocomputing. 2013, 121, 298-308. [32] Shang, J.; Chen, M. Y.; Ji, H. Q.; Zhou, D. H. Recursive transformed component statistical analysis for incipient fault detection. Automatica. 2017, 80, 313-327.

ACS Paragon Plus Environment

Page 10 of 11

Page 11 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54

Industrial & Engineering Chemistry Research

Off-line learning

On-line process monitoring

Training data

Testing data in the moving window

Build HMM

Calculate Pt

Calculate Threshold Pt*

Pt < Threshold

Yes

Unknown mode

No Mode recognition

Stable mode

Monitoring index LLP

Monitoring index MGD

Fault detection MGDLLP ACS Paragon Plus Environment

Transition mode

Monitoring index MGD

Fault detection MGD