Publication number | US6459939 B1 |

Publication type | Grant |

Application number | US 09/340,531 |

Publication date | Oct 1, 2002 |

Filing date | Jun 29, 1999 |

Priority date | Jun 29, 1999 |

Fee status | Lapsed |

Publication number | 09340531, 340531, US 6459939 B1, US 6459939B1, US-B1-6459939, US6459939 B1, US6459939B1 |

Inventors | Alan J. Hugo |

Original Assignee | Alan J. Hugo |

Export Citation | BiBTeX, EndNote, RefMan |

Patent Citations (8), Non-Patent Citations (4), Referenced by (24), Classifications (6), Legal Events (3) | |

External Links: USPTO, USPTO Assignment, Espacenet | |

US 6459939 B1

Abstract

A method for determining the performance of model predictive controllers requiring only closed-loop data and an estimate of the process deadtime. A numerical measure is calculated which is a comparison of the actual variation of the process to the variation that would be expected if the process were controlled by a model predictive controller with no move suppression and no model mismatch. The measure accounts for the fact that the disturbance model employed in the Model Predictive Controller may be different from the actual process disturbance.

Claims(1)

1. A method of determining the performance of Model Predictive Controllers comprising:

a) first means for estimating a closed-loop response from a given set of closed-loop data;

b) second, means for estimating a process disturbance model;

c) the said process disturbance model is estimated from the said closed-loop response,

d) the said process disturbance model is estimated from said closed-loop response for time periods less than a given process deadtime,

e) third, means for determining a process response under model predictive control with no move suppression and no model mismatch from said process disturbance model and said closed-loop data;

f) fourth, means for calculating an error between the said closed-loop data and said process response under model predictive control with no move suppression and no model mismatch; and

g) fifth means for comparing an actual process variation to said error.

Description

Patents:

U.S. Pat. No. 5,838,561 Automatic control loop monitoring and diagnostics

U.S. Pat. No. 5,796,606 Process information system for distributed control systems

U.S. Pat. No. 6,208,949 Method and apparatus for dynamical system analysis

This invention relates to a technique to accurately assess the performance of Model Predictive Controllers. Specifically, the technique accounts for the known disturbance model of the controller in assessing the controller performance.

Model Predictive Controllers (also termed Dynamic Matrix Controllers, Multi-Variable Controllers, etc.) are widespread in industry and account for significant economic benefits. It is necessary to evaluate the performance of these controllers for the following reasons:

1. Evaluate current benefits of the Model Predictive Controller (MPC) application.

2. Explore future control opportunities.

3. Analyze current controller shortcomings.

Performance is largely defined it two ways: 1) How well does the controller maintain the process measurement at setpoint; and 2) How effectively does the controller determine the optimum operating conditions. Methods for evaluating the first question are the subject of this patent.

The simplest and most direct technique of determining controller performance is to calculate the variance between the measurement (PV) and the setpoint (SP), i.e.,

While simple to calculate, this measure has the extreme disadvantage that it is dependent on the level of disturbances or setpoint changes (i.e., the only reason why the SP would not equal the PV is due to disturbance/setpoint changes), and is thus more an indication of the disturbance/setpoint spectrums than the capabilities of the controller.

For this reason Harris et al. (1992) (Harris, T., and Desborough, L., Performance Assessment Measures for Univariate Feedback Control, Cdn J. of Chem. Eng., 70, pp 1186-1197, 1992) defined a Performance Index measure (η) that is independent of the disturbance/setpoint change spectrum and can be readily used to determine the actual performance of the controller. This measure compares the performance of the controller to a theoretical Minimum Variance controller (i.e., the best physically realizable linear controller), and can be calculated from routine operating data. This latter property is a strong advantage, as no plant tests are required to determine the measure. This measure is common in industry and is often referred to as the Harris Performance Index. This technology described in this patent is a modification or extension to the Harris Performance Index to account for one of the limitations of Model Predictive Controllers.

