Unknown Inputs Nonlinear Observer for an Activated Sludge Process

This paper deals with the jointly estimation problem of unknown inputs and nonmeasured states of one altering aerated activated sludge process (ASP). In order to provide accurate and economic concentration measures during aerobic and anoxic phases, a cascade high gain observer (HGO) approach is developed. Only two concentrations are available; the other process’s states are assumed unavailable. The observer converges asymptotically and it leads to a good estimation of the unavailable states which are the ammonia and substrate concentration, as well as a quite reconstruction of the unknown inputs, which are the influent ammonia and the influent substrate concentrations. To highlight the efficiency of the proposed HGO with this MIMO system’s dynamics, simulation results are validated with experimental data.


Introduction
The environmental conservation and the biological treatment of waste water are required to preserve the ecosystems.Therefore, modeling an activated sludge (AS) Waste Water Treatment Process (WWTP) has been one important area of research for three decades.This AS wastewater treatment process is physically, chemically, and biologically very complex.Furthermore the AS is nonlinear and unsteady system; this is caused by the flow rate variations in a waste water and its composition, joined to reactions varying in time within mingled microorganisms.In order to model the biological plant of the activated sludge process, a lot of models have been suggested: Activated Sludge Process Model N1 (ASM(1)) in [1], the ASM(2) in [2], ASM(2d) in [3], and ASM(3) in [4].Due to the complexity of these models, various types of an activated sludge-plant reduced model have been put forward in the literature [5][6][7][8][9][10].In [11] the authors presented one new model to alternate activated sludge models.In fact the mean continuous model would capture the switching dynamics governing behavior (i.e., aerobic-anaerobic phases), without any consideration of any switching behavior.Recently, [12] proposed a generalized state-identification classification modeling.As a matter of fact, image processing and analysis have been utilized on the basis of linear regression modeling.
More generally, waste water treatment processes have a few measurement and control equipment.Under these circumstances various models are utilized in controller design.Among control methods applied to the bioprocess we quote robust multimodel control using the Quantitative Feedback Theory (QFT) techniques in [13], adaptive control in [14,15], and model predictive control (MPC) in [16,17] and recently networked control in [18].Nonetheless, the state variables of the ASM are not all available which present some difficulties in the eventual implementation.
In practice, for activated sludge process control and/or supervision, there are some nononline measured inputs, states, or concentrations (substrate concentration).Alternatively, these latter can also be hardly measured (ammonia and nitrate concentration): indeed the main difficulty towards the establishment of control strategies originates in general from the absence online, cheap, and reliable instrumentation.Even the existing sensors can be in many cases noisy and may require high maintenance and implementation cost.In the objective of solving such a problem, various observation methods are proposed.Among the observation method already studied we quote the following: [19] proposed for the nonlinear ASP model an extended Kalman filter, which works within the alternating aerobic-anoxic phase.The filter has been utilized for estimating states and nonstationary disturbances; [20] provided a brief summary of some results concerning state and parameter estimation approaches in relation to the chemical and biochemical processes.In [21] the authors resorted to a linear matrix inequality (LMI) to determine the values of the gain matrix of the asymptotic and classic observer, those are applied to a linearized bioprocess model.In the same framework a nonlinear observer which actually utilized both LMI and Lyapunov function was studied in [22].In [11], the authors applied one classical Luenberger observer which was dedicated to continuous nonlinear models after linearization and it allows the estimation of nitrate and ammonia concentrations.Furthermore, to better developing the control strategies in presence of unknown inputs, several researches are devoted to such observation.
Studying dynamical systems when having UIs has highly given an incentive for research activities in state estimation, fault reconstruction, and control theory.By means of one Unknown Input Observer (UIO), simultaneously estimating the system state and the UIs is achievable.In [23], the authors extended the MIMO systems triangular class.In actual fact, these systems involved nonlinear inputs.In addition, they proposed one full order high gain observer in order to conjointly estimate nonmeasured state and UIs.The authors in [24] gave an essential and sufficient condition for UIO existence.In [25], the authors worked out UIO for one nonlinear system's canonical observable form in presence of uncertainties in LMI terms.The author in [26] designed a nonlinear differential algebraic systems' observation scheme with UI.The authors in [27] extended a standard high gain observer (HGO) with the aid of a sliding-mode-based term.This latter would follow the UI vector.Extending the work of [28,29] designed a set of cascade high gain observers for many nonlinear MIMO systems.Each observer would provide the estimation of just an UI-vector component, with the exception of the last one which would really give an adjustment of full state variables.Given its utility, UIO has received big attention and has been applied to ASP, as in [30] where the authors suggested estimated states and the UIs as regards a minimized ASP nonlinear model through the use of the Extended Kalman Filter (EKF).Reference [31] treated the problems of Waste Water Treatment Plants (WWTP) sensor fault detection with the basic UIO (HGO).In fact, faults are just considered as unknown inputs.In a recent work, [32], an UIO is used to estimate jointly the states and an UI of the bioprocess linearized model.
The purpose of this paper is to estimate UIs and unavailable states of ASP.This contribution deals with the use of high gain UIO, which is relatively simple to design, can provide global or semiglobal stability results for a large class of nonlinear systems, robust to modeling uncertainty and external disturbances, and has a single adjustment parameter.The concentrations of ammonia and substrate in the inputs (  and   4  ) are very influential in the solution of the system model [33].So the continuous varying inputs may seriously affect the system control; such problem can be solved by the use of an unknown inputs observer.The paper is planned as follows.Initially we introduce a nonlinear ASP model.This model belongs to one specific nonlinear system class for which the development of a corresponding UIO is then detailed in section three; some preliminaries on the nonlinear systems class and state transformations are given in this section.Thereafter, in section four, and to facilitate applying the observation algorithm, we have proposed an appropriate repartition of ASP states.Simulation results are presented in the aim of highlighting the effectiveness of the proposed observer in curing one systematic routine so as to reconstruct jointly UIs and state variables for the ASP.A final conclusion ends the paper.

