Open access peer-reviewed chapter

Fuzzy Logic Modeling and Observers Applied to Estimate Compositions in Batch Distillation Columns

Written By

Mario Heras-Cervantes, Gerardo Marx Chávez-Campos, Héctor Javier Vergara Hernández, Adriana del Carmen Téllez-Anguiano, Juan Anzurez-Marin and Elisa Espinosa-Juárez

Submitted: September 29th, 2018 Reviewed: December 11th, 2018 Published: January 26th, 2019

DOI: 10.5772/intechopen.83479

Chapter metrics overview

812 Chapter Downloads

View Full Metrics


In this chapter, the analysis and design of a fuzzy observer based on a Takagi-Sugeno model of a batch distillation column are presented. The observer estimates the molar compositions and temperatures of the light component in the distillation column considering a binary mixture. This estimation aims to allow monitoring the physical variables in the process to improve the quality of the distillated product as well as to detect failures that could affect the system performance. The Takagi-Sugeno fuzzy model is based on eight linear subsystems determined by three premise variables: the opening percentage of the reflux valve and the liquid molar composition of the light element of the binary mixture in the boiler and in the condenser. The stability analysis and the observer gains are obtained by linear matrix inequalities (LMIs). The observer is validated by MATLAB® simulations using real data obtained from a distillation column to verify the observer’s convergence and analyze its response under system disturbances.


  • Takagi-Sugeno modeling
  • fuzzy observers
  • composition estimation
  • distillation column

1. Introduction

Distillation is the process to separate chemical substances most used in the industry, with the petrochemical (petroleum products) [1] and food (alcoholic beverages production) industries being the most important due to the current people lifestyle in which the daily use of petroleum fuels is essential for both transport and energy generation.

Fractional distillation is used to separate homogeneous liquid mixtures in which the difference between the boiling points of the components is less than 25°C. Each of the separate components is called fractions. In general, there are two operating modes: continuous and batch. In the continuous mode, the feeding of the liquid mixture and the extraction of the distilled product are carried out continuously. In the batch distillation, the mixture is initially deposited in the boiler; at the end of the process, the distillate and bottom product are extracted.

The batch operation is mainly used to separate small amounts of mixture, to obtain different qualities of the distilled product from the same mixture, or to separate multicomponent mixtures.

A batch distillation column is not operated using constant parameters; the control actions are continuously adjusted according to the state of the distillation; therefore, monitoring and controlling all the variables of the process are essential to improve the quality and quantity of the distilled product, as well as to guarantee the safety of the process and the operators. To fulfill or facilitate this objective, it is necessary to implement control techniques such as models, observers, and controllers.

In the literature, modeling and control techniques such as estimators, observers, fault detection systems, and control systems are applied to distillation columns in order to obtain a better analysis and understanding of the dynamics of the process, improving the quality of the distilled product and enhancing the user safety, among other tasks.

Distillation column simplified models present the basic principles of the process and its operation taking into account several considerations to describe the dynamics of the system in a simpler but understandable form. Authors in [2] present a simplified model of a binary distillation column, based on the liquid-vapor equilibrium of the binary mixture and the mass balance considering all the elements of the distillation column as plates.

Authors in [3] design a model based on the existence of the liquid and vapor molar fluids that vary in each column plate; the compositions of the bottom product and the distilled product are estimated using a dynamic model based on the mass and component balances. In [4] a low-order model of an ideal multicomponent distillation using the theory of nonlinear wave propagation is presented. Authors in [5] present a low-order model for a reactive multicomponent distillation column, in addition to designing a predictive control to obtain the best quality of the distilled product.

Rigorous models are more complete because they represent plate by plate the element balance of phases in each element of the distillation column (boiler, condenser, and plates). In these models, the mathematic expressions are determined by a series of differential equations given by mass, light component, or energy balances depending on the application, the control strategy, or the operation type. An important advantage of the rigorous modeling is the high resolution of the dynamics, having the disadvantage of combining a greater number of variables and expressions that make difficult the design, simulation, and implementation of controllers.

