DESCRIPTION TITRE : 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 : o 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é, o 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, o 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, o 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 : o pour chaque créneau, sélection aléatoire d’un deuxième nombre de quantiles, o 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 : o un radiateur, o un climatiseur, o un ballon d’eau chaude, o un chauffe-eau, o un réfrigérateur, et o une pompe à chaleur, chaque dispositif thermique étant de préférence du même type, - le procédé comprend en outre une étape d’affichage de l’indicateur collectif sur un afficheur du dispositif électronique d’estimation, - lors de l’étape d’obtention le modèle thermique comprend en outre un coefficient de gain thermique obtenu pour chaque dispositif, lors de l’étape de détermination, une durée individuelle de retour à un équilibre thermique après un effacement est déterminée pour chaque dispositif et pour chaque instant, et - l’étape d’obtention du modèle thermique comprend, pour chaque dispositif thermique, les sous-étapes suivantes : o 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é, o pour chaque instant de chaque séquence de pilotage, détermination d’un coefficient instantané de gain thermique à partir de la ou des températures mesurées par le dispositif thermique audit instant, à partir de la puissance, et à partir d’une durée entre deux instants successifs, o pour chaque séquence de pilotage, calcul d’un coefficient de gain thermique de séquence à partir des coefficients instantanés de gain thermique de la séquence de pilotage, o calcul d’un coefficient de gain thermique de dispositif à partir des coefficients de gain thermique de séquence de chaque séquence de pilotage le modèle comprenant chaque coefficient de gain thermique de dispositif. 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 : - [Fig.1] la figure 1 est une vue schématique d’un système de contrôle énergétique comprenant un dispositif électronique d’estimation selon l’invention ; - [Fig.2] la figure 2 est une représentation d’une interface présentée par le dispositif électronique d’estimation de la figure 1 ; - [Fig. 3] la figure 3 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 figure 1 ; - [Fig. 4] la figure 3 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 figure 1 ; - [Fig.5] la figure 5 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 figure 3 ; - [Fig.6] la figure 6 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 figure 4 ; et - [Fig.7] la figure 7 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 figure 1, 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 Iindiv(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 Iindiv(t), un module 70 de conversion des indicateurs individuels de potentiel de réduction de consommation Iindiv(t) pour une journée dite extrême, un module 75 de calcul d’une pluralité de quantiles Qi, un module 80 d’estimation d’un indicateur collectif de réduction de consommation Icol, un module 85 d’affichage de l’indicateur collectif Icol, 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 figure 1, 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 Tint(t), Text(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. [ ^{MATH 1]} 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 Text 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 Tint 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 Tc 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 Tint 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 : [ ^{MATH 2]} ^ _{^ ^^ ^^ ^^( ^^) ≥ ^^ ^^( ^^) − ^^ ^^} 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 Text. 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 Text 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 Text 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 Text 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. En variante, un plus grand nombre de types de journées est considéré. Par exemple, les types de journée sont formés par découpage des moyennes de températures extérieures en intervalles de 2°C. Ainsi, à titre d’exemples les types de journée sont : [-8 ; - 6[, [-6 ; -4[, [-4 ; -2[, [-2 ; 0[, [0 ; 2[, [2 ; 4[, [4 ; 6[, [6 ; 8[, [8 ; 10[, [10 ; 12[, [12 ; 14[, [14 ; 16[, [16 ; 18[, et [18 ; 20[. 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 figure 1. 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 figure 2, 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 Icol doit être interprété. En particulier, le tutoriel indique que les sorties suivantes sont obtenues pour chaque créneau th, 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). Selon une variante non-représentée, l’utilisateur du dispositif d’estimation 20 interagit avec le dispositif d’estimation 20 pour renseigner en outre : la puissance moyenne des dispositifs thermiques 25. 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 Iindiv(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. [ ^{MATH 3]} 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. [ ^{MATH 4] où t est l’instant auquel est calculé l’évolution de la température intérieure Tint(t), et} _{^^ ^^ est une durée inférieure à la durée d’effacement souhaitée Δ ^^.} L ^{e 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 Tint(t+ ^^ ^^) s’écarte de la température de} consigne Tc(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 : [ ^{MATH 5]} ^ _{^ ^^ − ^^ ^^ ^^ ^^( ^^ + ^^ ^^) = ^^ ^^} 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 ti, ti+1 dans la durée d’effacement effective, l’énergie ^^ _{^^ ^^ ^^ ^^} ⁽ ^^ _^^ → ^^ _^^+1 ⁾ 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 th, 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 [ ^{MATH 6]} % _{^^ ^^ ^^( ^^ → ^^ + Δ ^^) = 100 ∗} ^{^^ ^^ ^^ ^^( ^^−> ^^ + Δ ^^) ∑Δ ^^} _^^ ^{^^ = ^^ ^^ ^^ ^^ ^^ ^^( ^^ ^^−> ^^ ^^+1)} ^ _{^ ^^ ^^ ^^ ^^ = 100 ∗ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ∗ Δ ^^} où ^^ _{^^ ^^ ^^ ^^} ( ^^ _^^ −> ^^ _^^+1 ) est égal à zéro si la température intérieure Tint(t) s’éloigne de la ^{température de consigne Tc(t) de plus que de l’écart maximal souhaité ^^ ^^, et ^^ ^^ ^^ ^^ ^^( ^^ ^^−>} _{^^ ^^+1) respecte l’équation 3, et ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ est la puissance installée du dispositif thermique 25.} Un créneau th est un créneau horaire délimitant une durée entre deux instants choisis. Tous les créneaux th 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 Iindiv(t), par exemple comme étant égal au pourcentage de puissance effaçable % _{^^ ^^ ^^} ( ^^ → ^^ + Δ ^^), respectant ainsi l’équation suivante. [ ^{MATH 7]} 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 Iindiv(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. [ ^{MATH 8]} 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 ^^ − 7 une semaine avant le jour ^^, Textrê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-temperatur e-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 : ^{[MATH 9]} % _{^^ ^^ ^^ ^^ ^^ ^^ ^^ê ^^ ^^( ^^ → ^^ + ^^ ^^) = 100} 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 th 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 th, une distribution de quantiles Qi(th). La distribution de quantiles Qi(th) forme alors, pour chaque créneau th, une distribution du pourcentage %eff(th) de la puissance installée effaçable pendant le créneau th. Le module d’estimation 80 est configuré pour estimer les percentiles des quantiles ordonnés pour chaque créneau Perc(%eff(th)). En complément facultatif, le module d’estimation 80 est en outre configuré pour calculer, pour chaque créneau th et pour chaque moyenne de quantile, un facteur de surdimensionnement ^^( ^^ _ℎ ) dont chaque percentile respecte l’équation suivante. [ ^{MATH 10]} où ^^ ^^ ^^ ^^ est la fonction percentile. De même, le module d’estimation 80 est par exemple configuré pour calculer, pour chaque créneau th, 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 % _{^^ ^^ ^^} (t _h ), calculés par le module d’estimation 80. Préférentiellement, l’indicateur collectif comprend en outre, pour chaque créneau th, 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 Icol(th) 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(th) 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 Tint(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 figure 7 illustrant un organigramme d’un procédé d’estimation de l’indicateur collectif Icol(th) 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 50 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 Tint(t), Text(t) ont été acquises lors de l’étape d’obtention 110, l’indicateur individuel Iindiv(t) propre à chaque dispositif. A cet effet, le module de détermination 60 détermine ces indicateurs individuels Iindiv(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 Iindiv(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 Text moyenne pendant la journée à laquelle l’indicateur individuel Iindiv(t) est déterminé. La température extérieure moyenne Text 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êmes 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 QI(th) d’indicateurs individuels Iindiv(t) à la pluralité de créneaux th comme expliqué précédemment. Lors d’une étape d’estimation 160, le module d’estimation 80 estime pour chaque créneau th, l’indicateur collectif Icol(th). 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 Iindiv(t) en fonction de la journée type qui leur est associée. Le module d’estimation 80 sélectionne, pour chaque créneau th, uniquement les quantiles Qi(th) respectifs des indicateurs individuels Iindiv(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 Qi(th). Lors de la troisième sous-étape de calcul 165, le module d’estimation 80 calcule la moyenne des quantiles Qi(th) 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 th les moyennes des quantiles Qi(th) calculés. La distribution de quantiles Qi(th) forme alors, pour chaque créneau th, la distribution de pourcentages %eff(th) 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 Icol(th) comprend le pourcentage de la puissance installée propre à être économisée %eff(th), 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 % _{^^ ^^ ^^ ^^} . Selon une variante, le fenêtre d’interface utilisateur présente en outre un champ permettant à un utilisateur de saisir la puissance moyenne des dispositifs thermiques 25 considérés. 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 Tint(t) est supérieure à la température extérieure Text(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 [ ^{MATH 13]} ^ _{^ ^^ ^^ ^^( ^^) ≤ ^^ ^^( ^^) + ^^ ^^} et l’équation 5 devient ^{[MATH 14]} ^ _{^ ^^ ^^ ^^( ^^ + ^^ ^^) − ^^ ^^ = ^^ ^^} 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. En complément facultatif, le modèle thermique comprend en outre, pour chaque dispositif thermique 25, un coefficient ^^ _{^^ ^^ ^^ ^^ ^^ ^^} de gain thermique associé au dispositif thermique 25 et qui sera décrit ci-après. Le module d’obtention 50 est en outre configuré pour déterminer, pour chaque instant t de chaque séquence de pilotage, un coefficient ^^( ^^) instantané de gain 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, de la puissance de chaque dispositif thermique ^^( ^^), et à partir du pas de temps D entre deux instants t successifs. Le module d’obtention 50 est par exemple configuré pour appliquer l’équation suivante. Le coefficient instantané de gain ^^( ^^) représente la variation de température Tint(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 pilotage. Le module d’obtention 50 est en outre configuré pour calculer, pour chaque séquence de pilotage, un coefficient ^^ _{^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^} de gain 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 pilotage correspondante. En particulier, le module d’obtention 50 est configuré pour calculer, pour chaque séquence de pilotage, le coefficient de séquence ^^ _{^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^} comme étant la médiane des coefficients instantanés ^^( ^^) durant ladite séquence de pilotage. Le module d’obtention 50 est en outre configuré pour calculer, pour chaque dispositif thermique 25, un coefficient ^^ _{^^ ^^ ^^ ^^ ^^ ^^} de gain thermique de dispositif 25. Par exemple, le module d’obtention 50 est configuré pour calculer le coefficient ^^ _{^^ ^^ ^^ ^^ ^^ ^^} du dispositif 25 comme étant la médiane des coefficients ^^ _{^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^} des séquences de pilotage. On comprend alors que le coefficient de gain 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), de la puissance P(t) et du coefficient de déperdition thermique ^^ _{^^ ^^ ^^ ^^ ^^ ^^} , lorsque le dispositif thermique 25 pilote ladite température intérieure Tint(t). En outre, selon ce complément facultatif, le module de détermination 60 détermine en outre une durée individuelle ^^ _{^^ ^^ ^^ ^^ ^^} de retour à l’équilibre suite à un effacement. A cet ^{effet, le module de détermination 60 est configuré pour déterminer la durée individuelle} _{^^ ^^ ^^ ^^ ^^ ^^ à partir du coefficient de déperdition thermique de dispositif ^^ ^^ ^^ ^^ ^^ ^^ ^^, du coefficient de} gain thermique de dispositif ^^ _{^^ ^^ ^^ ^^ ^^ ^^} , de la puissance du dispositif ^^ _{^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^} , et d’une température estimée à la fin de l’effacement ^^ _{^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^} ( ^^ _{^^ ^^ ^^ ^^ ^^ ^^} ), où ^^ _{^^ ^^ ^^ ^^ ^^ ^^} désigne l’instant final de l’effacement. On comprend que l’instant final de l’effacement ^^ _{^^ ^^ ^^ ^^ ^^ ^^} est : - soit égal à l’instant ^^ pour lequel la température intérieure ^^ _{^^ ^^ ^^} ( ^^ _{^^ ^^ ^^ ^^ ^^ ^^} ) a atteint l ^{a valeur ^^ ^ − ^^ ^^, la température estimée à la fin de l’effacement} _{alors égale à ^^ ^^ − ^^ ^^,} - soit égal à la somme entre l’instant ^^ et la durée d’effacement souhaité Δt, la température estimée à la fin de l’effacement ^^ _{^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^} ( ^^ _{^^ ^^ ^^ ^^ ^^ ^^} ) étant alors égale à ^^ _{^^ ^^ ^^} ( ^^) + ^^ _{^^ ^^ ^^ ^^ ^^ ^^} ∗ Δ ^^ ∗ ( ^^ _{^^ ^^ ^^} ( ^^) − ^^ _{^^ ^^ ^^} ( ^^)) d’après l’équation [MATH 4]. Le module de détermination 60 est par exemple configuré pour calculer, pour chaque dispositif 25 et pour chaque instant ^^, l’évolution de la température intérieure ^^ _{^^ ^^ ^^} , à une pluralité de pas ^^ successifs depuis l’instant final de l’effacement selon l’équation suivante. Le module de détermination 60 est configuré pour déterminer la durée individuelle de retour à l’équilibre ^^ _{^^ ^^ ^^ ^^ ^^} comme étant la durée ^^ ∗ ^^ au bout de laquelle la température ^{intérieure ^^ ^^ ^^ ^^ ( ^^ ^^ ^^ ^^ ^^ ^^ ^^ + ^^ ∗ ^^) devient proche de la température de consigne ^^ ^^ ( ^^ ^^ ^^ ^^ ^^ ^^ ^^ +} _{^^ ∗ ^^) à une précision ^^ ^^. C’est-à-dire la durée individuelle ^^ ^^ ^^ ^^ ^^ ^^ = ^^ ∗ ^^ pour laquelle la} ^{condition suivante est vérifiée :} La précision ^^ _^^ est par exemple égale à 0.1°C. On comprend donc que la durée individuelle ^^ _{^^ ^^ ^^ ^^ ^^} de retour à l’équilibre calculée pour une dispositif 25 et un ^^ correspond à la durée de retour à l’équilibre lorsqu’un effacement a été réalisé à l’instant ^^. Selon ce complément facultatif, le module d’association 65 est configuré pour associer, à chaque durée individuelle ^^ _{^^ ^^ ^^ ^^ ^^} un type de journée parmi les types de journée précités. Selon ce complément facultatif, le module de calcul 75 est configuré pour calculer, à partir de toutes les durées individuelles déterminées ^^ _{^^ ^^ ^^ ^^} , et pour chaque type de journée, des quantiles ^^ _{^^ ^^ ^^ ^^ ^^} ( ^^ _ℎ ) de durée individuelle ^^ _{^^ ^^ ^^ ^^ ^^} de manière analogue aux quantiles QI(th) d’indicateurs individuels Iindiv. Selon ce complément facultatif, le module d’estimation 80 est configuré pour ensuite effectuer, à partir des quantiles ^^ _{^^ ^^ ^^ ^^ ^^} ( ^^ _ℎ ) de durée individuelle ^^ _{^^ ^^ ^^ ^^ ^^} , les mêmes actions que celles définies précédemment à partir des quantiles QI(th) d’indicateurs individuels Iindiv, pour estimer des percentiles des quantiles ordonnés pour chaque créneau Perc( ^^(t _h )). Selon ce complément facultatif, de manière non représentée, le module d’affichage 85 est configuré pour afficher, par exemple sur l’afficheur 30, en outre, les percentiles Perc( ^^(th)) des quantiles ( ^^ _ℎ ) de durée de retour à l’équilibre. Selon ce complément facultatif et de manière optionnelle, 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é, uniquement pour les créneaux pour lesquels les percentiles Perc( ^^(t _h )) des quantiles ( ^^ _ℎ ) de durée de retour à l’équilibre sont inférieurs à un seuil prédéfini. Selon ce complément facultatif, en fonctionnement, le procédé est modifié comme suit. Lors de l’étape d’obtention 110, le modèle thermique obtenu par le module d’obtention comprend en outre le coefficient ^^ _{^^ ^^ ^^ ^^ ^^ ^^} de chaque dispositif 25 comme expliqué précédemment. Pour cela, préférentiellement, lors de la sous-étape 112, la puissance instantanée ^^( ^^) de chaque dispositif 25 est mesurée à chaque instant t. Lors de l’étape de détermination 130, le module de détermination 60 détermine en outre chaque durée individuelle de retour à l’équilibre ^^ _{^^ ^^ ^^ ^^ ^^} , comme expliqué précédemment. Lors de l’étape d’association 140, le module d’association 65 associe à chaque durée individuelle ^^ _{^^ ^^ ^^ ^^ ^^} , un type de journée. L ^{ors de l’étape de calcul 150, le module de calcul 75 calcule chaque quantile} _{durée de retour à l’équilibre.} Lors de l’étape d’estimation 160, le module d’estimation 80 estime les percentiles Perc( ^^(th)) des quantiles de durée de retour à l’équilibre similairement aux percentiles des quantiles d’indicateurs individuels. Lors de l’étape d’affichage 170, le module d’affichage 85 affiche en outre les percentiles Perc( ^^(t _h )) des quantiles de durée de retour à l’équilibre. Optionnellement, lors de l’étape de pilotage 180, le module de pilotage 90 envoie une instruction d’arrêt du pilotage de la température intérieure Tint(t) au(x) dispositif(s) thermique(s) 25 en fonction de l’effacement souhaité, uniquement pour les créneaux pour lesquelles les percentiles Perc( ^^(th)) des quantiles de durée de retour à l’équilibre sont inférieurs à un seuil prédéfini. L’évaluation du temps de retour à l’équilibre permet de mieux rendre compte des conséquences d’un effacement sur les dispositifs thermiques et leur environnement.