Printer Friendly

Stability Monitoring of Batch Processes with Iterative Learning Control.

1. Introduction

Batch processes are the main way to produce the high value-added products such as fine chemicals industry and biopharming. However, in batch processes, sensor variables will exhibit strong within-batch autocorrelation as well as nontrivial batch-to-batch correlations [1]. Essentially, batch processes with dynamics along both time and batch directions are referred to as a two-dimensional (2D) system. The design of controller in batch processes is different from that in continuous ones. There are many ways to make sure of the stabilization of a feedback control loop when the process is operated continuously. Since the existence of 2D dynamic, the conventional feedback controller for continuous processes may not be effective for batch processes [2].

Recently, the iterative learning control (ILC) has received tremendous attention in batch processes [2]. It is a general technique used to improve the control performance of a system by refining control input from cycle to cycle. ILC of batch processes improves the quality of products by considering the use of previous batch information to control the current one. It not only does make sure of the system stabilization of every batch operation by using feedback controller, but also improves the tracking performance by using feed forward controller.

Since the controller performance is closely related to the quality and product yield of industrial products, some scholars have analyzed the ILC controller performance of batch processes. Chen and Kong introduced an optimal ILC index to evaluate the controller performance based on minimum variance [2]. Wei et al. proposed a new controller performance analysis (CPA) method based on linear quadratic Gaussian ILC under the condition that the system model is unknown [3]. However, the above methods focused mainly on the performance assessment of ILC controller. They did not study the stability of batch processes. Stability is a preoccupation of batch processes when the ILC is applied. The degradation of controller performance, actuator breakdown, or the unexpected interference all can make system stabilization worsen. Shi et al. defined the concept of batch convergence with considering the unknown disturbance and initial value [4]. Wang et al. introduced an ILC-PI controller and analyzed the robust asymptotically stable conditions of batch processes [5]. However, those stability research papers focused mainly on analyzing the concept of batch convergence. They did not show how to monitor the stability of the batch operation. Because of the intrinsically dynamic operations of batch processes, some stability problems still go undiagnosed for a prolong period of time. In this paper, stability monitoring of batch operation system with ILC is developed. Because of the large variation of batch processes in the operation condition during a batch run, the stability index during a batch run came into place for online stability monitoring. Siljak came up with a 2D system stability test method based on the Schwarz model and innerwise matrix [6]. Based on their methods, the paper proposed two kinds of stability indices along both time and batch directions.

Based on the 2D system theory, the closed-loop batch processes with ILC can be expressed 2D time series model [7]. The proposed method adopted a general 2D-ARMA model to express the ILC batch process with considering the uncertainty of the model and disturbance. However, it is not easy to identify that time series under close-loop conditions. Since there are many terms available in a 2D-ARMA, an effective fault monitoring method can be developed only if accurate 2D model identification is available. Aksasse and Radouane adopted the 2D Akaike information criterion (AIC) to identify the autoregressive model [8]. Yao and Gao proposed an automatic determination method based on the iterative stepwise regressions (SWR) to identify the 2D batch processes model and obtained a better identify result [9]. Recently, an adaptive least absolute shrinkage and selection operator (LASSO) had become one of the most popular approaches for automatic variable selection in linear regression [10]. Wang et al. compared the advantage and disadvantage of adaptive LASSO and SWR algorithm [11]. They indicated that the adaptive LASSO procedure in terms of order estimation is much better than the SWR and the SWR is slightly superior in terms of accuracy in values of true coefficients. The paper considered the advantage of those two methods. An integrated method of adaptive LASSO and SWR was adopted to identify the orders and coefficients of that 2D time series. Compared with the existing methods, the integrated method has better identification performance. Given the coefficients and orders of the identified model, the stability indices can be easily computed and monitored by SPC chart.

The rest of this paper is organized as follows: Section 2 will introduce the 2D-ARMA model and give the estimation of white noise. Section 3 will introduce our methodology including the derivation of the 2D-ARMA model of the batch process with ILC, the construction of innerwise matrix, and the definition of the stability indices. The effectiveness of the proposed approach is demonstrated in Section 4 by some simulations. The last section gives our conclusions.

2. 2D-ARMA Model

The 2D-ARMA model of a batch process with considering the influence of disturbance can be given by the following model [12]:

[mathematical expression not reproducible] (1)

