Réduction du modèle ASM 1 pour la commande optimale des petites stations d'épuration à boues activées Reduction of the ASM 1 model for optimal control of small-size activated sludge treatment plants

Résumé

L'adoption par l'Union Européenne de normes de rejets plus contraignantes implique une meilleure gestion des stations d'épuration. L'utilisation de modèles de simulation dynamique dans des schémas de commande en boucle fermée constitue une alternative intéressante pour répondre à ce problème.

Sur la base du modèle ASM 1, un modèle réduit est ici élaboré pour le procédé à boues activées en aération séquentielle, en vue de la commande optimale du système d'aération. Les simplifications considérées sont de deux types : (i) les dynamiques lentes du système sont identifiées au moyen d'une méthode d'homotopie, puis éliminées du modèle ; (ii) des simplifications plus heuristiques consistant à prendre en compte un composé organique unique et à éliminer la concentration des composés organiques azotés sont ensuite appliquées. Elles conduisent à un modèle simplifié de 5 variables. L'application d'une procédure d'identification paramétrique permet alors de démontrer que le comportement dynamique du modèle simplifié est en bonne adéquation avec celui du modèle ASM 1 sur un horizon de prédiction de plusieurs heures, même lorsque les concentrations de l'influent ne sont pas connues. Il est également vérifié que le modèle proposé est observable et structurellement identifiable, sous des conditions d'aérobiose et d'anoxie, à partir des mesures en ligne des concentrations en oxygène dissous, ammoniaque et nitrate.

Le modèle simplifié développé présente ainsi toutes les propriétés requises pour une future utilisation au sein de schémas de commande en boucle fermée, en vue de la commande optimale des petites stations d'épuration à boues activées.

Mots-clés:

• Traitement des eaux usées,
• procédé à boues activées,
• élimination des nutriments

Abstract

In order to meet the stricter wastewater effluent guidelines adopted by the European Union, wastewater treatment plants require better management strategies. Wastewater treatment process models have become a major tool to design closed-loop control schemes. However, the dynamic models that are currently used in the simulation of activated sludge treatment plants (ASM 1, ASM 2 and, more recently, ASM 3 models) are highly dimensional and are not appropriate for on-line implementation (e.g., for model predictive control or optimal control). It is therefore important to develop reduced models that could be used for this purpose.

A reduced model was developed to describe the behaviour of alternating activated sludge treatment plants, with the aim of applying it to the optimal control of an aeration system. The reduction scheme was based on appropriate simplifications to the ASM 1 model (which is more appropriate for open-loop control). The objective was to verify if accurate predictions could be made time periods of several hours (about 8 h).

The present results are related to an existing small-size wastewater treatment plant. This plant was designed for 15,000 population-equivalents (p.e.) and consists of a primary treatment stage (screening, grit removal, primary sedimentation), followed by a secondary treatment stage (biological treatment). The latter consists of a single aeration tank of about 2,050 m³ equipped with 3 turbines which are operated cyclically to create alternating aerobic and anoxic conditions. Ammonia is converted into nitrate during air-on periods (nitrification step) and nitrate is subsequently removed during air-off periods (denitrification step). It is important to note that a dynamic model, based on the ASM 1 model and calibrated from a set of input/output measurements over a one-day period (Chachuat, 2001), was used here as a reference to perform model reduction.

The following two-level simplification procedure was applied :

· A homotopy method was first used to establish relationships between the states and the dynamics of the system, via an eigenvalue decomposition. The components that are associated with the slowest dynamics are then assumed constant to reduce the state space dimension. Heterotrophic (X_B,H) and autotrophic (X_B,A) biomass and inert particulate organic compounds (X_I) were detected as the slow state variables. It was found that the short-term predictions of the dynamic model were not affected by assuming that X_I, X_B,H and X_B,A concentrations were constant. Eliminating these 3 state variables, along with the concentrations of soluble inert organic compounds (S_I), resulted in a 7-dimensional dynamic model.

· However, further simplifications were required to enable the on-line optimisation of the bioreactor aeration profiles with reasonable computational times. These simplifications consisted of taking into account the process specifications in order to reduce the state space dimension to 4 or 5, and were therefore based on more heuristic considerations. Both organic and nitrogenous compounds are under consideration: (i) a single organic compound (denoted as X_DCO) is formed by adding soluble and particulate organic compound concentrations, and (ii) the mathematical expression that describes the organic nitrogen hydrolysis process is simplified so that the dynamics with respect to soluble and particulate organic nitrogen are independent.

The two previous simplification steps produced a reduced 5-dimensional dynamic model with state variables X_DCO, S_NO, S_NH, S_ND and S_O. It should also be noted that the resulting model involved the parameters Y_H, _iNBM, K_S, K_NO, K_O,H, K_NH,A, η_NO,g and η_NO,h that are identical to those defined in the original ASM 1 model by Henze et al. (1987). In addition, 7 specific parameters were defined defined (θ₁, θ₂, θ₃, θ₄, θ₅, K_DCO, K_ND). These new parameters exhibited rather slow temporal variation, thus agreeing with the general ASM 1 model for short time periods.

Afterwards, a two-step procedure was applied to calibrate the model. This procedure first consisted of determining a reduced set of identifiable parameters by the use of both sensitivity and principal component analyses. Note that the inlet concentrations of organic compounds, ammonia nitrogen and soluble organic nitrogen may be considered as additional parameters since they are generally not measured on-line. The selected parameters (θ₁, θ₂, θ₃) and inlet concentrations (Xⁱⁿ_DCO, Sⁱⁿ_NH) were then estimated by the application of a local gradient search method (successive quadratic programming, SQP). Comparisons between the dynamic behaviour of both reduced and ASM 1 models show that accurate predictions can be obtained over time periods of several hours (8 h). It was also shown that the reduced model was observable and structurally identifiable under aerobic and anoxic conditions from dissolved oxygen, ammonia and nitrate concentration measurements. These results therefore demonstrate the ability of the reduced model to be embedded into closed-loop control schemes.

The conclusions from this work are twofold: (i) The reduced model can be used as a basis to construct an on-line observer to estimate the unmeasured state variables, the unknown (most sensitive) parameters and inlet concentrations; (ii) Non-linear model predictive control (NMPC) schemes can then be implemented to operate the aeration system so that the nitrogen discharge or the energy consumption are minimised (optimal control).

The initial results demonstrate that the application of NMPC strategies is likely to give large reductions of nitrogen discharge with respect to usual operating strategies (e.g., oxygen or redox control). Such closed-loop control schemes are particularly efficient in dealing with large influent variations (inlet flow rate, concentration and composition) resulting from both human activities and climatic conditions, and inherent modelling uncertainties. However, an experimental validation of this control strategy on a pilot scale or an industrial scale is required to confirm these results.

Keywords:

• Wastewater treatment,
• alternating activated sludge process,
• nutrient removal,
• ASM 1 model,
• model reduction,
• optimal control