The Minimum Variance controller referenced by Harris (1992) employs an open-loop model of the process in the following form:

This is standard time series notation for process transfer functions, and says that the output y_{t }is a function of the input u_{t }and an independent white noise input a_{t}. The input passes through a linear discrete process model, while the noise input passes through a linear discrete disturbance model Some model predictive controllers use a slightly different formulation of the process model, but the results are invariant to this different formulation.

The minimum variance controller requires full specification of the disturbance model (the second term on the right hand side (RHS) of equation 1). In contrast, the majority of model predictive controllers (MPC) use the following model in their design:

Note the simplified disturbance model term. The reason for this simplified model is that it does not require parameter identification from plant tests (generally very difficult as the input sequence a_{t }is by definition unscheduled and unmeasured), and this form of the disturbance model results in robust control.

This is problematic for the Harris Performance Index. The actual disturbance entering the process is often more accurately of the form:

For processes under MPC control and with the above disturbance entering the process, the Harris Performance Index will indicate poor control, even if the process model is exact and there are no restrictions in the input movement. While this is correct (the MPC controller is in fact sluggish under these circumstances compared to a minimum variance controller), it is of not much use to the practicing control engineer.

This is because the control engineer needs to determine if the sluggishness is a result of the two factors under his control, namely controller tuning and model fidelity. The control engineer has no control over the structure of the disturbance model in the controller. As the Harris Performance Index does not differentiate between the causes of poor control, the Harris Performance Index can be of limited use, or even misleading.

While the minimum variance controller used as a reference controller by Harris et al., requires specification of a full open-loop process and disturbance model, an estimate of the performance of a minimum variance controller does not require specifying or determining any process or disturbance model. In contrast, while a Model Predictive Controller does not require specifying a disturbance model, determining the performance of a Model Predictive Controller does require that the disturbance model be specified.

There are other less sophisticated performance measurements used by industry. A common one is an autocorrelation graph of the controller errors (i.e., PV-SP), but this suffers from the same deficiency as the Harris Performance Index. A disturbance of the form of Equation 3 will result in significant autocorrelation even if the process transfer function is perfect and there are no input move restrictions.

Of concern in the present invention is to determine the performance of Model Predictive Controllers using normal closed-loop operating data Also of concern in the present invention is for the performance measure to reflect the structure of assumed disturbance model in the MPC controller.

It is, therefore, a feature of the present invention to provide a method to determine the performance of a Model Predictive Controller.

It is, also, a feature of the present invention that it only requires normal operating data to determine the controller performance.

Yet another feature of the present invention is that the performance measure adequately reflects the structure of the assumed disturbance model in the Model Predictive Controller.

It is, also, a feature of the present invention that it has the same range and interpretation as the Harris Controller Performance Index.

Additional features and advantages of the invention will be set forth in part in the description which follows, and will in part be apparent from the description, or may be learned from practice of the invention. The features and advantages of the invention may be realized by means of the combinations and steps pointed out in the appended claims.

To achieve the objectives, features, and advantages, in accordance with the purpose of the invention as embodied and broadly described herein, a method to determine the performance of a Model Predictive Controller is provided. This performance measure accounts for the fixed disturbance model of Model Predictive Controllers and may be calculated using only on-line data.

Accordingly, besides the objects and advantages of the Model Predictive Controller Performance Index described above, several objects and advantages of the present invention are:

a) to assess the performance of a Model Predictive controller using only normal operating data;

b) to obtain a measure that has the same interpretation and range of the Harris controller performance index;

c) to obtain a measure that accounts for the fixed disturbance model structure of the controller;

d) to obtain a measure that is not indicative of the actual disturbance model of the process, but indicates the process model fidelity and tuning of the controller;

e) to estimate disturbance model parameters using only closed-loop data and an estimate of the deadtime.

The algorithm for Performance Assessment of Model Predictive Controllers closely follows the work by Harris et al. (1992). The derivation of the index and notation is largely the same, with a few important exceptions to handle the limited disturbance model of the Model Predictive Controller. The first 2 steps, which are unique to this patent application, involve estimating the disturbance model from closed-loop data; the last two steps describe the modifications to the Harris Performance Index to account for the disturbance model.

1. Calculate closed-loop impulse weights (ψ_{i}), such that:

*{tilde over (y)}* _{t}=(1+ψ_{1} *B+ψ* _{2} *B* ^{2}+ . . . +ψ_{m} *B* ^{m})*a* _{t} (4)

Here B is the backwards shift operator (i.e., B{tilde over (y)}_{t}={tilde over (y)}_{t−1}), {tilde over (y)}_{t }are the mean-corrected deviations from setpoint, and a_{t }is a normally-distributed input sequence. Since these cannot be estimated directly from normal operating data, it is necessary to estimate the inverse model, i.e.,

*{tilde over (y)}* _{t}=(*a* _{1} *B+a* _{2} *B* ^{2} *+ . . . +a* _{m} *B* ^{m})*{tilde over (y)}* _{t} *+a* _{t} (5)

Normal least squares may be used to estimate the α parameters from operating data. Equation (5) may be expressed as a function of a_{t }only, i.e.,

The impulse weights vy may be determined from this equation by the formula:

The impulse weights may also be determined by other methods, for instance, non-linear least squares.

2. Determine the disturbance noise model parameters from the impulse weights. These parameters are in the general disturbance model form:

For an autoregressive-integrated-moving average (ARIMA)(1,1,0) model, the autoregressive parameter can be approximated by:

The value of s is generally the interval when the impulse parameters are greater than some confidence value, and is always less than or equal to the deadtime. The key and novel observation here is that, since the controller has no affect of the process before the deadtime, then impulse weights before the deadtime are a function of the disturbance only (and not the controller or model mismatch). Therefore, these terms represent an unbiased estimate of the disturbance impulse weights. While Harris et al. (1992) employed equations similar to Equation 6 to determine the closed-loop response before the deadtime, employing these parameters to calculate the open-loop process disturbance model has not been proposed before.

It may be shown that, for a model predictive controller (a special case of the minimum variance controller) with no move suppression or model error and a step disturbance model, that an approximation to the closed loop response of the controller would be:

The contribution of the Model Predictive Controller at and after the deadtime is:

_{3}(*B*)*a* _{t}=(ψ_{3,b} *B* ^{b}+ψ_{3,b+1} *B* ^{b+1}+ . . . +ψ_{3,m} *B* ^{m})*a* _{t} (11)

Equations 10 and 11 are novel and unique to this patent application.

3. The Model Predictive Controller (MPC) Performance Index is defined as below. First, assume the closed loop impulse weights can be divided into 3 terms

Here b is the process deadtime, and

_{1}(*B*)=1+ψ_{1,1} *B+ψ* _{1,2} *B* ^{2}+ . . . +ψ_{1,b−1} *B* ^{b−1} (13)

_{2}(*B*)=1+ψ_{2,1} *B+ψ* _{2,2} *B* ^{2}+ . . . +ψ_{2,m} *B* ^{m} (14)

Equation (13) represents the response of the process before the deadtime, and is invariant to the controller. Equations (14) and (11) represent the response after the deadtime, which is affected by the controller. Equation (14) is the contribution due to model mismatch and detuning, and Equation (11) represents the response due to the controller having the wrong disturbance model. Note that Equation (12) is similar to Equation (5) in Harris et al, (1992), with the difference that the response of the process after the deadtime is separated into two parts.

Equation (12) may be rearranged to give:

All terms on the left hand side (LHS) of this equation are known. As in Harris et al. (1992), the first m α parameters in Equation 16 may be estimated using least squares and normal operating data. In matrix form, Equation 16 is:

where

The parameter estimates are found by solving the set of equations:

As in Harris et al. (1992), the residual mean square error is given by:

And the least squares estimate for the normalized performance index η (b) is

Additional Embodiments