where i represents the within-batch time index, k represents the batch-to-batch index, [] is the coefficient of the autoregressive (AR) part, and [b.sub.pq] is the coefficient of the moving average (MA) part. e(i,k) ~ N(0,[[sigma].sup.2]) is an 2D unobserved random white noise. M, N, P, Q are the orders of the 2DARMA model, respectively.

To estimate the value of [epsilon](i, k), we adopt a 2D long AR model to express the output y(i, k) of the batch system [11]:

[mathematical expression not reproducible] (2)

where M', N' [much greater than] M,N and the coefficient [c.sub.m'n'] and [??](i, k) of the above formulation can be obtained by linear regression.

The 2D-ARMA model (1) can be also considered as a discrete system with transfer function [13]:

[mathematical expression not reproducible] (3)

where [z.sup.-1.sub.1] ([z.sup.-1.sub.1] y(i, k) = y(i - 1, k)) is the back- shift operator for the within-batch time index i and [z.sup.-1.sub.2] ([z.sup.-1.sub.2] 1y(i,k) = y(i,k1)) is the back-shift operator for the batch-to-batch index.

3. The Stability Monitoring Scheme

3.1. The Batch Process with ILC. Consider a SISO system with performing repetitively an operation in a limited time T, called a batch run k. The feedback control structure of the batch system is shown in Figure 1. From Figure 1, the output of the batch process can be expressed as follows [2]:

y (i, k) = [G.sub.p] ([z.sup.-1.sub.1],k) u (i, k) + [G.sub.[epsilon]] ([z.sup.- 1.sub.1]) [epsilon] (i, k) i = 1,2, ... T; k = 1, 2, ... (4)

where T is the running time of each batch, [G.sub.p]([z.sup.-1.sub.1],k) and [G.sub.[epsilon]] ([z.sup.-1.sub.1], k) are the process and disturbance transfer function, respectively, and u(i, k) is the controller output.

According to the feedback control structure of each batch k, the controller output u(i, k) can be expressed as

??(??, ??) = [??.sub.??]([??.sup.?1.sub.1], ??) ([??.sup.sp] (??, ??) ? ??(??, ??)), (5)

where [x.sup.sp] (i,k) is the reference signal and [G.sub.c] ([z.sup.- 1.sub.1],k) is the feedback controller.

Substituting (5) into (4), the output can be expressed as

[mathematical expression not reproducible] (6)

The principle of ILC adopts the previous batch output y(i,k-1) and the reference signal [x.sup.sp](i, k-1) to control during the current batch k under the disturbance input [epsilon](i, k) [2]. The specific ILC structure is shown in Figure 2, and the update law can be expressed as follows:

[mathematical expression not reproducible] (7)

where [y.sup.sp](i) is the expected output of each batch and [G.sub.f]([z.sup.-1.sub.1], k) is the feedforward controller and adjusts the reference signal [x.sup.sp](i,k) by the output errors [epsilon](i,k - 1) = [y.sup.sp](i)-y(i,k-1).

Since [z.sup.-1.sub.2]y(i,k) = y(i,k - 1), the modified desired trajectory [x.sup.sp] (i, k) is given as

[mathematical expression not reproducible] (8)

Substituting (8) into (6), the following equation can be obtained:

[mathematical expression not reproducible] (9)

By adding and subtracting the set point [y.sup.sp](i), (9) can also be written as

[mathematical expression not reproducible] (10)

So the error e(i, k) can be written as a function of the previous batch error e(i, k - 1):

[mathematical expression not reproducible] (11)

According to Box et al.'s study [12], the above relation can be transformed into the following 2D-ARMA form:

[mathematical expression not reproducible] (12)

where r, s, v are the orders of G([z.sup.-1.sub.1]), H([z.sup.-1.sub.1]), N([z.sup.-1.sub.1]), respectively.

By moving H([z.sup.-1.sub.1])e(i, k - 1) term to the left, (12) can be written as

[mathematical expression not reproducible] (13)

And the ILC batch process can be considered as a discrete system with transfer function:

[mathematical expression not reproducible] (14)

3.2. Stability Analysis of Batch Process with ILC. According to the definition of bounded-input bounded-output (BIBO) stable, the 1D system is BIBO stable if all zeros A ([z.sup.-1.sub.1], 0) are not more than the unit circle. Consequently, the 2D polynomial A([z.sup.-1.sub.1], [z.sup.-1.sub.2]) is said to be BIBO stable if [6]

