Dynamic Data Reconciliation for Enhancing Performance of Minimum

Sep 8, 2016 - ABSTRACT: Healthy controllers are required in order for the control systems to maintain a high level of performance. In past research, m...
0 downloads 10 Views 1MB Size
Subscriber access provided by CORNELL UNIVERSITY LIBRARY

Article

Dynamic Data Reconciliation for Enhancing Performance of Minimum Variance Control in Univariate and Multivariate Systems Zhengjiang Zhang, and Junghui Chen Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.6b02532 • Publication Date (Web): 08 Sep 2016 Downloaded from http://pubs.acs.org on September 13, 2016

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 57

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

Dynamic Data Reconciliation for Enhancing Performance of Minimum Variance Control in Univariate and Multivariate Systems

Zhengjiang Zhanga and Junghui Chenb*

a

College of Physics and Electronic Information Engineering, Wenzhou University, Wenzhou 325035, People’s Republic of China b

Department of Chemical Engineering, Chung-Yuan Christian University, Chung-Li,Taoyuan, Taiwan, 320, Republic of China

Revision Submitted to Industrial & Engineering Chemistry Research

August, 2016

                                                              

To whom all correspondence should be addressed. Tel.: +886-3-2654107. Fax: +886-3-2654199. Email: [email protected]  1 

 

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

ABSTRACT: Healthy controllers are required in order for the control systems to maintain a high level of performance. In the past research, minimum variance control (MVC) played a crucial role as a benchmark in performance monitoring because of the attractive theoretical and computational properties associated with it. Since the influence of measurement errors has not been explicitly considered in the MVC theory in stochastic control systems, this paper first analyzes the influence of measurement errors on the control performance of MVC in both univariate and multivariate systems. And then the dynamic data reconciliation (DDR) methods are proposed and combined in the procedure of MVC/multivariate MVC (MMVC) to decrease the influence of measurement errors and to enhance the control performance. Considering both random measurement errors and gross errors, the effectiveness of MVC/MMVC combined with DDR on the variances of process outputs is illustrated by two simulated cases, including both univariate and multivariate systems. An actual pilot scale experiment with two inputs and two outputs is used to demonstrate the effectiveness of MMVC combined with DDR. Results of simulation and experiment show that MVC/MMVC combined with DDR can efficiently enhance the control performance.

KEYWORDS: dynamic data reconciliation; measurement noise; minimum variance control; stochastic system

1. INTRODUCTION In order to keep the plant operation stable, safe and up to the environmental standard, 2   

ACS Paragon Plus Environment

Page 2 of 57

Page 3 of 57

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

robust performance is required in control systems. Process control engineers and instrumentation

technicians

have

the

daunting

task

of

maintaining

and

trouble-shooting control systems. However, a typical modern industrial process operation consists of hundreds or even thousands of control loops. The poorer performing control loops may go unnoticed indefinitely or for a long period of time. A survey conducted by Desborough and Miller (2002) showed that only 16% of the controllers operate at their best to give the optimal performance and about 36% of the loops operate in a manual mode.1 It is not feasible for the operating personnel to manually monitor all the control loops in a plant. An automatic performance monitoring system should be used to continuously evaluate product quality and the production levels of a plant. The performance monitoring system is built upon measurable plant data that are reworked (usually statistically) to provide a clear and simple picture of the operating states of a plant. Then, one needs a means of comparison to evaluate the current performance.

The research on the control system performance assessment can be traced back to the 1960s and 1970s. Astrom (1967) and DeVries et al. (1978) carried out the fundamental research.2, 3 In the literature of control loop performance assessments, minimum variance control (MVC) had been widely used as a reference bound on the achievable performance since Harris (1989) proposed it.4 MVC is a simple, practical control strategy resulting from the application of the linear stochastic control theory.5 The basic idea of MVC is to express the minimum variance prediction of the plant output in the terms of the control input so that the predicted output can be driven to follow a desired output by solving a simple linear equation for appropriate control action.6 This groundbreaking research has received increasing attention in the process 3   

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

control field. Then several research efforts on the evaluation of the control performance have been made. The use of MVC has been extended to multivariate control systems based on the multivariate interpretation of the delay term, known as the interactor matrix, by Harris et al. (1996), Huang et al. (1997), Yan et al. (2015).7-9 Ko and Edgar (2001) extended MVC to the issue of the cascade control system performance assessment.10 They also investigated the minimum variance performance based PID controller in single loop case and developed an iterative solution for the calculation of the best achievable (minimum variance) PID control performance later.11. In recent years, MVC has been also developed for time-variant process systems,12,13 model predictive controllers,14,

15

PCA-based minimum variance

performance,16 decentralized controllers,17 adaptive optics systems,18 stochastic and deterministic control performance of batch processes,19 and nonlinear multivariate systems.20, 21 Shardt et al. (2012) gave a good review of the current techniques for performance assessment.22

One issue with using a minimum variance controller as a standard is that there are only a few industrial applications of MVC. Hugo (2006) concluded that there are three reasons why MVC is not always applicable in practice23: a) MVC contains no move suppression, so the resultant controller may be overly aggressive; b) MVC assumes that the process is perfectly known and is linear; c) Minimum variance controllers explicitly employ a full ARIMA model of the disturbance, as opposed to the limited disturbance model of industrial controllers.

As described by Hugo (2006), 23 the first two points are not significant deterrents to 4   

ACS Paragon Plus Environment

Page 4 of 57

Page 5 of 57

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

MVC while the last reason is more problematic, because the ARIMA model of the disturbance does not explicitly describe the influence of measurement errors in the sensor devices.

In most of the research papers, the measurement noise (also called random errors), which are small fluctuations of the measured variable from the true value and can be electrically induced or process-induced2, has been included in the original MVC through the concept of innovation process. However, the MVC theory for stochastic control systems does not explicitly consider the influence of measurement errors (containing 2 types of errors: random errors and gross errors) presented in actual sensor devices. When measurements are corrupted by random variations of the environment, they are said to be affected by random measurement errors. Gross errors in measurements usually occur for many different reasons, such as human errors, instrumental errors, fraudulent behavior and faults in systems. Gross errors can be quite large. Both errors have a detrimental effect on process control which requires accurate estimated states and measurements. In that situation, the MVC benchmark which is an indication of the inherent optimum would not be trustable. In the actual process system, the output signals are usually measured by sensor devices and the measured signals are then the feedback of the input signals of the controller for comparison. Since the sensor devices are not ideal and error-free, it is necessary to consider the influence of measurement errors on the control performance and to develop some corresponding dynamic data reconciliation (DDR) methods to enhance the control performance. The DDR methods have not been completely developed in the area of the control performance assessment yet.

5   

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

Page 6 of 57

Shown in this paper are several developments on MVC as follows: a) With additional consideration of measurement errors in the sensor devices, the influence of measurement errors on the control performance of MVC in both univariate and multivariate systems is analyzed. b) In order to decrease the influence of measurement errors on MVC and enhance the control performance, DDR methods are proposed and combined in the procedure of MVC. c) Considering both random measurement errors and gross errors, the effect of MVC combined with DDR on the variances of the process output is illustrated.

The remainder of this paper is organized as follows: MVC in stochastic univariate and multivariate systems is revisited in the next section. The influence of measurement noise on the control performance of MVC in univariate and multivariate systems is then analyzed in Section 3. The DDR methods are proposed and combined in the procedure of MVC to enhance the control performance in Section 4. In Section 5, the proposed MVC with DDR methods are used in three case studies to illustrate its effectiveness. Those case studies contain both univariate and multivariate systems. Finally, conclusions are drawn.

2. REVISIT

OF

MVC

IN

STOCHASTIC

UNIVARIATE

AND

MULTIVARIATE SYSTEMS Performance benchmarks are necessary to determine what the optimum operating state of the process is. This optimum operating state contains the inherent variance that cannot be compensated for by controllers. The operating optimum is then used as a benchmark to do comparative evaluation of the current operation. There are various 6   

ACS Paragon Plus Environment

Page 7 of 57

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

methods and most of them originate from the minimum variance theory developed by Harris (1989).4 The advantage of MVC is that no extra perturbations are necessary to determine the performance; only routine closed loop operating data are used. In this section, the MVC benchmark concept in univariate and multivariate systems would be reviewed.   (Figure 1 goes here.)   2.1 MVC in stochastic univariate system Consider the stochastic univariate feedback control system shown in Figure 1. In order to determine what the controller algorithm should be, consider a general controlled autoregressive moving average model of the closed loop response which can be described by:

