FR 3077119 A1 20190726 FR 1850513 A 20180123 PROCEDE ET DISPOSITIF DE DETERMINATION PREDICTIVE D'UN INDICATEUR DE CONFORT THERMIQUE D'UN LOCAL La présente invention a trait à un procédé et un dispositif de détermination prédictive d'un indicateur de confort thermique d'un local, pour des valeurs prédictives de conditions du local, en particulier pour des valeurs prédictives quelconques de la température de l'air à l'intérieur du local et de la température de l'air à l'extérieur du local. Il s'agit d'une évaluation prédictive du confort thermique du local car, avec la présente invention, l'indicateur de confort thermique peut être déterminé pour toute température de l'air à l'intérieur du local et toute température de l'air à l'extérieur du local, c'est-à-dire quelles que soient les conditions d'occupation du local et les conditions climatiques extérieures. L'invention vise en particulier à établir une cartographie prédictive du confort thermique dans un local, en déterminant un ou plusieurs indicateurs de confort thermique, à une ou plusieurs positions dans le local, pour des valeurs quelconques de la température de l'air à l'intérieur du local et de la température de l'air à l'extérieur du local. L'invention peut être appliquée pour prédire le confort thermique dans tout type de local, notamment dans un bâtiment ou dans un véhicule. Au sens de l'invention, on entend par "local" tout espace délimité par une enveloppe, plus particulièrement un espace d'habitat. Il peut s'agir d'un espace d'habitat fixe, tel qu'une maison individuelle ou un bâtiment, notamment à usage d'habitation ou tertiaire, ou une partie d'un tel bâtiment, par exemple un appartement dans un immeuble à plusieurs étages. Il peut également s'agir d'un espace d'habitat transportable, tel qu'une cabine dans un véhicule, en particulier dans une voiture, un camion, un train, un navire, un sous-marin, un avion, un vaisseau spatial. Dans un local, un occupant qui estime que son environnement thermique est trop chaud ou trop froid, même localement en un endroit précis du local, va améliorer son confort thermique en activant un dispositif de climatisation ou de chauffage. Cela tend à générer des surconsommations 2 d'énergie. Dans certains cas, il n'est pas possible de restaurer rapidement un environnement thermique confortable dans un local, ce qui peut entraîner en plus des surconsommations importantes d'énergie des risques pour la santé des occupants. Il est donc important d'évaluer de manière prédictive le 5 niveau de confort thermique d'un local pour identifier d'éventuels problèmes, en fonction des conditions d'utilisation du local et des conditions climatiques extérieures, et procéder si besoin à des travaux d'amélioration de l'isolation thermique du local. Pour évaluer le niveau de confort thermique d'un bâtiment, il existe 10 plusieurs indicateurs définis dans des normes internationales, parmi lesquels on peut citer à titre d'exemples non limitatifs : - la température opérative to (ou température de confort ressentie), qui est définie dans la norme ISO 7726; - la différence verticale de la température de l'air Jta,,, qui est définie 15 dans la norme ISO 7730; - l'asymétrie de température de rayonnement Jtpr, verticale ou horizontale, qui est définie dans la norme ISO 7730. Pour accéder de manière prédictive à ces indicateurs de confort thermique, il est nécessaire de connaître, là aussi de manière prédictive, 20 plusieurs grandeurs physiques du bâtiment, en particulier la température de l'air en différents points à l'intérieur du bâtiment et la température de surface des parois à l'intérieur du bâtiment (murs, vitrages, sols, plafonds, etc.). Les méthodes disponibles actuellement ne permettent pas de prédire de manière simple et rapide des températures d'air et de surfaces à 25 l'intérieur d'un bâtiment quelles que soient les conditions d'occupation du bâtiment et les conditions climatiques extérieures. Par exemple, une première méthode connue consiste à simuler le comportement thermique du bâtiment en utilisant un modèle numérique du bâtiment et plusieurs approximations. Cela permet de prédire les valeurs des 30 températures d'air et de surfaces à l'intérieur du bâtiment pour différentes conditions intérieures et extérieures au bâtiment, qui sont ensuite utilisées 3 pour calculer les indicateurs de confort thermique définis par les normes. Cette première méthode, qui est utilisée notamment par les bureaux d'étude thermique lors de la conception des bâtiments, nécessite du temps, des compétences en thermique et des ressources informatiques. De plus, cette 5 méthode repose sur un calcul dépendant des matériaux renseignés dans le modèle numérique. Or, dans le cas d'un logement neuf à réception, il est possible que la construction réellement achevée ne respecte pas les hypothèses prises durant la phase de conception. Dans le cas de logements anciens où il n'existe pas de connaissance des matériaux constituant 10 l'enveloppe, il est difficile de définir précisément un modèle du bâtiment. De plus, cette méthode ne permet pas de prédire la stratification de l'air et donc l'indicateur de différence verticale de la température de l'air Jta,v. Pour cela, il faut faire des simulations numériques des fluides ("Computational Fluid Dynamics", ou CFD), qui sont coûteuses en temps et ne sont pas précises. 15 Dans une deuxième méthode connue, les températures d'air et de surfaces à l'intérieur du bâtiment sont mesurées directement à l'aide de capteurs de température répartis dans le bâtiment. Les mesures sont réalisées sur un temps suffisamment long, en particulier supérieur à une année, et on estime alors le niveau de confort thermique futur du bâtiment en 20 supposant que les mêmes zones d'inconfort reviennent chaque année. La mise en oeuvre de cette deuxième méthode est longue, puisque les mesures doivent être réalisées sur au moins une année. De plus, cette méthode fournit des résultats approximatifs et ne permet pas d'évaluer le confort lors d'épisodes climatiques inhabituels tels que des canicules ou des épisodes de 25 grand froid. C'est à ces inconvénients qu'entend plus particulièrement remédier l'invention en proposant un procédé et un dispositif permettant de déterminer de manière prédictive un indicateur de confort thermique d'un local de manière simple et rapide, avec un coût modéré et une précision raisonnable, 30 quelles que soient les conditions d'occupation du local et les conditions climatiques à l'extérieur du local. 4 A cet effet, l'invention a pour objet un procédé de détermination prédictive d'un indicateur de confort thermique d'un local, pour des valeurs prédictives quelconques de la température de l'air à l'intérieur du local et de la température de l'air à l'extérieur du local, caractérisé en ce qu'il comprend 5 des étapes dans lesquelles : - on détermine expérimentalement au moins un paramètre d'un modèle thermique descriptif du local, par ajustement entre le modèle thermique et des mesures réalisées in situ dans le local, où les mesures comprennent des mesures de la température de l'air à l'intérieur du local et 10 des mesures soit d'au moins une température supplémentaire à l'intérieur du local, soit du flux thermique à travers des parois du local correspondant à l'indicateur de confort thermique ; - on calcule une ou des températures de confort prédictives intervenant pour la détermination de l'indicateur de confort thermique, pour 15 lesdites valeurs prédictives quelconques de la température de l'air à l'intérieur du local et de la température de l'air à l'extérieur du local. L'invention repose sur la détermination expérimentale d'au moins un paramètre d'un modèle thermique descriptif du local, par ajustement entre le modèle thermique et des mesures réalisées in situ dans le local, puis 20 l'utilisation de ce ou ces paramètres pour calculer des températures de confort prédictives, pour lesdites valeurs prédictives quelconques de la température de l'air à l'intérieur du local et de la température de l'air à l'extérieur du local, permettant d'accéder à la valeur prédictive de l'indicateur de confort thermique. 25 En pratique, des exemples de paramètres de modèles thermiques descriptifs du local susceptibles d'être déterminés dans le cadre de l'invention comprennent, de manière non limitative : - le coefficient de déperdition thermique HLC du local (en W.K1) ; - un coefficient de transmission effectif Heff (en W.K1), qui est 30 déterminé en un ou plusieurs points à l'intérieur du local choisis pour correspondre à l'indicateur de confort thermique ; par exemple, lorsque 5 l'indicateur de confort thermique est la température opérative to en un point à l'intérieur du local, on peut déterminer la valeur du coefficient de transmission effectif Heff en ce point à partir de mesures de la température moyenne radiante à l'intérieur du local en ce point à l'aide d'un capteur de température 5 globe noir ou, en variante, on peut déterminer les valeurs du coefficient de transmission effectif Heff des différents éléments de construction des parois du local vus en ce point à partir de mesures des températures de surface intérieure de ces éléments de construction ; selon un autre exemple, lorsque l'indicateur de confort thermique est la différence verticale de la température 10 de l'air Jta,,, on peut déterminer les valeurs du coefficient de transmission effectif Heff en deux points situés respectivement à 0,10 m et 1,70 m selon une verticale du local, à partir de mesures de la température de l'air à l'intérieur du local en ces deux points ; de manière générale, on peut définir un coefficient de transmission effectif Heff en chaque point de l'espace et/ou 15 pour chaque type de mesure de température (température de surface, température de l'air, température moyenne radiante) comme l'inverse de la résistance thermique entre ce point et le noeud de température extérieure, ce coefficient étant alors estimé par inversion de modèle ou identification des paramètres d'un modèle thermique descriptif du local et des mesures 20 expérimentales réalisées dans le local ; - le coefficient de transmission thermique U (en W.m-2.K1), qui est déterminé pour chaque élément de construction des parois du local correspondant à l'indicateur de confort thermique ; par exemple, lorsque l'indicateur de confort thermique est la température opérative to en un point à 25 l'intérieur du local, on peut déterminer le coefficient de transmission thermique U pour chaque élément de construction des parois du local vu en ce point, notamment à partir de mesures du flux thermique à travers ces éléments de construction ou de la température de surface de ces éléments de construction ; selon un autre exemple, lorsque l'indicateur de confort 30 thermique est l'asymétrie de température de rayonnement verticale Jtp,,,, on peut déterminer le coefficient de transmission thermique U pour chaque 6 élément de construction des parois horizontales du local en vis-vis, notamment à partir de mesures du flux thermique à travers ces éléments de construction ou de la température de surface de ces éléments de construction. 5 D'autres paramètres de modèles thermiques peuvent également être déterminés dans le cadre de l'invention, tels qu'un facteur solaire dans le cas où le modèle thermique du local prend en compte la contribution du rayonnement solaire au chauffage du local. Dans le cadre de l'invention, pour une paroi du local, on distingue 10 différents "éléments de construction" de la paroi lorsque la paroi est formée de plusieurs éléments ayant des propriétés d'isolation thermique significativement différentes. Ainsi, par exemple, pour une paroi du local formée par un mur muni d'un vitrage, la paroi comprend deux éléments de construction qui sont, respectivement, le mur et le vitrage. Les paramètres 15 tels que le coefficient de transmission effectif Heff ou le coefficient de transmission thermique U sont alors déterminés pour chaque élément de construction de la paroi. Selon un aspect de l'invention, le ou les paramètres du modèle thermique sont déterminés en utilisant les variations transitoires de la 20 température à l'intérieur du local lorsque le local est soumis à des sollicitations internes contrôlées et dans un environnement externe mesuré. L'analyse quantitative de la variation de températures à l'intérieur du local et/ou du flux thermique à travers les parois du local permet alors de déterminer quantitativement le ou les paramètres du modèle thermique sur 25 une période courte, s'étendant notamment sur quelques heures, en limitant l'influence des paramètres susceptibles de modifier le comportement dynamique du local. En particulier, la brièveté des mesures permet de s'affranchir de l'influence des conditions d'utilisation du local et des variations des conditions climatiques extérieures. 30 Dans un mode particulier de réalisation, le procédé selon l'invention comprend des étapes dans lesquelles : 7 - sur au moins deux périodes de temps Di, successives correspondant à des puissances Pic distinctes de chauffe appliquées à l'intérieur du local, on procède à une campagne de mesures de la température de l'air à l'intérieur du local et d'au moins une température à l'intérieur du local correspondant à 5 l'indicateur de confort thermique à intervalles de temps rapprochés, ainsi qu'à la détermination de la température de l'air à l'extérieur du local à intervalles de temps rapprochés ; - on détermine la valeur du coefficient de déperdition thermique HLC du local, par ajustement entre : 10 d'une part, un modèle thermique exprimant la variation temporelle de la température de l'air à l'intérieur du local en fonction de la puissance de chauffe appliquée dans le local et de la température de l'air à l'extérieur du local, et d'autre part, l'évolution mesurée de la température de l'air à l'intérieur 15 du local en fonction du temps ; - on détermine la valeur d'un coefficient de transmission effectif Heff en chaque point de mesure de la température à l'intérieur du local correspondant à l'indicateur de confort thermique, par ajustement entre : d'une part, un modèle thermique exprimant la variation temporelle de la 20 température à l'intérieur du local correspondant à l'indicateur de confort thermique en fonction de la puissance de chauffe appliquée dans le local et de la température de l'air à l'extérieur du local, et d'autre part, l'évolution mesurée de la température à l'intérieur du local correspondant à l'indicateur de confort thermique en fonction du temps ; 25 - on calcule la ou les températures de confort prédictives T'nfort pred intervenant pour la détermination de l'indicateur de confort thermique, à partir des valeurs préalablement déterminées du coefficient de déperdition thermique HLC du local et du coefficient de transmission effectif Hely, pour lesdites valeurs prédictives quelconques de la température de l'air à 30 l'intérieur du local T air in pred et de la température de l'air à l'extérieur du local Tair out pred selon la relation : 8 HLC Tconfort pred = Tair out pred u Tair in pred Tair out pred) iieff Dans ce mode de réalisation, à titre d'exemples : - lorsque l'indicateur de confort thermique est la température opérative to en un point à l'intérieur du local, on peut déterminer la valeur du coefficient de transmission effectif Heff en ce point à partir de mesures de l'évolution 5 temporelle de la température moyenne radiante à l'intérieur du local Trad in k conf en ce point et/ou on peut déterminer la valeur du coefficient de transmission effectif Heff pour chaque élément de construction des parois du local vu en ce point à partir de mesures de l'évolution temporelle de la température de surface à l'intérieur du local Tsurfin k conf de l'élément de 10 construction ; - lorsque l'indicateur de confort thermique est la différence verticale de la température de l'air Jta,, selon une verticale du local, on peut déterminer les valeurs du coefficient de transmission effectif Heff en deux points de hauteurs différentes selon cette verticale à l'intérieur du local à partir de 15 mesures de l'évolution temporelle de la température de l'air à l'intérieur du local Tair in k eeuf en ces deux points ; - lorsque l'indicateur de confort thermique est l'asymétrie de température de rayonnement Jtpr verticale, respectivement horizontale, on peut déterminer la valeur du coefficient de transmission effectif Heff de 20 chaque élément de construction des parois horizontales du local en vis-vis, respectivement des parois verticales du local en vis-vis, à partir de mesures de l'évolution temporelle de la température de surface à l'intérieur du local Tsurf in k conf Dans un autre mode de réalisation, le procédé selon l'invention 25 comprend des étapes dans lesquelles : - on détermine expérimentalement le coefficient de transmission thermique U de chaque paroi du local correspondant à l'indicateur de confort thermique, par ajustement entre un modèle thermique et des mesures réalisées in situ dans le local, où les mesures comprennent des mesures de de chacune de ces parois. 9 la température de l'air à l'intérieur du local et des mesures soit du flux thermique à travers la paroi, soit de la température de surface de la paroi à l'intérieur du local ; - on calcule une température de surface prédictive T'rf. p'd de 5 chaque paroi du local correspondant à l'indicateur de confort thermique, pour lesdites valeurs prédictives quelconques de la température de l'air à l'intérieur du local Tair in pred et de la température de l'air à l'extérieur du local Tair out pred, selon la relation : U Tsurf pred = Tair in predhin Tair in pred Tair out pred) avec hin le coefficient d'échange convectif de l'air à l'intérieur du local. 10 Dans ce mode de réalisation, à titre d'exemples : - lorsque l'indicateur de confort thermique est la température opérative to en un point à l'intérieur du local, on peut déterminer la valeur du coefficient de transmission thermique U de chaque élément de construction des parois du local vu en ce point, à partir de mesures du flux thermique à travers ces 15 éléments de construction ou de la température de surface de ces éléments de construction ; - lorsque l'indicateur de confort thermique est l'asymétrie de température de rayonnement Jtpr verticale, respectivement horizontale, on peut déterminer la valeur du coefficient de transmission thermique U de 20 chaque élément de construction des parois horizontales, respectivement verticales, du local en vis-vis, à partir de mesures du flux thermique à travers ces éléments de construction ou de la température de surface de ces éléments de construction. Selon une première variante, pour déterminer expérimentalement le 25 coefficient de transmission thermique U de chaque paroi du local correspondant à l'indicateur de confort thermique, on procède aux étapes suivantes : - sur au moins deux périodes de temps Di, successives correspondant à des puissances Pk distinctes de chauffe appliquées dans le local, on 10 procède à une campagne de mesures de la température de l'air à l'intérieur du local Tai, in I, et du flux thermique à travers la paroi qk à intervalles de temps rapprochés, ainsi qu'à la détermination de la température de l'air à l'extérieur du local Tai, 't k à intervalles de temps rapprochés ; 5 - on détermine la valeur du coefficient de transmission thermique U de la paroi par ajustement entre : d'une part, un modèle thermique exprimant la variation temporelle de la température de l'air à l'intérieur du local en fonction du flux thermique à travers la paroi et de la température de l'air à l'extérieur du local, et 10 d'autre part, l'évolution mesurée de la température de l'air à l'intérieur du local Tai, k ou du flux thermique à travers la paroi qk en fonction du temps. Selon une deuxième variante, pour déterminer expérimentalement le coefficient de transmission thermique U de chaque paroi du local 15 correspondant à l'indicateur de confort thermique, on procède aux étapes suivantes : - sur au moins deux périodes de temps Di, successives correspondant à des puissances Pk distinctes de chauffe appliquées dans le local, on procède à une campagne de mesures de la température de l'air à l'intérieur 20 du local Tai, k et de la température de surface de la paroi à l'intérieur du local T'rf. k à intervalles de temps rapprochés, ainsi qu'à la détermination de la température de l'air à l'extérieur du local Tai, 't k à intervalles de temps rapprochés ; - on détermine la valeur de la grandeur U Ihin, où hi, est le coefficient 25 d'échange convectif de l'air dans le local au voisinage de la paroi, par ajustement entre : d'une part, un modèle thermique exprimant la variation temporelle de la température de l'air à l'intérieur du local en fonction de la température de surface de la paroi à l'intérieur du local et de la température de l'air à 30 l'extérieur du local, et 11 d'autre part, l'évolution mesurée de la température de l'air à l'intérieur du local Tai, in I, ou de la température de surface de la paroi à l'intérieur du local T'rj in k en fonction du temps. On accède ainsi à la valeur de la grandeur U 1 hm, et on en déduit la 5 valeur du coefficient de transmission thermique U de la paroi à partir d'une valeur standard et répertoriée du coefficient d'échange convectif de l'air hi, pour le type de paroi considéré. A titre d'exemple, pour une paroi verticale de bâtiment, une valeur standard du coefficient d'échange convectif de l'air hi, est 8 ± 2 W.m-2.K1. 10 Dans le cadre de l'invention, on entend par "puissance de chauffe du local" toute condition opératoire générant une variation de la température à l'intérieur du local, pour des conditions de température données à l'extérieur du local. Il est entendu que la puissance de chauffe peut être positive, nulle ou négative. Une puissance de chauffe positive correspond à un apport de 15 chaleur dans le local, qui peut être obtenu à l'aide d'un appareil de chauffage. Une puissance de chauffe négative correspond à un apport de froid dans le local, qui peut être obtenu à l'aide d'un appareil de climatisation. Dans le cas d'une puissance de chauffe nulle ou sensiblement nulle appliquée dans le local, la variation de la température dans le local peut 20 résulter d'un déséquilibre entre les températures initiales de l'air à l'intérieur du local et à l'extérieur du local. Par souci de simplicité, dans le cadre de cette description, les termes "chauffe" ou "chauffage" peuvent désigner tout autant un apport de chaleur qu'un apport de froid. Les périodes de temps Di, peuvent être soit disjointes, soit 25 immédiatement successives les unes aux autres. Dans ce dernier cas, on peut considérer que le procédé est réalisé dans sa globalité sur une période de temps continue, formée par la succession des périodes de temps Di,. Selon une réalisation avantageuse, les mesures in situ du procédé selon l'invention sont mises en oeuvre sur deux périodes de temps D1 et D2 30 successives correspondant à deux consignes de puissances distinctes de chauffe P1 et P2 appliquées dans le local. 12 Selon un aspect de l'invention, la première puissance de chauffe P1 appliquée sur la première période de temps D1 et la deuxième puissance de chauffe P2 appliquée sur la deuxième période de temps D2 ont des valeurs éloignées l'une de l'autre. 5 Dans un exemple particulier de réalisation, la première puissance de chauffe P1 appliquée sur la première période de temps D1 est strictement positive ou strictement négative, de manière à atteindre une température de l'air à l'intérieur du local strictement supérieure ou strictement inférieure à la température de l'air à l'extérieur du local, puis la deuxième puissance de 10 chauffe P2 appliquée sur la deuxième période de temps D2 est nulle ou sensiblement nulle, de manière à obtenir une évolution libre du local résultant du déséquilibre entre la température de l'air à l'intérieur du local et la température de l'air à l'extérieur du local. De préférence, en vue de limiter le temps de mise en oeuvre du 15 procédé, tout en réduisant les contributions non contrôlées résultant par exemple de l'utilisation du local ou des conditions climatiques, les mesures in situ du procédé selon l'invention sont mises en oeuvre : - alors que le local est inoccupé ; et/ou - sur des périodes de temps Di, pour lesquelles le rayonnement 20 solaire est faible, de préférence nul. Dans le cadre de l'invention, le modèle thermique utilisé peut être de tout type connu de l'homme du métier. Dans le cadre de l'invention, on effectue un ajustement paramétrique entre un modèle thermique descriptif du local et des mesures réalisées in situ 25 dans le local. Les étapes du procédé selon l'invention, en particulier pour l'ajustement entre le modèle thermique et les mesures in situ, peuvent être mises en oeuvre à l'aide d'une unité de calcul appropriée. Il peut s'agir notamment d'un terminal informatique, par exemple un ordinateur portable ou non, une tablette, un téléphone intelligent ("smartphone" en anglais), 30 comportant un module d'acquisition pour acquérir les mesures requises par 13 le procédé et un module de traitement pour exécuter tout ou partie des étapes de calcul du procédé à partir des mesures acquises. Le modèle thermique utilisé peut être de tout type connu de l'homme du métier. 5 Selon un aspect de l'invention, le modèle thermique utilisé est un modèle d'identification paramétrique (ou modèle d'identification de système), tel que les modèles décrits dans l'ouvrage "Time Series Analysis" de Henrik Madsen, Chapman & Hall, 2008 (ISBN-13: 978-1420059670 0). En particulier, le modèle thermique utilisé peut être un modèle 10 autorégressif, notamment un modèle ARX (Auto Regressive model with eXternal inputs). En variante, d'autres modèles autorégressifs peuvent être utilisés comme par exemple des modèles ARMAX ou ARMA. De façon connue, un modèle ARX est un modèle autorégressif définissant une sortie y(t) en fonction d'une ou de plusieurs entrées u(t), v(t) 15 et d'un résidu aléatoire de modélisation caractérisé par un bruit blanc de moyenne nulle e(t), t désignant l'instant d'échantillonnage considéré. Le modèle ARX s'écrit : (1 + A(p-1))yk = B(p-1)uk + D(p-1)vk + ek avec : 1+ A(p-1) = 1+ aip-1 + + anp-n B(p1) = b1p-1 + + bnp-n D(p1) = d1p-1 + + dnp-n et : P-luk = uk-i P-lvk = vk-i 20 p-1 est appelé opérateur de retard et prend en compte les états passés influençant le système au temps présent. n est l'ordre du modèle ARX. Concernant le choix de l'ordre n du modèle ARX, la valeur de n doit être suffisamment élevée pour prendre en compte l'inertie du système, mais suffisamment faible pour éviter un sur-paramétrage du modèle. 14 Selon une caractéristique avantageuse, l'ordre n du modèle ARX est sélectionné de manière itérative, par incrémentation, en partant de n = 1 et jusqu'à ce qu'il n'y ait pas d'autocorrélation ou de corrélation croisée entre le résidu de modélisation e du modèle ARX et les évolutions des mesures in 5 situ en fonction du temps, déterminées sur les périodes de temps Dic, correspondant au modèle. Les étapes d'identification d'un modèle sont connues en soi et ne sont donc pas détaillées plus avant ici. Les coefficients à identifier par ajustement paramétrique (ou "fitting") du modèle ARX sur l'évolution des mesures in situ 10 sont les coefficients a, , b, di, j = 1, , n, auxquels est relié le paramètre recherché du modèle thermique. Selon un autre aspect de l'invention, le modèle thermique utilisé est un modèle R-C avec un nombre adapté de résistances et de capacités. De préférence, le modèle thermique est un modèle R-C simple avec une 15 résistance et une capacité. En variante, le modèle thermique peut être un modèle R-C plus complexe que le modèle R-C simple, de type mRnC avec m et n deux nombres entiers, tel qu'un modèle dit "2R2C" avec deux résistances et deux capacités, ou encore un modèle dit "3R2C" avec trois résistances et deux capacités. 20 En particulier, dans le cadre de l'invention, on peut utiliser un modèle R-C simple pour décrire un local, avec deux noeuds de température homogène, l'un à l'intérieur du local et l'autre à l'extérieur du local, qui sont séparés par une résistance représentant le coefficient de déperdition thermique global HLC du local et décrivant la perte par transmission et 25 infiltration à travers l'enveloppe du local. Le noeud de température à l'intérieur du local est relié à un condensateur qui représente la masse thermique ou capacité thermique effective C du local. La puissance injectée dans le local est compensée par la perte de chaleur à travers l'enveloppe et la chaleur emmagasinée dans la structure de l'enveloppe, ce qui est décrit par 30 l'équation : 15 P = HLC (Tai, in I,- Tai, 't k) + C dt où P est la puissance totale injectée dans le local, Tai, in k CITair in k et Tair out k sont respectivement la température de l'air à l'intérieur du local et la température de l'air à l'extérieur du local, HLC est le coefficient de déperdition thermique global du local et C est la capacité thermique effective du local. 5 On suppose que la réponse en température du local est une simple exponentielle décroissante et que sa constante de temps est le produit du coefficient de déperdition thermique global HLC du local et de la capacité thermique effective C du local. En réalité, la réponse du local est plus complexe et est la superposition d'un grand nombre d'exponentielles 10 décroissantes, mais il a été validé expérimentalement qu'en adaptant les conditions du test, en particulier la durée du test et la valeur de la puissance de chauffe appliquée dans le local, seule la constante de temps la plus grande joue un rôle et le modèle décrit précédemment est valide. En appliquant deux puissances P1 et P2 de chauffe du local de valeurs 15 différentes sur deux périodes de temps D1 et D2, il est alors possible de déterminer le coefficient de déperdition thermique global HLC du local selon l'équation : ai P2 - a2P1 alAT2in a2ATlin' où (a ) k=1 ou 2 est la pente de la tangente à l'évolution Tai,. in k(t) de la température de l'air à l'intérieur du local en fonction du temps, sur un 20 intervalle de temps Atk à l'intérieur de la période D k pour lequel l'évolution mesurée Tai,. in k(t) est sensiblement linéaire, et (AT ic=1 ou 2 est la différence entre la température moyenne de l'air à l'intérieur du local et la température moyenne de l'air à l'extérieur du local sur l'intervalle de temps Atk. 25 Par analogie, il est possible de déterminer le coefficient de transmission thermique U d'une paroi de séparation entre l'intérieur et l'extérieur du local selon la formule : HLC = 16 q2 - a2 alAT2rn a231Tirn' où (a ) Ic=1 ou 2 est la pente de la tangente à l'évolution T - air in k (t) de la température de l'air à l'intérieur du local en fonction du temps, sur un intervalle de temps Atk à l'intérieur de la période Dk pour lequel l'évolution mesurée Tair in k(t) est sensiblement linéaire, (AT ) Ic=1 ou 2 est la différence 5 entre la température moyenne de l'air à l'intérieur du local et la température moyenne de l'air à l'extérieur du local sur l'intervalle de temps Atk, et (q0k=1 ou 2 est le flux thermique moyen à travers la paroi pris sur la période de temps Dk ou, de préférence et pour plus de précision, pris sur l'intervalle de temps Atk. 10 Dans le cadre de l'invention, l'ajustement entre le modèle thermique et les mesures in situ peut se faire sur toute l'étendue de chaque période de temps Dk, ou sur un ou plusieurs intervalles de temps compris dans chaque période de temps Dk. De préférence, dans le cas où le modèle thermique est un modèle ARX, l'ajustement paramétrique est réalisé sur l'ensemble des 15 périodes de temps Dk. Selon une caractéristique avantageuse, pour chaque période de temps Dk, la puissance Pk de chauffe du local est appliquée au moyen d'au moins un appareil de chauffage ayant une source de puissance contrôlée. De préférence, le ou les appareils de chauffage sont peu inertes de manière 20 à assurer un chauffage rapide de l'air à l'intérieur du local. Au moins un appareil de chauffage à source de puissance contrôlée pour le chauffage du local peut être un équipement fixe du local, c'est-à-dire un appareil de chauffage installé dans le local indépendamment de la mise en oeuvre du procédé. Il peut notamment s'agir d'une pompe à chaleur dont 25 le coefficient de performance (COP) est connu. En variante, au moins un appareil de chauffage à source de puissance contrôlée pour le chauffage du local peut être un appareil rapporté dans le local spécifiquement pour la mise en oeuvre du procédé. Selon une autre variante, le chauffage du local sur chaque période de 30 temps Dk peut être mis en oeuvre à l'aide d'une combinaison d'au moins un U= 17 appareil de chauffage qui équipe le local de manière fixe, indépendamment de la mise en oeuvre du procédé, et d'au moins un appareil de chauffage rapporté dans le local spécifiquement pour la mise en oeuvre du procédé. Les appareils de chauffage du local peuvent être de type convectif, 5 radiatif ou par conduction, ou combiner plusieurs de ces technologies. De préférence, les appareils de chauffage du local sont des appareils électriques, ce qui permet de déterminer la puissance de chauffe de manière directe et précise. Des exemples d'appareils de chauffage électriques comprennent notamment des appareils de type convectif mettant en jeu le 10 soufflage d'air chauffé au moyen de résistances électriques ; des tapis ou des films chauffants ; des parasols radiants. En variante, les appareils de chauffage du local peuvent être des appareils fonctionnant au gaz ou au fioul, pour autant que les rendements des brûleurs et les débits de combustible puissent être estimés de manière suffisamment précise pour 15 accéder à la puissance de chauffe. Dans un mode de réalisation, les appareils de chauffage du local sont des petits convecteurs électriques répartis dans le local. Selon une variante, les appareils de chauffage du local sont des tapis chauffants électriques, que l'on répartit dans le local en les positionnant verticalement et enroulés sur 20 eux-mêmes, de sorte que toute la puissance thermique est dissipée dans l'air dans le local. Cet agencement permet un chauffage rapide et homogène de l'air à l'intérieur du local. Si le procédé de l'invention est mis en oeuvre avec un local comportant des cloisons internes qui délimitent plusieurs pièces ou zones du 25 local, on peut mesurer la température dans plusieurs pièces ou zones du local et considérer que la température dans le local à chaque temps t est la moyenne des mesures de température obtenues au temps t dans les différentes pièces ou zones du local pondérées chacune par le volume de la pièce ou zone. 30 Selon un aspect de l'invention, la puissance de chauffe délivrée dans le local est mesurée à l'aide d'au moins un capteur de puissance. Le ou chaque capteur de puissance peut être un capteur de tension (voltmètre) 18 et/ou un capteur de courant (ampèremètre). De préférence, le ou chaque capteur de puissance est un wattmètre, muni à la fois d'un capteur de tension et d'un capteur de courant. Cela permet une mesure précise de la puissance dans le local, en s'affranchissant d'éventuelles fluctuations de la 5 tension de secteur ou de la détermination de la résistance du ou de chaque appareil de chauffage. Dans le cadre de l'invention, on peut déterminer la température de l'air à l'extérieur du local par le biais d'une campagne de mesures à intervalles de temps rapprochés. En variante, la détermination de la température de l'air à 10 l'extérieur du local à intervalles de temps rapprochés peut être obtenue par interpolation de données météorologiques. De préférence, le procédé selon l'invention est mis en oeuvre sur des périodes de temps Di, pour lesquelles la température de l'air à l'extérieur du local est stable. 15 Selon un aspect avantageux de l'invention, on détermine l'indicateur de confort thermique à plusieurs positions dans le local et on en déduit une cartographie du confort thermique dans le local. Il est ainsi possible d'identifier des zones d'inconfort dans le local et de cibler les mesures qu'il convient de prendre pour améliorer le confort thermique. 20 Dans un mode de réalisation, au moins une partie des étapes du procédé selon l'invention sont déterminées par des instructions de programmes d'ordinateurs. En conséquence, l'invention a également pour objet un programme d'ordinateur sur un support d'enregistrement, ce programme étant 25 susceptible d'être mis en oeuvre dans un terminal informatique, ou plus généralement dans un ordinateur, ce programme comportant des instructions adaptées à l'exécution de tout ou partie des étapes d'un procédé tel que décrit ci-dessus. Ce programme peut utiliser n'importe quel langage de programmation, 30 et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée. 19 L'invention a aussi pour objet un support d'enregistrement lisible par un ordinateur, et comportant des instructions d'un programme d'ordinateur tel que mentionné ci-dessus. Le support d'enregistrement peut être n'importe quelle entité ou 5 dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une mémoire morte, une mémoire non volatile réinscriptible, par exemple une clé USB, une carte SD, une EEPROM, ou encore un moyen d'enregistrement magnétique, par exemple un disque dur. 10 Le support d'enregistrement peut aussi être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé. Le support d'enregistrement peut être un support transmissible tel qu'un signal électrique ou optique, qui peut être acheminé via un câble 15 électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau informatique. Un autre objet de l'invention est un dispositif de détermination prédictive d'un indicateur de confort thermique d'un local, pour des valeurs prédictives quelconques de la température de l'air à l'intérieur du local et de 20 la température de l'air à l'extérieur du local, ce dispositif étant caractérisé en ce qu'il comprend : - au moins un appareil de chauffage comportant une source de puissance contrôlée, configuré pour appliquer, sur au moins deux périodes de temps Di, successives, des puissances Pic distinctes de chauffe dans le 25 local ; - au moins un capteur de température configuré pour mesurer la température de l'air à l'intérieur du local Tai, mi, à intervalles de temps rapprochés sur les périodes de temps Di, ; - au moins un capteur de mesure configuré pour mesurer soit au 30 moins une température à l'intérieur du local, soit le flux thermique à travers des parois du local correspondant à l'indicateur de confort thermique, à intervalles de temps rapprochés sur les périodes de temps Di,; 20 - un terminal comportant un module de traitement configuré pour déterminer au moins un paramètre d'un modèle thermique descriptif du local, par ajustement entre le modèle thermique et les mesures réalisées dans le local à l'aide des capteurs du dispositif. 5 De manière avantageuse, le terminal comporte également un module d'acquisition pour acquérir les mesures requises par le procédé, le module de traitement étant configuré pour exécuter tout ou partie des étapes de calcul du procédé à partir des mesures acquises. Dans un mode de réalisation, le dispositif comprend également au 10 moins un capteur de température configuré pour mesurer la température de l'air à l'extérieur du local Tai, 'tic à intervalles de temps rapprochés sur les périodes de temps Di,. Le ou les capteurs de mesure de la température de l'air à l'intérieur du local, respectivement à l'extérieur du local, peuvent comprendre au moins un 15 capteur de température ambiante destiné à être positionné dans le volume d'air à l'intérieur du local, respectivement à l'extérieur du local, tel qu'un thermocouple, par exemple un thermocouple de type K ou de type T, ou un thermomètre à résistance, par exemple une sonde Pt100. Le ou les capteurs de mesure d'une température à l'intérieur du local 20 correspondant à l'indicateur de confort thermique peuvent comprendre, à titre d'exemples : - lorsque l'indicateur de confort thermique est la température opérative to en un point à l'intérieur du local, au moins un capteur de température moyenne radiante à l'intérieur du local destiné à être positionné en ce point, 25 tel qu'un thermomètre à globe noir, et/ou une pluralité de capteurs de température de surface à l'intérieur du local destinés à mesurer la température de surface de chaque élément de construction des parois du local vu en ce point ; - lorsque l'indicateur de confort thermique est la différence verticale de 30 la température de l'air Jta,, selon une verticale du local, au moins deux capteurs de température de l'air à l'intérieur du local destinés à être 21 positionnés en deux points de hauteurs différentes selon cette verticale à l'intérieur du local ; - lorsque l'indicateur de confort thermique est l'asymétrie de température de rayonnement Mi), verticale, respectivement horizontale, une 5 pluralité de capteurs de température de surface à l'intérieur du local destinés à mesurer la température de surface de chaque élément de construction des parois horizontales du local en vis-vis, respectivement des parois verticales du local en vis-vis. Le ou les capteurs de température de surface peuvent comprendre au 10 moins un capteur de température de surface destiné à être positionné sur la surface de l'élément de construction à l'intérieur du local, tel qu'un thermocouple fin ou un thermomètre à résistance plat. De manière alternative ou complémentaire, le ou les capteurs de température de surface peuvent comprendre au moins un capteur de rayonnement thermique destiné 15 à être positionné en regard de l'élément de construction, notamment un capteur de mesure ponctuelle, tels qu'un pyromètre, ou un imageur, tel qu'une caméra thermique, en particulier une caméra thermique infrarouge. Par ailleurs, le ou les capteurs de mesure du flux thermique à travers des parois du local correspondant à l'indicateur de confort thermique peuvent 20 comprendre une pluralité de capteurs de flux thermique destinés à mesurer le flux thermique à travers chaque élément de construction du local qui a un intérêt pour la détermination de l'indicateur de confort thermique, à titre d'exemples : - lorsque l'indicateur de confort thermique est la température opérative 25 to en un point à l'intérieur du local, une pluralité de capteurs de flux thermique destinés à mesurer le flux thermique à travers chaque élément de construction des parois du local vu en ce point ; - lorsque l'indicateur de confort thermique est l'asymétrie de température de rayonnement Mi), verticale, respectivement horizontale, une 30 pluralité de capteurs de flux thermique destinés à mesurer le flux thermique à 22 travers chaque élément de construction des parois horizontales, respectivement verticales, du local en vis-vis. Chaque capteur de flux thermique est destiné à être positionné sur une face de l'élément de construction ou de la paroi. Le capteur de flux 5 thermique peut être un fluxmètre ou un calorimètre. De manière avantageuse, le capteur de flux thermique est un fluxmètre compatible avec la norme ISO 9869:1994, notamment un fluxmètre à gradient. Selon un aspect de l'invention, la résistance au flux due au fluxmètre est prise en compte et une correction est appliquée au flux thermique 10 mesuré, de manière à obtenir le flux thermique concernant uniquement l'élément de construction ou la paroi dont on souhaite déterminer le coefficient de transmission thermique U. Cela est particulièrement important pour des parois à faible résistance thermique, telles que des simples vitrages. 15 Selon une caractéristique avantageuse, le ou chaque appareil de chauffage chauffe l'air à l'intérieur du local. Cela permet un chauffage rapide du local. Tel est le cas notamment avec une pluralité de convecteurs électriques répartis dans le local, ou avec des tapis chauffants électriques comme décrits précédemment, qui sont disposés verticalement dans le local 20 et enroulés sur eux-mêmes, de sorte que toute la puissance thermique est dissipée dans l'air. Selon une caractéristique avantageuse, le terminal comporte un module de pilotage de la source de puissance du ou de chaque appareil de chauffage. 25 Selon une autre caractéristique avantageuse, le dispositif comprend des moyens de liaison, notamment sans fil, entre le ou chaque capteur du dispositif et le terminal. Les caractéristiques et avantages de l'invention apparaîtront dans la description qui va suivre de deux exemples de réalisation d'un procédé et 30 d'un dispositif selon l'invention, donnée uniquement à titre d'exemple et faite en se référant aux figures annexées dans lesquelles : 23 - la figure 1 est une vue de dessus schématique de l'intérieur d'une chambre d'une maison, dont l'enveloppe comprend plusieurs éléments de construction, à savoir un sol, un plafond et quatre murs parmi lesquels celui orienté au sud est muni d'un vitrage, et où est installé un dispositif selon 5 l'invention pour la détermination d'une cartographie prédictive de la température opérative to, en tant qu'indicateur de confort thermique, pour des valeurs prédictives quelconques de la température de l'air à l'intérieur de la chambre et de la température de l'air à l'extérieur de la maison, conformément à l'invention ; 10 - la figure 2 est un graphe montrant l'évolution temporelle de la température de l'air à l'intérieur de la chambre Tai, ink(t), lors d'un test comprenant l'application d'une première puissance de chauffe /31 à l'intérieur de la chambre sur une première période de temps D1, suivie de l'application d'une deuxième puissance de chauffe P2 sensiblement nulle à l'intérieur de la 15 chambre sur une deuxième période de temps D2, de manière à laisser la chambre en refroidissement libre ; - la figure 3 est un graphe montrant l'évolution temporelle de la température de l'air à l'extérieur de la maison Tai, 'tic(t), lors du test décrit précédemment ; 20 - la figure 4 est un graphe montrant l'évolution temporelle de la température de surface à l'intérieur de la chambre du mur orienté à l'est, lors du test décrit précédemment ; - la figure 5 est une cartographie prédictive de la température opérative to, pour des valeurs prédictives quelconques de la température de 25 l'air à l'intérieur de la chambre et de la température de l'air à l'extérieur de la maison, obtenue conformément à l'invention à partir de la détermination du coefficient de déperdition thermique HLC de la chambre et du coefficient de transmission effectif Heff au niveau de chaque point de mesure de la température de surface à l'intérieur de la chambre ; 24 - la figure 6 est un graphe montrant l'évolution temporelle du flux thermique à travers le mur de la chambre orienté à l'est, lors du test décrit précédemment ; et - la figure 7 est une cartographie prédictive de la température 5 opérative to, pour des valeurs prédictives quelconques de la température de l'air à l'intérieur de la chambre et de la température de l'air à l'extérieur de la maison, obtenue conformément à l'invention à partir de la détermination du coefficient de transmission thermique U de chaque élément de construction de la chambre formant une séparation entre l'intérieur et l'extérieur de la 10 maison. PROTOCOLE EXPERIMENTAL Le procédé selon l'invention est mis en oeuvre pour la détermination d'une cartographie prédictive de la température opérative to à l'intérieur de la 15 chambre 1 d'une maison montrée sur la figure 1, pour des valeurs prédictives quelconques de la température de l'air à l'intérieur de la chambre et de la température de l'air à l'extérieur de la maison. Pour obtenir cette cartographie prédictive de la température opérative to, l'invention est mise en oeuvre selon deux exemples de réalisation, dont 20 les résultats sont visibles, respectivement, sur la figure 5 pour l'exemple 1 et sur la figure 7 pour l'exemple 2. Les mesures in situ dans la chambre 1 sont réalisées pour les deux exemples au cours d'une même sollicitation thermique de la chambre 1, telle que décrite ci-dessous. L'enveloppe de la chambre 1 comprend un sol, un plafond, un mur 13 25 orienté à l'ouest, un mur 14 orienté au nord, un mur 15 orienté à l'est, et un mur 16 orienté au sud, ce dernier étant muni d'un vitrage 17 qui est un double vitrage. La chambre 1 présente une surface au sol de 13,2 m2, une surface de vitrage de 1,6 m2, une hauteur intérieure de 2,3 m, un volume de 30,4 m3 et une surface totale d'enveloppe en contact avec l'extérieur de 30 30 m2. Le sol, le plafond et les murs 13 et 14 de la chambre 1 sont des parois internes de la maison, tandis que les murs 15 et 16 et le vitrage 17 sont des 25 éléments de construction formant une séparation entre l'intérieur et l'extérieur de la maison. Dans ce qui suit, on considère que les parois de la chambre 1 correspondant à des parois internes de la maison sont à l'équilibre, de sorte qu'au cours du test leur température de surface à 5 l'intérieur de la chambre est égale à la température de l'air à l'intérieur de la chambre et le flux thermique qui les traverse est nul. Pour les deux exemples, la sollicitation thermique de la chambre 1 est mise en oeuvre alors que la chambre 1 est inoccupée. De plus, la sollicitation thermique de la chambre 1 est mise en oeuvre en continu dans sa globalité 10 sur une seule période de temps nocturne, de manière à s'affranchir de la contribution du rayonnement solaire au chauffage de la chambre 1. Pour les deux exemples, la sollicitation thermique de la chambre 1 comprend une première phase de chauffage de la chambre sur une première période de temps D1, de 00h30 à 04h30, ce qui correspond à l'application 15 d'une première puissance de chauffe P1 strictement positive, puis une deuxième phase refroidissement libre de la chambre sur une deuxième période de temps D2 de 04h30 à 08h30, ce qui correspond à l'application d'une deuxième puissance de chauffe P2 sensiblement nulle. La deuxième période de temps D2 est immédiatement consécutive à la première période 20 de temps D1. L'exemple 1 illustre l'obtention de la cartographie de température opérative to à l'intérieur de la chambre 1 à partir de la détermination du coefficient de déperdition thermique global HLC de la chambre et du coefficient de transmission effectif Heff au niveau de chaque élément de 25 construction 15, 16, 17 de la chambre formant une séparation entre l'intérieur et l'extérieur de la maison. Pour ce faire, un thermocouple 50 de type K est positionné sur la face intérieure de chacun des éléments de construction 15, 16, 17 de manière à mesurer l'évolution temporelle de la température de surface de l'élément à l'intérieur de la chambre au cours du test. 30 L'exemple 2 illustre l'obtention de la cartographie de température opérative to à l'intérieur de la chambre 1 à partir de la détermination du 26 coefficient de transmission thermique U de chaque élément de construction 15, 16, 17 de la chambre formant une séparation entre l'intérieur et l'extérieur de la maison, en utilisant des mesures du flux thermique à travers ces éléments de construction. Pour ce faire, un fluxmètre 60 est positionné sur la 5 face intérieure de chacun des éléments de construction 15, 16, 17 de manière à mesurer l'évolution temporelle du flux thermique à travers l'élément au cours du test. Plus précisément, pour les deux exemples, le procédé selon l'invention est mis en oeuvre à l'aide du dispositif montré sur la figure 1, qui 10 comprend : - une pluralité de convecteurs électriques 20, qui sont positionnés au voisinage du centre de la chambre 1 pour le chauffage de l'air à l'intérieur de la chambre ; - un thermocouple 30 de type K, positionné dans l'air ambiant dans la 15 chambre 1, au milieu du volume d'air à 110 cm de hauteur, configuré pour mesurer la température de l'air à l'intérieur de la chambre ; - un thermocouple 40 de type K, placé dans l'air à l'extérieur de la maison, configuré pour mesurer la température de l'air à l'extérieur de la maison ; - les trois thermocouples 50 de type K, positionnés chacun sur la face 20 intérieure de l'un des éléments de construction 15, 16, 17 de la chambre ; - les trois fluxmètres 60, qui sont des fluxmètres à gradient de type HFP01 commercialisés par la société Hukseflux, positionnés chacun sur la face intérieure de l'un des éléments de construction 15, 16, 17 de la chambre ; - un terminal 70, par exemple un ordinateur de type PC portable, qui 25 comporte un module de traitement configuré pour réaliser un ajustement entre un modèle thermique descriptif de la chambre et les évolutions temporelles des températures et flux thermiques mesurées par les capteurs 30, 40, 50, 60. De manière avantageuse, le terminal 70 comprend un module de 30 pilotage des convecteurs électriques 20, pour appliquer la modulation temporelle de la puissance de chauffe de l'air à l'intérieur de la chambre 1 sur les périodes de temps D1 et D2. 27 Pour chaque période de temps Dic, la puissance Pic appliquée est sensiblement égale à la puissance de chauffe imposée par les convecteurs électriques 20, aux puissances résiduelles près, provenant notamment du matériel de mesure et de calcul présent dans la chambre pour la mise en 5 oeuvre du procédé. Des capteurs de puissance, sous la forme de boucles ampèremétriques, peuvent être prévus pour mesurer la puissance délivrée dans la chambre lors de la mise en oeuvre du procédé. Au cours des deux périodes de temps D1 et D2, des mesures sont réalisées toutes les dix secondes à l'aide des capteurs 30, 40, 50, 60. Les 10 résultats de mesures sont visibles sur la figure 2 pour le thermocouple 30 de mesure de la température de l'air à l'intérieur de la chambre ; sur la figure 3 pour le thermocouple 40 de mesure de la température de l'air à l'extérieur de la maison ; sur la figure 4 pour le thermocouple 50 de mesure de la température de surface à l'intérieur de la chambre du mur 15 orienté à l'est 15 (montré à titre d'exemple, étant entendu que les résultats de mesures des thermocouples 50 positionnés sur le mur 16 et le vitrage 17 ont aussi été obtenus) ; sur la figure 6 pour le fluxmètre 60 de mesure du flux thermique à travers le mur 15 de la chambre orienté à l'est (montré à titre d'exemple, étant entendu que les résultats de mesures des fluxmètres 60 positionnés 20 sur le mur 16 et le vitrage 17 ont aussi été obtenus). Les données de mesures brutes obtenues à l'aide des capteurs 30, 40, 50, 60 sont acquises par un module d'acquisition du terminal 70 et traitées par le module de traitement du terminal 70, qui exécute des instructions d'un programme d'ordinateur conforme à l'invention installé dans 25 le terminal pour réaliser l'ajustement paramétrique entre un modèle thermique descriptif de la chambre et les évolutions temporelles mesurées des températures et flux thermiques. EXEMPLE 1 30 Le coefficient de déperdition thermique global HLC de la chambre 1 est déterminé par ajustement ("fitting"), sur chaque période de temps D1 et D2, entre d'une part un modèle R-C simple avec une résistance et une 28 capacité, exprimant la variation temporelle de la température de l'air à l'intérieur de la chambre en fonction de la puissance de chauffe appliquée dans la chambre et de la température de l'air à l'extérieur de la maison, et d'autre part l'évolution mesurée de la température de l'air à l'intérieur de la 5 chambre en fonction du temps Tai, in k(t). Comme expliqué précédemment, la valeur du coefficient de déperdition thermique global HLC de la chambre 1 est obtenue selon la formule : ai P2 - a2P1 alAT2in a2Anin' où (a ) k=1 ou 2 est la pente de la tangente à l'évolution Tai, k(t) de la 10 température de l'air à l'intérieur de la chambre en fonction du temps, sur un intervalle de temps Atk à l'intérieur de la période Dk pour lequel l'évolution mesurée Tai, k(t) est sensiblement linéaire, et (AT km)1c=1 ou 2 est la différence entre la température moyenne de l'air à l'intérieur du local et la température moyenne de l'air à l'extérieur du local sur l'intervalle de temps 15 Atk, et (-Pk)k.i ou 2 est la puissance de chauffe appliquée dans la chambre sur la période de temps Dk. La figure 2 montre les périodes de temps D, et D2, ainsi que les intervalles de temps At, et At2 sur lesquels a été réalisée l'approximation linéaire, c'est-à-dire l'ajustement avec le modèle R-C. Le coefficient de transmission effectif Heff au niveau de chaque point 20 de mesure de la température de surface à l'intérieur de la chambre 1 est également déterminé, pour chaque point de mesure, par ajustement ("fitting"), sur chaque période de temps D, et D2, entre d'une part un modèle R-C simple avec une résistance et une capacité, exprimant la variation temporelle de la température de surface à l'intérieur de la chambre en 25 fonction de la puissance de chauffe appliquée dans la chambre et de la température de l'air à l'extérieur de la maison, et d'autre part l'évolution mesurée de la température de surface à l'intérieur de la chambre en fonction du temps Tsurf. HLC = in k(t)- 29 Comme précédemment, par analogie, la valeur du coefficient de transmission effectif Heff au niveau de chaque élément de construction de la chambre 1 formant une séparation entre l'intérieur et l'extérieur de la maison est obtenue selon la formule : H = ef f A wiTt-^ 2 sin a2ATism' 5 où (a ) k=1 ou 2 est la pente de la tangente à l'évolution T - sur f in k (t) de la température de surface de l'élément de construction à l'intérieur de la chambre en fonction du temps, sur un intervalle de temps Atk à l'intérieur de la période Dk pour lequel l'évolution mesurée TsurfinI c ( t) est sensiblement linéaire, et (ATksin)k=lou 2 est la différence entre la température moyenne de 10 surface à l'intérieur de la chambre et la température moyenne à l'extérieur de la maison sur l'intervalle de temps Atk, et (P k=1 ou 2 est la puissance de chauffe appliquée dans la chambre sur la période de temps Dk. La figure 4 montre les périodes de temps D1 et D2, ainsi que les intervalles de temps At, et At2 sur lesquels a été réalisée l'approximation linéaire, c'est-à-dire 15 l'ajustement avec le modèle R-C. Les résultats sont donnés dans le tableau 1 ci-dessous. Point de mesure P1 al ATi(s)m P2 a2 AT2(s)n H L C OU (W) (°C.jourl) (°c) (W) (°C.jourl) (°c) Heff (W.K1) Air 110 cm 6756 11,1 24,8 132 -16 14,1 198 Mur 15 6756 9,1 16,5 132 -12,9 11,8 276 Mur 16 6756 7,7 16,7 132 -12,4 11,9 285 Vitrage 17 6756 3,7 8,8 132 -6,9 5 596 Tableau 1 20 ai P2 - a2P1 Pour chaque élément de construction des parois de la chambre, la température de surface prédictive Tsurfpred de l'élément est alors 30 déterminée, pour des valeurs prédictives quelconques de la température de l'air à l'intérieur de la chambre Tairinpred et de la température de l'air à l'extérieur de la maison Tair 't pred, selon la relation : HLC Tsurf pred = Tair out pred u (Tair in pred Tair out pred) ieff La température opérative to à l'intérieur de la chambre est ensuite 5 déterminée, par itération pour des points successifs situés par exemple tous les 50 cm dans la chambre 1, selon une relation classique, bien connue de l'homme du métier, exprimant la température opérative to en fonction de la température de l'air à l'intérieur de la chambre et de la température de surface de chaque élément de construction vu en ce point de la chambre, 10 telle que par exemple la relation : 1 1 t0 = -2Tair in pred -2Tr avec i chaque élément de construction des parois de la chambre, Fi le facteur de vue de l'élément de construction i en chaque point considéré à l'intérieur de la maison, Tsurfpreai la température de surface prédictive de l'élément de construction j. 15 Dans cet exemple (qui est un exemple non limitatif, d'autres formules de facteur de vue peuvent bien entendu être utilisées dans le cadre de l'invention), on utilise pour le facteur de vue la formule donnée sur le site http://www.thermalradiation.net/sectionc/C-62. htm I, correspondant à un facteur de vue de l'extérieur d'un cylindre infiniment long vis-à-vis de plaques 20 parallèles infiniment longues, avec des cotes a, b1, b2 telles que montrées à titre d'exemple sur la figure 1 pour un occupant (élément 1 pour la formule, représenté par « x » sur la figure 1) vis-à-vis du mur 16 de la chambre 1 exposé au sud (élément 2 pour la formule), où r le rayon du cylindre correspondant à l'occupant est tel que r REVENDICATIONS1. Procédé de détermination prédictive d'un indicateur de confort thermique (to ; Jta,, ; Jtpr) d'un local, pour des valeurs prédictives quelconques de la température de l'air à l'intérieur du local (Tair in pred) et de la température de l'air à l'extérieur du local (Tair 't pred), caractérisé en ce qu'il comprend des étapes dans lesquelles : - on détermine expérimentalement au moins un paramètre (HLC, Hely, U) d'un modèle thermique descriptif du local, par ajustement entre le modèle thermique et des mesures réalisées in situ dans le local, où les mesures comprennent des mesures de la température de l'air à l'intérieur du local (Tair in k) et des mesures soit d'au moins une température supplémentaire à l'intérieur du local (Trad in k conf ; Tair in k conf ; Tsurf in k conf), soit du flux thermique (qk) à travers des parois du local correspondant à l'indicateur de confort thermique (to ; Jta,, ; Jtpr) ; - on calcule une ou des températures de confort prédictives (Tconfortpred ; Tsurf pred) intervenant pour la détermination de l'indicateur de confort thermique (to ; Jta,, ; Jtpr), pour lesdites valeurs prédictives quelconques de la température de l'air à l'intérieur du local (Tair in pred) et de la température de l'air à l'extérieur du local (Tair out pred). 2. Procédé selon la revendication 1, caractérisé en ce qu'il comprend des étapes dans lesquelles : - sur au moins deux périodes de temps Dk successives correspondant à des puissances Pk distinctes de chauffe appliquées à l'intérieur du local, on procède à une campagne de mesures de la température de l'air à l'intérieur du local (Tair in k) et d'au moins une température à l'intérieur du local (Trad in k conf ; Tair in k conf ; Tsurf in k conf) correspondant à l'indicateur de confort thermique (to ; Jta,, ; Jtpr) à intervalles de temps rapprochés, ainsi qu'à la détermination de la température de l'air à l'extérieur du local (Tair out k) à intervalles de temps rapprochés ; 36 - on détermine la valeur du coefficient de déperdition thermique HLC du local, par ajustement entre : d'une part, un modèle thermique exprimant la variation temporelle de la température de l'air à l'intérieur du local en fonction de la puissance de 5 chauffe appliquée dans le local et de la température de l'air à l'extérieur du local, et d'autre part, l'évolution mesurée de la température de l'air à l'intérieur du local (Tair in k) en fonction du temps ; - on détermine la valeur d'un coefficient de transmission effectif Heff 10 en chaque point de mesure de la température à l'intérieur du local correspondant à l'indicateur de confort thermique (to ; Jta,, ; Jtpr), par ajustement entre : d'une part, un modèle thermique exprimant la variation temporelle de la température à l'intérieur du local correspondant à l'indicateur de confort 15 thermique (t0; Jta,, ; Jtpr) en fonction de la puissance de chauffe appliquée dans le local et de la température de l'air à l'extérieur du local, et d'autre part, l'évolution mesurée de la température à l'intérieur du local (Trad in k conf ; Tair in k conf ; Tsurf in k conf) correspondant à l'indicateur de confort thermique (to ; Jta,, ; Jtpr) en fonction du temps ; 20 - on calcule la ou les températures de confort prédictives Teenfort pred intervenant pour la détermination de l'indicateur de confort thermique (t0; Jta,, ; Jtpr), à partir des valeurs préalablement déterminées du coefficient de déperdition thermique HLC du local et du coefficient de transmission effectif Hely, pour lesdites valeurs prédictives quelconques de la température de l'air 25 à l'intérieur du local Tair inpred et de la température de l'air à l'extérieur du local Tair out pred, selon la relation : HLC Tconfort pred = Tair out pred u Tair in pred Tair out pred) iief f 3. Procédé selon la revendication 2, caractérisé en ce que l'indicateur de confort thermique est la température opérative to en un point à l'intérieur 37 du local, et on détermine la valeur du coefficient de transmission effectif Heff en ce point à partir de mesures de l'évolution temporelle de la température moyenne radiante à l'intérieur du local (T'din k conf) en ce point et/ou on détermine la valeur du coefficient de transmission effectif Heff pour chaque 5 élément de construction des parois du local vu en ce point à partir de mesures de l'évolution temporelle de la température de surface à l'intérieur du local (T'rf in k conf) de chacun de ces éléments de construction. 4. Procédé selon la revendication 2, caractérisé en ce que l'indicateur de confort thermique est la différence verticale de la température de l'air Jta,, 10 selon une verticale du local, et on détermine les valeurs du coefficient de transmission effectif Heff en deux points de hauteurs différentes selon cette verticale à l'intérieur du local à partir de mesures de l'évolution temporelle de la température de l'air à l'intérieur du local (l'air in k conf.) en ces deux points. 5. Procédé selon la revendication 2, caractérisé en ce que l'indicateur 15 de confort thermique est l'asymétrie de température de rayonnement Jtpr verticale, respectivement horizontale, et on détermine la valeur du coefficient de transmission effectif Heff de chaque élément de construction des parois horizontales du local en vis-vis, respectivement des parois verticales du local en vis-vis, à partir de mesures de l'évolution temporelle de la température de 20 surface à l'intérieur du local (T'rf in k conf) de chacune de ces parois. 6. Procédé selon la revendication 1, caractérisé en ce qu'il comprend des étapes dans lesquelles : - on détermine expérimentalement le coefficient de transmission thermique U de chaque paroi du local correspondant à l'indicateur de confort 25 thermique (to ; Jta,, ; Jtpr), par ajustement entre un modèle thermique et des mesures réalisées in situ dans le local, où les mesures comprennent des mesures de la température de l'air à l'intérieur du local (Tai, in k) et des mesures soit du flux thermique (qk) à travers la paroi, soit de la température de surface de la paroi à l'intérieur du local (T'rj in k) ; 38 - on calcule une température de surface prédictive T'rf. p'd de chaque paroi du local correspondant à l'indicateur de confort thermique (to ; Jta,, ; Jtpr), pour lesdites valeurs prédictives quelconques de la température de l'air à l'intérieur du local Tair inpred et de la température de l'air à 5 l'extérieur du local Tair 't prod, selon la relation : U Tsurf pred = Tair in pred Tair in pred Tair out pred) avec hin le coefficient d'échange convectif de l'air à l'intérieur du local. 7. Procédé selon la revendication 6, caractérisé en ce que, pour déterminer expérimentalement le coefficient de transmission thermique U de chaque paroi du local correspondant à l'indicateur de confort thermique (to ; 10 Jta,, ; Jtpr), on procède aux étapes suivantes : - sur au moins deux périodes de temps Di, successives correspondant à des puissances Pic distinctes de chauffe appliquées dans le local, on procède à une campagne de mesures de la température de l'air à l'intérieur du local Tair in k et du flux thermique à travers la paroi qk à intervalles de 15 temps rapprochés, ainsi qu'à la détermination de la température de l'air à l'extérieur du local Tair out k à intervalles de temps rapprochés ; - on détermine la valeur du coefficient de transmission thermique U de la paroi par ajustement entre : d'une part, un modèle thermique exprimant la variation temporelle de la 20 température de l'air à l'intérieur du local en fonction du flux thermique à travers la paroi et de la température de l'air à l'extérieur du local, et d'autre part, l'évolution mesurée de la température de l'air à l'intérieur du local Tair in k ou du flux thermique à travers la paroi qk en fonction du temps. 25 8. Procédé selon la revendication 6, caractérisé en ce que, pour déterminer expérimentalement le coefficient de transmission thermique U de chaque paroi du local correspondant à l'indicateur de confort thermique (to ; Jta,, ; Jtpr), on procède aux étapes suivantes : 39 - sur au moins deux périodes de temps Di, successives correspondant à des puissances /3k distinctes de chauffe appliquées dans le local, on procède à une campagne de mesures de la température de l'air à l'intérieur du local Tai, in I, et de la température de surface de la paroi à l'intérieur du 5 local T'rf. k à intervalles de temps rapprochés, ainsi qu'à la détermination de la température de l'air à l'extérieur du local Tai, 't I, à intervalles de temps rapprochés ; - on détermine la valeur de la grandeur U I hm, où hi, est le coefficient d'échange convectif de l'air dans le local au voisinage de la paroi, par 10 ajustement entre : d'une part, un modèle thermique exprimant la variation temporelle de la température de l'air à l'intérieur du local en fonction de la température de surface de la paroi à l'intérieur du local et de la température de l'air à l'extérieur du local, et 15 d'autre part, l'évolution mesurée de la température de l'air à l'intérieur du local Tai, in I, ou de la température de surface de la paroi à l'intérieur du local T'rj in k en fonction du temps. 9. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que les mesures in situ sont mises en oeuvre sur deux 20 périodes de temps D1 et D2 successives correspondant à deux consignes de puissances distinctes de chauffe P1 et P2 appliquées dans le local. 10. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le modèle thermique est un modèle d'identification paramétrique, notamment un modèle ARX. 25 11. Procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce que le modèle thermique est un modèle R-C. 12. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on détermine l'indicateur de confort thermique à plusieurs positions dans le local et on déduit une cartographie du confort 30 thermique dans le local. 40 13. Programme d'ordinateur comportant des instructions pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 12 lorsque ledit programme est exécuté par un ordinateur. 14. Support d'enregistrement lisible par un ordinateur sur lequel est 5 enregistré un programme d'ordinateur comprenant des instructions pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 12. 15. Dispositif de détermination prédictive d'un indicateur de confort thermique (to ; Jta,, ; Jtpr) d'un local, pour des valeurs prédictives 10 quelconques de la température de l'air à l'intérieur du local (Tai,. in pred) et de la température de l'air à l'extérieur du local (Tati. 't pred), ce dispositif étant caractérisé en ce qu'il comprend : - au moins un appareil de chauffage (20) comportant une source de puissance contrôlée, configuré pour appliquer, sur au moins deux périodes 15 de temps Di, successives, des puissances Pk distinctes de chauffe dans le local ; - au moins un capteur de température (30) configuré pour mesurer la température de l'air à l'intérieur du local Tair in k à intervalles de temps rapprochés sur les périodes de temps Di, ; 20 - au moins un capteur de mesure (50; 60) configuré pour mesurer soit au moins une température supplémentaire à l'intérieur du local (Trad in k conf ; Tair in k conf ; Tsurf in k conf), soit le flux thermique (qk) à travers des parois du local correspondant à l'indicateur de confort thermique (to ; Jta,, ; Jtpr), à intervalles de temps rapprochés sur les périodes de temps 25 Dk ; - un terminal (70) comportant un module de traitement configuré pour déterminer au moins un paramètre (HLC, Hely, U) d'un modèle thermique descriptif du local, par ajustement entre le modèle thermique et les mesures réalisées dans le local à l'aide desdits capteurs (30, 50, 60) du dispositif. 30 16. Dispositif selon la revendication 15, caractérisé en ce que le ou chaque capteur (50) de mesure d'au moins une température supplémentaire 41 à l'intérieur du local (T'd in k con f ; Tair in k conf ; Tsurf in k conf) correspondant à l'indicateur de confort thermique (to ; Jta,, ; Jtpr) est un capteur de température moyenne radiante, un capteur de température d'air ou un capteur de température de surface, sélectionné en fonction de l'indicateur de 5 confort thermique (to ; Jta,, ; Jtpr).