FR 3087055 A1 20200410 FR 1859197 A 20181004 Estimation d'une grandeur physique par un système de mesure distribué La présente invention concerne de manière générale l'estimation d'une 5 grandeur physique variable dans le temps par un système de mesure distribué comportant une plateforme centralisée et un ensemble composé d'une pluralité d'équipements de mesure distribués, aptes chacun à mesurer une donnée locale représentative de ladite grandeur physique et à communiquer à distance avec ladite plateforme centralisée par l'intermédiaire d'un réseau de 10 communication. Un domaine d'application visé tout particulièrement, bien que non exclusivement, par l'invention est celui de la gestion et/ou de la surveillance de la consommation en électricité d'appareils électriques, par exemple des appareils électriques à usage domestique, situés chez un nombre 15 particulièrement élevé de consommateurs et reliés à un réseau de distribution d'électricité. Dans de tels systèmes, il est classique d'installer au moins un équipement de mesure sur une pluralité de sites (domiciles, bâtiments d'entreprises, collectivités...), relié localement à au moins un appareil 20 électrique sous surveillance de façon à pouvoir effectuer, de préférence en continu, des mesures de consommation électrique de cet appareil électrique. Une plateforme centralisée dédiée, sous forme d'un ou de plusieurs serveurs, obtient alors, de préférence en temps réel, des échantillons de mesures de consommation électrique communiqués par chaque équipement de mesure 25 local du système. La plateforme centralisée peut ainsi déterminer les courbes de charge élémentaires, c'est-à-dire la variation temporelle de la consommation relative à chacune des charges ou appareils électriques sous surveillance, et en déduire, par agrégation de ces courbes de charge élémentaires, une courbe de charge globale. 30 Pour un système dédié à la surveillance seule de la consommation en électricité, les équipements de mesure distribués sur les différents sites peuvent être typiquement des compteurs électriques dits intelligents, aptes à mesurer la puissance électrique consommée par un ou plusieurs appareils 2 électriques de manière à délivrer périodiquement des échantillons de mesure de la consommation électrique. La Demanderesse propose également un système distribué de modulation de la consommation d'un ensemble composé d'un nombre très élevé d'appareils électriques (tels que des chauffe-eau, des radiateurs électriques, des climatiseurs...) situés chez un grand nombre d'utilisateurs, dans lequel les différentes courbes de charge obtenues par des équipements de mesure distribués permettent notamment de proposer un service dit d'effacement diffus, par lequel il est possible de sélectionner en temps réel, dans l'ensemble des équipements de mesure reliés localement aux appareils électriques, un sous-ensemble d'équipements de mesure auxquels le système, en particulier la plateforme centralisée, va envoyer des ordres de coupure momentanée de fonctionnement, de manière à réduire la consommation d'énergie de l'ensemble d'une valeur de consigne donnée. Un système de ce type est par exemple dans la demande W02008/017754, ou dans la demande W02012/172242. Ce système, en étant également apte à moduler, via des ordres de coupure et de redémarrage, sélectivement l'alimentation de certains appareils électriques par analyse de la consommation d'un grand nombre d'utilisateurs, permet d'adapter la consommation électrique à la production d'électricité disponible à un instant donné. Lorsque que les systèmes de mesure sont massivement distribués, avec un ensemble comportant un nombre très important d'équipements de mesure distribués, la connaissance de la grandeur physique agrégée à partir des échantillons de mesure prélevés par les différents équipements de mesure requiert la transmission d'un nombre très élevé, typiquement de centaines de milliers (voire plus), de mesures par les équipements de mesure. Quel que soit le réseau de communication utilisé pour la transmission des données, cela génère une occupation importante de la bande passante du réseau. En outre, pour certains types de réseau de communication tels que des réseaux cellulaires, les volumes de données nécessaires à cette transmission peuvent impliquer des coûts significatifs. Ce problème constitue en conséquence une limitation à l'augmentation du périmètre de ces systèmes, c'est-à-dire du nombre d'équipements de mesure distribués. 3 Le problème précédent est accru dans les systèmes de mesure de consommation électrique massivement distribués, tels que celui décrit dans les demandes W02008/017754 et W02012/172242 précitées, pour lesquels il est important de connaître en temps réel la courbe de charge agrégée exacte à 5 partir de la transmission en temps réel d'un grand nombre de mesures de consommation. Dans de tels systèmes, la plateforme centralisée obtient des échantillons de mesure de consommation électrique à une période d'échantillonnage généralement de l'ordre de la dizaine de minutes. Une 10 courbe de charge agrégée plus précise pourrait être obtenue en réduisant la période d'échantillonnage, par exemple en utilisant une période d'échantillonnage fortement réduite, typiquement de l'ordre de la dizaine de secondes, voire de la seconde. Le nombre des mesures qui doivent être non seulement collectées mais également agrégées sur une période 15 d'échantillonnage par la plateforme centralisée devient dans ce cas trop important. Pour résoudre le problème lié au nombre important de mesures collectées, directement dépendant du nombre d'équipements de mesure compris dans l'ensemble et/ou de la fréquence d'échantillonnage des mesures, 20 il est connu de pratiquer à un nombre limité de mesures collectées au sein d'un sous-ensemble, dit de sondage, puis d'estimer la grandeur physique, par exemple la courbe de charge globale, par extrapolation des mesures collectées à partir de ce sous-ensemble. La précision de l'estimation dépend très fortement de la façon dont le sous-ensemble de sondage a été constitué. 25 Or, à ce jour, les solutions connues reposent sur des techniques de sondage conventionnelles pour lesquelles la constitution d'un échantillon représentatif, c'est-à-dire ici, la sélection du sous-ensemble d'équipements de mesure sondés, se fait classiquement, une fois pour toutes et à l'initiative du sondeur, par pur aléa, par stratification, par grappe ou par la méthode des 30 quotas. Ces types de sondage génèrent donc la pré-sélection d'une partie seulement des instruments de mesure constituant l'ensemble dans le système, dont l'appartenance à l'échantillon considéré comme représentatif est due à 4 l'observation d'un comportement passé (consommation lors de la dernière heure, de la dernière journée, etc.) ou d'une configuration déterminée (équipement mesuré, tarification de l'abonné, localisation, etc.). Cet échantillon ou sous-ensemble représentatif, quelle que soit la 5 méthode de constitution employée, est statique, ou quasi-statique, c'est-à-dire que sa composition varie faiblement par rapport à la fréquence de remontée des mesures. Or, au-delà de la difficulté notoire à sélectionner un sous-ensemble représentatif sur une période d'observation, il n'est pas garanti que cette 10 représentativité observée dans le passé prédise convenablement des mesures à venir. L'expérience montre que, pour qu'un sous-ensemble soit représentatif sur des observations à venir, il faut que le sous-ensemble sondé soit d'une très grande taille, ce qui limite considérablement les avantages du sondage. 15 La présente invention a pour but de pallier les inconvénients des systèmes utilisant les techniques de sondage conventionnelles précitées. Pour ce faire, l'invention a pour objet un procédé d'estimation d'une grandeur physique variable dans le temps par un système de mesure distribué comportant une plateforme centralisée et un ensemble composé d'une 20 pluralité d'équipements de mesure distribués aptes chacun à communiquer à distance avec ladite plateforme centralisée par l'intermédiaire d'un réseau de communication, le procédé comportant une étape de mesure, par chaque équipement de mesure dudit ensemble, d'une donnée locale représentative de ladite grandeur physique, une étape de collecte, sur une période 25 d'échantillonnage prédéterminée, d'échantillons de mesure provenant d'un sous-ensemble dudit ensemble dans laquelle chaque équipement de mesure appartenant audit sous-ensemble transmet à ladite plateforme centralisée un échantillon de mesure de la donnée locale via ledit réseau de communication, et une étape d'estimation de ladite grandeur physique sur la période 30 d'échantillonnage à partir d'une agrégation des échantillons de mesure collectés, le procédé étant caractérisé en ce qu'il comporte, à chaque période d'échantillonnage, une phase d'exécution, au niveau de chaque équipement de 5 mesure dudit ensemble, d'un algorithme local de décision apte à décider aléatoirement, sur la base d'une probabilité d'envoi prédéfinie pour chaque équipement de mesure, si l'équipement de mesure correspondant fait partie dudit sous-ensemble, et en ce que ladite grandeur physique est estimée à 5 partir des échantillons de mesure collectés et de la probabilité d'envoi prédéfinie pour chaque équipement de mesure appartenant audit sousensem ble. Ainsi, la présente invention permet de créer dynamiquement et aléatoirement, à chaque période d'échantillonnage, le sous-ensemble 10 d'équipements de mesure qui va participer au sondage et permettre par agrégation une estimation précise et rapide de la grandeur physique. En outre, la création du sous-ensemble s'effectue par le biais seul des équipements de mesure qui décident ou non, de façon aléatoire, de participer au sondage à chaque période d'échantillonnage. 15 Selon un mode d'implémentation possible du procédé selon l'invention, la probabilité d'envoi prédéfinie pour au moins un équipement de mesure dudit ensemble a une valeur variable en fonction de la valeur de l'échantillon de mesure effectuée et/ou de conditions de réalisation de la mesure et/ou de l'instant de la mesure; dans ce cas, si ledit au moins un équipement de 20 mesure fait partie dudit sous-ensemble, ledit au moins un équipement de mesure transmet la probabilité d'envoi à la plateforme centralisée via ledit réseau de communication pendant l'étape de collecte, en même temps que l'échantillon de mesure de la donnée locale. Selon un autre mode d'implémentation possible du procédé selon 25 l'invention, la probabilité d'envoi prédéfinie est transmise à la plateforme centralisée par chaque équipement de mesure appartenant audit sous-ensemble via ledit réseau de communication pendant l'étape de collecte, en même temps que l'échantillon de mesure de la donnée locale. Dans une implémentation possible, le procédé peut comporter en outre 30 une étape de mise à jour de l'algorithme local de décision et/ou de la probabilité d'envoi prédéfinie pour au moins un équipement de mesure dudit ensemble. 6 La grandeur physique globale peut être estimée par calcul selon la relation suivante dans laquelle N' est le nombre d'équipements de mesure du sous-ensemble, h est la valeur de la donnée locale mesurée par un équipement i de mesure, 5 Ki est la probabilité d'envoi prédéfinie pour un équipement i de mesure et T est la période d'échantillonnage. Dans certains cas, la période d'échantillonnage peut être inférieure ou égale à 10 minutes. La période d'échantillonnage peut notamment être avantageusement 10 sensiblement de l'ordre de la seconde ou de la dizaine de secondes, ce qui est particulièrement utile pour un système devant estimer en temps réel la grandeur physique sur une période d'échantillonnage. L'invention a également pour objet un système de mesure distribué pour estimer une grandeur physique variable dans le temps, le système 15 comportant une plateforme centralisée et un ensemble composé d'une pluralité d'équipements de mesure distribués aptes chacun à communiquer à distance avec ladite plateforme centralisée par l'intermédiaire d'un réseau de communication, le système étant configuré de sorte que, lors d'une étape de mesure, chaque équipement de mesure dudit ensemble mesure une donnée 20 locale représentative de ladite grandeur physique, lors d'une étape de collecte, chaque équipement de mesure appartenant à un sous-ensemble dudit ensemble transmet à ladite plateforme centralisée un échantillon de mesure de la donnée locale via ledit réseau de communication sur une période d'échantillonnage prédéterminée, et lors d'une étape d'estimation, la 25 plateforme centralisée estime ladite grandeur physique sur la période d'échantillonnage à partir d'une agrégation des échantillons de mesure collectés, caractérisé en ce que chaque équipement de mesure dudit ensemble comporte un module logiciel configuré pour exécuter, à chaque période d'échantillonnage, un algorithme local de décision apte à décider 30 aléatoirement, sur la base d'une probabilité d'envoi prédéfinie pour chaque 7 équipement de mesure, si l'équipement de mesure correspondant fait partie dudit sous-ensemble, et en ce que la plateforme centralisée est configurée pour estimer la grandeur physique à partir des échantillons de mesure collectés et de la probabilité d'envoi prédéfinie pour chaque équipement de 5 mesure appartenant audit sous-ensemble. La donnée locale est par exemple une puissance électrique soutirée par une charge locale comprenant un ou plusieurs appareils électriques domestiques raccordés à un réseau de distribution d'électricité, et la grandeur physique est alors la puissance électrique totale consommée. 10 L'invention sera mieux comprise à la lecture de la description qui suit, faite en référence aux figures annexées, dans lesquelles : - la figure 1 donne un exemple d'architecture d'un système de mesure distribué apte à implémenter l'invention; - la figure 2 illustre des étapes susceptibles d'être mises en oeuvre 15 dans un procédé d'estimation d'une grandeur physique conforme à l'invention, par exemple par le système de la figure 1. Dans la suite de l'exposé, l'invention sera décrite dans le cadre de la détermination d'une courbe de charge globale comprenant une agrégation de courbes de charge élémentaires correspondant à l'évolution temporelle de la 20 consommation électrique d'une pluralité d'équipements électriques. La grandeur physique à estimer est donc ici la puissance électrique totale consommée par une pluralité d'appareils électriques. Comme indiqué ci-avant, le principe de l'invention peut être néanmoins appliqué à d'autres grandeurs physiques, par exemple à d'autres types d'énergies consommées. 25 La figure 1 illustre l'architecture d'un système 1 de mesure distribué permettant l'estimation de la puissance électrique consommée par des appareils électriques (non représentés), par exemple à usage domestique, connectés à un réseau de distribution en électricité (non représenté) et situés sur un nombre N de sites. Le système comporte une plateforme centralisée 1, 30 hébergée par un ou plusieurs serveurs Internet, et un ensemble composé ici de N équipements de mesure 21,...2;,...2N distribués sur les N sites. Dans la 8 suite de l'exposé, l'indice i est utilisé pour toute référence à un équipement de mesure 2, situé sur le site i, i étant un entier variant de 1 à au plus N. Chaque équipement de mesure 2, est apte à communiquer à distance avec la plateforme centralisée 1 par l'intermédiaire d'un réseau 3 de 5 communication. A cet effet, chaque équipement de mesure 2, comporte par exemple un modem 20, de communication sans fil permettant d'établir avec un modem 10 de communication correspondant au niveau de la plateforme centralisée 1 une connexion de type téléphonie par paquet, tel que le GPRS, la 3G ou la 4G. Alternativement, la connexion peut s'effectuer via une liaison de 10 type Ethernet. Chaque équipement de mesure 2, est en outre apte à mesurer la puissance électrique consommée localement par les appareils électriques auxquels il est relié, par l'intermédiaire d'un module 22, de mesure des tensions et des courants. Les équipements de mesure 2, sont par exemple des 15 compteurs électriques intelligents installés sur chaque site i. Dans une autre architecture possible ou le système peut effectuer une modulation de la puissance électrique en proposant un service d'effacement diffus, chaque équipement de mesure 2, est composé d'un ou plusieurs boîtiers tels que des boîtiers pilote et boîtiers modulateurs décrits dans les demandes 20 W02008/017754 et W02012/172242 précitées, incorporées par référence. A la différence d'un compteur, ces équipements sont en outre aptes à exécuter des ordres de coupure/redémarrage de fonctionnement d'un ou plusieurs appareils électriques auxquels ils sont reliés, ordres qu'ils reçoivent de la plateforme centralisée pour la gestion d'effacements ponctuels. 25 Dans tous les cas, afin d'éviter une transmission importante de données, le système est configuré pour estimer la puissance totale consommée à partir d'échantillons de mesure effectués par un sous-ensemble d'équipements de mesure, ou sous-ensemble de sondage considéré comme représentatif. Plus précisément, le système est configuré pour implémenter les 30 étapes suivantes : 9 - une étape de mesure, lors de laquelle chaque équipement de mesure 2; de l'ensemble de N équipements mesure la puissance électrique Pi consommée localement ; - une étape de collecte, lors de laquelle seul un nombre N' inférieur à N 5 d'équipements de mesure 2. correspondant au sous-ensemble de sondage, transmet à la plateforme centralisée 1 un échantillon de mesure de la puissance électrique P; consommée localement via le réseau 3 de communication sur une période d'échantillonnage T prédéterminée ; et - une étape d'estimation, lors de laquelle la plateforme centralisée 1 10 estime la puissance consommée totale sur la période d'échantillonnage T à partir d'une agrégation des échantillons de mesure collectés. A la différence des systèmes utilisant des techniques conventionnelles de sondage avec un sous-ensemble représentatif fixe, le système selon l'invention prévoit de changer, à chaque période d'échantillonnage, le 15 périmètre du sous-ensemble de sondage, c'est-à-dire le choix des équipements de mesures qui vont effectivement transmettre leur mesure lors de l'étape de collecte. Pour ce faire, conformément à l'invention, chacun des N équipements de mesure 2; exécute, à chaque période d'échantillonnage, un algorithme local de 20 décision mémorisé dans un module logiciel 21; et apte à décider aléatoirement, sur la base d'une probabilité d'envoi mi prédéfinie pour chaque équipement 2; de mesure, si l'équipement de mesure correspondant va envoyer ou non son échantillon de mesure P; à la plateforme centralisée 1. En d'autres termes, chaque équipement de mesure 2; est apte à 25 décider, de façon aléatoire, s'il va ou non faire partie du sous-ensemble de sondage transmettant les échantillons de mesure sur la période d'échantillonnage T. Dans un mode de réalisation possible, tous les équipements de mesure 2; comportent la même probabilité d'envoi rri prédéfinie et le même algorithme 30 local de décision. L'algorithme local de décision peut suivre une loi de Bernoulli et implémenter par exemple les étapes suivantes : 10 - générer un nombre aléatoire, de valeur comprise par exemple entre 0 et 100 ; - si le nombre aléatoire généré est inférieur à la valeur affectée au pourcentage de la probabilité d'envoi rri, par exemple inférieure à 10 si la 5 probabilité d'envoi Tri est égale à 10%, l'équipement de mesure 2; transmet l'échantillon de mesure Pi à la plateforme centralisée 1 (en d'autres termes l'équipement de mesure fait alors partie du sous-ensemble de sondage). - sinon, l'équipement de mesure 2. ne transmet pas l'échantillon de mesure 9 à la plateforme centralisée 1. 10 Dans un autre mode de réalisation possible, la probabilité d'envoi n-i peut avoir une valeur fixe prédéfinie, mais différente d'un équipement de mesure à l'autre. Dans un autre mode de réalisation possible, la probabilité d'envoi n-i pour un ou plusieurs équipements 2; de mesure dudit ensemble (voire pour la 15 totalité des équipements de mesure) a une valeur variable en fonction de la valeur de l'échantillon 9 de mesure effectuée et/ou de conditions de réalisation de la mesure et/ou de l'instant de la mesure. Ainsi, à titre d'exemple non limitatif, on peut prévoir que la probabilité d'envoi Ili pour un équipement 2; de mesure soit égale à 0% si la valeur de l'échantillon 9 de 20 mesure effectuée sur la période d'échantillonnage considérée est nulle. En d'autres termes, cet équipement 2; de mesure n'envoie pas les mesures nulles, et ne fera pas partie du sous-ensemble de sondage. Si l'équipement 2; de mesure fait partie d'un système de modulation proposant le service d'effacement diffus, on peut également prévoir, à titre d'exemple, que la 25 probabilité d'envoi 7CÉ pour cet équipement 2; de mesure soit égale à 100% si la valeur de l'échantillon 9 de mesure effectuée sur la période d'échantillonnage considérée n'est pas nulle alors qu'un ordre de coupure des appareils électriques reliés à cet équipement devrait être en cours d'exécution. En d'autres termes, un équipement 2; de mesure peut choisir d'envoyer 30 systématiquement sa mesure dès lors qu'il juge d'une situation anormale. La plateforme centralisée doit, comme il sera apparent dans la suite, connaître la probabilité d'envoi zi correspondant à tout échantillon de mesure 11 qu'elle collecte. Aussi, on prévoit que les équipements 2; de mesure appartenant au sous-ensemble de sondage et dont la probabilité d'envoi varie, quelle qu'en soit la raison, transmettent leur probabilité d'envoi Tri à la plateforme centralisée 1 via ledit réseau 3 de communication pendant l'étape 5 120 de collecte, en même temps que leur échantillon de mesure P. Dans un autre mode de réalisation possible, on peut prévoir que la probabilité d'envoi Tri prédéfinie soit transmise systématiquement, qu'elle soit variable ou non, à la plateforme centralisée 1 par chaque équipement 2; de mesure appartenant audit sous-ensemble via ledit réseau 3 de communication 10 pendant l'étape 120 de collecte, en même temps que l'échantillon de mesure Il est possible de prévoir qu'au moins un ou plusieurs, voire la totalité des équipements de mesure de l'ensemble puisse réaliser des mises à jour de l'algorithme local de décision et/ou de la probabilité d'envoi Tri prédéfinie. Les 15 mises à jour peuvent être transmises par la plateforme centralisée 1. En variante, les mises à jour sont réalisées sur la base d'informations calculées ou mesurées en local par les équipements de mesure. Une fois que la plateforme centralisée 1 a collecté les différents échantillons de mesure IP transmis par les N' équipements de mesure du sous- 20 ensemble de sondage, et, le cas échéant, les probabilités d'envoi Tri correspondantes, elle va pouvoir traiter ces échantillons au niveau d'un module de traitement 11, et notamment estimer la puissance totale P(T) consommée sur la période d'échantillonnage T à partir des échantillons de mesure P; collectés et de la probabilité d'envoi Mi prédéfinie pour chaque 25 équipement de mesure appartenant audit sous-ensemble de sondage. On montre notamment qu'une bonne estimation de la puissance totale P(T) consommée sur une période d'échantillonnage T peut être obtenue par la relation suivante : pcn =v Pi(T) z1(T) i=1 12 laquelle ne prend en compte que les N' équipements de mesure du sous-ensemble de sondage. La relation précédente part bien entendu de l'hypothèse que les mesures transmises par les N' équipements sont toutes effectivement reçues. 5 Dans le cas contraire, la relation s'applique au nombre de mesures réellement collectées. L'expérience montre que dans le cas d'une valeur à estimer sur la durée (suite temporelle), la méthode de génération dynamique d'un sous-ensemble de sondage à chaque période d'échantillonnage T est bien plus précise que les 10 méthodes dites statiques ou quasi statiques. Cette précision peut s'expliquer par le fait que le biais de sélection s'apparente à des fluctuations dues à la constitution dynamique du sous-ensemble de sondage, et que la moyenne des erreurs s'approche alors de zéro. La figure 2 résume des étapes susceptibles d'être mises en oeuvre dans 15 un procédé d'estimation d'une grandeur physique conforme à l'invention, par exemple par le système de la figure 1 comportant N équipements de mesure 2; distribués, et pour lequel la grandeur physique est la puissance électrique totale consommée sur une période d'échantillonnage T : Le procédé d'estimation comporte une première étape 100 de mesure 20 d'une donnée locale représentative de la grandeur physique à estimer par chacun des N équipements 2; de mesure constituant l'ensemble de départ. Lors d'une étape 110, chacun des N équipements 2; de mesure de l'ensemble exécute son algorithme local de décision, afin de décider aléatoirement s'il va transmettre ou non un échantillon de mesure P; sur la 25 période d'échantillonnage T considérée. Cette étape 110 permet de déterminer dynamiquement le sous-ensemble de sondage constitué de N' équipements 2; de mesure qui ont chacun décidé, à l'issue de l'exécution de leur algorithme local de décision, de transmettre leur échantillon de mesure P; selon leur probabilité d'envoi ir1. Dans le cas où l'équipement de mesure considéré 30 possède, comme mentionné précédemment, plusieurs valeurs possibles de probabilité d'envoi 1Ti, cette étape comporte en outre la sélection (non 13 représentée) de la valeur de probabilité d'envoi qui sera effectivement associée à la mesure. Une fois le sous-ensemble de sondage déterminé pour la période d'échantillonnage T considérée, le procédé se poursuit par une étape 120 de 5 collecte, par la plateforme centralisée 1, des échantillons de mesure Pi transmis par les seuls équipements 2i de mesure du sous-ensemble de sondage. Comme on l'a vu précédemment, la probabilité d'envoi peut ou doit dans certains cas être également transmise avec l'échantillon de mesure Pi correspondant. 10 La dernière étape 130 comprend l'estimation de la grandeur physique P(T) sur la période d'échantillonnage T à partir des au plus N' échantillons P1 de mesure collectés à l'étape 120 par la plateforme centralisée 1 et des probabilités d'envoi Tri associées, de préférence par application de la relation P(T); P(T) = (T) 1=1 Les différentes étapes sont réitérées à chaque période 15 d'échantillonnage, de façon à obtenir à chaque fois un sous-ensemble de sondage différent. Cela permet notamment de faire participer l'intégralité des équipements de mesure et de conserver la représentativité sur la durée. La période d'échantillonnage T peut être de l'ordre de 10 minutes. 20 Néanmoins, l'utilisation d'un sous-ensemble de sondage créé dynamiquement permet de réduire considérablement cette période d'échantillonnage. On a pu montrer notamment qu'il était possible d'estimer précisément une courbe de charge agrégée pour un système massivement distribué comportant une centaine de milliers de sites avec une période d'échantillonnage T de l'ordre de 25 10 secondes, voire de la seconde, sans engendrer de coûts conséquents en termes de communication des données. La méthode d'estimation selon l'invention permet de réduire considérablement le nombre de sites devant envoyer des mesures (on passe de 80% pour une méthode de sondage classique avec sous-ensemble statique 14 à environ 10% grâce à l'invention) tout en conservant les critères de précisions demandés aux solutions distribuées pour la participation aux services système. Par rapport à une méthode de sondage classique, il n'est pas nécessaire 5 ici de faire une analyse préalable du parc d'équipements de mesure pour sélectionner un sous-ensemble représentatif. Il n'est pas non plus nécessaire de posséder une connaissance préalable du parc (localisation, type de tarification, typologie de site, etc.), connaissance pas toujours accessible et dont la collecte est chronophage et coûteuse. De plus, la configuration du parc 10 est sujette à évoluer (déménagements, travaux, etc.) et maintenir à jour une telle connaissance nécessite une allocation de ressources permanente. En outre, il n'est pas non plus nécessaire de surveiller la possible défaillance des équipements de mesure sélectionnés dans le sous-ensemble de sondage puisque le sous-ensemble de sondage change à chaque période 15 d'échantillonnage. Chaque équipement de mesure fait partie à un moment ou un autre du sous-ensemble de sondage, ce qui permet d'améliorer la représentativité et par suite la qualité de l'estimation. Dans le cas où les équipements de mesure ont une probabilité d'envoi 20 variable, on augmente en outre l'information REVENDICATIONS1. Procédé d'estimation d'une grandeur physique variable dans le temps par un système de mesure distribué comportant une plateforme centralisée (1) et un ensemble composé d'une pluralité d'équipements (2) de mesure distribués aptes chacun à communiquer à distance avec ladite plateforme centralisée (1) par l'intermédiaire d'un réseau (3) de communication, le procédé comportant une étape (100) de mesure, par chaque équipement (2,) de mesure dudit ensemble, d'une donnée locale représentative de ladite grandeur physique, une étape (120) de collecte, sur une période d'échantillonnage prédéterminée, d'échantillons de mesure provenant d'un sous-ensemble dudit ensemble dans laquelle chaque équipement de mesure appartenant audit sous-ensemble transmet à ladite plateforme centralisée (1) un échantillon de mesure (Pi) de la donnée locale via ledit réseau (3) de communication, et une étape (130) d'estimation de ladite grandeur physique sur la période d'échantillonnage à partir d'une agrégation des échantillons de mesure (Pi) collectés, le procédé étant caractérisé en ce qu'il comporte, à chaque période d'échantillonnage, une phase (110) d'exécution, au niveau de chaque équipement (2,) de mesure dudit ensemble, d'un algorithme (22,) local de décision apte à décider aléatoirement, sur la base d'une probabilité d'envoi (rci) prédéfinie pour chaque équipement (2) de mesure, si l'équipement de mesure correspondant fait partie dudit sous-ensemble, et en ce que ladite grandeur physique est estimée à partir des échantillons de mesure (Pi) collectés et de la probabilité d'envoi (ffi) prédéfinie pour chaque équipement de mesure appartenant audit sous-ensemble. 2. Procédé selon la revendication 1, caractérisé en ce que la probabilité d'envoi (rci) prédéfinie pour au moins un équipement (2,) de mesure dudit ensemble a une valeur variable en fonction de la valeur de l'échantillon (Pi) de mesure effectuée et/ou de conditions de réalisation de la mesure et/ou de l'instant de la mesure, et en ce que si ledit au 16 moins un équipement (2;) de mesure fait partie dudit sous-ensemble, ledit au moins un équipement de mesure transmet la probabilité d'envoi à la plateforme centralisée (1) via ledit réseau (3) de communication pendant l'étape (120) de collecte, en même temps que l'échantillon de 5 mesure (P;) de la donnée locale. 3. Procédé selon la revendication 1, caractérisé en ce que la probabilité d'envoi (n-t) prédéfinie est transmise à la plateforme centralisée (1) par chaque équipement (2;) de mesure appartenant audit sous-ensemble via ledit réseau (3) de communication pendant l'étape (120) de collecte, 10 en même temps que l'échantillon de mesure (P,) de la donnée locale. 4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte en outre une étape de mise à jour de l'algorithme local de décision et/ou de la probabilité d'envoi (n-e) prédéfinie pour au moins un équipement (2;) de mesure dudit ensemble. 15 5. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la grandeur physique globale est estimée par calcul selon la relation suivante dans laquelle N' est le nombre d'équipements de mesure du sous- ensemble, 9 est la valeur de la donnée locale mesurée par un 20 équipement i de mesure, Tri est la probabilité d'envoi prédéfinie pour un équipement i de mesure, et T est la période d'échantillonnage. 6. Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que la période d'échantillonnage est inférieure ou égale à 10 minutes. 25 7. Procédé selon la revendication 6, caractérisé en ce que la période d'échantillonnage est de l'ordre de la seconde ou de la dizaine de secondes. 17 S. Système de mesure distribué pour estimer une grandeur physique variable dans le temps, le système comportant une plateforme centralisée (1) et un ensemble composé d'une pluralité d'équipements 5 (2,) de mesure distribués aptes chacun à communiquer à distance avec ladite plateforme centralisée (1) par l'intermédiaire d'un réseau (3) de communication, le système étant configuré de sorte que, lors d'une étape (100) de mesure, chaque équipement (2,) de mesure dudit ensemble mesure une donnée locale représentative de ladite grandeur 10 physique, lors d'une étape (120) de collecte, chaque équipement de mesure appartenant à un sous-ensemble dudit ensemble transmet à ladite plateforme centralisée (1) un échantillon de mesure (Pi) de la donnée locale via ledit réseau (3) de communication sur une période d'échantillonnage prédéterminée, et lors d'une étape (130) d'estimation, 15 la plateforme centralisée (1) estime ladite grandeur physique sur la période d'échantillonnage à partir d'une agrégation des échantillons de mesure collectés, caractérisé en ce que chaque équipement (2,) de mesure dudit ensemble comporte un module logiciel (22,) configuré pour exécuter, à chaque période d'échantillonnage, un algorithme local 20 de décision apte à décider aléatoirement, sur la base d'une probabilité d'envoi (ff prédéfinie pour chaque équipement (2,) de mesure, si l'équipement de mesure correspondant fait partie dudit sous-ensemble, et en ce que la plateforme centralisée (1) est configurée pour estimer la grandeur physique à partir des échantillons de mesure (PO collectés 25 et de la probabilité d'envoi (Ire) prédéfinie pour chaque équipement de mesure appartenant audit sous-ensemble. 9. Système selon la revendication 8, dans lequel la donnée locale est une puissance électrique soutirée par une charge locale comprenant un ou plusieurs appareils électriques domestiques raccordés à un réseau de 30 distribution d'électricité, et la grandeur physique est la puissance électrique totale consommée.