Model Predictive Control of Industrial Processes
Model Predictive Control of Industrial Processes
processes
Vitali Vansovitš
Contents
2 (33)
Industrial process (Iru Power Plant)
3 (33)
Industrial process control
Two cascade control loop
Cascade control loop
Control loop
fast
slow
slowest
4 (33)
Industrial process control (Iru Power Plant)
C3: kp = 0,25; ti = 20
C2: kp = 0,9; ti = 240
C1: kp = 0,5; ti = 180
5 (33)
Industrial process (Iru Power Plant)
84
83
82
P. output SP
81
80 P. output ME
79
78
0:00 1:12 2:24 3:36 4:48 6:00 7:12
130
125
120
Setpoint
115
110 Measurement
105
0:00 1:12 2:24 3:36 4:48 6:00 7:12
3000
2800
2600
2400 Water flow
2200
2000
0:00 1:12 2:24 3:36 4:48 6:00 7:12
6 (33)
Neural network model
Recursive prediction
126
124
122
Temperature, deg
120
118
116
Temperature, deg
114 Predicted temperature, deg
112
0 100 200 300 400 500 600 700 800 900
Time, min
8 (33)
Neural network first results
5 minutes predictions were good
Simulation was not possible
After 10 cycles model output was not close enough to real
process output
After 100 cycles model output was already somewhere in
outer space
So, it was decided to start from simple things – linear
model identification with conventional methods.
9 (33)
Process identification
Process is simple, but nonlinear:
Q r m fuel r m fuel
T
Q c mwater T c mwater
10 (33)
Polynomial model identification
11 (33)
State-space model
12 (33)
Process model
State space model with 4 states (ss2) was selected
as:
it is simpler than oe221 (6 states in state space representation)
shows better results then ss1 (1 state)
13 (33)
MPC (Model predictive controller)
Model identification
To avoid perfect “theoretical results” model for MPC was
identified from the data acquired one year later at the
time when boiler load was different (~1/8 of previous
example) and nonconstant.
14 (33)
MPC structure
15 (33)
MPC observer
xˆ k | k xˆ k | k 1 L y k yˆ k | k 1
xˆ k 1 | k Axˆ k | k Bu k
yˆ k | k 1 Cxˆ k | k 1
xˆ k 1 | k AI L C xˆ k | k 1 Bu k AL y k
ek xk xˆ k | k 1 ek 1 A LC e k
16 (33)
MPC design (1)
2 2
There is a cost function: V k Z k T k Q
U k R
where
Z(k) – predicted process output
T(k) – reference trajectory
U(k) – process input changes
Q, R – weight matrices
2
a A
a T Aa
Z k xk uk 1 U k
k k xk uk 1
where G 2 T
QE k
T
H Q R
18 (33)
MPC constraints
Constraints are in the form:
U k U k Z k
E 0, F 0, G 0
1 1 1
T
minimize U k H U k G T U k
subject to inequality constraint.
Quadratic programming problem
19 (33)
MPC for water boiler
1 manipulated process input (gas flow)
2 measured process disturbances (water & air flows)
1 measured process output
20 (33)
MPC for water boiler (Simulink)
21 (33)
Outputs comparison
22 (33)
Control quality
23 (33)
Gas consumption
24 (33)
Earn some money
Same control quality, but less raw materials
25 (33)
Setpoint decreased 0,6 (1)
Output comparison
26 (33)
Setpoint decreased 0,6 (2)
Gas consumption difference
~8000 m3
27 (33)
Setpoint decreased 0,6 (3)
Sum of errors
28 (33)
Very rough estimation
Gas price 400 € / 1000 m3
Difference 8000 m3
Period 16 hours
Season length 50 days
Economy
= 50 * (24/16) * 8000 * 400 / 1000
= 240000 € / season
29 (33)
Some results
30 (33)
Future plans
31 (33)
References
Oliver Nelles, Nonlinear System Identification,
Springer, Berlin, 2001
J.M.Maciejowski, Predictive Control with Constraints,
pdf book, 2000
Real life
32 (33)
Thank you!
Questions???