FR 3038716 A1 20170113 FR 1556383 A 20150706 PROCEDE DE DETERMINATION D'UN DIAGNOSTIC THERMIQUE D'UN BÂTIMENT Domaine technique L'invention concerne un procédé de détermination d'un diagnostic thermique par la mesure d'une température intérieure d'un bâtiment et un dispositif pour la mise en oeuvre d'un tel procédé. Etat de la technique La connaissance du comportement thermique d'un bâtiment permet d'évaluer les besoins de rénovation énergétique éventuels, mais aussi d'améliorer le contrôle de la température à l'intérieur de ce bâtiment, notamment pour améliorer le confort des occupants et/ou réduire la consommation énergétique. Le comportement thermique d'un bâtiment est classiquement évalué par des mesures sur le bâtiment. Ces mesures ne permettent cependant qu'une évaluation très approximative, ne prenant pas forcément en compte toutes les réalités du terrain, et ne peuvent donc conduire qu'à des préconisations très générales. Il existe donc un besoin pour un procédé permettant de fournir des préconisations techniques plus fiables et précises. La présente invention vise à satisfaire, au moins partiellement, ce besoin. Résumé de l'invention L'invention propose un procédé de détermination d'un diagnostic thermique par la mesure d'une température intérieure d'un bâtiment, le procédé comportant les étapes successives suivantes : a) choix d'un modèle thermique initial établissant, au moyen d'un ensemble de paramètres, une relation théorique entre ladite température intérieure et un ensemble de facteurs d'influence de ladite température intérieure ; b) détermination d'un ensemble de valeurs réelles de ladite température intérieure et desdits facteurs d'influence ; 3038716 2 c) optimisation des paramètres du modèle thermique initial par analyse, au moyen d'une méthode d'analyse, dudit ensemble de valeurs réelles, optionnellement filtrées au moyen d'un filtre ; d) en fonction d'une fonction d'évaluation, modification dudit ensemble de valeurs 5 réelles et/ou dudit modèle thermique initial et/ou dudit filtre et/ou de ladite méthode d'analyse, et reprise à l'étape c) et/ou poursuite à l'étape e) ; e) utilisation du modèle thermique pour réaliser un diagnostic thermique et, de préférence, déterminer des préconisations techniques, en particulier pour la gestion du chauffage, les habitudes d'usage et/ou la rénovation énergétique du bâtiment. 10 Comme on le verra plus en détail dans la suite de la description, grâce à l'étape d), un procédé selon l'invention est évolutif et permet d'effectuer un diagnostic précis et fiable, et de préférence de rechercher des préconisations techniques mieux adaptées à la situation réelle. Un procédé selon l'invention peut encore comporter une ou plusieurs des caractéristiques optionnelles préférées suivantes : 15 la modification à l'étape d) dépend - de la valeur et/ou de l'évolution temporelle d'un ou plusieurs paramètres ; et/ou - d'un score représentant la correspondance des valeurs des paramètres du modèle initial optimisé à l'étape c) avec lesdites valeurs réelles ; et/ou - de l'ensemble de valeurs réelles utilisé à l'étape c) ; et/ou 20 - du modèle thermique initial. le modèle thermique initial est choisi dans une base de modèles types en fonction d'une ou plusieurs caractéristiques thermiques dudit bâtiment, par exemple en fonction du nombre et/ou du type d'ouvertures et/ou de l'âge du bâtiment ; - à l'étape b) et/ou à l'étape d), les valeurs réelles de la température intérieure du bâtiment et 25 des facteurs d'influence dudit ensemble sont mesurées dans une période d'analyse dont le début et la fin sont déterminés en fonction de l'évolution des facteurs d'influence ; - à l'étape c), les valeurs réelles de la température intérieure du bâtiment et des facteurs d'influence dudit ensemble sont filtrées ; le filtre appliqué dépend de l'objet auquel il est appliqué ; le filtre appliqué dépend de son efficacité à améliorer les valeurs des paramètres ; la modification du filtre dépend de la valeur d'un ou plusieurs paramètres ; 3038716 3 - à l'étape b) et/ou à l'étape d), les valeurs réelles de la température intérieure du bâtiment et des facteurs d'influence dudit ensemble sont mesurées dans une période d'analyse dont la durée est supérieure 1 heure, de préférence supérieure à 3 heures, et/ou inférieure à 24 heures, de préférence inférieure à 12 heures ; 5 - à l'étape b) et/ou à l'étape d), les valeurs réelles de la température intérieure du bâtiment et des facteurs d'influence dudit ensemble sont mesurées dans une période d'analyse pendant laquelle tous les facteurs d'influence sauf un, voire plusieurs, dit « facteur d'influence actif », sont sensiblement constants ; l'étape c) est exécutée plusieurs fois, chaque exécution étant réalisée à partir d'un 10 ensemble de valeurs réelles de la température intérieure du bâtiment et des facteurs d'influence mesurées dans une période d'analyse pendant laquelle tous les facteurs d'influence sauf un, voire plusieurs, dit « facteurs d'influence actifs », sont sensiblement constants, les facteurs d'influence actifs ayant été successivement analysés du fait de l'exécution répétée de l'étape c) comportant au moins 2, de préférence au moins 3, de 15 préférence au moins 4 facteurs d'influence différents, voire tous les facteurs d'influence ; les facteurs d'influence actifs comportent au moins 2, de préférence au moins 3, de préférence au moins 4 facteurs d'influence différents parmi la température extérieure, la vitesse du vent, l'orientation du vent, l'ensoleillement, la pression atmosphérique, l'humidité, la présence des occupants dans le bâtiment, la puissance et l'état de 20 fonctionnement du système de chauffage ; - à l'étape c), notamment pour tenir compte de l'inertie thermique du bâtiment, on compare la température intérieure réelle mesurée à un instant tl avec une température intérieure théorique obtenue avec des valeurs de facteurs d'influence mesurées à un instant t2 antérieur à tl, le décalage entre tl et t2 étant variable selon le facteur d'influence 25 considéré ; - à l'étape b), au moins une partie des facteurs d'influence est de préférence évaluée au moyen de réponses à au moins un questionnaire soumis à des occupants dudit bâtiment, un questionnaire étant de préférence soumis accompagné d'informations résultant d'une mise en oeuvre antérieure dudit procédé ; 30 - à l'étape c), on met en oeuvre une méthode d'analyse choisie dans le groupe constitué par l'identification, les statistiques, la méthode des moindres carrés, linéaire ou non, le filtre de Kalman, le filtre particulaire, les réseaux de neurones artificiels et le problème inverse. 3038716 4 - à l'étape e), on génère un rapport comportant un diagnostic exposant - l'environnement du bâtiment, en particulier l'environnement urbain ou l'exposition au vent ou au soleil ; et/ou - une évaluation de la résistance thermique globale et/ou de l'efficacité des menuiseries 5 et/ou des volets et/ou du chauffage et/ou de la ventilation et/ou des apports de chaleur interne et/ou du thermostat ; et/ou des préconisations pour : - isoler ou renforcer l'isolation en des emplacements déterminés du bâtiment, par exemple changer les volets ; et/ou 10 - modifier le chauffage, par exemple en préconisant des réglages de la puissance de la chaudière et/ou des durées et/ou des plages horaires de chauffage ; et/ou - préconiser des changements dans la gestion, par exemple de la ventilation et/ou de l'ouverture des volets et/ou de l'aération naturelle, et/ou des modifications dans l'emplacement du thermostat ; et/ou 15 - préconiser l'installation de capteurs et/ou de régulateurs pour affiner le modèle ; et/ou - préconiser l'installation de systèmes de contrôle, notamment de contrôle automatique, pour gérer plus efficacement les consommateurs d'énergie et/ou les moyens isolant, comme les volets ; et/ou - fournir des informations facilitant la mise en oeuvre des préconisations, comme des 20 adresses de prestataires et/ou des fiches techniques et/ou des informations sur les coûts à prévoir. L'invention concerne également : - un programme d'ordinateur comprenant des instructions de code de programme pour l'exécution d'une ou plusieurs, de préférence toutes les étapes d'un procédé selon 25 l'invention, lorsque ledit programme est exécuté par un ordinateur, - un support informatique sur lequel est enregistré un tel programme, par exemple une mémoire ou un CD-ROM, et - un module de traitement dans lequel est chargé un tel programme et apte à exécuter ledit programme. 30 3038716 5 Brève description des figures D'autres caractéristiques et avantages de l'invention apparaitront encore à la lecture de la description détaillée qui va suivre, et à l'examen du dessin annexé dans lequel : - la figure 1 représente, schématiquement, un exemple de dispositif selon l'invention ; 5 - la figure 2 représente l'évolution temporelle de différents facteurs d'influence F1-F4, à savoir le type de temps (pluie, nuageux, clair), la vitesse du vent, l'ensoleillement et la température extérieure, respectivement, et de la température intérieure Ti, t désignant le temps et J1-J5 désignant cinq jours consécutifs ; - la figure 3 représente le pourcentage du temps de chauffe (Pchauffe) en fonction de la 10 différence entre la température intérieure et la température extérieure (Ti-Te). Définitions - Un « bâtiment » désigne toute construction dans laquelle on souhaite effectuer un diagnostic thermique. En particulier, un bâtiment peut être un bâtiment habité, comme un appartement ou une maison, un bureau, un entrepôt ou tout autre bâtiment industriel. 15 - Un « modèle thermique » d'un bâtiment est un modèle qui, au moyen d'un ensemble de « paramètres », établit une relation entre une température intérieure au bâtiment et un ensemble de « facteurs d'influence » de la dite température intérieure. La température intérieure peut être en particulier une température mesurée au centre du bâtiment. 20 - Par « facteur d'influence », on désigne un facteur, variable dans le temps, susceptible d'agir sur la température intérieure. - Les « facteurs d'influence » comprennent les « facteurs météorologiques » et les « facteurs d'usage » du bâtiment. Les « facteurs météorologiques » sont les facteurs associés à l'environnement auquel le 25 bâtiment est soumis, et en particulier comprennent la température extérieure, la vitesse du vent, la direction du vent, le niveau d'ensoleillement, la pluviométrie, la couverture nuageuse et la pression atmosphérique. Les « facteurs d'usage » sont les facteurs liés à l'utilisation du bâtiment. Ils comportent en particulier le nombre d'occupants dans le bâtiment, la nature desdits occupants (femmes, 30 hommes, enfants, personnes âgées, animaux), et l'énergie apportée par le chauffage 3038716 6 (réglage de la chaudière) ou les autres sources de chaleur interne (cuisson, appareils ménagers, consommations de veille des appareils électroniques, par exemple). Les « caractéristiques thermiques » du bâtiment sont les caractéristiques sensiblement constantes dans le temps (à l'échelle des variations de la température intérieure) et qui 5 modifient l'action des facteurs d'influence sur la température intérieure du bâtiment. Ils peuvent être liés à la structure du bâtiment, comme la nature, la taille, l'orientation et le nombre d'ouvertures (portes, fenêtres simple vitrage, fenêtres double vitrage), la présence d'une Ventilation Mécanique Contrôlée (VMC), le nombre et le volume des pièces, les surfaces et les épaisseurs des murs, les épaisseurs d'isolant sur les murs et/ou sur la toiture, 10 l'existence d'un vide sanitaire ou d'une cave, le nombre d'étages, la hauteur des plafonds, le nombre de radiateurs, les positions des radiateurs. Ils peuvent être également liés à l'environnement physique du bâtiment, comme l'existence d'un autre bâtiment à proximité immédiate, voire accolé au bâtiment étudié, les caractéristiques de cet autre bâtiment (dimensions, position par rapport au bâtiment étudié), 15 l'environnement urbain ou rural, la présence de végétation et l'ampleur de cette végétation, le relief du terrain, la présence d'une route à proximité, les caractéristiques de cette route, comme le trafic en particulier. - Sauf indications contraires, par « comportant un », on entend « comprenant au moins un ». Description détaillée 20 Dispositif La figure 1 représente un bâtiment 10 comportant une chaudière 12, et un dispositif selon l'invention pour la mise en oeuvre d'un procédé selon l'invention. Le dispositif comporte : un thermomètre 16 disposé pour mesurer, à tout instant t, une température intérieure 25 Ti(t) à l'intérieur du bâtiment 10, un ensemble de capteurs 18' destinés chacun à mesurer, à tout instant t, un facteur d'influence de ladite température intérieure, un module de traitement 22. L'indice « i » est utilisé pour distinguer chacun des capteurs 18, représentés, ainsi que pour 30 désigner les données correspondantes. 3038716 7 Les capteurs peuvent être dans le bâtiment (capteur 181) ou à proximité du bâtiment (capteur 182), ou être à distance, c'est-à-dire typiquement à plus de 10 m du bâtiment. De préférence, un centre météorologique 24 traite les mesures fournies par les capteurs à distance (capteurs 183, 184, 185). 5 Les capteurs pour mesurer les facteurs d'influence sont connus et comprennent, par exemple, pour les facteurs météorologiques, un anémomètre, un thermomètre, ou un baromètre, ou, pour les facteurs d'usage, un détecteur de présence, un compteur d'évaluation de la puissance de chauffage ou un capteur 186 de mesure du fonctionnement de la chaudière 12. Dans un mode de réalisation, les facteurs d'usage peuvent être mesurés au moyen d'un 10 téléphone portable de type « smartphone », en particulier pourvu de moyens de géolocalisation ou d'un accéléromètre, ou par analyse d'un emploi du temps des occupants, communiqué par exemple par synchronisation d'un agenda électronique. Le module de traitement 22 peut être un ordinateur ou tout moyen équivalent apte à recevoir de l'information et à la traiter, par exemple un serveur. Le module de traitement 22 comporte 15 classiquement un processeur programmé pour la mise en oeuvre d'un procédé selon invention, une mémoire électronique, notamment pour stocker les données reçues du thermomètre 16 et des capteurs 18' et une interface de communication permettant notamment de communiquer avec le thermomètre 16 et les capteurs 18,. Le module de traitement est préférentiellement placé à distance du bâtiment 10, par exemple à plus 100 m du bâtiment 10. 20 Le module de traitement 22 est en communication avec le thermomètre 16, qui lui fournit des mesures de la température intérieure Ti(t), et avec les capteurs 18' qui lui fournissent, directement ou par l'intermédiaire du centre météorologique, des mesures Me) sur les facteurs d'influence, respectivement. La communication peut être unidirectionnelle. Elle peut être également bidirectionnelle, 25 notamment avec la chaudière 12, afin que le module de traitement puisse contrôler la chaudière 12. La communication peut être filaire ou hertzienne, de préférence hertzienne. Dans un mode de réalisation préféré, un questionnaire 26 est également utilisé pour interroger les occupants, notamment pour évaluer les caractéristiques thermiques du bâtiment et/ou les 30 facteurs d'usage. Les réponses sont entrées dans le module de traitement. 3038716 8 Procédé Le module de traitement est de préférence programmé pour mettre en oeuvre toutes les étapes d'un procédé selon l'invention et fournir un diagnostic et/ou des préconisations 28. A l'étape a), on choisit un modèle thermique initial qui, de manière générique, définit une 5 relation entre une température intérieure d'un bâtiment et des facteurs d'influence, sans que ce modèle corresponde à un bâtiment spécifique. Dans un mode de réalisation, le modèle thermique initial est basé sur la loi de Fourier de la conduction thermique. Le modèle thermique initial peut être linéaire ou non. Compte tenu de l'inertie thermique du bâtiment, les facteurs d'influence agissent avec retard 10 sur la température intérieure. Pour déterminer la température intérieure à un instant déterminé, un modèle thermique utilise donc classiquement des données mesurées antérieurement audit instant. De préférence, pour déterminer la température intérieure à un instant tl, le modèle thermique initial tient compte des valeurs des facteurs d'influence à plusieurs instants dans une période 15 d'analyse antérieure à l'instant tl. Le nombre de facteurs d'influence pris en compte par le modèle thermique initial peut être faible, par exemple inférieur à 5, pour accélérer le procédé, ou élevé, par exemple supérieur à 10, pour améliorer la précision du modèle thermique obtenu. Dans un mode de réalisation préféré, le modèle thermique initial est choisi dans une base de 20 modèles types, un modèle type correspondant de préférence à un type de bâtiment, c'est-à-dire à un bâtiment présentant des caractéristiques thermiques particulières. Par exemple, un modèle type peut être un modèle qui permet d'approximer grossièrement le comportement thermique d'une maison de moins de dix ans, à un étage, d'une surface comprise entre 150 et 200 m2. 25 Le choix du modèle type peut être guidé par les réponses apportées au questionnaire 26. A titre d'exemple, le modèle initial suivant, particulièrement simple, peut être utilisé : kl *Pchauffe k2*Psoieii = k3* (Ti - Te) + k4*dTi/dt dans lequel : - Pchauffe désigne la puissance de chauffage, 30 - Psoieji désigne la puissance dégagée par l'ensoleillement, - - Te désigne la température extérieure, - Ti désigne la température intérieure, 3038716 9 - t désigne le temps, - k, désigne le paramètre i. Les paramètres k, dépendent des caractéristiques thermiques du bâtiment. Parfois cette dépendance est connue. Par exemple, k3 est égal à environ ((1/R) + 1/Ra), R désignant le 5 coefficient de résistance thermique du bâtiment et Rire, désignant le coefficient lié aux infiltrations par les menuiseries, et k4 est égal à environ C, C désignant le coefficient caractéristique de l'inertie thermique du bâtiment. Les paramètres peuvent également être des matrices à plusieurs dimensions. k2 peut être par exemple une matrice avec les valeurs différentes suivant l'orientation et la hauteur du soleil. 10 Des mesures sur le bâtiment ou sur son environnement stable (route, bâtiment adjacent,...) et/ou les réponses au questionnaire permettent ainsi une première évaluation des valeurs de ces paramètres. Plus précisément, la détermination des valeurs réelles des caractéristiques thermiques du bâtiment peut résulter de mesures. Ces mesures peuvent cependant être difficiles à réaliser, 15 gênantes pour les occupants et/ou coûteuses. Alternativement, ou en complément, il est possible de transmettre un questionnaire 26 aux occupants afin qu'ils transmettent des informations permettant de déterminer les valeurs réelles des caractéristiques thermiques du bâtiment. La compréhension des questionnaires, de nature technique, peut être cependant difficile. Les 20 questionnaires peuvent être également partiellement ou mal remplis. Dans un mode de réalisation préféré, un ou plusieurs cycles c)-d) sont effectués avec - les valeurs des caractéristiques thermiques du bâtiment recueillies avec le questionnaire et/ou mesurées, et - les valeurs réelles de moins de 5, moins de 3, de préférence d'un unique facteur 25 d'influence, par exemple l'énergie apportée par chauffage, mesurée pendant une période d'analyse de préférence inférieure à un mois, 2 semaines, voire 1 semaine. Ces premiers cycles permettent avantageusement d'obtenir un modèle thermique grossier, et de générer un premier diagnostic thermique. Le modèle grossier permet également d'évaluer les caractéristiques thermiques du bâtiment 30 sans mesure ni questionnaire. De préférence, un deuxième questionnaire, identique ou différent du premier questionnaire et présentant tout ou partie du premier diagnostic thermique est fourni aux occupants du 3038716 10 bâtiment. Avantageusement, les inventeurs ont constaté que le deuxième questionnaire est mieux accepté et, de manière générale, mieux rempli que le premier questionnaire. Le deuxième questionnaire permet avantageusement d'affiner le modèle thermique. Tous les paramètres peuvent cependant être évalués par l'étude de l'impact des facteurs 5 d'influence, comme décrit ci-après. Le modèle initial susmentionné utilise trois facteurs d'influence : l'ensoleillement, la température extérieure et la puissance du chauffage. Des informations relatives aux caractéristiques du bâtiment, par exemple mesurées ou fournies par le questionnaire, permettent d'estimer R, C, voire Rail, et donc d'accélérer la 10 détermination des valeurs des paramètres k3 et k4. A l'étape b), on détermine un ensemble de valeurs réelles de la température intérieure du bâtiment 10 et de facteurs d'influence qui seront utilisées pour l'étape c) d'optimisation. L'étape b) est indépendante de l'étape a), c'est-à-dire peut être réalisée avant, pendant ou après l'étape a). Elle peut être réalisée exclusivement avant l'étape c). Elle peut être 15 poursuivie au moins en partie pendant l'étape c). Les valeurs réelles des différents facteurs d'influence et de la température intérieure sont mesurées par les capteurs 18, pendant une période de mesure d'une durée At de préférence au moins égale ou supérieure à six semaines. La période de mesure At peut correspondre à l'ensemble des valeurs réelles mesurées, appelé « historique ». 20 De préférence, si les caractéristiques thermiques du bâtiment sont modifiées, par exemple du fait du renforcement de l'isolation thermique, il en est de même de la période de mesure afin que, pendant ladite période de mesure, les caractéristiques thermiques du bâtiment soient constantes. La poursuite du procédé conduit avantageusement à une adaptation du modèle thermique aux nouvelles caractéristiques thermiques du bâtiment. 25 La fréquence de mesure pour un capteur, de préférence pour chaque capteur, est de préférence au moins égale à 1 mesure par 30 minutes, en particulier pour les facteurs météorologiques. Elle est de préférence supérieure ou égale à 1 mesure par minute pour les facteurs d'usage. De préférence, la fréquence de mesure est déterminée en fonction d'une évaluation du bruit sur lesdites valeurs réelles. Par exemple, la fréquence de mesure de la température intérieure 30 est d'autant plus élevée que le bruit sur les valeurs réelles de la température intérieure est petit. Avantageusement, les valeurs réelles de la température intérieure ainsi filtrées 3038716 11 permettent une appréciation plus précise de la variation de température et évite le « parasitage » par les petites variations locales. On appelle "période d'analyse" l'intervalle de temps pendant lequel les valeurs réelles de la température intérieure du bâtiment et des facteurs d'influence dudit ensemble sont mesurées. 5 Le début et la fin de la période d'analyse sont déterminés en fonction de l'évolution des facteurs d'influence et de la température intérieure. En particulier, une période d'analyse peut être déterminée pour que, pendant ladite période d'analyse, dite "première période d'analyse", seul un facteur d'influence, dit "actif', varie sensiblement. La description de l'étape d) et la figure 2 fournissent des exemples de périodes d'analyse ainsi 10 déterminées. L'étape c) consiste en une optimisation permettant de transformer le modèle thermique initial en un modèle thermique correspondant au bâtiment 10 particulier pour lequel on recherche un modèle thermique, dit « bâtiment étudié ». L'optimisation consiste à modifier les valeurs des paramètres du modèle thermique initial 15 pour les faire correspondre, le plus possible, à la situation réelle. La situation réelle est définie par les valeurs réelles de la température intérieure et des facteurs d'influence de l'ensemble déterminé à l'étape b). Lesdites valeurs réelles peuvent être préalablement filtrées. Le filtrage comporte une opération d'identification du « bruit » (anomalies à filtrer), puis une opération d'extraction de 20 ce bruit de l'ensemble des valeurs réelles. Par exemple, pendant certaines périodes d'analyse, le modèle peut théoriquement se résumer à une relation, par exemple linéaire, entre l'évolution de la température intérieure et un facteur d'influence, l'action des autres facteurs d'influence étant connue, par exemple constante. L'analyse des valeurs réelles à différents instants devrait confirmer cette relation. Les écarts de l'évolution de la température intérieure 25 mesurée par rapport à l'évolution linéaire de la température intérieure calculée par ladite relation peuvent être considérés comme du bruit et filtrés. L'origine du bruit n'est généralement pas connue. Le bruit peut par exemple résulter de l'action de facteurs d'influence mineurs, non pris en compte dans le modèle thermique initial, par exemple l'état d'ouverture ou de fermeture d'une porte. 30 De préférence, on recherche le filtre dans une base de filtres en testant une pluralité de filtres de ladite base de filtres. 3038716 12 La correspondance entre les valeurs des paramètres du modèle thermique initial et la situation réelle est classiquement évaluée au moyen d'un score. Par exemple, avec le modèle susmentionné, le score peut être égal à P kl *- chauffe + k2*Psoleil - k3* (Ti - Te) - k4*dTi/dt. Ce score, appliqué aux valeurs réelles doit être le plus proche 5 possible de zéro pour que la correspondance soit maximale. Le score peut être aussi, par exemple, la différence entre les températures intérieures théorique, calculée avec le modèle, et réelle. Toutes les techniques d'optimisation connues peuvent être envisagées. Dans un mode de réalisation, l'étape c) comporte les étapes suivantes : 10 cl) définition d'un modèle à tester comme étant le modèle thermique initial ; c2) test du modèle à tester avec lesdites valeurs réelles et détermination, par traitement des valeurs réelles au moyen d'une méthode d'analyse, d'un score représentant la correspondance du modèle à tester avec lesdites valeurs réelles ; c3) en fonction du score, modification des valeurs des paramètres du modèle à 15 tester de manière à améliorer le score et reprise à l'étape c2) ou poursuite à l'étape d). Au début d'un cycle, les valeurs des paramètres et le modèle sont dits « à tester ». Le test consiste à appliquer lesdites valeurs réelles au modèle à tester, de manière à déterminer un score, puis à modifier les valeurs des paramètres pour améliorer ce score. Dans 20 l'exemple ci-dessus, l'amélioration du score consiste à le faire tendre vers zéro. L'évaluation du score permet de décider si le cyclage sur les étapes c2) - c3) doit être poursuivi ou stoppé. Par exemple, il peut être décidé de poursuivre le cyclage si la différence entre les températures intérieures théorique et réelle dépasse un premier seuil, par exemple si la différence est supérieure à 5°C et que le premier seuil est, par exemple, de 2°C. 25 La décision peut également dépendre de scores obtenus lors de cycles précédents et/ou de paramètres aléatoires et/ou du nombre de cycles déjà effectués, par exemple. En particulier, il est possible qu'en dépit de la répétition des cycles, on ne parvienne pas à un score acceptable. La décision peut être alors de quitter le cyclage. Cette décision peut résulter, par exemple, d'un nombre de cycles supérieur à un nombre maximal prédéterminé. 30 Un paramètre aléatoire peut également autoriser la poursuite du procédé à l'étape d), bien que le résultat apparaisse insatisfaisant. 3038716 13 Si la décision est la poursuite du cyclage, les valeurs des paramètres testées sont modifiées, ce qui définit un nouveau modèle à tester. Un nouveau cycle est repris à l'étape c2). La méthode d'analyse utilisée pour modifier les valeurs des paramètres testées peut être une quelconque des méthodes connues : identification, statistiques, méthode des moindres carrés 5 linéaire ou non, filtre de Kalman, filtre particulaire, réseaux de neurones artificiels, problème inverse. Par exemple, la méthode des moindres carrés. Des méthodes d'analyse sont notamment décrites dans « Handbook-Mathematics » de -NBronshtein et al, édité par Springer-Verlag Berlin and Heidelberg GmbH & Co. K, « BeyondKalman-Filter-Particle-Applications » de Branko Ristic et al, publié par Artech Print on 10 Demand (January 31, 2004) ou dans « Inverse ProblemTheory and Methods for Model parameter estimation » de Albert Tarantola, édité par Society for Industrial and Applied Mathematics. A l'étape d), selon la valeur d'une fonction d'évaluation, le traitement est poursuivi, après modification des conditions d'optimisation, par une nouvelle étape c), ou le procédé se 15 poursuit à l'étape e). La fonction d'évaluation permet de décider si le cyclage sur les étapes c) - d) doit être poursuivi ou stoppé. La décision peut dépendre du score atteint. Par exemple, il peut être décidé de poursuivre le cyclage si la différence entre les températures intérieures théorique et réelle dépasse un deuxième seuil, par exemple si la différence est supérieure à 5°C et que le 20 deuxième seuil est, par exemple, de 2°C. La décision peut également dépendre de paramètres aléatoires et/ou du nombre de cycles c) d) déjà effectués, par exemple. En particulier, il est possible qu'en dépit de la répétition des cycles, on ne parvienne pas à un score acceptable. La décision peut être alors de quitter le cyclage. Cette décision peut résulter, par exemple, d'un nombre de cycles supérieur à un 25 nombre maximal prédéterminé. Un paramètre aléatoire peut également autoriser la poursuite du procédé à l'étape e), bien que le résultat apparaisse insatisfaisant. Si la décision est la poursuite du cyclage, une modification peut en particulier porter sur : (a) l'ensemble de valeurs réelles analysé, et/ou 30 (b) le modèle thermique initial, et/ou (c) ledit filtre, et/ou 3038716 14 (d) ladite méthode d'analyse. Une première modification, portant notamment sur le choix du modèle thermique initial et le choix de la méthode d'analyse des paramètres dudit modèle thermique, peut être notamment effectuée après analyse du modèle grossier décrit ci-dessus. 5 La modification à l'étape d) dépend de préférence : - de la valeur et/ou de l'évolution temporelle d'un ou plusieurs paramètres ; et/ou - d'un score représentant la correspondance du modèle initial optimisé à l'étape c) avec lesdites valeurs réelles ; et/ou - de l'ensemble de valeurs réelles utilisé à l'étape c) ; et/ou 10 - du modèle thermique initial. (a) De préférence, l'étape c) est mise en oeuvre successivement avec les ensembles de valeurs réelles de ladite température intérieure et de facteurs d'influence mesurées dans des périodes d'analyse différentes, voire disjointes, c'est-à-dire qui ne se superposent pas, même partiellement. 15 Par exemple, avec le modèle thermique initial fourni ci-dessus à titre d'exemple, une ou de préférence plusieurs premières périodes d'analyse peuvent être définies par des périodes pendant lesquelles la chaudière est éteinte, et donc Pchauffe est nul, et il n'y a pas de soleil, par exemple la nuit. Les facteurs d'influence que sont l'ensoleillement et le chauffage sont ainsi « neutralisés ». La température extérieure est alors le seul facteur d'influence qui agit sur la 20 température intérieure. L'analyse de telles premières périodes permet ainsi d'établir des relations entre les valeurs des paramètres k3 et k4. Une ou de préférence plusieurs premières périodes d'analyse peuvent être définies par des périodes pendant lesquelles la température extérieure est sensiblement constante et il n'y a pas de soleil, par exemple la nuit. Les facteurs d'influence que sont la température extérieure et 25 l'ensoleillement sont « neutralisés ». Le chauffage est alors le seul facteur d'influence qui agit sur la température intérieure. L'analyse de telles premières périodes permet ainsi d'établir des relations entre les valeurs des paramètres kl, k3 et k4. Une ou de préférence plusieurs premières périodes d'analyse peuvent être définies par des périodes pendant lesquelles la chaudière est éteinte, et donc Pchauffe est nul, et la température 30 extérieure est sensiblement constante. Les facteurs d'influence que sont le chauffage et la température extérieure sont « neutralisés ». L'ensoleillement est alors le seul facteur 3038716 15 d'influence qui agit sur la température intérieure. L'analyse de telles premières périodes permet ainsi d'établir des relations entre les valeurs des paramètres k2, k3 et k4. Une ou de préférence plusieurs premières périodes d'analyse peuvent être définies par des périodes pendant lesquelles la température intérieure est sensiblement constante et 5 l'ensoleillement est nul, par exemple parce que les mesures sont effectuées la nuit. Le modèle se résume alors à Pchauffe = (k3/k1)* (Ti - Te) La figure 3 représente un ensemble de mesures effectuées pendant de telles premières périodes. Une régression linéaire permet de déterminer une droite dont la pente est égale à 10 (k3/k1). Le processus peut se poursuivre avec d'autres premières périodes d'analyse. L'ensemble des relations entre les valeurs des paramètres ainsi établies permet de déterminer finement ces valeurs. L'utilisation de dites « premières périodes d'analyse » pour évaluer les valeurs des paramètres 15 est particulièrement efficace. En effet, elle permet de filtrer des périodes pendant lesquelles les actions des différents facteurs d'influence sont complexes et difficiles, voire impossibles à modéliser de manière satisfaisante. Par exemple, si toute la période de mesure était utilisée comme période d'analyse, la période d'analyse comporterait de nombreuses périodes dans lesquelles les valeurs réelles ne permettraient pas d'expliquer l'action combinée des facteurs 20 d'influence actifs, notamment parce que des mesures appropriées n'ont pas été effectuées. Par exemple les mesures ne permettent généralement pas de déterminer si une fenêtre est ouverte, ou avec quel degré d'ouverture. Les « premières périodes d'analyse » sont avantageusement des périodes définies pour éliminer les facteurs d'influence non maîtrisés. De préférence, la durée de toute période d'analyse est supérieure à 1 heure, de préférence 25 supérieure à 3 h, et/ou inférieure à 24 h, de préférence inférieure à 12 h. La durée d'une « première période d'analyse » peut en particulier dépendre du facteur d'influence actif et, de préférence, de son évolution. Par exemple, si le facteur d'influence actif est le vent. La période d'analyse peut être initialement d'une durée de 5 heures. Si la température intérieure varie beaucoup, peut-être à 30 cause de facteurs d'influence non identifiés, la durée est de préférence raccourcie, par exemple à 2 heures, pour pouvoir mieux identifier l'impact spécifique du vent. Les autres facteurs d'influence variant alors peu. 3038716 16 Au contraire si la température intérieure a des variations très lentes, la période d'analyse est de préférence allongée, par exemple à 8 heures, afin de disposer d'un ensemble de valeurs réelles plus grand. Ainsi, la durée d'une période d'analyse, de préférence de toute période d'analyse est de 5 préférence déterminée en fonction d'une analyse statistique des valeurs réelles, et en particulier, de préférence, en fonction de l'évolution temporelle d'un ou plusieurs paramètres. De préférence, pendant une première période d'analyse , tous les facteurs d'influence sauf le facteur d'influence actif, varient dans une plage de variation inférieure à une fraction de leur amplitude maximale de variation pendant une période déterminée, par exemple pendant les 12 10 heures précédant le début de la période d'analyse envisagée, ou pendant la période de mesure. Ladite fraction, c'est-à-dire la tolérance sur la variation, est de préférence inférieure à 20%. Elle est de préférence déterminée en fonction du facteur d'influence considéré. Par exemple, une variation de moins de 20% de la vitesse du vent peut être considérée comme correspondant à une période pendant laquelle la vitesse du vent est sensiblement constante. 15 Mais une variation de plus de 10% de la température extérieure peut au contraire être considérée comme significative (température extérieure non constante). De préférence ladite fraction est déterminée par analyse de l'influence des variations du facteur d'influence considéré sur les valeurs des paramètres du modèle. L'analyse peut être en particulier une analyse statistique des valeurs réelles. 20 Par exemple, le début d'une « première période d'analyse » pour l'ensoleillement (facteur d'influence actif) peut correspondre à un instant to à partir duquel la vitesse du vent varie de moins de 20% et la température extérieure varie de moins de 5% par rapport à leur amplitude maximale respective pendant les 12 heures précédant ledit instant. La fin de cette première période d'analyse peut correspondre au premier instant à partir duquel la vitesse du vent varie 25 de plus de 20% ou la température extérieure varie de plus de 5% par rapport à leur amplitude maximale respective pendant les 12 heures précédant ledit instant to. Comme représenté sur la figure 2, une première période d'analyse peut par exemple correspondre à un moment où l'ensoleillement était nul (nuit), où le vent était nul, où le ciel était dégagé, sans pluie, et où la température extérieure variait sensiblement (période Dl). 30 Une première période d'analyse peut également, par exemple, correspondre à un moment pendant lequel la température extérieure était relativement stable, l'ensoleillement nul et le vent augmentait rapidement (période D2), ou pendant lequel le vent était sensiblement nul, la 3038716 17 température extérieure était stable, le temps était clair, et l'ensoleillement augmentait rapidement (période D3). Pendant les périodes Dl, D2 et D3, les facteurs d'influence actifs étaient donc la température extérieure, la force du vent et l'ensoleillement, respectivement. Pendant une première période d'analyse, l'évolution de la température intérieure peut ainsi 5 être attribuée au facteur d'influence actif. La valeur du ou des paramètres du modèle associés au facteur d'influence actif peut être ainsi facilement estimée. La période d'analyse peut être la même quel que soit le facteur d'influence considéré. Elle peut être identique pour les facteurs d'influence et la température intérieure. Cependant, un facteur d'influence actif n'agit pas immédiatement sur la température intérieure du fait de l'inertie 10 thermique du bâtiment. De préférence, la période d'analyse pour un facteur d'influence est donc décalée temporellement, en fonction d'une évaluation de ladite inertie thermique, par rapport à la période d'analyse pour la température intérieure, la période d'analyse pour le facteur d'influence commençant avant la période d'analyse pour la température intérieure. De préférence, l'évaluation de l'inertie thermique, et donc ledit décalage temporel sont affinés à 15 mesure que le modèle thermique est optimisé. Le décalage temporel est de préférence déterminé en fonction du facteur d'influence considéré. Par exemple, le décalage pour le chauffage sera plus faible que le décalage pour la vitesse du vent, le chauffage agissant plus rapidement que la vitesse du vent sur la température intérieure. 20 De préférence, le décalage est déterminé par analyse de l'influence des variations du facteur d'influence considéré sur la température intérieure. L'analyse peut être en particulier une analyse statistique des valeurs réelles. L'étape d) permet de changer la période d'analyse de manière à rendre actifs différents facteurs d'influence lors des exécutions successives de l'étape c). De préférence, cette 25 succession d'exécutions de l'étape c) permet, à l'issue de l'étape d), d'avoir analysé au moins 2, de préférence au moins 3, de préférence au moins 4 facteurs d'influence différents, de préférence tous les facteurs d'influence. Dans un mode de réalisation préféré, les exécutions successives de l'étape c) se font d'abord à partir de « premières périodes d'analyse ». De préférence, au moins une première période 30 d'analyse est déterminée pour chacun des facteurs d'influence, et en particulier pour chacun des facteurs d'influence météorologiques. L'optimisation de l'étape c), répétée à partir des ensembles de valeurs réelles de ladite température intérieure et desdits facteurs d'influence de 3038716 18 chacune de ces premières périodes d'analyse, permet avantageusement d'obtenir rapidement un modèle thermique de qualité correcte. Certains facteurs d'influence peuvent être corrélés, de sorte qu'il n'existe pas de première période d'analyse dans laquelle un seul de ces deux facteurs d'influence évolue. 5 De préférence, le cycle des étapes c) et d) est poursuivi en déterminant, à l'étape d), des périodes d'analyse définies de manière que, pendant une période d'analyse, dite « deuxième période d'analyse », tous les facteurs d'influence sauf deux, dits « facteurs d'influence actifs », soient sensiblement constants. L'optimisation de l'étape c), répétée à partir des ensembles de valeurs correspondant permet d'améliorer encore la qualité du modèle 10 thermique. Le cycle des étapes c) et d) est de préférence ainsi poursuivi en déterminant, à l'étape d), des périodes d'analyse comportant un nombre croissant de facteurs d'influence actifs, de préférence en déterminant des troisièmes, quatrièmes, cinquièmes, etc. périodes d'analyse dans lesquelles tous les facteurs d'influence sauf trois, quatre, cinq etc. sont sensiblement 15 constants. De préférence, le nombre de facteurs d'influence actifs augmente linéairement, de préférence d'une unité à chaque changement. Un changement du nombre de facteurs d'influence actifs peut être décidé en particulier lorsque le cyclage des étapes c) et d) ne permet plus d'améliorer les résultats, c'est-à-dire de réduire l'écart entre les températures intérieures 20 théorique et réelle. La modification des périodes d'analyse à l'étape d) contribue considérablement à l'obtention rapide d'un modèle thermique précis. La modification de l'ensemble de valeurs réelles peut consister en une modification de la 25 période d'analyse, mais aussi, alternativement ou en complément, en une modification de la fréquence de mesure (échantillonnage). (b) A l'étape d), le modèle thermique initial peut être également modifié. Un changement de modèle thermique initial est de préférence décidé lorsque les résultats obtenus avec un premier modèle thermique initial sont estimés insuffisants et que cette 30 insuffisance est potentiellement attribuable audit modèle thermique initial. 3038716 19 De préférence, un changement de modèle thermique initial est décidé lorsqu'une ou plusieurs occurrences de l'étape c) ont abouti à des valeurs de paramètres typiques d'un autre modèle thermique initial. Par exemple, le procédé a pu conduire à des valeurs de paramètres correspondant à certaines caractéristiques thermiques caractéristiques d'un bâtiment type, par 5 exemple d'un bâtiment présentant un vide-sanitaire. On remplace alors le modèle thermique initial par le modèle thermique initial de ce bâtiment type. De préférence, on réutilise cependant les valeurs de paramètres déterminées précédemment pour le nouveau modèle thermique initial. (c) A l'étape d), le filtre peut être également modifié. 10 De préférence, le filtre peut être modifié afin de rechercher le filtre le plus adapté. De préférence, on recherche le filtre dans une base de filtres en testant une pluralité de filtres de ladite base de filtres. Le filtrage par rapport à des valeurs théoriques fournies par un modèle peut cependant conduire à filtrer de manière inappropriée les valeurs réelles. En particulier, les écarts peuvent 15 résulter d'imperfections dans le modèle. Par exemple, pendant certaines périodes d'analyse, le modèle peut théoriquement se résumer à une relation linéaire, entre l'évolution de la température intérieure et un facteur d'influence, mais cette relation susmentionnée peut ne pas être exactement linéaire dans la réalité. La recherche du filtre est de préférence renouvelée lorsque la période d'analyse à laquelle 20 appartiennent les valeurs réelles analysées n'a jamais été utilisée. Par exemple, un filtre peut être déterminé au regard de valeurs réelles d'une période d'analyse afin d'améliorer le modèle lors de la prise en compte de ces valeurs réelles. Cependant, si l'application de ce filtre à des valeurs réelles d'autres périodes d'analyse similaires conduit à dégrader le modèle, le filtre peut être modifié, voire supprimé. De préférence, le filtre appliqué dépend ainsi de son 25 efficacité à améliorer les valeurs des paramètres. La modification du filtre pour rechercher le filtre le plus adapté conduit avantageusement à un filtrage qui agit différemment en fonction de l'objet auquel il est appliqué, c'est-à-dire selon qu'il s'applique à la température intérieure ou à un facteur d'influence, et différent selon ledit facteur d'influence. 30 La modification du filtre peut également dépendre de la valeur d'un ou plusieurs paramètres. 3038716 20 Par exemple, le bruit sur la mesure de température intérieure correspond à des petites variations non corrélées à l'évolution générale. Ces petites variations peuvent résulter, par exemple, d'une ouverture de fenêtre, de porte, ou d'une personne qui touche le capteur de température. Une première évaluation des valeurs des paramètres peut amener à conclure que 5 l'inertie thermique du bâtiment est élevée. Le filtre peut alors être choisi pour filtrer particulièrement la température intérieure, ce qui permet d'éliminer efficacement ce bruit. (d) A l'étape d), la méthode d'analyse des valeurs réelles peut être également modifiée. La méthode d'optimisation du modèle est alors évolutive. Notamment, l'invention permet des ruptures pour rechercher une méthode d'optimisation plus adaptée et permettant d'aboutir à 10 un modèle thermique plus proche de la réalité. Les changements de méthode d'analyse pour évaluer les modifications à apporter aux valeurs des paramètres du modèle, à l'étape c), peuvent être réalisés à intervalles réguliers ou irréguliers, voire de manière aléatoire. Dans un mode de réalisation, lesdits changements de méthode d'analyse dépendent des 15 progrès effectués entre deux cycles successifs. En particulier, un changement est avantageux lorsque l'écart entre les températures théorique et réelle ne parvient plus à être réduit par modification des paramètres du modèle. De préférence, la méthode d'analyse mise en oeuvre dépend de l'ensemble de valeurs réelles analysé, c'est-à-dire de la période d'analyse. De préférence, la méthode d'analyse dépend de 20 la nature et/ou du nombre de dits « facteurs d'influence actifs », c'est-à-dire des facteurs d'influence qui varient sensiblement pendant la période d'analyse. Dans un mode de réalisation, des changements de la méthode d'analyse sont donc effectués suite à une modification d'une période d'analyse. Par exemple, si un seul facteur d'influence Fi est variable pendant la période d'analyse, 25 l'ensemble des valeurs réelles mesurées sera constitué de couples de valeurs réelles (Fi, Ti), mesurés à différents instants pendant la période d'analyse. Une analyse du nuage de points correspondant par la méthode des moindres carrés permettra d'établir une valeur pour un paramètre reliant Fi et Ti. Si seuls deux facteurs d'influence sont variables pendant une période d'analyse, la méthode d'analyse sera par exemple la méthode du problème inverse. 3038716 21 De préférence, la méthode d'analyse mise en oeuvre dépend des caractéristiques thermiques du bâtiment et/ou du modèle thermique initial, et en particulier des facteurs d'influence pris en compte dans le modèle thermique initial. Par exemple, lors d'une première occurrence de l'étape c), la méthode d'analyse dite « des 5 moindres carrés » peut être utilisée pour trouver une valeur d'un paramètre mettant en relation deux facteurs d'influence. Lors d'une deuxième occurrence, la méthode d'analyse dite « du problème inverse » peut être mise en oeuvre. Les méthodes mises en oeuvre pour modifier les valeurs des paramètres sont de préférence choisies parmi les méthodes suivantes : identification, statistiques, méthode des moindres 10 carrés linéaire ou non, filtre de Kalman, filtre particulaire, réseaux de neurones artificiels et problème inverse. De préférence, la méthode d'analyse mise en oeuvre dépend de l'ensemble de valeurs réelles déterminé à l'étape b) et/ou du modèle thermique initial et/ou du score permettant d'évaluer l'efficacité de l'optimisation. 15 A l'étape e), le modèle est utilisé pour analyser le comportement du bâtiment en réaction à des évolutions de facteurs d'influence et/ou de caractéristiques thermiques du bâtiment, par exemple résultant d'une modification de l'isolation, et/ou pour évaluer lesdites caractéristiques thermiques. L'étape e) peut être réalisée à tout instant après l'analyse d'au moins un dit ensemble de 20 valeurs réelles. En particulier, elle peut être réalisée alors que l'optimisation est encore en cours. Le modèle thermique utilisé est alors le meilleur modèle obtenu jusque-là. De préférence, un rapport exposant les résultats de l'étape e) est généré. Un « moteur de règles » peut être mis en oeuvre à cet effet. Le rapport peut comporter un diagnostic exposant 25 - l'environnement du bâtiment, en particulier l'environnement urbain ou l'exposition au vent ou au soleil ; - une évaluation de la résistance thermique globale et/ou de l'efficacité des menuiseries et/ou des volets et/ou du chauffage et/ou de la ventilation et/ou des apports de chaleur interne (cuisson, appareils ménagers, consommations de veille des appareils 30 électroniques, par exemple) et/ou du thermostat. Le rapport peut comporter des préconisations, notamment pour : 3038716 22 - isoler ou renforcer l'isolation en des emplacements déterminés du bâtiment, par exemple changer les volets ; - modifier le chauffage, par exemple en préconisant des réglages de la puissance de la chaudière et/ou des durées et/ou des moments du chauffage ; 5 - préconiser des changements de comportement, par exemple sur la ventilation et/ou l'ouverture des volets et/ou l'aération naturelle, et/ou des modifications dans l'emplacement du thermostat ; - préconiser l'installation de capteurs et/ou de régulateurs pour affiner le modèle ; - préconiser l'installation de systèmes de contrôle, notamment de contrôle automatique, 10 pour gérer plus efficacement les consommateurs d'énergie et/ou les moyens isolant, comme les volets ; - fournir des informations facilitant la mise en oeuvre des préconisations, comme des adresses de prestataire et/ou des fiches techniques et/ou des coûts à prévoir ; - préconiser des comportements en cohérence avec les variations des facteurs 15 d'influence, par exemple adapter sa tenue vestimentaire à la température extérieure. Dans un mode de réalisation, le module de traitement peut commander la chaudière (commande C sur la figure 1) et/ou les volets et/ou la VMC en conséquence. En fonction de la complexité du modèle, les préconisations peuvent être plus ou moins précises. Par exemple, si le modèle utilise la direction du vent comme facteur d'influence, il 20 est possible de préconiser des travaux d'isolation spécifiquement sur une façade déterminée. Comme cela apparaît clairement à présent, l'invention fournit un procédé permettant d'obtenir un modèle thermique précis et faible du comportement thermique d'un bâtiment. En particulier, les modifications à l'étape d), en particulier les modifications des périodes d'analyse et les changements de méthode d'analyse pour évaluer les modifications à apporter 25 aux valeurs des paramètres du modèle, autorisent des « ruptures » dans l'optimisation. Ces ruptures permettent avantageusement d'améliorer le modèle thermique. Il en résulte un diagnostic fiable et précis, et, de préférence, des préconisations bien adaptées à ce bâtiment. Bien entendu la présente invention n'est pas limitée aux modes de réalisation décrits, fournis 30 à titre d'exemples illustratifs et non limitatifs. 3038716 23 En particulier, dans un mode de réalisation, l'optimisation n'est jamais arrêtée, de sorte que le modèle thermique s'améliore en permanence et/ou s'adapte à toute évolution du bâtiment ou de son environnement, par exemple du fait d'un changement dans l'isolation thermique du bâtiment. REVENDICATIONSI. Procédé de détermination d'un diagnostic thermique par la mesure d'une température intérieure (Ti) d'un bâtiment (10), le procédé comportant les étapes successives suivantes : a) choix d'un modèle thermique initial établissant, au moyen d'un ensemble de paramètres, une relation théorique entre ladite température intérieure et un ensemble de facteurs d'influence (Fi) de ladite température intérieure ; b) détermination d'un ensemble de valeurs réelles de ladite température intérieure et desdits facteurs d'influence ; c) optimisation des paramètres du modèle thermique initial par analyse, au moyen d'une méthode d'analyse, dudit ensemble de valeurs réelles, optionnellement filtrées au moyen d'un filtre ; d) en fonction d'une fonction d'évaluation, modification dudit ensemble de valeurs réelles et/ou dudit modèle thermique initial et/ou dudit filtre et/ou de ladite méthode d'analyse, et reprise à l'étape e) et/ou poursuite à l'étape e) ; e) utilisation du modèle thermique obtenu pour réaliser un diagnostic thermique. 2. Procédé selon la revendication précédente, dans lequel la modification à l'étape d) dépend - de la valeur et/ou de l'évolution temporelle de la valeur d'un ou plusieurs paramètres ; et/ou - d'un score représentant la correspondance des valeurs des paramètres du modèle initial optimisé à l'étape c) avec lesdites valeurs réelles ; et/ou - de l'ensemble de valeurs réelles utilisé à l'étape c) ; et/ou - du modèle thermique initial. 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel, à l'étape a), le modèle thermique initial est choisi, en fonction d'un ou plusieurs facteurs d'influence propres audit bâtiment, dans une base de modèles types. 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel, à l'étape b) et/ou à l'étape d), les valeurs réelles de la température intérieure du bâtiment et des facteurs d'influence dudit ensemble sont mesurées dans une période d'analyse dont le 3038716 25 début et la fin sont déterminés en fonction de l'évolution temporelle des facteurs d'influence. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel, à l'étape b) et/ou à l'étape d), les valeurs réelles de la température intérieure du bâtiment et de 5 facteurs d'influence dudit ensemble sont mesurées dans une période d'analyse dont la durée est supérieure à 1 heure et inférieure à 24 heures. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel, à l'étape b) et/ou à l'étape d), les valeurs réelles de la température intérieure du bâtiment et de facteurs d'influence dudit ensemble sont mesurées dans une période d'analyse pendant 10 laquelle tous les facteurs d'influence sauf un, dit « facteur d'influence actif », sont sensiblement constants. 7. Procédé selon la revendication immédiatement précédente, dans lequel l'étape c) est exécutée plusieurs fois, chaque exécution étant réalisée à partir d'un ensemble de valeurs réelles de la température intérieure du bâtiment et des facteurs d'influence mesurées dans 15 une période d'analyse pendant laquelle tous les facteurs d'influence sauf un, dit « facteur d'influence actif », sont sensiblement constants pendant les douze heures précédant le début de ladite période d'analyse, les facteurs d'influence actifs ayant été successivement analysés du fait de l'exécution répétée de l'étape c) comportant au moins deux facteurs d'influence différents. 20 8. Procédé selon la revendication immédiatement précédente, les facteurs d'influence actifs ayant été successivement analysés du fait de l'exécution répétée de l'étape c) comportant au moins deux facteurs d'influence différents parmi la température extérieure, la vitesse du vent, l'orientation du vent, l'ensoleillement, la pression atmosphérique, l'humidité, la présence d'occupants dans le bâtiment, la puissance de chauffage et l'état de 25 fonctionnement du système de chauffage. 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel, à l'étape d), la modification de la méthode d'analyse est un choix dans le groupe constitué par l'identification, les statistiques, la méthode des moindres carrés, linéaire ou non, le filtre de Kalman, le filtre particulaire, les réseaux de neurones artificiels et le problème inverse. 30 10. Procédé selon l'une quelconque des revendications précédentes, dans lequel, à l'étape e), on génère un rapport comportant un diagnostic exposant 3038716 26 - l'environnement du bâtiment, en particulier l'environnement urbain ou l'exposition au vent ou au soleil ; et/ou - une évaluation de la résistance thermique globale et/ou de l'efficacité des menuiseries et/ou des volets et/ou du chauffage et/ou de la ventilation et/ou des apports de chaleur 5 internes et/ou du thermostat ; et/ou des préconisations pour : - isoler ou renforcer l'isolation en des emplacements déterminés du bâtiment; et/ou - modifier le chauffage; et/ou - préconiser des changements de comportement et/ou des modifications dans 10 l'emplacement du thermostat ; et/ou - préconiser l'installation de capteurs et/ou de régulateurs pour affiner le modèle ; et/ou - préconiser l'installation de systèmes de contrôle, pour gérer plus efficacement les consommateurs d'énergie et/ou les moyens isolants; et/ou - fournir des informations facilitant la mise en oeuvre des préconisations. 15 11. Dispositif de détermination d'un diagnostic thermique par la mesure d'une température intérieure d'un bâtiment, ledit dispositif comportant : un thermomètre (16) disposé pour mesurer, à tout instant t, une température intérieure Ti(t) à l'intérieur d'un bâtiment (10), - un ensemble de capteurs (18i) destinés chacun à mesurer, à tout instant t, un facteur 20 d'influence de ladite température intérieure, - un module de traitement (22) comportant un programme d'ordinateur comprenant des instructions de code de programme pour l'exécution d'un procédé selon l'une quelconque des revendications précédentes, lorsque ledit programme est exécuté par un ordinateur. 25