In [6], a model based on neural networks is presented in order to optimize the energy efficiency in a binary distillation column. Authors in [7] present a model of a binary distillation column based on neural networks. The neural network training and validation are performed using real data from a nine-plate pilot plant for a mixture of methanol and water. Authors in [8] present the simulation and optimization of a rigorous model for a batch reactive distillation column. Authors in [9] present the design and simulation of a discrete Kalman filter to estimate the molar compositions of the light component in a batch distillation column.

Generally, the light component composition measurement is performed offline using expensive instruments, so the implementation of state observers to estimate online this composition has become a frequent and important task. Authors in [10, 11, 12] present high-gain observers to estimate the light component composition in all the distillation column plates from the measurement of the temperature of some plates and the column actual inputs.

Due to the different distillation types and their mathematical representation, there are different types of observers for different applications. In [13] the authors present a discrete-time D-LPV observer to estimate sensors and actuators states and faults, using the H∞ approach applied to the estimation error.

In [14], the authors present a pair of extended Luenberger observers (complete and reduced order) to estimate the compositions of a multicomponent mixture from temperature measurements of the distillation column plates. The observers’ gains are calculated from the location of the closed-loop eigenvalues using a mathematical software.

In [15], a full-order nonlinear observer is presented to estimate the composition and temperatures of a distillation column. A nonlinear model obtained by the mass balance in each plate of the column is used, resulting in a set of high-order differential equations with nonlinear terms. The observer is validated in simulation to demonstrate his behavior and his robustness. The parametric representation or identification is another methodology used to estimate certain variables of the distillation columns, as presented in [16, 17].

The difficulty of designing and implementing the observers lies mainly in the nonlinear dynamics of the distillation column; thus, having a linear system would facilitate the design of observers and controllers to implement control strategies such as fault detection and diagnosis systems and automatic control and tolerant control in order to improve the performance and safety of the process, as well as the quality of the distilled product.

The Takagi-Sugeno fuzzy modeling is a tool to model and control complex systems using a nonlinear system decomposition in a multi-model structure formed by linear and not necessarily independent and fuzzy logic models [18, 19], where the representation of the nonlinear system is achieved by a weighted summation of the whole subsystems. The Takagi-Sugeno representation provides a solution to solve the problems in the design and implementation of control strategies for nonlinear systems.

In [20], the authors propose a methodology to design control techniques for systems represented in the Takagi-Sugeno form. In [21] the identification of a model of a binary distillation column, based on fuzzy models, taking into account 6 system inputs and 2 outputs for 64 rules is presented. The model is simulated using real data to validate its performance.

Authors in [22] present a controller of the molar composition of the distilled and bottom products for a binary distillation column using neural networks and fuzzy logic (ANFIS) based on a 2 × 2 MIMO system. In [23] an adaptive PID controller based on Takagi-Sugeno modeling to control the distilled and bottom products of a binary distillation column is presented.

Due to the close relationship between the fuzzy representation of nonlinear systems and the theory of linear matrix inequalities, different works based on both techniques have been developed, allowing to find solutions to the calculations corresponding to the observer and controller gains and the Lyapunov stability analysis. In [24], a methodology to design observers and controllers for a fuzzy system is proposed.

The main contribution of this work is the design of a fuzzy observer based on a Takagi-Sugeno model to estimate the molar compositions and temperatures of the light component in each plate of a binary distillation column. The observer performance is validated for applications such as system monitoring and fault detection.


2. Takagi-Sugeno fuzzy model for a batch distillation column

The objective of a batch distillation process is to separate two or more elements from a mixture, where the most volatile element is obtained as the distilled product. The equipment used to carry out the distillation process is the distillation columns or the distillation pilot plants, as shown in Figure 1.

Figure 1.

Distillation column.

A distillation column is composed by a boiler that provides the heat necessary to evaporate the mixture to be distilled; the body of the column is formed by n2plates, where a partial separation of the mixture is performed due to the phase equilibrium and a condenser that provides the necessary cooling to convert the distilled product to liquid [11].

Due to the high cost of the distillation columns and the amount of energy required in a distillation process, it is necessary to design adequate control strategies to monitor the main process variables and detect sensor and actuator faults to improve the quality of the distilled product and the safety of the process and the user.

