Robust High-Gain Observers Based Liquid Levels and Leakage Flow Rate Estimation

The paper aims to solve the problem of liquid level and leakage flow rate estimations for a state coupled four-tank process, that is why an UIO is developed to simultaneously estimate the unmeasured state variables and the perturbations considered as unknown inputs. We have proposed a state repartition that allows putting the model of the quadruple tank system to the canonical form for which the design of the observer is more easier. The observation scheme that uses a combination of high-gain observers and sliding mode observers allows improving robustness in the state estimation quality and a perfect reconstruction of the disturbance waveforms.


Introduction
The unknown inputs (UI) estimation issue and its mathematical formulation have received considerable interest over the last two decades in several domains (secure communication [1], civil engineering [2], biomedical domain [3], chemistry [4], etc.).In fact, the estimation of the UI is required in many engineering applications and scientific studies, especially where the plants embody unknown disturbances, faults, parameters mismatch, etc.Through an unknown input observer (UIO) the aim to simultaneously estimate the unmeasured system state and the UI can be achieved.The methods used in the literature can be classified into two types: decoupling of the UI or the estimation of UI by extending the system state vector.
First alternatives mainly rely on decoupling the UI through nonlinear transformation [5,6].Consequently, strict conditions are imposed on the disturbance distribution matrix and the UI structures.The framework given in [7] discusses the UI estimation subject in different conditions.Such assumption is recently released in [8], where the authors propose a systematic design methodology for state observers for a large class of nonlinear systems with bounded exogenous inputs.Sliding mode observers-(SMO-) based UI estimation is investigated in several works [9,10].The chattering phenomena which represent the main hindrance of such approach have recently been overcome by the higher-order SMO [11].The high-gain observer-(HGO-) based approaches have been also successfully used in the conjoint estimation of state variables and UI or faults.In [12], under some global Lipschitz assumptions, a cascade HGO for a large class of nonlinear MIMO systems is designed in such a way that each subobserver provides an estimation of one component of the UI vector except the last one which achieves a reconstruction of the whole state variables.
For the second alternative, the UIs are considered as a part of the system state under the condition that their variations are relatively slow with respect to the state dynamics (constant, time polynomial, etc.) [10].The corresponding observers are then constructed in such a way to estimate both the state vector and the UIs.Furthermore, in [13,14], the authors also estimate the states and disturbance using an extended state observer (ESO).Some authors have sought to extend conventional observation algorithms such as Lunberger and EKF, so that the state vector includes the UI ( [2,15,16]).Others, however, transform the UI identification problem into a constrained optimization problem which can be easily solved by adopting the linear matrix inequalities (LMIs) formalism [17,18].
To add robustness to the quality of estimation, a combination between the HGO and the SMO algorithms is performed in some recent works.In [19] auxiliary outputs are generated using high-gain approximate differentiators and then used in the design of SMO without the match requirement for linear MIMO systems.For SISO nonlinear Lipschitz systems with nonmatching uncertainty, a hybrid observer structure that combines a HGO with higher-order sliding mode term is proposed in [20].An extension of the HGO by a sliding mode term that follows the disturbance vector for affine input nonlinear MIMO class is investigated in [21].Conjoint estimation of state and UI is required in many industrial applications, namely, in the processes where liquid level control intervenes, e.g., in food processing, water treatment systems, breeding, and pharmaceutical and petrochemical industries.The simultaneous observation of the liquid levels and the leakage flow rates is a real problem in these industrial processes where liquids are pumped, stored in tanks, and then pumped to other tanks.Besides the intrinsic nonlinearities shown in the corresponding models and the strong coupling of its states, the presence of UI such as valves perturbations or unknown flow rates present more challenges whether in the control or in diagnostic objectives.
This paper aims to solve the problem of state estimation of a quadruple tank system in presence of the UI (Figure 1).A HGO is used together with the SMO so as to improve the quality of the liquid levels estimation and reconstruct the leakage flow rates of the underlying system.A change of the state coordinates is used for transforming the original system into the canonical observable form.We introduce multiple sliding modes to handle the disturbance inputs.Under a structural assumption for the disturbance distribution matrix, the multiple SMO is designed in order to guarantee the complete observability of the system with respect to the UI.The proposed method relies only on the output estimation error in the sliding surface.The paper is organized as follows: the next section introduces the quadruple tank system.Section 3 gives some preliminaries on the nonlinear systems class of study and the state transformation.The main results on the design and analysis of a robust nonlinear observer that combines the HGO and SMO are presented in Section 4. Section 5 is devoted to the simulation results with a comparison between the Robust High-Gain Observer (R-HGO), our proposed method, and the Extended Kalman Filter (EKF) algorithm.Finally, a conclusion and perspectives are drawn.