1) Extension to the multivariable (multiple coupled inputs and outputs) case may be done. The multivariable case for minimum variance control has been formulated (Harris, et. al., *Performance Assessment of Multivariable Controllers, *Automatica, 32, No 11, pp 1505-1518, 1996), and extending this for the MPC disturbance models would largely follow the modifications detailed in this patent.

2) Other controllers (such as Proportional-Integrating-Derivative (PID)) may also be evaluated using small modifications to this method.

3) Other controller disturbance models (for instance, a first order disturbance model) may be used rather than the step disturbance model. This would involve mainly a modification to Equation 10.

4) Feedforward controllers and setpoint change models may be incorporated into the analysis.

5) The performance index (Equation 21) may be defined differently by suitable algebraic manipulations.

6) The performance measurement may be extended to the time and frequency domain.

Conclusions, Ramifications, and Scope

Accordingly, the reader will see that the method described herein is able to measure the performance of a Model Predictive Controller so that it may be determined if it is detuned or has model mismatch. This method requires only on-line data, and is insensitive to the process disturbance or setpoint spectrums.

Although the description above contains many specificities, these should not be construed as limiting the scope of the invention but as merely providing illustrations of some of the preferred embodiments of the invention.

Thus the scope of the invention should be determined by the appended claims and their legal equivalents, rather than the derivations given.

Patent Citations

Cited Patent | Filing date | Publication date | Applicant | Title |
---|---|---|---|---|

US5347446 * | Feb 10, 1992 | Sep 13, 1994 | Kabushiki Kaisha Toshiba | Model predictive control apparatus |

US5535117 * | Mar 22, 1994 | Jul 9, 1996 | Kabushiki Kaisha Toshiba | Method and apparatus for controlling a process having a control loop using feedback control |

US5568378 * | Oct 24, 1994 | Oct 22, 1996 | Fisher-Rosemount Systems, Inc. | Variable horizon predictor for controlling dead time dominant processes, multivariable interactive processes, and processes with time variant dynamics |

US5740033 * | Oct 13, 1992 | Apr 14, 1998 | The Dow Chemical Company | Model predictive controller |

US5796606 | Feb 7, 1996 | Aug 18, 1998 | Noranda Inc. | Process information and maintenance system for distributed control systems |

US5838561 | Sep 23, 1996 | Nov 17, 1998 | Pulp And Paper Research Institute Of Canada | Automatic control loop monitoring and diagnostics |

US6208949 * | Jul 1, 1998 | Mar 27, 2001 | Adaptive Audio, Inc. | Method and apparatus for dynamical system analysis |

US6381505 * | Mar 11, 1999 | Apr 30, 2002 | Aspen Technology, Inc. | Robust steady-state target calculation for model predictive control |

Non-Patent Citations

Reference | ||
---|---|---|

1 | Cdn. J. of ChemEng 70, 1992 pp 1186-1197. | |

2 | Desborough, L.D., & Harris, T.L., 1992, Performance Assessment for Univariate Feedback Controllers. | |

3 | Harris, T.J., Seppala, C.T., & Desborough, L.D., 1999, A Review of Performance Monitoring and Assessment. | |

4 | J. of Process Control 9, 1999, pp 1-17, Techniques For Univariate and Multivariate Control Systems. |

Referenced by

Citing Patent | Filing date | Publication date | Applicant | Title |
---|---|---|---|---|

US6546358 * | Feb 9, 2000 | Apr 8, 2003 | Alan J. Hugo | Performance assessment of non-deadtime compensated controllers |

US6594593 * | Jun 16, 2000 | Jul 15, 2003 | Alan J. Hugo | Performance assessment of controllers applied to integrating processes |

US6618691 * | Aug 28, 2000 | Sep 9, 2003 | Alan J Hugo | Evaluation of alarm settings |

US6643554 * | Dec 14, 2000 | Nov 4, 2003 | Universal Dynamics Technologies Inc. | Method and apparatus for adaptive control of marginally stable systems |