Process Description
There is a various model of ASP design.Precisely an ASP is composed of three principal components.The first component is an aeration tank.The latter will serve as a bioreactor.The second component is a settling tank (final clarifier).It was used to separate AS solids and treated waste water.The third one is a Return Activated Sludge (RAS) which is a tool for transferring settled AS from the clarifier into the aeration tank influent, (Figure 1).Added to that, the atmospheric air, or pure oxygen, will be put into a combination of primary treated sewage (waste water) which will be mixed with organisms in the purpose of developing a biological floc.In general, such a mixture is named as Mixed Liquor.Frequently, the concentration of the dry solids of mixed liquor suspended solid ranges between 3 and 6 g/L.Thereby, we control any removal efficiency by means of multiple operative conditions, for example, the aeration tank hydraulic residence time.To make it clear, the latter is determined by the division of the volume of the aeration tank over the flow rate.It is worth mentioning other factors like influent load (BOD5: Biological Oxygen Demand within 5 days, COD: Chemical Oxygen Demand, Nitrogen, etc.).The last factor is related to the AS solids which are existing in the aeration tank.We can mention here two other factors: oxygen supply and temperature.Mixed liquor will be evacuated, at the aeration tank discharge, into settling tanks.In addition to that the supernatant (treated waste water) ought to or should run off in the purpose of being discharged to a natural water.It can as well undergo further treatment before discharge.Because of biological growth, excess sludge has an eventual accumulation beyond the desired mixed liquor suspended solid concentration in the aeration tank.We remove such a solid amount from the treatment process in the target of keeping the biomass ratio to supplied food in balance.This solid amount, which is called the waste activated sludge, should be kept aside in storage tanks and should be further treated by digestion.This can be done prior to disposal under anaerobic or aerobic conditions.

Materiel and Method.
The process considered in this work is the pilot ASWWTP unit in the Engineering Laboratory of Environmental Processes (ELEP) of the National Institution of Applied Sciences (NIAS) in Toulouse, France.It comprises a unique aeration tank with a volume   = 0.03 3 equipped with an aeration surface that provides the necessary oxygen in the reactor in the goal of creating nitrification and denitrification conditions.It is obvious that the settler is a tank in which, in general, biomass is recirculated within the tank.In [34] mass balance was applied for each soluble (substrate) or particulate component (biomass).This led to an eleven state model that was subsequently reduced to four state variables.The earlier nonlinear models like ASM(1) were unattractive because of the high complex scheme systems.As a consequence, our focus is on reduced order models, which are on the bases of some biochemical considerations, giving the description of the nonlinear process behavior [33].As a result, we distinguish aerobic and anoxic phases through the use of state variables associated with four substrate concentrations.As a matter of fact, aerobic phase is the aeration period.In the latter the air is injected within the reactor in big quantities in the objective of converting the pollution as ammonium nitrate.After that, in order to transform nitrate in nitrogen, in the anoxic phase, the aeration is stopped and also an optional extraneous carbon source is added inside the mixer.The operation of altering phases can be caused by the   oxygen transfer coefficient value which conspicuously varies from 0 (concerning the anoxic phase) to higher values (in what concerns the aerobic phase) (Figure 2).

System
Model.The downsized system model is made up of four state variables which have just one alternating phase operation (alternating aeration).The latter is caused by the   oxygen transfer coefficient value changing from 0 to higher values.Tables 1 and 3 provide the process model's variable and parameters, while Table 2 gives the state variables.The process nonlinear model differential equations are expressed in the following [19]:   with 2 ,   3 ,   4 , and   are the concentrations of the dissolved oxygen, the nitrate, the ammonia, and the readily biodegradable substrate, respectively.
is the concentration of substrate soluble in water and   4 is the concentration of the ammoniacal nitrogen entering the reactor.   and   4 are considered the IUs of the system.   is an oxygenous carbon source.  and   are the dilution rates determined as where   is the input flow.In addition to that,   is an external carbon source flow.  (for i=1..5) is a simpler form of the process kinetics   .It is directly linked to the standard mode ASM (1).Moreover,   is a specific parameter of the downsized nonlinear model.