2.1 Binary distillation column nonlinear model

The mathematical model of a distillation column consists of a set of differential equations that represent the dynamics of each plate of the column. Usually, the model of a distillation column is based on the balance of the light component in all the plates [10]. The model of a distillation column that considers a binary mixture as well as constant vapor and liquid flows through all columns is shown in Eq. (1):


where Vis the vapor molar flow, Lis the liquid molar flow, Mis the retained mass, xiis the liquid composition in plate i,and yiis the vapor composition in plate i.

The dynamic representation of a distillation column is based on submodels of the light component balance in each element of the column.

For the condenser, numbered as plate 1, Eq. (2) shows:


For the intermediate plates of the body of the distillation column, Eq. (3) shows:


where i=2,3n1.

For the boiler, numbered as plate 1, Eq. (4) shows:


where Dis the distilled product and nis the total number of plates.

2.1.1 Distillation column state-space model

The state-space dynamics of the distillation column, considering as states the light component composition in all the plates and as control inputs to the heating power, Qf,and the reflux, Rf, is expressed in Eq. (5):


where Aand Bmatrixes are defined by


where the molar fluids of liquid, L, and vapor, V, are determined by Eqs. (6) and (7):


where Hivapis the light component vapor enthalpy, Hjvapis the heavy component vapor enthalpy, and xnis the light component composition in the boiler.

The Gxiαifunction is determined by the light component composition xi, the relative volatility αi, the saturation pressure Pisat, the pressure in the column PT, the vapor molar V, and the activity coefficient γi, as expressed in Eq. (8):


2.2 Distillation column Takagi-Sugeno model

The Takagi-Sugeno model is the fuzzy representation of a nonlinear model obtained from the linear subsystems interpolation according to fuzzy rules having the form (Eq. 9):


where zitzptis the premise variables, Mi1Mp1is the fuzzy sets, ris the number of linear subsystems, xtis the state vector, utis the input vector, Aiis the linear submodel state matrices, and Biis the input vector for each subsystem.

Each consecutive linear equation represented by Aixt+Biutis called a subsystem and represents an operating point of the nonlinear system.

Given the pair xtut, the complete fuzzy model is obtained by using a singleton fuzifier, product-type inference mechanism, and center of gravity as a defuser, as described in Eq. (10):


where the vector for ppremise variables ztis defined by Eq. (11):


In addition, the calculated weight wiztfor each irule from the membership functions is defined by Eq. (12):


and the normalized weight hifor each rule is defined by Eq. (13):


2.3 Fuzzy observer

According to the structure of the fuzzy observer [24, 25, 26] expressed in Eq. (14),


The estimation error is determined by Eq. (15):


The fuzzy observer stability is demonstrated if each Ai, Cipair is observable and Pcomplies with the Lyapunov equation expressed in Eq. (16):




In [25], the demonstration that the observer is stable is presented as long as a positive definite matrix Pthat satisfies the system of linear matrix inequalities (LMIs) is found, as shown in Eq. (17):


The observer gains are defined by the LMI systems solution defined in Eq. (18):


3. Case of study: EDF-1000 distillation pilot plant

As a case of study, an EDF-1000 distillation pilot plant is used (Figure 2), consisting of 11 perforated plates, having 7 RTD (PT100) temperature sensors located in the condenser (plate 1); in plates 2, 4, 6, 8, and 10; and in the boiler (plate 11).

Figure 2.

EDF-1000 distillation pilot plant.

The boiler is composed of a 2-L tank to contain the mixture and a side tank to contain a 300-watt heating resistance.

3.1 Case of study: state-space nonlinear model

The following assumptions are considered in the EDF-1000 distillation pilot plant to simplify the designing state without significantly affecting the dynamics and precision of the model:

  • Constant pressure in the column.

  • Liquid output flows in the column.

  • No vapor retention.

  • Vapor and liquid balance in all the column plates.

  • Adiabatic distillation column.

  • Batch operation type.

Eq. (19) shows the mathematical model of a batch-type 11-plate EDF-1000 distillation pilot plant for an ethanol-water mixture, considering the compositions in each plate x1x2x10x11, where the condenser composition is x1and the boiler composition is x11; in addition, the control inputs are the heating power Qfand the reflux Rf:


