Next Article in Journal
A Hybrid Inverse Problem Approach to Model-Based Fault Diagnosis of a Distillation Column
Next Article in Special Issue
A Hybrid Framework for Simultaneous Process and Solvent Optimization of Continuous Anti-Solvent Crystallization with Distillation for Solvent Recycling
Previous Article in Journal
Thermal Isolation of a Clean Alloy from Waste Slag and Polymeric Residue of Electronic Waste
Previous Article in Special Issue
Modeling and Simulation of Reaction and Fractionation Systems for the Industrial Residue Hydrotreating Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Triple-Mode Model Predictive Control Using Future Target Information

National Center for International Research on Quality-Targeted Process Optimization and Control, College of Control Science and Engineering, Zhejiang University, Hangzhou 310027, China
*
Author to whom correspondence should be addressed.
Submission received: 9 December 2019 / Revised: 27 December 2019 / Accepted: 30 December 2019 / Published: 2 January 2020
(This article belongs to the Special Issue Process Optimization and Control)

Abstract

:
In this paper, we propose a triple-mode model predictive control (MPC) algorithm that uses future target information to improve tracking performance. To explicitly take into account the future target information in the MPC optimization, the proposed triple-mode control law encompasses three parts: (i) the future target information feedforward, (ii) the output feedback, and (iii) the extra degrees of freedom for constraint satisfaction. The first two parts of the control law are off-line designed through unconstrained MPC, and the optimal future trajectory horizon is obtained by golden section search based on the integral of squared error (ISE) criterion. The final part is calculated by the on-line MPC algorithm aiming to satisfy constraints. Furthermore, we analyze the feasibility and convergence properties of the proposed algorithm. The method is demonstrated by the simulation of the shell fundamental control problem and also tested on the coordinated control problem in the power plant. The test results show that the proposed algorithm can increase tracking performance dramatically due to the proper selection of future trajectory horizon.

1. Introduction

Concern about multivariable constrained system control has become a central issue in the control domain. As a result, much research in recent years has focused on the development of advanced control algorithms that are process friendly. Model predictive control (MPC) has been widely and successfully applied in the process industry, primarily due to its superior performance and the capability of dealing with a multivariate constrained problem [1,2,3]. Since the step-response model is more convenient to obtain and more suitable for the process industry, this model has been successfully applied to numerous industrial processes by Aspen Technology, Honeywell Hi-Spec, and others. There has been an increase in the adoption of dynamic matrix control in MPC commercial software packages such as DMC+, SMC, RMPCT, HIECON [1].
The key idea of model predictive control is to use the historical information and model prediction to calculate the optimal input trajectory by minimizing the deviation between the process output and the desired target in a receding horizon way. There exists a rich supporting theory [4,5,6] for its success in industry. Moreover, the applications of model predictive control have already gone beyond petrochemical or chemical industries and extended to power plants [7], traffic, robots [8] and any other fields.
The research to date has tended to focus on the regulation problem (steer the system to the origin) [4] for facilitating the analysis and design. The regulation problem assumes the hypothesis that the future setpoint is fixed or changes instantaneously, thus future target information is rarely discussed in the current literature. Considering the setpoint changes, Limón et al. [9] adopted a modified cost function and a stabilizing extended terminal constraint in MPC optimization to obtain the feasible setpoints in case that the setpoints are not feasible. Their work is still preserving no consideration of using the future target information. The MPC optimization have the capability of considering the future target information [10], which may be helpful in control since it allows the controller to react to the setpoint changes in advance [11]. Carrasco and Goodwin [12] proposed a feedforward model predictive control, in which the feedforward was designed individually before the feedback design. However, how much future target information (so-called “preview” in their work) should be used in feedforward design is not discussed in their work. Rossiter and Grinnell [13] developed an extended input horizon generalized predictive controller (EIHGPC) algorithm, in which they parameterized the degrees of freedom to spread the control increments over the prediction horizon and adopted a heuristic algorithm to select a proper future target horizon, hence improves the tracking performance. Their method only pertains to the unconstrained case. In [14], the authors proposed a dual-mode constrained MPC algorithm, which was composed of a linear state feedback and extra degrees of freedom. Under such a control law, the feedforward could not be designed individually and the future target horizon was determined through a trial and error process.
However, the inappropriate use of future target information may cause the deterioration of control performance [13,15]. There remains a need for an effective method that can determine a proper future target horizon. Moreover, the current implementations of MPC fail to separately design the feedforward, feedback and constraint satisfaction, which could not yield the desirable performance from each part. In this study, a triple-mode constrained model predictive control algorithm is developed, which has the merit of independent design for each mode of the control law. The main ingredients of the triple-mode control law are: (i) the feedforward design using future target information, (ii) the output feedback, and (iii) the extra degrees of freedom for constraint satisfaction. In addition, a golden section search method is proposed to optimize the future trajectory horizon, which can be used to design the feedforward to improve the tracking performance.
The outline of this article is as follows. First, we describe the typical dynamic matrix control. Next, we give the future target information feedforward structure derived from the unconstrained control law and propose an optimization method for the selection of future trajectory horizon. We then introduce the triple-mode model predictive control using future target information. Finally, we provide two illustrative simulation examples.

2. Preliminary

Mathematical notations used throughout this paper are defined as follows. Given a vector x R n , the symbol | | x | | denotes the Euclidean 2-norm and | | x | | P 2 denotes the quadratic form of x. e.g., | | x | | P 2 = x T P x ; A T denotes the transpose of A; Given two integers, l m , we define the set I l : m = { l , l + 1 , , m 1 , m } . The superscript ∗ denotes the optimal solution or cost according to the context.

Dynamic Matrix Control Algorithm