Quadruple Tank Process Modeling
The so-called quadruple tank system, introduced in [22], has recently attracted much attention as it exhibits characteristics of interest in both control research and education [23,24].In our case, the quadruple tank process, shown in Figure 1, is a slightly modified version compared to the design given in [22].This system consists of a liquid basin, two pumps, four tanks having the same area with orifices, and level sensors at the bottom of each tank.In this experimental setup, Pump 1 and Pump 2 provide, respectively, in-feed to tanks 3 and 4 and the outflows of tank 3 and tank 4 become in-feed to tank 1 and tank 2 as shown in Figure 1.The outflows of tank 1 and tank 2 are emptied into the liquid basin.The dynamic equations for the liquid level in the four tanks issued from the Bernoulli's law are as follows: where ℎ  (),   ,    ()     () are respectively the liquid level, the cross-sectional area, the inflow rate, and the outflow rate, for the  ℎ tank.Note that both pumps are identical.So, the inflow rates into the two top tanks 3 and 4 are given by where   is the pump's constant ( 3  −1 /).The outflow rate from the orifice at the bottom of each tank is Then, the outflow rate for each top tank is given by where  is the gravitational acceleration and   denotes the cross-sectional areas of the outflow orifice at the bottom of the  ℎ tank into the  ℎ tank and for each bottom tank.
where   denotes the cross-sectional area of the outflow orifice at the bottom of the  ℎ tank into the basin.Finally, note that for the four-tank system the following equation should be respected: We have considered that the four tanks have the same crosssectional area   =  for  = 1, .
Then, we can rewrite the quadruple tank model as follows: where  1 and  2 are the unknown leakage flow rates, respectively, from the bottom tanks 1 and 2.
The main objective of this work is to simultaneously estimate the missing liquid levels in both upper tanks ℎ 3 and ℎ 4 and the UI waveforms  1 and  2 with only the measurements of the liquid in both bottom tanks ℎ 1 and ℎ 2 .

Nonlinear Class of Study and State Transformation
We consider the following class of the affine input nonlinear MIMO systems to which the quadruple tank model belongs: where  ∈ M ⊂ R  ,  ∞ connected manifold of dimension , and we assume the state space of interest M to be compact; () and   (),  = 1, . . ., , are smooth vector fields on M; ℎ  (),  = 1, . . .,  are smooth functions from M  R,  =∈ M ⊂ R  , (, ) is a vector field on M; and the disturbance vector is represented by () = [ 1 (), . . .,   ()]  with   () denoting the disturbance signals that affect the system, and we assume that each   () is bounded.The traditional nonlinear transformation uses the structural properties of the system to decouple the known/UI by transforming the system into another domain.In order to design the nonlinear UIO, the system distribution vectors  1 (), . . .,   () must satisfy the involutive property [5].The outputs should also have vector relative degree corresponding to (, ) at each point  0 ∈ M.These assumptions in general are conserved.Instead of decoupling the UIs, we shall deal with the disturbances directly in our design.In the context of our paper, the relative degree of the system is defined with respect to the UI as follows.
Assumption 1. From the  outputs, there are at least  ≥  outputs with relative degrees   = 1 with respect to the UI,  = 1, . . ., .
Remark.Among the methods that allow constructing the state transformation, one can proceed as follows.
For each output   , we define the following transformation: Thereafter, let the transformation matrix be as follows: Consequently, model ( 9) can be transformed into the new coordinates with transformation (10), so that For the subsystems under the transformations  1 , . . .,   , the following structure can be obtained: Journal of Control Science and Engineering with Each subsystem in form ( 13) can be rewritten in a condensed form as follows: where ] , So, the whole system is given in  coordinates by where and Before the synthesis of the observer, some other assumptions are required as follows.
The following assumption is the key requirement that guarantees the reconstruction of all the UI from the multiple sliding mode.Assumption 6.The dynamics of states that are measured as outputs of  subsystems have the following structure: where   1 (x) ̸ = 0,  = 1, . . ., .

