FR 2871153 A1 20051209 FR 0405970 A 20040602 Procédé de traitement d'eaux à l'aide d'un réacteur biologique, dans lequel la vitesse d'air injecté 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 oeuvre 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. Cependant, le niveau de cette aération doit être correctement maîtrisé afin de faire coexister dans le même réacteur deux réactions antagonistes: l'une nécessitant de l'oxygène, la nitrification (N), et l'autre son absence, la dénitrification (DN). En effet, 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 sortir de celuici. 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. II 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 FR2 764 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 la consigne en oxygène dissous. 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 oeuvre un réacteur biologique intégrant une biomasse aérée par injection d'air dans lequel la commande de l'injection d'air est optimisée par rapport aux techniques de l'art antérieur. 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 oeuvre. 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, à 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 la vitesse dudit air injecté, caractérisé en ce qu'il comprend une mesure en continu de ladite charge entrante en N-NH4 (CvEDD) 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 air à injecter en fonction du temps. On note que la charge entrante CvEDD peut s'exprimer de la façon suivante: CvEDD = Q x CNH4 X 24 1000x1000xVéé où : Q est le débit d'alimentation (exprimé ici en litre par heure), CNH4 la concentration mesurée en entrée en ammonium (mg N-NH4 par litre) et Vaéré le volume du réacteur aéré (m3). 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 Vaéré =nf x S x h,,,at et S = nf x S, (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. 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. 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) = a(H(t) * CVEDD (t) CVCONSIGNE) + R, dans laquelle CvCONSIGNE 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/m3 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: Vaér(t)=a ÉCve(t)+ I3, avec _ Qa;r Vair nf x Su 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 10 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-NI-14 constant. Ainsi, la concentration d'ammonium de sortie s'exprime alors en fonction de la concentration d'entrée diluée: [NH4] (t) [NH4Je(t-Tg)-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 (Cve) est égale à la différence de charge entre l'entrée (CVEDD) convolué par H(t) (transformé de Laplace inverse d'une fonction de retard de phase H(s)) et la sortie (Cvs). 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) = a(H(t) * CVEDD (t) CVCONSIGNE) + R 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.02NH4 = 4,57 x Masse de N-NH4 éliminée Dénitrification: B.02NO3 = -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 a et (3 (a est égale à 25.7 par la théorie et 24.1 par la mesure, (3 est égale à -4.2 par la théorie et -3.6 par la mesure), R2 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) = a(H(t) * CVEDD (t) CVCONSIGNE(t + At)) + ( t, De cette façon, le terme Cvconsigne(t+At) 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: \n H(s) = où :1 1+ Vs nQ j 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 15 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 (Cvs) et/ou d'une concentration en ammonium dissous contenues dans lesdites eaux traitées. Dans ce cas, ladite mesure de ladite charge de sortie (Cvs) e s t préférentiellement exploitée pour ajuster ladite consigne par application de l'équation suivante: (1000 x 1000 x hmat x S)' le signal e(t) étant défini par:, e(t) _ LNH4ICONSIGNE (t)-LNH4Is(t) hn,at é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 hmat par S représentant 10 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] consigne(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, 15 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 25 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. CVCONSIGNE(t) - Cvs= erreur(t), 24xQxe(t) 5 où erreur(t) = Selon un premier mode de réalisation, ladite mesure d'une charge de sortie (Cvs) 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 5 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: sie>0,f(e)=e; si e s 0, f(e) = 1 exp(-k.e), où k > O. 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, +oo [, 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=1. Ainsi, au prix d'une faible déstabilisation du système, entraînée par 25 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)=x2, 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 oeuvre du procédé de traitement d'eaux tel que décrit précédemment, comprenant au moins un réacteur intégrant une biomasse aéré 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 (CvCONSIGNE) 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. De tels moyens de calcul pourront comprendre la mise en oeuvre de moyens logiciels. 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 (Cvs) 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é, 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 NNH4/m3/J. La vitesse en eau d'alimentation V. moyenne correspondante est de 1,2 mlh 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. II n'existe pas de zone de sous ou sur-aération, avant et après le 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. 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 CVEDD. Le régulateur reçoit de plus une information concernant la consigne de charge de sortie CvCONSIGNEÉ Un autre régulateur 54, dit Feedback , permet de recueillir des informations, et notamment la charge de sortie Cvs mesurée en continu. Le régulateur 54 reçoit également l'information concernant la consigne de charge de sortie CVCONSIGNEÉ Le régulateur utilisé pour la boucle de rétroaction peut être du type PID (régulateur par action proportionnelle) ou PFC (Predictive 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 Vair, en appliquant la loi de commande de la vitesse de l'air en fonction du temps: Vair (t)= a x (H (t)* CVEDD (t) CVCons,gne (t))+ N + VairRétroaction Les coefficients a 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: / \4 1+Vs 4Q / Par ailleurs, la commande issue du régulateur 54 est obtenue par le calcul correspondant aux instructions suivantes: CvCONSIGNE(t) - Cvs = erreur(t), 24xQxe(t) où erreur(t) = (1000 x 1000 x h at x S)' le signal e(t) étant défini par:, e(t) = LNH4]CONSIGNE(t) [NH41s(t), Avec les paramètres déjà énoncés (Veau=1,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) H(s) _ 18 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 la vitesse dudit air injecté, caractérisé en ce qu'il comprend une mesure en continu de ladite charge entrante en N-NH4 (CVEDD) 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. 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) * CVEDD (t) CVCONSIGNE) + R, dans laquelle CvCONSIGNE 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). 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) = a(H(t) x CVEDD (t) CVCONSIGNE(t + At)) + f3, 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: H(s) = 1 V 1+ s nQ j n est un paramètre d'ajustement de la diffusion dans ledit ou lesdits réacteurs; n où : - 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. 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 (Cvs) et/ou d'une concentration en ammonium dissous contenues dans lesdites eaux traitées. 6. Procédé de traitement d'eaux selon la revendication 5, caractérisé en ce que ladite mesure de ladite charge de sortie (Cvs) est exploitée pour ajuster ladite 10 consigne par application de l'équation suivante: CVCONSIGNE(t) - Cvs= erreur(t), 24xQxe(t) où erreur(t) = (1000 x 1000 x h t x S)' le signal e(t) étant défini par:, e(t) = [NHals (t)-[NH 4IcoNSJGNE(t), hmat é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 hmat. par S représentant le volume aéré dudit réacteur biologique. 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 (Cvs) 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) 25 réalisée de telle sorte que: sie>0,f(e)=e; - si e s 0, f(e) = 1 exp(-k.e), où k > O. 9. Dispositif pour la mise en oeuvre 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é 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 (CvCONS[GNE) 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 (Cvs) 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.