EP1756012A1 - Method for treating waters by means of a biological reactor in which the speed of continuously injectable air is adjustable and corresponding device - Google Patents

Abstract

The invention relates to a method for treating waters with the aid of a biological reactor for fighting against nitrogenous pollution which is called input load and contained in said waters, wherein the reactor comprises a biomass ventilate by injected air. The inventive method involves at least one injected air speed adjustment stage and consists in continuously measuring said input N-NH<sub

Description

Procédé de traitement d'eaux à l'aide d'un réacteur biologique, dans lequel la vitesse d'air injecté en continu dans le réacteur est régulée, et dispositif correspondant. La présente invention concerne le domaine du traitement de l'eau. Plus précisément, l'invention se rapporte aux traitements des effluents urbains ou industriels par un procédé d'élimination de la pollution azotée et carbonée grâce à une biomasse en culture libre ou fixée sur un matériau solide. Une technique classique de dépollution des eaux met en œuvre un réacteur biologique tel qu'un biofiltre ou une boue activée. Selon cette technique, le réacteur est aéré pour assurer le traitement de la pollution. L'art antérieur repose le plus souvent sur deux phases de traitement distinctes, l'une nécessitant de l'oxygène, la nitrification (N), et l'autre son absence, la dénitrification (DN). Ces deux phases du processus peuvent être réalisées soit dans un seul réacteur avec des périodes d'aération puis de non aération (aération séquencée), soit dans un réacteur muni de deux compartiments spécifiques (l'un aéré en continu, l'autre jamais). L'invention consiste à réaliser ces deux réactions antagonistes (Nitrification/Dénitrification) simultanément dans le même réacteur sans compartimentation spécifique dédiée à l'une ou l'autre des réactions et avec une aération permanente. Cependant, le niveau de cette aération doit être correctement maîtrisé car plus l'oxygène disponible dans le réacteur est en excès, plus la dénitrification est inhibée et vice et versa. Ainsi le rendement global du procédé de dépollution dépend directement de l'ajustement de l'aération. Dans les procédés de l'art antérieur, afin d'ajuster l'aération de la biomasse en fonction de la pollution à traiter, on effectue généralement des mesures de paramètres donnés dans le réacteur, ou en sortie de celui-ci. Ils ont pour but de déterminer l'état du système à un instant précis. A partir de cet état, une commande de débit de gaz est calculée puis appliquée. Il s'agit, le plus souvent d'une mesure d'ammonium, de nitrates, de potentiel d'oxydoréduction ou de l'oxygène dissous... D'autres paramètres, que l'on peut nommer « composites » sont aussi utilisés. Il s'agit de combinaisons linéaires des variables mesurées. Ainsi, dans le cadre de la conduite de procédé en traitement des eaux résiduaires urbaines, l'indicateur de charge polluante (masse de N-NH4 en kilogramme par mètre cube de réacteur aéré et par jour) a déjà été validé comme une donnée utile à la conduite. En effet, la demande de brevet publiée sous le numéro WO 01/02306 Al définit un estimateur de cette charge polluante à partir des mesures de conductivité et de turbidité. Cependant, il s'agit d'une estimation et non d'une mesure directe. Or dans le cas d'une régulation qui nécessite un degré de précision élevée, il n'est pas possible d'utiliser de tels outils. Une des justifications de l'utilisation de tels estimateurs tient au coût d'achat et d'entretien des analyseurs spécifiques. On comprend donc que ce faible coût est accompagné de performances très limitées. Actuellement, on distingue dans l'état de la technique différents types de stratégies de régulation. Selon l'une de ces stratégies de régulation, on calcule le rendement d'utilisation par la biomasse de l'oxygène distribué « Ct » pour réguler le débit de gaz. Ceci est décrit par le brevet publié sous le numéro FR2764 817. Cette technique nécessite l'emploi d'un Ct de référence qui lui-même dépend de la pollution à traiter. Or, la difficulté réside dans le fait que la charge polluante varie dans le temps et qu'il faut donc changer ce Ct de référence régulièrement afin qu'il soit toujours cohérent, ce qui en exclut l'application dans le cas d'effluent fortement variable. D'autres techniques utilisent des combinaisons linéaires des variables d'entrées du système (ammonium, nitrates...) pour calculer la commande en débit d'air à appliquer, tel que décrit notamment dans la demande internationale publiée sous le numéro WO 93/07089. Cependant, ces techniques s'appuient sur une modélisation empirique ou semi-empirique. Les fonctions de commande utilisées sont basées essentiellement sur les données issues de l'expérience des fonctionnements antérieurs. Si le modèle est désadapté du fait de l'évolution de la biomasse ou des performances du filtre, l'aération n'est alors plus optimale. L'art antérieur a également proposé des systèmes de contrôle utilisant une boucle de régulation dite « Feedback/Feedforward » (ce qui peut se traduire par boucle de contre réaction/boucle de prise en tendance). Selon une approche de ces systèmes, certains auteurs proposent des stratégies visant à piloter la concentration en oxygène dissous dans le réacteur biologique par l'utilisation d'une mesure d'ammonium en entrée. Selon une autre approche, le concept de charge en ammonium à traiter est utilisé pour prédire l'évolution nécessaire de Ia consigne en oxygène dissous. La présente invention n'utilise pas la mesure d'oxygène dissous qui représente uniquement l'oxygène disponible mais pas celui nécessaire à la biomasse. Selon encore d'autres techniques, la régulation du volume aéré s'appuie sur une compartimentation du réacteur (et non pas sur la régulation d'un débit ou d'une vitesse en gaz). Cependant, le caractère foncièrement discret (variation discontinue du volume aéré) de ce type de stratégie, incite les auteurs à toujours compléter leur système par une régulation de l'oxygène dissous. Quoiqu'il en soit, la majeure partie des publications s'appuie sur des résultats issus de simulation (Activated Sludge Model) ou d'essais batch avec des eaux synthétiques. En pratique, peu d'expériences à l'échelle pilote ou sur site ont été réalisées. En outre, les expériences conduites concernent quasi exclusivement des réacteurs biologiques à boues activées. L'invention a notamment pour objectif de pallier les inconvénients et/ou insuffisances de l'art antérieur. Plus précisément, l'invention a pour objectif de proposer un procédé de traitement des eaux mettant en œuvre un réacteur biologique, non compartimenté, intégrant une biomasse aérée par une injection d'air continue dans lequel la commande de l'injection d'air est optimisée par rapport aux techniques de l'art antérieur. En ce sens, l'invention a pour objectif de réaliser simultanément les étapes de nitrification et de dénitrification ainsi que l'abattement du carbone, dans une même enceinte et dans le même temps. L'invention a également pour objectif de fournir un tel procédé qui soit efficace en toute circonstance, y compris lorsque la charge polluante des eaux à traiter présente des variations importantes au fil du temps. L'invention a aussi pour objectif de fournir un tel procédé qui permette d'ajuster la commande d'air injecté en fonction des performances de la biomasse. Un autre objectif de l'invention est de fournir un tel procédé qui permette d'envisager des retours sur investissement relativement rapides et des gains en exploitation. En ce sens, l'invention a pour objectif de proposer un tel procédé qui permette d'éviter les surcoûts d'exploitation engendrés par des pics de charges polluantes journaliers. L'invention a encore pour objectif de fournir un tel procédé qui permette la commande d'air injecté sans qu'il soit nécessaire de procéder à des mesures en oxygène dissous. Encore un autre objectif de l'invention est de fournir un tel procédé qui soit facile à mettre en œuvre. Ces objectifs, ainsi que d'autres qui apparaîtront par la suite, sont atteints grâce à l'invention qui a pour objet un procédé de traitement d'eaux incluant les processus de nitrification et de dénitrification au sein d'une même enceinte d'un réacteur biologique en vue d'abattre la pollution azotée, dite charge entrante, contenue dans lesdites eaux, ledit réacteur intégrant une biomasse aérée par injection de gaz oxygéné, ledit procédé incluant au moins une étape de régulation de la vitesse dudit gaz injecté, caractérisé en ce que ladite injection de gaz est continue et en ce que lesdits processus de nitrification et de dénitrification sont essentiellement simultanés, ledit procédé comprenant une mesure en continu de ladite charge entrante en N-NH₄ (Cv_EDD) contenue dans lesdites eaux, ladite charge entrante mesurée étant pondérée par un décalage temporel pour en déduire une commande de la vitesse dudit gaz à injecter en fonction du temps. On note qu'on entend par injection continue de gaz oxygéné une injection permanente quoique variable. On note également que la charge entrante CV_E00 peut s'exprimer de la façon suivante : CVEDD = ^• 1000 x 1000 x V_** où : Q est le débit d'alimentation (exprimé ici en litre par heure), C_NH4 la concentration mesurée en entrée en ammonium (mg N-NH4 par litre) et V_aéré le volume du réacteur aéré (m³). Dans le cas où il existe plusieurs réacteurs identiques en parallèle le volume aéré est égal au produit du nombre de réacteur aéré en fonctionnement (nf) par la surface unitaire (S₁₁) et par la hauteur (hmat) unitaire de chaque réacteur soit V_aéIé =nf x S₁₁ x h_mat et S = nf x S_u (S surface totale du réacteur aéré). Ainsi, le procédé selon l'invention s'appuie sur une mesure directe de la charge comme paramètre d'entrée d'un modèle mathématique de prédiction des besoins en air (ou plus généralement des besoins en oxygène). De cette façon, comme cela va apparaître plus clairement par la suite, on obtient un procédé de traitement des eaux qui donne une commande plus efficace et plus précise que les procédés de l'art antérieur. Cette commande permet une injection de gaz continu mais variable dans le temps correspondant aux besoins de la biomasse et autorisant simultanément la nitrification et la dénitrification. Ce décalage temporel (qui peut être évolutif) par lequel on pondère la mesure de la charge entrante permet de prendre en compte notamment le décalage lié au transit du point de mesure au lieu de dégradation. Ce décalage lié au transit correspond sensiblement au décalage entre l'instant réel ou l'eau pénètre dans le réacteur et l'instant où l'ammonium commence à être dégradé. On tient compte également grâce à ce décalage du temps de passage dans le réacteur ou du temps de séjour moyen. Comme cela va apparaître plus clairement par la suite, un résultat particulièrement intéressant d'un tel procédé est qu'il permet de limiter les sur ou sous aérations lors des pics de charge polluante journaliers qui engendrent des surcoûts d'exploitation. Aussi, la présente invention, grâce à des gains de performances importants, permet un retour rapide sur investissement : de 1 à 2 ans. Par ailleurs, le procédé selon l'invention peut être appliqué aux biofiltres ainsi qu'à la plupart des procédés de traitement d'eaux, tels que notamment les boues activées y compris les bioréacteurs à membranes, les procédés à cultures fixées comme les biofiltres, les lits fluidisés, et les procédés à cultures mixtes. Selon une solution préférée, ladite ou lesdites étapes de régulation sont basées sur la loi suivante de commande de la vitesse de l'air en fonction du temps : Vair (t) = α(H(t) * Cv_EDD (t) - Cv_CONSIGNE) + β, dans laquelle Cv_CONSIGNE est une consigne de charge de sortie et H(t) est la transformée de Laplace inverse d'une fonction de retard de phase H(S). On note que le signe « * » indique un produit de convolution. Ainsi, on peut corriger en permanence la commande d'injection d'air en fonction des variations de la charge entrante dans la mesure réalisée en continu, ce qui permet de conserver une grande efficacité du traitement sur des effluents différents, et notamment sur des effluents à fortes variations de charge au cours du temps. Le raisonnement conduisant à l'établissement d'une telle loi de commande de prise en tendance est explicité de façon détaillée ci-après. L'objectif de la loi de commande est de relier la charge de N-NH4 entrante (en kg N-NH4/m³ de réacteur aéré/J) dans le réacteur biologique au débit d'air puis à la vitesse en air, la charge de N-NH4 entrante étant la charge réellement appliquée au réacteur, c'est-à-dire celle évaluée à partir du mélange de l'eau décantée et de l'eau recirculée. Cette variable permet de prendre en compte simultanément les variations de débit et de concentration. De plus, c'est un paramètre de dimensionnement des sites industriels. Or, on observe que la charge réelle éliminée est constante pour un débit d'air donné sur un horizon de 24 heures. Cette hypothèse est vérifiée par l'analyse de plusieurs essais en eau tamponnée à débit d'air constant. Lors d'expériences, on a fait varier la concentration en ammonium et le débit d'alimentation de sorte que la charge appliquée demeure constante. Ces expériences ont conduit à observer, après la phase transitoire, un retour à un niveau égal de la charge en N-NH4 éliminée avec pour seul paramètre invariant le débit d'air. Ceci est clairement attesté par le graphique de la figure 1 qui montre l'évolution au cours du temps de : - la charge entrante 11 ; - la charge éliminée 12 ; - la charge de sortie 13. Une relation linéaire entre la charge éliminée (Cve) en N-NH4 et la vitesse en air, en régime permanent a donc été établie : V_olr(ή= a -Cv_e(ή₊ β , avecProcess for treating water with the aid of a biological reactor, in which the air velocity continuously injected into the reactor is regulated, and corresponding device. The present invention relates to the field of water treatment. More precisely, the invention relates to the treatment of urban or industrial effluents by a process for eliminating nitrogen and carbon pollution by means of a biomass in free culture or fixed on a solid material. A conventional water depollution technique uses a biological reactor such as a biofilter or activated sludge. According to this technique, the reactor is ventilated to ensure the treatment of the pollution. The prior art is most often based on two distinct treatment phases, one requiring oxygen, nitrification (N), and the other its absence, denitrification (DN). These two phases of the process can be carried out either in a single reactor with periods of aeration and then without aeration (sequenced aeration), or in a reactor equipped with two specific compartments (one aerated continuously, the other never) . The invention consists in carrying out these two antagonistic reactions (nitrification / denitrification) simultaneously in the same reactor without specific compartmentation dedicated to one or other of the reactions and with permanent aeration. However, the level of this aeration must be properly controlled because the more oxygen available in the reactor is in excess, the more denitrification is inhibited and vice versa. Thus, the overall efficiency of the depollution process depends directly on the adjustment of the aeration. In the processes of the prior art, in order to adjust the aeration of the biomass as a function of the pollution to be treated, measurements of given parameters are generally carried out in the reactor, or at the outlet thereof. They aim to determine the state of the system at a precise moment. From this state, a gas flow control is calculated and then applied. It is most often a measure of ammonium, nitrates, oxidation-reduction potential or dissolved oxygen ... Other parameters, which we can call "composites" are also used. These are linear combinations of the measured variables. Thus, as part of the urban wastewater treatment process, the pollutant load indicator (mass of N-NH4 in kilograms per cubic meter of aerated reactor per day) has already been validated as a useful data source. the driving. Indeed, the patent application published under the number WO 01/02306 A1 defines an estimator of this pollutant load from the conductivity and turbidity measurements. However, it is an estimate and not a direct measure. But in the case of a regulation that requires a high degree of precision, it is not possible to use such tools. One of the justifications for using such estimators is the cost of purchasing and maintaining specific analyzers. We understand that this low cost is accompanied by very limited performance. Currently, there are different types of control strategies in the state of the art. According to one of these control strategies, the utilization efficiency by the "Ct" distributed oxygen biomass is calculated to regulate the gas flow rate. This is described by the patent published under the number FR2764 817. This technique requires the use of a reference Ct which itself depends on the pollution to be treated. However, the difficulty lies in the fact that the pollutant load varies over time and it is therefore necessary to change this reference Ct regularly so that it is always consistent, which excludes the application in the case of strongly effluent. variable. Other techniques use linear combinations of the input variables of the system (ammonium, nitrates, etc.) to calculate the airflow control to be applied, as described in particular in the international application published under the number WO 93 / 07089. However, these techniques rely on empirical or semi-empirical modeling. The control functions used are based essentially on the data resulting from the experience of previous operations. If the model is mismatched due to changes in biomass or filter performance, then aeration is no longer optimal. The prior art has also proposed control systems using a regulation loop called "Feedback / Feedforward" (which can be translated by feedback loop / loop tending). According to an approach of these systems, some authors propose strategies to control the concentration of dissolved oxygen in the biological reactor by the use of an input ammonium measurement. According to another approach, the concept of ammonium feedstock to be treated is used to predict the necessary evolution of the dissolved oxygen setpoint. The present invention does not use dissolved oxygen measurement which represents only available oxygen but not that required for biomass. According to still other techniques, the regulation of the aerated volume is based on compartmentalisation of the reactor (and not on the regulation of a flow rate or a gas velocity). However, the fundamentally discrete nature (discontinuous variation of the aerated volume) of this type of strategy, encourages the authors to always complete their system by a regulation of the dissolved oxygen. Nevertheless, the majority of publications are based on results from simulation (Activated Sludge Model) or batch tests with synthetic waters. In practice, few pilot-scale or on-site experiments have been conducted. In addition, the experiments conducted concern almost exclusively activated sludge biological reactors. The invention particularly aims to overcome the disadvantages and / or shortcomings of the prior art. More specifically, the invention aims to propose a water treatment method using a biological reactor, not compartmentalized, incorporating a biomass aerated by a continuous air injection in which the control of the air injection is optimized over the techniques of the prior art. In this sense, the objective of the invention is to simultaneously carry out the nitrification and denitrification steps as well as the carbon abatement, in the same enclosure and at the same time. The invention also aims to provide such a method that is effective in all circumstances, including when the polluting load of the water to be treated has significant variations over time. Another object of the invention is to provide such a method which makes it possible to adjust the injected air control as a function of the biomass performances. Another objective of the invention is to provide such a method which makes it possible to envisage relatively rapid returns on investment and gains in operation. In this sense, the object of the invention is to propose such a method which makes it possible to avoid the additional operating costs generated by daily pollutant load peaks. It is another object of the invention to provide such a method which allows the control of injected air without the need for dissolved oxygen measurements. Yet another object of the invention is to provide such a method which is easy to implement. These objectives, as well as others that will appear later, are achieved thanks to the invention which relates to a water treatment process including nitrification and denitrification processes within the same enclosure of a biological reactor for reducing the nitrogen pollution, said incoming charge, contained in said water, said reactor incorporating a biomass aerated by oxygenated gas injection, said method including at least one step of regulating the speed of said injected gas, characterized in said gas injection is continuous and in that said nitrification and denitrification processes are essentially simultaneous, said method comprising a continuous measurement of said incoming charge of N-NH ₄ (Cv _EDD ) contained in said waters, said measured incoming charge being weighted by a time shift to deduce a control of the speed of said gas to be injected as a function of time. It is noted that continuous injection of oxygen gas means a permanent injection, although variable. It is also noted that the incoming load CV _E00 can be expressed as follows: CVEDD = ^• 1000 x 1000 x V _** where: Q is the feed rate (expressed here in liters per hour), C _{NH4 is} the measured input concentration of ammonium (mg N-NH4 per liter) and V _aerated the volume of the aerated reactor (m ³ ). In the case where there are several identical reactors in parallel, the aerated volume is equal to the product of the number of aerated reactor operating (nf) by the unit area (S ₁₁ ) and the unit height (hmat) of each reactor is V _aéIé = nf x S ₁₁ xh _mat and S = nf x S _u (S total area of the aerated reactor). Thus, the method according to the invention is based on a direct measurement of the load as an input parameter of a mathematical model for predicting air requirements (or more generally oxygen requirements). In this way, as will become clearer later, a water treatment process is obtained which gives a more efficient and precise control than the prior art processes. This control allows a continuous but variable injection of gas corresponding to the needs of the biomass and simultaneously allowing nitrification and denitrification. This temporal shift (which can be scalable) by which the measurement of the incoming charge is weighted makes it possible to take into account in particular the offset linked to the transit of the measurement point instead of degradation. This shift related to the transit substantially corresponds to the shift between the actual moment when the water enters the reactor and the moment when the ammonium begins to be degraded. This shift in the reactor passage time or the average residence time is also taken into account. As will become clearer later, a particularly interesting result of such a method is that it allows to limit over or under aeration during daily pollutant load peaks that generate additional operating costs. Also, the present invention, thanks to significant performance gains, allows a quick return on investment: from 1 to 2 years. Furthermore, the process according to the invention can be applied to biofilters as well as to most water treatment processes, such as in particular activated sludge including membrane bioreactors, fixed culture processes such as biofilters, fluidized beds, and mixed-culture processes. According to a preferred solution, said one or more regulation steps are based on the following law of control of the speed of the air as a function of time: Vair (t) = α (H (t) * Cv _EDD (t) - Cv _SETPOINT ) + β, in which Cv _SETPOINT is an output load setpoint and H (t) is the inverse Laplace transform of a H (S) phase delay function. Note that the sign "*" indicates a product of convolution. Thus, it is possible to permanently correct the air injection control as a function of the variations in the incoming charge in the measurement carried out continuously, which makes it possible to maintain a high efficiency of the treatment on different effluents, and especially on effluents. with strong load variations over time. The reasoning leading to the establishment of such a trend control law is explained in detail below. The objective of the control law is to connect the feedstock of incoming N-NH4 (in kg N-NH4 / m ³ of aerated reactor / J) in the biological reactor to the air flow rate and then to the air velocity. incoming charge of N-NH4 being the charge actually applied to the reactor, that is to say the one evaluated from the mixture of the settled water and recirculated water. This variable makes it possible to simultaneously take into account variations in flow and concentration. In addition, it is a sizing parameter for industrial sites. However, it is observed that the actual charge eliminated is constant for a given air flow over a 24-hour horizon. This assumption is verified by the analysis of several constant air flow buffered water tests. In experiments, the ammonium concentration and the feed rate were varied so that the applied load remained constant. These experiments led to observe, after the transient phase, a return to an equal level of the charge in N-NH4 eliminated with only parameter invariant the air flow. This is clearly attested by the graph of FIG. 1 which shows the evolution over time of: the incoming charge 11; - the eliminated charge 12; - The output charge 13. A linear relationship between the eliminated charge (Cve) in N-NH4 and the steady-state air velocity has therefore been established: V _olr (ή = a -Cv _e (ή ₊ β, with

^{a>r =} nf xS_u ^{a> r =} nf xS _u

Cependant la relation obtenue ne prend pas encore en compte les phénomènes transitoires. Une des particularités de l'invention réside donc dans la prise en compte des phénomènes transitoires pour l'obtention des paramètres de la relation précédente. En effet, pour relier la charge éliminée à la vitesse en air, il faut tenir compte du temps de passage dans le réacteur ou du temps de séjour moyen Tg (temps de transit moyen d'une particule fluide dans un réacteur considéré). Le paramètre Tg est obtenu pratiquement en étudiant le coefficient de cross- corrélation reliant la concentration d'eau décantée diluée (eau décantée + eau recirculée) à la mesure de sortie du procédé. Or, d'après les observations précédentes, à débit d'air fixé la charge éliminée est constante, et si l'on fixe le débit d'alimentation, cela revient à considérer un abattement (différence des concentrations Entrée/Sortie) en N-NH4 constant. Ainsi, la concentration d'ammonium de sortie s'exprime alors en fonction de la concentration d'entrée diluée :However, the relation obtained does not yet take into account the transient phenomena. One of the peculiarities of the invention therefore lies in taking transient phenomena into account in order to obtain the parameters of the preceding relation. Indeed, to connect the removed charge to the air velocity, it is necessary to take into account the passage time in the reactor or the average residence time Tg (average transit time of a fluid particle in a reactor considered). The Tg parameter is obtained practically by studying the cross-correlation coefficient connecting the diluted decanted water concentration (decanted water + recirculated water) to the output measurement of the process. However, according to the above observations, at a fixed air flow rate, the eliminated charge is constant, and if the feed rate is set, this amounts to considering an abatement (difference in input / output concentrations) in N- Constant NH4. Thus, the output ammonium concentration is then expressed as a function of the diluted inlet concentration:

[NH4]_s(t) s[NH4]_e(t-T_g)-Cst où Cst est une constante que l'on va définir par exemple par la mesure. De plus, le coefficient de cross-corrélation est maximum pour le pas n correspondant au décalage temporel Tg≈nTe (où Te est la période d'échantillonnage). Par extension, la charge éliminée est calculée en faisant la différence entre la charge d'entrée et la charge de sortie à laquelle on applique un retard égale à Tg. On peut dès lors exprimer la vitesse en air à l'instant t en fonction de la charge éliminée calculée au même instant (voir figure 2, courbe F(t)). Néanmoins, cette méthode ne tient pas compte de l'instant réel où l'eau décantée diluée pénètre dans le réacteur et où l'ammonium commence à être dégradé, décalage lié au transit du point de mesure au lieu de dégradation. De même, elle ne tient pas compte de la variation d'aération qui n'a probablement pas un effet immédiat et de l'hydraulique du système. L'ensemble de ces phénomènes induit une dispersion des valeurs (comme le montre la figure 2 dans le cas d'un biofiltre) et donc une incertitude sur la vitesse en air à appliquer jusqu'à ±15% autour de 7 Nm/h. Une autre particularité de l'invention tient donc à l'application d'une fonction de retard de phase H(s) qui permet de corriger ces omissions et d'obtenir la linéarité attendue (comme le montre la courbe H(t)*F(t) sur la figure 2). Ainsi, la charge éliminée au cours du temps (Cv_e) est égale à la différence de charge entre l'entrée (Cv_EDD) convolué par H(t) (transformé de Laplace inverse d'une fonction de retard de phase H(s)) et la sortie (Cv₈). Cette méthodologie permet d'obtenir la relation linéaire entre la vitesse en air et la charge éliminée calculée précédemment. On obtient donc la loi de commande : Vair (t) = α(H(t) * Cv₁₃₀₀ (t) - Cv_C0NSIGNE) + β Par ailleurs, il est connu que l'utilisation de modèles avec des paramètres empiriques présente des inconvénients en terme de déploiement sur sites. En effet, l'obtention sur une usine des valeurs nécessaires aux réglages n'est pas toujours compatible avec les priorités d'exploitation. La présente invention offre donc la possibilité d'obtenir des préréglages de manière théorique. Ceci a été mis en évidence par une réification du modèle, consistant à rapprocher la loi empirique observée d'une loi théorique afin de s'affranchir des problèmes de dimensionnement et de transfert d'échelle. L'explication de la loi empirique obtenue vient du calcul des besoins théoriques en oxygène correspondant à la nitrification et à l'élimination du carbone. Les besoins en 02 sont déterminés par les relations connues suivantes : Nitrification : B.O2NH4 = 4,57 x Masse de N-NH4 éliminée Dénitrification : B.O2NO3 = -2.86 x (Masse de N-NH4 éliminée-Masse de NO3 produite) Carbone : B.O2DCOs = 0,90 x Masse de DCOs éliminée Le besoin en oxygène pour le réacteur biologique s'exprime donc comme la somme des besoins précédents. Les besoins en air découlent de ces besoins en oxygène. En effet, une relation les lie avec le rendement de transfert, dont la valeur diminue avec la vitesse en air. La figure 3 est un graphique permettant de comparer les demandes empirique 32 et théorique 31 en air pour un biofiltre. Ce graphique montre que les besoins théoriques en air ainsi calculés concordent bien avec la loi empirique. Il existe donc une base théorique pour l'obtention des paramètres recherchés sur site On précise que les coefficients des deux régressions linéaires correspondent aux valeurs de α et β (α est égale à 25.7 par la théorie et 24.1 par la mesure, β est égale à -4.2 par la théorie et -3.6 par la mesure), R² représentant le coefficient de confiance sur 100% associé à chaque régression linéaire (les coefficients obtenus d'un côté par la mesure et de l'autre par la théorie sont quasiment identiques). Préférentiellement, on exprime la loi de commande en boucle ouverte, de la façon suivante : Vair (t) = α(H(t) * Cv_EDD (t) - Cv_C0NSIGNE(t + Δt)) + β, De cette façon, le terme Cvconsigne(t+Δt) permet d'anticiper une variation future de la consigne de la charge de sortie. Cette solution peut offrir une souplesse supplémentaire non négligeable pour l'optimisation de l'exploitation d'une station. Selon un mode de réalisation préférentiel, ladite fonction de retard de phase est du type :[NH4] _s (t) s [NH4] _e (tT _g ) -Cst where Cst is a constant that will be defined for example by the measurement. In addition, the cross-correlation coefficient is maximum for step n corresponding to the time offset Tg≈nTe (where Te is the sampling period). By extension, the eliminated charge is calculated by differentiating between the input load and the output load to which a delay equal to Tg is applied. It is then possible to express the air speed at time t as a function of the eliminated charge calculated at the same time (see Figure 2, curve F (t)). However, this method does not take into account the actual moment when the diluted decanted water enters the reactor and where the ammonium begins to be degraded, offset due to the transit of the measurement point instead of degradation. Likewise, it does not take into account the aeration variation that probably does not have an immediate effect and the hydraulics of the system. All of these phenomena induce a dispersion of the values (as shown in Figure 2 in the case of a biofilter) and therefore an uncertainty on the air velocity to be applied up to ± 15% around 7 Nm / h. Another peculiarity of the invention is therefore the application of a phase delay function H (s) which makes it possible to correct these omissions and to obtain the expected linearity (as shown by the curve H (t) * F (t) in Figure 2). Thus, the charge eliminated over time (Cv _e ) is equal to the charge difference between the input (Cv _EDD ) convoluted by H (t) (inverse Laplace transform of a phase delay function H (s) )) and the output (Cv ₈ ). This methodology makes it possible to obtain the linear relationship between the air speed and the eliminated load calculated previously. We thus obtain the control law: Vair (t) = α (H (t) * Cv ₁₃₀₀ (t) - Cv _C0NSIGN ) + β Furthermore, it is known that the use of models with empirical parameters has disadvantages in terms of deployment on sites. In fact, obtaining on a plant the values necessary for the adjustments is not always compatible with the operating priorities. The present invention thus offers the possibility of obtaining presets in a theoretical manner. This has been demonstrated by a reification of the model, consisting of bringing the empirical law observed closer to a theoretical law in order to overcome the problems of sizing and scaling. The explanation of the empirical law obtained comes from the calculation of the theoretical oxygen requirements corresponding to the nitrification and the elimination of the carbon. The O 2 requirements are determined by the following known relationships: Nitrification: B.O2NH4 = 4.57 x Mass of N-NH4 removed Denitrification: B.O2NO3 = -2.86 x (Mass of N-NH4 removed-Mass of NO3 produced) Carbon: B.O2DCOs = 0.90 x Mass of CODs eliminated The oxygen requirement for the biological reactor is therefore expressed as the sum of the previous needs. Air requirements arise from these oxygen requirements. Indeed, a relationship links them with the transfer efficiency, whose value decreases with air speed. Figure 3 is a graph for comparing empirical 32 and theoretical 31 air demand for a biofilter. This graph shows that the theoretical air requirements thus calculated agree well with the empirical law. There is therefore a theoretical basis for obtaining the parameters sought on site It is specified that the coefficients of the two linear regressions correspond to the values of α and β (α is equal to 25.7 by theory and 24.1 by measurement, β is equal to -4.2 by the theory and -3.6 by the measure), R ² representing the 100% confidence coefficient associated with each linear regression (the coefficients obtained on one side by the measurement and on the other by the theory are almost identical ). Preferably, the open loop control law is expressed in the following way: Vair (t) = α (H (t) * Cv _EDD (t) - Cv _C0NSIGN (t + Δt)) + β, In this way, the term Cv (t + Δt) makes it possible to anticipate a future variation of the setpoint of the output load. This solution can offer significant additional flexibility for optimizing the operation of a station. According to a preferred embodiment, said phase delay function is of the type:

ou : - n est un paramètre d'ajustement de la diffusion dans ledit ou lesdits réacteurs ; - V est le volume apparent dudit ou desdits réacteurs biologiques ; - Q est le débit d'alimentation desdites eaux à traiter ; - s est la transformée de Laplace de la variable t. Cette fonction correspond à la fonction de transfert dans l'espace de Laplace d'une série de n réacteurs parfaitement agités (RPA), n étant un nombre entier positif et non nul. Cette fonction de retard de phase revêt ici la forme mathématique d'une série de réacteurs parfaitement agités (RPA) dont les paramètres sont le débit et le volume. Le débit est asservi à celui de filtration de la colonne car seul le volume apparent V est utilisé pour ajuster la fonction. Plus généralement, toutes les expressions mathématiques utilisant des combinaisons de fonction de transfert de réacteurs idéaux peuvent être utilisées. La fonction de retard de phase peut donc être de différents types selon d'autres modes de réalisation envisageables. Par exemple, cette fonction peut être exprimée sous la forme d'une fonction de transfert représentant une série de quatre RPA, avec n=4. Selon une solution avantageuse, le procédé comprend au moins une étape de mesure d'une charge de sortie (Cv₅) et/ou d'une concentration en ammonium dissous contenues dans lesdites eaux traitées. Dans ce cas, ladite mesure de ladite charge de sortie (Cv₅) est préférentiellement exploitée pour ajuster ladite consigne par application de l'équation suivante : Cv_CONSi_GNE(t) - Cv₈= erreur(t), , . 24 x Qx e(t) ou erreuήt) = - ±¹ r, ^W (1000 x 1000 x A₁^ x 5) le signal e(t) étant défini par :, e(t) = [NH₄]_CONSlGNE(t)-[NH₄]_s(t) h_mat étant la hauteur de ladite biomasse dans ledit réacteur et S étant la surface dudit ou desdits réacteurs, le volume résultant du produit h_mat par S représentant le volume aéré dudit réacteur biologique. Cette formule est issue du calcul de charge, e(t) dans cette formule représentant la différence de concentration entre la consigne [NH4]_consîgne(t) et la mesure de sortie [NH4]_s(t). La présente invention, par l'intermédiaire d'un tel système de rétroaction, apporte une correction continue au modèle de prédiction de sorte à toujours être en adéquation avec la réalité. Ainsi, on peut compenser les erreurs de modèle de prise en tendance ou les perturbations non-mesurables, ceci en adjoignant un terme de rétroaction. En effet, une boucle fermée avec prise en tendance permet d'obtenir de bien meilleurs résultats qu'un régulateur « FeedForward » ou « FeedBack » seul. Une des spécificités de la rétroaction réside donc dans l'évaluation de l'erreur. En effet, il ne s'agit pas uniquement d'une différence entre la mesure de l'ammonium dissous en sortie et de la consigne (noté e(t)), mais d'une différence entre une charge de sortie et une charge de consigne (noté erreur(t)). Bien qu'il existe une relation entre e(t) et erreur(t), la différence entre ces deux paramètres est fondamentale. Selon un premier mode de réalisation, ladite mesure d'une charge de sortie (Cv_s) et/ou d'une concentration en ammonium dissous contenues dans lesdites eaux traitées, est réalisée en continu. Le régulateur prend implicitement en compte les variations du débit d'alimentation pour le calcul de l'action à appliquer. Selon un mode de réalisation préférentiel, le procédé comprend une étape de transformation dudit signal e(t) réalisée de telle sorte que : - si e > 0, f(e) = e ; - si e ≤ 0, f(e) = 1 - exp(-k.e), où k > 0. On note que k permet de graduer l'intensité de la transformation du signal. On apporte ainsi une amélioration grâce à une transformée du signal e(t). En effet, la concentration en ammonium de sortie n'est pas une fonction strictement linéaire du débit de gaz puisqu'elle est bornée par zéro. Par exemple, lorsque la consigne est de 2 mg-N-NH4/L, l'erreur mesure/consigne ne pourra prendre que des valeurs comprises dans l'intervalle [-2 , +∞ [, dans ce cas, il y a une dissymétrie néfaste de l'intervalle de régulation. Une telle transformée permet donc de prendre en compte cette dissymétrie. Cette transformation permet d'amplifier la réaction du régulateur lorsqu'il s'approche de la limite inférieure. Dans le cas de la consigne égale à 2mg-N-NH4/L, l'erreur prend alors des valeurs comprises dans l'intervalle [-6.4 , +∞ [ si k=l. Ainsi, au prix d'une faible déstabilisation du système, entraînée par l'augmentation du gain dans ce régime de fonctionnement, on limite l'effet de la saturation. On peut imaginer d'autres transformations telles que f(x)=x², qui permettrait d'avoir une variation de commande douce lorsque l'on est proche de la valeur de consigne et de plus en plus forte à mesure que l'on s'en éloigne. Enfin, l'utilisation de boucle de régulation Feedback et Feedforward offre une sécurité accrue car la commande est une composition de plusieurs signaux. Lorsqu'une information est absente, les autres permettent tout de même une régulation même si elle est moins performante. L'invention concerne également un dispositif pour la mise en œuvre du procédé de traitement d'eaux incluant des phase de nitrification et de dénitrification au sein d'une même enceinte d'un réacteur biologique tel que décrit précédemment, ledit réacteur intégrant une biomasse aérée par injection d'air et des moyens de régulation de la vitesse dudit air injecté, caractérisé en ce qu'il comprend : - des moyens de mesure en continu de ladite charge entrante ; - des moyens de paramétrage d'une consigne de charge de sortie et/ou d'une concentration de sortie (Cv_C0NSIGNE) ; - des moyens de calcul destinés à agir sur lesdits moyens de régulation sur la base d'une loi de commande dans laquelle ladite charge entrante mesurée est notamment pondérée par un décalage temporel pour en déduire une commande de la vitesse dudit air à injecter en fonction du temps, ladite injection d'air étant continue et lesdits processus de nitrification et de dénitrification étant essentiellement simultanés. De tels moyens de calcul pourront comprendre la mise en œuvre de moyens logiciels. En outre, on peut noter que la présente invention n'utilise pas de mesure ou de consigne en oxygène dissous pour atteindre ces objectifs. Selon un mode de réalisation préféré, le dispositif comprend une boucle de rétroaction comprenant des moyens de mesure d'une charge de sortie (Cv_s) et/ou d'une concentration en ammonium dissous contenues dans lesdites eaux traitées, et des moyens de comparaison de ladite charge de sortie avec ladite consigne de charge de sortie. Dans ce cas, lesdits moyens de comparaison sont préférentiellement reliés auxdits moyens de calcul, en vue d'ajuster ladite consigne prise en compte dans ladite loi de commande. D'autres caractéristiques et avantages de l'invention apparaîtront plus clairement à la lecture de la description suivante d'un mode de réalisation préférentiel de l'invention, donné à titre d'exemple illustratif et non limitatif, et des dessins annexés parmi lesquels : - la figure 1 est un graphique montrant les variations de la concentration en ammonium à charge constante, en eau tamponnée ; - la figure 2 montre deux relevés de vitesses d'air en fonction de la charge éliminée calculée pour un biofiltre ; - la figure 3 est un graphique montrant une comparaison des demandes empirique et théorique en air pour un biofiltre ; - la figure 4 est une vue schématique d'un unité pilote selon l'invention ; - la figure 5 est une représentation schématique d'un bloc de régulation «Feedback/Feedformward » selon l'invention ; - la figure 6 montre les courbes de relevés obtenus avec le procédé selon l'invention, à débit d'alimentation et taux de recirculation fixes ; - la figure 7 montre les courbes de relevés obtenus avec le procédé selon l'invention, à débit d'alimentation variable et taux de recirculation fixes ; - la figure 8 montre les courbes de relevés obtenus avec le procédé selon l'invention, à débit d'alimentation et taux de recirculation variables. L'efficacité du procédé selon l'invention est démontrée ci-après par l'intermédiaire d'essais conduits à l'aide d'une unité pilote telle que représentée sur la figure 4. Tel que cela apparaît, cette unité est composée de deux colonnes de plexiglas d'une hauteur de 5 m pour un diamètre intérieur de 29 cm. La hauteur de matériau (h mat) utilisé pour les essais en Nitrification/Dénitrification (NDN) simultanée est de 2,75 m et le diamètre moyen des billes qui le constitue est de 3,34 mm ± 0.19 mm. L'effluent traité est issu d'un réseau d'eau urbaine ; il subit une décantation primaire (décanteur lamellaire) avant d'être acheminé gravitairement dans un bac tampon 42 de 30 litres agité en permanence. L'influent est ensuite relevé par l'intermédiaire de deux pompes SEEPEX afin de remplir les deux colonnes de mise en charge. Celles-ci offrent une perte de charge maximum disponible de 2,40 mCe. Une fraction 43 de l'eau traitée est réutilisée dans le cadre de la nitrification / dénitrification pour alimenter le pilote. Cette eau est mélangée avec l'eau d'alimentation dans les colonnes de mise en charge. Là aussi, deux pompes Seepex permettent de recirculer les débits souhaités. L'eau traitée non recirculée 44 issue des deux colonnes se mélange dans un bac commun de 10 litres d'où sont prélevés les échantillons à analyser. Dans cette configuration, le pilote fonctionne à l'image d'un site comprenant deux cellules de filtration. On homogénéise ainsi le traitement et on dispose surtout d'une redondance des mesures en ligne. Deux rampes d'air (non représentées) situées à 20 cm du fond de chaque colonne permettent d'injecter l'air procédé de manière continue mais variable dans tout le réacteur, et deux autres orifices situés au fond de chacune d'elles permettent l'injection de l'air lavage. Dans les deux cas, la production d'air est assurée par un réseau d'air comprimé. La charge appliquée durant les essais, calculée sur la totalité du matériau, est comprise en 0,3 et 0,6 kg N-NH4/m3/J pour une moyenne de 0,45 kg N- NH4/m3/J. La vitesse en eau d'alimentation V₆₃₁₁ moyenne correspondante est de 1,2 m/h pour un taux de recirculation moyen de 125%. Ce taux diffère de la NDN classique (c'est-à-dire avec un réacteur avec une zone aérée et une zone non aérée) puisque les rendements mis en jeu sont différents. Trois cas d'études illustrent les performances de l'invention : le premier (figure 6) démontre l'efficacité de la régulation à débit et taux de recirculation fixes. Le second (figure 7) illustre les performances à débit variable et taux de recirculation fixe. Et le dernier (figure 8), montre les résultats obtenus à débit et taux de recirculation variable. La figure 6 illustre donc qu'il est possible d'éliminer le pic de charge journalier par une anticipation des besoins en air. Il n'existe pas de zone de sous ou sur-aération, avant et après Ie pic de charge, contrairement à l'état de l'art précédent. La dénitrification n'est donc pas pénalisée pour atteindre l'objectif de nitrification puisque la commande d'aération permet de maintenir une concentration d'ammonium identique tout au long de la journée. On constate que la régulation objet de l'invention fonctionne correctement et anticipe le pic de charge, dans les conditions de débit d'alimentation variable (figure 7). Excepté pendant une courte période où la charge est telle que la saturation en vitesse en air a été atteinte. En effet, l'aération a été bridée pour que celle-ci ne dépasse pas les 15 Nm/h afin de rester représentatif de la majeure partie des installations. Enfin, dans le cas d'une utilisation avec des débits d'alimentation et taux de recirculation variables (figure 8), ce qui n'est pas le cas représentatif des installations industrielles, les résultats, bien que dégradé, restent nettement supérieurs à l'état de l'art antérieur. Ces résultats sont obtenus à l'aide d'un bloc de régulation tel que celui illustré par la figure 5 qui correspond à un mode de réalisation préférentiel du procédé selon l'invention. Tel que représenté, une eau à traiter 51 est dirigée vers un réacteur biologique 52, les données de débit d'alimentation Q et de charge entrante étant associées à cette eau à traiter. On rappelle que, selon le principe de l'invention, les processus de nitrification et de dénitrification sont réalisés au sein d'une même enceinte du réacteur (le réacteur présentant en l'occurrence un compartiment unique) et que ces processus de nitrification et de dénitrification sont essentiellement simultanés. Des mesures sont effectuées concernant cette eau à traiter à l'aide d'un premier régulateur 53, dit régulateur « FeedForward » qui effectue notamment une mesure en continu de la charge entre C_VEDD. Le régulateur reçoit de plus une information concernant la consigne de charge de sortie Cv_C0NSIGNE. Un autre régulateur 54, dit « Feedback », permet de recueillir des informations, et notamment la charge de sortie C_vs mesurée en continu. Le régulateur 54 reçoit également l'information concernant la consigne de charge de sortie Cv_C0NSIGNE. Le régulateur utilisé pour la boucle de rétroaction peut être du type PID (régulateur par action proportionnelle) ou PFC (Prédictive Functional Control). Ces régulateurs sont paramétrés pour émettre chacun une commande. Les commandes correspondantes sont traitées de façon à agir sur la vitesse de l'air V_3J1., en appliquant la loi de commande de la vitesse de l'air en fonction du temps : Vair(t)= a x (H(O* Cv_EDD (ή- Cv_Consigne (*)> β + Vair_Mtmaction Les coefficients α et β peuvent être obtenus par des mesures ou par calcul, tel qu'indiqué précédemment. Vair rétroaction correspond à la variation de débit d'air nécessaire calculée par le régulateur de la boucle de rétroaction. De plus, dans le cas de l'unité pilote décrite précédemment, la fonction de retard de phase H(s) s'exprime de la façon suivante : or: n is an adjustment parameter of the diffusion in said reactor or said reactors; V is the apparent volume of said one or more biological reactors; - Q is the feed rate of said water to be treated; - s is the Laplace transform of the variable t. This function corresponds to the transfer function in the Laplace space of a series of n perfectly stirred reactors (RPA), n being a positive integer and not zero. This phase delay function here takes the mathematical form of a series of perfectly stirred reactors (RPA) whose parameters are flow and volume. The flow rate is slaved to that of filtration of the column because only the apparent volume V is used to adjust the function. More generally, all mathematical expressions using ideal reactor transfer function combinations can be used. The phase delay function can therefore be of different types according to other possible embodiments. For example, this function can be expressed as a transfer function representing a series of four RPAs, with n = 4. According to an advantageous solution, the process comprises at least one step of measuring an outlet charge (CV ₅ ) and / or a dissolved ammonium concentration contained in said treated water. In this case, said measurement of said output load (Cv ₅ ) is preferably used to adjust said setpoint by applying the following equation: Cv _CONS i _GNE (t) - Cv ₈ = error (t),,. 24 x Qx e (t) or erreuήt) = - ± ¹ r, ^W (1000 x 1000 x A ₁ ^ x 5) the signal e (t) being defined by _:, e (t) = [NH ₄ ] _CONSlGNE ( t) - [NH ₄ ] _s (t) h _mat being the height of said biomass in said reactor and S being the surface of said reactor (s), the resultant volume of the product _mat by S representing the aerated volume of said biological reactor. This formula is derived from the design load, e (t) in the formula representing the difference in concentration between the set [NH4] _Set (t) and the measurement output [NH4] _s (t). The present invention, through such a feedback system, provides a continuous correction to the prediction model so as to always be in adequacy with reality. Thus, one can compensate for errors in the trend model or non-measurable disturbances by adding a feedback term. In fact, a closed loop with a tendency to take a trend makes it possible to obtain much better results than a "FeedForward" or "FeedBack" regulator alone. One of the specificities of feedback is therefore the evaluation of the error. Indeed, it is not only a difference between the measurement of dissolved ammonium output and the setpoint (denoted e (t)), but of a difference between an output charge and a charge of setpoint (noted error (t)). Although there is a relationship between e (t) and error (t), the difference between these two parameters is fundamental. According to a first embodiment, said measurement of an outlet charge (Cv _s ) and / or a dissolved ammonium concentration contained in said treated water, is carried out continuously. The controller implicitly takes into account the variations of the feed rate for the calculation of the action to be applied. According to a preferred embodiment, the method comprises a step of transforming said signal e (t) carried out such that: - if e> 0, f (e) = e; if e ≤ 0, f (e) = 1 - exp (-ke), where k> 0. Note that k makes it possible to graduate the intensity of the transformation of the signal. This provides an improvement through a signal transformation e (t). Indeed, the output ammonium concentration is not a strictly linear function of the gas flow since it is bounded by zero. For example, when the setpoint is 2 mg-N-NH4 / L, the measurement / setpoint error can not be take only values in the interval [-2, + ∞ [, in this case, there is a dissymmetry harmful to the control interval. Such a transform therefore makes it possible to take this asymmetry into account. This transformation amplifies the regulator's reaction as it approaches the lower limit. In the case of the setpoint equal to 2mg-N-NH4 / L, the error then takes values in the range [-6.4, + ∞ [if k = l. Thus, at the price of a low destabilization of the system, driven by the increase in gain in this operating regime, the effect of saturation is limited. We can imagine other transformations such that f (x) = x ² , which would have a soft control variation when we are close to the setpoint and stronger and stronger as we away from it. Finally, the use of Feedback and Feedforward control loop offers increased security because the control is a composition of several signals. When information is absent, the others still allow regulation even if it is less efficient. The invention also relates to a device for implementing the water treatment method including nitrification and denitrification phase within the same enclosure of a biological reactor as described above, said reactor incorporating an aerated biomass by air injection and means for regulating the speed of said injected air, characterized in that it comprises: means for continuously measuring said incoming charge; means for setting an output load setpoint and / or an output concentration (Cv _C0NSIGN ); calculation means intended to act on said regulation means on the basis of a control law in which said measured incoming charge is weighted in particular by a time shift for deriving a control of the speed of said air to be injected as a function of time, said air injection being continuous and said nitrification and denitrification processes being essentially simultaneous. Such calculation means may include the implementation of software means. In addition, it may be noted that the present invention does not use a dissolved oxygen measurement or setpoint to achieve these objectives. According to a preferred embodiment, the device comprises a feedback loop comprising means for measuring an exit charge (Cv _s ) and / or a dissolved ammonium concentration contained in said treated water, and means of comparison. of said output load with said output load setpoint. In this case, said comparison means are preferably connected to said calculation means, in order to adjust said setpoint taken into account in said control law. Other characteristics and advantages of the invention will emerge more clearly on reading the following description of a preferred embodiment of the invention, given by way of illustrative and nonlimiting example, and the appended drawings among which: FIG. 1 is a graph showing the variations of the ammonium concentration at constant load, in buffered water; FIG. 2 shows two air velocity readings as a function of the eliminated charge calculated for a biofilter; FIG. 3 is a graph showing a comparison of empirical and theoretical air demands for a biofilter; FIG. 4 is a schematic view of a pilot unit according to the invention; FIG. 5 is a schematic representation of a "Feedback / Feedformward" control block according to the invention; FIG. 6 shows the curves of readings obtained with the method according to the invention, with feed rate and fixed recirculation rates; FIG. 7 shows the curves of readings obtained with the method according to the invention, with variable feed rate and fixed recirculation rates; - Figure 8 shows the survey curves obtained with the method according to the invention, variable feed rate and recirculation rate. The efficiency of the method according to the invention is demonstrated below by means of tests carried out using a pilot unit as shown in FIG. 4. As it appears, this unit is composed of two Plexiglas columns with a height of 5 m and an inside diameter of 29 cm. The material height (h matt) used for the simultaneous Nitrification / Denitrification (NDN) tests is 2.75 m and the average diameter of the balls which constitute it is 3.34 mm ± 0.19 mm. The treated effluent comes from an urban water network; it undergoes a primary settling (lamellar decanter) before being conveyed gravitarily in a buffer tank 42 of 30 liters stirred continuously. The influent is then raised via two SEEPEX pumps to fill the two loading columns. These offer a maximum available pressure drop of 2.40 mCe. A fraction 43 of the treated water is reused as part of the nitrification / denitrification to feed the pilot. This water is mixed with the feed water in the loading columns. Here again, two Seepex pumps make it possible to recirculate the desired flow rates. The non-recirculated treated water 44 from the two columns is mixed in a common tank of 10 liters from which the samples to be analyzed are taken. In this configuration, the driver works like a site with two filtration cells. This homogenizes the processing and is mainly redundant online measures. Two air ramps (not shown) located 20 cm from the bottom of each column can inject the process air continuously but variable throughout the reactor, and two other holes located at the bottom of each of them allow injection of air wash. In both cases, the production of air is provided by a compressed air network. The load applied during the tests, calculated on the totality of the material, is between 0.3 and 0.6 kg N-NH4 / m3 / J for an average of 0.45 kg N-NH4 / m3 / J. The corresponding mean feed velocity V ₆₃₁₁ is 1.2 m / h for an average recirculation rate of 125%. This rate differs from the conventional NDN (that is to say with a reactor with an aerated zone and a non-aerated zone) since the yields involved are different. Three cases of studies illustrate the performances of the invention: the first one (FIG. 6) demonstrates the efficiency of the regulation with rate and fixed recirculation rates. The second (Figure 7) illustrates the performance at variable rate and fixed recirculation rate. And the last (Figure 8), shows the results obtained at variable flow rate and recirculation. FIG. 6 thus illustrates that it is possible to eliminate the peak of daily load by anticipating the air requirements. There is no under or over-aeration zone, before and after the load peak, contrary to the state of the art above. Denitrification is therefore not penalized to achieve the nitrification objective since the aeration control makes it possible to maintain an identical ammonium concentration throughout the day. It is found that the regulation object of the invention operates properly and anticipates the peak load, under the conditions of variable feed rate (Figure 7). Except for a short period when the load is such that saturation in air velocity has been reached. Indeed, the aeration was restrained so that it does not exceed the 15 Nm / h to remain representative of the majority of the installations. Finally, in the case of use with variable feed rates and recirculation rates (Figure 8), which is not the representative case of industrial facilities, the results, although degraded, remain significantly higher than the state of the art. These results are obtained using a control block such as that illustrated in Figure 5 which corresponds to a preferred embodiment of the method according to the invention. As shown, a water to be treated 51 is directed to a biological reactor 52, the feed rate data Q and incoming load being associated with this water to be treated. It is recalled that, according to the principle of the invention, the nitrification and denitrification processes are carried out within the same chamber of the reactor (the reactor having in this case a single compartment) and that these processes of nitrification and denitrification are essentially simultaneous. Measurements are made concerning this water to be treated using a first regulator 53, called "Feedforward" regulator, which notably performs a continuous measurement of the charge between C _VEDD . The controller also receives information about the output load setpoint Cv _C0NSIGNE . Another regulator 54, called "Feedback", makes it possible to collect information, and in particular the output load C _vs measured continuously. The regulator 54 also receives the information concerning the output load setpoint Cv _C0NSIGNE . The controller used for the feedback loop can be of the PID (Proportional Proportional Controller) or PFC (Predictive Functional Control) type. These controllers are set to each issue a command. The corresponding commands are processed so as to act on the velocity of the air V _3J1 ., Applying the law of control of the velocity of the air as a function of time: Vair (t) = ax (H (O * Cv _EDD (ή- Cv _Setpoint (*)> β + Vair _Mtmaction The coefficients α and β can be obtained by measurements or by calculation, as indicated above. The feedback is the required change in airflow calculated by the feedback loop controller. In addition, in the case of the pilot unit described above, the phase delay function H (s) is expressed as follows:

Par ailleurs, la commande issue du régulateur 54 est obtenue par le calcul correspondant aux instructions suivantes : Cv_C0Nsi_GNE(t) - Cv₅= erreur(t), , v TA x Q x eit) ou erreur(t) = - ^¹ -, ^W (1000 x 1000 x h_mat x S) Furthermore, the command from the regulator 54 is obtained by the calculation corresponding to the following instructions: Cv _C0N if _GNE (t) - Cv ₅ = error (t) _,, v TA x Q x eit) or error (t) = - ^ ¹ -, ^W (1000 x 1000 xh _mat x S)

le signal e(t) étant défini par :, e(t) = [NH₄]_comGNE(t)-[NH₄]_s(t),the signal e (t) being defined by _:, e (t) = [NH ₄ ] _comGNE (t) - [NH ₄ ] _s (t),

Avec les paramètres déjà énoncés (Veau=l,2 m/h et hmat = 2,75 m), on définit erreur(t) de la façon suivante : erreur(t) = 1,05.10^"2e(t) With the parameters already stated (Cal = l, 2 m / h and hmat = 2.75 m), one defines error (t) in the following way: error (t) = 1.05.10 ^"2 e (t)

Claims

REVENDICATIONS

1. Procédé de traitement d'eaux, à l'aide d'un réacteur biologique en vue d'abattre la pollution azotée, dite charge entrante, contenue dans lesdites eaux, ledit réacteur intégrant une biomasse aérée par injection d'air, ledit procédé incluant au moins une étape de régulation de Ia vitesse dudit air injecté, caractérisé en ce qu'il comprend une mesure en continu de ladite charge entrante en N-NH₄ (Cv_EDD) contenue dans lesdites eaux, ladite charge entrante mesurée étant pondérée par un décalage temporel pour en déduire une commande de la vitesse Vair dudit air à injecter en fonction du temps. 1. A method of water treatment, using a biological reactor to reduce nitrogen pollution, said incoming charge, contained in said water, said reactor incorporating a biomass aerated by air injection, said method including at least one step of regulating the speed of said injected air, characterized in that it comprises a continuous measurement of said incoming charge of N-NH ₄ (Cv _EDD ) contained in said waters, said measured incoming charge being weighted by a time offset to deduce a control of the air velocity of said air to be injected as a function of time.

2. Procédé de traitement d'eaux selon la revendication 1, caractérisé en ce que ladite ou lesdites étapes de régulation sont basées sur la loi suivante de commande de la vitesse de l'air en fonction du temps : Vair (t) = α(H(t) * Cv_EDD (t) - Cv_C0NSIGNE) + β, dans laquelle Cv_CONSIGNE est une consigne de charge de sortie et H(t) est la transformée de Laplace inverse d'une fonction de retard de phase H(S). 2. Water treatment method according to claim 1, characterized in that said one or more regulation steps are based on the following law of control of the speed of the air as a function of time: Vair (t) = α ( H (t) * Cv _EDD (t) - Cv _C0NSIGN ) + β, in which Cv _SETPOINT is an output load setpoint and H (t) is the inverse Laplace transform of a phase delay function H (S ).

3. Procédé de traitement d'eaux selon l'une des revendications 1 et 2, caractérisé en ce que ladite ou lesdites étapes de régulation sont basées sur une loi de commande s 'exprimant de la façon suivante : Vair (t) = α(H(t) x Cv_EDD (t) - Cv_C0NSIGNE(t + Δt)) + β, 3. Water treatment method according to one of claims 1 and 2, characterized in that said one or more regulation steps are based on a control law s' expressing as follows: Vair (t) = α ( H (t) x Cv _EDD (t) - Cv _C0NSIGN (t + Δt)) + β,

4. Procédé de traitement d'eaux selon l'une des revendications 2 et 3, caractérisé en ce que ladite fonction de retard de phase est du type :4. Water treatment method according to one of claims 2 and 3, characterized in that said phase delay function is of the type:

où : - n est un paramètre d'ajustement de la diffusion dans ledit ou lesdits réacteurs ; - V est le volume apparent dudit ou desdits réacteurs biologiques ; - Q est le débit d'alimentation desdites eaux à traiter ; - s est la transformée de Laplace de la variable t. where: - n is an adjustment parameter of the diffusion in said reactor (s); V is the apparent volume of said one or more biological reactors; - Q is the feed rate of said water to be treated; - s is the Laplace transform of the variable t.

5. Procédé de traitement d'eaux selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'il comprend au moins une étape de mesure d'une charge de sortie (Cv_s) et/ou d'une concentration en ammonium dissous contenues dans lesdites eaux traitées. 5. Water treatment method according to any one of claims 1 to 4, characterized in that it comprises at least one step of measuring an output load (Cv _s ) and / or a concentration of dissolved ammonium contained in said treated water.

6. Procédé de traitement d'eaux selon la revendication 5, caractérisé en ce que ladite mesure de ladite charge de sortie (Cv₅) est exploitée pour ajuster ladite consigne par application de l'équation suivante : Cv_CONSi_GNE(t) - Cv₅= erreur(t), , v 24 x Qx e(t) ou erreur(t) = ±-¹ -, ^w (1000 x 1000 x h_mat x S) le signal e(t) étant défini par :, e(t) = [NH₄]_s(t)-[NH₄]_CONSIGNE(t),6. A method of water treatment according to claim 5, characterized in that said measurement of said output load (Cv ₅ ) is used to adjust said setpoint by applying the following equation: Cv _CONS i _GN E (t) - Cv ₅ = error (t),, v 24 x Qx e (t) or error (t) = ± - ¹ -, ^w (1000 x 1000 x h _mat x S) the signal e (t) being defined by:, e (t) = [NH ₄ ] _s (t) - [NH ₄ ] _SETPOINT (t),

h_mat étant la hauteur de ladite biomasse dans ledit réacteur et S étant la surface dudit ou desdits réacteurs, , le volume résultant du produit h_mat. par S représentant le volume aéré dudit réacteur biologique. h _mat being the height of said biomass in said reactor and S being the surface of said one or more reactors, the resulting volume of the product h _mat . by S representing the aerated volume of said biological reactor.

7. Procédé de traitement d'eaux selon l'une des revendications 5 et 6, caractérisé en ce que ladite mesure d'une charge de sortie (Cv₈) et/ou d'une concentration en ammonium dissous contenues dans lesdites eaux traitées, est réalisée en continu 8. Procédé de traitement d'eaux selon l'une des revendications 5 et 6, caractérisé en ce qu'il comprend une étape de transformation dudit signal e(t) réalisée de telle sorte que : - si e > 0, f(e) = e ; - si e ≤ 0, f(e) = 1 - exp(-k.e), où k > 0. 9. Dispositif pour la mise en œuvre du procédé de traitement d'eaux selon l'une quelconque des revendications 1 à 8, comprenant au moins un réacteur intégrant une biomasse aérée par injection d'air et des moyens de régulation de la vitesse dudit air injecté, caractérisé en ce qu'il comprend : - des moyens de mesure en continu de ladite charge entrante ; - des moyens de paramétrage d'une consigne de charge de sortie et/ou d'une concentration de sortie (Cv_C0NSIGNE) ; - des moyens de calcul destinés à agir sur lesdits moyens de régulation sur la base d'une loi de commande dans laquelle ladite charge entrante mesurée est notamment pondérée par un décalage temporel pour en déduire une commande de la vitesse dudit air à injecter en fonction du temps. 10. Dispositif selon la revendication 9, caractérisé en ce qu'il comprend une boucle de rétroaction comprenant des moyens de mesure d'une charge de sortie (Cv₈) et/ou d'une concentration en ammonium dissous contenues dans lesdites eaux traitées, et des moyens de comparaison de ladite charge de sortie avec ladite consigne de charge de sortie. 11. Dispositif selon la revendication 10, caractérisé en ce que lesdits moyens de comparaison sont reliés auxdits moyens de calcul, en vue d'ajuster ladite consigne prise en compte dans ladite loi de commande. 7. Water treatment method according to one of claims 5 and 6, characterized in that said measurement of an output charge (Cv ₈ ) and / or a dissolved ammonium concentration contained in said treated water, is carried out continuously 8. Water treatment method according to one of claims 5 and 6, characterized in that it comprises a step of transforming said signal e (t) carried out so that: - if e> 0 , f (e) = e; if e ≤ 0, f (e) = 1 - exp (-ke), where k> 0. 9. Device for implementing the water treatment method according to any one of claims 1 to 8, comprising at least one reactor integrating a biomass aerated by air injection and means for regulating the speed of said injected air, characterized in that it comprises: - means for continuously measuring said incoming charge; means for setting an output load setpoint and / or an output concentration (Cv _C0NSIGN ); calculating means intended to act on said regulating means on the basis of a control law in which said measured incoming charge is weighted in particular by a time offset to deduce a command of the speed of said air to be injected as a function of the time. 10. Device according to claim 9, characterized in that it comprises a feedback loop comprising means for measuring an output charge (Cv ₈ ) and / or a dissolved ammonium concentration contained in said treated water, and means for comparing said output load with said output load setpoint. 11. Device according to claim 10, characterized in that said comparing means are connected to said calculating means, for adjusting said instruction taken into account in said control law.