Robust High-Gain Observer Design
4.1.Observer Synthesis.Our objective consists in synthesizing an observer to simultaneously estimate the unmeasured state and the UI without assuming any model for the latter.According to [21], for the subsystem (11) satisfying Assumptions 4-6, the proposed observer can be designed as follows: where for all  = 1, . . ., .V  () is a scalar-valued robust term given by the sliding mode estimation: (for  = 1, . . ., ) is the sliding mode estimation gain and   is the boundary layer design parameter.In summary, the whole proposed observer of the system (17) is given by where with is a definite positive solution of the following algebraic Lyapunov equation: It can be explicitly given as follows: with  V  = !/( − V)!V! and  > 1 is the sole design parameter.
The proof of the error convergence is detailed in [21].

Observer Form in the Original State Coordinates.
Under Assumption 1 and using (11), we can write then the observer can be written in the original state coordinates as follows: where   is given by The UIs can be reconstructed from their respectively equivalent control signals as follows: for  = 1, . . ., , where  is a small positive scalar.
The UI estimation relies on the output estimation error and hence the estimation can be performed online together with state estimation.

Application of the UIO on the Four-Tank System
To show the effectiveness of the proposed nonlinear UIO previously described, we consider the following intuitive state transformation: In the  coordinate, system (8) can be written as follows: with z1 1 = ẑ1 1 −  1 1 ; and z2 1 = ẑ2 1 −  1 1 2. The UIs can be estimated from the robust term through the multiple sliding modes (32) as follows:  By means of inverse transformation (29), the observer can be given the original coordinates as (30).
Both trapezoidal profiles for the leakage flow rates are imposed as disturbance inputs to the plant (Figure 4).The time evolution of the liquid levels ℎ  for ( = 1 ⋅ ⋅ ⋅ 4) issued from the model simulation is compared to their respective estimates provided by the observer ℎ  for ( = 1 ⋅ ⋅ ⋅ 4) (Figures 2 and 3).Notice that, with a choice of the synthesis parameter  = 1 and  1 =  2 = 30, we remark that all the estimates need less than 2 to track well their true value at the transient.For both leakage flow rate estimations   ( = 1, 2), as shown in Figure 4, their reconstruction is quite precise when they are constant, whereas a little bounded error (less than 10% ) is recorded when   varies linearly with time.

Comparison between the Robust High-Gain Observer (R-HGO) and the Extended Kalman Filter (EKF).
In order to highlight the features of the R-HGO design, besides its privilege in time computation and in the number of synthesis parameters, it is compared with the standard EKF algorithm which is one of the most industrial diffused observer [26,27].After several attempts to adjust the design parameters of this last ((0), (0), (0)), we have choose the following values: (0) = 10 −9  6 ;, (0) = 5 ⋅ 10 −8  2 , (0) = 5 ⋅ 10 −6  6 .Simulation results for conjoint state and UIs estimation of both techniques are illustrated in Figures 2-4.It is shown that when the abrupt disturbance  2 () is applied, only the estimates that arise from the R-HGO remain rallied around the trues states.This is can be explicated by the local nature of the EKF which approximates the nonlinear model only around some small neighborhood of the operating point.Moreover, through the zoom of the  1 () estimation, we remark that the EKF induces a biased reconstruction, whereas, in spite of the ripples arising from the sliding mode term in the R-HGO, the mean value of the estimated signal is more near to its truth waveform.

Conclusion
A combination of HGO and SMO is used for a conjoint estimation of the state variables and the UIs.The robust terms designed from the sliding surfaces allow preserving a little bound estimation error when the disturbance occurs.Besides, it contributes to the reconstruction of the UI waveforms.As an application, both liquid levels of the upper tanks and leakage flow rates in both bottom tanks are conjointly estimated for a quadruple tank process.Simulation results demonstrate a good estimation performance especially when the UIs are constant or vary relatively slow.
In the majority of engineering applications, measurements are collected only at the sampling instants.So, a logical extension of our work consists in the development of a continuous-discrete time UIO for the quadruple tank process.

Figure 3 :
Figure 3: Estimation of the non-measured liquid levels.
. ., 4. So, let us define the stationary parameter   as follows: { 1 , . . .,  7 ,  9 } ,  31 ,  41 ,  2 ,  32 ,  42 ,  31 +  32 ,  41 +  42 ]  8 =  10 = The objective is to reconstruct the liquid level of tanks 3 and 4 and the UIs  1 and  2 .Only the measurements of ℎ 1 and ℎ 2 are considered available.The previous assumptions are not very restrictive and they can be verified for a large class of MIMO nonlinear systems.The appropriate state observer can be rewritten in the following form: