• No results found

LS-SVM Identification of a Distillation Column

N/A
N/A
Protected

Academic year: 2021

Share "LS-SVM Identification of a Distillation Column"

Copied!
1
0
0

Bezig met laden.... (Bekijk nu de volledige tekst)

Hele tekst

(1)

LS-SVM Identification of a Distillation Column

Bart Huyck, Jos De Brabanter

KaHo Sint Lieven - Department Industrieel Ingenieur

Email: bart.huyck@kahosl.be jos.debrabanter@kahosl.be

Filip Logist, Jan Van Impe

K.U.Leuven - Department of Chemical Engineering (CIT)

Email: filip.logist@cit.kuleuven.be jan.vanimpe@cit.kuleuven.be

Kris De Brabanter, Bart De Moor

K.U.Leuven - Department of Electrical Engineering (ESAT - SCD)

Email: kris.debrabanter@esat.kuleuven.be bart.demoor@esat.kuleuven.be

1 Introduction

In a world where economic and environmental issues be-come more and more important, efficient knowledge of the behavior of a process has become indispensable. Mathemat-ical models are heavily exploited for the predicting of pro-cess behaviour, e.g., in view propro-cess monitoring and control. In the case of control, prediction and simulation is mostly done by linear models [1]. In the academic world, however, an evolution towards nonlinear models can be observed. For nonlinear systems, a variety of possible model structures and techniques exist , e.g., neural networks, wavelets, fuzzy models and Least Squares Support Vector Machines[2]. In this paper [3], we focus on the applicability of LS-SVMs for black-box system identification of a pilot scale binary distil-lation column. The LS-SVM models are compared to stan-dard linear techniques as transfer function models, subspace state-space models and the Box-Jenkins type models.

2 Experimental set-up

In the pilot-scale distillation set-up, four variables can be manipulated (MV): the reboiler duty, the duty of the feed heater and the reflux and feed flow rate. The controlled vari-ables (CV) of the system are the temperatures in the reboiler and at the top of the column. Only these variables are em-ployed in de modeling of the column as they are strongly related to the quality of the final products. Other measure-ments are available, but are currently omitted for control.

3 Results

Comparision of the 10-fold crossvalidation is performed on an estimation dataset for both the reboiler and the top tem-perature. For the reboiler temperature, a MISO model with lag 16 is selected based on the clear minimum of the mean squared error (MSE) of cross-validation (CV-MSE). Despite the clear minimum in CV-MSE, calculation of a fit value as defined in [4] and the MSE on a completely new vali-dation dataset, a model with lag 28 fits the measured data

tighter. The fit value changes from 80% for lag 16 to 86% for lag 28. For the top temperature a clear minimum can-not be found. This is probaly caused by an inproper model class. From lag 30 on, the CV-MSE hardly changes, so any model with a higher lag can be chosen. Based on the Akaike Information Criterion and a fit value, the model with lag 35 is selected. Comparison of the LS-SVM models with linear transfer function, state-space and polynomial models mod-els demonstrate that there is always a slight inprovement for the simulation of both temperatures for the LS-SVM model compaired to the linear models but the difference is not more than a 3 to 5% for the fit value.

4 Conclusion

The use of LS-SVMs to simulate a measured temperature compared to some well-known linear models types for a bi-nary distillation column is discussed. In this real life exam-ple, there is always a linear model describing the measured temperature very accurately for both the top as well as the reboiler temperature. The LS-SVM always compete with the best linear model, but is only slightly better.

5 Acknowledgements

Work supported in part by Projects OT/09/025/TBA, PFV/10/002 (Center-of-Excellence Optimization in Engineering), IOF-SCORES4CHEM (KP/09/005) of the Research Council of the Katholieke Universiteit Leuven and IUAP VI/4. J. Van Impe holds the chair Safety Engineering sponsored by the Belgian chemistry and life sciences federation essenscia.

References

[1] S. J. Qin and T. A. Badgwell. “A Survey of Industrial Model Pre-dictive Control Technology,” Contr Eng Pract, 11:733764, 2003

[2] J. A. K. Suykens, T. Van Gestel, J. De Brabanter, B. De Moor and J. Vandewalle “Least Squares Support Vector Machines” World Scientific Pub. Co., Singapore, 2002 (ISBN 981-238-151-1)

[3] B. Huyck, K. De Brabanter, F. Logist, J. De Brabanter, J. Van Impe and B. De Moor. “Identifcation of a Pilot Scale Distillation Column: A Kernel Based Approach” submitted to IFAC World congress Milano (2011) [4] L. Ljung. System Identification: Theory for the User, Second Edi-tion. Prentice Hall, Upper Saddle River, New Jersey, 1999.

Referenties

GERELATEERDE DOCUMENTEN

Future work for the presented method includes the exten- sion of the method to other block oriented structures like Wiener-Hammerstein systems where, after the identification of

In this paper the BLA approach is used to model the linear block and these results are used to help LS-SVM modeling the nonlinear part.. For the proposed method it will be shown

Methods dealing with the MIMO case include for instance: In [13] basis functions are used to represent both the linear and nonlinear parts of Hammerstein models; in [14], through

In this paper a new approach based on Least Squares Support Vector Machines (LS-SVMs) is proposed for solving delay differential equations (DDEs) with single-delay.. The proposed

This paper introduces a software tool SYM-LS-SVM-SOLVER written in Maple to derive the dual system and the dual model represen- tation of LS-SVM based models,

We exploit the properties of ObSP in order to decompose the output of the obtained regression model as a sum of the partial nonlinear contributions and interaction effects of the

In this paper a new approach based on Least Squares Support Vector Machines (LS-SVMs) is proposed for solving delay differential equations (DDEs) with single-delay.. The proposed

Future work for the presented method includes the exten- sion of the method to other block oriented structures like Wiener-Hammerstein systems where, after the identification of