y (t )  z

k

B( z 1 ) D( z 1 ) u (t )  d (t ) A( z 1 ) A( z 1 )

(1)

where u(t ) is the manipulated input, and d (t ) is a random external disturbance 2 which is assumed to be zero mean and Gaussian white noise of variance  . y

refers to the controlled variable or the process output. The scaling should be chosen so that y  0 when the process is at the desired set-point or steady state. Control will not be completely free from variance seeing that the manipulated variable takes k th time instant to affect the process output, y . Thus, MVC is used to minimize the variance of the output at time t  k given all the information up to time t ,

7   

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

 y2  E  y(t  k )2 

Page 8 of 57

(2)

E indicates the statistical expectation. The cost function in Eq.(2) is the equation that describes the deviation of the output from its desired value, so it contains the variance to be minimized by the controller.

Since z  k in Eq.(1) represents a k -step delay in the control signal, the disturbance term D ( z 1 ) can be decomposed into two parts, one related to the past time and the other related to the future time.

D( z 1 )  A( z 1 ) F ( z 1 )  z  k G( z 1 )

(3)

Then the output is simplified as a controllable variance part of the output signal and an uncontrollable part caused by time delay of the process.

 B( z 1 ) F ( z 1 )  G( z 1 ) y (t  k )  Fd (t  k )   u (t )  y(t )  1 1 D( z )  D( z ) 

(4)

The task is now to minimize the cost function (Eq.(2)) by reducing the feedback dependent term represented by the terms in the square bracket of Eq. (4). The terms in the square bracket need to be zero by changing the manipulated input ( u ).

B( z 1 ) F ( z 1 ) G ( z 1 )  u ( t ) y (t )  0 D( z 1 ) D( z 1 )

  8   

ACS Paragon Plus Environment

(5)

Page 9 of 57

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

Now the objective function in Eq.(2) to be minimized for the value of the output at the

k th time instant in the future can be represented by

 

2 MVC y

 E  ( Fd (t  k )) 2 

(6)

  One can determine what variance of the uncontrollable part in the output is and the variance will be the theoretical lower bound performance benchmark. The benchmark is an indication of the inherent optimum set by the process design and equipment. With the inherent optimum capability defined in a control loop, the detection of the poor control performance is a relatively simple task. What one needs to do is just to compare the normal variance of the output signal with the benchmark.

2.2 MVC in stochastic multivariate system The stochastic feedback control system with multiple inputs and multiple outputs (MIMO) is shown in Figure 1. The MIMO linear time-invariant stationary stochastic process model can be described by

y (t )  A( z 1 )u(t )  B( z 1 )d(t )

(7)

where A ( z 1 ) and B( z 1 ) are proper, rational transfer function matrices for the process plant and disturbances respectively; y(t ) and u(t ) are the output vector and the input vector separately; d(t ) represents a white noise vector of disturbances with zero mean and covariance matrix Σ d .

9   

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

Page 10 of 57

The difficulty in the multivariate controller performance assessment is the factorization of the time delay matrix, which is known as the interactor matrix. Assume that A ( z 1 ) is a full rank transfer function matrix. There is a unique, non-singular, unitary interactor matrix ( D ), and

 ( z 1 )  K lim D( z 1 ) A ( z 1 )  lim A 1

z 1  0

z 0

(8)

 ( z 1 ) is the delay-free transfer function matrix of A ( z 1 ) , and K is a full where A

rank constant matrix. Premultiplying both sides of Eq.(7) by z  k D( z 1 ) , where k is the maximum delay order of interactor matrix D( z 1 ) , yields

 ( z 1 )u(t )  N( z 1 )d(t ) y (t )  z  k A

(9)

where y (t )  z  k D( z 1 )y (t ) and N( z 1 )  z  k D( z 1 )B( z 1 ) . Huang and Shah (1999) have found that a unitary interactor matrix is an optimal factorization of time delays for multivariate systems in terms of minimum variance control and controller performance assessment.24 The disturbance term N ( z 1 ) can be decomposed into two parts, one related to the past time and the other related to the future time,

N( z 1 )  F( z 1 )  z  k R( z 1 )

(10)

Like the univariate system, the output is simplified as an uncontrollable part of the output signals and a controllable part because there is time delay of the process.

10   

ACS Paragon Plus Environment

Page 11 of 57

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

y (t )  F( z 1 )d(t )  F( z 1 )B( z 1 )1 A( z 1 )u(t )  z  k R( z 1 )B( z 1 )1 y(t )

(11)

After the feedback dependent, the square bracket term in Eq.(11), is reduced to zero, minimum variance control of the closed-loop response can be obtained.

 

2 MMVC y

 E  y (t )T y (t )   trace{Cov  Fd (t ) }

(12)

where trace is the sum of the diagonal elements of a matrix. It is shown that multiplying a unitary interactor matrix to the output variables does not change the quadratic measure or the variance of the corresponding output variables. Thus, the benchmark of the multivariate control system is obtained. These existing results with minimum variance control as the benchmark did not explicitly consider the influence of measurement errors. However, according to practical experiences, there is no sensor measurement without measurement noise. Measurement noise is the variation in the sensor reading that does not correspond to changes in the process and can be caused by background electrical interference, mechanical vibration and process fluctuations. For cases in which there is a relatively large amount of noise, the impact of the measurement noise on the performance of control systems will be discussed in the next section.

3. INFLUENCE OF MEASUREMENT NOISE ON MVC Measured process data inevitably contain some inaccurate information because measurements are obtained by imperfect instruments, but the bound of the control loop performance assessment discussed before was estimated from the data without 11   

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

Page 12 of 57

considering the influence of measurement noise. With the influence of measurement noise, the feedback signals are no long the actual output signals ( y(t )   for univariate systems and y(t ) for multivariate systems), and the measured signals contain the inherent variability associated with the measurement noise. The influence of measurement errors that degrades the control performance of MVC in univariate systems and of MMVC in multivariate systems would be separately discussed in this section.

(Figure 2 goes here.)

3.1 MVC in the stochastic univariate system with measurement noise Under the consideration of the measurement noise for the univariate system shown in Figure 2, the raw process measurement ( ym (t ) ) can be expressed as

ym (t )  y(t )   (t )

(13)

where  (t ) is the measurement noise assumed to be normally distributed (  (t ) ~ N (0,  2 ) ) and  2 is the variance of measurement noise.

After

substituting

Eq.(13)

into

Eq.(4)

and

regrouping

the

terms

into

feedback-invariant and feedback-dependent terms, the process model can be described as

12   

ACS Paragon Plus Environment

Page 13 of 57

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

y (t  k ) 

B( z 1 ) F ( z 1 ) G( z 1 ) u ( t ) ( ym (t )   (t ))  Fd (t  k )  D( z 1 ) D( z 1 )

 B( z 1 ) F ( z 1 )  G( z 1 ) G( z 1 )   =Fd (t  k )  ( t ) u ( t ) ym (t )    1 1 1 D( z ) D( z )  D( z ) 

(14)

Therefore, the output variance of the minimum variance controller satisfies the following equality

 

2 MVC y

2   G ( z 1 )    E  ( Fd (t  k ))   E   (t )   1  D ( z )   2

(15)

  Equality in Eq.(15) is obtained when the minimum variance controller is used. The minimum variance controller for this case can be obtained by solving the following equation,

B( z 1 ) F ( z 1 ) G ( z 1 )  u ( t ) ym (t )  0 D( z 1 ) D( z 1 )

(16)

Note that the first term on the right-hand side of Eq.(15) is generated by the process disturbance, and the second term is generated by the measurement noise. Compared with Eq.(6), Eq.(15) indicates that the output variance under MVC in the univariate system is increased because of the influence of measurement noise.

3.2 MVC in the stochastic multivariate system with measurement noise Like the univariate system, the measured signals y m (t ) in the multivariate system shown in Figure 2 can also be expressed as 13   

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

y m (t )  y (t )  ε(t )

Page 14 of 57

(17)

where ε(t) represents an uncorrelated random noise vector at the point t with zeros mean and the covariance matrix Σ ε . After Eq.(17) is substituted into Eq.(11), feedback-invariant and feedback-dependent terms can be yielded. The multivariate process model can be described as