The design of the MPC algorithm in this paper follows the standard DMC approach proposed by Garcia and Morshedi [2]. Without loss of generality, we consider a MIMO linear dynamic system with n y R controlled variables and n u R manipulated variables. The DMC algorithm adopts a non-parametric step response model as the predictor, thus the ith output y i ( t ) R at time t can be formulated by the following step response model
y i ( t ) = j = 1 n u k = 1 a i , j , k Δ u j ( t k ) , ( 1 j n y )
in which Δ u j ( t ) R denotes the jth input increment, i.e., Δ u j ( t ) u j ( t ) u j ( t 1 ) . a i , j , k R denotes the kth element of the corresponding step response, and the step response sequence of the ith output to the jth input is { a i , j , 1 , a i , j , 2 , , a i , j , n , } . For a stable plant this sequence will asymptotically reach a constant value, i.e., a i , j , n + m a i , j , n , ( m > 0 ).
Based on the step response model, the output predictions over n future steps can be expressed as follows
Y ( t + 1 ) = Λ Y ( t ) + A Δ u ( t )
where
Y ( t ) = [ y T ( t | t 1 ) , y T ( t + 1 | t 1 ) , , y T ( t + n 1 | t 1 ) ] T y ( t + k | t 1 ) = [ y 1 ( t + k | t 1 ) , y 2 ( t + k | t 1 ) , , y n y ( t + k | t 1 ) ] T , k I 0 : n 1 Δ u ( t ) = [ Δ u 1 ( t ) , Δ u 2 ( t ) , , Δ u n u ( t ) ] T Λ = 0 I n y 0 0 0 0 0 I n y 0 0 0 0 0 I n y 0 0 0 0 0 I n y 0 0 0 0 I n y , A = A 1 A 2 A n A t = a 1 , 1 , t a 1 , 2 , t a 1 , n u , t a 2 , 1 , t a 2 , 2 , t a 2 , n u , t a n y , 1 , t a n y , 2 , t a n y , n u , t , t I 1 : n
Δ u ( t ) R n u is a vector of the future control moves at time t. Y ( t ) R n y · n is the dynamic states of the system, and y ( t + k | t 1 ) R n y is the future output vector at time t + k assuming the input remain constant starting at time t 1 , which means Δ u ( t + j ) = 0 for j 0 . The evolution of predictions can be obtained by translation of one-step prediction plus the contribution made by the latest control move Δ u ( t ) .
To achieve the goal of offset-free tracking, a bias is introduced into the model to incorporate the feedback. It’s assumed that the model prediction and the measured output at the current time is due to a constant step disturbance at the output. The constant output step disturbance assumption is standard in industrial applications [16]. The model bias e ( t ) R n y is based on the latest output measurement y ˜ ( t ) R n y and the predicted outputs y ^ ( t | t 1 ) R n y
e ( t ) = y ˜ ( t ) y ^ ( t | t 1 )
in which y ^ ( t | t 1 ) refers to the prediction of output at time t based on the measurements up to time t 1 . The bias is added to compensate the output prediction
Y ^ ( t | t ) = Y ^ ( t | t 1 ) + H e ( t )
where Y ^ ( t | t 1 ) denotes the estimate of Y ( t ) at time t before the latest measurements are obtained. Y ^ ( t | t ) is the revised estimate of Y ( t ) on the basis of the feedback, and H : = d i a g { I n y , , I n y } refers to the weighting matrix of the bias term.
For open-loop stable systems, the future control moves Δ U M ( t ) R M · n u over the control horizon (M) is optimized to drive the model predicted outputs Y ^ P ( t + 1 | t ) R P · n y as closely as possible to a desired future trajectory R ( t + 1 | t ) R P · n y over prediction horizon (P). The DMC optimization problem is defined as
min Δ u M ( t ) J ( t ) = R ( t + 1 | t ) Y ^ P ( t + 1 | t ) Q 2 + Δ U M ( t ) W 2
s . t . Y ^ P ( t + 1 | t ) = Λ P Y ^ ( t | t ) + A P M Δ U M ( t )
Δ u m i n Δ U M ( t ) Δ u m a x
u m i n U M ( t ) u m a x
where u m a x , u m i n , Δ u m a x , Δ u m i n are input constraints and increment constraints, respectively. The matrix A P M is the dynamic matrix with P × M block matrix and Λ P is the shift matrix with appropriate dimension. Q , W are the weighting matrices for outputs and the controller moves. The above notations read
Y ^ P ( t + 1 | t ) = y ^ T ( t + 1 | t ) , , y ^ T ( t + P | t ) T , Δ U M ( t ) = Δ u T ( t | t ) , , Δ u T ( t + M 1 | t ) T , R ( t + 1 | t ) = r T ( t + 1 | t ) , , r T ( t + P | t ) T , r T ( t + k | t ) = r 1 ( t + k | t ) , , r n y ( t + k | t ) T , k I 1 : P A P M = A 1 0 0 A 2 A 1 0 A M A M 1 A 1 A P A P 1 A P M + 1 , Λ P = I P · n y × P · n y 0 × Λ , Q = d i a g ( Q 1 , , Q P ) , Q 1 = d i a g { q 1 , , q n y } W = d i a g ( W 1 , , W M ) , W 1 = d i a g { w 1 , , w n u } .
At every sample time, controller evaluates its future input trajectories by solving the optimization problem (5) and sends the first step of its entire future input trajectory to its actuators. However, although the future target information R ( t + 1 | t ) of the controlled variable were embedded in the optimization problem of DMC algorithm, the default choice in DMC algorithm was r ( t + k | t ) = r ( t + 1 | t ) , k I 2 : P . Little attention has been paid to the selection of an appropriate future trajectory horizon and the use of future target information.