where Aand Bmatrixes are defined by


where the output matrix is defined by an 11 × 11 identity matrix as shown in Eq. (20).


3.2 Case of study: Takagi-Sugeno fuzzy model

The state-space model of the case of study expressed in Eq. (10) contains nonlinearities in the A and B matrices since both are dependent on the states; in addition, the reflux input (Rf) disturbs all the states of the system (compositions) when it is activated.

The premise variables considered in the fuzzy model are the reflux valve opening percentage (z3=Rf) due to its effect in all the states of the system; the light component composition in the condenser (z2=x1), where the distillate product is obtained; and the light component composition in the boiler z1=x11, where the original mixture is located; besides, the boiler is directly related to the heating power (one of the system input variables) that provides the energy required to evaporate the mixture.

According to the steady-state dynamics of the distillation column, trapezoidal-type membership functions are chosen with two rules for each one, as expressed mathematically in Eqs. (21) and (22). This type of function is selected because it allows a greater range in the universe of discourse, where the belonging degree of belonging is 1, which prevents an oscillation when the states stabilize:


The number of subsystems of the Takagi-Sugeno fuzzy model is dependent on the number of combinations that the membership functions have; for the case of study considering three premise variables (z1=x11, z2=x1, and z3=Rf), each one with two rules (zimaxand zimin), the number of subsystems is 23=8.

Ifz1tisM1,z2tis andM3z2tisM6
Ifz1tisM1,z2tis andM4z2tisM6
Ifz1tisM2,z2tis andM3z2tisM6
Ifz1tisM2,z2tis andM4z2tisM6

where M1and M2are fuzzy sets for z1t, M3andM4are fuzzy sets for z2t, and M5and M6are fuzzy sets for z3t.

According to the fuzzy rules, the Takagi-Sugeno fuzzy model of the EDF-100 distillation pilot plant is expressed in Eq. (23):


3.3 Case of study: fuzzy observer

From the general structure of the fuzzy observer presented in Eq. (14) and the fuzzy Takagi-Sugeno model expressed in Eq. (9), the fuzzy observer for the 11-plate distillation column is expressed in Eq. (24):


The fuzzy observer scheme proposed for the distillation column of the case study is shown in Figure 3, taking as inputs the heating power and the reflux action; the estimated states are the molar compositions in all the plates, and the measured outputs are the temperatures in the condenser; in plates 2, 4, 6, 8, and 10; and in the boiler.

Figure 3.

Observer scheme.

Using the LMI system for eight rules, the LMI system for the fuzzy observer is obtained. The LMI’s characteristics of each function are expressed in Eq. (25):


The LMIs that represent the membership functions overlaps are expressed in Eq. (26):


where Pis positive definite diagonal matrix P>0of 11 × 11 dimension, and Nis an auxiliary matrix dependent on the number of states (11) and the measured outputs (7), resulting in a 7 × 11 dimension.

The resulting system of 36 LMIs is solved using the MATLAB® lmiedittool. Once the LMI system is solved, with P>0to guarantee the closed-loop stability of each subsystem, the Kigains are calculated in Eq. (27):


4. Results and discussion

The fuzzy observer validation is performed using real data from a distillation process; the main characteristics are shown in Tables 1 and 2.

EtOHvolume in the boiler1000ml
H2Ovolume in the boiler1000ml
Process total pressure662mmHg

Table 1.

Process parameters.

Density (pi)0.7891kg/m3
Molecular weight (Wi)46.06918.0528g/mol
Boiling temperature (Tbi)78.4100°C

Table 2.

Mixture parameters.

The process operates during 50 min, taking the initial compositions x0=[0.8555, 0.8525, 0.8480, 0.8412, 0.8309, 0.8148, 0.7896, 0.7483, 0.6767, 0.5369, 0.2300] in steady state.

The system input signals used in the validation stage are the reflux signal, considering an opening percentage of 20%, starting at the 5 min and lasting 30 min, and the heating power.

Figure 4 shows the light component compositions calculated by the nonlinear model in the 11 plates of the distillation column; this simulation is carried out considering no disturbances in the system.