y (t )  F( z 1 )d(t )  z  d R( z 1 )F( z 1 )1 z  d D( z 1 )ε(t )  ( z 1 )u(t )  R( z 1 )F( z 1 )1 z  d D( z 1 )y (t )   z  d  A m 

(18)

Therefore, the output variance of the minimum variance controller satisfies the following equality.

 

2 MMVC y

 trace{Cov  Fd (t ) }  trace{Cov  R ( z 1 )F ( z 1 ) 1 D( z 1 )ε (t ) }

(19)

  Equality in Eq.(19) is obtained when the minimum variance controller is used. Note that the first term on the right-hand side of Eq.(19) is generated by the process disturbances, and the second term is caused by the measurement noise. Compared with Eq.(12), Eq.(19) indicates that the variances of outputs under MMVC are also increased because of the influence of measurement noise.

14   

ACS Paragon Plus Environment

Page 15 of 57

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

4. MVC

COMBINED

WITH

DDR

TO

ENHANCE

CONTROL

PERFORMANCE When the measurements containing measurement errors are used by the controller, the performance of MVC would be affected. Bai demonstrated that the DDR was an effective tool to decrease propagation of measurement errors inside control loops.25 DDR can be used to decrease measurement errors inside control loops. With the uncertainty ranges of measurements determined by the measurement errors of measuring instruments as well as instrument degradation, 26 the DDR uses dynamic process models to reconcile the process measurements. Therefore, it can be combined with MVC to enhance the control performance. In this section, the effectiveness of the DDR based MVC (DDR-MVC) is proposed to overcome the degradation of the control performance in the presence of significant measurement noise.

4.1.DDR-MVC in a stochastic univariate system For univariate systems, at the time point t   the process output in (Eq.(1)) can be decomposed into two parts. One is the predicted part ( yˆ (t ) ) from the past input measurements; the other, the unpredicted part (  (t ) ) which is the inherent variability associated with the process disturbance,

y (t ) 

B( z 1 ) D( z 1 ) u t  k  d (t ) ( ) A( z 1 ) A( z 1 )   

(20)

 (t )

yˆ ( t )

where  (t )   represents the model prediction error with the Gaussian distribution (  (t ) ~ N (0, 2 ) ) and  2 is the corresponding variance. 15   

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

Page 16 of 57

(Figure 3 goes here.)

The stochastic univariate feedback control system with DDR is shown in Figure 3. According to Bayes’ rule, the posterior distribution of the reconciled output ( yr (t ) ) given the measurements ( ym (t ) ) and the predicted output ( yˆ (t ) ) is proportional to the product of two likelihood functions of the measured output and the predicted output,

p ( yr (t ) ym (t ), yˆ (t ))  L ( ym (t ) yr (t )) L ( yˆ (t ) yr (t )) p ( yr (t ))

(21)

where L ( ym (t ) yr (t )) and L ( yˆ (t ) yr (t )) are the likelihood functions. Based on the assumptions about the measurement noise and the model predictions, the likelihood functions can be shown as

L( ym (t ) yr (t )) 

 1 ( ym (t )  yr (t ))2  exp    2 2  2 

(22)

L( yˆ (t ) yr (t )) 

 1 ( yˆ (t )  yr (t )) 2  exp    2 2  2 

(23)

1

1

where  2 and  2 are the variances of the measured error and the predicted error, respectively. p( yr (t )) is a priori probability for the reconciled signal yr (t ) and it is a constant. Thus, the posterior distribution becomes

p( yr (t ) ym (t ), yˆ (t )) 

 1  ( ym (t )  yr (t ))2 ( yˆ (t )  yr (t ))2   1  exp     (24) 2 2 (2 )2   2    16 

 

ACS Paragon Plus Environment

Page 17 of 57

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

The estimation of yr (t ) based on a principle similar to the maximum likelihood is the maximum posterior distribution,

yr (t )  arg max p ( yr (t ) ym (t ), yˆ (t ))

(25)

p ( yr (t ) ym (t ), yˆ (t )) is proportional to L( ym (t ) yr (t )) L ( yˆ (t ) yr (t )) ,

y r (t )

yr ( t )

Since

minimizes the augmented likelihood L( ym (t ) yr (t )) L( yˆ (t ) yr (t )) . The reconciled signal yr (t ) is

yr (t )  yˆ (t )  (1   2 2 ) 1 ( ym (t )  yˆ (t ))

(26)

The estimation error of DDR can be given

 (t )  yr (t )  y (t )

(27)

Thus, with DDR, the model output is

B( z 1 ) F ( z 1 ) G( z 1 ) y (t  k )  u (t )  ( yr (t )   (t ))  Fd (t  k ) D( z 1 ) D( z 1 )  B( z 1 ) F ( z 1 )  G( z 1 ) G ( z 1 ) ( t ) u ( t ) yr (t )    Fd (t  k )     1 1 1 D( z ) D( z )  D( z ) 

(28)

Like Eq.(15), the output variance of the minimum variance controller satisfies the following equality 17   

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

2  G( z 1 )    (t )    E ( Fd (t  k ))   E  1  D( z )  

 

2 MVC y

2

Page 18 of 57

(29)

Using Eq.(26), the variance of  (t ) can be calculated by

var( (t ))  (  2   2 )1

(30)

The detailed derivations of the equations for the SISO system and for the MIMO system are similar; to save space, only equations for the MIMO system are derived in Appendix A. Eq.(30) indicates that the variance (  (t ) ) of the reconciled error yr (t ) is less than those of measured error  (t ) and model prediction error  (t ) . Thus, when Eq.(29) is compared with Eq.(15), the output variance under DDR-MVC is decreased since the variance of  (t ) is less than that of the measurement noise  (t ) . When the variance of the measurement noise (  ) increases significantly, the variance of the model prediction error ( ) will give the main contribution to the variance of

 (t ) . Thus, the influence of the process disturbance and measurement instrument bias, outliers, and slow drifts can be reduced efficiently by DDR. In the SISO case, since only one signal enters the DDR block, DDR can be considered as a kind of digital filter.

4.2. DDR-MVC in a stochastic multivariate system

In a multivariate system, at the time point t   the process output in (Eq.(9)) can be decomposed into two parts, 18   

ACS Paragon Plus Environment

Page 19 of 57

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

y (t )  A( z 1 )u(t )  B( z 1 )d(t )   yˆ ( t )

(31)

δ (t )

where δ(t ) is  normally distributed ( N (0, Σ δ ) ). Σ δ

is the corresponding

covariance matrix, which can be obtained by 

Σδ  Cov(B( z 1 )d(t ))

(32)

The stochastic feedback control MIMO system with DDR can be shown in Figure 3. According to Bayes’ rule, DDR uses both the measurement and the process model information to derive the reconciled signals y r (t ) ,

p(y r (t ) y m (t ), yˆ (t ))  L(y m (t ) y r (t )) L(yˆ (t ) y r (t )) p(y r (t ))

(33)

where p ( y r (t ) y m (t ), yˆ (t )) represents the posterior knowledge based on both the measurement information y m (t ) and the process model information yˆ (t ) . L ( y m (t ) y r (t )) and L ( yˆ (t ) y r (t )) are the likelihood functions. Based on the

assumptions about the measurement noise and the model predictions, the likelihood functions follow that

L(y m (t ) y r (t )) 

1 (2 )

N /2

Σε

1/2

1 exp( (y r (t )  y m (t ))T Σε 1 (y r (t )  y m (t ))) 2

19   

ACS Paragon Plus Environment

(34)

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

L(yˆ (t ) y r (t )) 

1 (2 ) N /2 Σδ

1/2

Page 20 of 57

1 exp( (y r (t )  yˆ (t ))T Σδ 1 (y r (t )  yˆ (t ))) 2

(35)

where Σ ε and Σ δ are the covariances of the measured error and the predicted error, respectively. p(y r (t )) is a priori probability for the reconciled signals y r (t ) and it is a constant. Thus, the posterior distribution becomes,

p( y r (t ) y m (t ), yˆ (t )) 

1 (2 ) Σε Σδ N

1/2

 1 exp  ( y r (t )  y m (t ))T Σε 1 (y r (t )  y m (t ))  2

(36)



