FR 3085767 A1 20200313 FR 1858027 A 20180906 Procédé et système de pilotage et de régulation d'un système de chauffage d'un local [1] L'invention concerne le domaine de la domotique. Plus particulièrement, l'invention concerne le domaine de la régulation et du pilotage d'un système de chauffage dans un local, à partir d'une consigne de température. [2] Actuellement, un grand nombre de fabricants proposent des systèmes de pilotage du chauffage. La majorité d'entre eux offrent à l'usager la possibilité de piloter le chauffage en to fonction du confort thermique souhaité ou d'une consommation maximale souhaitée. Notamment, l'utilisateur peut définir trois températures de consigne à respecter, comprenant une température dite de confort pendant les phases de présence, une température dite de réduit, pendant les phases d'absence et une température dite de nuit, pendant les phases de présence en sommeil. 15 [3] On connaît notamment le brevet EP3045999 qui décrit un tel système. Plus précisément, ce document décrit un système de régulation et de pilotage d'un système de chauffage à partir de la connaissance d'une donnée de consommation de chauffage. Cette donnée est notamment obtenue par un dispositif de comptage d'énergie mis en place sur le 20 système de chauffage, difficile à mettre en place sur des systèmes de chauffage déjà existants. En variante, la donnée de consommation de chauffage est obtenue par une méthode de désagrégation de courbe de charge pour extraire la part chauffage de la consommation totale du local. Toutefois, ce type de méthode est source d'erreurs. 25 [004] La présente invention vient améliorer la situation. [005] A cet effet, elle propose un procédé, mis en oeuvre par des moyens informatiques, de pilotage et de régulation d'un système de chauffage d'un local, sur une période de temps, à partir d'une consigne reçue de température souhaitée définissant une température de consigne, 30 le procédé comportant les étapes : - élaborer un modèle de prévision d'une consommation d'énergie globale du local permettant de maintenir une température intérieure du local égale à la température de consigne sur ladite période de temps, ledit modèle s'exprimant comme fonction de la température de consigne par une équation - Gglobat® = GV (P ote'3,1 - Scui Talon 35 2 où : o Ce job représente la prévision de consommation d'énergie globale du local sur la période de temps, 5 o GV représente un coefficient de déperditions thermiques du local, o YQrateaH eprésente un écart positif entre la température intérieure estimée du local, et une température extérieure, ladite température intérieure estimée dépendant de la température de consigne, et ledit écart positif étant calculé à chaque pas de temps sur la période de temps, o Serti représente un seuil de température intérieure au-dessus de laquelle un système de chauffage du local est mis en état de désactivation, et oTao' représente une consommation d'énergie totale en usages non thermosensibles du local, obtenue à partir de données statistiques ; - puis, définir une valeur maximale de la consommation d'énergie globale sur la période de temps, et, à partir du modèle de prévision de la consommation d'énergie globale du local : o obtenir en sortie une nouvelle température de consigne à atteindre dans le local de sorte que la valeur de la consommation d'énergie globale du local est inférieure ou égale à ladite valeur maximale à l'issue de la période de temps. [6] Grâce à ces dispositions, le procédé selon l'invention permet d'effectuer une corrélation entre une température de consigne souhaitée et une consommation maximale souhaitée. De plus, le procédé ne nécessite pas l'utilisation de dispositif intrusif pour déterminer la consommation en chauffage d'un local. [7] Selon une réalisation, la construction du modèle comprend une étape préalable de : - une estimation de la température intérieure du local, l'estimation étant effectuée itérativement à chaque pas de temps At, ladite estimation étant exprimée par : - K2 ((Tin c( où : (i) est la température intérieure estimée à l'itération i, 3 o sont des paramètres thermiques du local, représentant respectivement une remontée en température et une chute en température du local, .K1, K, étant obtenus à partir de données historiques du local, 5 o Étatchaufftzge est un booléen représentant un état d'activation ou de désactivation du système de chauffage, o mars représente la température extérieure du bâtiment obtenue à partir de données météorologiques prévisionnelles ; avec une première itération où une valeur de la température initiale 7l est choisie parmi l'une et/ou l'autre de : o une température par défaut, o la température de consigne, o une température mesurée ; et 15 - une détermination d'un écart positif entre la température intérieure estimée du local et la température extérieure à chaque pas de temps est donnée par : où : représente l'écart positif entre la température intérieure estimée du 20 local et la température extérieure sur toute la période de temps ; o t(i) et T,TOD représentent respectivement ladite température intérieure estimée et ladite température extérieure à chaque pas de temps. [8] Ainsi, le modèle de prévision de la consommation d'énergie globale du local prend en 25 compte la consommation en chauffage mais aussi d'autres facteurs influant sur la température intérieure tels que la température extérieure. [9] Selon une réalisation, à chaque pas de temps, une comparaison entre la valeur de la température intérieure estimée et la valeur de la température de consigne est effectuée et, si la 30 valeur de la température estimée à ce pas de temps est supérieure ou égale à la valeur de la température de consigne : 4 - le booléen État., , est forcé à zéro. [10] Ainsi, lorsque la température intérieure atteint la température de consigne, le système de 5 chauffage est éteint. Le modèle de consommation tient compte de l'état éteint ou actif du système de chauffage de manière à obtenir précisément la consommation en énergie du local. [11] Selon une réalisation, à chaque pas de temps, une comparaison entre la valeur de la température intérieure estimée et la valeur de la température de consigne est effectuée et, si la 10 valeur de la température estimée à ce pas de temps est inférieure à la valeur de la température de consigne : - le booléen Êtat,,,ufic est forcé à 1. [0 1 2] De la même manière, lorsque la température intérieure est inférieure à la température de 15 consigne, le système de chauffe est à activer. Le modèle prend en compte l'état du système de chauffage afin de déterminer précisément la consommation en énergie du local. [013] Selon une réalisation, le procédé comprend une étape d'initialisation des valeurs du coefficient de déperditions thermiques GV, du seuil de température intérieure Seuil et de la 20 consommation d'énergie totale en usages non thermosensibles Talon au moyen d'une simulation thermique du local, permettant de générer des données de consommation en énergie estimées, lesdites valeurs initialisées des paramètres 6V, Seuil et Talon étant injectées dans ledit modèle de prévision. 25 [014] Ainsi, le procédé est opérationnel immédiatement et ne nécessite aucune phase d'initialisation lors de laquelle la consommation en énergie du local ne peut être estimée. Le procédé proposé présente donc un grand intérêt pour l'utilisateur, du fait de son efficacité immédiate. 30 [015] Selon une réalisation, au cours de la période de temps, les valeurs du coefficient de déperditions thermiques GV, du seuil de température intérieure Seuli et de la consommation d'énergie totale en usages non thermosensibles Talon sont recalculés, le recalcul comprenant : 5 - obtenir des données historiques de la température intérieure mesurée, de la température extérieure et de la consommation globale en énergie de chauffage, - calculer l'écart positif Potawientre la température intérieure mesurée du local 5 et la température extérieure, à chaque pas de temps pour chaque jour passé de la période de temps, et - si un nombre de jours passés, pour lesquels l'écart positif Potein est supérieur à un écart maximum et pour lesquels la consommation globale en énergie est supérieure à 0, est supérieur à un nombre de jours maximum, le 10 coefficient de déperditions thermiques 6V, le seuil de température intérieure Seuil et la consommation d'énergie totale en usages non thermosensibles Talon sont calculés selon un premier mode, ou - si un nombre de jours passés, pour lesquels l'écart positif Po !H est inférieur à un écart minimum, est supérieur à un nombre de jours minimum, le 15 coefficient de déperditions thermiques GY, le seuil de température intérieure Seuil et la consommation d'énergie totale en usages non thermosensibles Talon sont calculés selon un deuxième mode, - autrement, le coefficient de déperditions thermiques Gin, le seuil de température intérieure Sevî et la consommation d'énergie totale en usages non 20 thermosensibles Tai' conservent leurs valeurs respectives déterminées lors de l'étape d'initialisation. [016] Ainsi, les valeurs du coefficient de déperditions thermiques du seuil de température intérieure Seuil et de la consommation d'énergie totale en usages non thermosensibles Tai. 25 sont déterminées précisément en fonction du local et des données historiques en cours de fonctionnement. Cela permet de ne pas nécessiter de phase d'apprentissage préalable lors de laquelle le procédé n'est pas opérationnel. D'autre part, cela permet de calibrer les valeurs de ces coefficients de manière précise de sorte que la consommation en énergie peut être déterminée spécifiquement au local et à l'utilisateur du local. 30 6 [017] Selon une réalisation, le procédé comprend : - définir plusieurs températures de consigne respectivement associées à une pièce du local parmi une pluralité de pièces ; 5 - appliquer le modèle de prévision d'une consommation d'énergie globale aux pièces du local, de manière à ce que la température intérieure de chaque pièce atteigne la température de consigne associée à ladite pièce ; - puis, définir une valeur maximale de la consommation d'énergie globale sur la période de temps, et, à partir du modèle de prévision de la consommation 10 d'énergie globale du local : - obtenir en sortie de nouvelles températures de consigne à atteindre dans les pièces du local respectivement associées audites températures de consignes, de sorte que la valeur de la consommation d'énergie globale du local est inférieure ou égale à ladite valeur maximale à l'issue de la période de temps. 15 [18] Ainsi, le procédé permet également d'effectuer un pilotage et une régulation en chauffage pièce par pièce. Cela a l'avantage d'augmenter encore la précision du pilotage et de la régulation puisque la consommation en énergie de chaque pièce est précisément calculée. D'autre part, cela permet également de réduire la consommation en énergie du local puisque 20 toutes les pièces peuvent être à des températures différentes. Notamment, des températures de consigne peu élevées peuvent être associées à des pièces non utilisées ou moins utilisées. [19] Selon une réalisation, les valeurs du seuil de température intérieure Seuil et de la consommation d'énergie totale en usages non thermosensibles Tala conservent leurs 25 valeurs respectives déterminées pour le local, le procédé comprenant en outre, au cours de la période de temps, un calcul d'un coefficient de déperditions thermiques Glipt.èc' pour chacune des pièces du local , ledit calcul comprenant : - obtenir des données historiques de la température extérieure ; et, pour chaque pièce du local : 30 - obtenir des données historiques de la température intérieure mesurée et de la consommation en énergie de chauffage de ladite pièce du local ; - calculer un écart positif PotaiPëëcs entre la température intérieure mesurée de ladite pièce et la température extérieure à chaque pas de temps pour chaque jour passé de ladite période de temps, et 7 - si un nombre de jours passés, pour lesquels l'écart positif Fotdoeipiècg est supérieur à un écart maximum et pour lesquels la consommation globale en énergie de chauffage est supérieure à 0, est supérieur à un nombre de jours 5 maximum, le coefficient de déperditions thermiques de ladite pièce C.lisa est calculé par régression linéaire de l'équation : Ccel"Leia119(pièce) CV (PotakebipÉce - - autrement, le coeffic ent de déperditions thermiques GlAà', est calculé au moyen de l'équation : to GVrni,g = Ratiopiic, où Ratios est choisi parmi l'un et/ou l'autre - le ratio d'un nombre de radiateurs dans ladite pièce du local par le nombre de radiateurs total du local ; - le ratio d'une puissance de chauffage dans ladite pièce par une puissance de 15 chauffage totale du logement ; - le ratio d'une surface de ladite pièce par une surface totale du logement. [20] Ainsi, le recalcul du coefficient de déperdition thermique permet d'ajuster le modèle à 20 chacune des pièces à piloter. Cela permet d'augmenter la précision et l'efficacité du procédé selon l'invention puisque la construction du modèle est faite pièce par pièce. [21] Selon une réalisation, au moins deux températures de consigne sont définies, et de préférence au moins trois températures de consigne, chaque température de consigne étant 25 associée à une plage horaire récurrente d'un jour de la période de temps, les étapes de régulation du procédé étant effectuée pour chacune des températures de consigne, en fonction de la plage horaire dans laquelle les étapes du procédé sont mises en oeuvre. [22] Ainsi, le procédé est particulièrement bien adapté aux cas où plusieurs températures de 30 consigne sont reçues, que ce soit pour un pilotage du local dans sa globalité ou pièce par pièce. De plus, cela permet à un utilisateur de diminuer la valeur de la température de consigne pour une plage horaire lors de laquelle il est absent du local, par exemple. Cela entraine donc une réduction de la consommation en énergie du local. 8 [23] La présente invention vise aussi un programme informatique caractérisé en ce qu'il comporte des instructions pour la mise en oeuvre du procédé selon l'invention, lorsque ce 5 programme est exécuté par un processeur. [24] La présente invention vise aussi un système informatique de pilotage et de régulation d'un système de chauffage d'un local, sur une période de temps, à partir d'une consigne reçue de température souhaitée définissant une température de consigne, comportant un circuit de 10 traitement informatique pour la mise en oeuvre du procédé selon l'invention. [25] D'autres avantages et caractéristiques de l'invention apparaitront à la lecture de la description détaillée ci-après d'exemples de réalisations de l'invention, et à l'examen des dessins annexés sur lesquels 15 [26] La figure 1 représente schématiquement un système de pilotage et de régulation d'un système de chauffage d'un local selon un exemple de réalisation, [27] La figure 2 est un ordinogramme des principales étapes du procédé de pilotage et de 20 régulation d'un système de chauffage d'un local selon un exemple de réalisation, [28] La figure 3 est un ordinogramme des principales étapes de détermination du cas à utiliser pour calculer le coefficient de déperditions thermiques GV, le seuil de température intérieure Seüff et la consommation d'énergie totale en usages non thermosensibles Talon , 25 [29] La figure 4 illustre le calcul du coefficient de déperditions thermiques GV, du seuil de température intérieure et du consommation d'énergie totale en usages non thermosensibles Tait selon le cas d'initialisation, 30 [030] La figure 5 illustre le calcul du coefficient de déperditions thermiques G1â, du seuil de température intérieure Seuil et du consommation d'énergie totale en usages non thermosensibles Talaiet selon un mode, 9 [031] La figure 6 illustre le calcul du coefficient de déperditions thermiques 6V, du seuil de température intérieure Se-ull et du consommation d'énergie totale en usages non thermosensibles Talon selon un autre mode, 5 [32] La figure 7 est un ordinogramme illustrant les principales étapes de l'estimation de la température intérieure du local, [33] La figure 8 est un ordinogramme des principales étapes du procédé de pilotage et de 10 régulation d'un système de chauffage pour chaque pièce d'un local. [34] La figure 1 représente schématiquement un système SYST de pilotage et de régulation d'un système de chauffage d'un local, ou système SYST de pilotage et de régulation ou système SYST selon un exemple de réalisation. Par « local » on entend tout type d'habitation, 15 de bureaux, d'entreprises et plus généralement toute forme de bâtiment comprenant un système de chauffage. [35] Le système SYST de pilotage et de régulation comprend un dispositif DISP de pilotage et de régulation, ou dispositif DISP, utilisé dans un local. Le dispositif DISP comprend 20 notamment un processeur PROC sur lequel sont chargées des instructions permettant de mettre en oeuvre les étapes du procédé décrites ci-après. Le dispositif DISP comprend en outre une mémoire MEM. Le dispositif DISP peut également comprendre une interface utilisateur (non représentée). L'interface utilisateur peut notamment comprendre un écran, un clavier etc. 25 [036] En variante, le dispositif DISP est connecté à une interface utilisateur distincte, par tout moyen de connexion tel qu'une connexion sans-fil, internet, etc. L'interface utilisateur peut alors être un ordinateur, une télécommande ou encore un téléphone intelligent ou un boîtier dédié. 30 [037] Le local comprend au moins un ensemble de chauffage EC. L'ensemble de chauffage EC comprend tous les dispositifs consommant de l'énergie, qui peuvent être commandés et qui servent à chauffer le local. Notamment, l'ensemble de chauffage EC comprend au moins un système de chauffage SC, par exemple un radiateur, une chaudière à gaz, à fuel ou à bois. 10 [38] Selon une réalisation particulière, l'ensemble de chauffage EC peut être relié à un dispositif de mesure de la consommation d'énergie, par exemple un compteur COMPT. Le compteur COMPT peut être un compteur intelligent. 5 [39] L'ensemble de chauffage EC peut également être relié à un thermostat THERM permettant, d'une part, de mesurer la température intérieure du local et, d'autre part, d'activer ou de désactiver l'un des systèmes de chauffage SC de l'ensemble de chauffage EC. 10 [040] Le compteur COMPT peut être connecté au dispositif DISP. Notamment, le compteur COMPT transmet les données de consommation en énergie du local au dispositif DISP. La connexion peut être effectuée par une connexion sans-fil, par radio, par internet ou par une connexion filaire. 15 [041] Le thermostat THERM peut être connecté au dispositif DISP. Le thermostat THERM transmet les données de température intérieure du local, et éventuellement les états d'activation ou de désactivation des systèmes de chauffage, au dispositif DISP. Le dispositif DISP est apte à communiquer en retour avec le thermostat THERM de manière à piloter les systèmes de chauffage SC. La connexion peut être effectuée par une connexion sans-fil, par 20 radio, par internet ou par une connexion filaire. [42] Le dispositif DISP peut en outre être relié à un serveur SERV, par exemple par une connexion sans-fil de type connexion internet. Notamment, le serveur SERV peut servir à stocker les historiques de consommation en énergie du local, les températures intérieures 25 mesurées et estimées, des prévisions météorologiques comprenant la température extérieure, etc. [43] La figure 2 illustre les principales étapes du pilotage et de la régulation de systèmes de chauffage d'un local dans sa globalité sur une période de temps. 30 [44] A l'étape SI, une température de consigne Teo', est reçue. La température de consigne Too's peut être entrée par un utilisateur, via l'interface utilisateur par exemple. Cette température de consigne est à atteindre dans le local. [045] Selon une variante de réalisation non illustrée, l'utilisateur renseigne un planning horaire sur la période. Plus précisément, l'utilisateur renseigne au moins deux plages horaires 11 et de préférence au moins trois plages horaires. Les plages horaires peuvent être définies de manière journalière, de sorte que chaque jour de la période donnée est découpé selon ces plages horaires. En variante, les plages horaires peuvent être définies hebdomadairement, de 5 sorte que chaque semaine de la période donnée est découpée selon ces plages horaires, ou encore mensuellement. A titre indicatif, la période de temps donnée peut s'étendre de quelques semaines, à quelques mois voire quelques années. [46] Les plages horaires peuvent correspondre à une plage de présence, lors de laquelle 10 l'utilisateur est présent dans le local, une plage d'absence, lors de laquelle l'utilisateur est absent du local ou une plage de réduit, par exemple lorsque l'utilisateur dort. [47] L'utilisateur peut associer une température de consigne à chacune de ces plages horaires. Les températures de consigne peuvent différer selon les plages horaires. Le système 15 SYST est alors à même de construire un planning dans lequel chaque plage horaire est associée à une température de consigne à respecter. [48] La figure 2 présente le cas où une plage horaire n'a pas été définie, la régulation et le pilotage se font sur la base d'une unique température de consigne. Toutefois, il est entendu 20 que les étapes décrites ci-après s'appliquent de la même façon dans le cas où un planning a été construit. [49] Les étapes S2 à S4 illustrent l'estimation de la température intérieure sur la période donnée. Pour calculer la prévision de température intérieure, le procédé utilise une équation 25 d'évolution de la température intérieure en fonction des caractéristiques du bâtiment, de la température de consigne demandée et de données météorologiques (notamment la température extérieure) pouvant s'écrire comme suit Tmtti) = rmrif - 1) -- (Ki -État,h,,,if où : 30 o Tînt est la température intérieure estimée à l'itération t, o Vit, K, sont des paramètres thermiques du local, représentant respectivement une remontée en température et une chute en température du local, Ki, K; étant obtenus à partir de données historiques du local, 12 o Étatchnft est un booléen représentant un état d'activation ou de désactivation du système de chauffage, représente la température extérieure du bâtiment obtenue à partir de données 5 météorologiques prévisionnelle. [50] L'équation ci-dessus indique que la température intérieure évolue au cours du temps, en fonction de sa valeur au pas de temps précédent, de l'état du système de chauffage (en marche ou arrêté) et de l'écart entre températures intérieure et extérieure. [51] L'estimation de la température intérieure Tint comprend une première itération où ) = Trnriari )-Text(1 )))31 est égale à la température de consigne. Ainsi, l'estimation de la température intérieure Tint est directement corrélée à la température de consigne Lens reçue à l'étape SI [052] La température intérieure Tm( est estimée à chaque pas de temps I. Selon une réalisation, le pas de temps peut être égal à 15 minutes. La prévision de la température intérieure Tin, dans le local est faite sur la période donnée. Ainsi la température intérieure Tint est estimée à chaque pas de temps jusqu'à la fin de la période donnée où 1=1. Les valeurs de la 20 température intérieure Tmt sont par exemple stockées sur le serveur SERV. [53] Le chauffage du local est principalement lié à l'écart entre la température intérieure Tint du local et la température extérieure Tett, obtenue si besoin à un pas de temps égal à celui de T,et. Plus précisément, la température extérieure peut être obtenue au moyen de données 25 météorologiques stockées sur le serveur SERV ou encore au moyen d'une sonde de température placée à l'extérieur du local et connectée au dispositif DISP. [54] A l'étape S5, l'écart positif Potd,gif entre la température intérieure Tint et la température extérieure Text est calculé sur l'ensemble de la période comme suit : 30 10 15 où : 13 PVtdogg représente l'écart positif entre la température intérieure estimée du local et la température extérieure sur toute la période de temps; Tr(zi (1) représentent respectivement ladite température intérieure estimée et 5 ladite température extérieure à chaque pas de temps. [055] Cela permet d'obtenir, à l'étape S6 la consommation globale en énergie du local sur la période de temps en fonction de la température de consigne Te., reçue à l'étape Si, définie par la formule : sk,aait = G V Potan Talon 10 où : représente la prévision de consommation d'énergie globale du local sur la période de temps, représente un coefficient de déperditions thermiques du local, 15 représente un seuil de température intérieure au-dessus de laquelle un système de chauffage du local est mis en état de désactivation, et oTalon représente une consommation d'énergie totale en usages non thermosensibles du local, obtenue à partir de données statistiques. 20 [056] La construction du modèle permettant de prévoir la consommation en énergie du local par rapport à une température de consigne permet, à l'inverse, d'obtenir une valeur de température de consigne Toens' en fonction d'une consommation globale maximale en énergie, souhaitée par l'utilisateur. 25 [057] Les étapes S2 à S6 peuvent être réalisées sur le serveur SERV. [058] Ainsi, à l'étape S7, l'utilisateur peut entrer une consommation globale en énergie maximale. 30 [059] A partir de cette consommation globale maximale et du modèle construit, le procédé retourne en sortie une ou plusieurs valeurs de températures de consignes T0,','. L'utilisateur peut choisir la valeur de la température de consigne Tcons' qu'il souhaite atteindre dans le local. 14 [60] En variante, la ou les températures de consignes Teo',' présentées à l'utilisateur à l'étape S8 ont été filtrées afin de présenter à l'utilisateur les valeurs se rapprochant le plus de la 5 température de consigne Teo', initialement reçue à l'étape Si. [61] Dans le cas où un planning a été créé à partir de plusieurs valeurs de consigne Tc.., chacune associée à une plage horaire distincte, autant de valeurs de consigne Ternis' sont retournées à l'étape S8. Ainsi, si l'utilisateur a défini trois plages horaires associées respectivement à une température de consigne Teolls, un triplet de température Te'; est retourné à l'étape S8. L'utilisateur peut alors faire son choix entre les différents triplets de températures. [62] La détermination de la température de consigne Tcons' peut être effectuée en local. Par 15 exemple, le dispositif DISP détermine cette température de consigne. [63] Ainsi, la construction du modèle de prévision de la consommation globale en énergie est réalisée sur le serveur SERV et la détermination d'une température de consigne Lon; est réalisée en local sur le dispositif DISP. Cela permet de réaliser les opérations nécessitant 20 davantage de ressources en calculs à distance, de manière à ne pas surcharger le dispositif DISP, tandis que les opérations plus simples sont réalisées en local de manière à diminuer le temps de latence dû aux échanges de données entre le serveur SERV et le dispositif DISP. [64] La ou les températures de consigne Toms' retenues par l'utilisateur sont envoyées au 25 thermostat THERM. Le dispositif DISP se charge de piloter et réguler les systèmes de chauffage SC afin que la température intérieure Tint du local atteigne la température de consigne Tconi. [65] De manière plus précise, la réussite du modèle de prévision de la consommation globale 30 en énergie repose sur une connaissance précise et personnalisée, en fonction du local notamment, des paramètre thermiques K1, K2 (appelés paramètres dynamiques dans la suite de la description) ainsi que du coefficient de déperditions thermiques du seuil de température intérieure Seuil et de la consommation d'énergie totale en usages non thermosensibles Talon (appelés paramètres statiques dans la suite de la description). 15 [66] Les paramètres dynamiquesK1, K, permettent d'estimer, à partir d'une température de consigne Teons et d'une température extérieure Text, l'évolution dynamique de la température 5 intérieure Tint. [67] Ces paramètres sont nommés « remontée en température » et « chute en température » et peuvent être calculés par simple régression linéaire sur la base de données historique du local. Notamment, les données historiques du local peuvent correspondre aux estimations de 10 la température intérieure Tint, la température de consigne TCons et l'évolution de la température extérieure Te,, sur au moins une partie de la période de temps. Ces données historiques peuvent être stockées dans le serveur SERV. Ainsi, et pour assurer la cohérence du modèle, il apparaît que les équations utilisées pour estimer la température intérieure Tint et les paramètres dynamiques Kr, K sont les mêmes. 15 [68] Lorsqu'aucun historique n'est disponible, on peut fixer les paramètres à des valeurs par défaut, qui correspondent à un comportement d'un bâtiment moyen. Dans ce mode de réalisation, on prend par exemple comme valeurs par défaut K1=1.5 et K2=-0.05. 20 [069] Les paramètres statiques permettent d'estimer une consommation globale d'énergie nécessaire au respect d'une température de consigne connue. Comme dit plus haut, les paramètres statiques comprennent : le niveau de déperditions du local GV, la consommation en énergie des usages non thermosensibles Talon, s'exprimant en kWh, et la température intérieure (ou l'écart entre températures intérieure et extérieure) Seuil au-dessous de laquelle 25 les besoins de chauffage sont nuls ou négligeables. Plus la valeur de GV est élevée, plus il faudra de l'énergie pour maintenir le local à une température choisie. [070] Les paramètres statiques peuvent être calculés à partir des données d'historique du local, pouvant être stockées dans le serveur SERV. Plus précisément, le calcul des paramètres 30 statiques peut être effectué par régression linéaire de l'équation = GV (Potd,en - Sei ti + Talon et en utilisant les données d'historiques. 16 [071] Si les données d'historique ne sont pas suffisantes, des valeurs par défaut peuvent être utilisées. Ces valeurs sont par exemple calculées à partir d'un logiciel de simulation thermique et/ou d'un questionnaire renseigné par l'utilisateur sur le local. 5 [72] La figure 3 est un ordinogramme illustrant les principales étapes d'un procédé permettant de déterminer quelles valeurs des paramètres statiques sont à utiliser. [73] Aux étapes S9 à S11, les données historiques de la consommation globale en énergie, io de la température intérieure Tint mesurée du local et de la température extérieure Te,/ sont obtenues. Ces données historiques sont obtenues pour chaque jour passé d'au moins une partie de la période de temps. [74] A partir de ces données, l'écart positif Fo aentre la température intérieure Tint du 15 local et la température extérieure Text., à chaque pas de temps pour chaque jour passé de la période de temps est calculé à l'étape S12. [75] A l'étape S13, le nombre de jours où l'écart positif Faa gest strictement supérieur à un écart positif maximal, et où la consommation globale en énergie est supérieure à zéro est 20 dénombré. Si ce nombre de jours est supérieur à un nombre de jours maximal, prédéfini, les paramètres statiques sont calculés selon un premier mode, appelé « mode hiver » à l'étape S14. [76] Autrement, il est déterminé à l'étape S15 si un nombre de jours où l'écart positif 25 Pardsaest strictement inférieur à un écart minimum est supérieure à un nombre de jours minimum, prédéfini. Si c'est le cas, les paramètres statiques sont calculés selon un deuxième mode, appelé « mode été » à l'étape S17. Si ce n'est pas le cas, les paramètres statiques sont calculés selon un mode « initialisation ». A l'étape S18, les valeurs des paramètres statiques GV, Seuil et Talon sont retournés. 30 [77] Selon un exemple de réalisation, l'écart minimum est égale à 48°C.h, l'écart maximal est égal à 192°C.h, le nombre de jours minimum est égal à 5jours et le nombre de jours maximal est égal à 20jours. 17 [078] Les figures 4 à 6 sont des ordinogrammes illustrant les étapes permettant de calculer les paramètres statiques respectivement selon le mode « initialisation », « été » et « hiver ». 5 [079] Pour le mode « initialisation » des paramètres statiques, le procédé utilise par exemple un « profil utilisateur », qui permet de déterminer à gros traits les caractéristiques du local et de sa consommation en énergie. Ce profil peut être renseigné par l'utilisateur au moyen d'une série de questions (par exemple l'année de construction du logement, la surface, le nombre d'occupants, le type de système de chauffage, le type de système de production d'Eau Chaude 10 Sanitaire). Le profil peut être plus ou moins détaillé, pour améliorer la précision des estimations. [80] Selon une première variante, décrite à la figure 4, le profil client est injecté dans un logiciel de simulation thermique. [81] La valeur initialisée de Talon est calculée à partir des données historiques annuelles de la consommation globale en énergie et de la consommation en chauffage du local obtenues à l'étape S20. La différence entre ces deux consommations divisée par le nombre de jours sur l'année (365) permet d'obtenir la valeur de Talon à l'étape S21. [82] Aux étapes S22 et S23, les données historiques d'une consommation de chauffage hebdomadaire et d'une température extérieure moyenne pour N sous-périodes, chacune d'une durée d'une semaine, sont obtenues. 25 [083] Puis, pour chaque sous-période N d'une semaine, l'écart positif Potd,gaest calculé à l'étape S24, en prenant comme valeur de la température intérieure 20°C et un pas de temps égale à 24h, de sorte que l'écart positif Potd,,au est obtenu pour un jour représentatif de chaque sous-période N d'une semaine. 30 [084] Puis, la consommation en chauffage journalière moyenne est obtenue à partir de la consommation en chauffage hebdomadaire à l'étape S25. Aux étapes S26 et S27, seules les valeurs de la consommation en chauffage hebdomadaire supérieures à 0.001 et les valeurs de l'écart positif potg supérieures à 12 sont conservées de manière à ne pas fausser les résultats. 15 20 18 [85] A l'étape 528, les valeurs du niveau de déperditions thermiques GV et de Seuil sont obtenues par régression linéaires de l'équation C Mau! flaur = A Fttogni + où GV = A et Seuil = -B/A. [86] Les valeurs des paramètres statiques GV, Seuil et Talon du mode « initialisation » sont obtenues à l'étape S29. [87] Selon une deuxième variante de réalisation non représentée, les valeurs des paramètres statiques GV, Seuil et Talon du mode « initialisation » sont obtenues au moyen d'une base de données issue d'une simulation thermique ou d'une analyse statistiques des données de l'utilisateur. [88] La figure 5 illustre le calcul des paramètres statiques dans le mode « été ». [89] Aux étapes S30 et S31, les données historiques de la température intérieure T11,1 mesurée, de la température extérieure Text et de la consommation globale journalière sont 20 obtenues. [90] L'écart entre la température intérieure Tan et la température extérieure Tee est calculé à l'étape S32. Puis, à l'étape S33, le calcul de la médiane C des consommations globales journalières pour les jours pour lesquels l'écart entre la température intérieure Tint et la 25 température extérieure Text est inférieur à l'écart minimum est effectué. L'erreur quadratique moyenne RMSE est calculée pour ces mêmes jours selon la formule : Où N est le nombre de jours pour lesquels l'écart entre la température intérieure Tau et la 30 température extérieure Tee est inférieur à l'écart minimum. 5 10 15 [091] A l'étape S35, les valeurs de l'erreur RMSE et de la médiane C associées à l'écart minimum sont retournées. 19 [92] Les étapes S33 à S35 peuvent être répétées pour différentes valeurs de l'écart minimum. Par exemple, à chaque nouveau calcul, l'écart minimum est augmenté ou diminué d'un pas de 12°C.h, avec un écart minimum compris entre 48 et 192°C.h. De cette manière, une pluralité 5 de triplets de valeurs (ANISE, C, écart minimum) est obtenue à l'étape S35. [93] A l'étape S36, le paramètre statique Talon prend la valeur de la médiane C du triplet pour lequel l'erreur RMSE est la plus petite. Les paramètres statiques GV et Seuil prennent les valeurs des paramètres statiques initialisés, calculés à la figure 4. 10 [94] La figure 6 illustre calcul des paramètres statiques dans le mode « hiver » [95] A l'étape S40, l'écart positif Rota,atH journalier est calculé en fonction des données historiques des températures intérieures et extérieures. A l'étape S41, les données historiques 15 journalières de la consommation globale en énergie Caïc sont obtenues. [96] A l'étape S42, il est déterminé si le nombre de jours pour lesquels l'écart positif PotdEnest inférieur à l'écart maximal est supérieur à un nombre de jours maximum. Si c'est le cas, le coefficient C prend la valeur de la médiane de la consommation globale en énergie 20 pour les jours pour lesquels l'écart positif Porgegmest inférieur à l'écart maximal, comme illustré à l'étape S43. Autrement, le coefficient C prend la valeur du paramètre statique Talon initialisé (étape S44). [97] A l'étape S45 a lieu le calcul d'un coefficient A par régression linéaire de l'équation ("Ce:, - c) = A ' (Potp,,e) -Seuil) 25 Pour les jours où l'écart positif [Iotas/lest supérieur à l'écart minimum et où la consommation globale est supérieure à zéro. 30 [098] A l'étape 546 a lieu le calcul d'une consommation en énergie journalière estimée C,stim, pour chaque jour disponible des données historiques du local, telle que : 20 Cesti-7n = [099] A l'étape S47, l'erreur RIVISE est calculée pour chaque jour disponible telle que : RMSE - 5 [0100] Les valeurs du quadruplet (écart minimum, A, C et RNISE) sont retournées. [0101] Les étapes S42 à S47 décrites ci-avant peuvent être répétées pour différentes valeurs de l'écart minimum. Par exemple, à chaque itération de ces étapes la valeur de l'écart 10 minimum augmente ou diminue d'un pas de 12 °C.h, avec l'écart minimum compris entre 48°C.h et 192°C.h. Ainsi, une pluralité de quadruplets (écart minimum, A, C et RMSE) est retournée à l'étape S47. [0102] A l'étape S48, il est identifié si, parmi la pluralité de quadruplets retournés à l'étape 15 S47, il existe au moins un quadruplet donc une valeur du coefficient A est positive. Si c'est le cas, il est déterminé quelles sont les valeurs des coefficients A, C et de l'écart minimum pour lesquelles l'erreur RMSE est la plus petite. A l'étape S50, les paramètres statiques prennent respectivement les valeurs : GV=A 20 Talon = C, Seuil = écart minimum. [0103] Si ce n'est pas le cas, les paramètres statiques GV, Talon et Seuil conservent les valeurs du mode « initialisation » décrit en référence à la figure 4 (étape S49). [0104] Une fois calculés, les paramètres statiques GV, Talon et Seuil peuvent ne plus l'être. Par exemple, on peut prévoir de calculer à nouveau ces paramètres lors d'un changement d'utilisateur du local, ou encore lors d'un changement de mode de consommation d'énergie par exemple. 25 30 [0105] Le calcul de ces paramètres est particulièrement avantageux puisqu'il est fait à partir de données historiques. Les valeurs des paramètres statiques sont donc directement corrélées à 21 la consommation en énergie réelle de l'utilisateur. D'autre part, 1 initialisation de ces valeurs permet au procédé d'être directement opérationnel. 5 [0106] La figure 7 illustre de manière plus précise l'estimation de la température intérieure Tint à partir d'au moins une température de consigne. [0107] La prévision de température intérieure Tint permet de passer d'une température de consigne souhaitée à une température intérieure respectant cette température de consigne. to Dans le cas où l'utilisateur a défini un planning de contraintes dans lequel plusieurs plages horaires sont respectivement associées à une température de consigne, la prévision de la température intérieure Tint permet de respecter ces températures de consigne durant la plage horaire associée. 15 [0108] On décrit ci-après le cas où une unique température de consigne est reçue. Toutefois, les étapes décrites ci-après s'appliquent au cas où plusieurs températures de consigne sont reçues. [0109] Afin d'effectuer l'estimation de la température intérieure Tint, la température de 20 consigne Te', est reçue. Les prévisions de la température extérieure T,t sur la période de temps sont également reçues. Par exemple, ces prévisions sont reçues à partir du serveur SERV qui a accès aux prévisions météorologiques de la période donnée. Ces données sont par exemple reçues à chaque pas de temps, par exemple égal à 15 minutes. 25 [0110] Les valeurs des paramètres dynamiquesKi, K2 sont également reçues. Enfin, une température intérieure initiale peut être reçue. Cette valeur correspond par exemple à une température intérieure mesurée par le thermostat. En variante, il est possible de prendre une température intérieure par défaut, par exemple égale à 20°C. En variante, il est également possible de prendre une température intérieure initiale égale à la température de consigne 30 Tcons [0111] L'estimation comprend le calcul successif, à chaque pas de temps, d'une valeur Tint, qui représente l'estimation de la température intérieure, et une valeur Etatchaifflage, qui représente l'état du système de chauffage (1 si Marche ou 0 si Arrêt) à chaque pas de temps. 35 L'estimation de la valeur de Tint à chaque pas de temps i s'effectue à partir des données Tint et 22 Etate,''frage du pas de temps précédent (H), des données de la température de consigne Te'', et de la température extérieure Te', au pas de temps i et des valeurs des paramètres dynamiques K 5 [0112] Ainsi, l'étape S51 illustre l'initialisation des différents paramètres utiles à l'estimation de la température intérieure Tint. La température intérieure Tint au pas de temps i=0 est égale à la température intérieure initiale, et le booléen Etatch''ff'g, est forcé à 1 si la température intérieure initiale est inférieure à la température de consigne Toens et à 0 sinon. 10 [0113] A l'étape S52, le pas de temps est incrémenté de 1. L'étape S53 comprend alors l'estimation de la température intérieure Tait au pas de temps suivant i+1. Puis à l'étape 554, il est déterminé si la température intérieure au pas de temps i+1 est inférieure à la température de consigne. Si tel est le cas, le booléen Etatchauirage est forcé à 1, sinon il est forcé à 0. Tant 15 que le pas de temps est inférieur au pas de temps illustrant la fin de la période de temps, les étapes sont répétées en incrémentant systématiquement de 1 le pas de temps. A la fin de la période de temps (étape S56), le procédé est stoppé. La température intérieure estimée Tint à chaque pas de temps est retournée en sortie. Grâce à cette estimation de la température intérieure sur toute la période de temps, il est possible de construire le modèle de 20 consommation d'énergie tel que décrit en référence à la figure 2. [0114] Le procédé selon l'invention permet en outre de piloter et de réguler un système de chauffage de pièces d'un local, indépendamment les unes des autres. 25 [0115] Dans ce mode de réalisation, plusieurs températures de consigne l'00', respectivement associées à une pièce du local sont définies. Le modèle de prévision de la consommation d'énergie globale peut être appliqué à chacune des pièces du local, en fonction de la température de consigne Tc. associée à cette pièce. Une valeur maximale de la consommation d'énergie globale peut être définie, de sorte qu'à partir du modèle de prévision, 30 le procédé permet de retourner en sortie de nouvelles températures de consigne l'00', respectivement associées à une pièce du local. [0116] Afin de maintenir un degré de précision élevée, les paramètres statiques GV et Seuil peuvent être recalculés localement, c'est-à-dire pour chaque pièce. Le paramètre statique 23 Talon n'est pas à recalculer puisqu'il représente la consommation en énergie du local des usages non-thermosensibles. 5 [0117] Selon une réalisation, le paramètre Seuilpiè' de chaque pièce est égal au paramètre Seuil calculé pour le local. [0118] Selon une autre réalisation, le paramètre Seuilpiece de chaque pièce est obtenu par régression linéaire, tel que décrit ci-avant. Dans ce cas, le paramètre Talon est fixé à zéro. 10 [0119] Le calcul du paramètre GVp,,' est décrit ci-après pour une seule pièce. Il est entendu que ce calcul est fait pour chacune des pièces du local pouvant être pilotée au moyen du dispositif DISP. 15 [0120] Selon une première variante de réalisation, le paramètre Ratiopièce est calculé. Selon une réalisation, ce peut être le ratio du nombre de systèmes de chauffage (par exemple le nombre de radiateurs) présent dans la pièce par le nombre de systèmes de chauffage total dans le local. 20 [0121] En variante, ce peut être le ratio de la surface de la pièce par rapport à la surface totale du local. [0122] Selon une autre variante, ce peut être le ratio de la puissance du chauffage installé dans la pièce par rapport à la puissance du chauffage installé dans le local. 25 [0123] Ce paramètre est un indicateur sommaire des besoins de chauffage de la pièce, relativement à tout le local. [0124] Le paramètre statique GVpià' est calculé comme étant le produit du paramètre statique 30 GV du local par le paramètre Ratiopie'. [0125] Selon une deuxième variante de réalisation, le paramètre statique GVpièce de la pièce est calculé par régression linéaire sur les données historiques journalières selon l'équation : C; Mau f fPtieR = G V - SU, 24 Pour les jours où l'écart positif journalier Pu gHP pour la pièce est supérieur au paramètre statique Seuil du local. 5 [0126] Les deux variantes de réalisation sont l'une ou l'autre mises en oeuvre en fonction de certaines conditions. [0127] Ainsi, pour déterminer quelle variante de réalisation mettre en oeuvre, un dénombrement du nombre de jours pour lesquels l'écart positif Potace du logement est 10 strictement supérieur à l'écart maximum, par exemple égal à 192°C.h et pour lesquels la consommation globale du local est supérieure à 0. Si le nombre de jours est supérieurs à un nombre de jours maximal, par exemple égal à 20 jours, alors la deuxième variante de réalisation est utilisée. Autrement, la première variante de réalisation est utilisée. 15 [0128] La figure 8 est un ordinogramme illustrant les principales étapes du procédé de pilotage et de régulation pièce par pièce, dans le cas où le local comprend des pièces pilotables et des pièces non pilotables. [0129] Par « non-pilotables » on entend des pièces munies d'un système de chauffage ne 20 pouvant pas être piloté au moyen du dispositif DISP. Afin de résoudre ce problème, il est possible de définir deux sous-locaux fictifs. [0130] Un premier sous-local est défini comme comprenant une température intérieure identique pour toutes les pièces. La température intérieure peut être calculée comme la 25 moyenne des températures intérieures mesurées par le thermostat sur une sous-période de temps par exemple égale à une semaine. Ce premier sous-local est le local non-pilotable. [0131] Le deuxième sous-local est un local piloté pour lequel la température intérieure de chacune des pièces est directement liée à la température de consigne Tcoits souhaitée. 30 [0132] Ci-après on décrit un exemple de réalisation dans lequel un local comprend N pièces dont seulement P pièces sont pilotables, avec P égal ou inférieur à N. [0133] Pour chaque pièce pilotable p, il est possible de calculer la consommation de 35 chauffage pour la pièce p, sachant qu'elle est pilotée. Dans ce cas la température intérieure 25 peut varier en fonction de la température de consigne Te'' souhaitée dans cette pièce. Il est également possible de calculer la consommation de chauffage pour la pièce p, en supposant qu'elle n'est pas pilotée. Dans ce cas, la température intérieure est constante et identique pour 5 toutes les pièces. Enfin, il est possible de calculer la consommation de chauffage pour tout le local, en supposant qu'aucune pièce n'est pilotée. Cela est possible grâce aux coefficients statiques GVpréce et Seuilpièce de chaque pièce, ainsi qu'aux coefficients statiques du local GV et Seuil. fo [0134] Il est donc possible d'écrire : pilot , n piloté cimmtrapg (Lx ) ncinzu = non p lett ch t-rdaQ (pièce 'jPZ Et (Pièces non pilotées) 15 [0135] Soit et 20 [0136] Ainsi, en substituant le même terme dans les deux équations on obtient : 26 [0137] Afin de déterminer la consommation de chauffage du local la même approche que celle décrite en référence aux figures précédentes est utilisée. 5 [0138] Aux étapes S60, S61 et S62, les paramètres GV, Seuil, Talon, Potes GVp,e', Seuilp,,' et Fürd,efripïèr, sont reçus. L'écart positif Pot Ndu local est calculé quotidiennement, en prenant comme valeur de la température intérieure Tm/ une valeur unique, choisie parmi la moyenne des températures intérieures mesurées par un capteur dans 10 les pièces équipées d'un tel capteur sur une sous-période de temps, par exemple une semaine, ou une valeur par défaut, par exemple 20°C, ou encore la moyenne de la température intérieure mesurée dans une seule pièce sur une sous-période de temps, par exemple une semaine. 15 [0139] L'écart positif par pièce Pot agfips, est calculé pour chaque pièce, à partir de la température intérieure estimée, basée sur la température de consigne Toons Pour une pièce p donnée, l'écart positif Potannteest calculé en récupérant la température de consigne Teolls de cette pièce puis en estimant la température intérieure répondant à cette consigne. 20 [0140] A partir des valeurs des écarts positifs Patdgemnè et Pot ainsi que des paramètres statiques GV, GVpiece, Seuil, Seuilpie' et Talon, il est possible d'estimer les consommations de chauffage pour un jour j selon les formules suivantes C'unafl'fLoccti (tan ptioté, j) GV /rua (0;(Pettd, (r) - Seuil)) = GV - 25 = G Vota, CO- Seldienè,,, [0141] Ces formules correspondent respectivement aux étapes S64, S65 et S66. 30 27 [0142] Afin d'obtenir les consommations en chauffage rCdlagniagq L,,paGti, pour toutes les pièces, la somme de chacune de ces respectivement aux étapes S66 et S68. consommations en chauffage est effectuée L'étape S69, la consommation en chauffage global du logement est obtenue selon la formule Cth'ff's« (Local() REVENDICATIONS1. Procédé, mis en oeuvre par des moyens informatiques, de pilotage et de régulation d'un système de chauffage d'un local, sur une période de temps, à partir d'une consigne reçue de température souhaitée définissant une température de consigne, le procédé comportant les étapes élaborer un modèle de prévision d'une consommation d'énergie globale du local permettant de maintenir une température intérieure du local égale à la température to de consigne sur ladite période de temps, ledit modèle s'exprimant comme fonction de la température de consigne par une équation : Ce [ ter= GV m (P otsei - Seuil) -f- Talon où : représente la prévision de consommation d'énergie globale du local 15 sur la période de temps, o 6V représente un coefficient de déperditions thermiques du local, o Pot représente un écart positif entre la température intérieure estimée du local, et une température extérieure, ladite température intérieure estimée dépendant de la température de consigne, et ledit écart positif étant calculé à 20 chaque pas de temps sur la période de temps, o eu= d représente un seuil de température intérieure au-dessus de laquelle un système de chauffage du local est mis en état de désactivation, et o Taiem représente une consommation d'énergie totale en usages non thermosensibles du local, obtenue à partir de données statistiques ; 25 - puis, définir une valeur maximale de la consommation d'énergie globale sur la période de temps, et, à partir du modèle de prévision de la consommation d'énergie globale du local : o obtenir en sortie une nouvelle température de consigne à atteindre dans le local de sorte que la valeur de la consommation d'énergie globale du local 30 est inférieure ou égale à ladite valeur maximale à l'issue de la période de temps. 29 2. Procédé selon la revendication 1, dans lequel la construction du modèle comprend une étape préalable de : - une estimation de la température intérieure du local, l'estimation étant effectuée 5 itérativement à chaque pas de temps At, ladite estimation étant exprimée par : (i) = ', 1) -I- (K1 x Êtatcix¢r,;ja g K2 où : o T.u est la température intérieure estimée à l'itération i, o> Ki, K2. sont des paramètres thermiques du local, représentant respectivement une remontée en température et une chute en température du local, K1, étant obtenus à partir de données historiques du local, o Étatts, est un booléen représentant un état d'activation ou de désactivation du système de chauffage, o T,ar, représente la température extérieure du bâtiment obtenue à partir de données météorologiques prévisionnelles ; avec une première itération égale où une valeur de la température initiale TnitOÏ est choisie parmi l'une et/ou l'autre de : o une température par défaut, o la température de consigne, o une température mesurée ; et - une détermination d'un écart positif entre la température intérieure estimée du local et la température extérieure à chaque pas de temps est donnée par où : o Potd représente l'écart positif entre la température intérieure estimée du local et la température extérieure sur toute la période de temps; o (i) et T ) représentent respectivement ladite température intérieure estimée et ladite température extérieure à chaque pas de temps. 30 3. Procédé selon la revendication 2, dans lequel à chaque pas de temps, une comparaison entre la valeur de la température intérieure estimée et la valeur de la température de consigne est effectuée et, si la valeur de la température estimée à ce pas de temps est 5 supérieure ou égale à la valeur de la température de consigne - le booléen trat,hauffu est forcé à zéro. 4. Procédé selon l'une des revendications 2 à 3, dans lequel à chaque pas de temps, une comparaison entre la valeur de la température intérieure estimée et la valeur de la 10 température de consigne est effectuée et, si la valeur de la température estimée à ce pas de temps est inférieure à la valeur de la température de consigne : - le booléen Étaijc,a, est forcé à 1. 5 Procédé selon l'une des revendications 2 à 4, dans lequel le procédé comprend une 15 étape d'initialisation des valeurs du coefficient de déperditions thermiques (3V, du seuil de température intérieure Seau et de la consommation d'énergie totale en usages non thermosensibles Talor au moyen d'une simulation thermique du local permettant de générer des données de consommation en énergie estimées, lesdites valeurs initialisées des paramètres 6V, Seuil et 'Talon étant injectées dans ledit modèle de 20 prévision. 6. Procédé selon la revendication 4, caractérisé en ce que, au cours de la période de temps, les valeurs du coefficient de déperditions thermiques 6V, du seuil de température intérieure Seuil et de la consommation d'énergie totale en usages non thermosensibles 25 sont recalculés, le recalcul comprenant : - obtenir des données historiques de b température intérieure mesurée, de la température extérieure et de la consommation globale en énergie de chauffage, - calculer l'écart positif Pa entre la température intérieure mesurée du local et la température extérieure, à chaque pas de temps pour chaque jour passé de 30 la période de temps, et 31 - si un nombre de jours passés, pour lesquels l'écart positif P est supérieur à un écart maximum et pour lesquels la consommation globale en énergie est supérieure à 0, est supérieur à un nombre de jours maximum, le 5 coefficient de déperditions thermiques 6V, le seuil de température intérieure et la consommation d'énergie totale en usages non thermosensibles Tain., sont calculés selon un premier mode, ou - si un nombre de jours passés, pour lesquels l'écart positif FatA,4a est inférieur à un écart minimum, est supérieur à un nombre de jours minimum, le 10 coefficient de déperditions thermiques gir, le seuil de température intérieure Siet la consommation d'énergie totale en usages non thermosensibles TG ot sont calculés selon un deuxième mode, - autrement, le coefficient de déperditions thermiques GV, le seuil de température intérieure Sezel et la consommation d'énergie totale en usages non 15 thermosensibles T'aurai conservent leurs valeurs respectives déterminées lors de l'étape d'initialisation. 7 Procédé selon l'une des revendications 1 à 6, caractérisé en ce qu'il comprend : - définir plusieurs températures de consigne respectivement associées à une 20 pièce du local parmi une pluralité de pièces ; - appliquer le modèle de prévision d'une consommation d'énergie globale aux pièces du local, de manière à ce que la température intérieure de chaque pièce atteigne la température de consigne associée à ladite pièce ; - puis, définir une valeur maximale de la consommation d'énergie globale sur la 25 période de temps, et, à partir du modèle de prévision de la consommation d'énergie globale du local : - obtenir en sortie de nouvelles températures de consigne à atteindre dans les pièces du local respectivement associées audites températures de consignes, de sorte que la valeur de la consommation d'énergie globale du local est 30 inférieure ou égale à ladite valeur maximale à l'issue de la période de temps. 32 8. Procédé selon h revendication 6 en combinaison avec la revendication 7, dans lequel les valeurs du seuil de température intérieure Seuil et de la consommation d'énergie totale en usages non thermosensibles rctfori conservent leurs valeurs respectives déterminées 5 pour le local, le procédé comprenant en outre, au cours de la période de temps, un calcul d'un coefficient de déperditions thermiques Gir,:'. pour chacune des pièces du local , ledit calcul comprenant : - obtenir des données historiques de la température extérieure , et, pour chaque pièce du local : 10 - obtenir des données historiques de la température intérieure mesuréeet de la consommation en énergie de chauffage de ladite pièce du local ; - calculer un écart positif Patssunc, entre la température intérieure mesuréede ladite pièce et la température extérieure à chaque pas de temps pour chaque jour passé de ladite période de temps, et 15 - si un nombre de jours passés, pour lesquels l'écart positif Potat,gup est supérieur à un écart maximum et pour lesquels la consommation globale en énergie de chauffage est supérieure à 0, est supérieur à un nombre de jours maximum, le coefficient de déperditions thermiques de ladite pièce fi zE est calculé par régression linéaire de l'équation : Cce'ffa'gqpièce) = CV 4, (Pot erra tj,- Seuil} 20 - autrement, le coefficient de déperditions thermiques G ,ëce est calculé au moyen de l'équation : Ratiopitc- où Ratio - est choisi parmi l'un et/ou l'autre : mets 25 - le ratio d'un nombre de radiateurs dans ladite pièce du local par le nombre de radiateurs total du local ; - le ratio d'une puissance de chauffage installée dans ladite pièce par une puissance de chauffage totale du logement , - le ratio d'une surface de ladite pièce par une surface totale du logement. 30 33 9. Procédé selon l'une des revendications 1 à 8, dans lequel au moins deux températures de consigne sont définies, et de préférence au moins trois températures de consigne, chaque température de consigne étant associée à une plage horaire récurrente d'un jour 5 de la période de temps, les étapes de régulation du procédé étant effectuée pour chacune des températures de consigne, en fonction de la plage horaire dans laquelle les étapes du procédé sont mises en oeuvre. 10. Programme informatique caractérisé en ce qu'il comporte des instructions pour la mise fil) en oeuvre du procédé selon l'une des revendications précédentes, lorsque ce programme est exécuté par un processeur. 11. Système informatique de pilotage et de régulation d'un système de chauffage d'un local, sur une période de temps, à partir d'une consigne reçue de température souhaitée 15 définissant une température de consigne, comportant un circuit de traitement informatique (DISP) pour la mise en oeuvre du procédé selon l'une des revendications 1 à 9. 20