FR 2911960 A1 20080801 FR 0752882 A 20070125 La présente invention concerne le domaine de l'analyse et du contrôle en ligne de la composition des effluents liquides, telles par exemple l'eau potable, les eaux urbaines ou industrielles. La présente invention vise plus particulièrement un procédé de détection de la variabilité de la composition d'un tel effluent. Dans les procédés de détection déjà connus, les mesures sont effectuées en ligne, par exemple dans une conduite disposée en amont d'une station de traitement d'eau usée. Un tel procédé permet de détecter une anormalité dans la composition de l'effluent, notamment physico-chimique ou biologique, susceptible de provoquer un incident dans le procédé de traitement situé en aval du lieu où sont effectuées les mesures. L'acquisition en ligne de ces mesures et leur traitement mathématique permettent d'alerter l'opérateur de l'arrivée d'un effluent ayant une composition différente de celle de l'effluent alimentant habituellement la station de traitement. Le cas échéant, l'opérateur prendra des mesures appropriées, par exemple en déviant l'effluent différent vers une autre unité de traitement et/ou en prélevant un échantillon de cet effluent qui sera analysé a posteriori. Parmi les procédés de détection de la variabilité en ligne d'un effluent, certains mettent en oeuvre la spectrophotométrie Ultra-Violet visible (ci-après dénommée UV-VIS ) permettant à partir d'un ensemble de signaux, de représenter la variabilité qualitative de l'effluent. Pour les aspects qualitatifs, le spectre UV-VIS fait office d'empreinte digitale de l'effluent. En d'autres termes, un effluent donné possède un spectre UV-VIS propre comprenant des points particuliers (dénommés points isobestiques) qui, en toute théorie, ne varient pas au cours du temps. Dès lors que l'on détecte un changement dans l'allure du spectre 30 UV mesuré en ligne, on comprend que ce changement d'allure traduit une modification de la composition de l'effluent. L'information obtenue indique à l'opérateur s'il y a conservation globale ou non de la qualité de l'effluent mais ne permet pas d'évaluer la dangerosité de l'effluent en question. Enfin, ces points n'existent pas 35 forcément pour tous les effluents, ce qui rend inefficace cette méthode dans ce cas. Un autre mode connu de traitement des spectres UV, tel que décrit dans les documents FR 2 783 322 et FR 2 787 883, consiste à utiliser la déconvolution des spectres UV, ce traitement se basant sur l'hypothèse selon laquelle chaque spectre UV d'un effluent peut être considéré comme une combinaison linéaire d'un petit nombre de spectres UV connus, appelés spectres de référence, et notamment relatifs à des composés particuliers de l'effluent étudié. Cette méthode suppose donc d'avoir une connaissance a priori de la composition de l'effluent ou des composés susceptibles de s'y trouver, ce qui impose de réaliser des analyses en laboratoires préalables. En outre, une bibliothèque de spectres de composés ou de compositions prédéterminés doit être préalablement réalisée. Par ailleurs, le document EP 1 070 954 mentionne l'utilisation de la spectroscopie dans le spectre UV-visible pour la détermination par déconvolution de spectres de paramètres globaux. Ce document propose une surveillance temporelle de l'effluent de manière à déterminer le moment précis où une pollution éventuelle commence à se manifester. Pour chacun des sites sur lequel on réalisera les mesures, cette méthode nécessite donc une caractérisation préalable de l'effluent via des méthodes conventionnelles, ainsi que la calibration et l'identification préalable de spectres de référence spécifiques à l'effluent analysé. Il en résulte que les procédés de qualification connus sont totalement liés à l'environnement dans lequel ils sont destinés à être mis en oeuvre. En particulier, dans le cas de la spectrophotométrie UV-VIS, cela suppose que l'effluent absorbe dans la gamme de longueur d'onde considérée. Un autre inconvénient des procédés de l'art antérieur est qu'ils ne sont pas capables de traiter d'autres types de mesures, tels notamment le pH, la conductivité, ou la température. Un but de la présente invention est de proposer un procédé de qualification de la variabilité de la composition d'un effluent, dans lequel on réalise une succession de mesures au cours du temps d'au moins un premier et un deuxième paramètres de cet effluent, lequel procédé est indépendant de l'environnement dans lequel il est mis en oeuvre tout en étant capable de traiter des mesures tant paramétriques que non-paramétriques. L'invention atteint son but par le fait que, à chaque pas de temps : on détermine une première dérivée pour chacun des premier et deuxième paramètres ; on détermine au moins une deuxième dérivée entre les premier et deuxième paramètres ; on définit, pour la première dérivée de chacun des paramètres, des premiers domaines logiques comprenant au moins un premier domaine logique normal correspondant à une variabilité normale de ladite première dérivée, et un premier domaine logique anormal correspondant à une variabilité anormale de ladite première dérivée ; on attribue des premières probabilités d'appartenance à chacun des premiers domaines logiques pour la première dérivée de chacun des paramètres ; on définit des deuxièmes domaines logiques comprenant au moins un deuxième domaine logique normal correspondant à une variabilité normale de ladite deuxième dérivée, et un deuxième domaine logique anormal correspondant à une variabilité anormale de ladite deuxième dérivée ; - on attribue des deuxièmes probabilités d'appartenance de la deuxième dérivée à chacun des deuxièmes domaines logiques ; on définit des domaines logiques globaux à partir des premiers et deuxièmes domaines logiques, les domaines logiques globaux comprenant au moins un domaine logique global normal et un domaine logique global anormal ; on attribue des probabilités d'appartenance globales de l'ensemble constitué par les premières dérivées des paramètres et au moins une deuxième dérivée à chacun des domaines logiques globaux ; on qualifie la variabilité de la composition de l'effluent à partir de ces probabilités d'appartenance globales. Ainsi, grâce à l'invention, il n'est pas nécessaire de connaître a priori la composition habituelle de l'effluent pour pouvoir détecter une anomalie dans la variabilité de sa composition. En outre, grâce à l'invention, il n'est pas non plus nécessaire de procéder à des analyses préalables afin d'initialiser le procédé. Tout au contraire, le procédé selon l'invention ne nécessite qu'une succession temporelle de mesures de paramètres pour s'adapter à l'environnement dans lequel il est mis en oeuvre, ledit environnement pouvant être par exemple un réseau d'assainissement ou un réseau d'eau potable ou bien encore un site industriel. Cela signifie que le procédé selon la présente invention peut être mis en oeuvre de façon simple et rapide, sans nécessiter les étapes préalables de calibration ou d'analyse en laboratoire telles que celles effectuées dans les procédés connus. Par ailleurs, le procédé selon l'invention peut avantageusement se baser tant sur l'analyse de la variabilité de mesures de spectres, que sur l'analyse de la variabilité de mesures de grandeurs, tels par exemple mais non exclusivement le pH, le potentiel Rédox, la conductivité ou la température. Un des principaux objectifs de l'invention est de renseigner immédiatement et en temps réel l'opérateur de l'arrivée éventuelle d'un effluent ayant une composition à risque. Pour ce faire, on analyse l'évolution des variations de chacun des paramètres pris séparément d'une part, ainsi que l'évolution de la variation relative entre les premier et deuxième paramètres d'autre part. Au sens de la présente invention, les premières dérivées sont prises par rapport à une même grandeur, de préférence le temps, tandis que les deuxièmes dérivées dénotent la variation d'un des paramètres par rapport à l'autre, de préférence au cours du temps. Dans le cas particulier où le paramètre dépend d'une variable non temporelle, la deuxième dérivée, pourra être par exemple la dérivée du paramètre par rapport à ladite variable, pour un pas de temps donné. Autrement dit, on analyse d'une part, à l'aide des premières dérivées, comment chacun des paramètres évolue, de préférence dans le temps. On parlera ici de variances individuelles des paramètres, la première dérivée étant avantageusement une dérivée temporelle. D'autre part, on analyse comment les paramètres évoluent l'un par rapport à l'autre à l'aide des deuxièmes dérivées. On parlera ici de variances croisées desdits paramètres. Au demeurant, les variances croisées permettent notamment de déterminer si les évolutions des paramètres sont similaires ou à tout le moins corrélées. De manière préférentielle, le nombre de paramètres est 5 sensiblement plus grand que deux. De manière préférentielle, il existe autant de deuxièmes dérivées que de combinaisons entre les paramètres. L'idée sous jacente est de comparer les variances individuelles et croisées de manière à établir une indication de la variance globale de la 10 composition de l'effluent. Le traitement mathématique permettant la mise en oeuvre du procédé selon la présente invention repose sur la théorie de la logique floue, bien connue par ailleurs, si bien que la comparaison mentionnée ci-dessus sera obtenue par une opération d'inférence, bien connue par 15 ailleurs. En fonction de l'indication de variance globale, on pourra qualifier la variation de la composition de l'effluent de normal ou d'anormal, ou par tout autre qualificatif approprié. De manière plus précise, l'attribution des premières probabilités 20 d'appartenance à chacun des premiers domaines logiques pour la première dérivée de chacun des paramètres permet de caractériser la variance individuelle des paramètres, tandis que l'attribution des deuxièmes probabilités d'appartenance de la deuxième dérivée à chacun des deuxièmes domaines logiques permet quant à elle de caractériser la 25 variance croisée des paramètres. Enfin, on comprend que l'étape d'attribution des probabilités d'appartenance globales permet de caractériser la variance globale des paramètres. Le traitement mathématique appliqué dans la présente invention 30 utilise préférentiellement au moins une base de données glissante dont la longueur est variable, la base de données étant constituée par les valeurs des paramètres à des pas de temps précédents, grâce à quoi le procédé s'adapte seul aux évolutions naturelles de l'effluent. Par glissante , on entend que la base de données contient un 35 nombre prédéterminé de valeurs, une nouvelle entrée dans la base s'accompagnant de la suppression de la plus ancienne valeur stockée. De manière particulièrement intéressante, cela implique que l'opérateur n'a pas à fixer de seuils pour déterminer si l'effluent varie, ni d'adapter un tel seuil à l'évolution naturelle de l'effluent. Au contraire, grâce à l'invention, de tels seuils sont déterminés 5 automatiquement par le procédé, en fonction des données enregistrées dans la base de données. Comme on l'a déjà mentionné, le procédé selon l'invention ne nécessite aucune mesure, préalable ou non, de l'effluent en laboratoire. Par ailleurs, aucune calibration préalable de l'outil mathématique utilisé 10 n'est nécessaire. Ainsi s'affranchit-on de mesures ciblées sur certaines caractéristiques de l'effluent, tant et si bien que le procédé selon l'invention offre une approche globale du risque de dangerosité lié à la variabilité de la composition de l'effluent. 15 Enfin, on pourra choisir un critère de qualification simple, par exemple un signal rouge signifiant que l'effluent a fortement varié, ceci impliquant un risque fort de dysfonctionnement du traitement, associé par exemple à un signal vert signifiant que l'effluent n'a pas ou peu varié. 20 Avantageusement, pour la première dérivée de chacun des paramètres, les premiers domaines logiques sont définis à partir d'un multiple de l'écart-type évalué à partir d'un ensemble constitué par des valeurs de la première dérivée dudit paramètre déterminées à des pas de temps précédents. 25 De manière préférentielle, pour la première dérivée de chacun des paramètres, les premiers domaines logiques sont également définis à partir de la moyenne ou de la médiane d'un ensemble constitué par des valeurs de la première dérivée dudit paramètre déterminées à des pas de temps précédents. 30 Avantageusement, les deuxièmes domaines logiques sont définis à partir d'un multiple de l'écart type évalué à partir d'un ensemble constitué par des valeurs de ladite au moins une deuxième dérivée déterminées à des pas de temps précédents. Avantageusement, les domaines logiques globaux sont définis à 35 partir d'un multiple de l'écart-type évalué à partir d'un ensemble constitué par les valeurs des premières dérivées des paramètres et des valeurs de ladite au moins une deuxième dérivée déterminées à des pas de temps précédents. Les valeurs déterminées à des pas de temps précédents mentionnées ci-dessus sont avantageusement stockées dans au moins une base de données glissante, telle que celle mentionnée ci-dessus. En vue d'affiner la qualification de la variabilité de la composition de l'effluent, on définit en outre, pour la première dérivée de chacun des paramètres, un premier domaine logique intermédiaire correspondant à une variabilité peu normale de ladite première dérivée. On définit en outre, pour ladite au moins une deuxième dérivée, un deuxième domaine logique intermédiaire correspondant à une variabilité peu normale de ladite deuxième dérivée. De surcroît, on peut également définir un ou plusieurs domaines logiques globaux intermédiaires, correspondant à une variation peu normale de la composition de l'effluent. En pratique, on pourra utiliser un signal orange complétant les signaux rouge et vert mentionnés ci-dessus, ledit signal orange traduisant une variation de l'effluent dans un domaine acceptable, le risque étant estimé comme étant modéré mais nécessitant une vigilance accrue (risque potentiel d'évolution vers un signal rouge). Selon un premier mode de réalisation, les probabilités d'appartenance globales sont calculées à partir d'une somme pondérée des premières et deuxièmes probabilités d'appartenance. Selon un second mode de réalisation, on définit en outre : - des troisièmes domaines logiques comprenant au moins un troisième domaine logique normal et un troisième domaine logique anormal, les troisièmes domaines logiques étant définis à partir des premiers domaines logiques ; des quatrièmes domaines logiques comprenant au moins un quatrième domaine logique normal et un quatrième domaine logique anormal, les quatrièmes domaines logiques étant définis à partir des deuxièmes domaines logiques ; on attribue en outre des troisièmes probabilités d'appartenance d'un ensemble constitué par les premières dérivées des paramètres à chacun des troisièmes domaines logiques ; - on attribue en outre des quatrièmes probabilités d'appartenance d'un ensemble constitué par au moins ladite au moins une deuxième dérivée à chacun des quatrièmes domaines logiques ; - les domaines logiques globaux étant définis à partir des troisièmes et quatrièmes domaines logiques. Dans ce second mode de réalisation, on insère des domaines logiques supplémentaires, à savoir les troisièmes et quatrièmes domaines logiques de manière à affiner la qualification de la variabilité de la composition de l'effluent et à éviter les fausses alertes. De manière préférentielle, le troisième domaine logique normal correspond à l'union des premiers domaines logiques normaux , le troisième domaine logique anormal correspond à l'union des premiers domaines logiques anormaux , le quatrième domaine logique normal correspond à l'union des deuxièmes domaines logiques normaux , tandis que le quatrième domaine logique anormal correspond à l'union des deuxièmes domaines logiques anormaux . Avantageusement, on définit en outre un troisième domaine logique peu normal qui correspond de préférence à l'union des premiers domaines logiques peu normaux , et un quatrième domaine logique peu normal qui correspond de préférence à l'union des deuxièmes domaines logiques peu normaux . De manière préférentielle, on définit cinq domaines logiques globaux. Le premier domaine logique global correspond de préférence à l'union du troisième domaine logique normal et du quatrième domaine logique normal . Le deuxième domaine logique global correspond de préférence à l'union de l'union du troisième domaine logique normal et du quatrième domaine logique peu normal et de l'union du troisième domaine logique peu normal et du quatrième domaine logique normal . Le troisième domaine logique global correspond de préférence à l'union de l'union du troisième domaine logique normal et du quatrième domaine logique anormal avec l'union du troisième domaine logique peu normal et du quatrième domaine logique peu normal , et de l'union du troisième domaine logique anormal et du quatrième domaine logique normal . Le quatrième domaine logique global correspond de préférence à l'union de l'union du troisième domaine logique peu normal et du quatrième domaine logique anormal , et de l'union du troisième domaine logique anormal et du quatrième domaine logique peu normal . Le cinquième domaine logique global correspond de préférence à l'union du troisième domaine logique anormal et du quatrième 10 domaine logique anormal . Enfin, après les avoir calculées, on attribue les probabilités d'appartenance globales de l'ensemble constitué par les premières dérivées des paramètres et les secondes dérivées aux domaines logiques globaux. 15 Avantageusement, les probabilités d'appartenance globales sont calculées à partir d'une somme pondérée des troisièmes et quatrièmes probabilités d'appartenance. De manière préférentielle, ladite somme pondérée fait intervenir au moins un coefficient de pondération dynamique, par exemple l'écart-type 20 évalué à partir d'un ensemble constitué par des valeurs de la première dérivée déterminées à des pas de temps précédents. Selon une première variante, le domaine logique global qualifiant la variabilité de la composition de l'effluent est celui pour lequel la probabilité d'appartenance globale est maximale. 25 De manière préférentielle, seul le cinquième domaine logique global est utilisé pour déterminer la variation de la composition de l'effluent. On comprend que ce cinquième domaine logique global constitue celui pour lequel la probabilité d'une variabilité particulièrement anormale de la composition de l'effluent est la plus forte. 30 Selon une deuxième variante, le domaine logique global qualifiant la variabilité de la composition de l'effluent est le domaine logique le plus anormal pour lequel la probabilité d'appartenance globale est non nulle. Selon une troisième variante, le domaine logique global qualifiant la 35 variabilité de la composition de l'effluent est le domaine logique le plus anormal pour lequel la probabilité d'appartenance globale est supérieure à un seuil prédéterminé. Selon un mode de réalisation préférentiel de l'invention, le premier paramètre est l'absorbance de l'effluent mesuré pour une première longueur d'onde. De manière préférentielle, chacun des paramètres correspond à l'absorbance de l'effluent pour une longueur d'onde donnée. Préférentiellement, on mesure à chaque pas de temps l'absorbance de l'effluent pour une pluralité de longueurs d'onde comprises de 10 préférence entre 200 nm et 700 nm. Avantageusement, la première dérivée est une dérivée de l'absorbance par rapport au temps pour une longueur d'onde donnée, tandis que la deuxième dérivée est une dérivée de l'absorbance par rapport à la longueur d'onde calculée à un pas de temps donné. 15 La présente invention concerne également un dispositif de suivi de la variabilité de la composition d'un effluent mettant en oeuvre le procédé selon la présente invention. D'autres caractéristiques et avantages de l'invention apparaîtront mieux à la lecture de la description qui suit d'un mode de réalisation de 20 l'invention donné à titre d'exemple non limitatif. La description se réfère aux figures annexées sur lesquelles : - la figure 1 est une représentation schématique d'un dispositif de suivi de la variabilité de la composition d'un effluent mettant en oeuvre le procédé selon la présente invention ; 25 la figure 2 est un diagramme représentant les étapes du mode de mise en oeuvre préféré de la présente invention ; - la figure 3 est un graphique représentant des fonctions d'appartenance en triangle permettant de déterminer des probabilités d'appartenance des premières dérivées aux 30 premiers domaines logiques ; et la figure 4 est un graphique représentant des fonctions d'appartenance en triangle permettant de déterminer des probabilités d'appartenance des deuxièmes dérivées aux deuxièmes domaines logiques. 35 A l'aide des figures 1 à 4, on va maintenant décrire plus en détail un mode préféré de mise en oeuvre du procédé selon la présente invention dont l'objectif est d'évaluer la variabilité temporelle de la composition d'un effluent de l'industrie chimique. La figure 1 représente un dispositif 10 conforme à l'invention mettant en oeuvre le procédé selon l'invention. Comme on le voit sur cette figure, le dispositif 10 comprend des moyens de mesure 12 destinés à mesurer l'absorbance de l'effluent s'écoulant dans une conduite 14, cette dernière étant préférentiellement mais non exclusivement disposée en amont d'une station de traitement 15 de l'effluent. Les valeurs mesurées par les moyens de mesure 12 sont stockées dans des moyens de mémoire 16, tandis qu'un microprocesseur 18 est prévu pour effectuer des opérations de calcul conduisant à la qualification de la variabilité de la composition de l'effluent. Le résultat des opérations de calcul est préférentiellement donné par une interface visuelle 20 comportant trois feux de couleur, par exemple vert, orange et rouge caractérisant le niveau de normalité de la variabilité de l'effluent. Le procédé selon l'invention va maintenant être expliqué plus en détail à l'aide du diagramme représenté sur la figure 2. Dans ce mode de mise en oeuvre, on réalise au cours du temps une succession de mesures d'une pluralité de paramètres, à savoir l'absorbance de l'effluent pour plusieurs longueurs d'onde. Autrement dit, à chaque pas de temps, on mesure, au cours d'une étape S101, l'absorbance de l'effluent pour plusieurs longueurs d'onde, celles-ci étant comprises entre 200 et 700 nm, de préférence entre 230 et 500 nm. De préférence, les mesures sont effectuées en faisant varier la longueur d'onde par pas de 1 nm, tandis qu'un pas de temps t dure deux heures. Le diagramme de la figure 2 représente les étapes qui sont réalisées à chaque pas de temps. Ainsi dispose-t-on, dans cet exemple de mise en oeuvre où l'on mesure l'absorbance entre 230 nm et 500 nm, de 271 paramètres qui ont été mesurés simultanément toutes les deux heures sur l'effluent. On comprend donc qu'au sens de l'invention, le premier paramètre est l'absorbance à la longueur d'onde de 230 nm, le deuxième paramètre est l'absorbance à la longueur d'onde de 231 nm, ainsi de suite jusqu'au 271ème paramètre, qui est l'absorbance à 500 nm. Evidemment, le nombre de paramètres n'est pas limité à 271 et l'on peut tout à fait choisir un autre nombre de paramètres sans sortir du cadre de la présente invention. De manière préférentielle mais non nécessairement, les moyens de mémoire comportent une base de données stockant les mesures de paramètres pour une pluralité de pas de temps t. Au cours d'une seconde étape S102, on détermine une première dérivée pour chacun des 271 paramètres en calculant la différence entre deux valeurs d'absorbance consécutives dans le temps à une même longueur d'onde et en divisant le résultat obtenu par un écart de temps At existant entre deux mesures. Ainsi, si l'on considère le premier paramètre, à savoir l'absorbance 15 Abs de l'effluent à 230 nm, la première dérivée du 1eC paramètre d1,1(t) a pour expression : dl l (t) Abs23o (t) ù Abs23o (t ù At) = At 20 Il en est de même pour la première dérivée du 2ème paramètre d1,2(t), qui a pour expression : d 1 2 (t) = AbS231(t) ù AbS231(t ù At) Et ainsi de suite jusqu'au 271ème paramètre. 25 Au cours d'une troisième étape S103, on détermine au moins une deuxième dérivée entre les paramètres. Plus précisément, on définit, dans cet exemple, 270 deuxièmes dérivées des 271 paramètres. La lère deuxième dérivée d2,1(t) est calculée comme étant la différence entre les valeurs d'absorbance à 231 nm et à 30 230 nm, soit : d2,1(t) Abs 231 (t) ù Abs 23o (t 231 ù 230 At On comprend donc que la deuxième dérivée d2,1(t) dénote bien l'écart relatif entre les paramètres Abs231(t) et Abs230(t) et que la fonction t-> d2,1(t) fournit la variation temporelle de l'écart relatif entre ces deux paramètres. Il en est de même pour la 2nde deuxième dérivée d2,2(t), qui a pour expression : AbS232 (t)ù AbS231 (t) 232 ù 231 Et ainsi de suite jusqu'à la 270ème deuxième dérivée. Dans une quatrième étape S104, on définit pour chaque jème première dérivée (j=1 à 271), c'est-à-dire pour la première dérivée de chacun des paramètres, trois premiers domaines logiques (DL,(j), i=1 à 3, j=1 à 271) déterminés à partir des calculs suivants : a) on calcule une moyenne d1,1(t)et un écart-type 6d,,J (t) à partir des 50 valeurs calculées antérieurement de chaque jème première dérivée (j=1à271): 1 49 d1,>(t) = S0ld1,./ (tùitt) 1 49 1 (d .(tù'At)(t)) b) on calcule pour chaque jème première dérivée (j = 1 à 271) l'écart absolu E1,j(t) à la moyenne des premières dérivées : 1,> (t) = d1,1(t)ù d1,1(tl 25 c) Les premiers domaines logiques DEI(j,t)(avec i=1 à 3, j=1 à 271) sont alors définis au pas de temps tcomme suit : - le 1er des premiers domaines logiques de chaque jème première dérivée variation normale de la première dérivé , DLl (j,t), est caractérisé en ce que : c1,;(t) > 261,j(t). Au cours d'une cinquième étape S105, on attribue des probabilités d'appartenance à l'instant t des 271 premières dérivées à chacun des premiers domaines logiques, P(DLI(j,t)) 0=1 à 3, j=1 à 271), en 10 effectuant les sous-étapes suivantes : a) on norme entre -1 et 1 l'écart absolu obtenu pour chaque jème première dérivée, avec la correspondance suivante : Si e,,,(t) = 0 alors e1,, (tl -1; Si e1 (t) >_ 30, (t).d,,1(t) alors 1 e1,, ( 15 - Sinon e1(tll=û1+ 2'e, f(t) 361,, (t) b) on procède à la fuzzyfication (i.e. à l'attribution des probabilités d'appartenance de chaque jème première dérivée aux premiers domaines logiques D4(j,t) de chaque valeur normée e1 1(t)I avec des fonctions d'appartenance en triangle telles que représentées 20 sur la figure 3. Sur la figure 3, on a représenté le cas où, pour la 50ème première dérivée avec une valeur de IIe1,50 (t0 )1= 0.3 à l'instant to, on obtient P(DL,(50,t0))=o, P(DLi(50,t0))=0.7, P(DL; (50,t0))=0.3, qui 25 correspondent aux probabilités d'appartenance de la 50ème première dérivée aux premiers domaines logiques. Au cours d'une sixième étape S106, on définit des deuxièmes domaines logiques, à savoir, pour chaque kème deuxième dérivée trois 30 deuxièmes domaines logiques (DL2(k,t), i=1 à 3, k=1 à 270), déterminés à partir des calculs suivants : a) on calcule la moyenne et l'écart-type des 50 dernières valeurs de chaque deuxième dérivée : 49 d2k(t)=ù1d2k(tùiot) 50 i=o 49 2 ù (d2,1 (t ù iAt) ù d2,k (t )) 49 i=o b) on calcule pour chaque kème deuxième dérivée l'écart absolu E2,k(t) à la moyenne des deuxièmes dérivées : e2,k (t) d2,k (t) - d2,k (t) c) les deuxièmes domaines logiques DLZ(k,t) sont alors définis comme suit : - le 1er des deuxièmes domaines logiques variation normale de la kème deuxième dérivée , DLZ (k, t) , est caractérisé en ce que E2,k(t) 262,k(t). Au cours d'une septième étape S107, on attribue les probabilités d'appartenance des kème deuxièmes dérivées aux deuxièmes domaines logiques, P(DLZ (k, t)) (i=1 à 3, k=1 à 270), en effectuant les sous-étapes suivantes : a) on norme entre -1 et 1 l'écart absolu obtenu, avec la correspondance suivante : Si e2k(t)=0 alors Ile2,k(t)l=-1; Si e'k(t)%362k(t).d2k(t) alors Ile2,k(II =+1; Sinon Ile2k(t)l=ù1i 2 x e2 k (t) 3o-2,k (t) b) on procède à la fuzzyfication (i.e. à l'attribution des probabilités d'appartenance de chaque kème deuxièmes dérivées aux deuxièmes domaines logiques DLZ (k,t) de chaque valeur normée Ile, k (t)I avec et : 6d2k (t)= des fonctions d'appartenance en triangle telles que représentées sur la figure 4. Sur la figure 4, on a représenté le cas où, pour la 50ème deuxième dérivée avec une valeur de IIe250 (t0 )1= 0.4 à l'instant to, on obtient P(DLZ (50, t0 )) = 0 , P(DL22 (50, t0 )) = 0.6 , P(DL32 (50, t0 )) = 0.4 , qui correspondent aux probabilités d'appartenance de la 50ème deuxième dérivée aux deuxièmes domaines logiques. Au cours d'une huitième étape S108, préférentielle mais non nécessaire, on définit des troisièmes domaines logiques à savoir, trois troisièmes domaines logiques ( DL3(t) , i=1 à 3) pour les premières dérivées, ceux-ci correspondant à : a) Pour le 1' des troisièmes domaines logiques DL`3 (t) à l'union des 1ers premiers domaines logiques normaux ; b) Pour le 2ème des troisièmes domaines IogiquesDL233(t), à l'union des 2èmes premiers domaines logiques peu normaux ; c) Pour le 3ème des troisièmes domaines IogiquesDL3(t), à l'union des 3èmes premiers domaines logiques anormaux ; Au cours d'une neuvième étape S109, également préférentielle mais non nécessaire, on calcule la probabilité d'appartenance aux troisièmes domaines logiques de l'ensemble des 271 premières dérivées en faisant la somme des probabilités d'appartenance de chaque jième première dérivée aux premiers domaines logiques pondérée par l'écart-type des premières dérivées : 271 P(DL3 (t)) _ cs1 1 (t) x P(DLi (j, t)), avec i=1 à 3 =1 Au cours d'une dixième étape S110, également préférentielle mais non nécessaire, on définit des quatrièmes domaines logiques, à savoir trois quatrièmes domaines logiques (D4 (t), i=1 à 3) pour les deuxièmes dérivées, ceux-ci correspondant à : a) Pour le 1er des quatrièmes domaines logiques, DL4 (t), à l'union des fers deuxièmes domaines logiques normaux ; b) Pour le 2ème des quatrièmes domaines logiques, DLS (t), à l'union des 2èmes deuxièmes domaines logiques peu normaux ; c) Pour le 3ème des quatrièmes domaines logiques, DL, (t), à l'union des 3èmes deuxièmes domaines logiques anormaux ; Au cours d'une onzième étape S111, également préférentielle mais non nécessaire, on calcule la probabilité d'appartenance de chaque 270 deuxièmes dérivées aux trois quatrièmes domaines logiques en faisant la somme des probabilité d'appartenance de chaque kème deuxièmes dérivées aux deuxièmes domaines logiques pondérée par l'écart-type des deuxièmes dérivées : 270 P(DL4(t)) = a-z(t)xP(DL2(k,t)), avec i=1 à 3 k=1 Au cours d'une douzième étape S112, on définit des domaines logiques globaux à partir desquels on définira la variabilité globale de l'effluent, de préférence cinq domaines logiques globaux ( DL5 (t), avec 1=1 à 5) à partir des troisièmes et quatrièmes domaines logiques. Ils correspondent à : a) Pour le 1er des domaines logiques globaux, DL5 (t), à l'union du 1er 20 troisième domaine logique normal et du 1' quatrième domaine logique normal ; b) Pour le 2ème des domaines logiques globaux, DL5 (t), à l'union de l'union du 1er troisième domaine logique normal et du 2ème quatrième domaine logique peu normal et de l'union du 2ème 25 troisième domaine logique peu normal et du 1eC quatrième domaine logique normal ; c) Pour le 3ème des domaines logiques globaux, DL5 (t), à l'union de l'union du ter troisième domaine logique normal et du 3ème quatrième domaine logique anormal , avec l'union du 2ème 30 troisième domaine logique peu normal et du 2ème quatrième domaine logique peu normal , et de l'union du 3ème troisième domaine logique anormal et du 1er quatrième domaine logique normal ; d) Pour le 4ème des domaines logiques globaux, DL5 (t), à l'union de 35 l'union du 2ème troisième domaine logique peu normal et du 3ème quatrième domaine logique anormal et l'union du 3ème troisième domaine logique anormal , et du 2ème quatrième domaine logique peu normal ; e) Pour le Sème des domaines logiques globaux, DLS (t), à l'union du 3ème troisième domaine logique anormal et du 3ème quatrième domaine logique anormal ; Au cours d'une treizième étape S113, on attribue des probabilités d'appartenance globales (P(DL5(t)), avec 1=1 à 5) de l'ensemble constitué par les premières dérivées des paramètres et les secondes dérivées à chacun des domaines logiques globaux ( DL5(t) , avec 1=1 à 5), de préférence en calculant ll lesdites probabilités comme suit : P(DL5 (t// = P(DL3 (t))x P(DL4 (t)) P(DL5 (t)) [P(DL3 (t))x P(DL'4 (t))j+ [P(DL3 (t))x P(DL4 (t))] P(DL5 (t)) [P(DL3 (t ))x P(DL4 (t))j+ [P(DL3 (t))x P(DL4 (t))j+ 1P(DL3 (t))x P(DL4 (t))] P(DL5 (t)) [P(DL3 (t))x P(DL24 (t))j+ [P(DL3 (t))x P(DL4 (t))] P(DL5 (t)) P(DL3 (t))x P(DL34 (t)) Au cours d'une quatorzième étape S114, on norme le résultat 20 obtenu en divisant chaque probabilité d'appartenance par la somme des probabilités d'appartenance aux 5 domaines logiques globaux. Puis, au cours d'une quinzième étape S115, on qualifie la variabilité de la composition de l'effluent à partir des probabilités d'appartenance globales en utilisant, de préférence, uniquement le 5ème domaine logique 25 global. On effectue les tests suivants (donnés à titres illustratifs): Si la probabilité d'appartenance au 5ème domaine logique global est supérieure à un premier seuil, par exemple 66%, alors cela signifie que la composition de l'effluent a fortement varié. Un feu rouge est affiché par l'interface 20. 30 Si la probabilité d'appartenance au 5ème domaine logique global est comprise par exemple entre 33% et 66% alors cela signifie que la composition de l'effluent a varié. Un feu orange est affiché par l'interface 20. Si la probabilité d'appartenance au Sème domaine logique global 35 est inférieure à un second seuil, par exemple, 33% alors cela signifie que la composition de l'effluent n'a pas varié. Un feu vert est affiché par l'interface 20. Grâce à ce mode de mise en oeuvre, non limitatif, les fausses alertes sont avantageusement minimisées. En effet, il faut qu'une majorité des premières et deuxièmes dérivées présente des probabilités élevées d'appartenance aux domaines logiques anormaux pour que le 5ème domaine logique global présente une probabilité d'appartenance forte et déclenche une alarme. Comme on l'a mentionné ci-dessus, les étapes S107 à S111 ne sont pas indispensables à la mise en oeuvre de l'invention. En effet, une variante particulière conforme à l'invention, consiste à définir les domaines logiques globaux directement à partir des premiers et deuxièmes domaines logiques sans passer par des troisièmes et quatrièmes domaines logiques, à attribuer des probabilités d'appartenance globales de l'ensemble constitué par les premières dérivées des paramètres et les secondes dérivées à chacun des domaines logiques globaux, puis à qualifier la variabilité de la composition de l'effluent à partir des probabilités globales de manière similaire à ce qui est effectué au cours des étapes S112 à S115. REVENDICATIONS 1. Procédé de qualification de la variabilité de la composition d'un effluent, dans lequel on réalise une succession de mesures au cours du temps d'au moins un premier et un deuxième paramètres physico- chimiques de cet effluent, ledit procédé étant caractérisé en ce que, à chaque pas de temps : on détermine une première dérivée pour chacun des premier et deuxième paramètres ; - on détermine au moins une deuxième dérivée entre les premier et deuxième paramètres ; - on définit, pour la première dérivée de chacun des paramètres, des premiers domaines logiques comprenant au moins un premier domaine logique normal correspondant à une variabilité normale de ladite première dérivée, et un premier domaine logique anormal correspondant à une variabilité anormale de ladite première dérivée ; on attribue des premières probabilités d'appartenance à chacun des premiers domaines logiques pour la première dérivée de chacun des paramètres ; on définit des deuxièmes domaines logiques comprenant au moins un deuxième domaine logique normal correspondant à une variabilité normale de ladite deuxième dérivée, et un deuxième domaine logique anormal correspondant à une variabilité anormale de ladite deuxième dérivée ; - on attribue des deuxièmes probabilités d'appartenance de la deuxième dérivée à chacun des deuxièmes domaines logiques ; - on définit des domaines logiques globaux à partir des premiers et deuxièmes domaines logiques, les domaines logiques globaux comprenant au moins un domaine logique global normal et un domaine logique global anormal ; on attribue des probabilités d'appartenance globales de l'ensemble constitué par les premières dérivées des paramètres et au moins une deuxième dérivée à chacun des domaines logiques globaux ; - on qualifie la variabilité de la composition de l'effluent à partir de ces probabilités d'appartenance globales. 2. Procédé de qualification selon la revendication 1, dans lequel, pour la première dérivée de chacun des paramètres, les premiers domaines logiques sont définis à partir d'un multiple de l'écart type évalué à partir d'un ensemble constitué par des valeurs de la première dérivée dudit paramètre déterminées à des pas de temps précédents. 3. Procédé de qualification selon la revendication 1 ou 2, dans lequel, pour la première dérivée de chacun des paramètres, les premiers domaines logiques sont également définis à partir de la moyenne ou de la médiane d'un ensemble constitué par des valeurs de la première dérivée dudit paramètre déterminées à des pas de temps précédents. 4. Procédé de qualification selon l'une quelconque des revendications 1 à 3, dans lequel, les deuxièmes domaines logiques sont définis à partir d'un multiple de l'écart type évalué à partir d'un ensemble constitué par des valeurs de ladite au moins une deuxième dérivée déterminées à des pas de temps précédents. 5. Procédé de qualification selon l'une quelconque des revendications 1 à 4, dans lequel, les domaines logiques globaux sont définis à partir d'un multiple de l'écart type évalué à partir d'un ensemble constitué par des valeurs des premières dérivées des paramètres et des valeurs de ladite au moins une deuxième dérivée déterminées à des pas de temps précédents. 6. Procédé de qualification selon l'une quelconque des revendications 1 à 5, dans lequel on définit en outre, pour la première dérivée de chacun des paramètres, un premier domaine logique intermédiaire correspondant à une variabilité peu normale de ladite première dérivée. 7. Procédé de qualification selon l'une quelconque des revendications 1 à 6, dans lequel on définit en outre, pour ladite au moins unedeuxième dérivée, un deuxième domaine logique intermédiaire correspondant à une variabilité peu normale de ladite deuxième dérivée. 8. Procédé de qualification selon l'une quelconque des revendications 1 à 7, dans lequel on définit en outre un domaine logique global intermédiaire. 9. Procédé de qualification selon l'une quelconque des revendications 1 à 8, dans lequel les probabilités d'appartenance globales sont calculées à partir d'une somme pondérée des premières et deuxièmes probabilités d'appartenance. 10. Procédé de qualification selon l'une quelconque des revendications 1 à 8, dans lequel on définit en outre : des troisièmes domaines logiques comprenant au moins un troisième domaine logique normal et un troisième domaine logique anormal, les troisièmes domaines logiques étant définis à partir des premiers domaines logiques ; des quatrièmes domaines logiques comprenant au moins un quatrième domaine logique normal et un quatrième domaine logique anormal, les quatrièmes domaines logiques étant définis à partir des deuxièmes domaines logiques ; procédé dans lequel : - on attribue en outre des troisièmes probabilités d'appartenance d'un ensemble constitué par les premières dérivées des paramètres à chacun des troisièmes domaines logiques ; on attribue en outre des quatrièmes probabilités d'appartenance d'un ensemble constitué par au moins ladite au moins une deuxième dérivée à chacun des quatrièmes domaines logiques ; procédé dans lequel : - les domaines logiques globaux sont définis à partir des troisièmes et quatrièmes domaines logiques. 11. Procédé de qualification selon la revendication 10, dans lequel les probabilités d'appartenance globales sont calculées à partir d'une somme pondérée des troisièmes et quatrièmes probabilités d'appartenance. 12. Procédé de qualification selon la revendication 11, dans lequel ladite somme pondérée fait intervenir au moins un coefficient de pondération dynamique. 10 13. Procédé de qualification selon l'une quelconque des revendications 1 à 12, dans lequel le domaine logique global qualifiant la variabilité de la composition de l'effluent est celui pour lequel la probabilité d'appartenance globale est maximale. 15 14. Procédé de qualification selon l'une quelconque des revendications 1 à 12, dans lequel le domaine logique global qualifiant la variabilité de la composition de l'effluent est le domaine logique le plus anormal pour lequel la probabilité d'appartenance globale est non nulle. 20 15. Procédé de qualification selon l'une quelconque des revendications 1 à 12, dans lequel le domaine logique global qualifiant la variabilité de la composition de l'effluent est le domaine logique le plus anormal pour lequel la probabilité d'appartenance globale est 25 supérieure à un seuil prédéterminé. 16. Procédé de qualification selon l'une quelconque des revendications 1 à 15, dans lequel la première dérivée est une dérivée temporelle. 30 17. Procédé de qualification selon l'une quelconque des revendications 1 à 16, dans lequel le premier paramètre est l'absorbance de l'effluent mesuré pour une première longueur d'onde. 18. Procédé de qualification selon l'une quelconque des revendications 35 1 à 17, dans lequel chacun des paramètres est l'absorbance de l'effluent pour une longueur d'onde donnée.5 19. Procédé de qualification selon la revendication 18, dans lequel la première dérivée est une dérivée de l'absorbance par rapport au temps pour une longueur d'onde donnée, tandis que la deuxième dérivée est une dérivée de l'absorbance par rapport à la longueur d'onde calculée à un pas de temps donné. 20. Dispositif (10) de suivi de la variabilité de la composition d'un effluent mettant en oeuvre le procédé selon l'une quelconque des revendications 1 à 19 et comportant des moyens de mesure (12), des moyens de mémoire (16) et un microprocesseur (18).