FR 3024254 A1 20160129 FR 1457222 A 20140725 Procédé de détection d'anomalies dans un réseau de distribution, en particulier distribution d'eau » La présente invention concerne un procédé de détection d'anomalies dans un système d'approvisionnement en fluide, en particulier un réseau de distribution de fluides newtoniens, encore plus particulièrement pour la distribution d'eau, mais aussi par exemple dans les réseaux urbains de chauffage et/ou de climatisation. La performance opérationnelle est au coeur de la gestion de tels réseaux . Le niveau de performance peut être sensiblement amélioré grâce à des outils de détection et localisation d'anomalies hydrauliques, par exemple sur les réseaux de transport et de distribution d'eau potable. Dans de tels réseaux, différents types d'incidents peuvent survenir et avoir des conséquences graves pour l'exploitant. Il peut s'agir de fuites, d'interruption de fourniture aux consommateurs, défauts de fonctionnement d'un composant, erreur de manipulation lors d'une intervention ou suite à une intervention, malveillance, comportement anormal ou dysfonctionnement chez un consommateur. Sur des réseaux très étendus et partiellement inaccessibles (enterrés) la surveillance directe est peu efficace ou souvent impossible. La surveillance par des instruments de mesure n'est pas non plus suffisante, car une anomalie locale peut avoir des effets sur une zone étendue du réseau ou même sur le réseau tout entier, ce qui complique le diagnostic. On connaît des méthodes de détection basées sur des techniques statistiques. Elles ne permettent ni de caractériser ni de localiser l'anomalie. D'autres méthodes de détection connues sont basées sur des modèles hydrauliques utilisant des algorithmes évolutionnaires (algorithmes génétiques) pour détecter et/ou localiser des anomalies. Ces méthodes évolutionnaires sont très gourmandes en ressources de calcul, ce qui les rend difficilement applicables à grande échelle. 3024254 2 Le but de la présente invention est de proposer un procédé efficace pour détecter et/ou localiser et/ou caractériser et/ou quantifier une anomalie dans un réseau de distribution pouvant le cas échéant être vaste, comme l'est un réseau de distribution d'eau potable. 5 Suivant l'invention, le procédé de détection d'anomalies dans un réseau de distribution, comprend : - établissement pour le réseau d'un modèle direct comprenant les éléments suivants : 10 o des données structurelles du réseau, o des données opérationnelles relatives à des paramètres choisis comme paramètres d'entrée primaires décrivant un scénario opérationnel, o des lois reliant lesdites données structurelles, les paramètres 15 d'entrée primaires et des paramètres de sortie primaires (simuler le modèle); - Exécution du modèle direct pour déterminer des valeurs théoriques des paramètres de sortie primaires (résoudre le problème inverse); - obtention de mesures des paramètres de sortie primaires ; 20 - comparaison des valeurs théoriques et des mesures ; - dans des cas où la comparaison révèle au moins un paramètre de sortie primaire anormal du fait d'un écart important entre valeur théorique et mesure, établissement d'un modèle inverse qui est dérivé du modèle direct en introduisant l'au moins un paramètre de 25 sortie primaire anormal comme un paramètre d'entrée secondaire du modèle inverse, le modèle inverse comportant de plus au moins un paramètre de sortie secondaire en plus des paramètres de sortie primaires conservés ; - exécution du modèle inverse. 30 Les données structurelles sont typiquement relatives à la topographie du réseau, à la nature et aux caractéristiques de ses composants etc. On appelle « paramètre d'entrée primaire » un paramètre de fonctionnement. 35 Les paramètres de fonctionnement sont typiquement variables dans le temps. Par exemple dans un réseau de distribution d'eau un débit, une pression etc., 3024254 3 dont on utilise la valeur comme donnée d'entrée pour le modèle direct. Cette valeur peut être fournie par un instrument de mesure où être connue d'une autre manière, comme par exemple une statistique de consommation en fonction de la date et de l'heure. . En général les paramètres d'entrée primaires 5 comprennent des paramètres descriptifs des conditions aux limites du réseau. Par exemple, dans un réseau hydraulique, le débit d'alimentation du réseau et les débits fournis aux consommateurs. D'une manière générale, l'adjectif « primaire » signifie « relatif au modèle 10 direct » alors que l'adjectif « secondaire » signifiera « relatif au modèle inverse ». Un « scénario » est un cas de fonctionnement du réseau. Par exemple, dans un réseau hydraulique, le scénario commence avec un certain niveau de 15 remplissage dans un réservoir d'alimentation, et au cours du déroulement du scénario le réservoir est réapprovisionné avec un certain débit, tandis qu'il fournit de l'eau au réseau en fonction des débits prélevés aux points de consommation. Les différents débits évoluent selon des chronogrammes respectifs qui font partie du scénario. 20 Les lois sont en particulier les lois de la physique applicables. Par exemple, dans un réseau hydraulique, la perte de charge (chute de pression) entre deux points du réseau en fonction de la vitesse d'écoulement et des données structurelles du réseau entre les deux points. 25 Les « paramètres de sortie primaires » sont des paramètres calculables par exécution du modèle direct pour obtenir leur valeur théorique, et qui sont en même temps mesurables grâce à l'instrumentation du réseau qui fournit une mesure de chacun des paramètres de sortie primaires, ou connus par exemple 30 d'après des statistiques antérieures. Dans le procédé selon l'invention, on compare la valeur théorique et la valeur réelle des paramètres de sortie primaires et on considère qu'il y a anomalie dans le réseau en cas de discordance entre les valeurs théoriques et mesurées. 35 3024254 4 Lorsqu'un paramètre de sortie primaire est « anormal » en ce sens qu'il y a entre ses valeurs théoriques et de mesure un écart dépassant une marge d'erreur acceptable, on passe à la résolution (exécution) d'un modèle inverse pour lequel le paramètre de sortie primaire anormal devient un paramètre 5 d'entrée secondaire auquel on donne comme valeur la valeur de mesure. Notamment pour que le modèle inverse ait les degrés de liberté mathématiques nécessaires, on introduit dans le modèle inverse au moins une variable supplémentaire constituée par un nouveau paramètre de sortie secondaire. De préférence on choisit judicieusement ces nouveaux paramètres de sortie 10 secondaires en fonction de l'anomalie de résultat de la comparaison avec les mesures suite à l'exécution du modèle direct. L'exécution du modèle inverse fournit des résultats qui sont soit directement informatifs, par exemple les nouveaux paramètres de sortie ont une valeur 15 théorique qui permet d'identifier l'anomalie avec une certitude suffisante, soit permettent de converger vers une telle identification par voie itérative comme il sera exposé plus loin. Ainsi, l'invention exploite finement la cohérence d'ensemble des variables 20 (paramètres de fonctionnement) liées par les « lois », typiquement les équations de comportement du réseau, en particulier hydraulique. L'au moins un paramètre de sortie secondaire peut être un paramètre ajouté au modèle hydraulique, par rapport au modèle direct. 25 Par exemple, si l'on soupçonne une fuite en un point du réseau, on peut ajouter un débit de fuite en tant que paramètre de sortie secondaire que l'on cherche à déterminer à l'aide du modèle inverse. L'au moins un paramètre de sortie secondaire peut aussi correspondre à au 30 moins une donnée supprimée, par rapport au modèle direct. Par exemple si un paramètre d'entrée primaire indique une vanne en position ouverte, mais que l'anomalie révélée par l'écart entre la valeur théorique et la valeur mesurée pour au moins un paramètre d'entrée primaire suggère que cette vanne pourrait être fermée ou partiellement fermée, le paramètre d'entrée primaire 35 en question est choisi comme paramètre de sortie secondaire, et supprimé des paramètres d'entrée pour l'exécution du modèle inverse. 3024254 5 Dans un mode de réalisation préféré, on sélectionne l'au moins un paramètre de sortie secondaire en appliquant au moins un critère révélateur d'une vraisemblance que l'au moins un nouveau paramètre de sortie secondaire soit 5 en cause dans les écarts constatés entre mesures et valeurs théoriques. En reprenant l'exemple précédent, une valeur de débit mesuré qui est anormale par rapport à la valeur théorique en un point donné du réseau rend par exemple vraisemblable soit la fermeture d'une vanne qui était considérée 10 ouverte, soit l'existence d'une fuite importante, soit un débit à contre-courant anormal en aval (liste non-limitative). Selon une particularité très préférentielle de l'invention, certaines au moins des données opérationnelles, des valeurs théoriques et des mesures comprennent 15 des séries temporelles, c'est-à-dire des séries de valeurs associées chacune à un horaire. Ainsi, le modèle direct et le modèle inverse sont opérés non seulement de façon statique, mais également de façon dynamique. Ceci permet au modèle de faire 20 des contrôles de cohérence supplémentaires. Par exemple le modèle peut comporter une loi qui relie la variation du niveau d'un réservoir à une somme de débits mesurés en différents points du réseau. En outre, en cas d'anomalie, on obtient non seulement la valeur de l'écart à un instant du fonctionnement, mais une fonction d'écart, donnant l'écart en fonction de tout ou partie des 25 paramètres primaires (exécution du modèle direct) ou secondaires (exécution du modèle secondaire), qu'ils soient d'entrée ou de sortie. Les mesures, qu'elles soient utilisées pour obtenir certaines données fonctionnelles, ou qu'elles concernent les paramètres de sortie pour obtenir les 30 valeurs de mesure à comparer aux valeurs théoriques, sont alors également obtenues sous forme de séries temporelles. Au lieu de les utiliser à l'état brut, on leur fait de préférence subir un lissage avant de les comparer aux valeurs théoriques. Ainsi les fonctions d'écart sont elles-mêmes automatiquement lissées. 35 Dans une étape d'initialisation au début de la mise en oeuvre d'un scénario, certaines données fonctionnelles peuvent être difficilement accessibles. C'est 3024254 6 typiquement le cas des consommations. Par exemple, s'agissant d'un réseau de distribution d'eau, il est difficile de connaître la consommation heure par heure alimentée par le réseau ou par chaque sous-réseau de ce réseau. On prévoit selon l'invention de remplacer de telles données difficilement accessibles en 5 temps réel par des données connues historiquement ou par expérience. Par exemple, on connaît par des études antérieures la consommation estimée d'un sous-réseau donné à chaque heure du jour et de la nuit, assortie d'information sur l'incertitude de ces valeurs. On peut, au moins au stade de l'initialisation, utiliser ces séries temporelles comme données fonctionnelles. 10 Certaines données structurelles ou fonctionnelles peuvent être inconnues ou incertaines. Par exemple, dans un réseau de distribution d'eau, l'altitude des points de mesure de pression est importante, puisque chaque mètre d'altitude supplémentaire correspond à une chute de pression de 10 kPa environ, qu'il est 15 nécessaire de prendre en compte pour la justesse des résultats obtenus avec les équations régissant la dynamique des écoulements. Des diamètres de canalisation, des rugosités de canalisation, et autres paramètres influant sur les pertes de charge dans l'écoulement, peuvent être mal connus, notamment dans des installations anciennes. 20 Notamment pour un tel cas, il peut être judicieux, au cours d'une étape d'initialisation, d'utiliser un modèle inverse pour déterminer certaines données structurelles ou opérationnelles impossibles ou difficiles à obtenir directement, le modèle inverse étant exécuté en tenant pour acquises les valeurs de certains 25 paramètres de sortie primaires, d'après les valeurs connues ou mesurées pour ces paramètres. Suivant un aspect important de l'invention, on intègre dans l'exécution du modèle direct des calculs de sensibilité de certains au moins des paramètres de 30 sortie primaires, aux variations de certains au moins des autres paramètres, et pour le modèle inverse on choisit comme paramètres de sortie secondaires, des paramètres d'entrée primaires auxquels le paramètre de sortie primaire présentant l'anomalie est plus particulièrement sensible. 35 La sensibilité peut être obtenue théoriquement comme la dérivée d'un paramètre par rapport à chacun des autres, connaissant les lois qui relient ces 3024254 7 paramètres entre eux, ou encore expérimentalement en comparant les variations de la mesure du paramètre de sortie avec les variations des autres paramètres pour toutes les valeurs de la série. La sensibilité peut être différente dans différentes plages de valeurs du paramètre dont la sensibilité 5 est évaluée, et dans différentes plages de valeurs du paramètre par rapport auquel la sensibilité est évaluée. La sensibilité peut dès lors s'exprimer comme une fonction ou comme un tableau de valeurs. La sensibilité peut être positive ou négative, selon que le paramètre de sortie évolue dans le sens, ou respectivement en sens inverse du paramètre par rapport auquel la sensibilité 10 du paramètre de sortie est calculée. Dans un mode de mise en oeuvre préféré : - par le modèle direct, on obtient une fonction d'écart de la valeur de l'écart entre valeur théorique et mesure pour le paramètre de sortie 15 primaire anormal en fonction des paramètres de sortie secondaires ; - on calcule d'après les sensibilités le gradient de la fonction d'écart par rapport aux paramètres ; - par l'application du modèle inverse on détermine un ensemble de nouvelles données, modifiant les données initiales, et pour lesquelles 20 le gradient global de la fonction d'écart présente une valeur nulle, indicative de sensibilités positives et négatives qui se neutralisent par rapport aux paramètres ; - on sélectionne parmi ces paramètres un sous-ensemble formé de ceux pour lesquels la sensibilité est de même signe ; 25 - on recommence le procédé de façon itérative en restreignant les paramètres de sortie secondaires à ceux faisant partie du sous-ensemble ; - et ainsi de suite jusqu'à localisation de l'anomalie à un seul paramètre ou un petit nombre de paramètres. 30 Pour constituer le sous-ensemble, on choisit de préférence le signe de sensibilité correspondant aux paramètres les plus vraisemblablement en cause, d'après leur nature, la valeur réaliste qui résulterait d'une variation dans le sens ramenant la valeur théorique anormale vers la valeur de mesure 35 correspondante. 3024254 8 Typiquement, lorsqu'un seul paramètre de sortie primaire présente une anomalie de la valeur mesurée par rapport à la valeur théorique, on privilégie l'hypothèse d'une erreur de mesure, surtout si le paramètre présentant l'anomalie n'est relié aux limites topographiques du réseau qu'à travers 5 d'autres points de mesure. Dans ce cas, on peut commencer par vérifier l'appareil de mesure avant application du modèle inverse. Mais on peut appliquer le modèle inverse de façon itérative pour conclure à l'erreur de mesure en constatant qu'aucune combinaison de valeurs n'est plausible pour les autres paramètres si la mesure anormale est prise comme hypothèse. Si la 10 mesure provient d'un point de mesure en limite de réseau, l'anomalie a plus de chances de provenir d'un élément exceptionnel au-delà des limites topographiques du réseau. Par exemple un consommateur injecte de l'eau de puits dans le réseau, au lieu de soutirer de l'eau fournie par le réseau. Il faut peut-être alors remettre en cause les données opérationnelles sur les points de 15 consommation. Dans un mode de mise en oeuvre élaboré, on exécute cycliquement le modèle direct et la comparaison des valeurs théoriques et mesurées, et l'exécution du modèle inverse est automatiquement activée lorsqu'au moins deux paramètres 20 de sortie primaires présentent une anomalie de la valeur mesurée par rapport à la valeur théorique. De préférence, les éléments de modèle comprennent en outre des indications d'incertitudes sur les données et les paramètres, et au cours de la comparaison 25 ces incertitudes sont prises en compte pour considérer une mesure de paramètre de sortie comme étant conforme ou au contraire anormale. DESCRIPTION DE MODES DE RÉALISATION 30 D'autres particularités et avantages de l'invention ressortiront encore de la description ci-après, relative à des exemples non-limitatifs relatifs à un réseau de distribution d'eau potable, également appelé « réseau hydraulique ». 35 3024254 9 Observation préliminaire La description qui suit vaut description de toute particularité qu'elle contient, qu'elle soit prise isolément des autres particularités même faisant partie du 5 même paragraphe ou de la même phrase, et vaut description de toute combinaison de telles particularités, dès lors qu'une telle particularité ou combinaison de particularités est distinctive de l'état de la technique et offre un effet technique, que ce soit seule ou en combinaison avec des concepts présentés ci-dessus, et que cette particularité soit exprimée dans les termes 10 même de cette description ou dans des termes plus ou moins généralisés. Vocabulaire On utilisera le vocabulaire suivant : 15 Série temporelle : séquence finie de données scalaires indexées par le temps, généralement espacées d'une durée constante. Modèle hydraulique : sans autre précision, désigne le modèle hydraulique d'un 20 réseau de transport et/ou de distribution d'eau potable : ensemble de données décrivant la topologie (graphe du réseau), la topographie (altitude), les caractéristiques des tronçons définissant leur comportement hydraulique (longueur, diamètre intérieur, matériau et/ou rugosité, coefficient de perte de charge singulière, singularités, etc.), les caractéristiques des équipements 25 définissant leur comportement hydraulique (diamètre, coefficient de perte de charge, valeur de consigne, courbe de pompe, etc.) consommations nominales associées aux noeuds du graphe, ainsi qu'un scénario décrivant les conditions de fonctionnement. 30 Scénario : Ensemble des données conjoncturelles ou opérationnelles décrivant les conditions aux limites utilisées avec un modèle hydraulique pour simuler le comportement d'un système hydraulique dans des conditions particulières, telles que les niveaux initiaux des réservoirs, les séries temporelles de modulation des consommations par secteur et/ou par type de consommateur, 35 les états des pompes, vannes, consignes d'équipements actifs pilotés et/ou les règles d'asservissement. 3024254 10 Mesure : série temporelle issue d'un système d'acquisition de données. Une mesure est associée à un composant élémentaire particulier du modèle hydraulique (arc ou tronçon pour une mesure de débit, noeud pour une mesure 5 de pression, réservoir pour une mesure de niveau). Modèle direct : La simulation du modèle hydraulique consiste à résoudre numériquement les équations de l'écoulement stationnaire décrites et paramétrées par le modèle hydraulique, et son évolution dynamique, en 10 utilisant les conditions aux limites décrites dans le scénario. Les résultats sont des débits, vitesses, perte de charge, états, pressions en tout point du système et à tout instant, ces résultats étant habituellement présentés sous forme de série temporelle ou sous forme de carte. 15 Modèle inverse : (résolution du problème inverse) consiste à résoudre numériquement les équations de l'écoulement stationnaire et leur évolution dans le temps. Elles sont décrites et paramétrées par le modèle hydraulique et un scénario, par rapport à un ou plusieurs paramètres du modèle ou du scénario, en utilisant les mesures de certains débits (ou vitesse), et/ou 20 pressions et/ou niveaux de réservoir. Données de télé-relève : série d'index de consommation pour un compteur, mesurés à une périodicité de 1 à 24 heures, et télé transmise en général au moins une fois par jour. 25 Caractéristique : valeur scalaire ou série temporelle contenant une information significative pour le processus étudié. Généralités 30 Le mode de mise en oeuvre du procédé selon l'invention utilise conjointement des données opérationnelles issues d'un système de télésurveillance, permettant la construction d'un scénario, un modèle hydraulique du système du réseau hydraulique, et un système d'analyse des résultats de l'exécution du 35 modèle. 3024254 11 Le programme de modélisation hydraulique intègre un calcul explicite des sensibilités des résultats (mesures) par rapport aux données d'entrée (consommations aux noeuds et résistances hydraulique aux arcs), permettant le calcul du gradient de l'erreur de simulation par rapport aux paramètres, puis 5 la résolution du problème inverse. L'utilisation du modèle inverse permet aussi de pallier l'absence de certaines mesures nécessaires à la construction d'un scénario. 10 L'utilisation du modèle direct permet de détecter des anomalies de comportement du système. L'utilisation du modèle inverse permet de localiser, dans l'espace et éventuellement dans le temps, l'origine possible de l'anomalie. 15 L'utilisation optionnelle de données de télé relève permet d'actualiser la répartition spatiale des consommations décrites dans le modèle hydraulique. Description détaillée 20 Les éléments qui sont fournis en entrée pour l'exécution du modèle, la comparaison avec les mesures des paramètres de sortie primaires, et le cas échéant l'exécution du modèle inverse comprennent : - les données structurelles du réseau, qui peuvent comprendre des consignes et des données d'état des équipements ; 25 - les paramètres sélectionnés pour constituer les paramètres d'entrée primaires, ainsi que les séries temporelles de leurs valeurs scalaires, obtenues par mesure ou par des données de télé-relève disponibles pour tout ou partie des abonnés, ou encore par des statistiques existantes ; - les paramètres sélectionnés pour constituer les paramètres de sortie 30 primaires ainsi que les séries temporelles des mesures de ces paramètres ; - des informations sur la précision et les incertitudes des mesures et des paramètres, leurs pondérations relatives, l'activation des paramètres sélectionnés pour être des paramètres d'entrée secondaires et des paramètres de sortie secondaires lors de la résolution du problème inverse. 35 3024254 12 Opérations élémentaires - Construction du scénario à partir des mesures et des séries d'état 5 Cette opération consiste à transcrire les mesures dans une forme utilisable par le logiciel de modélisation hydraulique. Par exemple, la somme algébrique des séries temporelles de débits d'entrée/sortie d'un secteur hydraulique est convertie en série de modulation des consommations nominales du secteur. 10 - Résolution du problème direct, calcul des caractéristiques Cette opération consiste à lancer le moteur de simulation hydraulique, avec comme données d'entrée le modèle hydraulique du système et le scénario créé par l'opération précédente. Les sorties exploitées sont les séries temporelles 15 simulées (valeurs théoriques des paramètres de sortie primaires) comparables à des mesures disponibles pour ces paramètres. - Calcul des caractéristiques 20 Pour chaque point de mesure (paramètre de sortie primaire), on calcule la caractéristique sur la base de la série des résidus, définie comme la différence entre valeur simulée (valeur théorique) et valeur mesurée (lissée ou non), calculée au pas d'acquisition de la mesure. La caractéristique indique s'il y a ou non anomalie. 25 - Résolution du problème inverse (cas général) La résolution du problème inverse est une opération itérative qui consiste à minimiser à l'aide d'un algorithme d'optimisation, une fonction d'écart entre 30 mesure et simulation (valeur théorique) au point de mesure (algorithme opérant par la méthode de la somme des carrés ou la méthode des valeurs absolues des séries de résidus), par rapport à un ensemble de paramètres configurable. Dans la fonction d'écart, chaque point de chaque série est pondéré par un coefficient, calculé à partir de l'incertitude a priori et de la 35 variance empirique locale des mesures. 3024254 13 Le calcul se fait en utilisant des mesures horodatées dans une fenêtre temporelle prédéterminée (par exemple, de 4 à 24 heures). Chaque itération nécessite le calcul de la fonction d'écart en s'appuyant sur le 5 calcul du problème direct. La méthode d'optimisation utilisée fait appel au calcul du gradient de la fonction d'écart par rapport aux paramètres. La valeur de ce gradient est calculée explicitement à partir des sensibilités (dérivées) des mesures par rapport aux paramètres du modèle hydraulique, qui est exécuté explicitement de manière imbriquée avec la résolution du problème direct. 10 L'un des aspects de l'innovation est d'une part le calcul explicite des sensibilités en même temps que la résolution du problème direct, et leur utilisation pour calculer le gradient ; d'autre part la possibilité de choisir par rapport à quels paramètres le problème inverse est résolu, par exemple : 15 - Données structurelles (caractéristiques hydrauliques des tronçons et équipements telles que la rugosité ou les coefficients de pertes de charge singulières (PCS)) - Côte (altitude) des capteurs de pression et/ou de niveau - Répartition spatiale des consommations 20 - Modulation temporelle d'un ou plusieurs types de consommation de répartition spatiale prédéfinie. Un autre élément important de l'innovation réside dans l'utilisation des sensibilités pour localiser l'anomalie, lorsque le paramètre identifié est associé à 25 plusieurs composants du modèle. Dans ce cas, lorsque la solution est obtenue, les sensibilités prennent des valeurs positives pour certains composants, négatives pour d'autres, car le gradient d'ensemble prend une valeur nulle à l'optimum. Il suffit alors de répéter la méthode en résolvant par rapport à l'un des deux sous-ensembles de composants dont le gradient est de même signe 30 pour déterminer la localisation de l'anomalie. La séquence ainsi formée de sous-ensembles emboîtés de composants localise l'anomalie de plus en plus précisément. La précision de la localisation est d'autant plus importante que la densité de 35 mesures est grande, et leur place judicieuse pour la détermination du type de paramètre sélectionné. 3024254 14 - Problème inverse appliqué à la reconstitution d'une mesure essentielle au cours d'une étape d'initialisation 5 Lorsqu'une mesure essentielle à la constitution d'un scénario (par exemple un débit d'entrée ou de sortie d'un secteur) est indisponible, la résolution du modèle inverse permet de reconstituer cette mesure manquante, dès lors que certaines conditions d'observabilité sont réunies, en l'occurrence que d'autres mesures existent dont les valeurs sont en relation avec la mesure manquante 10 (par exemple des mesures de niveau de réservoir et/ou de pression). Un scénario initial est alors constitué en utilisant une série prédéterminée, et le problème inverse résolu par rapport au paramètre constitué par la série manquante. 15 - Exécution du modèle inverse appliquée à la détection d'anomalies Appliquée à l'identification de données structurelles telles que rugosités ou PCS (perte de charge singulière), la résolution du problème inverse (exécution du modèle inverse) permet de caler le modèle hydraulique, de repérer des 20 anomalies de résistance hydraulique telles que des erreurs de saisie de diamètre ou de longueur de conduites, de localiser des vannes fermées oubliées après une intervention. Appliquée à l'identification de paramètres structurels tels que les altitudes des 25 capteurs ou les consignes d'équipements, la résolution du problème inverse permet de diagnostiquer des pertes de représentativité du modèle hydraulique après une intervention qui a entraîné la modification de l'un des paramètres du réseau. 30 Appliquée à l'identification de la répartition spatiale des consommations ou la série des coefficients de modulation d'un ou plusieurs types de consommateurs, la résolution du problème inverse permet de localiser les anomalies de consommation. L'écart de chaque paramètre entre ses valeurs théoriques et de mesure constitue un indicateur quantitatif ou une caractéristique utilisable pour 35 la détection d'anomalie. 3024254 15 - Interprétation des caractéristiques et anomalies Indicateurs de conformité : pour chaque caractéristique, on calcule un ou plusieurs indicateurs quantitatifs : critère de Nash, « Index of Agreement », 5 taux de dépassement d'un seuil d'erreur. Anomalies : pour chaque mesure on calcule un indicateur qualitatif (bonne / assez bonne / médiocre / mauvaise) à partir de la position de ses indicateurs de conformité par rapport à des seuils fixés opérationnellement. Les valeurs 10 classées « médiocre » ou « mauvaise » sont identifiées comme anomalies. Une anomalie isolée sur un secteur est classée « anomalie de mesure » ; l'observation de plusieurs anomalies sur un secteur déclenche une « anomalie réseau ». 15 Avantages L'ensemble constitué par les composants exposés plus haut peut être connecté au système d'information technique de l'opérateur d'un système d'approvisionnement d'eau potable. L'ensemble est construit de telle manière 20 qu'il peut être activé régulièrement, à une fréquence comprise entre la fréquence d'acquisition des données mesurées, et une fois par jour. Dans ces conditions la détection et la caractérisation de l'anomalie sont beaucoup plus rapides qu'avec les méthodes habituellement utilisées et l'efficacité opérationnelle en est considérablement améliorée. 25 Dans la phase de mise en place d'un modèle hydraulique pour un réseau, l'invention permet de réduire fortement le temps de calage du modèle. Premier exemple de réalisation 30 L'invention est configurée pour modéliser un secteur hydraulique alimenté par gravité par deux réservoirs de tête A et B. Le modèle hydraulique est constitué d'un fichier au format 35 PICCOLO/GANESSA_SI M 3024254 16 Le générateur de scénario utilisé est GANESSA_GS. Le moteur de simulation hydraulique incorporant la résolution de problème (exécution de modèle) direct ou inverse est GANESSA_SIM. 5 Les éléments alimentant le modèle sont le débit d'entrée de chaque réservoir, le niveau de chaque réservoir, les débits de sortie de chaque réservoir et la pression en un point éloigné, à la fréquence d'une fois par jour. Le débitmètre de sortie du réservoir A présente une anomalie consistant à sous-compter d'environ 10%. 10 La modélisation hydraulique directe sur 24 heures produit, entre autres, le cumul de la somme algébrique des débits d'entrée et de sortie des réservoirs. La caractéristique « niveau calculé moins niveau mesuré » renseigne sur la justesse des mesures de débit, un écart positif en fin de période indiquant un 15 sous-comptage des débits de sortie, ou un sur comptage des débits d'entrée. L'indicateur calculé à partir de cette caractéristique met en évidence une anomalie ; le signe de la caractéristique pour chacune des mesures de débit de sortie indique le sens de l'anomalie. 20 Le modèle inverse configuré avec comme paramètres de sortie secondaires, une consommation supplémentaire fixe au point de sortie de chaque réservoir et une série temporelle de coefficients, calcule l'ajustement des deux séries d'intensité pour chaque point de sortie de réservoir. La série associée au point de sortie du réservoir A présente un profil proportionnel au débit de sortie de A, 25 correspondant à une consommation de l'ordre de 10%, tandis que la série associée à B présente des intensités négligeables. Second exemple 30 L'invention est configurée pour modéliser un secteur hydraulique alimenté par une conduite maîtresse en deux points A et B. Le modèle hydraulique est constitué d'un fichier au format EPANET. Le générateur de scénario utilisé est GANESSA_GS. Le moteur de simulation hydraulique incorporant la résolution du problème 35 inverse est GANEssA_SIM. 3024254 17 Les mesures sont les débits et pression d'entrée (points A et B), et les mesures de pression obtenues à six points additionnels répartis sur le réseau de distribution. 5 Dans un premier temps l'invention a été utilisée hors ligne pour ajuster le calage hydraulique du modèle : le modèle inverse a permis d'identifier l'altitude précise des points de mesure de pression, la résistance hydraulique des deux vannes d'entrée, et la rugosité moyenne des 4 principales classes de matériaux dont sont constituées les conduites du secteur. 10 Dans un second temps l'invention est mise en oeuvre à une fréquence quotidienne. Une faible consommation additionnelle (par exemple 1 1/s) est répartie sur l'ensemble des noeuds du modèle, par exemple au prorata des longueurs de conduites qui leur sont raccordées, et attribuée à un code 15 consommateur particulier (par exemple 'LEAK'), et une série temporelle est associée à ce code 'LEAK', initialisée avec des valeurs nulles ou unitaires. Le problème inverse est alors résolu par rapport aux coefficients de la série temporelle des coefficients de 'LEAK'. La caractéristique correspondante est alors analysée : si la consommation correspondante est inférieure à un seuil 20 prédéterminé, alors la situation est considérée comme normale. Dans le cas contraire, une anomalie de consommation est présumée. L'invention est alors à nouveau mise en oeuvre après avoir réalloué la consommation additionnelle aux seuls noeuds pour lesquels la sensibilité calculée à l'étape précédente est négative, et ce à plusieurs reprises, jusqu'à ce que le nombre de noeuds 25 restant soit inférieur à un seuil prédéterminé. Les derniers ensembles de noeuds correspondants sont alors reportés comme localisation de l'anomalie, et la valeur moyenne de la consommation additionnelle (associée à 'LEAK') comme son intensité. 30 Bien entendu l'invention n'est pas limitée aux exemples décrits. REVENDICATIONS1. Procédé de détection d'anomalies dans un réseau de distribution, comprenant les étapes suivantes : établir pour le réseau de distribution un modèle direct comprenant les éléments suivants o des données structurelles du réseau de distribution relatives à la topographie dudit réseau de distribution et comprenant des consignes ou des données d'états d'équipement dudit réseau de distribution, o des données opérationnelles relatives à des paramètres choisis comme paramètres d'entrée primaires décrivant un scénario opérationnel, o des lois reliant lesdites données structurelles, les paramètres d'entrée primaires et des paramètres de sortie primaires - exécuter le modèle direct par un moteur de simulation hydraulique pour déterminer des valeurs théoriques des paramètres de sortie primaires avant ou après tout ou partie de ce qui précède, obtenir des mesures issues de systèmes d'acquisition de données des paramètres de sortie primaires ; - comparer les valeurs théoriques avec les mesures ; - dans des cas où la comparaison révèle au moins un paramètre de sortie primaire anormal du fait eun écart important entre valeur théorique et mesure, établir un modèle inverse qui est dérivé du modèle direct en introduisant l'au moins un paramètre de sortie primaire anormal comme un paramètre d'entrée secondaire du modèle inverse, le modèle inverse comportant de plus au moins un paramètre de sortie secondaire en plus de ceux des paramètres de sortie primaires qui sont conservés - exécuter par un moteur de simulation hydraulique le modèle inverse ; - localiser une anomalie du réseau de distribution à partir de l'exécution du modèle inverse. 2. Procédé selon la revendication 1, caractérisé en ce que au moins un paramètre de sortie secondairè est un paramètre ajouté au modèle, par rapport au modèle direct. Procé 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que au moins un paramètre de sortie secondaire correspond à au moins une donnée supprimée, par rapport au modèle direct. 4. Procédé selon l'une des revendications 1 à 3, caractérisé en ce qu'on sélectionne au moins un paramètre de sortie secondaire en appliquant des critères indiquant une vraisemblance que au moins un nouveau paramètre de sortie secondaire soit en cause dans les écarts constatés entre mesures et valeurs théoriques. 5. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que certaines au moins des données opérationnelles, des valeurs théoriques et des mesures comprennent des séries temporelles c'est-à-dire chacune une série de valeurs associées chacune à un horaire. 6. Procédé selon la revendication 5, caractérisé en ce que les séries temporelles formées de valeurs obtenues à partir de mesures sont lissées suivant une durée caractéristique configurable. 7. Procédé selon l'une des revendications 1 à 6, caractérisé en ce que dans une étape d'initialisation on utilise un modèle inverse pour déterminer certaines données structurelles ou opérationnelles impossibles ou difficiles à obtenir directement, le modèle inverse étant exécuté en tenant pour acquises les valeurs de certains paramètres de sortie primaires. 8. Procédé selon l'une des revendications 1 à 7, caractérisé en ce qu'on intègre dans l'exécution du modèle direct des calculs de sensibilité de certains au moins des paramètres de sortie primaires aux variations de certains au moins des autres paramètres, et en ce que pour le modèle inverse on choisit comme paramètres de sortie secondaires, des paramètres d'entrée primaires auxquels le paramètre de sortie primaire anormal est plus particulièrement sensible. 9. Procédé selon la revendication 8, caractérisé en ce que : - par le modèle direct, on obtient une fonction d'écart de la valeur de l'écart entre valeur théorique et mesure pour le paramètre de sortie primaire anormal, en fonction des paramètres de sortie secondaires ; - -aux calculs d'après les sensibilités le gradient de la fonction d'écart par rapport aux paramètres ; par l'application du modèle inverse on détermine un ensemble de nouvelles données, modifiant les données initiales, et pour lesquels le gradient global de la fonction d'écart présente une valeur nulle, indicative de sensibilités positives et négatives qui se neutralisent ; - on sélectionne parmi les paramètres par rapport auxquels les sensibilités sont calculées, un sous-ensemble formé de ceux pour lesquels la sensibilité est de même signe ; - on recommence le procédé de façon itérative en restreignant les paramètres de sortie secondaires à ceux faisant partie du sous- ensemble ; et ainsi de suite jusqu'à localisation de l'anomalie à un seul paramètre ou un petit nombre de paramètres. 10. Procédé selon l'une des revendications 1 à 6, caractérisé en ce que lorsqu'un seul paramètre de sortie primaire présente une anomalie de la valeur mesurée par rapport à la valeur théorique, on privilégie l'hypothèse d'une erreur de mesure. 11. Procédé selon l'une des revendications 1 à 10, caractérisé en ce qu'on exécute cycliquement le modèle direct et la comparaison des valeurs théoriques et mesurées, et en ce que l'exécution du modèle inverse est automatiquement activée lorsqu'au mois deux paramètres de sortie primaires présentent une anomalie de la valeur mesurée par rapport à la valeur théorique. 12.Procédé selon l'une des revendications 1 à 11, caractérisé en ce que les éléments de modèle comprennent en outre des indications d'incertitudes -21- 3024254 sur les données et les paramètres, et au cours de la comparaison ces incertitudes sont prises en compte pour considérer une mesure de paramètre de sortie comme étant conforme ou au contraire anormale. 5 13.Utilisation d'un procédé selon l'une quelconque des revendications 1 à 12, comprenant en outre une amélioration d'un niveau de performance d'un réseau de distribution.