Figure 4.

Light component composition, nonlinear model.

Figure 5 shows the composition estimated by the observer with initial conditions different from the nonlinear model. The composition initial conditions of the observer are ẋ̂0=[0.8555, 0.8525, 0.8480, 0.8412, 0.8309, 0.8148, 0.7896, 0.7483, 0.6767, 0.5369, 0.2300]. The convergence time of the observer is 48 sec.

Figure 5.

Light component estimated by the observer. Different initial conditions.

Figure 6 shows the result of the fuzzy observer simulation considering a disturbance in the plate 6 composition.

Figure 6.

Plate 6 light component estimation under disturbance.

Figure 7 shows the plate 6 composition estimated by the fuzzy observer and the comparison with the nonlinear model to verify the observer’s convergence.

Figure 7.

Plate 6 light component estimation comparison.

Figure 8 shows the simulation of the compositions estimated by the fuzzy observer under seven perturbations in the composition of plates 4, 6, and 8.

Figure 8.

Plates 4, 6, and 8 light component estimations under disturbances.

In Figure 9, a comparison between the light component compositions in plates 4, 6, and 8 is shown to verify the convergence of the fuzzy observer with the nonlinear model.

Figure 9.

Plates 4, 6, and 8 light component estimation comparison.

It can be noted that the light component composition estimation has a minimum error compared with the composition obtained from the nonlinear model using real data; besides, the convergence time is suitable for an online failure detection system or different control tasks.


5. Conclusions

In this work a fuzzy observer to estimate the molar composition of the light component in each plate of a binary distillation column is presented. The observer is based on a Takagi-Sugeno fuzzy model of eight rules, taking as premise variables the light component composition in the condenser and in the boiler, as well as the reflux signal.

The gains are calculated by means of LMIs in order to guarantee the stability for each of the closed loop linear subsystems.

The fuzzy observer is validated in simulation by using real data from an 11-plate binary distillation column for an ethanol-water mixture, considering seven RTD temperature sensors and a 300-watt heating resistor. In order to validate the fuzzy observer convergence, different tests were carried out with ideal and different initial conditions between the nonlinear system and the fuzzy observer, as well as disturbances in the light component composition in the nonlinear system.

The observer under different initial conditions has a convergence time of 45 sec. It is also shown that the observer’s convergence time under disturbances in the composition of the plate 1 is 1 sec.