3. Main Results

In this section, we first derive the explicit control law of the unconstrained MPC, which includes two parts: (i) the future target information feedforward, and (ii) the output feedback. Based on the future target information feedforward, we demonstrate that a poor future trajectory horizon has quite an impact on the system’s performance. Then we propose an optimization method to determine the future trajectory horizon and the tripe-mode model predictive control using future target information. The triple-mode MPC is composed of the designed unconstrained control law and the extra degrees of freedom. The future target information is taken into account explicitly to improve tracking performance.

3.1. Illustrations of Poor Future Trajectory Horizon

Now we show that the unconstrained MPC control law will provide possibilities in the explicit expression of feedforward and feedback. Consider the performance objective J ( t ) in problem (5), one can derive the unconstrained law by minimizing it w.r.t the future control moves, according the extremum condition d ( J ( t ) ) / d ( Δ U M ( t ) ) = 0 , we have
Δ U M ( t ) = ( ( A P M ) T Q A P M + W ) 1 ( A P M ) T Q ( R ( t + 1 | t ) Λ P Y ^ ( t | t ) )
and making the following substitutions
F r ( ( A P M ) T Q A P M + W ) 1 ( A P M ) T Q K ( ( A P M ) T Q A P M + W ) 1 ( A P M ) T Q Λ P
then we can get
Δ u M ( t ) = F r R ( t + 1 | t ) K Y ^ ( t | t )
It can be seen that the future target information is involved in the first term of the control law: F r R ( t + 1 | t ) , which can be seen as the feedforward in the unconstrained control law. In addition, the second part of the control law can be seen as the output feedback based on the current output prediction.
For convenience, we define the future trajectory horizon in feedforward term for each output as f i R , i I 1 : n y , which indicates the length of future target information considered in the optimization. Any future targets beyond the future trajectory horizon will not be taken into account in the optimization problem, even if known. For a given f i , i I 1 : n y , the future target information adopted in MPC can be expressed as
r i ( t + f i + k | t ) = r i ( t + f i | t ) , k I 1 : P f i
Furthermore, by combining the specific future trajectory horizon, the unconstrained law can be represented as follows:
F r R ( t + 1 | t ) = F 1 r 1 + + F n y r n y = F ˜ 1 r ˜ 1 + + F ˜ n y r ˜ n y
where F i r i , i I 1 : n y refers to the ith feedforward contribution, and F ˜ i r ˜ i , i I 1 : n y is the reconstructed future target information feedforward component based on the designed future trajectory horizon f i .
F i = F i , F n y + i , , F ( P 1 ) · n y + i , F ˜ i = F i , F n y + i , , F ( f i 2 ) · n y + i , F ( f i 1 ) · n y + i + + F ( P 1 ) · n y + i , r i = r i ( t + 1 | t ) r i ( t + 2 | t ) r i ( t + P | t ) , r ˜ i = r i ( t + 1 | t ) r i ( t + 2 | t ) r i ( t + f i | t ) , i I 1 : n y
and F j , j I 1 : P · n y denotes the jth column vector of the feedforward matrix F r .
It indicates that, for the ith output, the original P-dimension target vector r i has been replaced by a designed f i -dimension vector r ˜ i . A major current focus in improving the control performance is how to select a proper f i , i I 1 : n y . It is known that the typical DMC algorithm considers no setpoint changes during the optimizaiton, which is equivalent to set the future trajectory horizon as f i = 1 , i I 1 : n y .
Next, we examine the closed-loop performances of the MPC algorithm adopting different choices of future trajectory horizon. Consider a single-input, single-output system described by the following transfer function
G ( s ) = 1 ( 10 s + 1 ) ( 5 s + 1 )
The sampling time is chosen as 5 s, a step change of setpoint happens at sampling instant 15 and the unconstrained MPC algorithm is applied to drive the system to its target. The arbitrary choices of future trajectory horizon are f = 1 , 6 , 15 . The tuning parameters of DMC are: P = 15 , M = 3 and Q = I , R = 5 I in which I is the identity matrix.
As illustrated by Figure 1, with a default choice of future trajectory horizon in MPC ( f = 1 ), the system response is slugglish when the setpoint changes, in contrast, setting the future trajectory horizon as the prediction horizon could result in an aggressive response. A proper selection f = 6 deserves a desirable tracking performance.
The integral of squared error (ISE) criterion is widely used to evaluate the system performance in control engineering, the discrete form of ISE reads
I S E = t = 0 e r r o r ( t ) 2
in which e r r o r ( t ) refers to the difference between the expected value (target value) and the actual response value at time instant t.
The ISE values shown in Table 1 demonstrate that the standard MPC algorithm with no anticipation ( f = 1 ) owns large ISE value during the transition from one target to another. On the other hand, with large number of future trajectory horizon ( f = 15 ), the system would react far in advance and lead to bad performance too.

3.2. Future Trajectory Horizon Optimization Algorithm