( y r (t )  yˆ (t ))T Σδ 1 ( y r (t )  yˆ (t )) 

The estimation of y r (t ) , based on a principle similar to maximum likelihood is the maximum posterior distribution,

y r (t )  arg max p ( y r (t ) y m (t ), yˆ (t )) y r (t )

(37)

Since the p ( y r (t ) y m (t ), yˆ (t )) is proportional to L ( y m (t ) y r (t )) L ( yˆ (t ) y r (t )) , y r (t ) minimizes the augmented likelihood L ( y m (t ) y r (t )) L ( yˆ (t ) y r (t )) . The reconciled signals y r (t ) are

y r (t )  yˆ (t )  ( Σ ε 1  Σδ 1 ) 1 Σε 1 ( y m (t )  yˆ (t ))  yˆ (t )  K ( y m (t )  yˆ (t ))

where K  ( Σ ε 1  Σ δ 1 ) 1 Σ ε 1 . The estimation error of DDR can be given by 20   

ACS Paragon Plus Environment

(38)

Page 21 of 57

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

γ (t )  y r (t )  y (t )  δ(t )  K  ε (t )  δ(t ) 

(39)

Thus, under DDR, the model output is

y (t )  F( z 1 )d(t )  z  d R( z 1 )F( z 1 )1 z  d D( z 1 ) γ(t )  ( z 1 )u(t )  R( z 1 )F( z 1 )1 z  d D( z 1 )y (t )   z  d  A r 

(40)

Like Eq.(19), the output variance of the minimum variance controller satisfies the following equality

 

2 MVC y

 trace{Cov  Fd(t )}  trace{Cov G( z 1 )F( z 1 )1 D( z 1 )γ (t ) }

(41)

Using Eq.(39), the variance of γ(t ) can be calculated by

Vr  Cov ( γ (t ))  Cov (δ(t )  K (ε (t )  δ(t )))  ( Σε 1  Σδ 1 ) 1

(42)

The detailed derivations of Eqs.(38), (39) and (42) can be found in Appendix A. Eq.(42) indicates that the covariance matrix of the reconciled errors ( γ(t ) ) is a function of Σ ε (the covariance matrix of the measured error) and Σ δ (the covariance matrix of the model prediction error). Since both Σ ε and Σ δ are positive definite symmetric matrices, the trace of covariance matrix Vr

are less

than that in Σ ε and in Σ δ , which means that the reconciled signals are more precisely estimated than the measured signals. If MMVC uses the reconciled signals 21   

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

as the feedback signals instead of the measured signals, the variances of outputs will be decreased. Therefore, the DDR based MMVC (DDR-MMVC) can efficiently enhance the control performance.

The control structure in Figure 1 is widely used for the benchmarking purpose (like MVC). It is a typical feedback control structure used in the actual process. The traditional MVC does not explicitly consider the influence of measurement errors presented in actual sensor devices. The measurement errors would affect the variances of process outputs. In that situation, the MVC benchmark which is an indication of the inherent optimum would not be trustable. In this paper, the influence of measurement errors on the control performance of MVC is considered, so the DDR method combined in the procedure of MVC is proposed to reduce the influence of measurement errors. An additional DDR filter shown in Figure 3 is used to improve the control performance, but the feedback control structure is still the same as the original MVC benchmark of the feedback control structure. This means that the tight control performance of the proposed DDR-MVC would be gotten. The control structure which contains an additional DDR filter may be far from the common practice, but modern controllers on microprocessors and digital computers can be implemented easily.

5. CASE STUDIES To illustrate the effectiveness of the proposed methods developed in Section 4, three case studies are presented. For easy explanation, two simulation cases that consist of the given dynamic system with the measurement noise are used because they illustrate how to calculate the DDR-MVC and DDR-MMVC performance assessment indices 22   

ACS Paragon Plus Environment

Page 22 of 57

Page 23 of 57

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

and how to enhance the indices when the controlled system contains the measurement noise. Then in the third case, the data from the real experiment of the feedback control loop are presented to help readers compare the DDR-MMVC based performance assessment method with the conventional MVC based performance assessment method in terms of the performance improvement.

Case study 1: MVC in a stochastic univariate feedback control system The case study based on Hugo (2006)23 is selected to explain the performance assessment procedure developed in Section 4.1. The stochastic univariate feedback control system considered is:

0.08z 3 1 y(t )  u(t )  d (t ) 1 1  0.92 z 1  z 1 0.08 z 3 (1  z 1 ) 1  0.92 z 1 u t d (t ) ( )   (1  0.92 z 1 )(1  z 1 ) (1  0.92 z 1 )(1  z 1 )

(43)

where d (t ) is the zero mean Gaussian random variable, which has a standard deviation of   0.01 . The parameters of the PI controller ( be k1 

k1  k2 z 1 ) are selected to 1  z 1

1 0.92 and k 2  , both of which are derived from the control law of 3*0.08 3* 0.08

MVC.23

Two sets of conditions are used to demonstrate the performance

improvement of this SISO control system based on the proposed DDR-MVC:

Condition 1: Measurements with two different sizes of measurement noise In the first condition, the measurement noise (  (t ) ) is assumed to be normally distributed (  (t ) ~ N (0,  2 ) ). Two different sizes of measurement noise with 23   

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

Page 24 of 57

  0.025 and   0.05 are compared.

From Eq. (43), the predicted part ( yˆ (t ) ) and the unpredicted part (  (t ) ) can be derived as follows,

(1  0.92 z 1 )(1  z 1 ) y(t )  0.08z 3u (t )  (1  0.92 z 1 )d (t ) y(t )  1.92 z 1 y(t )  0.92 z 2 y(t )  0.08 z 3u(t )  (1  0.92 z 1 )d (t )

(44)

y(t )  1.92 y(t  1)  0.92 y(t  2)  0.08u(t  3)  (1  0.92 z 1 )d (t )     (t )

yˆ ( t )

If the standard deviation of the disturbance is given, the standard deviation of the model prediction error ( ) can be calculated by

 2  var (1  0.92 z 1 )d (t )  = 1  0.922   2

(45)

(Figure 4 goes here.)

Thus, from Eq.((26)), the reconciled output under MVC and the corresponding estimation error of DDR can be calculated (Eq.((27))). When   0.025 , The results of MVC based on measured signal and reconciled signal are shown in Figure 4. The standard deviations of the process outputs with MVC based on the measured signal and the reconciled signal are 0.0252 and 0.0224, respectively. The results demonstrate that the variance of the process outputs decreases using DDR.

When   0.05 , the results of MVC based on measured signal and reconciled signal are shown in Figure 5. The standard deviations of the process outputs with MVC 24   

ACS Paragon Plus Environment

Page 25 of 57

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 measured signal and the reconciled signal are 0.0367, and 0.0277, respectively. The results show that the variance of the process output with MVC based on measured signal is sensitive to the measurement noise. While the standard deviation of measurement noise is increased to be 0.05, the variance of the process output with MVC based on measured signal increases obviously, from 0.0252 to 0.0367. However, the variance of the process output with MVC based on reconciled signal increases a little, from 0.0224 to 0.0277, which demonstrates that the performance of MVC can be enhanced by DDR.

(Figure 5 goes here.)

Condition 2: Measurements with noise and gross errors The outliers are added to the measurements at time step 100 and 300. The magnitude of outliers is equal to 3.0. The bias is introduced into the measurements from time step 100 to 200. The magnitude of bias is equal to 0.05. Slow linear drift happens from time step 500 to 600 with a gradient of 0.0005. Like the previous discussion in Condition 1, the results of MVC based on measured signal and reconciled signal are shown in Figure 6. The standard deviations of the process outputs with MVC based on measured signal and reconciled signal are 0.0325 and 0.0266, respectively. With the influence of measurement instrument bias, outliers, and slow drifts, the variance of process output with MVC based on measured signal increases obviously, from 0.0252 (Condition 1) to 0.0325 (Condition 2). The variance of process output with MVC based on reconciled signal increases a little, from 0.0224 to 0.0266. The results indicate that the influence of measurement noise can be reduced efficiently by DDR. And the performance of MVC can be improved by DDR and it is robust to the measurement instrument bias, outliers, and slow drifts. 25   

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

(Figure 6 goes here.)