US6993396 | Mar 20, 2003 | Jan 31, 2006 | John Peter Gerry | System for determining the health of process control feedback loops according to performance assessment criteria |

US7076411 * | Jan 15, 2002 | Jul 11, 2006 | National Instruments Corporation | System and method for performing a hardware-in-the-loop simulation using a plurality of graphical programs that share a single graphical user interface |

US8005575 | Jun 1, 2006 | Aug 23, 2011 | General Electric Company | Methods and apparatus for model predictive control in a real time controller |

US8527252 * | Jul 28, 2006 | Sep 3, 2013 | Emerson Process Management Power & Water Solutions, Inc. | Real-time synchronized control and simulation within a process plant |

US9163828 * | Oct 31, 2011 | Oct 20, 2015 | Emerson Process Management Power & Water Solutions, Inc. | Model-based load demand control |

US9217565 | Feb 7, 2011 | Dec 22, 2015 | Emerson Process Management Power & Water Solutions, Inc. | Dynamic matrix control of steam temperature with prevention of saturated steam entry into superheater |

US9235657 | Mar 13, 2013 | Jan 12, 2016 | Johnson Controls Technology Company | System identification and model development |

US9335042 | Aug 16, 2010 | May 10, 2016 | Emerson Process Management Power & Water Solutions, Inc. | Steam temperature control using dynamic matrix control |

US9436179 | Mar 13, 2013 | Sep 6, 2016 | Johnson Controls Technology Company | Systems and methods for energy cost optimization in a building system |

US9447963 | Feb 7, 2011 | Sep 20, 2016 | Emerson Process Management Power & Water Solutions, Inc. | Dynamic tuning of dynamic matrix control of steam temperature |

US20030132962 * | Jan 15, 2002 | Jul 17, 2003 | Santori Michael L. | System and method for performing a hardware-in-the-loop simulation using a plurality of graphical programs that share a single graphical user interface |

US20060047351 * | Aug 27, 2004 | Mar 2, 2006 | Alan Hugo | Process controller output and alarm setting evaluation |

US20080027704 * | Jul 28, 2006 | Jan 31, 2008 | Emerson Process Management Power & Water Solutions, Inc. | Real-time synchronized control and simulation within a process plant |

US20110288846 * | May 21, 2010 | Nov 24, 2011 | Honeywell International Inc. | Technique and tool for efficient testing of controllers in development (h-act project) |

US20120221124 * | May 7, 2012 | Aug 30, 2012 | Fisher-Rosemount Systems, Inc. | Using autocorrelation to detect model mismatch in a process controller |

US20130110298 * | Oct 31, 2011 | May 2, 2013 | Emerson Process Management Power & Water Solutions , Inc. | Model-based load demand control |

CN103092149B * | Oct 31, 2012 | Apr 12, 2017 | 爱默生过程管理电力和水解决方案公司 | 基于模型的发电量需求控制 |

CN104199296A * | Aug 15, 2014 | Dec 10, 2014 | 上海交通大学 | Linear regression performance evaluation method with forgetting factor |

CN105278527A * | Dec 1, 2015 | Jan 27, 2016 | 中国烟草总公司郑州烟草研究院 | Real-time performance evaluation method suitable for single loop control system in tobacco processing process |

WO2005042673A1 * | Oct 30, 2003 | May 12, 2005 | Council Of Scientific & Industrial Research | Method for predicting amount of gas produced during coal gasification |

Classifications

U.S. Classification | 700/44, 700/29, 702/182 |

International Classification | G05B13/04 |

Cooperative Classification | G05B13/048 |

European Classification | G05B13/04D |

Legal Events

Date | Code | Event | Description |
---|---|---|---|

Apr 19, 2006 | REMI | Maintenance fee reminder mailed | |

Oct 2, 2006 | LAPS | Lapse for failure to pay maintenance fees | |

Nov 28, 2006 | FP | Expired due to failure to pay maintenance fee | Effective date: 20061001 |

Rotate