FR 3141544 A1 20240503 FR 2211153 A 20221026 Procédé d’estimation d’un indicateur collectif de réduction de consommation énergétique pour un ensemble de dispositifs thermiques, dispositif électronique d’estimation et produit programme d’ordinateur associés La présente invention concerne un procédé d’estimation d’un indicateur collectif de réduction de consommation énergétique pour un ensemble de dispositifs thermiques. Chaque dispositif thermique est propre à piloter une température d’un fluide dans un logement. Le procédé étant mis en œuvre par un dispositif électronique d’estimation et comprend les étapes suivantes : obtention d’un modèle thermique de chaque dispositif thermique, acquisition d’une durée d’effacement souhaitée, détermination, pour chaque dispositif thermique, d’un indicateur individuel de potentiel de réduction de consommation pour une pluralité d’instants successifs à partir du modèle thermique, et de la durée d’effacement souhaitée, chaque indicateur individuel étant représentatif de l’énergie propre à être économisée par ledit dispositif thermique pendant une durée inférieure ou égale à la durée d’effacement souhaitée, et pour une pluralité de créneaux, estimation de l’indicateur collectif de réduction de consommation énergétique à partir des indicateurs individuels. Figure pour l'abrégé : figure 1 Procédé d’estimation d’un indicateur collectif de réduction de consommation énergétique pour un ensemble de dispositifs thermiques, dispositif électronique d’estimation et produit programme d’ordinateur associés La présente invention concerne un procédé d’estimation d’un indicateur collectif de réduction de consommation énergétique pour un ensemble de dispositifs thermiques. La présente invention concerne également un dispositif électronique d’estimation d’un indicateur collectif de réduction de consommation énergétique pour un ensemble de dispositifs thermiques. La présente invention concerne également un produit programme d’ordinateur propre à mettre en œuvre un tel procédé. L’invention concerne le domaine du contrôle de la consommation énergétique d’un ensemble de dispositifs thermiques, en particulier d’un ensemble de dispositifs thermiques à usage individuel. On entend par « usage individuel » un usage réalisé par des particuliers. La répartition énergétique à l’échelle d’une société telle qu’une ville, une région, un pays ou un continent souffre actuellement de deux problématiques principales. D’une part la demande en énergie est croissante, notamment en énergie électrique. En effet, de plus en plus de dispositifs de consommation électrique équipent chaque individu ainsi que leur habitation. Cette croissance s’explique également par la croissance démographique. D’autre part, la transition énergétique pousse au développement de nouvelles solutions de production d’énergie électrique telles que les solutions éoliennes, hydroliennes, ou solaires. Bien que ces technologies de production énergétiques présentent de nombreux avantages, chacune présente l’inconvénient d’être une source d’énergie intermittente. Autrement dit, ces sources d’énergie ne produisent pas continuellement la même quantité d’énergie. Ceci conduit à des périodes de surproduction d’énergie, i.e. lorsque la production surpasse la demande, ou à des périodes de productions insuffisantes, i.e. lorsque la demande est supérieure à la production. La deuxième situation est particulièrement problématique puisque tous les individus et infrastructures ne pourront alors pas bénéficier de l’énergie électrique demandée. Dans une telle situation, on connait des systèmes de réseau intelligent (de l’anglais, smart grid ) permettant une allocation préférentielle de l’énergie électrique aux utilisateurs et/ou aux infrastructures considérés comme prioritaires. Toutefois, cette solution est une solution du type tout ou rien. Autrement dit, la demande énergétique de certains est satisfaite, tandis que la demande énergétique d’autres est complètement insatisfaite. En particulier, les logements d’habitation ne sont généralement par considérés comme prioritaires face à d’autres infrastructures telles que les hôpitaux. En parallèle, le développement technologique, et notamment domotique, fournit de plus en plus de dispositifs connectés au sein des habitations. Ces dispositifs ont pour avantage de pouvoir être pilotés à distance, par exemple via une application de téléphone intelligent (de l’anglais, smartphone ). Parmi ces dispositifs, certains sont particulièrement consommateurs d’énergie tels que les radiateurs, les climatiseurs, les ballons d’eau chaude, les chauffe-eaux, les réfrigérateurs et les pompes à chaleur. L’invention vise donc à permettre une réduction intelligente de l’énergie électrique consommée en tirant profit des dispositifs thermiques connectés. À cet effet, l’invention a pour objet un procédé d’estimation d’un indicateur collectif de réduction de consommation énergétique pour un ensemble de dispositifs thermiques, chaque dispositif thermique étant propre à piloter une température d’un fluide dans un logement, le procédé étant mis en œuvre par un dispositif électronique d’estimation, le procédé comprenant les étapes suivantes : obtention d’un modèle thermique de chaque dispositif thermique, acquisition d’une durée d’effacement souhaitée, détermination, pour chaque dispositif thermique, d’un indicateur individuel de potentiel de réduction de consommation pour une pluralité d’instants successifs à partir du modèle thermique, et de la durée d’effacement souhaitée, chaque indicateur individuel étant représentatif de l’énergie propre à être économisée par ledit dispositif thermique pendant une durée inférieure ou égale à la durée d’effacement souhaitée, et pour une pluralité de créneaux, estimation de l’indicateur collectif de réduction de consommation énergétique à partir des indicateurs individuels. Avec le procédé selon l’invention, il est possible de quantifier une consommation énergétique qui pourrait être économisée par les dispositifs thermiques propres à piloter la température d’un fluide dans un logement. On entend par logement une zone propre à recevoir des utilisateurs humains. Le logement est par exemple un logement d’habitation, un bureau d’une entreprise ou un commerce. Le logement d’habitation est par exemple un appartement ou une maison. Selon des modes de réalisation particuliers, le procédé d’estimation comprend une ou plusieurs des caractéristiques suivantes, prise(s) isolément, ou suivant toutes les combinaisons techniquement possibles : l’indicateur collectif de réduction de consommation énergétique comprend un pourcentage de puissance installée de l’ensemble des dispositifs thermiques propre à être économisée lors d’un créneau, dit pourcentage de puissance installée effaçable, l’étape d’acquisition comprend en outre l’acquisition d’un nombre de dispositifs thermiques installés, l'indicateur collectif comprenant en outre, pour chaque créneau, la quantité d’énergie propre à être économisée, ladite quantité étant calculée à partir du nombre de dispositifs thermiques installés, du pourcentage de puissance installée effaçable et d’une puissance installée moyenne des dispositifs thermiques, le dispositif électronique d’estimation est propre à recevoir des données de température et de puissance provenant des dispositifs thermiques, et dans lequel, lors de l'étape d’obtention, le modèle thermique est déterminé à partir des données de température provenant des dispositifs thermiques, l’étape d’obtention du modèle thermique comprend, pour chaque dispositif thermique, les sous-étapes suivantes : mesure de température(s) et de puissance, par le dispositif thermique, à différents instants successifs, la pluralité d’instants formant des séquences de pilotage pendant lesquels le dispositif thermique pilote la température du fluide et des séquences de repos pendant lesquels aucun pilotage de la température n’est effectué, pour chaque instant de chaque séquence de repos, détermination d’un coefficient instantané de déperdition thermique à partir de la ou des températures mesurées par le dispositif thermique audit instant, et à partir d’une durée entre deux instants successifs, pour chaque séquence de repos, calcul d’un coefficient de déperdition thermique de séquence à partir des coefficients instantanés de déperdition thermique de la séquence de repos, calcul d’un coefficient de déperdition thermique de dispositif à partir des coefficients de déperdition thermique de séquence de chaque séquence respectant les contraintes suivantes : une durée de la séquence de repos est supérieure à un premier seuil prédéfini, et une température dans un environnement du dispositif thermique respecte une contrainte prédéfinie, le modèle comprenant chaque coefficient de déperdition thermique de dispositif, l’étape d’acquisition comprend en outre l’acquisition d’un écart maximal souhaité à une température de consigne propre à chaque dispositif thermique, lors de l’étape de détermination, chaque indicateur individuel étant déterminé en fonction de la température de consigne propre au dispositif thermique et de l’écart maximal, chaque indicateur individuel est déterminé en outre à partir du coefficient de déperdition thermique associé à chaque dispositif thermique, chaque indicateur individuel dépendant de l’énergie propre à être économisée en désactivant le dispositif sans qu’une température propre au dispositif thermique ne s’éloigne de la température de consigne de plus de l’écart maximal à la température de consigne, le procédé comprenant en outre, entre les étapes de détermination et d’estimation, une étape de calcul lors de laquelle, une pluralité de quantiles d’indicateurs individuels sont calculés à une pluralité de créneaux, l'étape d’estimation comprenant les sous étapes suivantes qui sont itérées un premier nombre de fois : pour chaque créneau, sélection aléatoire d’un deuxième nombre de quantiles, calcul d’une moyenne des quantiles sélectionnés, l’indicateur collectif de réduction de consommation énergétique comprenant, pour chaque créneau, les moyennes des quantiles calculées, l’étape d’acquisition comprend en outre l’acquisition d’un pourcentage de confiance souhaité, l’indicateur collectif de réduction de consommation énergétique comprenant, pour chaque créneau, le percentile des moyennes des quantiles, calculées selon le pourcentage de confiance souhaité, le procédé comprend, entre l’étape de détermination et l’étape d’estimation, une étape d’association à chaque indicateur individuel de potentiel réduction de consommation, d’un type de journée en fonction d’une valeur moyenne d’une température dans un environnement du dispositif thermique associé lors de la détermination de l’indicateur individuel, l’étape d’acquisition comprend en outre l’acquisition d’un type de journée choisi parmi une pluralité de types de journée, l’étape d’estimation comprenant une sous-étape de filtrage de quantiles d’indicateurs individuels en fonction de la journée type qui leur est associée, l'indicateur collectif étant estimé uniquement à partir des quantiles d’indicateurs individuels associés à un type de journée correspondant au type de journée reçu, chaque dispositif thermique est d’un type choisi dans le groupe consistant en : un radiateur, un climatiseur, un ballon d’eau chaude, un chauffe-eau, un réfrigérateur, et une pompe à chaleur, chaque dispositif thermique étant de préférence du même type, et le procédé comprend en outre une étape d’affichage de l’indicateur collectif sur un afficheur du dispositif électronique d’estimation. La présente invention concerne également un produit programme d’ordinateur comprenant des instructions logicielles qui, lorsqu’elles sont exécutées par un ordinateur, mettent en œuvre le procédé d’estimation tel que décrit ci-dessus. La présente invention concerne en outre un dispositif électronique d’estimation d’un indicateur collectif de réduction de consommation énergétique pour un ensemble de dispositifs thermiques, chaque dispositif thermique étant propre à piloter une température d’un fluide dans un logement, le dispositif d’estimation comprenant : un module d’obtention configuré pour obtenir un modèle thermique de chaque dispositif thermique, un module d’acquisition configuré pour acquérir une durée d’effacement souhaitée, un module de détermination configuré pour déterminer, pour chaque dispositif thermique, un indicateur individuel de potentiel de réduction de consommation pour une pluralité d’instants successifs à partir du modèle thermique, et de la durée d’effacement souhaitée, chaque indicateur individuel étant représentatif de l’énergie propre à être économisée par ledit dispositif thermique pendant une durée inférieure ou égale à la durée d’effacement souhaitée, et un module d’estimation configuré pour estimer, pour une pluralité de créneaux, l’indicateur collectif de réduction de consommation énergétique à partir des indicateurs individuels. Selon un mode de réalisation particulier, le dispositif d’estimation comprend la caractéristique suivante : le module d’acquisition est configuré pour acquérir en outre un écart maximal souhaité à une température de consigne, le module de détermination étant configuré pour déterminer l’indicateur individuel de potentiel de réduction de consommation en outre à partir de l’écart maximal souhaité. Ces caractéristiques et avantages de l’invention apparaitront plus clairement à la lecture de la description qui va suivre, donnée uniquement à titre d’exemple non limitatif et faite en référence aux dessins annexés, sur lesquels : la est une vue schématique d’un système de contrôle énergétique comprenant un dispositif électronique d’estimation selon l’invention ; la est une représentation d’une interface présentée par le dispositif électronique d’estimation de la ; la est une représentation schématique de deux premières grandeurs d’un indicateur collectif de réduction de consommation énergétique obtenu par le dispositif électronique d’estimation représenté sur la ; la est une représentation schématique de deux deuxièmes grandeurs d’un indicateur collectif de réduction de consommation énergétique obtenu par le dispositif électronique d’estimation représenté sur la ; la est une représentation schématique d’un percentile choisi des deux premières grandeurs de l’indicateur collectif de réduction de consommation énergétique de la ; la est une représentation schématique d’un percentile choisi des deux deuxièmes grandeurs de l’indicateur collectif de réduction de consommation énergétique de la ; et la est un organigramme d’un procédé d’estimation d’un indicateur collectif de réduction de consommation énergétique selon l’invention. Sur la , un système de contrôle énergétique 10 est représenté. Le système 10 comprend une pluralité de logements 15 par exemple répartis dans une ville, dans une région, dans un pays, ou sur un continent. Le système de contrôle 10 comprend en outre un dispositif électronique d’estimation 20. Le dispositif électronique d’estimation 20 est destiné à être utilisé par un exploitant du système de contrôle 10 pour contrôler, au cours du temps, une activation ou une désactivation du fonctionnement de dispositif(s) thermique(s) 25 dans les logements 15. Dans chacun des logements 15 est disposé un dispositif thermique 25 respectif propre à piloter la température d’un fluide à l’intérieur du logement 15. On entend par fluide, un gaz ou un liquide. Selon un exemple non représenté, un ou plusieurs logements 15 comprennent plusieurs dispositifs thermiques 25. Chaque dispositif thermique 25 est par exemple un radiateur. Le fluide dont la température est pilotée par le radiateur 25 est l’air à l’intérieur du logement 15. En particulier, chaque dispositif thermique 25 est un radiateur commandable à distance. Par exemple, chaque radiateur est commandable depuis un téléphone intelligent (de l’anglais, smartphone ) via une application telle qu’une application fonctionnant dans un environnement nuagique (de l’anglais, cloud network ). Ainsi, chaque radiateur 25 est apte à ce que son fonctionnement soit activé ou désactivé à distance. Lorsque le radiateur 25 est en fonctionnement activé, il est configuré pour piloter la température T int à l’intérieur du logement 15 pour qu’elle atteigne une température de consigne T c souhaitée. Le dispositif électronique d’estimation 20 est par exemple un ordinateur. Le dispositif d’estimation 20 comprend alors préférentiellement un écran d’affichage 30, aussi appelé moniteur ou afficheur, et une tour 35. La tour 35 est connectée à l’écran d’affichage 30 et propre à diffuser des images sur ledit écran 30. La tour 35 comprend de préférence un processeur 40 et une mémoire 45. De manière connue en soi, le processeur 40 est propre à implémenter des instructions logicielles contenues dans la mémoire 45. La mémoire 45 stocke une pluralité de modules logiciels. En particulier, la mémoire 45 stocke un module 50 d’obtention d’un modèle thermique de chaque dispositif thermique 25, un module 55 d’acquisition de paramètre(s), un module 60 de détermination d’indicateurs individuels de potentiel de réduction de consommation I indiv (t), et préférentiellement un module 65 d’association d’un type de journée à chaque indicateur individuel de potentiel de réduction de consommation I indiv (t), un module 70 de conversion des indicateurs individuels de potentiel de réduction de consommation I indiv (t) pour une journée dite extrême, un module 75 de calcul d’une pluralité de quantiles Q i , un module 80 d’estimation d’un indicateur collectif de réduction de consommation I col , un module 85 d’affichage de l’indicateur collectif I col , et un module 90 de pilotage d’une réduction de consommation. En variante non représentée, chacun desdits modules stockés dans la mémoire 45 est réalisé sous la forme d’un composant de logique programmable, tel qu’un FPGA (de l’anglais, Field Programmable Gate Array ) ou encore d’un circuit intégré, tel qu’un ASIC (de l’anglais, Application Specific Integrated Circuit ). Lorsque le dispositif électronique d’estimation 20 est réalisé sous la forme d’un ou plusieurs logiciels, comme représenté sur la , c’est-à-dire qu’il comprend un programme d’ordinateur, il est en outre apte à être enregistré sur un support lisible par ordinateur, non représenté. Le support lisible par ordinateur est par exemple un medium apte à mémoriser les instructions électroniques et à être couplé à un bus d’un système informatique. À titre d’exemple, le support lisible est un disque optique, un disque magnéto-optique, une mémoire ROM, une mémoire RAM, tout type de mémoire non-volatile (par exemple EPROM, EEPROM, FLASH, NVRAM), une carte magnétique ou une carte optique. Sur le support lisible est alors mémorisé sur un programme d’ordinateur comprenant les instructions logicielles. Le module d’obtention 50 est configuré pour obtenir un modèle thermique de chaque dispositif thermique 25. Le modèle thermique comprend par exemple pour chaque dispositif thermique 25, un coefficient de déperdition thermique associé au dispositif thermique 25 et qui sera décrit ci-après. Le module d’obtention 50 est par exemple configuré pour mesurer, ou obtenir des mesures, à une pluralité d’instants t successifs, de température T int (t) à l’intérieur du logement 15, de température T ext (t) à l’extérieur du logement 15, et de température de consigne T c (t) et de puissance P(t) de chaque dispositif thermique 25. Préférentiellement, la pluralité d’instants successifs t forme des séquences de pilotage pendant lesquelles le dispositif thermique 25 pilote la température T int (t) à l’intérieur du logement 15 pour qu’elle atteigne la température de consigne T c (t), et des séquences de repos pendant lesquelles aucun pilotage de la température T int (t) n’est effectué. Le module d’obtention 50 est par exemple en outre configuré pour déterminer, pour chaque instant t de chaque séquence de repos, un coefficient instantané de déperdition thermique. Plus particulièrement, le module d’obtention 50 est configuré pour déterminer ce coefficient instantané à partir des températures T int (t), T ext (t) mesurées auxdits instants t, et à partir d’une durée D entre deux instants successifs t. La durée D est appelée pas de temps entre deux instants t successifs. Le module d’obtention 50 est par exemple configuré pour appliquer l’équation suivante. Le coefficient instantané représente la variation de température T int (t+D) à l’intérieur du logement 15 pendant la durée D en fonction de la différence de température entre l’intérieur T int (t) et l’extérieur T ext (t) du logement 15, lors d’une séquence de repos. Le module d’obtention 50 est en outre configuré pour calculer, pour chaque séquence de repos, un coefficient de déperdition thermique de séquence. Le module d’obtention 50 est par exemple configuré pour calculer le coefficient de séquence à partir des coefficients instantanés de la séquence de repos correspondante. En particulier, le module d’obtention 50 est configuré pour calculer, pour chaque séquence de repos, le coefficient de séquence comme étant la médiane des coefficients instantanés durant ladite séquence de repos. Le module d’obtention 50 est en outre configuré pour calculer, pour chaque dispositif thermique 25, un coefficient de déperdition thermique de dispositif 25. À cet effet, le module d’obtention 50 est par exemple configuré pour sélectionner, parmi les séquences de repos, uniquement celles ayant une durée supérieure à un premier seuil prédéfini et pour lesquelles la température extérieure T ext (t) respecte une contrainte prédéfinie. Le premier seuil prédéfini est par exemple égal à une heure. La contrainte prédéfinie sur la température extérieure T ext (t) est par exemple que la température extérieure moyenne T ext pendant la séquence soit inférieure ou égale à 15°C. Le module d’obtention 50 est configuré pour alors calculer le coefficient de déperdition thermique du dispositif 25 à partir des coefficients de déperdition thermique de séquence correspondants aux séquences sélectionnées. Par exemple, le module d’obtention 50 est configuré pour calculer le coefficient du dispositif comme étant la médiane des coefficients des séquences sélectionnées. En complément facultatif, le module d’obtention 50 est configuré pour calculer le coefficient de déperdition thermique de dispositif 25, uniquement si le nombre de séquences de repos du dispositif thermique 25 respectant les contraintes de durée et de température extérieure T ext est supérieur à un deuxième seuil prédéfini, par exemple égal à 5. On comprend alors que le coefficient de déperdition thermique de dispositif représente la variation médiane de température intérieure T int (t) par unité de temps et en fonction de la température extérieure T ext (t), lorsque le dispositif thermique 25 ne pilote pas ladite température intérieure T int (t). Le module d’acquisition 55 est configuré pour acquérir une durée d’effacement souhaitée et une puissance installée de chaque dispositif thermique 25. La puissance installée de chaque dispositif thermique 25 est la puissance maximale que chaque dispositif thermique 25 est propre à déployer. On entend par « effacement » le fait de suspendre le pilotage de la température intérieure T int d’un ou plusieurs dispositifs thermiques 25. Durant un effacement, la température de consigne T c (t) n’est pas nécessairement respectée. La durée d’effacement souhaitée correspond à la durée maximale pendant laquelle le pilotage de la température intérieure T int de chaque dispositif thermique 25 cesseraient, alors que la température de consigne n’est pas nécessairement atteinte. Autrement dit, en l’absence d’effacement, durant ladite durée, chaque dispositif thermique 25 piloterait la température T int en activant le dispositif thermique 25. La durée d’effacement est par exemple égale à 30 minutes ou 1 heure. En complément facultatif, le module d’acquisition 55 est configuré pour acquérir en outre le nombre de dispositifs thermiques 25 dans le système de contrôle énergétique 10, un écart maximal souhaité à la température de consigne T c (t) pendant la durée d’effacement , un type de journée, et un pourcentage de confiance souhaité qui sera décrit ci-après. L’écart maximal souhaité est un écart maximal acceptable à la température de consigne T c donnée à chaque dispositif thermique 25. Ainsi, si cet écart maximal souhaité est par exemple égal à 0,5°C, alors uniquement des effacements pour lesquels la température T int dans l’environnement 15 ne s’éloigne pas de la température de consigne T c (t) de plus de l’écart maximal souhaité sont acceptables. Autrement dit, lors d’un effacement, la température intérieure T int doit respecter à tout instant t l’équation suivante : L’écart maximal souhaité est par exemple égal à 0,5°C, 0,75°C, 1°C, 1,5°C, 2°C, ou 3°C. On entend par « type de journée », une segmentation des journées en fonction de la température extérieure T ext . Par exemple, l’ensemble des journées de chaque dispositif thermique 25 est segmenté en au moins deux types de journées, préférentiellement de deux à cinq types de journées, plus préférentiellement cinq types de journées, respectivement nommées : journée chaude, journée tempérée, journée froide, journée très froide, et journée extrême. Une journée chaude est une journée pour laquelle la température extérieure T ext moyenne entre 7h et 15h et entre 18h et 20h, est supérieure à 15°C. Une journée tempérée est une journée pour laquelle la température extérieure moyenne T ext est comprise, pour la même période de temps, entre 10°C et 15°C. Une journée froide est une journée pour laquelle la température extérieure T ext moyenne est comprise, pour la même période de temps, entre 5°C et 10°C. Une journée très froide est une journée pour laquelle la température extérieure moyenne T ext , pour la même période de temps, est inférieure à 5°C. Une journée extrême est une journée étalon définie par une norme, par exemple dans les règles du mécanisme de capacité de RTE. Le pourcentage de confiance sera détaillé ci-après. Préférentiellement, le module d’acquisition 55 est configuré pour acquérir les paramètres précités, depuis un utilisateur du dispositif d’estimation 20, par exemple via des périphériques non-représentés sur la . Par exemple, le dispositif d’estimation 20 est configuré pour afficher sur l’afficheur 30, l’interface présentée sur la figure 2. Ainsi, l’utilisateur du dispositif d’estimation 20 interagit avec le dispositif d’estimation 20 pour renseigner les paramètres, i.e. le nombre de dispositifs thermiques 25, le type de journée, la durée d’effacement , l’écart maximal souhaité et le pourcentage de confiance . Ceci est fait par exemple via des boutons et/ou des menus déroulants. L’interface représentée sur la , indique par exemple en outre un tutoriel présentant la manière d’interagir avec cette interface. En particulier, lorsque les dispositifs thermiques 25 sont des radiateurs, ce tutoriel indique que les informations complémentaires suivantes sur les entrées à renseigner : Nombre de radiateurs : le nombre de radiateurs dans votre portefeuille, Type de journée : définie à partir de la température extérieure moyenne sur la journée entre les horaires suivants : 7h-15 et 18h-20h et à partir de la journée à températures extrême selon les règles de mécanisme de capacité RTE, Durée d’effacement : durée maximum de l’effacement ciblé, Ecart max à la TC de consigne : contrainte d’écart entre la température intérieure et la température de consigne, lors d’un effacement. Si la contrainte est atteinte, l’effacement s’arrête pour le radiateur en question avant la fin du créneau d’effacement, et % de confiance : pourcentage de chance que l’output atteigne une certaine valeur. Le tutoriel met en avant la présence d’un bouton « Run all » et indique que l’utilisateur doit cliquer dessus. Enfin, le tutoriel précise la manière dont l’indicateur collectif I col doit être interprété. En particulier, le tutoriel indique que les sorties suivantes sont obtenues pour chaque créneau t h , comme décrit ci-après : le pourcentage de la puissance installée qui est effaçable (en %), le facteur de surdimensionnement (l’inverse de la part de la puissance installée qui est effaçable), le volume d’énergie effaçable (en KWh), le capacité moyenne effaçable (en KW). Le module de détermination 60 est configuré pour déterminer, pour chaque dispositif thermique 25, un indicateur individuel de potentiel de réduction de consommation I indiv (t) à la pluralité d’instants t successifs. L’indicateur individuel de potentiel de réduction de consommation I indiv (t) est aussi appelé indicateur individuel dans la suite de la description. Pour cela, le module de détermination 60 est configuré pour déterminer chaque indicateur individuel I indiv (t) à partir du modèle thermique et de la durée d’effacement souhaitée . Préférentiellement, le module de détermination 60 détermine ledit indicateur individuel I indiv (t) en outre à partir de l’écart maximal souhaité . A cet effet, le module de détermination 60 est configuré pour calculer, à chaque instant t, l’énergie consommée par chaque dispositif thermique 25 entre l’instant t et l’instant successif t+D, i.e. à chaque pas de temps D. Préférentiellement, le module de détermination 60 est configuré pour calculer cette énergie à partir de la puissance P(t) du dispositif thermique 25 et du pas de temps D, par exemple selon l’équation suivante. Le module de détermination 60 est en outre configuré pour, à chaque instant t , simuler l’évolution de la température intérieure T int (t) en l’absence de pilotage de la température intérieure T int (t) par le dispositif thermique 25. Par exemple, l’évolution de la température suit l’équation suivante. où t est l’instant auquel est calculé l’évolution de la température intérieure T int (t), et est une durée inférieure à la durée d’effacement souhaitée . Le module de détermination 60 est configuré pour déterminer s’il existe une durée au bout de laquelle la température intérieure T int (t+ ) s’écarte de la température de consigne T c (t) d’une valeur égale à l’écart maximal souhaité . Autrement dit, le module de détermination 60 est configuré pour déterminer s’il existe une durée respectant l’équation suivante : Le module de détermination 60 est configuré pour, s’il existe une telle durée alors déterminer cette durée comme la durée effective d’effacement. Autrement dit, l’effacement n’aura lieu que pendant la durée , si la durée est inférieure ou égale à la durée d’effacement souhaitée , sinon il aura lieu pendant la durée d’effacement souhaitée . Ainsi, partant de l’instant t, l’effacement aura lieu à tous instants successifs inférieurs à , la durée étant par exemple différente d’un dispositif thermique 25 à un autre. Le module de détermination 60 est configuré pour déterminer, pour chaque instant t i , t i+ 1 dans la durée d’effacement effective, l’énergie consommée par le dispositif thermique 25. Pour chaque instant t et pour chaque dispositif thermique 25, le module de détermination 60 est configuré pour calculer un pourcentage de puissance installée de l’ensemble des dispositifs thermiques 25 propre à être économisé lors d’un créneau t h , aussi appelé pourcentage de puissance installée effaçable comme étant égal au rapport entre la puissance effaçable et la puissance qui aurait été consommée si le dispositif thermique 25 fonctionnait à pleine puissance, i.e. à la puissance installée durant la durée d’effacement souhaitée , par exemple selon l’équation suivante où est égal à zéro si la température intérieure T int (t) s’éloigne de la température de consigne T c (t) de plus que de l’écart maximal souhaité , et respecte l’équation 3, et est la puissance installée du dispositif thermique 25. Un créneau t h est un créneau horaire délimitant une durée entre deux instants choisis. Tous les créneaux t h ont par exemple une même durée. Le module de détermination 60 est configuré pour ensuite déterminer, pour chaque instant t et pour chaque dispositif thermique 25, un indicateur individuel respectif I indiv (t), par exemple comme étant égal au pourcentage de puissance effaçable , respectant ainsi l’équation suivante. L’indicateur individuel I indiv (t) représente donc le rapport entre l’énergie consommée effaçable pendant la durée d’effacement souhaitée et l’énergie qui aurait été consommée si le dispositif thermique 25 fonctionnait à pleine puissance. Le module d’association 65 est configuré pour associer, à chaque indicateur individuel I indiv (t), un type de journée parmi les types de journée précités. Préférentiellement, le module d’association 65 est configuré, pour associer, à chaque indicateur individuel I indiv (t), un type de journée parmi une journée chaude, une journée tempérée, une journée froide, ou une journée très froide. Le module d’association 65 est configuré pour associer le type de journée à chaque indicateur individuel I indiv (t) en fonction de la température extérieure moyenne de la journée associée. La température extérieure moyenne est par exemple calculée entre 7h et 15h et entre 18h et 20h, pour le dispositif thermique 25 en question. Le module de conversion 70 est configuré pour convertir chaque puissance effaçable en une puissance effaçable extrême ( , afin de déterminer un indicateur individuel I indiv (t) correspondant à un indicateur de référence associé à la journée de type extrême. À cet effet, le module de conversion 70 est configuré pour appliquer les équations suivantes. où est la puissance effaçable ) par un dispositif thermique 25 pendant un créneau à un jour , est la puissance effaçable ) par le même dispositif thermique 25 pendant le même créneau à un jour une semaine avant le jour , T extrême (t) est une température extrême représentant une chronique de 48 valeurs de températures, par pas de temps demi-horaire, telle que présentée dans l’Annexe 1.1.2 des règles du mécanisme de capacité de RTE, TFLS(t) est une température correspondant à une moyenne pondérée de températures demi-horaires de 32 stations météo dont le résultat est ensuite lissé, tel que disponible sur le site internet https://data.enedis.fr/explore/dataset/donnees-de-temperature-et-de-pseudo-rayonnement/information/, et représente une chronique de valeurs en Watt par degré Celsius à un pas demi−horaire pour chaque dispositif thermique 25, associée à un couple de deux chroniques, la première chronique étant une chronique de puissances à un pas demi−horaire et la deuxième chronique étant une chronique de Température France Lissée (TFLS), tel qu’expliqué en Annexe E des règles du mécanisme de capacité de RTE. Le module de conversion 70 est configuré pour déterminer l’indicateur I indiv (t) de référence correspondant à la journée de type extrême, comme étant égal au pourcentage de puissance effaçable extrême ( , par exemple selon l’équation suivante : Le module de calcul 75 est configuré pour calculer, à partir de tous les indicateurs individuels déterminés I indiv (t), et pour chaque type de journée, des quantiles Q I (t h ) d’indicateur individuel. Autrement dit, le module de calcul 75 est configuré pour calculer, pour une pluralité de créneaux t h , quel pourcentage de la puissance installée est au moins économisable par 50% des dispositifs thermique 25, par 20% des dispositifs thermiques 25, par 80% des dispositifs thermiques 25, etc. Les créneaux t h sont espacés les uns des autres d’une durée égale à la durée d’effacement souhaitée . Ainsi, le créneau t h commence à un instant et s’étend jusqu’à l’instant . Par exemple, si la durée d’effacement souhaitée est égale à 30 min, alors 48 créneaux t h sont considérés. Le module d’estimation 80 est configuré pour estimer, pour chaque créneau t h , l’indicateur collectif I col (t h ) de réduction de consommation énergétique à partir des indicateurs individuels I indiv (t). À cet effet, le module d’estimation 80 est configuré pour filtrer, les indicateurs individuels I indiv (t) associés au type de journée acquis par le module d’acquisition 55. Le module d’estimation 80 est configuré pour ensuite itérer un premier nombre de fois les actions suivantes. Pour chacun du premier nombre de fois, le module d’estimation 80 est configuré pour sélectionner, pour chaque créneau t h , aléatoirement un deuxième nombre de quantiles Q I (t h ), par exemple égal au nombre de dispositifs thermiques 25, parmi les quantiles Q I (t h ) correspondant aux indicateurs individuels filtrés. Toujours pour chacun du premier nombre de fois, le module d’estimation 80 est configuré pour calculer une moyenne des quantiles Q I (t h ) sélectionnés. Ainsi, le module d’estimation 80 est configuré pour obtenir, pour chaque créneau t h , un nombre de quantiles Q I (t h ) moyennés égal au premier nombre. À titre d’exemple, le premier nombre est égal à 1000 et le deuxième nombre est égal à 500. Ainsi, le module d’estimation 80 est configuré pour calculer, pour chaque créneau t h , 1000 quantiles Q I (t h ) moyennés, chacun obtenu en moyennant 500 quantiles Q I (t h ) sélectionnés aléatoirement parmi les quantiles Q i (t h ) correspondant aux indicateurs individuels I indiv (t) filtrés. Le module d’estimation 80 est ensuite configuré pour ordonner, pour chaque créneau t h , les moyennes des quantiles Q i (t h ) calculés, par exemple du plus petit au plus grand. Ainsi, le module d’estimation 80 est configuré pour obtenir, pour chaque créneau t h , une distribution de quantiles Q i (t h ). La distribution de quantiles Q i (t h ) forme alors, pour chaque créneau t h , une distribution du pourcentage % eff (t h ) de la puissance installée effaçable pendant le créneau t h . Le module d’estimation 80 est configuré pour estimer les percentiles des quantiles ordonnés pour chaque créneau Perc(% eff (t h )). En complément facultatif, le module d’estimation 80 est en outre configuré pour calculer, pour chaque créneau t h et pour chaque moyenne de quantile, un facteur de surdimensionnement dont chaque percentile respecte l’équation suivante. où est la fonction percentile. De même, le module d’estimation 80 est par exemple configuré pour calculer, pour chaque créneau t h , une quantité d’énergie effaçable , aussi appelé volume d’énergie effaçable, dont chaque percentile respecte l’équation suivante. : où est la puissance installée moyenne des dispositifs thermiques 25 dans le système 10, et est le nombre de dispositifs thermiques 25 dans le système 10. En outre, le module d’estimation 80 est configuré pour calculer, pour chaque créneau t h , une capacité moyenne effaçable , par exemple dont les percentiles vérifient l’équation suivante. L’indicateur collectif I col (t h ) comprend avantageusement, pour chaque créneau t h , les pourcentages de puissance installée propres à être économisés calculés par le module d’estimation 80. Préférentiellement, l’indicateur collectif comprend en outre, pour chaque créneau t h , les facteurs de surdimensionnement , les quantités d’énergie effaçables , et les capacités moyennes effaçables . Le module d’affichage 85 est configuré pour afficher, par exemple sur l’afficheur 30, l’indicateur collectif I col (t h ) obtenu par le module d’estimation 80. Sur la partie haute de la figure 3, la distribution d’un pourcentage % eff (t h ) de chaque créneau t h est représentée. Ceci correspond à un exemple d’affichage de l’indicateur collectif I col (t h ) obtenu par le module d’estimation 80. Sur la partie basse de la figure 3, la distribution de facteurs de surdimensionnement ) est représentée. Sur la partie haute de la figure 4, la distribution de quantités d’énergie effaçables est représentée. Sur la partie basse de la figure 4, la distribution de capacités moyennes effaçables , est représentée. En complément facultatif, lorsque le pourcentage de confiance est acquis par le module d’acquisition 55, le module d’affichage 85 est configuré pour alors n’afficher que le percentile correspondant au pourcentage de confiance dans la distribution d’un pourcentage % eff (t h ). Il en va de même pour les distributions de facteur de surdimensionnement , de quantité d’énergie effaçable , et de capacité moyenne effaçable . Sur la partie haute de la figure 5 le percentile de pourcentage effaçable % eff (t h ) correspondant à un pourcentage de confiance de 90% est représenté. Sur la partie basse de la figure 5, le percentile de facteur de surdimensionnement correspondant à un pourcentage de confiance de 90% est représenté. Sur la partie haute de la figure 6, le percentile de quantité d’énergie effaçable correspondant à un pourcentage de confiance de 90% est représenté. Sur la partie basse de la figure 6, le percentile de capacité moyenne effaçable correspondant à un pourcentage de confiance de 90% est représenté. Le module de pilotage 90 est configuré pour commander le pilotage de la stratégie d’effacement. Autrement dit, le module de pilotage 90 est configuré pour envoyer une instruction d’arrêt du pilotage de la température intérieure T int (t) au(x) dispositif(s) thermique(s) 25 en fonction de l’effacement souhaité. Le fonctionnement du système de contrôle énergétique 10 et plus particulièrement du dispositif électronique d’estimation 20 va maintenant être décrit en référence à la illustrant un organigramme d’un procédé d’estimation de l’indicateur collectif I col (t h ) de réduction de consommation, selon l’invention. Lors d’une étape d’obtention 110, le module d’obtention 50 obtient le modèle thermique de chaque dispositif thermique 25. Pour cela, l’étape d’obtention 110 comprend préférentiellement une sous-étape de mesure 112 lors de laquelle le module d’obtention 50 acquiert à la pluralité d’instants t, les mesures de température intérieure T int (t), de température extérieure T ext (t), et avantageusement les mesures de température de consigne T c (t) et de puissance P(t), depuis chaque dispositif thermique 25. Puis, l’étape d’obtention 110 comprend une sous-étape de détermination 114, lors de laquelle le module d’obtention 55 détermine, pour chaque instant t de chaque séquence de repos, le coefficient instantané de déperdition thermique , par exemple à partir des températures mesurées T int (t), T ext (t), et à partir d’une durée D entre deux instants t successifs. Pour cela, le module d’obtention 50 applique par exemple l’équation 1. Puis, l’étape d’obtention 110 comprend de préférence une première sous-étape de calcul 116 lors de laquelle le module d’obtention 50 calcule le coefficient de déperdition thermique de séquence à partir des coefficients instantanés de déperdition thermique de la séquence de repos, par exemple par calcul d’une médiane des coefficients instantanés . Ensuite, l’étape d’obtention 110 comprend une deuxième sous-étape de calcul 118 lors de laquelle le module d’obtention 50 calcule le coefficient de déperdition thermique du dispositif 25 à partir des coefficients de déperdition thermique de séquence de chaque séquence respectant les contraintes précitées, par exemple par calcul d’une médiane des coefficients de séquence . Lors d’une étape d’acquisition 120, le module d’acquisition 55 acquiert la durée d’effacement souhaitée . Préférentiellement, lors de l’étape d’acquisition 120, le module d’acquisition 55 acquiert également le nombre de dispositifs thermiques 25, le type de journée, l’écart maximal souhaité à la température de consigne T c (t) et optionnellement le pourcentage de confiance . Ensuite, lors d’une étape de détermination 130, le module de détermination 60 détermine pour chaque instant t auquel les températures T int (t), T ext (t) ont été acquises lors de l’étape d’obtention 110, l’indicateur individuel I indiv (t) propre à chaque dispositif. A cet effet, le module de détermination 60 détermine ces indicateurs individuels I indiv (t) comme expliqué précédemment, et notamment en appliquant les équations 3 à 7. Puis, lors d’une étape d’association 140, le module d’association 65 associe à chaque indicateur individuel I indiv (t), un type de journée, par exemple en fonction d’une valeur moyenne d’une température dans un environnement du dispositif thermique 25. Ladite valeur moyenne est par exemple la température extérieure T ext moyenne pendant la journée à laquelle l’indicateur individuel I indiv (t) est déterminé. La température extérieure moyenne T ext est par exemple calculée entre 7h et 15h et entre 18h et 20h, pour le dispositif thermique 25 en question. Lors d’une étape de conversion 145, le module de conversion 70 calcule l’indicateur individuel de référence associé aux journées extrême en convertissant chaque puissance effaçable P eff en la puissance effaçable extrême correspondante . Pour cela, le module de conversion 70 applique par exemple les équations 8 et 9. Lors d’une étape de calcul 150, le module de calcul 75 calcule les quantiles Q I (t h ) d’indicateurs individuels I indiv (t) à la pluralité de créneaux t h comme expliqué précédemment. Lors d’une étape d’estimation 160, le module d’estimation 80 estime pour chaque créneau t h , l’indicateur collectif I col (t h ). Pour cela, l’étape d’estimation 160 comprend préférentiellement une sous-étape de filtrage 162 lors de laquelle le module d’estimation 80 filtre les indicateurs individuels I indiv (t) en fonction de la journée type qui leur est associée. Le module d’estimation 80 sélectionne, pour chaque créneau t h , uniquement les quantiles Q i (t h ) respectifs des indicateurs individuels I indiv (t) associés au type de journée acquis pendant l’étape d’acquisition 120. L’étape d’estimation 160 comprend en outre une sous-étape de sélection 164 et une troisième sous-étape de calcul 165, qui sont itérées le premier nombre de fois. Lors de la sous-étape de sélection 164, le module d’estimation 80 sélectionne aléatoirement le deuxième nombre de quantiles Q i (t h ). Lors de la troisième sous-étape de calcul 165, le module d’estimation 80 calcule la moyenne des quantiles Q i (t h ) sélectionnés. L’étape d’estimation 160 comprend en outre une quatrième sous-étape de calcul 166, lors de laquelle le module d’estimation 80 ordonne, pour chaque créneau t h les moyennes des quantiles Q i (t h ) calculés. La distribution de quantiles Q i (t h ) forme alors, pour chaque créneau t h , la distribution de pourcentages % eff (t h ) de la puissance installée effaçable pour chaque créneau t h . Lors de la quatrième sous-étape de calcul 166, le module d’estimation 80 calcule préférentiellement en outre le facteur de surdimensionnement , la quantité d’énergie effaçable , et la capacité effaçable , par exemple à l’aide des équations 10 à 12. L’indicateur collectif I col (t h ) comprend le pourcentage de la puissance installée propre à être économisée % eff (t h ), et préférentiellement le facteur de surdimensionnement , la quantité d’énergie effaçable , et la capacité moyenne effaçable . Lors d’une étape d’affichage 170, le module d’affichage 85 affiche l’indicateur collectif I col (t h ) pour chaque créneau t h , par exemple sur l’afficheur 30. En complément facultatif, lorsqu’un pourcentage de confiance est acquis lors de l’étape d’acquisition 120, le module d’affichage 85 n’affiche que le percentile correspondant au pourcentage de confiance dans la distribution de pourcentages % eff (t h ), et avantageusement les facteurs de surdimensionnement , quantités d’énergie effaçables , et capacités moyennes effaçables . Lors d’une étape de pilotage 180, le module de pilotage 90 pilote la stratégie d’effacement comme décrit précédemment. Avantageusement, le module d’acquisition est propre à afficher sur l’écran d’affichage 30 une fenêtre d’interface utilisateur présentant un champ permettant à un utilisateur de saisir la durée d’effacement souhaitée , comme représenté sur la figure 2. Préférentiellement, la fenêtre d’interface présente des champs permettant aussi à l’utilisateur de saisir le nombre de dispositifs thermiques 25, le type de journée, l’écart maximal souhaité à la température de consigne T c (t) et le pourcentage de confiance . Ceci permet à un utilisateur de mettre en œuvre facilement l’étape d’acquisition 120 et donc de faire varier les paramètres d’entrée pour la mise en œuvre des étapes 130 à 170. Selon une variante, chaque dispositif thermique 25 est par exemple un climatiseur, un ballon d’eau chaude, un chauffe-eau, un réfrigérateur, ou une pompe à chaleur. Préférentiellement, chaque dispositif thermique 25 est du même type. Autrement dit, chaque dispositif thermique 25 est un climatiseur, chaque dispositif thermique 25 est un ballon d’eau chaude, chaque dispositif 25 est un chauffe-eau, chaque dispositif 25 est un réfrigérateur, ou chaque dispositif 25 est une pompe à chaleur. Le dispositif d’estimation 20 est analogue à celui présenté précédemment. Le fonctionnement du dispositif d’estimation 20 est analogue aux différences près listées ci-après. Lorsque le dispositif thermique 25 est un climatiseur, la température intérieure T int (t) est supérieure à la température extérieure T ext (t). Le coefficient est calculé de manière analogue mais quantifie le réchauffement du logement 15 par l’air à l’extérieur du logement 15. Ainsi, l’équation 2 devient et l’équation 5 devient Lorsque le dispositif thermique 25 est un ballon d’eau chaude, un chauffe-eau ou une pompe à chaleur, le fluide est de l’eau à l’intérieur d’un réservoir du dispositif thermique 25. La température intérieure T int (t) est la température de l’eau dans le réservoir. La température extérieure T ext (t) est la température dans le logement 15. Enfin, lorsque le dispositif thermique 25 est un réfrigérateur, le fluide est l’air dans une enceinte du réfrigérateur. La température intérieure T int (t) est la température dans le réfrigérateur. La température extérieure T ext (t) est la température dans le logement 15. Dans cet exemple la température intérieure T int (t) est plus faible que la température extérieure T ext (t). Ainsi, les équations 2 et 5 sont respectivement remplacées par les équations 13 et 14. Ainsi, le procédé selon l’invention permet d’estimer le pourcentage de la puissance installée qui pourrait être économisé en appliquant une stratégie d’effacement. Cette information permet, en elle-même, de pallier le déséquilibre du réseau électrique lors des jours de pointe, sans avoir recours à des moyens de production thermiques très polluants et coûteux. Ainsi, l’invention permet de limiter la production de gaz à effet de serre tel que le dioxyde de carbone en rationalisant l’utilisation de dispositifs thermiques. Procédé d’estimation d’un indicateur collectif (I col (t h )) de réduction de consommation énergétique pour un ensemble de dispositifs thermiques (25), chaque dispositif thermique (25) étant propre à piloter une température d’un fluide dans un logement (15), le procédé étant mis en œuvre par un dispositif électronique d’estimation (20), le procédé comprenant les étapes suivantes : obtention (110) d’un modèle thermique ( ) de chaque dispositif thermique (25), acquisition (120) d’une durée d’effacement souhaitée ( ), détermination (130), pour chaque dispositif thermique (25), d’un indicateur individuel (I indiv (t)) de potentiel de réduction de consommation pour une pluralité d’instants successifs (t) à partir du modèle thermique ( ) , et de la durée d’effacement souhaitée ( ), chaque indicateur individuel (I indiv (t)) étant représentatif de l’énergie propre à être économisée par ledit dispositif thermique (25) pendant une durée inférieure ou égale à la durée d’effacement souhaitée ( ), et pour une pluralité de créneaux (t h ), estimation (160) de l’indicateur collectif (I col (t h )) de réduction de consommation énergétique à partir des indicateurs individuels (I indiv (t)). Procédé selon la revendication 1, dans lequel l’indicateur collectif (I col (t h )) de réduction de consommation énergétique comprend un pourcentage de puissance installée (% eff (t h )) de l’ensemble des dispositifs thermiques (25) propre à être économisée lors d’un créneau (t h ), dit pourcentage de puissance installée effaçable. Procédé selon la revendication 2, dans lequel l’étape d’acquisition (120) comprend en outre l’acquisition d’un nombre ( ) de dispositifs thermiques (25) installés, l'indicateur collectif (I col (t h )) comprenant en outre, pour chaque créneau (t h ), la quantité d’énergie propre à être économisée ( ), ladite quantité ( ) étant calculée à partir du nombre ( ) de dispositifs thermiques (25) installés, du pourcentage de puissance installée effaçable ( ) et d’une puissance installée moyenne ( ) des dispositifs thermiques (25). Procédé selon l’une quelconque des revendications précédentes, dans lequel le dispositif électronique d’estimation (20) est propre à recevoir des données de température (T int (t), T ext (t), T c (t)) et de puissance (P(t)) provenant des dispositifs thermiques (25), et dans lequel, lors de l'étape d’obtention (110), le modèle thermique ( ) est déterminé à partir des données de température provenant des dispositifs thermiques (25). Procédé selon la revendication précédente, dans lequel l’étape (110) d’obtention du modèle thermique ( ) comprend, pour chaque dispositif thermique (25), les sous-étapes suivantes : mesure (112) de température(s) (T int (t), T ext (t), T c (t)) et de puissance (P(t)), par le dispositif thermique (25), à différents instants successifs (t), la pluralité d’instants (t) formant des séquences de pilotage pendant lesquels le dispositif thermique (25) pilote la température du fluide et des séquences de repos pendant lesquels aucun pilotage de la température n’est effectué, pour chaque instant (t) de chaque séquence de repos, détermination (114) d’un coefficient instantané de déperdition thermique ( ) à partir de la ou des températures (T int (t), T ext (t)) mesurées par le dispositif thermique (25) audit instant (t), et à partir d’une durée (D) entre deux instants successifs (t), pour chaque séquence de repos, calcul (116) d’un coefficient de déperdition thermique de séquence à partir des coefficients instantanés de déperdition thermique ( ) de la séquence de repos, calcul (118) d’un coefficient de déperdition thermique de dispositif ( ) à partir des coefficients de déperdition thermique de séquence de chaque séquence respectant les contraintes suivantes : une durée de la séquence de repos est supérieure à un premier seuil prédéfini, et une température (T ext (t)) dans un environnement du dispositif thermique (25) respecte une contrainte prédéfinie, le modèle comprenant chaque coefficient de déperdition thermique de dispositif ( ). Procédé selon l’une quelconque des revendications précédentes, dans lequel l’étape d’acquisition (120) comprend en outre l’acquisition d’un écart maximal souhaité ( ) à une température de consigne (T c (t)) propre à chaque dispositif thermique (25), lors de l’étape de détermination (130), chaque indicateur individuel (I indiv (t)) étant déterminé en fonction de la température de consigne (T c (t)) propre au dispositif thermique (25) et de l’écart maximal ( ). Procédé selon les revendications 5 et 6, dans lequel chaque indicateur individuel (I indiv (t)) est déterminé en outre à partir du coefficient ( ) de déperdition thermique associé à chaque dispositif thermique (25), chaque indicateur individuel (I indiv (t)) dépendant de l’énergie propre à être économisée en désactivant le dispositif sans qu’une température (T int (t)) propre au dispositif thermique (25) ne s’éloigne de la température de consigne (T c (t)) de plus de l’écart maximal ( ) à la température de consigne (T c (t)). Procédé selon la revendication précédente, comprenant en outre, entre les étapes de détermination (130) et d’estimation (160), une étape de calcul (150) lors de laquelle, une pluralité de quantiles (Q i (t h )) d’indicateurs individuels (I indiv (t)) sont calculés à une pluralité de créneaux (t h ), l'étape d’estimation (160) comprenant les sous étapes suivantes qui sont itérées un premier nombre de fois : pour chaque créneau, sélection (164) aléatoire d’un deuxième nombre de quantiles (Q i (t h )), calcul (165) d’une moyenne des quantiles (Q i (t h )) sélectionnés, l’indicateur collectif (I col (t h )) de réduction de consommation énergétique comprenant, pour chaque créneau, les moyennes des quantiles (Q i (t h )) calculées. Procédé selon la revendication précédente, dans lequel l’étape d’acquisition (120) comprend en outre l’acquisition d’un pourcentage de confiance souhaité ( ), l’indicateur collectif (I col (t h )) de réduction de consommation énergétique comprenant, pour chaque créneau (t h ), le percentile des moyennes des quantiles (Q i (t h )), calculées selon le pourcentage de confiance souhaité ( ). Procédé selon la revendication 9, dans lequel l’étape d’acquisition (120) comprend en outre l’acquisition d’un type de journée choisi parmi une pluralité de types de journée, l’étape d’estimation (160) comprenant une sous-étape de filtrage (162) de quantiles (Q i (t h )) d’indicateurs individuels en fonction de la journée type qui leur est associée, l'indicateur collectif (I col (t h )) étant estimé uniquement à partir des quantiles (Q i (t h )) d’indicateurs individuels associés à un type de journée correspondant au type de journée reçu. Procédé selon l’une quelconque des revendications précédentes, dans lequel le procédé comprend, entre l’étape de détermination (130) et l’étape d’estimation (160), une étape d’association (140) à chaque indicateur individuel (I indiv (t)) de potentiel réduction de consommation, d’un type de journée en fonction d’une valeur moyenne d’une température (T ext (t)) dans un environnement du dispositif thermique (25) associé lors de la détermination (140) de l’indicateur individuel (I indiv (t)). Procédé selon l’une quelconque des revendications précédentes, dans lequel chaque dispositif thermique (25) est d’un type choisi dans le groupe consistant en : un radiateur, un climatiseur, un ballon d’eau chaude, un chauffe-eau, un réfrigérateur, et une pompe à chaleur, chaque dispositif thermique (25) étant de préférence du même type. Procédé selon l’une quelconque des revendications précédentes, dans lequel le procédé comprend en outre une étape (170) d’affichage de l’indicateur collectif sur un afficheur (30) du dispositif électronique d’estimation (20). Produit programme d’ordinateur comprenant des instructions logicielles qui, lorsqu’elles sont exécutées par un ordinateur, mettent en œuvre le procédé d’estimation selon l’une quelconque des revendications précédentes. Dispositif (20) d’estimation d’un indicateur collectif (I col (t h )) de réduction de consommation énergétique pour un ensemble de dispositifs thermiques (25), chaque dispositif thermique (25) étant propre à piloter une température d’un fluide dans un logement (15), le dispositif d’estimation (20) comprenant : un module d’obtention (50) configuré pour obtenir un modèle thermique ( ) de chaque dispositif thermique (25), un module d’acquisition (55) configuré pour acquérir une durée d’effacement souhaitée ( ), un module de détermination (60) configuré pour déterminer, pour chaque dispositif thermique (25), un indicateur individuel (I indiv (t)) de potentiel de réduction de consommation pour une pluralité d’instants successifs (t) à partir du modèle thermique ( ) , et de la durée d’effacement souhaitée ( ), chaque indicateur individuel (I indiv (t)) étant représentatif de l’énergie propre à être économisée par ledit dispositif thermique (25) pendant une durée inférieure ou égale à la durée d’effacement souhaitée ( ), et un module d’estimation (80) configuré pour estimer, pour une pluralité de créneaux (t h ), l’indicateur collectif (I col (t h )) de réduction de consommation énergétique à partir des indicateurs individuels (I indiv (t)). Dispositif (20) selon la revendication 15, dans lequel le module d’acquisition (55) est configuré pour acquérir en outre un écart maximal ( ) souhaité à une température de consigne (T c (t)), le module de détermination (60) étant configuré pour déterminer l’indicateur individuel (I indiv (t)) de potentiel de réduction de consommation en outre à partir de l’écart maximal souhaité ( ).