Condition 3: Comparison results with Monte Carlo simulations In the last condition, Monte Carlo simulations are used to demonstrate the reliable comparison results. The exponential filter, which is by far the most commonly used first order filter in industrial applications27, is used for comparisons. The measurement with only random errors (measurement noise) is considered firstly, and then measurement with both noise and gross errors is considered as the second case of measurement. The simulation is run 500 times repetitively in each case of measurement. The expected standard deviations of measured signal, filtered signal and reconciled signal are provided in Table 1 to illustrate the effectiveness of DDR.

(Table 1 goes here.)

As can be seen in Table 1, the expected standard deviations of both filtered signal and reconciled signal are decreased compared with the expected standard deviation of measured signal. The exponential filter can decrease the expected standard deviation of measured signal by 12.00% (from 0.0342 to 0.0301) with measurement noise, while the DDR can decrease the expected standard deviation of measured signal by 16.08% (from 0.0342 to 0.0287) with measurement noise. When measurement with both noise and gross errors is considered, the exponential filter can decrease the expected standard deviation by 13.00% (from 0.0354 to 0.0308), while the DDR can decrease the expected standard deviation by 17.80% (from 0.0354 to 0.0291). The results indicate that the influence of measurement noise and gross errors can be 26   

ACS Paragon Plus Environment

Page 26 of 57

Page 27 of 57

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

reduced much more efficiently by DDR. The effectiveness of DDR is more obvious comparing with that of exponential filter.

Case study 2: MVC in a stochastic multivariate feedback control system The case study for the multivariate system comes from Huang and Shah (1997).28 The system involves two manipulated variables and two controlled variables, which are widely used in literature.5, 29 The process and disturbance transfer function matrices are defined as follows:

 z 1 1  0.4 z 1 A( z 1 )    0.3z 1  1  0.1z 1

4 z 2  1  0.1z 1  z 2  1  0.8 z 1 

1   1  0.5 z 1 B( z 1 )    0.5 1  0.5 z 1

0.6  1  0.5 z 1   1  1  0.5 z 1 

(46)

The set points are assumed to be zeros for both outputs. The disturbances are two-dimensional normal-distributed white noise with the covariance matrix Σd  I . The system has a general interactor matrix which can be determined as

0.9578z 0.2873z  D( z 1 )   2 2  0.2873z 0.9578z 

(47)

The MVC controllers are designed according to the two single loops without considering the interactions5, 29:

27   

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

 0.5  0.2 z 1  1  0.5 z 1 C( z 1 )    0  

   0.25  0.2 z 1  1 1  (1  0.5 z )(1  0.5z ) 

Page 28 of 57

0

(48)

The vector of the measurement noise ε(t) is assumed to be normally distributed with covariance matrix Σ ε  I . MMVC based on measured signals and based on reconciled signals are studied in this example. Two sets of conditions are used to demonstrate the performance improvement of this MIMO control system based on the proposed DDR-MMVC.

Condition 1: Measurements with measurement noise With the given covariance matrix of the disturbance and the model prediction error ( Σ δ in Eq.(32)), the standard deviation of the model prediction error ( Σ δ ) is given by

Σδ  Cov (B( z 1 )d(t ))

(49)

Thus, from Eq.(38), the reconciled output under MVC can be calculated and the corresponding estimation error of DDR can also be calculated (Eq.(39)). In Figure 7, the variances of the measured signals are obviously greater than the variances of the reconciled signals. In fact, the covariance of the process outputs for MMVC based on the measured signals is

2.1088 0.0705 ΣMMVC    0.0705 2.0780 28   

ACS Paragon Plus Environment

(50)

Page 29 of 57

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

The MMVC benchmark based on the measured signals can be calculated as

J MMVC  trace  Σ MMVC   4.1868 . The covariance of the process outputs for MMVC based on the reconciled signals after simulation is

1.8962 0.0332 ΣDDRMMVC    0.0332 1.7602

(51)

And the MMVC benchmark based on the reconciled signals can be calculated as

J DDR  MMVC  trace  Σ DDR  MMVC   3.6564 . If the measured signals are used as feedback signals, MMVC will result in larger variances of the process outputs. However, since the dynamic data reconciliation reduces the variances of the measured signals obviously; MMVC based on the reconciled signals will result in smaller variance of the process outputs.

(Figure 7 goes here.)

Condition 2: Measurements with measurement noise and gross errors The outliers are added to the first measurements (process output 1) at time step 100 and 300, and the second measurements (process output 2) at time step 600 and 800. The magnitude of outliers is equal to 3.0. The sign of outliers in the process output 1 is positive, and the sign of outliers in the process output 2 is negative. The biases are introduced into the first measurements from time step 100 to 200, and the second measurements from time step 300 to 400. The magnitude of biases is equal to 2.5. Slow linear drift happens at the same gradient of 0.05 from time step 500 to 600 for 29   

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

Page 30 of 57

the first measurement and from time step 600 to 700 for the second measurement. Based on the proposed DDR-MMVC, the results of the reconciled signals of the process outputs are shown in Figure 8. The variances of measured signals are obviously greater than the variances of the reconciled signals. From the simulation results, the covariance of the process outputs for MMVC based on measured signals is

(Figure 8 goes here.)

2.2983 0.0246 ΣMMVC    0.0246 4.0533

(52)

And the calculated MMVC benchmark based on the measured signals is

J MMVC  trace  Σ MMVC   6.3515 . After the simulation, the covariance of the process outputs for MMVC based on the reconciled signals is

1.9090 0.1388  ΣDDRMMVC    0.1388 2.5633

(53)

And the calculated MMVC benchmark based on the reconciled signals is

J DDR  MMVC  trace  Σ DDR  MMVC   4.4723 . Comparing the above results with those in the first condition, the variances of the process outputs with MMVC based on the measured signals increase obviously, from 2.1088 (Condition 1) to 2.2983 (Condition 2) for the process output 1, and from 2.0780 (Condition 1) to 4.0533 (Condition 2) for the process output 2. The variance of process output 2 especially increases over 95%. With the proposed DDR-MMVC, the variances of the reconciled signals increase only 30   

ACS Paragon Plus Environment

Page 31 of 57

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

a little, from 1.8962 (Condition 1) to 1.9090 (Condition 2) for process output 1, and from 1.7602 (Condition 1) to 2.563 (Condition 2) for process output 2. Moreover, the MMVC benchmarks based on the measured signals increase 51.70%, while the MMVC benchmarks based on the reconciled signals increase only 22.31%. The results show that the application of DDR can decrease the variances of the measured signals obviously. DDR is robust to the influence of the measurement instrument bias, outliers, and slow drifts. When the MMVC uses the reconciled signals as feedback signals instead of the measured signals, the variances of the outputs are decreased. DDR-MMVC can efficiently enhance the control performance.

Condition 3: Comparison results with Monte Carlo simulations Like Case study 1, in this condition, Monte Carlo simulations are also used to demonstrate the effectiveness of DDR in Case study 2 with the multivariable system. The exponential filter is also used for comparisons. The measurements containing only random errors (measurement noise) are considered firstly, and then measurements with both random errors and gross errors are considered as the second case of measurement. The simulation is run 500 times repetitively in each case of measurement. The expected standard deviations of measured signals, filtered signals and reconciled signals are provided in Table 2.

(Table 2 goes here.)

As can be seen in Table 2, the expected standard deviations of the reconciled signal are decreased more significantly compared with the results of the univariate system with Monte Carlo simulations. Considering measurement noise only, the DDR can 31   

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

decrease the expected standard deviation of measured signal by 37.29% (from 1.7432 to 1.0931) in process output 1 and by 38.94% (from 1.7652 to 1.0778) in process output 2, while the DDR can decrease the expected standard deviation only by 16.08% in the univariate system. Considering the system with measurement noise and gross error, the DDR can decrease the expected standard deviation of measured signal by 37.05% (from 1.9263 to 1.2125) in process output 1 and by 38.90% (from 2.1149 to 1.2922) in process output 2, while the DDR can decrease the expected standard deviation only by 17.80% in the univariate system. The results suggest that more efficient error reduction can be derived in the multivariate system because more redundant information of the process model is used by the DDR. The results also indicate that the influence of measurement noise and gross errors can be reduced much more efficiently by DDR comparing with the exponential filter.

