Advanced Process Control-HYSYS
Advanced Process Control-HYSYS
• chemicals • aerospace
• food processing • metallurgy
• automotive • pulp and paper
3-33
3-34 Advanced Control
MPC Theory
Currently most model predictive control techniques like Dynamic
Matrix Control (DMC) and Model Algorithmic Control (MAC) are
based on optimization of a quadratic objective function involving
the error between the set point and the predicted outputs. In
these cases, a discrete impulse response model can be used to
derive the objective function.
Let a0, a1, a2,...,aT represent the value of the unit step response
function obtained from a typical open loop process, as shown in
the figure below:
Figure 3.16
3-34
Control Theory 3-35
ˆ (3.28)
cn = cn
T
cˆ n + 1 = c 0 + ¦ ai 'mn + 1 – i (3.29)
i=1
Figure 3.17
T
ˆ
cn + 1 = c0 + ¦ hi mn + 1 – i (3.30)
i=1
3-35
3-36 Advanced Control
T
ˆ
cn = c0 + ¦ hi mn – i (3.31)
i=1
T
ˆ ˆ
cn + 1 = cn + ¦ hi 'mn + 1 – i (3.32)
i=1
ˆ ˆ (3.33)
c n* + 1 – c n + 1 = c n – c n
T
c n* + 1 = c n + ¦ hi 'mn + 1 – i (3.34)
i=1
3-36
Control Theory 3-37
T
c n* + j = c n* + j – 1 + ¦ hi 'mn + j – i (3.35)
i=1
where:
j = 1, 2, ..., V
c n* + 1 a1 0 0 } 0 'm n cn + P1
c n* + 2 a2 a1 0 0 'm n + 1 cn + P2
c n* + 3 a3 a2 a1 0 'm n + 2 cn + P3
= + (3.36)
. . . . .
. . . . .
. . . . .
c n* + V a V a V – 1 a V – 2 } a V – U + 1 'm n + U – 1 cn + PV
where:
i
ai = ¦ hj (3.37)
j=1
i
Pi = ¦ Sj for i = 1 , 2 , }, V (3.38)
j=1
T
Si = ¦ h i 'm n + j – i for j = 1 , 2 , }, V (3.39)
i = j+1
3-37
3-38 Advanced Control
Denoting the set point trajectory, (in other words, the desired
values of the set point V time steps into the future), as
r n + j j = 1 , 2 , } , V , Equation (3.36) can be written as:
ˆ ˆ (3.40)
E = – A 'm + E'
where:
A = the V u U triangular matrix
'm = the U u 1 vector of future control moves.
ˆ and E'
E ˆ = the closed loop and open loop predictions,
respectively, and are defined as follows:
r n + 1 – c n* + 1 En – P1
rn + 2 – c n* + 2 En – P2
ˆ = . ˆ .
E E' =
. .
. .
rn + V – cn + V En – PV
–1 ˆ
'm = A E' (3.41)
3-38
Control Theory 3-39
J > 'm @ = E
ˆ
ˆ TE (3.42)
–1 ˆ ˆ
'm = A T A A T E' = K c E' (3.43)
where:
–1
A T A A T = pseudo-inverse matrix
ˆ ˆ (3.44)
J > 'm @ = E T * u E + 'm T * y 'm
where:
* u and * y are positive-definite weighting matrices for
predicted errors and control moves, respectively. These
matrices allows you to specify different penalties to be
placed on the predicted errors resulting in a better
tuned controller.
–1 ˆ ˆ
'm = A T * u A + * y A T * u E' = K c E' (3.45)
3-39
3-40 General Guidelines
Capacity
The ability of a system to attenuate incoming disturbances is a
function of the capacitance of a system and the period of the
disturbances to the system. From Terminology section,
attenuation is defined as:
K
Attenuation = 1 – ---------------------------- (3.46)
ZW 2 + 1
3-40
Control Theory 3-41
Dead Time
The dead time has no effect on attenuating disturbances to open
loop systems. However, it does have a significant negative effect
on controllability. Dead time in a process system reduces the
amount of gain the controller can implement before
encountering instability. Because the controller is forced to
reduce the gain, the process is less able to attenuate
disturbances than the same process without dead time.
3-41
3-42 General Guidelines
Figure 3.18
3-42
Control Theory 3-43
The following rules are approximate. They help you obtain tight
control. You can adjust the tuning parameters further if the
closed loop response is not satisfactory. Tighter control and
better performance can be achieved by increasing the gain.
Decreasing the controller gain results in a slower, but more
stable response.
3-43
3-44 General Guidelines
Flow Control
Flow in a pipe is typically a fast responding process. The dead
time and capacity associated with a length of pipe is generally
small. It is therefore not unusual for the process to be limited by
the final control element (valve) dynamics. You can easily
incorporate valve dynamics in the Aspen HYSYS model by
modifying the valve parameters in the Actuator page of the
Dynamics tab.
Tuning a flow loop for PI control is a relatively easy task. For the
flow measurement to track the set point closely, the gain, Kc,
should be set between 0.4 and 0.65 and the integral time, Ti,
should be set between 0.05 and 0.25 minutes.
3-44
Control Theory 3-45
This ensures that the control valve is wide open for a level of
75% and completely shut when the level is 25% for a set point
level of 50%. If PI control is desired, the liquid level controller is
typically set to have a gain, Kc, between 2 and 10. The integral
time, Ti, should be set between 1 and 5 minutes.
3-45
3-46 General Guidelines
Temperature Control
Temperature dynamic responses are generally slow, so PID
control is used. Typically, the controller gain, Kc, should be set
between 2 and 10, the integral time, Ti, should set between 2
and 10 minutes, and the derivative time Td, should be set
between 0 and 5 minutes.
Tuning Methods
An effective means of determining controller tuning parameters
is to bring the closed loop system to the verge of instability. This
is achieved by attaching a P-only controller and increasing the
gain such that the closed loop response cycles with an amplitude
that neither falls nor rises over time. At a system’s stability
margins, there are two important system parameters, the
ultimate period and the ultimate gain, which allow the
calculation of the proportional, integral, and derivative gains.
3-46
Control Theory 3-47
Figure 3.19
3-47
3-48 General Guidelines
4h
K u = ------
Sa
Ultimate Period
P u = Period taken from limit cycle
Controller Gain
Ku
K c = -------
3.2
Controller Integral Time
T i = 2.2P u
3-48
Control Theory 3-49
3-49
3-50 General Guidelines
For more information on These specifications should be reserved for variables that
setting pressure-flow
physically remain constant in a plant. For example, you can
specifications in a
dynamic simulation case, specify the exit pressure of a pressure relief valve since the exit
see Chapter 1 - pressure typically remains constant in a plant.
Dynamic Theory.
3-50
Control Theory 3-51
Connections Tab
In the Connections tab, you can specify/select the variable
information entering and exiting the controller.
3-51
3-52 General Guidelines
Operating
Description
Variable
Physical Valve It is recommended that a physical valve be used as the
operating variable for a controller. The controller’s
output signal, OP, is the desired actuator position of
the physical valve. With this setup, a more realistic
analysis of the effect of the controller on the process is
possible. Material flow through the valve is calculated
from the frictional resistance equation of the valve and
the surrounding unit operations. Flow reversal
conditions are possible and valve dynamics can be
modeled if a physical valve is selected.
Material Stream If a material stream is selected as an operating
variable, the material stream’s flow becomes a P-F
specification in the dynamic simulation case. You must
specify the maximum and minimum flow of the
material stream by clicking the Control Valve button.
The actual flow of the material stream is calculated
from the formula:
OP %
Flow = ------------------ Flow max – Flow min + Flow min
100
Aspen HYSYS varies the flow specification of the
material stream according to the calculated controller
output, OP. (Therefore, a non-realistic situation can
arise in the dynamic case since material flow is not
dependent on the surrounding conditions.)
Energy Stream If an energy stream is selected as an operating
variable, you can select a Direct Q or a Utility Fluid
Duty Source by clicking the Control Valve button.
If the Direct Q option is selected, specify the maximum
and minimum energy flow of the energy stream. The
actual energy flow of the energy stream is calculated
similarly to the material flow:
OP %
Energy Flow = ------------------ Flow max – Flow min + Flow min
100
If the Utility Fluid option is chosen, you need to specify
the maximum and minimum flow of the utility fluid.
The heat flow is then calculated using the local overall
heat transfer coefficient, the inlet fluid conditions, and
the process conditions.
3-52
Control Theory 3-53
Parameters Tab
The direction of the controller, the controller’s PV range, and the
tuning parameters can be specified in the Parameters tab.
For more information Tuning parameters are specified in the tuning field.
about the choice of tuning
parameters for each
controller, see Section
3.5.3 - Choosing
Controller Tuning
Final Control Elements
Parameters.
Set the range on the control valve at roughly twice the steady
state flow you are controlling. This can be achieved by sizing the
valve with a pressure drop between and 15 and 30 kPa with a
valve percent opening of 50%. If the controller uses a material
or energy stream as an operating variable (OP), the range of the
stream’s flow can be specified explicitly in the FCV property view
of the material or energy stream. This property view is displayed
by clicking on the Control Valve button in the PID Controller
property view.
For more information The final control element can be characterized as a linear, equal
about the
characterization of final
percentage, or quick opening valve. Control valves also have
control elements in time constants which can be accounted for in Aspen HYSYS.
Aspen HYSYS, see
Modeling Hardware
It is suggested that a linear valve mode be used to characterize
Elements section from
Section 3.3.1 - the valve dynamics of final control elements. This causes the
Available Control actual valve position to move at a constant rate to the desired
Operations.
valve positions much like an actual valve in a plant. Since the
actual valve position does not move immediately to the OP% set
by the controller, the process is less affected by aggressive
controller tuning and can possibly become more stable.
3-53
3-54 General Guidelines
Figure 3.20
From the Strip Charts tab, add a new strip chart by clicking the
Add button and activate the variables to be displayed on the
strip chart. No more than six variables should be selected for
each strip chart to keep it readable.
Figure 3.21
Click on the Strip Chart button in the View group to see the strip
chart. Size as desired and then right-click on the strip chart.
Select Graph Control command from the Object Inspect menu.
3-54
Control Theory 3-55
There are six tabs, where you can manipulate the strip chart
display features, set the numerical ranges of the strip chart for
each variable, the nature of the lines for each variable, and how
the strip chart updates and plots the data.
Figure 3.22
3-55
3-56 General Guidelines
Changing the step size causes the model to run slowly, but
during the initial switch from steady state to Dynamics mode,
the smaller step sizes allow the system to initialize better and
enable close monitoring of the controllers to ensure that
everything was set up properly.
A smaller step size also increases the stability of the model since
the solver can more closely follow changes occurring in the
plant.
Figure 3.23
3-56
Control Theory 3-57
3-57
3-58 General Guidelines
Stability
It is shown that the stability of a closed loop process depends on
the controller gain. If the controller gain is increased, the closed
loop response is more likely to become unstable. The controller
gain, Kc, input in the PID Controller operation in Aspen HYSYS is
a unitless value defined in Equation (3.47).
OP% u PV Range
K c = -------------------------------------------- (3.47)
error
3-58
Control Theory 3-59
3.6 References
1
Svrcek, Bill. A Real Time Approach to Process Controls First Edition
(1997) p. 91
2 Svrcek, Bill. A Real Time Approach to Process Controls First Edition
(1997) p. 70
3 Svrcek, Bill. A Real Time Approach to Process Controls First Edition
(1997) p. 105-123
4
Ogunnaike, B.A. and W.H. Ray. Process Dynamics, Modelling, and
Control Oxford University Press, New York (1994) p. 531
5
Seborg, D. E., T. F. Edgar and D. A. Mellichamp. Process Dynamics
and Control John Wiley & Sons, Toronto (1989) p. 649-667
3-59