[mathematical expression not reproducible] (15)

Let us write A([z.sup.-1.sub.1], [z.sup.-1.sub.2]) as a polynomial in [z.sup.-1.sub.2] with coefficients that are polynomials in [z.sup.-1.sub.1]:

[mathematical expression not reproducible] (16)

In the pioneering paper [14], the concept of innerwise matrix has been used to analyze the stability of single dimensional systems. Kanellakis et al. extended it to 2D case and created the following innerwise matrix [[DELTA].sub.2N]([z.sup.-1.sub.1]) [13]:

[mathematical expression not reproducible], (17)

where [bar.d.sub.N] ([z.sup.-1.sub.1]) is the complex conjugate of [d.sub.N] ([z.sup.-1.sub.1]).

According to (16) and (17), the innerwise matrix [[DELTA].sub.2][z.sup.-1.sub.1] ([z.sup.-1]) of batch process with ILC can be constructed as

[mathematical expression not reproducible] (18)

where [??.sub.1]([??.sup.?1.sub.1]) [congruent] ??([??.sup.?1.sub.1]) = 1+[??.sub.1] [??.sup.?1.sub.1] + ... + [??.sub.??][??.sup.???.sub.1] and [??.sub.0]([??.sup.?1.sub.1]) ? ??([??.sup.?1.sub.1]) = [?.sub.0] + [?.sub.1][??.sub.?1.sub.1] + ... + [?.sub.??][??..sup.???.sub.1].

Siljak pointed that, for positivity checking of such an innerwise matrix, one requires the positivity checking of the matrix at one point, say at z = 1, and the positivity checking of the determinant for all [absolute value of z] = 1 [6]. Based on their study, we proposed two indices to analyze the stability of batch process with ILC:

R1 = max_positive_root_of [a([z.sup.-1.sub.1],u) = 0] [less than or equal to] 1, (19)

R2 = determination [[DELTA].sub.2] ([z.sup.-1.sub.1])] [greater than or equal to] 0 (20)

[for all] [z.sup.-1.sub.1],[absolute value of [z.sup.-1.sub.1]] = 1,

where R1 is the stability index in the within-batch time direction and R2 is the stability index in the batch-to-batch direction. The system will not be stable if the within-batch dynamic is not stable. Even the system is stable in the within-batch time direction if (19) is valid, the system can be unstable in the batch-to-batch if (19) and (20) are violated.

3.3. The Process of Stability Monitoring Scheme. In order to ensure the safe operation of batch processes with ILC, the paper proposed a real-time stability monitoring method, which contains two phases. Phase 1 is to establish control limits in normal period for each stability index. Phase 2 is to monitor each stability index online and to determine each stability index in control. The details are shown as follows.

Phase 1. (1) Collect the dynamic batch processes data in a normal period.

(2) Identify the 2D-ARMA model for each batch in the normal period by adopting an adaptive LASSO algorithm to identify the order of 2D-ARMA model and an SWR algorithm to determine the coefficients of 2D-ARMA model.

(3) Construct the innerwise matrix A2([z.sup.-1]) based on the identified coefficients of 2D-ARMA model and calculate the stability index R1 and R2.

(4) Establish Shewhart control charts for each stability indices.

Phase 2. (1) Collect the dynamic batch processes data at the current batch k.

(2) Identify the 2D-ARMA model for this batch k by adopting adaptive LASSO and SWR algorithm.

(3) Construct the innerwise matrix [[DELTA].sub.2]([z.sup.-1.sub.1]) and calculate the stability indices R1 and R2 for this batch k.

(4) Determine whether R1 and R2 for this batch k are in control.

4. Simulation

A batch process is given as follows:

[mathematical expression not reproducible] (21)