Case study 3: Experiments on air-pressure tank systems In the previous cases, the derivations for DDR-MVC or DDR-MMVC are based on the full ARIMA models of the disturbance and the measurement noise. The explicit models are usually unavailable in advance in many real processes. In this experiment, a pilot scale of a quadruple tank apparatus, shown in Figure 9, is used to illustrate the data driven based modeling strategy with all information retrieved from routine closed loop measurement data to enhance the current control performance.

(Figure 9 goes here.)

The MIMO experimental system shown in Figure 9 consists of four interconnected air tanks arranged in two parallel trains of two tanks in series built upon a steel framework. Supply air flows into the system through two air-actuated control valves 32   

ACS Paragon Plus Environment

Page 32 of 57

Page 33 of 57

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

which serve as the manipulated variables for the system. The air flows through tanks before exiting to the atmosphere. Specifically, the air flowing through control valve 1 (Cv1) proceeds into tank 1 and subsequently into tank 2 downstream before exiting the system. Additionally, a portion of the flow from the control valve can be routed into the downstream tank of the adjacent train (tank 4). Similarly, control valve 2 (Cv2) affects the pressure in tank 3 and 4, with cross-flow effect on tank 2. Valves V14 and V23 are directly responsible for the cross-train flow. In the system, the measured outputs are the pressures in tank 2 and tank 4, P2 and P4 (voltages). This arrangement makes both pressures in the two bottom tanks become a function of the two control valves. The tanks are equipped with differential-pressure-to-current (DP/I) transducers to provide a continuous measurement of the pressure. The computer is connected to a AD-622 analog/digital I/O expansion card from Advantech. The expansion board uses a 14-bit converter; therefore, the digital signals are 14-bit. The analog signals from the measured levels are amplified and conditioned using EDM35 (4-20mA/0-5 volts) modules. The visual programming control interface and the controller algorithms (MVC and PI) are developed using MATLAB for this apparatus.

From the derivations in Section 4.2, the reconciled signals ( y   t  ) rely on covariance matrices (  and  ) of the measurement noise and model prediction error (Eq.(38)). The variance or the covariance matrices needs to be estimated based on routine closed loop process measurement data. The covariance matrix  of the vector of measurement noise can be provided by sensor manufactures. The data based estimation for the covariance matrix (  ) from the autoregressive model is given in Appendix B. To get the reconciled signals ( y   t  ) in Eq.(38), the prediction outputs 33   

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

( yˆ (t ) ) should be calculated. In order to estimate the process model ( A( z 1 ) ) and the disturbance model ( B( z 1 ) ), the persistent excitation condition should be imposed on the test signals to guarantee that the estimation yields unique solutions. The pseudo random binary sequence (PRBS) test signals are used to derive a total of 1200 sets of data collected from the closed loop operation; then the conventional identification algorithms can be directly applied to identifying the process model and the disturbance model.

Now the operating performance of the controlled system using the proposed method is evaluated. Two different kinds of the control strategies (MVC and PI) are separately applied. Figures 10 and 11 show the two outputs of MVC and of PI respectively. In Figure 10, the variances of the process output with MVC based on the measurement signals and the reconciled signals in tank 2 are 7.6210-5 and 3.8610-5. The variances in tank 4 are 1.96110-4 and 1.21010-4. In Figure 11, the variances of the process output with PI based on the measurement signals and the reconciled signals in tank 2 are 1.76910-4 and 8.3710-5. The variances in tank 4 are 3.07110-4 and 1.79910-4. This result demonstrates that the variances of the process outputs decrease because the application of a DDR to the process measurement would allow the controllers to access more accurate and precise data. The corresponding control performance would be improved. (Figure 10 goes here.) (Figure 11 goes here.)

34   

ACS Paragon Plus Environment

Page 34 of 57

Page 35 of 57

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

6. CONCLUSIONS From the research reports, about 66% to 80% of the controllers of industrial processes are not performing as well as they are expected. This detrimental effect on the process profitability causes the product variance.  MVC plays a crucial role as a benchmark in performance monitoring. However, the influence of measurement errors is not explicitly considered in the theory for MVC of stochastic control systems. In fact, measurements such as flow rates, pressures, and temperatures from a chemical process are inherently inaccurate. They are contaminated by random errors and possibly gross errors. With additional consideration of measurement errors in the sensor devices, the variance of the output under MVC is obviously increased. Such influence of measurement errors degrades the control performance of MVC in both univariate and multivariate systems. A model-based approach, DDR, is considered to decrease measurement noise inside control loops. The derivations for MVC/MMVC combined with DDR demonstrate that the reconciled signal/signals by DDR is/are more precisely estimated than the measured signal/signals. The variances of outputs by DDR-MVC/DDR-MMVC will be decreased and the control performance can then be enhanced. The effectiveness of DDR-MVC/DDR-MMVC is illustrated by simulation cases and an actual pilot scale experiment, including both univariate and multivariate systems. Results show that both MVC and MMVC benchmarks can be improved

by

reconciled

signal/signals

as

feedback

signal/signals.

DDR-MVC/DDR-MMVC can efficiently enhance the control performance.

The conventional MVC approach would not be applied when the disturbance distribution is non-Gaussian. Entropy is a general measure of uncertainty. It has more physical relevance than the variance of the arbitrary random variables, so a more 35   

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

general approach by incorporating the entropy for the feedback control loop subjected to the unknown disturbance distribution is worth investigation in the future. In some practice applications when MVC is implemented, there are some problems, such as high gain, wide bandwidth and unrealistically large control signal variation. The actuators would result in saturation and excessive wear and tear. In the future, the extension work will be focused on DDR combined with minimal achievable variance of measured output, Such as generalized minimum variance control and linear quadratic Gaussian (LQG) control, which consider both control signals and output signals. Also, DDR is based on the information of operating data and process model, the process model need to be identified firstly based on the previous operating data. DDR combined with a data-driven based subspace model under LQG control to establish the benchmark will be further developed in the future.

The aim of using DDR in the procedure of MVC is to decrease the influence of measurement errors and to enhance the control performance. It can derive accurate approximate results of the actual outputs (y). However, the control structure considered in this paper is known as one degree of freedom control structure. This structure has the disadvantage as the feedback loop properties cannot be designed independently of the reference tracking transfer-function30. The generalized control structure can represent a wide class of controllers, including the one and two degrees of freedom configurations, as well as feedforward and estimation schemes and many others. Besides, it can also be used to formulate optimization problems for controller design and is widely used in the robust control theory31,32. Considering the disturbance, measurement noise or even model uncertainty and uncertain time delay, the use of the generalized control structure combined with the robust controller can 36   

ACS Paragon Plus Environment

Page 36 of 57

Page 37 of 57

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

provide more accurate results. The use of generalized control structure would be a good direction for our further research.



APPENDIX A. THE DETAILED DERIVATIONS OF THE COVARIANCE MATRIX OF THE RECONCILED ERRORS

To use a non-linear optimizer to find the values of the reconciled output yr (t ) for the univariate system and y r (t ) for the multivariate system), just minimize the posterior distribution (of Eq.(25) for the univariate system and of Eq.(37) for the multivariate system). Both derivations are similar, so only the posterior distribution for

the

multivariate

system

is

discussed.

The

partial

derivatives

of

p ( y r (t ) y m (t ), yˆ (t )) (Eq.(36)) with respect to y r (t ) can be expressed analytically

using the equation

p(y r (t ) y m (t ), yˆ (t ))  Σε 1 (y r  y m )  Σδ 1 (y r  yˆ ) y r (t )

(54)

 

By setting the derivative of p ( y r (t ) y m (t ), yˆ (t )) with respect to y r (t ) to zero, after properly reorganizing the equation, the value of y r (t ) can be obtained,

y r  yˆ  K (y m  yˆ )

(55)

where K  ( Σ ε 1  Σ δ 1 ) 1 Σ ε 1 . Thus the estimation error of the reconciled output ( γ ) is

37   

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

γ  yr  y  δ  K (ε  δ)

Page 38 of 57

(56)

The covariance matrix of γ is

Vr  Cov( γ )  Cov δ  K (ε  δ)   Cov  Kε  (I  K )δ   K T Σ ε K  ( I  K )T Σ δ ( I  K )

(57)

where K T Σ ε K and (I  K )T Σ δ (I  K ) can be rewritten as

