FR 3070495 A1 20190301 FR 1757919 A 20170828 Détection de défaillance d'un réseau de distribution électrique L'invention relève du domaine de la gestion d'un réseau de distribution électrique. H est de la responsabilité d'un gestionnaire d'un réseau de distribution électrique de surveiller le bon fonctionnement du réseau. Des mesures préventives, notamment l'entretien, permettent de limiter les risques de dysfonctionnement. Des mesures curatives, notamment les interventions d'opérateurs sur le terrain, permettent de réduire la durée et l' ampleur des disfonctionnements lorsqu'ils surviennent. Pour permettre des interventions rapides et efficaces, il convient de détecter la survenance d'un incident au plus tôt et de le qualifier le plus précisément possible. Les réseaux de transport et de répartition de l'énergie électrique correspondent à la partie amont de l'infrastructure électrique et acheminent l'énergie électrique sur de longues distances. De tels réseaux sont généralement équipés de systèmes de surveillance au moins en partie automatisés. Les réseaux de transport et de répartition comprennent notamment des capteurs aptes à transmettre à distance des données de surveillance de l'état du réseau. Les réseaux de distribution correspondent à la partie avale de l'infrastructure et acheminent l'énergie électrique sur de courtes distances mais jusque chaque noeud de consommation. Les réseaux de distribution sont donc plus denses et nombreux que les réseaux de transport et de répartition. Par conséquent, équiper le réseau de distribution de systèmes de surveillance et en assurer la maintenance présenteraient un coût élevé. Les réseaux de distribution sont généralement dépourvus de tels systèmes de surveillance. En pratique, un incident sur un réseau de distribution est généralement signalé au gestionnaire du réseau par un ou plusieurs utilisateurs constatant l'interruption de la fourniture d'énergie électrique. De tels signalements, généralement effectués par téléphone, sont peu précis et peu réactifs. Par exemple, il est difficile de distinguer une anomalie sur le réseau de distribution de la responsabilité du gestionnaire d'une anomalie sur le réseau privé (par exemple domestique) de la responsabilité de l'utilisateur lui-même. Il en résulte un nombre important de fausses alertes à traiter. Depuis peu, les noeuds consommateurs des réseaux de distribution, parfois appelé « Point De Livraison » (PDL), peuvent être équipés de dispositifs communicants. Par exemple, des compteurs communicants, connus sous la référence commerciale « Linky » sont en cours de 35 déploiement en France. De tels compteurs comprennent un module de communication apte à 2 transmettre des données à distance de manière au moins en partie automatisée. De tels modules sont agencés pour transmettre des informations relatives à la consommation en aval du noeud consommateur correspondant, notamment dans un but de facturation. De tels modules ne sont pas prévus pour transmettre des informations relatives au réseau de distribution en amont. Les 5 modules de communication des noeuds consommateurs ne permettent donc pas, en l'état, de détecter des anomalies sur le réseau de distribution. L'invention vient améliorer la situation. 10 H est proposé un procédé de détection d'anomalie sur un réseau, mis en oeuvre par ordinateur, incluant une pluralité de noeuds consommateurs, chacun des noeuds consommateurs étant apte à transmettre une confirmation via le réseau de distribution électrique en réponse à une requête reçue par l'intermédiaire du réseau de distribution électrique. Le procédé comprend : a) collecter les retours de requêtes transmises à chacun des noeuds consommateurs via le réseau 15 de distribution électrique à une pluralité d'instants successifs d'une période d'observation et à une pluralité d'instants successifs d'une période courante, b) pour chaque requête, établir une variable d'état de type booléen dont la valeur correspond respectivement à la réception d'une confirmation ou à l'absence de réception d'une confirmation, 20 c) pour chaque noeud consommateur et pour chacune des périodes, déduire une distribution statistique des transitions d'états possibles entre deux instants successifs, d) pour chaque noeud consommateur et pour chacune des périodes, calculer une vraisemblance de la distribution statistique, e) pour un ensemble de plusieurs noeuds consommateurs, calculer un indice d'écart à la 25 normalité en fonction des vraisemblances de la période d'observation et de la période courante des noeuds consommateurs, f) comparer l'indice obtenu à une valeur seuil prédéfinie, g) uniquement dans le cas où l'indice obtenu dépasse la valeur seuil, générer un signal d'alerte. 30 Un tel procédé permet de détecter avec un haut degré de fiabilité des anomalies sur le réseau de de distribution électrique. Le compromis entre une détection d'un maximum d'anomalies et un minimum de fausses alertes est amélioré. En particulier, le procédé permet de distinguer un disfonctionnement réel du réseau de distribution d'un défaut de communication des noeuds consommateurs. Ainsi, la surveillance du réseau peut être en grande partie automatisée et réalisée sensiblement en continu. En cas de disfonctionnement, l'anomalie peut être détectée en 3 temps réel et des mesures curatives peuvent être mises en oeuvre au plus tôt. Les désagréments causés par un disfonctionnement sont de courte durée. Selon un autre aspect, il est proposé un module de détection d'anomalie sur un réseau de 5 distribution électrique comprenant une unité de traitement, le réseau incluant une pluralité de noeuds consommateurs, l'unité de traitement étant reliée au réseau et chacun des noeuds consommateurs étant apte à transmettre une confirmation via le réseau de distribution électrique en réponse à une requête reçue par l'intermédiaire du réseau de distribution électrique. L'unité de traitement est agencée de manière à: 10 a) interroger chacun des noeuds consommateurs via le réseau de distribution électrique à une pluralité d'instants successifs d'une période d'observation et à une pluralité d'instants successifs d'une période courante, b) pour chaque interrogation, établir une variable d'état de type booléen dont la valeur correspond respectivement à la réception d'une confirmation ou à l'absence de réception d'une 15 confirmation, c) pour chaque noeud consommateur et pour chacune des périodes, déduire une distribution statistique des transitions d'états possibles entre deux instants successifs, d) pour chaque noeud consommateur et pour chacune des périodes, calculer une vraisemblance de la distribution statistique, 20 e) pour un ensemble de plusieurs noeuds consommateurs, calculer un indice d'écart à la normalité en fonction des vraisemblances de la période d'observation et de la période courante des noeuds consommateurs, f) comparer l'indice obtenu à une valeur seuil prédéfinie, g) uniquement dans le cas où l'indice obtenu dépasse la valeur seuil, générer (113) un signal 25 d'alerte. Selon un autre aspect, il est proposé un programme d'ordinateur apte à être mis en oeuvre en un procédé de détection d'anomalie sur un réseau tel que les procédés définis dans les présentes. Selon encore un autre aspect, il est proposé un support d'enregistrement non transitoire lisible 30 par un ordinateur sur lequel est enregistré un programme comprenant des instructions pour la mise en oeuvre d'un des procédés proposés dans les présentes. Les caractéristiques suivantes peuvent, optionnellement, être mises en oeuvre. Elles peuvent être mises en oeuvre indépendamment les unes des autres ou en combinaison les unes avec les 35 autres : 4 - L'indice d'écart à la normalité est équivalent à une somme ou une moyenne des différences de vraisemblances entre la période d'observation et la période courante de chaque noeud consommateur. Cela permet de détecter aisément et rapidement une anomalie touchant 5 simultanément un nombre important de noeuds consommateurs sans pour autant déclencher de fausses alertes en cas de disfonctionnement limité à un unique noeud (anomalie probablement indépendante du réseau lui-même) ou de disfonctionnement à durée très courte par rapport à la période d'observation (anomalie à très faible durée pour laquelle déclencher une alerte est probablement inutile). 10 - Les interrogations et les réceptions de confirmations sont effectuées par Courtant Porteur en Ligne sur le réseau de distribution électrique. Ce type de communication ne nécessite ni d'installer un réseau de communication filaire dédié, ni de prévoir des communications sans fil. Les équipements physiques permettant la mise en oeuvre de la communication par CPL sont 15 parfois déjà installés ou en cours de déploiement. Ainsi, un réseau électrique existant peut être mis en conformité avec le procédé proposé ici par des mises à jour logicielles plus rapides et moins coûteuses que des installations physiques sur le terrain. - La période d'observation est sélectionnée parmi plusieurs en fonction des conditions de 20 fonctionnement du réseau de distribution électrique au cours de la période courante. Ainsi, les périodes d'observation peuvent être choisies analogues à la période courante de manière à limiter les biais contextuels dans les comparaisons. Par exemple, on peut considérer qu'il est préférable d'utiliser une période d'observation correspondant à une soirée d'hiver (consommation électrique élevée sur le réseau) pour comparer les comportements d'une période 25 courante correspondant elle aussi à une soirée d'hiver. - La période d'observation s'étend sur une durée comprise entre 1 seconde et 100 jours ; la période courante s'étend sur une durée comprise entre 10 millisecondes et 7 jours ; les instants successifs de l'une au moins de la période d'observation et de la période courante sont 30 mutuellement espacés d'une durée comprise entre 1 milliseconde et 24 heures ; et/ou l'une au moins de la période d'observation et de la période courante comprend entre 2 et 500 instants successifs. D'après les essais de la demanderesse, de tels ordres de grandeur permettent de détecter la plupart des anomalies tout en présentant un faible taux de faux positifs. 5 - Les opérations sont mises en oeuvre de manière répétée de sorte que les variable d'état relatives à la période courante sont mises à jour sensiblement en temps réel. Une surveillance en temps réel permet de générer un signal d'alerte rapidement après qu'une anomalie est survenue. Les délais d'interventions pour remédier aux anomalies s'en trouvent réduits. 5 - Les opérations sont mises en oeuvre distinctement pour plusieurs ensembles de noeuds consommateurs d'un même réseau. Ceci permet, en cas d'alerte, non seulement une rapidité de la détection d'anomalie mais aussi une rapidité d'identification de l'origine de l'anomalie L'efficacité des équipes d'intervention sur le terrain s'en trouve amélioré. En outre, les réponses 10 à apporter aux anomalies peuvent être adaptées en fonction des noeuds consommateurs concernés, par exemple en fonction des activités humaines correspondant aux noeuds consommateurs (réseaux électriques domestiques, industriels ou d'importance vitale telle que les hôpitaux). 15 - La composition de chaque ensemble est sélectionnée en fonction des positions des noeuds consommateurs sur le réseau de sorte qu'un signal d'alerte est associé à une zone prédéfinie du réseau. Ceci permet aussi, en cas d'alerte, d'identifier rapidement l'origine de l'anomalie D'autres caractéristiques, détails et avantages de l'invention apparaîtront à la lecture de la 20 description détaillée ci-après, et à l'analyse des dessins annexés, sur lesquels : - la figure 1 est une représentation schématique d'un réseau à surveiller ; - la figure 2 représente les données brutes résultant d'une surveillance d'un ensemble de noeuds pendant une période commune ; - la figure 3 est un diagramme représentant différentes opérations mises en oeuvre selon 25 l'invention et sous la forme d'une boucle itérative ; - la figure 4 représente les données brutes résultant d'une surveillance de deux noeuds pendant une période commune ; - la figure 5 représente les données brutes résultant d'une surveillance d'un ensemble de noeuds pendant une période commune réelle ; 30 - la figure 6 est une représentation analogue à celle de la figure 5 d'une surveillance pendant laquelle des anomalies sont simulées ; - la figure 7 représente l'évolution comparée d'un indice calculé selon deux modes de réalisation à partir des données de la surveillance représentée en figure 6 ; et - la figure 8 représente l'évolution du taux de détection des anomalies en fonction du nombre de 35 fausses alarmes selon quatre méthodes distinctes. 6 Les dessins et la description ci-après contiennent, pour l'essentiel, des éléments de caractère certain. Ils pourront donc non seulement servir à mieux faire comprendre la présente invention, mais aussi contribuer à sa définition, le cas échéant. 5 En figure 1, un réseau 1 de distribution d'électricité est représenté. Le réseau 1 comprend un noeud source 3 et une pluralité de noeuds consommateurs 5' l'indice i permettant de distinguer les noeuds consommateurs 5 les uns des autres. Dans l'exemple décrit ici, le réseau 1 comprend neuf noeuds 5, répartis en trois groupes de trois noeuds. 10 Dans l'exemple décrit ici, un noeud source 3 correspond à un transformateur disposé entre le réseau de transport moyenne tension (ou haute tension type A HTA) et le réseau de distribution basse tension (BT) en aval. Le réseau Basse Tension (BT) choisi en exemple présente une structure arborescente. Le sommet de la structure arborescente correspond au noeud source 3, 15 matérialisé ici par un poste de transformation HTA/BT. Chaque noeud consommateur 5 forme un point de branchement le long d'un des tronçons du réseau 1. Les noeuds consommateurs 5 correspondent sensiblement aux « Point De Livraison (PDL) ». Chaque noeud consommateur 5 forme une interface entre le réseau 1 et un réseau électrique 20 secondaire, généralement privé. En pratique, le noeud consommateur 5 est matérialisé par un ensemble d'équipements, dont un compteur propre à comptabiliser la consommation électrique du réseau secondaire à des fins de facturation. Par abus de langage, les noeuds consommateurs 5 peuvent être appelés compteurs. 25 Dans le présent contexte, des anomalies 7, 9 se produisent sur le réseau 1. L'anomalie 7 a lieu immédiatement en aval du noeud source 3 et en amont des bifurcations du réseau 1. L'anomalie 7 entraîne une interruption de la fourniture d'énergie de tous les noeuds 5 du réseau 1. L'anomalie 9 a lieu en aval d'une des bifurcations du réseau 1. L'anomalie 9 entraîne une interruption de la fourniture d'énergie de certains seulement des noeuds 5 du réseau 1, ici les 30 noeuds 54, 55 et 56. Dans les modes de réalisation décrit ci-après, le procédé permet de détecter une anomalie 7, 9 ayant pour conséquence l'interruption de la fourniture d'électricité de certains au moins des noeuds consommateurs 5 du réseau 1. Le procédé permet de détecter la survenance d'au moins 35 une anomalie 7, 9 sans nécessairement en identifier la source ou la partie du réseau 1 concernée. 7 Dans des modes de réalisation particuliers, le procédé permet non seulement de détecter la survenance d'une anomalie 7, 9 mais permet en outre d'en identifier certaines au moins des caractéristiques, par exemple les parties du réseau 1 concernées. Dans ce dernier cas, les précisions de localisation permettent d'adapter en conséquence les mesures curatives à mener et 5 donc de réduire encore les durées d'interruption de la fourniture d'énergie pour les noeuds consommateurs 5 impactés. Sur le réseau 1, chacun des noeuds consommateurs 5 est équipé d'un module de communication, ou interface de communication, apte à transmettre des données vers un concentrateur commun 10 Chacun des noeuds consommateurs 5 est apte à transmettre, au moins, une confirmation via le réseau 1 de distribution électrique en réponse à une requête reçue par l'intermédiaire du réseau 1 de distribution électrique. Chaque noeud consommateur 5 est communicant. Par exemple, les compteurs intelligents ou AMM pour « Automated Meter Management » sont aptes à échanger des données avec le concentrateur. Le concentrateur collecte les données de plusieurs AMM. Le 15 concentrateur est, par exemple, situé au niveau du noeud source 3. Dans l'exemple décrit ici, le concentrateur transmet à son tour des informations à une unité de traitement distante et commune à plusieurs concentrateurs, par exemple via un réseau de téléphonie mobile. Ainsi, l'unité de traitement de données comprend un processeur couplé de 20 manière opérationnelle avec une mémoire pour mettre en oeuvre au moins en partie un procédé tel que défini dans les présentes. En variante, une partie au moins du traitement des données peut être mis en oeuvre par les concentrateurs et/ou localement par des équipements voisins des concentrateurs. 25 Parmi les AMM développés, certains comprennent des modules de communication fonctionnant par Courant Porteur en Ligne (CPL). Ainsi, les données sont échangées par le réseau électrique lui-même sans qu'il soit nécessaire de prévoir un réseau de communication dédié entre le concentrateur et les noeuds consommateurs 5. Le réseau de communication utilise physiquement le réseau câblé électrique. 30 Parmi les protocoles de communication CPL, certains se composent d'un terminal maître (ici le concentrateur) et de plusieurs terminaux esclaves (ici les modules de communication des noeuds consommateurs 5). L'ensemble des terminaux, ou noeuds, est parfois appelé « grappe CPL ». Contrairement à d'autres configurations connues, une telle configuration maître-esclave ne 35 permet pas aux noeuds consommateurs 5 de transmettre activement un signal d'alerte au 8 concentrateur pour indiquer une perte d'alimentation. Les modules de communication des noeuds consommateurs 5 sont généralement alimentés en énergie par le réseau de distribution lui-même. Par conséquent, les noeuds consommateurs 5 cessent de fonctionner (et de transmettre des données) lorsque la fourniture d'énergie est interrompue. 5 Dans de premiers essais, la demanderesse a tenté d'interpréter l'absence de communication de certains noeuds consommateurs 5 détectée depuis le concentrateur comme un indicateur d'une anomalie Cette approche est représentée en figure 2. 10 Le concentrateur interroge chacun des noeuds consommateurs 5 par une requête. Lorsqu'il est en mesure de la faire, le noeud consommateur 5 répond par une confirmation. Pour chaque requête, le concentrateur reçoit donc une confirmation ou ne reçoit pas de confirmation. La réception d'une confirmation est notée « OK » ci-après et représentée par un remplissage grisé sur les figures 2 et 4 à 6. L'absence de confirmation est notée « KO » ci-après et représentée par un 15 remplissage noir sur les figures 2 et 4 à 6. Dans les exemples décrits ici, la requête et la confirmation attendue, à l'image d'une commande « ping », sont dépourvues d'autre information. Ainsi, la consommation énergétique d'un tel fonctionnement et les risques d'anomalie de communication sont réduits. En variante, l'opération d'interrogation par le concentrateur et de réponse par chaque noeud consommateur 5 peut être combinée à d'autres 20 opérations de communication. Par exemple, une requête du concentrateur peut viser à relever une consommation sur le noeud consommateur 5. La réponse du noeud consommateur 5, respectivement l'absence de réponse, peut être interprétée comme une confirmation OK, respectivement une absence de confirmation KO. 25 Dans l'essai représenté en figure 2, chacun des noeuds consommateurs 5 d'un ensemble (ou grappe) est interrogé par un concentrateur via un réseau 1 commun de distribution électrique. Les requêtes/interrogations sont effectuées, pour chacun des noeuds consommateurs 5, à intervalles réguliers. En figure 2, l'abscisse correspond au temps tandis que chacune des lignes de la partie supérieure correspond à un noeud consommateur 5. Les confirmations reçues (OK) 30 sont représentées en grisé. Les absences de confirmation (KO) sont représentées en noir. Les absences d'information sont représentées en blanc. Dans un premier essai, il a été choisi arbitrairement que l'absence de confirmation (KO) d'au moins 20 noeuds consommateurs 5 par intervalle entraîne l'émission d'une alerte 11. En partie 35 inférieure de la figure 2, les alertes 11 émises au cours du temps sont représentées par du noir, 9 tandis que l'absence d'alerte est représentée par du blanc. Durant cet essai, aucune anomalie d'alimentation n'a eu lieu sur le réseau 1 pendant la période étudiée. Toutes les alertes sont donc, ici, de fausses alertes. 5 Dans cet exemple, les absences de confirmation KO sont dues à d'autres raisons que des interruptions d'alimentation des noeuds consommateurs 5. Les absences de confirmation KO sont manifestement dues à des problèmes propres à la communication. En effet, le réseau 1 est, par nature, conçu pour transporter de l'énergie électrique et pas des données. La technologie du CPL est imparfaite : des erreurs de communication sont fréquentes. Par exemple, les charges 10 connectées au réseau de distribution peuvent provoquer des atténuations du signal CPL ou générer des bruits parasites dans les bandes de fréquences utilisées pour les émissions CPL. La qualité des communications CPL peut évoluer dans le temps en fonction des conditions d'utilisation du réseau électrique. De plus, de telles perturbations sont souvent géographiquement hétérogènes : les communications avec les noeuds consommateurs 5 sont 15 perturbées différemment en fonction de l'emplacement des noeuds consommateurs 5 sur le réseau 1. En outre, les modules de communication présents au niveau des noeuds consommateurs 5 sont prévus initialement pour assurer des services de relève (pour la facturation), pour programmer les grilles tarifaires et pour mettre à jour les composants logiciels. Pour assurer de telles fonctions, il est généralement admis qu'un taux de réponse du 20 module de communication de 90% est suffisant. Les modules de communication installés ou en cours de déploiement sont donc prévus pour une telle performance Par conséquent, en fonctionnement normal, il est beaucoup plus probable qu'une absence de confirmation KO résulte d'une erreur de communication plutôt que d'une interruption de la fourniture d'énergie du noeud consommateur 5. 25 Ainsi, la méthode décrite jusqu'ici ne donne pas entièrement satisfaction, pour au moins les trois raisons suivantes : - la probabilité d'un échec de communication est bien plus élevée que la probabilité d'une interruption de l'alimentation électrique ; 30 - la qualité des communications varie d'un noeud consommateur 5 à l'autre ; et - la qualité des communications des noeuds consommateurs5 varie dans le temps. Il est maintenant fait référence à la figure 3. 10 Chacun des noeuds consommateurs 5 est interrogé via le réseau 1 de distribution électrique. Les interrogations sont répétées à une pluralité d'instants t successifs d'une période d'observation T (ou d'apprentissage) et à une pluralité d'instants t' successifs d'une période courante T'. Dans l'exemple décrit ici, le concentrateur émet des requêtes à destination de chacun des noeuds 5 consommateurs 5. Les interrogations et les réceptions de confirmation OK sont effectuées via le réseau 1 lui-même, ici au moyen de CPL. Dans l'exemple décrit ici, une requête spécifique est envoyée à chaque noeud consommateur 5. Un tel fonctionnement est, par exemple, compatible avec le protocole de communication par CPL normé et connu sous la référence « G1 ». 10 En variante, une requête générale à tous les noeuds peut être transmise (type « broadcast »). Le retour de chaque noeud 5, c'est-à-dire la réponse ou l'absence de réponse, est ensuite répertoriée. L'état des noeuds (OK ou KO) peut aussi être déduit en fonction de l'existence d'une réponse ou de l'absence de réponse à tout type de requêtes, notamment des requêtes non spécifiquement dédiées à la mise en oeuvre du présent procédé. 15 Les retours de plusieurs noeuds peuvent être au moins en partie regroupés dans des ensembles de données stockés dans une mémoire. Un tel fonctionnement est, par exemple, compatible avec le protocole de communication par CPL normé et connu sous la référence «G3 ». Selon ce protocole, les retours des noeuds peuvent être transmis indirectement de noeuds en noeuds, par 20 exemple au moyen d'ensemble de données de type « table de voisins ». En réponse à chaque requête, une réponse/confirmation est reçue ou rien n'est reçu. Dans une première opération 101, les retours des requêtes sont collectés et répertoriés. À réception de la confirmation ou à non réception de la confirmation, une deuxième opération 103 est mise en 25 oeuvre : il est établi une variable d'état Bu de type booléen dont la valeur correspond respectivement à la réception d'une confirmation (OK) ou à l'absence de réception d'une confirmation (KO). La période d'observation T est choisie comme une référence. La période d'observation T est 30 sélectionnée pour être représentative du fonctionnement normal des communications entre les noeuds consommateurs 5 et le concentrateur en l'absence d'anomalie d'alimentation sur le réseau 1. La période courante T' correspond à la période pour laquelle on souhaite identifier l'existence 35 éventuelle d'anomalies. Lorsque le procédé est mis en oeuvre dans le cadre d'une surveillance 11 en temps réel du réseau 1, la période courante T' correspond au temps réel de la surveillance tandis que la période d'observation T correspond à une période antérieure pour laquelle les résultats de la deuxième opération 103 ont été préalablement enregistrés. Par exemple, la première opération 101 et la deuxième opération 103 sur la période d'observation T peuvent 5 être réalisées avant que la première opération 101 et la deuxième opération 103 ne soient réalisées sur la période courante T'. Les résultats de la deuxième opération 103 menée sur une période d'observation T peuvent être mémorisés. Puis, dans un second temps, la première opération 101 et la deuxième opération 103 sont mises en oeuvre sur la période courante T'. 10 La période d'observation T peut être sélectionnée parmi plusieurs périodes d'observations en fonction des conditions de fonctionnement du réseau 1 au cours de la période courante T'. Par exemple, la période d'observation T est sélectionnée de sorte que les conditions d'utilisation du réseau 1 soient équivalentes à celles de la période courante T'. Si on souhaite surveiller le réseau 1 lorsque les consommations sont élevées et le réseau électrique fortement sollicité (par 15 exemple le soir en hiver), on peut sélectionner une période d'observation T correspondant à un soir d'hiver antérieur. Ainsi, l'influence des variations temporelles de la qualité des communications sur les résultats du procédé est réduite. Dans des modes de réalisation, par exemple pour une surveillance en temps réel, la période d'observation T peut être glissante. La période d'observation T peut être sélectionnée immédiatement antérieure à la période courante 20 T' de sorte les conditions d'utilisation du réseau restent similaires et comparables entre la période d'observation T et la période courante T'. Le procédé peut être mis en oeuvre de manière répétée, notamment en continu pour une surveillance en temps réel. Dans ces cas, les résultats obtenus sur une période courante T' 25 peuvent être mémorisés. Si la période courante T' ne comporte aucune anomalie d'alimentation électrique, elle peut devenir une période d'observation T pour une itération ultérieure du procédé. Pour faciliter la compréhension, les opérations sont présentées comme se déroulant les unes à la 30 suite des autres. En variante, les opérations de collecte 101 des retours des noeuds et d'attribution 103 des variables d'état B sont mises en oeuvre préalablement. Les données qui en résultent peuvent être stockées au moins temporairement. Puis, les opérations suivantes de traitement des données peuvent être mises en oeuvre dans un deuxième temps. 12 Dans une troisième opération 105, pour chaque noeud consommateur 5 et pour chacune des périodes T, T', il est déduit une distribution statistique P' P', des transitions d'états possibles (0K->OK, OK->K0, KO->OK, KO->K0) entre les instants t successifs de ladite période T, respectivement T'. 5 Au cours de la troisième opération 105, la succession des transitions entre deux variables d'état peut être modélisée sous la forme d'une chaîne de Markov à deux états OK, KO pour chaque noeud consommateur 5, en considérant que l'état à un instant est indépendant de l'état de l'instant précédent. Avec deux variables d'états OK et KO, il existe donc quatre transitions 10 possibles : (0K->OK), (0K->K0), (K0->OK), et (K0->K0). Par exemple, on considère la séquence suivante d'un noeud consommateur 5: OK -> KO -> OK -> OK-> KO -> KO -> OK -> KO -> KO -> KO -> KO 15 Dans cet exemple, le comptage des transitions (3,->By) possibles de variables d'état B, est le suivant : (Nb(OK OK) = 1 Nb(OK KO) = 3 Nb(K0 KO) = 4 Nb(K0 OK) = 2 La distribution statistique des transitions possibles peut aussi prendre la forme d'une matrice de 20 transition. Par exemple, les lignes correspondent à l'état à un instant t et les colonnes à l'état à l'instant suivant t+1. Les éléments de cette matrice représentent l'ensemble des probabilités P de transition d'un état vers un autre. Une telle matrice peut se calculer de la manière suivante : Nb(OK OK) P(OK OK) = Nb(OK OK)+ Nb(OK KO) Nb(OK KO) P(K0 KO) = Nb(K0 KO) + Nb(K0 OK) Nb(K0 OK) P(K0 OK) = Nb(K0 KO) + Nb(K0 OK) P(OK KO) = Nb(OK OK) + Nb(OK KO) Nb(K0 KO) 25 Avec l'exemple de séquence ci-avant, on obtient les probabilités de transition suivantes : 13 1 P(OK OK) = -4 3 P(OK KO) = -4 4 P(K0 KO) = -6 2 P(K0 OK)= -6 La matrice de transition prend donc la forme suivante : OK KO M = OK 0.25 0.75 KO 0.33 0.67 5 Dans une quatrième opération 107, pour chaque noeud consommateur 5, la vraisemblance L1, respectivement L,', des transitions de la période courante T' est calculé en fonction de la distribution statistique P, pendant la période d'observation T, respectivement en fonction de la distribution statistique P' pendant la période courante T'. 10 Le calcul de la vraisemblance L, par rapport à la distribution statistique P, pendant la période d'observation T permet de détecter un comportement inhabituel du noeud 5. Ainsi, en fonction du comportement habituel en communication du noeud 5, il est possible d'estimer la pertinence potentielle des défauts de réponse (KO) futurs pour déduire une anomalie d'alimentation. 15 Le calcul de la vraisemblance L', par rapport à la distribution statistique P', pendant la période courante T' permet de quantifier la variabilité des dernières transitions d'un noeud consommateur 5 au cours de la période courante T'. Le calcul de L,' permet de quantifier la prédictibilité, c'est-à-dire le caractère stochastique de la succession des dernières transitions présumée être un processus Markovien par comparaison avec le comportement du noeud 20 consommateur 5 durant cette même période courant T', et ce indépendamment du comportement des autres noeuds consommateurs 5. Ainsi, on peut distinguer, au sein de la période courante T', un comportement erratique (très bref donc probablement non représentatif) d'un comportement significatif sur la période courante T' (donc probablement représentatif). Cette quantification est indépendante du comportement des autres noeuds consommateurs 5. 25 14 Un exemple de calcul de la vraisemblance L,, L', de la distribution statistique 131, P', est décrit ci-après. Des calculs de vraisemblance différents peuvent être mis en oeuvre pour mesurer les mêmes aspects. 5 La vraisemblance L, par rapport à la période d'observation T est calculée selon l'équation suivante : Neval Li = br(Bx By) x log(Pi( Bx By)) x,y avec 10 i, l'indice identifiant le noeud consommateur 5, Neval le nombre total de transitions observées pendant une séquence d'évaluation de la période courante T', Nbr(Bx->By) le nombre de transitions de type Bx-*By observées pendant la période courante T', et 15 131(B, -> By) la probabilité d'observer la transition de type Bx-*By pendant la période d'observation T, l'ensemble des transitions de variables d'état possibles correspondant ici aux quatre transitions possibles : (0K->OK), (0K->K0), (K0->OK), et (K0->K0). 20 La vraisemblance L', par rapport à la période courante T' de la séquence d'évaluation est calculée selon l'équation suivante : L'1 = 1 NbT,(Bx By) x log(P' i(Bx By) Neval avec P',(B, -> By) la probabilité d'observer la transition de type Bx-*By pendant la période courante T'. 25 Dans cet exemple, le facteur 1/Nevai permet de moyenner la vraisemblance sur le nombre de transitions de la séquence d'évaluation de sorte que le résultat est rendu moins dépendant de la longueur de la séquence d'évaluation. 30 Dans une opération 109, pour un ensemble de plusieurs noeuds consommateurs 5, un indice S (ou score) d'écart à la normalité est calculé en fonction des vraisemblances L,, L' par rapport à 15 la période d'observation T et par rapport à la période courante T' de chaque noeud consommateur 5. L'indice S peut s'exprimer, par exemple, sous la forme d'une différence telle que dans l'équation suivante : S 5 Pour chaque noeud consommateur 5' - plus les transitions d'état de la séquence d'évaluation sont vraisemblables dans le contexte de la période courante, c'est-à-dire conformes au comportement récent du noeud consommateur 5' et 10 - moins les transitions d'état de la période d'observation T sont vraisemblables par rapport au comportement passé, c'est-à-dire peu conformes au comportement habituel du noeud consommateur 5' - plus l'indice S est élevé. 15 Ainsi, l'indice S est particulièrement élevé lorsqu'est observé une combinaison de : - un comportement non transitoire (ou stable) à l'échelle de la période d'évaluation (comparaison court terme), et - un comportement anormal par rapport à un comportement de référence observé lors d'une période comparable antérieure (comparaison à long terme). 20 Dans le cadre spécifique des communications des noeuds consommateurs 5, l'indice S peut subir une opération optionnelle de correction. Par exemple, l'indice S prend la forme d'une moyenne (ou somme) dont les éléments relatifs à des noeuds 5 particuliers sont ignorés, par exemple pour supprimer ou atténuer les effets dans la détection d'anomalies de noeuds présentant des 25 comportements aberrants ou dont les disfonctionnements sont connus et traités par ailleurs. Dans une opération 111, l'indice S ou l'indice corrigé S. obtenu est comparé à une valeur seuil S'f prédéfinie. 30 Dans une opération 113, un signal d'alerte est généré uniquement dans le cas où l'indice S (ou S.) obtenu dépasse la valeur seuil Sref- 16 La figure 4 représente les réponses (ou l'absence de réponse) au cours du temps de deux noeuds consommateurs 5' 5, disposés sur un réseau 1 commun Les deux noeuds consommateurs 5' 5, ont été sélectionnés de sorte qu'ils présentent un taux de réponse identique sur la période de temps observée. Autrement, dit la qualité de leur communication est, en valeur globale, 5 identique. Une unique anomalie d'alimentation électrique 7 a eu lieu pendant la période d'observation. Or, on constate sur la figure 4 que la succession d'absences de réponse KO (en noir) correspondant à l'anomalie 7 est plus aisément identifiable sur le chronogramme du noeud consommateur 5, que sur le chronogramme du noeud consommateur 5,. Autrement dit, bien que les deux noeuds consommateurs 5' 5, présentent des taux de réponse identiques, il est plus facile 10 de distinguer une anomalie d'alimentation à partir des erreurs de communication sur l'un des noeuds consommateurs que sur l'autre. Le noeud 5, de la figure 4 présente une répartition des erreurs de communication homogène dans le temps. L'alternance de réponses OK et d'absences de réponse KO est rapide. Ainsi, il semble 15 aisé de repérer une succession d'absences de réponse KO correspondant à une anomalie d'alimentation et présentant une longue durée parmi les brèves successions d'absences de réponse KO dues aux erreurs de communication. Par comparaison, le noeud 5, de la figure 4 présente une répartition des erreurs de communication hétérogène dans le temps, avec de longues successions d'absences de réponse KO. Ainsi, il semble difficile de distinguer une 20 succession d'absences de réponse KO correspondant à une anomalie d'alimentation d'une succession d'absences de réponse KO correspondant à des erreurs de communication. Les matrices de transition des noeuds 5, et 5 de la figure 4 sont reproduites ci-dessous : OK KO Mi = OK 0,45 0,55 KO 0,45 0,55 25 OK KO M = OK 0,9 0,1 KO 0,1 0,9 Les deux matrices amènent le constat suivant : les probabilités d'aller vers l'état KO (évènement se produisant en cas de coupure d'alimentation) diffèrent entre les deux noeuds consommateurs 5, et 5,. Pour le second noeud 5,, les probabilités d'atteindre l'un ou l'autre des 30 deux états est fortement dépendant de l'état précédent. Le procédé décrit ci-avant mesure cette dépendance au moyen de l'analyse statistique et l'utilise pour distinguer les noeuds en fonction 17 de leur pertinence. L'observation d'une transition est d'autant plus pertinente qu' elle est singulière (probabilité d'occurrence faible). L'exemple de la figure 4 illustre l'approche proposée par la demanderesse. 5 L'indice S tel que décrit ci-avant permet une analyse combinée du comportement statistique court terme et long terme de chaque noeud consommateur indépendamment les uns des autres. Le procédé permet de distinguer les comportements anormaux mais transitoires des comportements anormaux et durables à l'échelle de la période d'observation. Par construction 10 de l'indice S, les bruits statistiques sont éliminés de l'analyse, ce qui permet de limiter le taux de faux positifs (fausse alarme). En outre, l'analyse statistique est effectuée distinctement sur chaque noeud consommateur 5, ce qui permet de réduire les effets des différences intrinsèques des noeuds consommateurs 5 les uns 15 par rapport aux autres, notamment en terme de performance en communication. Le calcul final de l'indice S agrégé sur un ensemble de plusieurs noeuds consommateurs 5 permet de faire ressortir des anomalies faibles au niveau de chacun des noeuds consommateurs 5 mais simultanées. La remise à zéro en cas de réception de confirmation (OK) permet d'éliminer l'impact des séquences dépourvues de confirmations (KO) entrecoupées de confirmation (OK), 20 c'est-à-dire du bruit ne correspondant manifestement pas à une coupure de l'alimentation du noeud consommateur 5. La figure 5 montre un relevé de données dans des conditions réelles et en l'absence d'anomalie sur une pluralité de noeuds consommateurs 5, chaque noeud correspondant à une ligne, au cours 25 du temps représenté en abscisse. La figure 6 montre les mêmes données que pour la figure 5 à ceci près qu'une coupure d'alimentation est simulée sur certains noeuds consommateurs 5. La figure 7 représente l'évolution de l'indice mesuré selon le procédé décrit ci-avant au cours du temps à partir des données représentées en figures 6. La courbe 300 correspond à la valeur de l'indice S sans correction tandis que la courbe 400 correspond à la valeur de l'indice S. après 30 correction. On constate que chacun des deux procédés permet d'identifier sans ambiguïté la période concernée par une anomalie alors même qu'une faible part des noeuds 5 est concernée par la coupure d'alimentation. 35 18 En outre, l'application de la correction permet de faire disparaître certaines variations 301 non pertinentes et visibles sur la courbe 300. La figure 8 représente le taux de détection de coupures (en ordonnée et en pourcentage) en 5 fonction du nombre de faux positifs (en abscisse et en nombre par semaine d'un ensemble de noeuds test). Quatre méthodes sont comparées. Une première méthode, correspondant à la courbe référencée 4000, consiste à dénombrer à chaque instant t le nombre de noeuds consommateurs 5 pour lesquels une absence de 10 confirmation (KO) est observée et de déclencher une alerte au-dessus d'une valeur seuil. Cela correspond à la méthode présentée ci-avant en référence à la figure 2. Une deuxième méthode, correspondant à la courbe référencée 3000, consiste à dénombrer à chaque instant t le nombre d'absences de confirmation (KO) et de déclencher une alerte au-15 dessus d'une valeur seuil en excluant du compte les noeuds présentant un taux général de réponse trop faible (technique de « liste noire »). La courbe référencée 2000 correspond à un mode de réalisation du procédé, en l'absence de correction de l'indice S. 20 La courbe référencée 1000 correspond à un mode de réalisation du procédé, avec application de la correction de l'indice S.. Les essais comparatifs de la figure 8 ont été menés dans les conditions suivantes : 25 - panel de 103 noeuds consommateurs communiquant par CPL (sous protocole « G1 ») sélectionnés parmi des noeuds fonctionnant en conditions réelles dans le contexte du déploiement des compteurs « Linky » en France ; - les procédés correspondant aux courbes 1000 et 2000 ont été mis en oeuvre par ordinateur, ici par implémentation dans un logiciel d'un concentrateur ; 30 - les requêtes sont espacées d'un intervalle constant et supérieur à une seconde ; - les calculs sont effectués de manière itérative en temps réel ; - la période d'observation T correspond aux 15 jours précédents la période courante T' et aux données enregistrées distinctement pour chaque noeud. La période d'observation T est donc, ici, glissante dans le temps ; 19 - la séquence d'évaluation est définie comme la suite des 15 dernières transitions (Neva' = 15) pour chaque noeud ; - les coupures d'alimentation sont simulées par l'attribution artificielle de la valeur « KO » à la variable d'état l31,' d'un sous ensemble de noeuds pendant une période commune (passage des 5 données brutes réelles de la figure 5 à celles simulées de la figure 6). Les courbes de la figures 8 sont obtenues en faisant varier le seuil de déclenchement pour les deux premières méthodes, respectivement la valeur de référence S'f pour les deux procédés de traitement statistique. Les marges d'erreur du redressement sont calculées par inférence 10 statistique (« bootstrap ; ré-échantillonnage des ensembles de noeuds avec 1000 itérations). Un exemple de résultat comparé est donné dans le tableau suivant en fixant le nombre de fausses alarmes par semaine à dix : Méthode Taux de détection Comptage KO (courbe 4000) 4,60 % Comptage KO avec liste noire (courbe 3000) 63,8 % Indice statistique sans correction (courbe 2000) 77,1 % Indice statistique avec correction (courbe 1000) 88,1 % 15 Réciproquement, le tableau suivant donne un nombre de fausses alarmes (par grappe et par semaine) pour un taux de détection des anomalies fixé à 50% : Méthode Nombre de fausses alarmes Comptage KO 110 Comptage KO avec liste noire 3,76 Score statistique sans correction 0,625 Score statistique avec correction 0,00305 20 Le procédé du score statistique avec correction comparé à la méthode du simple comptage des absences de réponses (KO) permet donc de diminuer le nombre de fausses alarmes par un facteur 36 000 tout en conservant le même taux de détection. Quel que soit le compromis choisi entre taux de détection et fausse alarme par réglage du seuil 25 de déclenchement, les procédés proposés ici permettent d'améliorer significativement la 20 situation par rapport à une surveillance par comptage, avec ou sans application d'une liste noire. Les techniques de liste noire, bien que permettant de réduire le nombre de fausse alarme, conduisent souvent à ignorer une partie du réseau, ce qui n'est pas souhaitable. 5 Dans des variantes, le procédé peut être appliqué en parallèle sur plusieurs groupes de noeuds 5 d'un même réseau 1. Les groupes peuvent être complètement distincts ou bien comprendre des noeuds 5 communs. La composition de chaque ensemble peut être sélectionnée en fonction des positions des noeuds consommateurs 5 sur le réseau 1 de sorte qu'un signal d'alerte est associé à une zone prédéfinie du réseau 1. Ainsi, lorsqu'une alerte est générée pour un groupe, le signal 10 peut être associé à une zone du réseau 1. L'invention ne se limite pas aux exemples de procédé et de système décrits ci-avant, seulement à titre d'exemple, mais elle englobe toutes les variantes que pourra envisager l'homme de l'art REVENDICATIONS1. Procédé mis en oeuvre par ordinateur de détection d'anomalie sur un réseau (1) de distribution électrique incluant une pluralité de noeuds consommateurs (5), chacun des noeuds 5 consommateurs (5) étant apte à transmettre une confirmation via le réseau (1) de distribution électrique en réponse à une requête reçue par l'intermédiaire du réseau (1) de distribution électrique, le procédé comprenant : a) collecter (101) les retours de requêtes transmises à chacun des noeuds consommateurs (5) via le réseau (1) de distribution électrique à une pluralité d'instants (t) successifs d'une période 10 d'observation (T) et à une pluralité d'instants (t') successifs d'une période courante (T'), b) pour chaque requête, établir (103) une variable d'état (B') de type booléen dont la valeur correspond respectivement à la réception d'une confirmation (0K) ou à l'absence de réception d'une confirmation (KO), c) pour chaque noeud consommateur (5) et pour chacune des périodes (T, T'), déduire (105) une 15 distribution statistique (P,, P',) des transitions d'états possibles (0K->OK, OK->K0, KO->OK, KO->K0) entre deux instants successifs, d) pour chaque noeud consommateur (5) et pour chacune des périodes (T, T'), calculer (107) une vraisemblance (L,, L',) de la distribution statistique (P,, P',), e) pour un ensemble de plusieurs noeuds consommateurs (5), calculer (109) un indice (S; Scor) 20 d'écart à la normalité en fonction des vraisemblances (L' L',) de la période d'observation (T) et de la période courante (T') des noeuds consommateurs (5), f) comparer (111) l'indice (S ; Scor) obtenu à une valeur seuil (Sreô prédéfinie, g) uniquement dans le cas où l'indice (S; Scor) obtenu dépasse la valeur seuil (Sref), générer (113) un signal d'alerte. 25 2. Procédé selon la revendication 1, dans lequel l'indice (S; Scor) d'écart à la normalité est équivalent à une somme ou une moyenne des différences de vraisemblances (LI, L',) entre la période d'observation (T) et la période courante (T') de chaque noeud consommateur (5). 30 3. Procédé selon l'une des revendications précédentes, dans lequel les requêtes et les réceptions de confirmations (0K) sont effectuées par Courtant Porteur en Ligne sur le réseau (1) de distribution électrique. 22 4. Procédé selon l'une des revendications précédentes, dans lequel la période d'observation (T) est sélectionnée parmi plusieurs en fonction des conditions de fonctionnement du réseau (1) de distribution électrique au cours de la période courante (T'). 5. Procédé selon l'une des revendications précédentes, dans lequel : - la période d'observation (T) s'étend sur une durée comprise entre 1 seconde et 100 jours, - la période courante (T') s'étend sur une durée comprise entre 10 millisecondes et 7 jours, - les instants successifs (t, t+1) de l'une au moins de la période d'observation (T) et de la période courante (T') sont mutuellement espacés d'une durée comprise entre 1 milliseconde et 10 24 heures, et/ou - l'une au moins de la période d'observation (T) et de la période courante (T') comprend entre 2 et 500 instants successifs. 6. Procédé selon l'une des revendications précédentes, dans lequel les opérations sont mises 15 en oeuvre de manière répétée de sorte que les variable d'état (Bu) relatives à la période courante (T') sont mises à jour sensiblement en temps réel. 7. Procédé selon l'une des revendications précédentes, dans lequel les opérations sont mises en oeuvre distinctement pour plusieurs ensembles de noeuds consommateurs (5) d'un même 20 réseau (1). 8. Procédé selon la revendication 7, dans lequel la composition de chaque ensemble est sélectionnée en fonction des positions des noeuds consommateurs (5) sur le réseau (1) de sorte qu'un signal d'alerte est associé à une zone prédéfinie du réseau (1). 25 9. Module de détection d'anomalie sur un réseau (1) de distribution électrique comprenant une unité de traitement, le réseau (1) incluant une pluralité de noeuds consommateurs (5), l'unité de traitement (3) étant reliée au réseau (1) et chacun des noeuds consommateurs (5) étant apte à transmettre une confirmation (13) via le réseau (1) de distribution électrique en réponse à une 30 requête (11) reçue par l'intermédiaire du réseau (1) de distribution électrique, l'unité de traitement étant agencée de manière à : a) interroger (101) chacun des noeuds consommateurs (5) via le réseau (1) de distribution électrique à une pluralité d'instants (t) successifs d'une période d'observation (T) et à une pluralité d'instants (t') successifs d'une période courante (T'), 23 b) pour chaque interrogation, établir (103) une variable d'état (Bu) de type booléen dont la valeur correspond respectivement à la réception d'une confirmation (0K) ou à l'absence de réception d'une confirmation (KO), c) pour chaque noeud consommateur (5) et pour chacune des périodes (T, T'), déduire (105) une 5 distribution statistique (Pi, P',) des transitions d'états possibles (0K->OK, OK->K0, KO->OK, KO->K0) entre deux instants successifs, d) pour chaque noeud consommateur (5) et pour chacune des périodes (T, T'), calculer (107) une vraisemblance (Li, L'i) de la distribution statistique (Pi, P',), . e) pour un ensemble de plusieurs noeuds consommateurs (5), calculer (109) un indice (S ; SCOL 1 10 d'écart à la normalité en fonction des vraisemblances (Li, L'i) de la période d'observation (T) et de la période courante (T') des noeuds consommateurs (5), f) comparer (111) l'indice (S ; Scor) obtenu à une valeur seuil (Sref) prédéfinie, g) uniquement dans le cas où l'indice (S; Scor) obtenu dépasse la valeur seuil (Sref), générer (113) un signal d'alerte. 15 10. Programme informatique comportant des instructions pour la mise en oeuvre du procédé selon l'une des revendications 1 à 8, lorsque ce programme est exécuté par un processeur. 11. Support d'enregistrement non transitoire lisible par un ordinateur sur lequel est enregistré 20 un programme comprenant des instructions pour la mise en oeuvre du procédé selon l'une des revendications 1 à 8.