The feedback and feedforward controllers are [G.sub.c] = 1/[(1-0.9z.sup.-1] and [G.sub.f] = 0.5. Under those two controllers, the process can track the setpoint after seven batches, which is shown in Figure 3. In subsequent monitoring, a window with 30 batches, each batch has 1000 data samples and is used to estimate the 2D-ARMA model. Seventy moving windows are created by moving the window of 30 batches forward one batch at a time. The first 30 windows are used as "normal" data to estimate the control limits. The limits UCL = [mu] + 3 [sigma] and LCL = [mu] - 3[sigma] are the conventional control limits, which are widely used in SPC. Variations beyond these two limits indicate that there are substantial changes to the process. The system become unstable as variations exceeded the stability limits SL = 1 or SL = 0. The within-batch dynamic is unstable as the value of R1 is greater than 1 and the batch-to-batch dynamic is unstable as the value of R2 is smaller than 0. Two step changes [G.sub.p] = (0.6-0.3[z.sup.-1])/(1-0.7[z.sup.-1] [right arrow] [G.sub.p] = (0.45-0.3[z.sup.-1] 1)/(1-0.7[z.sup.-1] and Gp = (0.6-0.3[z.sup.-1])/(1-0.7[z.sup.-1] [right arrow] [G.sub.p] = (0.6[z.sup.-1] - 0.3[z.sup.-2])/(1 - 0.7[z.sup.-1] have been introduced at 31st window. Figures 4 and 5 show the output of different batches when those two step changes happened, respectively. It can be clearly see that the system is unstable when there is a process change Gp = (0.6-0.3[z.sup.-1])/(1-0.7[z.sup.-1] [right arrow] [G.sub.p] = (0.6[z.sup.-1] - 0.3[z.sup.-2])/(1 - 0.7[z.sup.-1].

Figure 6 shows the stability monitoring with a process change Gp = (0.6 - 0.3[z.sup.-1])/(1 - 0.7[z.sup.-1]) [right arrow] [G.sub.p] = (0.45 - 0.3[z.sup.-1])/(1 - 0.7[z.sup.-1]). From Figure 6, it can be clearly see that those two indices are not more than the stability limit although they are beyond the control limit at the 35th window, which indicated that the within-batch process and batch-to-batch dynamics are still stable. Figure 7 shows the stability monitoring with a process change [G.sub.p] = (0.6 0.3[z.sup.-1])/(1 - 0.7[z.sup.-1]) [right arrow] [G.sub.p] = (0.6[z.sup.-1] - 0.3z-2)/(1 - 0.7[z.sup.-1]). From Figure 7, R1 had become smaller when the process change occurred, which indicated that within-batch process is still stable. However, R2 was beyond the stability limit SL = 0, which indicated that the batch-to-batch dynamics suffer a substantial change leading to system instability.

5. Conclusion

In this paper, the batch process with ILC was derived as a 2DARMA model. An integrated technology of adaptive LASSO and SWR was adopted to identify the 2D-ARMA model. The paper analyzed the stability of ILC batch process based on innerwise matrix and proposed two stability indices. Finally, conventional SPC chart was used to monitor those stability indices. Some simulation results showed the effectiveness of the proposed method.

Competing Interests

The authors declare that they have no competing interests.


This work was supported by the National Natural Science Foundation of China under Grants nos. 61603347, 61603346, and 61603348 and the Doctor Startup Funds Foundation of Zhengzhou University of Light Industry under Grants nos. 2015BSJJ024 and 2015BSJJ026.


[1] N. Lu, Y. Yao, F. Gao, and F. Wang, "Two-dimensional dynamic PCA for batch process monitoring," AIChE Journal, vol. 51, no. 12, pp. 3300-3304, 2005.

[2] J. Chen and C.-K. Kong, "Performance assessment for iterative learning control of batch units," Journal of Process Control, vol. 19, no. 6, pp. 1043-1053, 2009.

[3] S. Wei, J. Cheng, and Y. Wang, "Data-driven two-dimensional LQG benchmark based performance assessment for batch processes under ILC," in Proceedings of the 9th IFAC Symposium on Advanced Control of Chemical Processes (ADCHEM '15), pp. 291-296, Whistler, Canada, June 2015.

[4] J. Shi, F. Gao, and T.-J. Wu, "Robust iterative learning control design for batch processes with uncertain perturbations and initialization," AIChE Journal, vol. 52, no. 6, pp. 2171-2187, 2006.

[5] Y. Wang, Y. Yang, and Z. Zhao, "Robust stability analysis for an enhanced ILC-based PI controller," Journal of Process Control, vol. 23, no. 2, pp. 201-214, 2013.

[6] D. D. Siljak, "Stability criteria for two-variable polynomials," Institute of Electrical and Electronics Engineers. Transactions on Circuits and Systems, vol. 22, no. 3, pp. 185-189, 1975.

[7] Y. Yao and F. R. Gao, "Batch process monitoring in score space of two dimensional dynamic Principal Component Analysis (PCA)," Industrial & Engineering Chemistry Research, vol. 46, no. 24, pp. 8033-8043, 2007.

[8] B. Aksasse and L. Radouane, "Two-dimensional autoregressive (2-D AR) model order estimation," IEEE Transactions on Signal Processing, vol. 47, no. 7, pp. 2072-2077, 1999.

[9] Y. Yao and F. R. Gao, "Statistical monitoring and fault diagnosis of batch processes using two-dimensional dynamic information," Industrial and Engineering Chemistry Research, vol. 49, no. 20, pp. 9961-9969, 2010.

[10] H. Zou, "The adaptive lasso and its oracle properties," Journal of the American Statistical Association, vol. 101, no. 476, pp. 1418-1429, 2006.

[11] Y. Wang, Y. Yao, Y. Zheng, and D. S. Wong, "Monitoring of within batch and batch-to-batch dynamics using adaptive LASSO," International Journal of System Control and Information Processing, vol. 1, no. 4, 2015.

[12] G. E. Box, G. M. Jenkins, G. Reinsel, and G. M. Ljung, Time Series Analysis: Forecasting and Control, Wiley Series in Probability and Statistics, John Wiley & Sons, Hoboken, NJ, USA, 5th edition, 2015.

[13] A. Kanellakis, S. Tzafestas, and N. Theodorou, "Stability tests for 2-D systems using the schwarz form and the inners determinants," IEEE Transactions on Circuits and Systems, vol. 38, no. 9, pp. 1071-1077, 1991.

[14] E. I. Jury, "Theory and applications of the inners," IEEE Transaction on Automatic Control, vol. 63, no. 7, pp. 1044-1069, 1971.

Yan Wang, Junwei Sun, Taishan Lou, and Lexiang Wang

College of Electronic and Information Engineering, Zhengzhou University of Light Industry, Zhengzhou 450002, China

Correspondence should be addressed to Yan Wang;

Received 14 November 2016; Accepted 5 January 2017; Published 31 January 2017

Academic Editor: Xavier Leoncini

Caption: Figure 1: The feedback control structure of each batch k.

Caption: Figure 2: The schematics of ILC with the feedback controller ([G.sub.c]([z.sup.-1.sub.1], k)) and the feedforward controller ([G.sub.f]( [z.sup.-1.sub.1], k)).

Caption: Figure 3: The controlled output.

Caption: Figure 4: The controlled output with a process change [G.sub.p] = (0.6- 0.3[z.sup.-1])/(1-0.7[z.sup.-1] [right arrow] [G.sub.p] = (0.45 - 03 [z.sup.-1])/{1 - 0.7[z.sup.-1]).

Caption: Figure 5: The controlled output with a process change [G.sub.p] = (0.6 0.3[z.sup.-1])/{1 - 0.7[z.sup.-1]) [right arrow][G.sub.p] = (0.6 [z.sup.-1] 0.3 [z.sup.-2])/(1 - 0.7[z.sup.-1]).

Caption: Figure 6: Stability monitoring with a process change [G.sub.p] = (0.6 0.3[z.sup.-1])/{1 - 0.7[z.sup.-1]) [right arrow][G.sub.p] = (0.45 - 0.3[z.sup.- 1])/{1 - 0.7[z.sup.-1]).

Caption: Figure 7: Stability monitoring with a process change [G.sub.p] = (0.6-0.3[z.sup.-1])/(1-0.7z [G.sup.f] = (0.6 [z.sup.-1] -0.3[z.sup.-2])/{1-0.7[z.sup.-1]).
COPYRIGHT 2017 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Wang, Yan; Sun, Junwei; Lou, Taishan; Wang, Lexiang
Publication:Advances in Mathematical Physics
Date:Jan 1, 2017
Previous Article:Application of Adjoint Data Assimilation Method to Atmospheric Aerosol Transport Problems.
Next Article:Quasi-Particles, Thermodynamic Consistency, and the Gap Equation.

Terms of use | Privacy policy | Copyright © 2020 Farlex, Inc. | Feedback | For webmasters