K T Σ ε K  K T Σ ε ( Σ ε 1  Σ δ 1 ) 1 Σ ε 1  K T ( Σ ε 1  Σ δ 1 ) 1

(I  K )T Σ δ (I  K )  K T ( Σ ε Σ δ 1 )T Σ δ Σ ε Σ δ 1K  K T ( Σ ε Σ δ  1 ) T ( Σ ε 1  Σ δ  1 ) 1

(58)

(59)

where the measurement noise are independent from each other; hence the covariance matrix of the measurement noise ( Σ ε ) is diagonal.

After Eqs.(58) and (59) are substituted into Eq.(57), the covariance matrix of γ can be obtained

Vr  K T ( Σ ε 1  Σδ 1 ) 1  K T ( Σε Σ δ 1 )T ( Σε 1  Σ δ 1 ) 1  [K T  K T ( Σε Σ δ 1 )T ]( Σ ε 1  Σ δ 1 ) 1  K T [I  ( Σε Σ δ 1 )T ]( Σ ε 1  Σ δ 1 ) 1  K T K  T ( Σ ε  1  Σ δ  1 )  1  ( Σ ε  1  Σ δ  1 ) 1

38   

ACS Paragon Plus Environment

(60)

Page 39 of 57

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



APPENDIX B. DATA BASED ESTIMATION OF THE COVARIANCE MATRIX OF THE MODEL PREDICTION ERRORS.

Data based estimation of the covariance matrix (  δ for the multivariate system) of the model prediction errors would be derived as below. Because the derivations for the estimation of the variance  2 for the univariate system and of the covariance  δ for the multivariate system are similar, so only the derivations for the multivariate system are discussed.

The multivariate system described (Eq.(7)) is operated by feedback-only controllers in the closed loop, and the controller is expressed in a transfer function form as

u (t )   C ( z  1 ) y m (t )

(61)

After substituting Eq.(61) into Eq.(7), the corresponding measured outputs of the system for the disturbance changes and the measurement noise are

y m (t )   A ( z 1 )C( z 1 ) y m (t )  B ( z 1 )d (t )  ε (t )

(62)

 

where p(t )  B( z 1 )d(t )  ε(t ) is defined as the unmeasured part. Thus, the first term on the right side of Eq.(62) can be estimated by a finite-length autoregressive model. Using the prediction error pˆ (t ) , the estimated covariance matrix of p(t ) can be calculated by

Σ p  cov pˆ (t )   cov[B ( z 1 )d (t )  ε (t )] 39   

ACS Paragon Plus Environment

(63)

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

Page 40 of 57

The disturbances to be considered in this paper are caused by the process equipment and only act on the process, the measurement noise is the variation in the sensor reading that does not correspond to changes in the process. Therefore, the disturbances and the measurement noise are assumed to be independent. the covariance matrix of p(t ) is the linear combination of the covariance matrix of the model prediction error ( Σδ ) and the covariance matrix of the measurement noise ( Σε ),

Σp  cov[B( z 1 )d(t )]  cov[ε(t )]  Σδ  Σε

(64)

Thus, the estimated covariance matrix of the model prediction error can be calculated by

Σδ  Σp  Σε

(65)

 



AUTHOR INFORMATION Corresponding Author *Tel.: +886-3-2654107. Fax: +886-3-2654199. E-mail: [email protected]. Notes The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors gratefully acknowledge Ministry of Science and Technology, R.O.C. (MOST 103-2221-E-033-068-MY3); National Natural Science Foundation of 40 

 

ACS Paragon Plus Environment

Page 41 of 57

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

China (No. 61374167); Natural Science Foundation of Zhejiang Province (No. LQ14F030006); and Science and Technology Planning Project of Zhejiang Province (No. 2015C31157; 2014C31074; 2014C31093) for financial support.



REFERENCES

(1). Desborough, L.; Miller, R. Increasing customer value of industrial control performance monitoring-Honeywell's experience. In AIChE symposium series (pp. 169-189). New York; American Institute of Chemical Engineers; 1998. (2). Aström, K. J. Computer control of a paper machine-an application of linear stochastic control theory. IBM Journal 1967, 11, 389-05. (3). DeVries, W. R.; Wu, S. M. Evaluation of process control effectiveness and diagnosis of variation in paper basis weight via multivariate time-series analysis. IEEE Transactions on Automatic Control 1978, 23, 702-708. (4). Harris, T. J. Assessment of control loop performance. The Canadian Journal of Chemical Engineering 1989, 67, 856-861. (5). Uduehi, D.; Ordys, A.; Grimble, M.; Majecki, P.; Xia, H. Controller benchmarking procedures—data-driven methods. In Process Control Performance Assessment (pp. 81-126). Springer London. 2007. (6). Li, Z.; Evans, R. J. Minimum-variance control of linear time-varying systems. Automatica 1997, 33, 1531-1537. (7). Harris, T. J.; Boudreau, F.; MacGregor, J. F. Performance assessment of multivariable feedback controllers. Automatica 1996, 32, 1505-1518. (8). Huang, B.; Shah, S. L.; Fujii, H. The unitary interactor matrix and its estimation using closed-loop data. Journal of Process Control 1997, 7, 195-207. (9). Yan, Z.; Chan, C. L.; Yao, Y. Multivariate control performance assessment and 41   

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

control

system

monitoring

via

hypothesis

test

on

Page 42 of 57

output

covariance

matrices. Industrial & Engineering Chemistry Research 2015, 54, 5261−5272. (10).

Ko, B. S.; Edgar, T. F. Performance assessment of multivariable feedback

control systems. Automatica 2001, 37, 899-905. (11).

Ko, B. S.; Edgar, T. F. PID control performance assessment: The single‐

loop case. AIChE Journal 2004, 50, 1211-1218. (12).

Huang, B. Minimum variance control and performance assessment of

time-variant processes. Journal of Process Control 2002, 12, 707-719. (13).

Olaleye, F. B.; Huang, B.; Tamayo, E. Feedforward and feedback controller

performance assessment of linear time-variant processes. Industrial & Engineering Chemistry Research 2004, 43, 589-596. (14).

Patwardhan, R. S.; Shah, S. L.; Qi, K. Z. Assessing the performance of model

predictive controllers. The Canadian Journal of Chemical Engineering 2002, 80, 954-966. (15).

Harrison, C. A.; Joe Qin, S. Minimum variance performance map for

constrained model predictive control. Journal of Process Control 2009, 19, 1199-1204. (16).

Chen, J.; Wang, W. Y. PCA− ARMA-Based Control Charts for Performance

Monitoring of Multivariable Feedback Control. Industrial & Engineering Chemistry Research 2010, 49, 2228-2241. (17).

Sendjaja, A. Y.; Kariwala, V. Minimum variance benchmark for performance

assessment of decentralized controllers. Industrial & Engineering Chemistry Research 2012, 51, 4288-4298. (18).

Kulcsár, C.; Raynaud, H. F.; Petit, C.; Conan, J. M. Minimum variance

prediction and control for adaptive optics. Automatica 2012, 48, 1939-1954. 42   

ACS Paragon Plus Environment

Page 43 of 57

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

(19).

Farasat, E.; Huang, B. Deterministic vs. stochastic performance assessment

of iterative learning control for batch processes. AIChE Journal 2013, 59, 457-464. (20).

Alipouri, Y.; Poshtan, J. A linear approach to generalized minimum variance

controller design for MIMO nonlinear systems. Nonlinear Dynamics 2014, 77, 935-949. (21).

Grimble, M. J.; Majecki, P. Non-linear generalised minimum variance

control using unstable state-dependent multivariable models. IET Control Theory & Applications 2013, 7, 551 - 564. (22).

Shardt, Y.; Zhao, Y.; Qi, F.; Lee, K.; Yu, X.; Huang, B.; Shah S.

Determining the state of a process control system: current trends and future challenges. Canadian Journal of Chemical Engineering 2012, 90, 217–245. (23).

Hugo,

A.

J.

Performance

assessment

of

single-loop

industrial

controllers. Journal of Process Control 2006, 16, 785-794. (24).

Huang, B. and Shah, S. L. Performance Assessment of Control Loops.

Springer-Verlag, London. 1999. (25).

Bai, S.; Thibault, J.; McLean, D. D. Dynamic data reconciliation: Alternative

to Kalman filter. Journal of Process Control 2006, 16, 485-498. (26).

