Publication number | US7930044 B2 |
Publication type | Grant |
Application number | US 11/799,782 |
Publication date | Apr 19, 2011 |
Filing date | May 3, 2007 |
Priority date | Sep 7, 2006 |
Fee status | Lapsed |
Also published as | CN101523308A, CN101523308B, US20080065242, WO2008030374A2, WO2008030374A3 |
Publication number | 11799782, 799782, US 7930044 B2, US 7930044B2, US-B2-7930044, US7930044 B2, US7930044B2 |
Inventors | Fakhruddin T Attarwala |
Original Assignee | Fakhruddin T Attarwala |
Export Citation | BiBTeX, EndNote, RefMan |
Patent Citations (22), Non-Patent Citations (7), Referenced by (10), Classifications (5), Legal Events (3) | |
External Links: USPTO, USPTO Assignment, Espacenet | |
Not Applicable
Not Applicable
1. Field of Invention
This invention relates to incorporating a correction for dynamic variance of manipulated variables and controlled variables in steady state optimization in order to minimize dynamic violation of the controlled variables.
2. Background of the Invention
Optimization of a process having a plurality of independently controlled, manipulated variables and at least one controlled variable which is dependent on the manipulated variables using deterministic steady state process model and optimization method has widely been practiced in various industry including refining, chemical and pulp manufacturing.
Typically, results of the steady state optimization are applied to the (dynamic) process system having dynamic responses to the changes made. The resulting dynamic responses of the process alter the predicted steady state values of the controlled variables. The updated predicted steady state value of the controlled variables is subsequently used in the optimization method to calculate new change in the manipulated variables. Depending on the scope and time horizon, in practice, a steady state optimization may include feedback of actual performance into updating of the results of the steady state optimization on a continuing basis to a varying degrees; at one extreme, applications such as real-time process control, this feedback is automatic and well defined whereas at the other extreme applications such as plant wide production planning and control the feedback correction is not well defined and even not present. In case of the latter, essentially the steady state optimization is used to generate a future plan of actions but with very limited form of feedback correction mechanism for corrections to the future plan of actions when the actual performance of the process system deviates from the planned targets. In which case, the whole process of production planning and control is essentially an ad-hoc mix of a very detailed and comprehensive plan generation with almost non-existent feedback correction consistent with the initial plan generation. Instead, most often a new plan is generated with a few changes based on exception.
The constraint limits (low/high limits) of the variables used in the steady state optimization are also to be honored by the process system performance during the transient from current state to the final desired optimal steady state. These low/high limits are required to be not violated or in the least the violation be kept to the minimum.
In practice due to reasons of measured and unmeasured disturbances to the process system and due to the inherent dynamic characteristics of the process, the controlled variables would unavoidably violate their constraint limits. These dynamic constraint violations become serious and cause degrading of controllability of the process when one or more manipulated variables saturate by way of hitting one of its limits or of its final control valve output limits. Furthermore, because of marked changed in the response of the process when a manipulated variable saturates, the controllability of the process degrades significantly. Therefore, in practice, it becomes a difficult problem in which on one hand the optimizer is required to push up against the manipulated variables limit for increased optimality whereas on the other hand, at the manipulated variables limit, controllability of the process is adversely affected. Therefore, it is not uncommon to find that at times just when the process seemed to be at or very close to the optimal steady state, the process system seemed to become “out of control”. To counter this problem, in practice, most model predictive control implementation provide for what is generally known in the prior art as “linearization of valves output” for saturating manipulated variable. Although this technique offers some relief from the “out of control” problem, it practice, it introduces its own problem in that the “linearization” itself could be erroneous and therefore, induce dynamic oscillation of its own.
When a manipulated variable saturates, it essentially becomes ineffective for affecting the process in the desired manner. For instance, when a fuel gas flow valve saturates, process outlet temperature is not controllable to decreasing inlet temperature by manipulation of the fuel gas flow set point. Therefore, when this happens the process outlet temperature kind of behaves as if it is icicles hanging. That is, what is otherwise should be a two sided process response to its set point, at the manipulated variable saturation, the process response becomes one sided, either looking like icicles or spikes.
One of the objects of the invention presented herein is to maintain a certain degree of controllability of the process at or near manipulated variable valve saturation or its limit that would eliminate and minimize the loss of controllability problem mentioned above.
Therefore, there is a need to improve the process of steady state optimization that would reflect the dynamic variations of the variables. Another object of the invention presented herein is to control dynamic violation of low/high limits of the controlled variables in general.
One area where steady state optimization has widely been applied is in model predictive control, MPC. For the purpose of exposition, this application type will be used hereon for this invention, however, the method is equally applicable to other types of application involving steady state optimization as it relates to a process system having a plurality of manipulated variables and controlled variables with dynamic response characteristics. MPC has been widely applied in the process industry for control and optimization of complex multivariable processes. Intrinsically MPC attempts to optimize future behavior of complex multivariable processes while minimizing dynamic violation of low and high constraints (limits) of the variables.
In practice a number of different approaches have been developed and commercialized in dealing with dynamic constraints violation. U.S. Pat. Nos. 4,349,869 and 4,616,308 describe an implementation of MPC control called Dynamic Matrix Control incorporating a set of different tuning weights for controlled variables within the dynamic controller. These weights are adjusted as dynamic process value approaches the constraints limits. In another U.S. Pat. No. 6,122,555, and its commercial implementation known as RMPCT, it incorporates a funnel around the violating limits that permits the dynamic controller to determine a suitable trajectory returning the dynamic process value back towards the limit.
However, all of these implementations of MPC have one thing in common when it comes to constraint violation of the controlled variables. They all rely on employing in one form or another dynamic tuning weights in the dynamic controller of a model predictive controller to dampen down dynamic violation. They do not include any measure of correction in the steady state optimizer to compensate for likely dynamic violations of the constraints. The steady state optimizer is strictly formulated with reference to predicted steady state and the constraints satisfaction based on the assumption that the dynamic value of the controlled variables will eventually settle down to the predicted steady state values. Inherently in all these implementations there is this built-in unresolved problem of dynamic constraint resolution. The dynamic constraints resolution problem is essentially solved by way of feedback of the process response and the dynamic tuning weights. Therefore, most of these implementations do not offer a robust and sustained quality of control in which dynamic constraints violations are well behaved and remain within reasonable bounds. Instead, when dynamic constraints violations become excessive the control engineer is required to adjust the tuning weights albeit by trail and error method, not an easy task at all even in a moderate size MPC.
The steady state optimizer of most MPC implementations does not incorporate any explicit measure that would prevent dynamic constraint violation. The steady state optimizer is purely based on predicted steady state value of the variables and the steady state model gains.
An improved and more effective method of dealing with dynamic constraint violation is needed for robust and sustained performance of model predictive controllers.
Referring to
In the interest of clarity, not all features of actual implementation of a MPC controller are described in this specification.
The dynamic MPC calculation has been studied extensively (see e.g., S. J. Qin and T. A. Badgwell, “An Overview of Industrial Model Predictive Control Technology”, in Fifth International Conference on Chemical Process Control, J. C. Kantor, C. E. Garcia, and B. Caranhan, Eds., No. 93 in AIChE Symposium Series 316, 1997, pp. 232-256). The goal of the steady state MPC calculation is to recalculate the targets for the dynamic controller 511, every time the MPC controller 513 executes, because disturbances entering the system or new input information from the operator may change the location of the optimal steady state. This separation of the MPC algorithm into a steady-state and dynamic calculation was alluded to, for example, by C. Cutler, A. Morshedi, and J. Haydel, “An Industrial Perspective on Advanced Control”, AICheE National Meeting Washington, D.C., 1983 and is now common in industrial MPC technology.
Briefly, the overall system depicted in
On the other hand, the plant 501's operation is symbolically represented in
As shown in
Typically, the steady state optimizer uses a steady-state version of the dynamic process model used for the dynamic controller move calculation. The recalculated optimal steady state is then passed to the dynamic controller. The steady-state target calculation for linear MPC, represented by block 512 in
Although the issue of dynamic violation of the constraints in the dynamic optimization has been known to be a difficult problem, little progress has been made to develop simple and effective methods of preventing much of the dynamic violations.
Accordingly, the present invention relates to a new formulation of the steady-state optimization that explicitly accounts for dynamic variations of the variables in the steady-state target calculations so as to minimize and prevent dynamic violations of the constraints of the variables.
The foregoing and other features and aspects of the subject invention will be best understood with reference to a detailed description of a specific embodiment of the invention, which follows, when read in conjunction with the accompanying drawings, in which:
Basically, the overall functioning of the apparatus and the system as proposed by the invention as per
In practice, there are a number of ways in which the proposed invention can be utilized. In its simplest form, the proposed invention can be implemented as an adjunct to an existing steady state optimization control system with the addition of block 416 and a correction term for each of the inequality constraints as per the invention stated below. The invention anticipates that the dynamic variance correction can be incorporated in a variety of manners as a part of a steady state optimization of linear or non-linear formulation relating to a dynamical system comprising of at least one manipulated variable and at least one controlled variable. The specific embodiment of the invention as disclosed herein relating to model predictive control is meant to be an exemplary application and not limiting to it.
The steady state optimization of a model predictive controller is typically formulated as follows. For the constrained steady state optimization, the optimization problem is stated as follows,
J (v, p) 1.1
Max J
Subject to
v≧0 1.2.1
v_{min}≦v≦v_{max} 1.2.2
Where
v_{min }is low limit of process variables, v
v_{max }is high limit of process variables, v
p is economic price associated with v
The variables v are further classified as independent and dependent variables. Independent variables are considered inputs to the process and dependent variables are outputs of the process. Thus,
[v]=[u, y] 1.3.1
In a general case, y is describable by a non-linear function of u as
y=G(u) 1.4
and the objective function J( ) in eqn 1.1 is describable as a non-linear function of input variables u and output variables y.
For the purpose of exposition of the invention, it is assumed that y is a linear function of u and the objective function is also a linear function. However, the invention is applicable to non-linear process system with non-linear objective function as what is proposed in the invention is independent of the process system model and the objective function as will be evident from the details to follow. For the sake exposition, the invention will be described in reference to its application in model predictive control system though not limited to it in anyway.
In Model Predictive Control system with linear models, objective function J is of the form
J=Σp_{i}v_{i} 1.1.1
Where v is set of Variables within the scope of the controller
p_{i }is economic price associated with v_{i }
A linear form of model is used to describe change in y to change in u as
Δy=g(Δu) 1.4.1
where g() is process based relationship defining steady state changes in y for changes in u. One particular form of g() which is commonly used is that of linear form, that is
Δy=GΔu 1.4.2
The matrix G is generally termed as process gains matrix. Eqn 1.4.2 is a linear characterization of u and y. In a non-linear application, g() above would be more involved function of the variables.
Using eqn 1.4.2, the dependent variable y new value is calculated as follows in terms of change in the independent variables values:
y _{j} =y _{j} *+Σg _{ji} Δu 1.4.3
y _{j} =y _{j} *+Σg _{ji}(u _{i} −u _{i}*) 1.4.3.1
Where,
In case of a non-linear form of g(Δu) in eqn 1.4.1, a non-linear method of solution for y would be used. For the purpose of exposition of the invention, linear from of g() as stated by eqn 1.4.3.1 will be used hereon. Those ordinarily skilled in the art would understand as to how a non-linear form of g(Δu) can be adequately treated within the framework of the invention proposed herein.
In order to provide a way of compliance between the steady state optimization solution and the dynamic controller move to follow later, in the prior art most model predictive control would include what is generally referred to as optimization step limitations as follows:
|u_{i}−u_{i}*|≦min(OptimizationStep_{i}, 10*MaxDynamicControllerMove_{i}) 1.4.4
where,
OptimizationStep_{i }is an operator set value for independent variable i,
MaxDynamicControllerMove_{i }is an operator set value for independent variable i,
Inclusion of constraints 1.4.4 is a way to ensure the resulting solution does not force the dynamic controller calculation to become infeasible albeit in an indirect way. However, in practice it is often not easy to set values in constraint set 1.4.4 consistently. Consequently, the constraint set 1.4.4 is not really useful and in reality becomes much harder to change and maintain.
As used herein, the term “dynamic violation of the constraints” as applied to manipulated variables and controlled variables means that the variables do not violate upper and lower constraints at any time.
Even with the inclusion of constraint set 1.4.4, basically the steady state optimization part of a model predictive control in the prior art does not include explicitly any constraints directly or indirectly relating to the dynamic value of the dependent variables (or controlled variables as generally referred to).
An equivalent formulation of the above stated constrained optimization problem using change in value of the variables could be used, what is generally known as “delta” form. However, for the purpose of exposition and without loss of generality, the absolute value formulation will be used hereon.
In the prior art, there are a number of formulations for calculation of dynamic move of the manipulated variables for a model predictive control. They differ in their details in terms of the manner of dealing with dynamic interactions of the controlled variables, dynamic constraint violation and adherence of the dynamic value of the controlled variables to their desired set point as calculated by the steady state optimizer. The details of the dynamic move calculation is presented herein for the sake of providing the context within which the invention presented is most relevant but not limited to it.
For the sake of exposition of the invention presented herein, it suffices to summarize the solution of calculation of dynamic move, ΔU as
Δu=f(e, Ĝ, u*, b, c) 2.0
Where
As practiced in the art, the weight vector B and C are set by the practitioner based on practical experience of operation of the process. Nevertheless, these weights are at best based on the understanding of the process, devoid of any structured basis. Typically, they are set initially based on a few simulation studies to ascertain expected outcomes of action in different scenarios of constraint violations and conflicts. However, they are not often updated easily. Therefore, for most part these tuning values remain unchanged except when a significant problem arises in that the controller actions do not jibe with the particular process conditions. In which case, these weights are adjusted only after the fact. They will be changed once again when something else goes wrong. These weights provide only a limited degree of control over the dynamic violation of the controlled variables. It does not possess any self-adaptive property. Consequently, tuning these weights becomes a problem in combinatorial values, which for even a moderately sized problem becomes a daunting task.
The proposed invention includes in the steady state optimization, information relating to dynamic variance of the variables so as to ensure that dynamic violation of the constraints of the variables are significantly reduced and eliminated in the sense of the stochastic measure. It proposes to modify the inequality constraints set 1.2.2 as follows separately for independent variables and dependent variables; for independent variables:
u _{i} ^{min} ≦u _{i} +n _{i}δ(u _{i})≦u _{i} ^{max} 1.2.3
0≦n _{i}δ(u _{i})≦0.5(u _{i} ^{max} −u _{i} ^{min}) 1.2.3.1
where,
δ^{2}(u_{i}) is current value of variance of dynamic value of independent variable i
n_{i }is a what is herein termed as standard deviation weight for independent variable i,
and n_{i}≧0
n_{i}δ(u_{i}) in eqn 1.2.3 and its use later is termed hereon as variance correction for independent variable i.
n_{i }can be an operator set value or can be determined as a part of a constrained optimization solution in which it is maximized subject to,
0≦n_{i}≦n_{i}* 1.2.3.2
where n_{i}* is an operator set maximum value of n_{i}.
Similarly for controlled variables, y an equivalent constraint formulation would be
y _{j} ^{min} ≦y _{j} +n _{j}δ(y _{j})≦y _{j} ^{max} 1.2.4
0≦n _{j}δ(y _{j})≦0.5(y _{j} ^{max} −y _{j} ^{min}) 1.2.4.1
where,
δ^{2}(y_{j}) is current value of variance of dynamic value of dependent variable j,
n_{j }is a what is herein termed as standard deviation weight for dependent variable j,
and n_{j}≧0
n_{j}δ(y_{j}) in eqn 1.2.4 and its use later is termed hereon as variance correction for dependent variable j.
n_{j }can be an operator set value or can be determined as a part of a constrained optimization solution in which it is maximized subject to,
0≦n_{j}≦n_{j} 1.2.4.2
where n_{j}* is an operator set maximum value of n_{j}.
Note: δ in eqn 1.2.3 and 1.2.4 is what is generally known as standard deviation and termed as “sigma” commonly in the art. We will refer to δ by its generic name as sigma. Therefore, n_{i }in eqn 1.2.3.1 and hereon will be referred to as standard deviation weight or alternatively as sigma weight. The proposed invention does not restrict sigma weight to be an integer value, though use of integer value would provide a more meaningful significance. Therefore, a value of 3 or 6 would impart the meaning of 3-sigma or 6-sigma variance control respectively. Furthermore, the variance is calculated over a time horizon consistent with the time to steady state of the process and/or individual process variable. For instance, a controlled variable, the time horizon can be the maximum time to steady state for changes in one or more manipulated variables. Those ordinarily skilled in the art would appreciate that this time horizon for variance can be chosen judiciously depending on the underlying dynamic response of the variables, its frequency and the magnitude of measured as well as unmeasured disturbances. For instance, for a process with time to steady state of 120 minutes, the time horizon for variance can be some multiple of 120 minutes for all variables. Too small a time horizon would cause the controller to be unnecessarily too responsive and too high a time horizon would cause the controller not to be sufficiently responsive to dynamic violations.
An alternate more general formulation of eqn 1.2.3 and 1.2.4 can be written for u and y as follows
u _{i} ^{min} +n _{i} ^{max}δ(u _{i})≦u _{i} ≦u _{i} ^{max} −n _{i} ^{max}δ(u _{i}) 1.2.5
u _{i} ^{min } +n _{i} ^{min}δ(u _{i})≦u _{i} ^{max} 1.2.5.1
u _{i} ^{min} ≦u _{i} ^{max} −n _{i} ^{max}δ(u _{i})_{i} 1.2.5.2
u _{i} ^{min} +n _{i} ^{min}δ(u _{i})≦u _{i} ^{max} −n _{i} ^{max}δ(u _{i}) 1.2.5.3
n_{i} ^{min}≧0 1.2.5.4
n_{i} ^{max}≧0 1.2.5.5
δ(u _{i})≧0 1.2.5.6
where,
n_{i} ^{min }is considered to be standard deviation weight for minimum constraint
n_{i} ^{max }is considered to be standard deviation weight for maximum constraint
i is index for vector [u] with i=1, number of manipulated variables
n_{i} ^{min }and n_{i} ^{max }permit a asymmetrical correction of the constraints which could be useful in many instances. For example, in controlling a reactor bed temperature with a quench flow, for all practical purposes n_{i} ^{max}=0. However, for the reason of generality, eqn 1.2.5-1.2.5.6 provide a comprehensive and flexible method of specifying the variance correction under varying situations. Constraint 1.2.5.3 is important in that to prevent the corrected high/low limits do not cross. n_{i} ^{min }and n_{i} ^{max }can be an operator set values or can be determined as part of a constrained optimization solution in which their values are not to exceed the operator set values.
In similar manner, eqn 1.2.4 is generalized in the same manner as eqns 1.2.5-1.2.5.6 as eqn 1.2.6-1.2.6.6.
y _{j} ^{min} +n _{j} ^{min}δ(y _{j})≦y _{j} ≦y _{j} ^{max} −n _{j} ^{max}δ(y _{j}) 1.2.6
y _{j} ^{min} +n _{j} ^{min}δ(y _{j})≦y _{j} 1.2.6.1
y _{j} ^{min} ≦y _{j} ^{max} −n _{j} ^{max}δ(y _{j}) 1.2.6.2
y _{j} ^{min} +n _{j} ^{min}δ(y _{j})≦y _{j} ^{max} −n _{j} ^{max}δ(y _{j}) 1.2.6.3
n_{j} ^{min}≧0 1.2.6.4
n_{j} ^{max}≧0 1.2.6.5
δ(y _{j})≧0 1.2.6.6
where,
n_{j} ^{min }is considered to be standard deviation weight for minimum constraint
n_{j} ^{max }is considered to be standard deviation weight for maximum constraint
j is index for vector [y] with j=1, number of controlled variables
n_{j} ^{min }and n_{j} ^{max }permit a asymmetrical correction of the constraints which could be useful in many instances. For example, in protecting a reactor bed temperature, for all practical purposes n_{j} ^{min}=0. However, for the reason of generality, eqn 1.2.6-1.2.6.6 provide a comprehensive and flexible method of specifying the variance correction under varying situations. Constraint 1.2.6.3 is important in that to prevent the corrected high/low limits do not cross. n_{j} ^{min }and n_{j} ^{max }can be an operator set values or can be determined as part of a constrained optimization solution in which their values are not to exceed the operator set values.
For the sake of simplicity, hereon eqn 1.2.3 and 1.2.4 will be used. However, the eqn 1.2.5 and 1.2.6 could be used just as well with the exception of additional standard deviation weights for each of the variables.
At first glance it would seem that use of variance correction for independent variables would not be appropriate, as it would unnecessarily restrict the scope of optimization leading to at times infeasible solution. However, when this problem is analyzed in an overall sense and wider context, it will be apparent that in fact, counter-intuitively the use of variance correction for independent variables has its own stabilizing effect. First of all, the reduced range of the independent variables (normally considered to be the manipulated variables) would permit the size of optimization step to be reduced. This in turn would reduce size of dynamic move by the controller. The reduced dynamic move by the controller would in turn keep dynamic oscillations in the dependent variables to minimum, thereby having a stabilizing effect on the overall optimization and control of the process. As for the reduced range in the independent variables causing infeasible solution, this would not pertain to real processes. Because, real processes have what is generally known as throughput variables such as rate of feed and/or rate of energy input such as rate of fuel gas etc that can always be reduced to achieve a feasible solution. Of course, in the penultimate case, this could lead to the process operating at its minimum throughput level or possibly shutting it down completely.
Furthermore, the standard deviation weight involved with variance correction can be set at different values by the operator permitting an appropriate extent of variance correction consistent with the nature of the variable and its effect on the process response. Hence, both the independent variables and dependent variables can selectively be set up to have variance correction at different extent.
For the purpose of description of the present invention, use of variance for correction of the inequality constraints as stated above and the manner of its use is termed as variance control in steady state optimization.
In the above mentioned equations relating to variance correction, for manipulated variables, the variance term δ(u_{i}) is calculated from the actual process value of the variables based on a time horizon and frequency appropriate to the control system design and its operation. It is conceivable that the variance correction term can be construed in many different ways depending on the dynamic behavior characteristics of the variables and its impact on the controller performance. For the sake of exposition, a simple form of the variance correction term is used in describing the present invention though not limited to it.
For controlled variables, basically the variance term, δ(y_{j}) can be calculated in more than one way. One simple way is to use actual process value (dynamically varying). That is to say δ(y_{j}) is δ(y_{j}(pv)), where y_{j}(pv) refers to process value of variable y_{j}. Another way is to use the difference of process value of the controlled variable and its set point/target as calculated by the steady state optimization. That is to say δ(y_{j}) is δ(y_{j}(pe)), where y_{j}(pe) is difference of process value and set point of variable y_{j}. Again, the variance term can be calculated variously depending on the control time horizon and controller frequency. Those skilled in the art can judiciously chose the appropriate method of calculation. For the sake of exposition, hereon, the variance term for the controlled variable is based on process value of the controlled variable.
Those skilled in the art would appreciate the general applicability of the inequality of constraints formulation stated above.
Thus, eqn 1.2.5-1.2.5.6 and eqns 1.2.6-1.2.6.6 define the general form of what is termed herein as “dynamic variance correction” in model predictive control. These eqns define a dynamically changing constraint space for the variables involved.
A graphical depiction of eqn 1.2.3 and 1.2.4 is shown in
In essence, the variance correction conduits (101,102), (201,202) set up an internal dynamic low/high limits in accordance with the short-term variations of the variables as determined by actual measurements. Thus, greater the amount of the variance, greater would be the contraction of the normal limits. In practice what this means is that the steady state optimization becomes automatically and adaptively less aggressive as the dynamic controller is attempting to attain the desired steady state. This in turn would have an overall stabilizing effect on the process as a whole. The expanding and contracting dynamic constraint range allows the dynamic controller to respond to changes in the steady state targets in the presence of disturbances and the process dynamics. Intrinsically, the proposed dynamic variance correction is a way of stabilizing the process under control without unduly pushing it to the constraints by way of optimization. It provides for the optimization to back off from pushing the process too hard whilst the process is moving towards the targets and/or responding to the measured/unmeasured disturbances including any model-mismatch error effects.
Basically, the variance correction leads to what is depicted in
The use of variance correction effectively provides for an adaptive method of internal steady state optimization constraints correction for actual measured variance of the process. By setting the standard deviation weight appropriately at the value of 3 and greater provides for what is generally known in the art as 3 sigma and higher stochastic control limits. Higher value of standard deviation weight for a manipulated variable naturally improves controllability of the process at or near its saturation. Obviously, greater controllability of the process affords greater degree of operability of the process, which in turn can protect the process against sudden large disturbances.
For example, for an exothermic reaction process, at or very near saturation of the quench valves, the process is at safety risk for normal process disturbances. As in the prior art, there is no effective means for providing a margin of safety for the process to operate for normal process disturbances. Whereas, the proposed use of variance as per the invention in the steady state optimization provides an immediate and direct feedback correction to the overall optimization and control following a change in the actual short term dynamic variations in any of the variables thus providing a timely margin of safety to the operability of the process.
Tighter control of product qualities is an important criterion for consistent improved process performance. Use of sigma weight of 3 and higher in respect of the controlled variables effectively provide for 3 sigma and higher degree of quality control. In fact, in reference to critical process variables for safe operation of the process, this would provide an adaptive safety of margin of operation in response to changing dynamic process variations.
Another more general embodiment of the invention of described above relates to use of steady state optimization for operations of dynamic process system in which various resources are spent and controlled. The use of dynamic variance correction within a steady state optimization method can be applied in a variety of physical systems. The optimization may be either linear or non-linear and either linear or non-linear models can describe the process system. Those ordinary skilled in the field would appreciate how the present invention can be incorporated within such other optimization situations.
An application of the proposed invention is demonstrated by its application to a DeEthanizer process unit as depicted in
In
Case One: As per the prior art, it does not include any variance control.
Case Two: As per the invention, it includes variance control with standard deviation weight of 3 for all inequality constraints of the manipulated variables and the controlled variables.
For the purpose of exposition, performance comparison of two quality variables namely, overhead product impurity and bottom product impurity are highlighted in
In case one, both product qualities are violated constantly at their optimum targets.
In case two, both product qualities are only occasionally violated at their optimum targets. Thus, for most of the time, both product qualities remain within their specification. Occasional violations of product qualities can be further removed by increasing the standard deviation weight if desired.
This comparison study clearly demonstrates the benefit of application of the present invention with much improved performance. As evident from comparison of the performance of the exemplary process system without and with the dynamic violation correction reveal the superior performance when the dynamic violation correction as proposed by the present invention is used with the steady state optimization. Without the dynamic violation correction, both products violate their high limits continuously, whereas, with dynamic violation correction both products have very minimal dynamic violations.
Thus, with the dynamic variance correction as proposed by the invention presented herein, on-spec optimal production is achieved within the stochastic measure as per the sigma value used, whereas without the dynamic variance correction, off-spec optimal production is achieved with no stochastic measure of product spec violations.
More interestingly, without the dynamic variance correction, the fuel gas manipulated variables saturated to its high limit for the entire time period indicating that it was not available for control at all. Whereas, with the dynamic variance correction, the fuel gas exhibit a variation indicating that it was participating in control.
Furthermore, the feed rate variations without the dynamic variance correction is higher than with the dynamic variance correction, indicating that with the dynamic variance correction, not only the product specification are better controlled but also the throughput is more stable.
The proposed invention is applicable in general to any steady state optimization pertaining to a dynamic system where as a result of the optimization; changes are applied to the dynamic system resulting in dynamic responses from the system. Given that in real world most systems are of dynamic in nature, the method of dynamic variance correction can be incorporated as part of any method of optimization that pertains to its steady state condition. As stated the proposed invention applies to the steady state optimization specifically independent of how the changes to the process are calculated and applied or to that matter independent of the nature of the variables and the time frequency of the system response and the control cycle. The dynamic system can be of any form and the variables involved can be of any form and shape. There are no restrictions as to the type of variables and the system structure except that the variables are measurable directly or indirectly.
Cited Patent | Filing date | Publication date | Applicant | Title |
---|---|---|---|---|
US4349869 * | Oct 1, 1979 | Sep 14, 1982 | Shell Oil Company | Dynamic matrix control method |
US4616308 * | Dec 2, 1985 | Oct 7, 1986 | Shell Oil Company | Dynamic process control |
US5602689 * | Jun 1, 1995 | Feb 11, 1997 | Rodime Plc | Digital servo control system for use in disk drives, including sample integrity tester for reducing effects of spurious sampled position values |
US5638230 * | Jun 1, 1995 | Jun 10, 1997 | Rodime Plc | Digital servo control system for use in disk drives |
US5677809 * | Jun 1, 1995 | Oct 14, 1997 | Rodime Plc | Digital servo control system for use in disk drives, having state space observer including integral compensator |
US5680272 * | Jun 1, 1995 | Oct 21, 1997 | Rodime Plc | Digital servo control system for use in disk drives |
US5684650 * | Feb 25, 1994 | Nov 4, 1997 | Rodime, Plc | Digital servo control system for use in disk drives |
US5914830 * | Nov 25, 1997 | Jun 22, 1999 | Rodime Plc | Digital servo control system including sample integrity tester for reducing effects of spurious sampled position values |
US6122555 * | May 5, 1997 | Sep 19, 2000 | Honeywell International Inc. | System and methods for globally optimizing a process facility |
US6347018 * | May 3, 1999 | Feb 12, 2002 | Rodime Plc | Digital servo control system for use in disk drives |
US7058617 * | Sep 14, 2000 | Jun 6, 2006 | Pavilion Technologies, Inc. | Method and apparatus for training a system model with gain constraints |
US7187989 * | Dec 20, 2004 | Mar 6, 2007 | Fakhruddin T Attarwala | Use of core process models in model predictive controller |
US7194318 * | Jan 31, 2004 | Mar 20, 2007 | Fakhruddin T Attarwala | Integrated optimization and control using modular model predictive controller |
US7213006 * | Nov 4, 2005 | May 1, 2007 | Pavilion Technologies, Inc. | Method and apparatus for training a system model including an integrated sigmoid function |
US7315846 * | Apr 3, 2006 | Jan 1, 2008 | Pavilion Technologies, Inc. | Method and apparatus for optimizing a system model with gain constraints using a non-linear programming optimizer |
US7337022 * | Jun 19, 2003 | Feb 26, 2008 | Fisher-Rosemount Systems, Inc. | Constraint and limit feasibility handling in a process control system optimizer |
US7433743 * | May 27, 2002 | Oct 7, 2008 | Imperial College Innovations, Ltd. | Process control using co-ordinate space |
US20040049295 * | Jun 19, 2003 | Mar 11, 2004 | Wilhelm Wojsznis | Constraint and limit feasibility handling in a process control system optimizer |
US20060224534 * | Apr 3, 2006 | Oct 5, 2006 | Hartman Eric J | Method and apparatus for training a system model with gain constraints using a non-linear programming optimizer |
US20070168057 * | Dec 5, 2006 | Jul 19, 2007 | Fisher-Rosemount Systems, Inc. | Multi-objective predictive process optimization with concurrent process simulation |
US20080140227 * | Dec 4, 2007 | Jun 12, 2008 | Attarwala Fakhruddin T | Dynamic model predictive control |
US20100057222 * | Nov 10, 2009 | Mar 4, 2010 | Aspen Technology, Inc. | Computer method and apparatus for constraining a non-linear approximator of an empirical process |
Reference | ||
---|---|---|
1 | Cutler, C.R., "An Industrial Perspective On Advanced Control." 1983 AIChE, Annual Meeting, Washington DC. | |
2 | Garcia, Carlos E. "Quadratic Programming Solution of Dynamic Matrix Control" 1986 pp. 73-87, Gordon & Breach Science Publishers S.A. vol. 46, Chem Eng. Commun. | |
3 | Kantor, Jeffrey et al. Chemical Process Control-V, Assessment & New Directions for Research, 1997, pp. 232-256, Cache, AIChE Symposium Series 316, vol. 93. | |
4 | Lam, David K., "Optimal steady-state transitions under constrained predictive control" 2006 Elsevier B.V., 1593-1598,16th European Symposium Computer Aided Process Engineering. | |
5 | Nikolaou, Michael, Model Predictive Controllers: A Critical Synthesis of Theory and Industrial Needs, 2001, Academic Press, Advances in Chemical Engineering Series. | |
6 | Qin, S. Joe, "A survey of industrial model predictive control technology." 2002, Elsevier Science Ltd., 733-764, Control Engineering Practice 11, Aspen Technology, Inc. | |
7 | Tyagunov, Andrey A., "High-Performance Model Predictive Control for Process Industry" 2004, Library Technische Universiteit Eindhoven. |
Citing Patent | Filing date | Publication date | Applicant | Title |
---|---|---|---|---|
US8036760 * | Sep 26, 2008 | Oct 11, 2011 | Fisher-Rosemount Systems, Inc. | Method and apparatus for intelligent control and monitoring in a process control system |
US8112163 * | Mar 23, 2009 | Feb 7, 2012 | Fakhruddin T Attarwala | Embedded dynamic alarm control system |
US8271104 * | Apr 22, 2010 | Sep 18, 2012 | Terndrup Consult & Associates Bv | Method and system for dynamic optimisation of industrial processes |
US8504175 * | Jun 2, 2010 | Aug 6, 2013 | Honeywell International Inc. | Using model predictive control to optimize variable trajectories and system control |
US8594828 * | Sep 30, 2008 | Nov 26, 2013 | Rockwell Automation Technologies, Inc. | System and method for optimizing a paper manufacturing process |
US8620461 | Sep 24, 2009 | Dec 31, 2013 | Honeywell International, Inc. | Method and system for updating tuning parameters of a controller |
US20100082120 * | Apr 1, 2010 | Rockwell Automation Technologies, Inc. | System and method for optimizing a paper manufacturing process | |
US20100274368 * | Oct 28, 2010 | Terndrup Henrik | Method and system for dynamic optimisation of industrial processes | |
US20120004813 * | Jan 5, 2012 | Joachim Baumgarten | Method for adjusting at least one working unit of a self-propelled harvesting machine | |
US20130218399 * | Apr 14, 2011 | Aug 22, 2013 | Snecma | Method and system for monitoring the level of oil contained in a tank of an aircraft engine |
U.S. Classification | 700/44, 700/28 |
International Classification | G05B13/02 |
Cooperative Classification | G05B13/048 |
European Classification | G05B13/04D |
Date | Code | Event | Description |
---|---|---|---|
Nov 28, 2014 | REMI | Maintenance fee reminder mailed | |
Apr 19, 2015 | LAPS | Lapse for failure to pay maintenance fees | |
Jun 9, 2015 | FP | Expired due to failure to pay maintenance fee | Effective date: 20150419 |