Applying from two to seven perturbations to the nonlinear system, it is demonstrated that the observer is robust under multiple perturbations, enabling its implementation as a residual generator in a fault detection system.


  1. 1. Ibrahim D, Jobson M, Li J, Guilln Goslbez G. Optimization based design of crude oil distillation units using surrogate column models and a support vector machine. Chemical Engineering Research and Design. 2018;134:212-225
  2. 2. Skogestad S. Dynamics and control of distillation columns: A tutorial. Chemical Engineering Research and Design. 1997;75(6):539-562
  3. 3. Cingara A, Jovanovic M, Mitrovic M. Analytical first-order dynamic model of binary distillation column. Chemical Engineering. 1990;45(12):3585-3592
  4. 4. Kienle A. Low-order dynamic models for ideal multicomponent distillation processes using nonlinear wave propagation theory. Chemical Engineering Science. 2000;55(10):1817-1828
  5. 5. Balasubramhanya LS, Doyle FJ Iii. Nonlinear model-based control of a batch reactive distillation column. Journal of Process Control. 2000;10(2):209-218
  6. 6. Osuolale FN, Zhang J. Energy efficiency optimisation for distillation column using artificial neural network models. Energy. 2016;106:562-578
  7. 7. Singh AK, Singh HP, Mishra S. Validation of ann-based model for binary distillation column. In: Proceeding of International Conference on Intelligent Communication, Control and Devices. Singapore: Springer; 2017. pp. 235-242
  8. 8. Lopez-Saucedo ES, Grossmann IE, Segovia-Hernández JG, Hernández S. Rigorous modeling, simulation and optimization of a conventional and nonconventional batch reactive distillation column: A comparative study of dynamic optimization approaches. Chemical Engineering Research and Design. 2016;111:83-99
  9. 9. Sree Latha Chopparapu VI, George IT, Bhat VS. Design and simulation of kalman filter for the estimation of tray temperatures in a binary distillation column. International Journal of Pure and Applied Mathematics. 2017;114(9):11-18
  10. 10. Hammouri H, Targui B, Armanet F. High gain observer based on a triangular structure. International Journal of Robust and Nonlinear Control. 2002;12(6):497-518
  11. 11. Cervantes MH, Anguiano ACT, Ramrez MCG, Zaragoza CMA, Archundia ER. Validación experimental de un observador de alta ganancia reconfigurable para una columna de destilación. Revista Iberoamericana de Automtica e Informtica Industrial RIAI. 2015;12(4):397-407
  12. 12. Astorga CM, Quintero-Marmol E, Vela LG, Torres L. Observador no lineal para la estimación de concentraciones en un proceso de destilación metanol/etanol. Información tecnológica. 2006;17(6):35-41
  13. 13. López-Estrada FR, Ponsart JC, Astorga-Zaragoza CM, Camas-Anzueto JL, Theilliol D. Robust sensor fault estimation for descriptor-LPV systems with unmeasurable gain scheduling functions: Application to an anaerobic bioreactor. International Journal of Applied Mathematics and Computer Science. 2015;25(2):233-244
  14. 14. Quintero-Marmol E, Luyben WL, Georgakis C. Application of an extended Luenberger observer to the control of multicomponent batch distillation. Industrial & Engineering Chemistry Research. 1991;30(8):1870-1880
  15. 15. Lang L, Gilles ED. Nonlinear observers for distillation columns. Computers & Chemical Engineering. 1990;14(11):1297-1301
  16. 16. Zhu Y. Distillation column identification for control using Wiener model. IEEE. 1999;5:3462-3466
  17. 17. Safdarnejad M, Gallacher JR, Hedengren JD. Dynamic parameter estimation and optimization for batch distillation. Computers and Chemical Engineering. 2016;86:18-32
  18. 18. Angelov PP, Filev DP. An approach to online identification of Takagi-Sugeno fuzzy models. IEEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics). 2004;34(1):484-498
  19. 19. Takagi T, Sugeno M. Fuzzy Identification of Systems and its Applications to Modeling and Control. IEEE transactions on systems, man, and cybernetics. IEEE. 1985;1:116-132
  20. 20. Wang HO, Tanaka K, Griffin M. Parallel distributed compensation of nonlinear systems by takagi-sugeno fuzzy model. IEEE. 1992;2:531-538
  21. 21. Sanandaji BM, Salahshoor K, Fatehi A. Multivariable GA-based identification of TS fuzzy models: MIMO distillation column model case study. In: de n Fuzzy Systems Conference, 2007. FUZZ-IEEE 2007;2007
  22. 22. Sivakumar R, Balu K. ANFIS based distillation column control. In International Journal of Computer Applications Special issue on Evolutionary Computation. 2010;2:67-73
  23. 23. Mishra P, Kumar V, Rana KPS. A fractional order fuzzy PID controller for binary distillation column control. Expert Systems with Applications. 2015;42(22):8533-8549
  24. 24. Tanaka K, Sugeno M. Stability analysis and design of fuzzy control systems. Fuzzy Sets and Systems. 2004;45(2):135-156
  25. 25. Tanaka K, Ikeda T, Wang HO. Fuzzy regulators and fuzzy observers: Relaxed stability conditions and LMI-based designs. IEEE Transactions on Fuzzy Systems. 1998;6(2):250-265
  26. 26. Cao Y-Y, Frank PM. Analysis and synthesis of nonlinear time-delay systems via fuzzy control approach. IEEE Transactions on Fuzzy Systems. 2000;8(2):200-211

Written By

Mario Heras-Cervantes, Gerardo Marx Chávez-Campos, Héctor Javier Vergara Hernández, Adriana del Carmen Téllez-Anguiano, Juan Anzurez-Marin and Elisa Espinosa-Juárez

Submitted: September 29th, 2018 Reviewed: December 11th, 2018 Published: January 26th, 2019