Jiang, X.; Liu, P.; Li, Z. Data reconciliation for steam turbine on-line

performance monitoring. Applied Thermal Engineering 2014, 70, 122-130. (27).

Narasimhan S.; Jordache C. Data reconciliation and gross error detection.

In Measurement Errors and Error Reduction Techniques (pp. 32-58). Gulf Publishing Company, Houston. 2000. (28).

Huang, B.; Shah, S. L.; Kwok, E. K. Good, bad or optimal? Performance

assessment of multivariable processes. Automatica 1997, 33, 1175-1183. 43   

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

(29).

Zhao Y.; Su H.; Chu J.; Gu Y. Multivariable Control Performance

Assessment Based on Generalized Minimum Variance Benchmark, Chinese Journal of Chemical Engineering 2010, 18, 86-94. (30).

Grimble, M. J. H2/ h∞ robust control design for a generalized control

structure and flight control application. International Journal of Systems Science 1995, 26, 2043-2068. (31).

Skogestad, S.; Postlethwaite, I. Multivariable Feedback Control: Analysis

and Design, John Wiley & Sons, Chichester. 1996. (32).

Zhou, K.; Doyle, J. C.; Glover, K. Essentials of Robust Control, Prentice hall,

Upper Saddle River, 1998.

44   

ACS Paragon Plus Environment

Page 44 of 57

Page 45 of 57

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

TOC graphic



45   

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

Table Caption Table 1 Comparison results for the univariate system with 500-repetition simulations. Table 2 Comparison results for the multivariate system with 500-repetition simulations

ACS Paragon Plus Environment

Page 46 of 57

Page 47 of 57

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

Table 1 Comparison results for the univariate system with 500-repetition simulations Expected standard deviation

Conditions of measurement

Measured signal

Filtered signal

Reconciled signal

  0.025

0.0342

0.0301

0.0287

Random errors with   0.025 and gross errors

0.0354

0.0308

0.0291

Random errors with at

Table 2 Comparison results for the multivariate system with 500-repetition simulations Conditions of measurements Random errors with at Σ ε  I Random errors with Σ ε  I and gross errors

Expected standard deviation Measured signals Filtered signals Reconciled signals output 1 output 2 output 1 output 2 output 1 output 2 1.7432 1.7652 1.2582 1.3299 1.0931 1.0778 1.9263

2.1149

1.4750

ACS Paragon Plus Environment

1.7526

1.2125

1.2922

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

Figure Captions Figure 1 The stochastic univariate (or multivariate) feedback control system. Figure 2 The stochastic univariate (or multivariate) system under the consideration of measurement noise. Figure 3 The stochastic univariate (or multivariate) feedback control system with DDR. Figure 4 Results of MVC for the univariate system based on the measured signal and the reconciled signal when there are measurement noise with   0.025 in the univariate system. Figure 5 Results of MVC for the univariate system based on the measured signal and the reconciled signal when there is measurement noise with   0.05 in the univariate system. Figure 6 Results of MVC for the univariate system based on the measured signal and the reconciled signal when there are measurement noise with   0.05 and the gross errors in the univariate system. Figure 7 Measured signals and reconciled signals of the process outputs when there are two measurement noises with Σε  I in the multivariate system. Figure 8 Measured signals and reconciled signals of the process outputs when there are two measurement noises with Σε  I and the gross errors in the multivariate system. Figure 9 Schematic of the four-tank air pressure system in Case Study 3. Figure 10 Process outputs with MVC based on the measurement signals and the reconciled signals in tank 2 and tank 4 in the air-pressure tank experiments. Figure 11 Process outputs with PI controllers based on the measurement signals and the reconciled signals in tank 2 and tank 4 in the air-pressure tank experiments.

ACS Paragon Plus Environment

Page 48 of 57

Page 49 of 57

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

d (t )/d(t )

D ( z 1 ) A( z 1 ) B ( z 1 )

r (t )

r (t )



e(t )

C ( z 1 )

u (t )

e(t )

C( z 1 )

u(t )

z k

B ( z 1 ) A( z 1 ) A ( z 1 )

y (t ) y (t )

  Figure 1 The stochastic univariate (or multivariate) feedback control system.

d (t )/d(t )

D ( z 1 ) A( z 1 ) B ( z 1 )

r (t ) r (t )



e(t )

C ( z 1 )

u (t )

e(t )

C( z 1 )

u(t )

z k

B ( z 1 ) A( z 1 ) A ( z 1 )

y (t ) y (t )

ym (t ) y m (t )

 (t ) / ε(t )

  Figure 2 The stochastic univariate (or multivariate) system under the consideration of measurement noise.

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

d (t )/d(t ) D ( z 1 ) A( z 1 ) B ( z 1 )

r (t ) r (t )



e(t )

C ( z 1 )

u (t )

e(t )

C( z 1 )

u(t )

z k

B ( z 1 ) A( z 1 ) A ( z 1 )

yr (t )

ym (t )

y r (t )

y m (t )

y (t )

y (t )

 (t ) / ε(t )

Figure 3 The stochastic univariate (or multivariate) feedback control system with DDR.

0.08

MVC based on measured signal MVC based on reconciled signal

0.06 0.04 Process output

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 50 of 57

0.02 0 -0.02 -0.04 -0.06 -0.08

0

100

200

300

400

500 600 Time step

700

800

900

1000

Figure 4 Results of MVC for the univariate system based on the measured signal and the reconciled signal when there are measurement noise with   0.025 in the univariate system.

ACS Paragon Plus Environment

Page 51 of 57

0.15

MVC based on measured signal MVC based on reconciled signal

0.1

0.05 Process output

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

0

-0.05

-0.1

-0.15

-0.2

0

100

200

300

400

500 600 Time step

700

800

900

1000

Figure 5 Results of MVC for the univariate system based on the measured signal and the reconciled signal when there is measurement noise with   0.05 in the univariate system.

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

0.1

MVC based on measured signal MVC based on reconciled signal

0.05

Process output

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 52 of 57

0

-0.05

-0.1

-0.15

0

100

200

300

400

500 600 Time step

700

800

900

1000

Figure 6 Results of MVC for the univariate system based on the measured signal and the reconciled signal when there are measurement noise with   0.05 and the gross errors in the univariate system.

ACS Paragon Plus Environment

Page 53 of 57

Process output 1

10

measured reconciled

5

0

-5

0

100

200

300

400

500 600 Time step

700

800

10 Process output 2

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

900

1000

measured reconcilied

5

0

-5

0

100

200

300

400

500 600 Time step

700

800

900

1000

Figure 7 Measured signals and reconciled signals of the process outputs when there are two measurement noises with Σε  I in the multivariate system.

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

Process output1

10

measured reconciled

5 0 -5 -10

0

100

200

300

400

500 600 Time step

700

800

10 Process output2

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 54 of 57

900

1000

measured reconciled

5 0 -5 -10

0

100

200

300

400

500 600 Time step

700

800

900

1000

Figure 8 Measured signals and reconciled signals of the process outputs when there are two measurement noises with Σε  I and the gross errors in the multivariate system.

ACS Paragon Plus Environment

Page 55 of 57

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

Cv2

Cv1 V14

Tank 1

V23

Tank 3 Controller

Controller

Tank 2 P2

Tank 4 PT P4

PT

Figure 9 Schematic of the four-tank air pressure system in Case Study 3.

ACS Paragon Plus Environment

Page 56 of 57

measured reconciled

2.84

2.82

2.8 600

Process output in tank 4

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

Process output in tank 2

Industrial & Engineering Chemistry Research

700

800

900 Time step

1000

1100

1200

measured reconciled

2.8

2.78

2.76 600

700

800

900 Time step

1000

1100

1200

Figure 10 Process outputs with MVC based on the measurement signals and the reconciled signals in tank 2 and tank 4 in the air-pressure tank experiments.

ACS Paragon Plus Environment

measured reconciled

2.84 2.82 2.8 2.78 600

Process output in tank 4

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

Process output in tank 2

Page 57 of 57

700

800

900 Time step

1000

1100

1200

measured reconciled

2.82 2.8 2.78 2.76 600

700

800

900 Time step

1000

1100

1200

Figure 11 Process outputs with PI controllers based on the measurement signals and the reconciled signals in tank 2 and tank 4 in the air-pressure tank experiments.

ACS Paragon Plus Environment