Since the closed-loop response with smaller ISE value is better, we advocate use of the ISE value to select the proper future trajectory horizon. The ISE value can be regarded as a function of future trajectory horizon, i.e., I S E i = g ( f i ) , i I 1 : n y . A better f i deserves lower function value g ( f i ) , we take the idea of golden section search algorithm [17] to optimize the future trajectory horizon. The algorithm can be summarized as follows:
According to Algorithm 1, we plot the ISE value subject to the iteration of a k , b k , λ k , μ k as Figure 2 shows. The optimal future trajectory horzion obtained is f = μ 6 = 6 . Which testifies the result of last section.
Algorithm 1: The optimization of Future trajectory horizon based on ISE value
  • For any controlled variable y i , i I 1 : n y ;
  • Step 1. Given the prediction horizon P, set a 1 = 1 , b 1 = P , then calculate the initial point:
    λ 1 = a 1 + 0.382 ( b 1 a 1 ) μ 1 = a 1 + 0.618 ( b 1 a 1 )
    and set k = 1 ;
  • Step 2. Give a specified target, govern the process by using the control law Equation (8) and evaluate the ISE value g ( · ) over the entire runtime;
  • Step 3. If a k = = λ k , set f i = μ k or if b k = = μ k set f i = λ k , and the optimal future trajectory horizon is obtained, break the optimization; Otherwise, if g ( λ k ) > g ( μ k ) , go to Step 4, else if g ( λ k ) g ( μ k ) , go to Step 5;
  • Step 4. Set
    a k + 1 = λ k b k + 1 = b k λ k + 1 = μ k μ k + 1 = a k + 1 + 0.618 b k + 1 a k + 1
    and go to Step 6;
  • Step 5. Set
    a k + 1 = a k b k + 1 = μ k μ k + 1 = λ k λ k + 1 = a k + 1 + 0.382 b k + 1 a k + 1
    and go to Step 6;
  • Step 6. Set k = k + 1 , go to Step 2.

3.3. Triple-Mode MPC Algorithm Using Future Target Information

For the conventional MPC algorithm, the optimization problem (5) is solved with n u · M degree of freedom Δ U M ( t ) to tackle constraints. However, the original optimization problem treats the input as decision variables, which breaks down the inherent future target information feedforward structure as described before. What’s more, the default choice of future trajectory horizon f i = 1 , i I 1 : n y in the MPC algorithm may lead to sluggish tracking performance. There exists no effective method to determine a proper future target horizon for MPC optimization.
With the proposed Algorithm 1 used, we next present a triple-mode MPC algorithm using future target information. Different from the typical DMC algorithm as described in preliminary section, the proposed triple-mode control law encompasses three parts: (i) the future target information feedforward, (ii) the output feedback, and (iii) the extra degrees of freedom. Each mode of the control law can be designed individually. The first two parts of the control law are off-line designed through unconstrained MPC, and the optimal future trajectory horizon is obtained by using Algorithm 1. The final part is calculated by the on-line MPC algorithm aiming to satisfy constraints.
Enhance the unconstrained MPC control law Equation (8) with an additional term C ( t ) R n u · M
Δ U M ( t ) = F r R ( t + 1 | t ) K Y ^ ( t | t ) + C ( t )
in which C ( t ) = [ c ( t | t ) , c ( t + 1 | t ) , , c ( t + M 1 | t ) ] T and c ( t + i | t ) , i I 0 : M 1 R N u representing degrees of freedom to get the constraints satisfied.
It indicates that the original degree of freedom Δ U M ( t ) in constrained DMC optimization (5) is replaced by C ( t ) , which not only provide the feedforward design from future target information but also guarantee the advantage of constraint satisfaction. At each sampling time, the controller solves the following optimization problem
min C ( t ) C ( t ) W 2
s . t . Y ^ P ( t + 1 | t ) = Λ P Y ^ ( t | t ) + A P M Δ U M ( t )
Δ U M ( t ) = F r R ( t + 1 | t ) K Y ^ ( t | t ) + C ( t )
Δ u m i n Δ U M ( t ) Δ u m a x
u m i n U M ( t ) u m a x
Y ^ P ( t + P | t ) = R ( t + P | t )
where Equation (17f) refers to the terminal constraint used to ensure closed-loop stability when model predictive control is employed. After optimization, the optimal solution C ( t ) = [ c ( t | t ) , c ( t + 1 | t ) , , c ( t + M 1 | t ) ] T can be obtained, and the optimal input increment is nothing but Δ u ( t ) = d T ( F r R ( t + 1 | t ) K Y ^ ( t | t ) ) + c ( t | t ) with d = [ I n u , 0 , , 0 ] T . Note that the proposed triple-mode MPC optimiazation problem (17) is a quadratic programming [18] with n u × M degrees of freedom. The degrees of freedom are the same as the original DMC optimization problem (5), thus the computation load of proposed triple-mode MPC is the same as which of DMC.
A detailed implementation procedure of the proposed triple-mode MPC algorithm (Algorithm 2) is summarized as follows.
Algorithm 2: Triple-mode MPC algorithm using future target information
  • Offline computation:
  • Carry out the Algorithm 1 to obtain the future trajectory horizon f i , i R n y , the first two parts of the control law are designed through Equations (7) and (10);
  • Online computation:
  • Step 1. At sampling time t, controller receives measurement from the sensors and evaluate the modified estimate of predictions according to Equation (4);
  • Step 2. Solve the optimization problem (17) to evaluate its future input trajectory;
  • Step 3. Injects the optimal input into the plant;
  • Step 4. Set t = t + 1 , and go to Step 1.
In this section, the triple-mode MPC using future target information is presented. As has been previously stated, this triple-mode control law is based on the usage of the perturbations on the fixed dual-mode unconstrained control law. To this end, the following assumption is considered.
Assumption A1.
The step response sequence a i , j , k , i I 1 : n y , j I 1 : n u , k N + is a Cauchy sequence: ϵ > 0 , N N , n , m , ( m , n > N | a i , j , n a i , j , m | < ϵ ) .
Assumption A2.
The weighting matrices Q R P · n y × P · n y and R R M · n u × M · n u are positive definite matrices. At least as many inputs as outputs is required: n u n y .
The future target is admissible if there exist a set of control sequence Δ U M ( t ) such that the system output Y ^ P ( t ) can reach it and the constraints (e.g., Δ u m i n Δ U M ( t ) Δ u m a x , u m i n U M ( t ) u m a x ) can be satisfied during the stage. The following theorem follows immediately.
Theorem 1.
Consider that Assumption 1 and 2 hold, and the future target R ( t + 1 | t ) is admissible. Then, for any feasible initial state Y ^ ( t | t ) , the proposed triple-mode controller Δ u ( t ) = d T ( F r R ( t + 1 | t ) K Y ^ ( t | t ) ) + c ( t | t ) asymptotically steers the system to R ( t + 1 | t ) .
The proof can be found in the Appendix A.