Studied Nonlinear Class and Design of Unknown Input Observer
It is worth considering the nonlinear multioutput systems.These latter can be as follows: ]  is an unknown function.The unknown input V = V 1 ∈ R  , the known input  ⊂  a compact set of R  , and the output  ∈ R  1 , (, ) ∈ R  with   (, ) ∈ R   . 1 has been related to the unknown inputs.One notates 1 ≤  1 ≤  such as System ( 1) is presented in its detailed form by the following equations: . . .
(A.2) The function Such assumption is a rank condition, which has to be opted for to estimate UIs.
Consequently, the next inequality is needed to design UIO: In particular, we need to have 3.1.State Transformation.Now, a set of coordinate transformations should be introduced so as to canonically form system (4) corresponding to the HGO synthesis.Just as it should be, an augmented system is made up of 2 blocks such that their outputs have the same values as those of the original system (4). 1 , . . .,   1 +1 and V 1 are related first-block variables.On the other hand, variables linked to the last one should be able to be simple copies of the full system states.
As a consequence, we interpret the studied immersion like an injective map  which is defined as follows: : R + → R  1 +⋅⋅⋅+  1 +1 ++ and 2 ) , with: Thus, the two subsystems in the  coordinates are given in the following: . . .
with  1 () is an unknown bounded function.
Using the adopted notation, system ( 12) can be written in the new coordinates  as follows: . . .
. . . with and 3.2.Observer Design.Our purpose is to synthesize an observer to estimate at the same time the vector of UIs as well as the nonmeasured state without the assumption of whatever model for UIs.According to Farza et al. ( 2004), a coordinate transformation is necessary for synthesizing such an observer.This may make one appropriate system form to synthesize the observer.Hence we should write systems ( 20)-( 21) as follows: with with and By referring to the framework given in [28], we opt for HGO synthesized for nonlinear MIMO systems, which are canonically nontriangular.The observer of systems ( 24)-( 25) is designed as follows: with > 0 is a real number that represent just one observer design parameter. 1 = 2 − 3;  2 = 1.
Converging the estimation error dynamic was demonstrated in detail in [28].

Application of UIO to ASP
with and and 4.2.Validation of the UIO with Experimental Data.We have taken advantage of the experiment done during the modeling and identification of the pilot unit [33] to highlight the performance of the studied observer.The states issued from the process simulation model are compared to the estimated states issued from our observer model as well as their measured values (Figures 4 and 5).The real data (measured) was carried out for the ASWWTP over 6 hours and with a sampling period of 20 minutes [33].In fact polluted water which comes from sanitary network in the city of Toulouse (France) feeds continuously the pilot unit in order to purify it.We obtained the measurement by taking samples every 20 minutes and analyze them in a biological laboratory to extract the three concentrations (oxygen, nitrate, and ammonia).
The real curves of the UIs are compared with their estimated ones.It is clear that the corresponding curves are almost superimposed.In addition, the UIs (the influent ammonia concentration   4 and the influent substrate concentration    ) are perfectly reconstructed by the UIO despite their variation.We have chosen a variable profile of the UIs to highlight the effectiveness of the UIO.Indeed, the proposed scenario reflects what is really happening in open-air treatment plants.For example, the increase in  concentration of the "influent ammonia" is due to the arrival of water loaded with heavy mass caused by industrial discharges, while decreasing in concentration of the "influential biodegradable substrate" is due to important (continuous) reversals of rain.The found results evidently show a thorough compromise between estimated and real curves.Moreover, we can show that the state variables signals of the model and the observer are rallied with the experimental measurements.This is quantified by the relative error (Table 4).The unknown inputs estimation error converges to zero.The advantages of such observer are that it is relatively simple to design and it has only one adjustment parameter.

Conclusion
A large class of nonlinear MIMO systems that involves UIs has been studied in this paper in order to synthesize UIO.This observer permits jointly estimating UIs and state variables.We have applied it to one complicated biochemical ASP.The latter is a reduced activated sludge model having switching anoxic and aerobic phases.By using only two measurements, we have achieved the reconstruction of two state variables (  and   4 ) and two UI signals (  4 and   ).Indeed, the obtained simulation results demonstrate quite estimation performances of the suggested observer compared to experimental data.Such observer should be combined with a state control law to ensure that the nitrogen concentration at the settler output does not exceed the standard European norm.

Figure 2 :
Figure 2: Combine aerobic and anoxic conditions to achieve nitrification and denitrification.

Figure 4 :
Figure 4: Estimation of dissolved oxygen, biodegradable substrate, and the influent biodegradable substrate concentration.