4. Applications

In this section, dynamic simulations are carried out to evaluate the performance of the proposed triple-mode controllers. This section presents two illustrative applications for the proposed methodology: (i) the shell fundamental control problem (case study I), and (ii) the coordinated control problem in the power plants (case study II). We aim at showing the better performance can be achieved by proper selection of future trajectory horizon. The simulation case study was performed on a PC with Intel Core i5 Quad 1.6GHz, 16GB memory and in the Matlab platform.

4.1. System I: A Heavy Oil Fractionator

In this example, consider a plant shown in Figure 3 presented in Prett and Morari [19]. in which a gaseous feed stream is separated by removing heat.
Here we will only consider the two-by-two subsystem, which consists of two inputs: the top draw u 1 and the side draw u 2 ; two outputs: the top end point y 1 and the side end point y 2 . The model for this system is:
y 1 y 2 = 4 . 05 e 30 s 50 s + 1 1.77 e 30 s 60 s + 1 5.39 e 20 s 50 s + 1 5.72 e 15 s 60 s + 1 u 1 u 2
The system will be represented using a 60-coefficient step response model with a 5-min sampling time. Set the initial operating point to be the origin and the upper bound of both output and input is 0.5 , the lower bound of both output and input is 0.5 . Three MPC algorithms are used to compare the system performance: the typical DMC, the proposed triple-mode MPC (Algorithm 2), and the DMC algorithm with an artificial selection of future target horizon (denoted as DMC-FF in the following). Tuning parameters of the MPC controllers are the same as follows: a prediction horizon of P = 60 and a control horizon of M = 8 . The weighting matrices are set as identity matrix. In the DMC-FF algorithm, the future target horizon is selected as f 1 = 40 , f 2 = 40 . In the triple-mode MPC algorithm, the optimal future target horizon can be obtained as f 1 = 13 , f 2 = 13 according to Algorithm 1. To compare the tracking performance, the evolution of setpoint is given as: the initial setpoint is coincident with the system’s initial operating point; at sample time 20, change the setpoint s p 1 of output y 1 to its upper bound and the setpoint s p 2 of output y 2 is also changed to its upper bound when sample time 40; at sample time 60, set the setpoint of output y 1 to the origin and the setpoint of output y 2 is also set to the origin at sample time 80.
Figure 4 and Figure 5 present the output response and input response of the controlled system.
As it can be seen in Figure 4, both two outputs are sluggish under the typical DMC algorithm when the setpoint changed. On the other hand, using an artificial selection of future target horizon may cause the system to response to the future target information far in advance. The performance is degraded because the typical DMC algorithm overlooks the information of future target and the DMC-FF algorithm uses too much future target information. In contrast, the triple-mode MPC algorithm which taking into account the future target information properly can respond to subtlety, the performance of the system is improved dramatically.
Next, Table 2 shows the sum of the ISE value of the controlled system for each control strategy throughout the simulation horizon. The DMC-FF using the future target information can outperform the typical DMC, and the triple-mode MPC using an optimal future target horizon achieves the best tracking performance among the three algorithms.

4.2. System II: A Power Plant Coordinated Control System

Coordinated control system (CCS) is the system that coordinates the steam turbine and the boiler at the same time to respond to the various load demands quickly [20]. The boiler and the turbine possess different dynamic characteristics. The dynamic from the combustion rate to the main steam pressure is a slow process with large delay and great inertia. In contrast, the response from the turbine governor valve to the unit output is a relatively fast process.
The simplified coordinated control system can be shown in Figure 6.
It can be seen that this system is a two-input, two-output control system, in which the manipulated variables are the fuel command u B ( % ) and the opening level of the turbine valve μ t ( % ) , the controlled variables are the main steam pressure p t ( M P a ) and the output power N E ( M W ) . Consider a 660 M W unit in Datang Panshan Power Plant, China, in which the CCS is composed of the boiler of Type H G 2023 / 17.6 Y M 4 boiler and the turbine of Type N 600 16.7 / 537 / 537 I . The nonlinear model from the load 50 % to load 100 % was identified through the step test [21]. When the load N E = 550 M W , the linearized model is:
Δ N E Δ p t = 2257 s 1 + 250 s 1 + 12 s 18 . 6 1 35 s 1 + 189 s 1 + 250 s 1 + 12 s e 23 s 0.279 1 + 250 s 0.571 1 35 s 1 + 189 s 1 + 250 s e 23 s Δ μ t Δ u B
In this case study, it was established that the control system must always obey the speed restriction on the manipulated variables | Δ μ t | 0.5 % , | Δ u B | 0.5 % , and the system is initially at N E = 550 MW, p t = 16.9 Mpa. The sampling time is T = 5 s. Three MPC algorithms are used to compare the system performance: the typical DMC, the proposed triple-mode MPC (Algorithm 2), and the DMC algorithm with an artificial selection of future target horizon (denoted as DMC-FF in the following). Tuning parameters of the three MPC algorithms are the same as follows: the prediction horizon is P = 300 and the control horizon is M = 20 ; the weighting matrices are set as: the weight for N E is q 1 = 1 and for p t is q 2 = 100 , the weight for inputs is W = 5 I . In the DMC-FF algorithm, the future trajectory horizon for outputs N E and p t are set to f 1 = 80 , f 2 = 120 respectively. In addition, the triple-mode MPC adopts f 1 = 22 , f 2 = 26 to design the feedforward based on Algorithm 1. The main purpose of the CCS is to keep the output power responding to load demands quickly when the unit running in either constant or sliding pressure operation. The following two scenarios are considered:

4.2.1. Constant Pressure Operation

When the unit is running in constant pressure operation mode, the main steam pressure is maintained at the rated value regardless of any changes in load demand.
Assuming that the evolution of load demand is given as: the initial load demand is coincident with the system’s initial operating point N E = 550 MW; at sample time 200, change the load demand to 600 MW then remain constant; at sample time 700, set the target load to 500 MW and then push it to the 550 MW at time instant 1500. In order to prevent the step load disturbance of the target load from impacting the entire control system of the unit, the setpoint of the output power is formed by limiting the load change rate after the target load is changed [22]. In general, the load change rate is determined by the operating personnel according to the operating conditions of the unit equipment. Here we set the load change rate as 5 MW/min.
Figure 7 presents the dynamics of the controlled outputs N E and p t . The figure shows the responses obtained by the typical DMC, the triple-mode MPC and the DMC-FF algorithms. The control actions implemented by the manipulated inputs μ t and u B that caused the responses in Figure 7, are shown in Figure 8.
It can be observed the output power governed by triple-mode MPC can track the load demand more quickly while maintaining steam pressure at the constant value. In contrast, the DMC-FF algorithm uses too much future target information and responses to the load demand far in advance, which leads to large oscillations in the main steam pressure. The main steam pressure controlled by the typical DMC algorithm also oscillates a lot when the load demand changes. In this case, the typical DMC algorithm even outperforms the DMC-FF algorithm due to the improper selection of future target horizon in the DMC-FF algorithm.
Table 3 reports the square errors between the output response and the desired setpoint value. In contrast, the proposed triple-mode MPC is expressed explicitly in terms of future information feedforward, and is hence expected to yield improved results for tracking performance.

4.2.2. Sliding Pressure Operation

We then proceed with the running of the unit concerning sliding pressure operation, and the optimal operating curve is plotted in Figure 9. The main steam pressure increases in line with the increase of the load when the load is less than 550 MW, and maintains constant when the load exceeds 550 MW. The load demands used for tracking are the same as before.
Figure 10 presents the dynamics of the controlled outputs N E and p t . The figure shows the responses obtained by the typical DMC, the proposed triple-mode MPC and the DMC-FF algorithms. The control actions implemented by the manipulated inputs μ t and u B that caused the responses in Figure 10, are shown in Figure 11.
Analyzing the results in the first panel of Figure 10, one can notice that for the trajectory of the controlled output N E under the proposed triple-mode MPC algorithm is closer to the desired trajectory than which under the typical DMC algorithm. In addition, the trajectory governed by the DMC-FF algorithm can also respond to the load demand more quickly than the typical DMC algorithm. (The curves of triple-mode MPC and DMC-FF almost overlap in the left magnified figure) However, the DMC-FF algorithm responses to the load demand far too early, which leads to the loss of tracking performace as we can see in the right magnified figure.
On the other hand, as shown in the second panel of Figure 10, when the setpoint of the main steam pressure changes with the load demand, the accuracy of tracking deteriorates a lot under either the typical DMC or the DMC-FF algorithms. In contrast, using the triple-mode MPC can drive the main steam pressure p t to track its setpoint without too much performance loss.
A quantitative comparison of square errors for the whole simulation is also shown in Table 4. In the table, the power output and pressure errors of the triple-mode MPC algorithm are dramatically reduced compared to those of the typical DMC and the DMC-FF algorithms.
Remark 1.
In many practical processes such as power plants [23], wind tunnels [24], there exists future target information that can be used in feedforward design. When the MPC algorithm is applied to real systems, the closed-loop performance is related to the accuracy of the model and the future target information. If the model is consistent with the plant, and the future target information is correct, the control performance obtained is as good as which obtained in the case study.

5. Conclusions

In this work, we propose a triple-mode model predictive control algorithm using future target information. In order to explicitly take into account the future target information in the MPC optimization, the proposed triple-mode control law includes the future target information feedforward term, and the future trajectory horizon can be optimized off-line to design the feedforward. Furthermore, the extra degrees of freedom are added to append the control law such that the constraints can be met when the MPC is running. Experimental results on the shell fundamental control problem, and the coordinated control problem in the power plant demonstrate that, by incorporating future target information feedforward explicitly and using an optimal future target horizon, the proposed triple-mode MPC algorithm could better hedge against the future target information, mitigate the sluggish response, and improve the tracking performance.

Author Contributions

Methodology, Z.X. and J.Z.; Writing—original draft, M.C.; Writing—review and editing, Z.X. and J.Z.; Perform the simulations, M.C.; Analysis of results, M.C., Z.X. and J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key Research and Development Program (No. 2017YFA0700300, No. 2017YFB0603703), NSFC-Zhejiang Joint Fund for the Integration of Industrialization and Informatization (No. U1809207).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of Theorem 1

Proof of Theorem 1.
A deviation model is obtained by computing an error model with respect to the future targets. To do so, consider a desired future trajectory described by R ( t + 1 | t ) , and define deviation variables
Y ¯ P ( t + 1 | t ) = Y ^ P ( t + 1 | t ) R ( t + 1 | t )
that satisfy the dynamic model
Y ¯ P ( t + 1 | t ) = Λ P Y ¯ ( t | t ) + A P M Δ U M ( t )
The zero regulation problem applied to the system in deviation variables finds Δ U M ( t ) that takes Y ¯ ( t | t ) to zero.
Feasibility: Assume that the state at the current sample time t is Y ¯ ( t | t ) . Also assume that the optimal solution is C ( t ) = [ c ( t | t ) , c ( t + 1 | t ) , , c ( t + M 1 | t ) ] T with an optimal cost V ( Y ¯ ( t | t ) , C ( t ) ) . Let Y ¯ ( t + 1 | t ) be the state at the next sampling time. Consider a control sequence
Δ U ˜ M ( t + 1 ) = F r R ( t + 2 | t + 1 ) K Y ^ ( t + 1 | t + 1 ) + C ˜ ( t + 1 )
with C ˜ ( t + 1 ) = [ c ( t + 1 | t ) , , c ( t + M 1 | t ) , 0 ] T . Then it is easy to see that Δ U ˜ M ( t + 1 ) is a feasible solution at time t + 1 due to the feasibility of the optimal solution at time t and enforced terminal constraint. Consequently, Y ^ P ( t + 1 + P | t + 1 ) = R ( t + 1 + P | t + 1 ) .
Convergence: Consider the feasible solution at time t + 1 previously presented. Following standard steps in the stability proofs of MPC (Mayne et al. [25]), we get that
V ( Y ¯ ( t + 1 | t + 1 ) , C ( t + 1 ) ) V ( Y ¯ ( t + 1 | t + 1 ) , C ˜ ( t + 1 ) ) V ( Y ¯ ( t | t ) , C ( t ) ) r ( t + 1 | t ) Y ^ ( t + 1 | t ) q 2 Δ u ( t | t ) w 2
Due to the definite positiveness of the optimal cost and its nonincreasing evolution, we infer that lim t Y ¯ P ( t + 1 | t ) Q 2 = 0 and lim t Δ U M ( t ) W 2 = 0 . Consequently, the system is steered to R ( t + 1 | t ) . □

References

  1. Qin, S.J.; Badgwell, T.A. A survey of industrial model predictive control technology. Control Eng. Pract. 2003, 11, 733–764. [Google Scholar] [CrossRef]
  2. Garcia, C.E.; Morshedi, A. Quadratic programming solution of dynamic matrix control (QDMC). Chem. Eng. Commun. 1986, 46, 73–87. [Google Scholar] [CrossRef]
  3. Rawlings, J.B. Tutorial overview of model predictive control. IEEE Control Syst. Mag. 2000, 20, 38–52. [Google Scholar]
  4. Rawlings, J.B.; Mayne, D.Q. Model Predictive Control: Theory and Design; Nob Hill Pub.: San Francisco, CA, USA, 2009. [Google Scholar]
  5. Maciejowski, J.M. Predictive Control: With Constraints; Pearson Education: London, UK, 2002. [Google Scholar]
  6. Goodwin, G.; Seron, M.M.; De Doná, J.A. Constrained Control and Estimation: An Optimisation Approach; Springer Science & Business Media: Berlin, Germany, 2006. [Google Scholar]
  7. Zhang, K.; Zhao, J.; Zhu, Y. MPC case study on a selective catalytic reduction in a power plant. J. Process Control 2018, 62, 1–10. [Google Scholar] [CrossRef]
  8. Di Carlo, J.; Wensing, P.M.; Katz, B.; Bledt, G.; Kim, S. Dynamic locomotion in the mit cheetah 3 through convex model-predictive control. In Proceedings of the 2018 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Madrid, Spain, 1–5 October 2018; pp. 1–9. [Google Scholar]
  9. Limón, D.; Alvarado, I.; Alamo, T.; Camacho, E.F. MPC for tracking piecewise constant references for constrained linear systems. Automatica 2008, 44, 2382–2387. [Google Scholar] [CrossRef]
  10. Clarke, D.W.; Mohtadi, C. Properties of generalized predictive control. Automatica 1989, 25, 859–875. [Google Scholar] [CrossRef]
  11. Middleton, R.H.; Chen, J.; Freudenberg, J.S. Tracking sensitivity and achievable performance in preview control. Automatica 2004, 40, 1297–1306. [Google Scholar] [CrossRef]
  12. Carrasco, D.S.; Goodwin, G.C. Feedforward model predictive control. Annu. Rev. Control 2011, 35, 199–206. [Google Scholar] [CrossRef]
  13. Rossiter, J.; Grinnell, B. Improving the tracking of generalized predictive control controllers. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 1996, 210, 169–182. [Google Scholar] [CrossRef]
  14. Dughman, S.; Rossiter, J. Systematic and effective embedding of feedforward of target information into MPC. Int. J. Control 2017, 1–15. [Google Scholar] [CrossRef] [Green Version]
  15. Valencia-Palomo, G.; Rossiter, J.; López-Estrada, F. Improving the feed-forward compensator in predictive control for setpoint tracking. ISA Trans. 2014, 53, 755–766. [Google Scholar] [CrossRef] [PubMed]
  16. Qin, S.J.; Badgwell, T.A. An Overview of Industrial Model Predictive Control Technology; AIChE Symposium Series; American Institute of Chemical Engineers: New York, NY, USA, 1997; Volume 93, pp. 232–256. [Google Scholar]
  17. Kiefer, J. Sequential minimax search for a maximum. Proc. Am. Math. Soc. 1953, 4, 502–506. [Google Scholar] [CrossRef]
  18. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  19. Prett, D.M.; Morari, M. The Shell Process Control Workshop; Elsevier: Amsterdam, The Netherlands, 2013. [Google Scholar]
  20. Moon, U.C.; Lee, Y.; Lee, K.Y. Practical dynamic matrix control for thermal power plant coordinated control. Control Eng. Pract. 2018, 71, 154–163. [Google Scholar] [CrossRef]
  21. Jizhen, L.; Liang, T.; Deliang, Z.; Xinping, L. Analysis on the Nonlinearity of Load-Pressure Characteristics of a 660MW Unit. Power Eng. 2005, 25, 533–536. [Google Scholar]
  22. Wood, A.J.; Wollenberg, B.F.; Sheblé, G.B. Power Generation, Operation, and Control; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  23. Flynn, D. Thermal Power Plant Simulation and Control; Number 43; IET: London, UK, 2003. [Google Scholar]
  24. Green, J.; Quest, J. A short history of the European Transonic Wind Tunnel ETW. Prog. Aerosp. Sci. 2011, 47, 319–368. [Google Scholar] [CrossRef]
  25. Mayne, D.Q.; Rawlings, J.B.; Rao, C.V.; Scokaert, P.O. Constrained model predictive control: Stability and optimality. Automatica 2000, 36, 789–814. [Google Scholar] [CrossRef]
Figure 1. The closed-loop output response with respect to different future trajectory horizon.
Figure 1. The closed-loop output response with respect to different future trajectory horizon.
Processes 08 00054 g001
Figure 2. The evolution of golden section search for future trajectory horizon of illustrated example by using Algorithm 1.
Figure 2. The evolution of golden section search for future trajectory horizon of illustrated example by using Algorithm 1.
Processes 08 00054 g002
Figure 3. Schematic diagram of “Shell” heavy oil fractionator [5]. Reproduced with permission from J.M.Maciejowski, predictive control with constraints; published by Pearson education, 2002.
Figure 3. Schematic diagram of “Shell” heavy oil fractionator [5]. Reproduced with permission from J.M.Maciejowski, predictive control with constraints; published by Pearson education, 2002.
Processes 08 00054 g003
Figure 4. Comparison of output response obtained by DMC, triple-mode MPC and DMC-FF. The top end point is contained in the first panel and the side end point is contained in the second panel.
Figure 4. Comparison of output response obtained by DMC, triple-mode MPC and DMC-FF. The top end point is contained in the first panel and the side end point is contained in the second panel.
Processes 08 00054 g004
Figure 5. Comparison of input response obtained by DMC, triple-mode MPC and DMC-FF. The top draw is contained in the first panel and the side draw is contained in the second panel.
Figure 5. Comparison of input response obtained by DMC, triple-mode MPC and DMC-FF. The top draw is contained in the first panel and the side draw is contained in the second panel.
Processes 08 00054 g005
Figure 6. Simplified block diagram of the coordinated control system.
Figure 6. Simplified block diagram of the coordinated control system.
Processes 08 00054 g006
Figure 7. Comparison of output response. The output power N E is contained in the first panel and the main steam pressure p t is contained in the second panel.
Figure 7. Comparison of output response. The output power N E is contained in the first panel and the main steam pressure p t is contained in the second panel.
Processes 08 00054 g007
Figure 8. Comparison of input response. The turbine valve μ t is contained in the first panel and the fuel command u B is contained in the second panel.
Figure 8. Comparison of input response. The turbine valve μ t is contained in the first panel and the fuel command u B is contained in the second panel.
Processes 08 00054 g008
Figure 9. Sliding pressure operation curve of the 660 MW unit.
Figure 9. Sliding pressure operation curve of the 660 MW unit.
Processes 08 00054 g009
Figure 10. Comparison of output response. The output power N E is contained in the first panel and the main steam pressure p t is contained in the second panel.
Figure 10. Comparison of output response. The output power N E is contained in the first panel and the main steam pressure p t is contained in the second panel.
Processes 08 00054 g010
Figure 11. Comparison of input response. The turbine valve μ t is contained in the first panel and the fuel command u B is contained in the second panel.
Figure 11. Comparison of input response. The turbine valve μ t is contained in the first panel and the fuel command u B is contained in the second panel.
Processes 08 00054 g011
Table 1. Closed-loop system performance with different choices of f.
Table 1. Closed-loop system performance with different choices of f.
f = 1 f = 6 f = 15
ISE0.80550.17970.3192
Table 2. Comparison of tracking error with respect to typical DMC, the triple-mode MPC and the DMC-FF algorithms.
Table 2. Comparison of tracking error with respect to typical DMC, the triple-mode MPC and the DMC-FF algorithms.
DMCTriple-Mode MPCDMC-FF
ISE6.37710.70880.9515
Table 3. Comparison of square errors under constant pressure operation.
Table 3. Comparison of square errors under constant pressure operation.
Error of N E (MW)Error of p t (MPa)
DMC431.70575.6357
Triple-mode MPC104.39300.0331
DMC-FF877.625132.3255
Table 4. Comparison of square errors under sliding pressure operation.
Table 4. Comparison of square errors under sliding pressure operation.
Error of N E (MW)Error of p t (MPa)
DMC501.757010.7935
Triple-mode MPC124.56500.1009
DMC-FF753.261122.1475

Share and Cite

MDPI and ACS Style

Chen, M.; Xu, Z.; Zhao, J. Triple-Mode Model Predictive Control Using Future Target Information. Processes 2020, 8, 54. https://0-doi-org.brum.beds.ac.uk/10.3390/pr8010054

AMA Style

Chen M, Xu Z, Zhao J. Triple-Mode Model Predictive Control Using Future Target Information. Processes. 2020; 8(1):54. https://0-doi-org.brum.beds.ac.uk/10.3390/pr8010054

Chicago/Turabian Style

Chen, Minghao, Zuhua Xu, and Jun Zhao. 2020. "Triple-Mode Model Predictive Control Using Future Target Information" Processes 8, no. 1: 54. https://0-doi-org.brum.beds.ac.uk/10.3390/pr8010054

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

Article Metrics

Back to TopTop