FR 3062723 A1 20180810 FR 1750903 A 20170203 1 PROCEDE DE DETERMINATION D'UNE GRANDEUR THERMIQUE D'UNE ZONE D'UN BATIMENT Domaine de l'invention L'invention se situe dans le domaine de la caractérisation thermique d'un bâtiment et concerne un procédé, un dispositif et un produit programme d'ordinateur pour la détermination d'une grandeur thermique d'une zone d'un bâtiment. Etat de la Technique Dans l'univers du bâtiment, il existe deux types d'approche de caractérisation des performances thermiques d'un bâtiment. Le premier type d'approche est basé sur les méthodes dites inverses et vise à déterminer les performances globales du bâtiment, par exemple la résistance thermique et/ou la capacité thermique du bâtiment. Cette approche a pour principe d'utiliser des mesures ou données générées à partir de logiciels de simulation pour identifier les paramètres thermiques du bâtiment à l'aide d'algorithmes dédiés. Elle s'appuie sur un modèle numérique du bâtiment, c'est-à-dire une description de l'enveloppe du bâtiment et de ses éventuels systèmes de chauffage et/ou refroidissement. Cette approche a pour objectif de fiabiliser les prédictions réalisées à l'aide de modèles (en particulier pour une utilisation dans le cadre de stratégies de contrôle-commande). Cette approche vise aussi à déterminer les grandeurs thermiques de l'enveloppe du bâtiment dans sa globalité et non pas localement. Le second type d'approche est une approche de détermination des performances thermiques de parois à l'aide d'une instrumentation dédiée. Cette approche vise à caractériser finement les propriétés thermophysiques de parois et fait appel à un protocole expérimental précis, comprenant dans certains cas des dispositifs de chauffage supplémentaires 3062723 2 pour créer artificiellement une différence entre les températures de surfaces de la paroi, et à des capteurs permettant de mesurer des variables liées aux grandeurs physiques locales de la paroi (par exemple températures de surface, flux de chaleur). Cette approche est fiable mais elle présente le 5 défaut de nécessiter une instrumentation avancée et intrusive. En effet, cette approche nécessite d'utiliser des capteurs d'utilisation peu courante (comme les thermocouples, les fluxmètres, ...) et de les mettre en oeuvre selon un protocole expérimental contraignant. De plus, pour certaines configurations de bâtiment, il peut être difficile voire impossible de déployer ces capteurs, 10 notamment dans le cas d'une mesure en étage élevé ou si le revêtement de façade rend la pose des capteurs difficile, comme c'est le cas sur un mur en briques. Il en ressort que lorsque l'objectif de l'approche est de fiabiliser les 15 modèles en vue de réaliser des prédictions, la vérification de la validité du modèle est réalisée sur la base d'une comparaison entre les prédictions réalisées à l'aide du modèle calibré et des mesures. La validité des paramètres thermiques de l'enveloppe en tant que tels ne fait pas partie des objectifs visés et n'est donc pas assurée. 20 Lorsque l'objectif est de caractériser les performances thermiques du bâtiment, la principale limite est liée au fait que ces méthodes se focalisent sur la caractérisation du bâtiment dans son entier. Elles nécessitent donc d'avoir la connaissance de la géométrie complète du 25 bâtiment, information qu'il est souvent difficile d'obtenir sur les bâtiments existants. Par ailleurs, ces méthodes ne permettent de calculer que des coefficients globaux (par exemple le coefficient de déperdition global). Ce type de méthode requiert généralement en outre de collecter, en plus de la température, des mesures supplémentaires (telles que le rayonnement et/ou l'humidité). Certaines de ces approches peuvent aussi nécessiter d'avoir une connaissance préalable des coefficients globaux de l'enveloppe du bâtiment, et se cantonnent à une répartition de cette valeur globale sur les couches internes modélisées de la paroi. 3062723 3 Résumé de l'invention L'invention vise à pallier tout ou partie des problèmes cités plus haut en proposant un procédé de détermination d'une valeur d'une grandeur thermique d'une zone d'un bâtiment, telle que la résistance thermique et/ou 5 la capacité thermique des parois et ouvrants, dans un temps limité et de façon fiable à partir de mesures simples, peu coûteuses et non intrusives. Il en résulte une facilité et robustesse d'utilisation du procédé car l'utilisateur peut déterminer la valeur d'une grandeur thermique d'une paroi ou d'un ouvrant d'une zone précise, c'est-à-dire localement, grâce à une 10 instrumentation simple et sans connaissances particulières du bâtiment. A cet effet, l'invention a pour objet un procédé, mis en oeuvre par ordinateur, de détermination d'une valeur d'au moins une grandeur thermique d'une zone d'un bâtiment soumise à au moins une donnée 15 thermique environnementale et ayant un comportement thermique réel, comprenant: - Une étape de construction d'un modèle thermique initial de la zone ayant un comportement thermique théorique, le modèle thermique initial utilisant la valeur de la au moins une grandeur thermique, la 20 valeur de la au moins une grandeur thermique étant initialisée à une valeur initiale, - Une étape d'entrée dans le modèle thermique initial de la zone d'au moins une donnée thermique environnementale issue d'une mesure, - Une étape d'exécution d'un algorithme d'identification au modèle 25 thermique initial modifiant la valeur initiale de la au moins une grandeur thermique jusqu'à obtenir un modèle thermique calibré ayant un comportement thermique sensiblement égal au comportement thermique réel, - Une étape d'extraction de la valeur de la au moins une grandeur 30 thermique du modèle thermique calibré. Selon un mode de réalisation, le procédé de détermination selon l'invention comprend préalablement à l'étape de construction du modèle thermique initial de la zone une étape de relevé des caractéristiques 35 géométriques de la zone du bâtiment, et l'étape de construction du modèle 3062723 4 thermique initial est basée sur le relevé des caractéristiques géométriques de la zone du bâtiment. Avantageusement, le procédé selon l'invention comprend après 5 l'étape d'extraction de la valeur de la au moins une grandeur thermique du modèle thermique calibré une étape d'affichage de la valeur de la au moins une grandeur thermique. Avantageusement, l'étape de construction du modèle thermique 10 initial est réalisée par analogie électrique. Avantageusement, le procédé selon l'invention comprend, préalablement à l'étape d'entrée de la au moins une donnée thermique environnementale issue d'une mesure, une étape de mesure de la au moins 15 une donnée thermique environnementale sur une période prédéterminée pour obtenir la au moins une donnée thermique environnementale issue d'une mesure. Avantageusement, l'étape de mesure de la au moins une donnée 20 thermique environnementale sur une période prédéterminée comprend la mesure d'au moins une température intérieure de la zone et une température extérieure à la zone. Avantageusement, l'étape de mesure de la au moins une donnée 25 thermique environnementale sur une période prédéterminée comprend la mesure du rayonnement global du bâtiment. Avantageusement, l'étape d'exécution d'un algorithme consiste à exécuter un algorithme de recherche d'une solution optimale avec une 30 approche de régression. Avantageusement, l'algorithme d'identification est le krigeage, ou un algorithme Particle Swarm Optimization (PSO) ou un algorithme Support Vector Machines (SVM) ou une combinaison d'un algorithme Particle Swarm 35 Optimization (PSO) et Support Vector Machines (SVM). 3062723 5 Selon un mode de réalisation, la au moins une grandeur thermique est la résistance thermique d'une paroi de la zone et/ou la capacité thermique de la zone. 5 L'invention concerne aussi un dispositif de détermination d'une valeur d'au moins une grandeur thermique d'une zone d'un bâtiment soumise à au moins une donnée thermique environnementale et ayant un comportement thermique réel, le dispositif comprenant des moyens pour mettre en oeuvre les étapes du procédé selon l'invention. 10 L'invention concerne aussi un produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes du procédé de détermination selon l'invention, lorsque ledit programme est exécuté sur un ordinateur. 15 Plusieurs applications peuvent bénéficier des avantages apportés par l'invention. Par exemple, l'invention peut servir de support à un processus de certification thermique de bâtiment, dans laquelle une sélection de zones thermiques peut être caractérisée afin de vérifier la conformité du 20 bâtiment à un standard ou un label donné. La rapidité et la simplicité de l'invention permettent en effet d'envisager, pour un coût limité, son application simultanée à plusieurs zones du bâtiment. Cette invention s'applique à tout type de bâtiment, sans distinction 25 d'usage, c'est-à-dire qu'il s'agisse d'un bâtiment résidentiel, tertiaire, collectif, industriel ou d'une maison individuelle. L'invention peut s'appliquer à des bâtiments neufs, nouvellement construits, et à des bâtiments existants, sans distinction d'époque. 30 L'invention concerne aussi un dispositif de détermination d'une valeur d'au moins une grandeur thermique d'une zone d'un bâtiment soumise à au moins une donnée thermique environnementale et ayant un comportement thermique réel, le dispositif comprenant des moyens pour mettre en oeuvre les étapes d'un tel procédé. 35 3062723 6 L'invention concerne aussi un produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes du procédé de construction, lorsque ledit programme est exécuté sur un ordinateur. 5 Description des figures L'invention sera mieux comprise et d'autres avantages apparaîtront à la lecture de la description détaillée d'un mode de réalisation 10 donné à titre d'exemple, description illustrée par le dessin joint dans lequel : - la figure 1 représente schématiquement le principe du procédé de détermination d'une valeur d'au moins une grandeur thermique de deux zones d'un bâtiment selon l'invention, 15 - la figure 2 représente schématiquement les étapes d'un mode de réalisation d'un procédé de détermination d'une valeur d'au moins une grandeur thermique d'une zone d'un bâtiment selon l'invention, - la figure 3 représente schématiquement un exemple d'un modèle thermique initial utilisant les grandeurs thermiques à déterminer 20 selon l'invention, - la figure 4 représente un exemple de zone à laquelle le procédé selon l'invention a été appliqué pour en déterminer la valeur des grandeurs thermiques. 25 Description détaillée de l'invention Par souci de clarté, les mêmes éléments porteront les mêmes repères dans les différentes figures. 30 Dans la description, l'invention est décrite pour une (ou plusieurs) zone(s). Le terme zone est à comprendre comme une zone thermique qui est un espace d'un bâtiment supposé thermiquement homogène. Une zone thermique peut consister en une seule pièce ou en plusieurs pièces adjacentes. Une zone thermique est ceinte de parois (sur les côtés, sur le 3062 723 7 haut et sur le bas) de même nature ou non, pouvant comprendre des ouvrants (porte(s) et/ou fenêtre(s)). Chaque paroi peut être opaque ou non, par exemple un mur traditionnel ou une paroi vitrée. La zone peut comprendre ou non des parois donnant sur l'extérieur du bâtiment. Elle peut 5 s'étendre horizontalement sur une ou plusieurs pièces adjacentes et/ou verticalement sur plusieurs étages du bâtiment. La configuration la plus courante, mais non limitative, pour l'application de l'invention est celle dans laquelle une zone thermique 10 correspond à une pièce du bâtiment dont une paroi comprend généralement un ou des ouvrant(s) et donne sur l'extérieur du bâtiment, les autres parois de la zone thermique donnant sur des pièces adjacentes. L'intérêt de l'invention réside dans la possibilité d'évaluer de façon 15 fiable notamment la résistance thermique réelle des parois et ouvrants d'une zone donnée d'un bâtiment existant et de répéter cette évaluation, de façon simultanée ou en séquence sur autant de zones que nécessaire. La figure 1 représente schématiquement le principe du procédé 20 de détermination d'une valeur d'au moins une grandeur thermique de deux zones 11, 12 d'un bâtiment 10 selon l'invention. La détermination d'une valeur d'au moins une grandeur thermique d'une zone 11, 12 du bâtiment 10 selon l'invention est représentée symboliquement par la référence 20. Chaque zone 11, 12 du bâtiment 10 est soumise à au moins une donnée 25 thermique environnementale 21, 22. Cette au moins une donnée thermique environnementale peut être la température extérieure, la température des zones adjacentes, le rayonnement solaire appliqué à la zone, etc. Chaque zone 11, 12 du bâtiment a un comportement thermique réel. Ce comportement thermique réel se traduit notamment par l'évolution de la 30 température intérieure de la zone 11, 12 en fonction de la ou des données thermiques d'environnement auxquelles elle est soumise. Le procédé selon l'invention est basé sur l'utilisation d'algorithmes d'identification capables de calculer la valeur des grandeurs thermiques des 35 parois de la zone 11, 12 à partir de données géométriques de la zone 3062723 8 considérée (longueur et hauteur des parois et ouvrants) et des mesures 21, 22 de températures d'ambiance (c'est-à-dire température intérieure de la zone) et extérieure et rayonnement global réalisées sur une période prédéterminée par exemple plusieurs jours. L'exécution d'un algorithme 5 d'identification permet de déterminer la valeur d'au moins une grandeur thermique 31, 32 de la zone considérée, comme expliqué en détail ci-après. La grandeur thermique peut être la résistance thermique d'une paroi (par exemple Rth1, Rth2 pour la zone 11 et Rth3, Rth4 pour la zone 12), la capacité thermique de la zone (par exemple C pour la zone 11 et C' pour la 10 zone 12),.... Une caractéristique importante du procédé selon l'invention est la détermination locale des grandeurs thermiques d'une zone. En effet, les mesures des données thermiques environnementales sont effectuées au 15 niveau d'une zone, en vue de qualifier la valeur de la grandeur thermique des parois et ouvrants de cette zone et non de déterminer la valeur de la grandeur thermique relative au bâtiment dans son ensemble. Une autre caractéristique importante du procédé selon l'invention est qu'il peut comprendre une étape de mesure de la au moins une donnée 20 thermique environnementale sur une période prédéterminée à partir d'une instrumentation simple et non intrusive. En effet, cette étape de mesure ne s'appuie que sur des mesures de températures et ne nécessite aucune mesure additionnelle (telle que la mesure d'un flux thermique par exemple). Une autre caractéristique importante du procédé selon l'invention 25 est que l'étape de mesure peut être réalisée sur une période prédéterminée assez courte, par exemple quelques jours, ou une semaine, et de préférence deux semaines. Par ailleurs, pendant la période prédéterminée sur laquelle a lieu l'étape de mesure, la zone considérée peut être soit passive (c'est-à-dire non occupée, avec une évolution non contrôlée de la température), soit faire 30 l'objet d'apports thermiques (chauffage ou refroidissement) à l'aide d'un système de chauffage ou de refroidissement contrôlé, auquel cas l'apport de puissance issu du chauffage doit être connu. La zone considérée peut également être occupée, et dans ce cas, afin de ne pas perdre en précision lors de la détermination des grandeurs thermiques, il est souhaitable de 35 connaître les apports (chauffage, aération, etc.). Il en résulte une flexibilité en 3062 723 9 termes d'occupation de la zone puisqu'il n'est pas nécessaire d'évacuer la zone. Enfin, l'étape d'exécution de l'algorithme d'identification du procédé selon l'invention nécessite peu de puissance de calcul puisqu'il peut 5 être mis en oeuvre sur un ordinateur de bureau et/ou à distance. La figure 2 représente schématiquement les étapes d'un mode de réalisation d'un procédé de détermination d'une valeur d'au moins une grandeur thermique d'une zone d'un bâtiment selon l'invention. 10 Procédé, mis en oeuvre par ordinateur, de détermination d'une valeur d'au moins une grandeur thermique d'une zone d'un bâtiment soumise à au moins une donnée thermique environnementale et ayant un comportement thermique réel. Selon l'invention, le procédé comprend une étape 203 de construction d'un modèle thermique initial 204 de la zone ayant un 15 comportement thermique théorique, le modèle thermique initial 204 utilisant la valeur de la au moins une grandeur thermique Rth1, Rth2, C, la valeur de la au moins une grandeur thermique étant initialisée à une valeur initiale Rth10, Rth20, Co. L'étape 203 de construction du modèle thermique initial peut être réalisée par analogie électrique, comme expliqué ultérieurement. 20 Préalablement à l'étape 203 de construction du modèle thermique initial de la zone, le procédé selon l'invention peut comprendre une étape 201 de relevé de caractéristiques géométriques basiques de la zone du bâtiment. Le relevé des caractéristiques géométriques de la zone concernée comprend notamment le relevé de la dimension des murs, du sol, du plafond 25 et des ouvrants. Les parois sont supposées planes. Il n'est pas nécessaire de connaître le positionnement du ou des ouvrants dans les parois ainsi que l'épaisseur des murs. Ce relevé peut être effectué sur site ou les caractéristiques géométriques de la zone peuvent être extraites d'un plan ou d'un modèle géométrique préexistant. Et l'étape 203 de construction du 30 modèle thermique initial est basée sur le relevé des caractéristiques géométriques de la zone du bâtiment. Le modèle est créé sur la base des caractéristiques géométriques relevées. Les valeurs des grandeurs thermiques, que l'on cherche donc à déterminer, sont initialisées de façon arbitraire. L'étape 203 de construction du modèle thermique initial peut 35 comprendre une étape de simplification visant à fusionner les parois afin de 3062 723 10 focaliser l'identification sur les paramètres d'intérêt et réduire le temps de calcul. Par exemple, dans le cas d'une zone à six parois latérales dont une seule paroi donne sur l'extérieur, si on souhaite déterminer les valeurs des grandeurs thermiques de la paroi donnant sur l'extérieur, il est possible de 5 fusionner dans le modèle thermique les cinq autres parois latérales intérieures, ainsi que le plafond et le sol, en une seule paroi intérieure. Le procédé selon l'invention comprend une étape 210 d'entrée dans le modèle thermique initial 204 de la zone d'au moins une donnée thermique environnementale issue d'une mesure 209. 10 Préalablement à l'étape 210 d'entrée de la au moins une donnée thermique environnementale issue d'une mesure, le procédé selon l'invention peut avantageusement comprendre une étape 209 de mesure de la au moins une donnée thermique environnementale sur une période prédéterminée pour obtenir la au moins une donnée thermique environnementale issue 15 d'une mesure. La au moins une donnée thermique environnementale peut être les températures et le rayonnement. Les températures comprennent la température d'ambiance de la zone considérée, la température d'ambiance des zones adjacentes à la zone considérée et la température extérieure. Ces données peuvent être mesurées à l'aide de l'instrumentation du bâtiment si 20 elle le permet, ou à l'aide de capteurs de températures et d'enregistreurs de données déployés spécifiquement pour le procédé de détermination. La durée de la mesure des données thermiques environnementales peut être de quelques jours et avantageusement supérieure ou égale à quinze jours. Le rayonnement est le rayonnement global, obtenu soit à partir de données 25 météorologiques de la zone géographique dans laquelle se situe le bâtiment, soit à partir des données mesurées par la station météorologique du bâtiment. La période d'échantillonnage des mesures de températures doit être horaire a minima, avantageusement de trente minutes et 30 préférentiellement de quinze minutes. Le procédé selon l'invention comprend une étape 205 d'exécution d'un algorithme d'identification au modèle thermique initial modifiant la valeur initiale de la au moins une grandeur thermique Rth10, Rth20, Co jusqu'à 35 obtenir un modèle thermique calibré ayant un comportement thermique 3062 723 11 sensiblement égal au comportement thermique réel. Autrement dit, l'étape 205 est itérative jusqu'à l'obtention d'un modèle thermique calibré qui, soumis aux données thermiques environnementales, résulte en un comportement thermique proche du comportement thermique réel de la 5 zone. Par exemple, il est souhaitable que le modèle thermique calibré présente une évolution de sa température intérieure de la zone en fonction des données thermiques d'environnement sensiblement identique à l'évolution constatée pendant l'étape 209 de mesure. Pour une précision maximale des résultats obtenus, on peut exiger que l'évolution thermique du 10 modèle thermique calibrée soit exactement la même que celle de la zone. Néanmoins, afin de réduire le temps de calcul, il est possible de tolérer une différence que l'utilisateur peut spécifier, comme il est usuel de le faire dans ce type de calculs. L'étape 205 d'exécution d'un algorithme consiste à exécuter un 15 algorithme de recherche d'une solution optimale avec par exemple une approche de régression. Enfin, le procédé selon l'invention comprend une étape 207 d'extraction de la valeur 208 de la au moins une grandeur thermique Rthl , 20 Rth2, C du modèle thermique calibré. La au moins une grandeur thermique peut être la résistance thermique d'une paroi de la zone et/ou la capacité thermique de la zone. Le procédé selon l'invention peut comprendre, après l'étape 207 25 d'extraction de la valeur 208 de la au moins une grandeur thermique Rthl , Rth2, C du modèle thermique calibré une étape d'affichage de la valeur de la au moins une grandeur thermique. L'affichage peut être fait sur l'écran de l'ordinateur. Il est ainsi aisé de comparer la valeur des grandeurs thermiques avec des valeurs normées pour vérifier la conformité à un standard ou label 30 donné de la zone ou plus généralement du bâtiment. Comme déjà mentionné, l'étape 209 de mesure de la au moins une donnée thermique environnementale sur une période prédéterminée peut comprendre la mesure d'au moins une température intérieure de la 35 zone et une température extérieure à la zone. L'étape 209 de mesure de la 3062 723 12 au moins une donnée thermique environnementale sur une période prédéterminée peut comprendre la mesure du rayonnement global du bâtiment. 5 La figure 3 représente schématiquement un exemple d'un modèle thermique initial utilisant les grandeurs thermiques à déterminer selon l'invention. L'exemple présenté sur la figure 3 est un modèle thermique par analogie électrique, qui sera désigné dans la suite par l'abréviation MTAE. Le modèle MTAE repose sur une analogie électrique permettant de représenter 10 tout ou partie d'un bâtiment sous la forme d'un circuit électrique Résistance-Capacité (RC). L'intérêt de ce modèle réside en particulier dans sa flexibilité : il permet de modéliser tout le bâtiment ou seulement une partie de ce dernier, mais également de faire varier le niveau de détails de la modélisation (nombre de parois modélisées, nombre de couches dans les 15 parois, etc.). Dans son principe, cette approche permet de représenter la topologie du bâtiment à l'aide d'un graphe non dirigé G où un noeud représente un point de mesure de la température. Par exemple le noeud 42 correspond à la température extérieure de la paroi extérieure 41 et le noeud 43 correspond à la température ambiante de la zone considérée. Ce graphe 20 est composé de noeuds N et d'arcs E. Chaque noeud iEN est associé à une température Ti et à une capacité Ci, et chaque arête a=(iEN,jEN)EE est associée à une résistance R(i,j) telle que R(i,j)=R(j,i). Ainsi, la résistance 44 est égale à la résistance 45 et correspond à la résistance thermique de la paroi extérieure 41. Les parois intérieures peuvent être définies séparément 25 ou bien fusionnées entre elles, toutes ou quelques unes seulement. Sur la figure 3, les parois intérieures ont été fusionnées en une paroi intérieure 46. A cette paroi intérieure 46 peut être associé un ouvrant 47. De même, un ouvrant extérieur 48 peut être défini. Il est aussi possible d'ajouter d'autres éléments dans le modèle, comme par exemple la ventilation 49. 30 La dynamique du circuit est décrite par un système d'équations différentielles du premier ordre de la forme : , a ri c.7t (t) = A i(t) + ET)>(t) 3062723 13 où ri(t) = [Mt), T2 (t), ... , Tn(t)r est le vecteur des températures au temps t, E> = [C1, C2, ..., Ce]' le vecteur des capacités, et cl) = [4,1,1)2, ...,(pnr le vecteur des flux (radiatif, solaire et net). Les éléments de la matrice A = (Aii, i, j E N) sont donnés par: 5 si i # j et (i, j) E E 1 si i # j et (i, j) E E si i = j M A.. - Ce modèle peut être, en fonction des objectifs visés et des caractéristiques de la zone, plus ou moins détaillé. L'approche MTAE offre en effet la possibilité de moduler le nombre de parois modélisées en 10 fusionnant certaines parois entre elles. L'intérêt de la fusion de parois, comme mentionné précédemment, est de réduire le nombre de paramètres à identifier et, ainsi, de garantir la convergence des algorithmes d'identification et de limiter le temps de calcul. La configuration la plus courante consiste à modéliser indépendamment la ou les parois donnant sur l'extérieur, mais à 15 fusionner l'ensemble des parois internes, pour être en mesure de déduire la résistance thermique de la ou des paroi(s) externe(s) et la capacité totale de la zone. Dans ce cas, la zone est alors modélisée sous la forme de deux parois (et ouvrants associés), une externe et une interne représentant la fusion des parois internes. Comme on peut le voir, sur la figure 3, les murs 20 sont supposés être monocouches et sont modélisés à l'aide d'un circuit 2R3C (2 résistances, 3 capacités). Il est à noter que cet exemple est donné uniquement à titre d'illustration : l'approche peut en effet s'appliquer avec des approches alternatives de modélisation (exemple : modélisation multicouche). 25 Une fois les étapes 201 (relevé des caractéristiques géométriques de la zone), 203 (construction du modèle thermique initial de la zone), 209 (mesure des données thermiques environnementales) effectuées, il faut ensuite réaliser l'étape 205 d'exécution d'un algorithme d'identification afin 30 d'obtenir la valeur des grandeurs thermiques recherchées. Ces dernières peuvent varier selon les données disponibles en entrée, mais a minima, il 3062723 14 s'agit du coefficient de transfert de chaleur de chaque paroi modélisée (U=1/RS, où R est la résistance de la paroi et S sa surface), de la capacité (C) de chaque paroi, du coefficient de transfert de chaleur (U) de chaque ouvrant. 5 L'invention peut se baser sur tout type d'algorithme d'identification : de nombreuses approches existent et, en fonction des objectifs visés en priorité, sont plus ou moins adaptées. Dans le cadre de la mise en oeuvre expérimentale de l'invention, trois approches ont été 10 étudiées : l'algorithme PSO (pour son acronyme anglo-saxon Particle Swarm Optimization), ainsi que deux approches basées sur des modèles de régression (une combinaison des Support Vector Machines (SVR) et de l'algorithme PSO, et le krigeage). Une étude comparative basée sur un cas d'étude simple ayant suggéré la supériorité du krigeage en termes de 15 robustesse et de performances, c'est ce dernier algorithme qui est privilégié. Le krigeage est une technique d'interpolation statistique qui permet d'approximer un modèle numérique donné y(x) à l'aide de la combinaison d'une fonction de régression et d'un processus stochastique : K(x) = F(x) + z(x) 20 où K est le modèle de krigeage (i.e. méta-modèle) approximant la réponse réelle K(x) E I': pour un vecteur d'entrée de dimension n:x ED c 11:11, F(x) la fonction de régression et z(x) le processus stochastique, D c II:11 représente les contraintes (bornes) s'appliquant à la variable x. Dans le 25 cadre de l'invention, le modèle numérique dont on souhaite approximer le comportement est l'erreur y(x) E ii: entre la température mesurée dans la zone visée et la température calculée à l'aide du modèle MTAE. Si l'on note TmTAE(x) la température calculée par le modèle MTAE, OÈIXEDC II:11 est le vecteur des paramètres à identifier dans le modèle MTAE (résistances, 30 capacités...) et D c Kri est l'espace des solutions possibles physiquement, et Tref la température mesurée de la zone visée, on cherche donc à approximer l'erreur y(x) telle que: Il y(x) = TmTAE(x) -Tref II à l'aide du modèle de krigeage K(x) tel que : 3062723 15 K (x) = F(x) + z(x) Dans le cadre du prototypage de l'invention, la fonction de régression F(x) est un modèle de régression polynomial consistant en une combinaison linéaire de p fonctions f(x) polynomiales combinant les 5 variables d'entrée : (X) = (x) OÙ (fli,...,flp) E P sont les coefficients de régression. La partie stochastique quant à elle est un processus gaussien centré caractérisé par sa fonction de covariance cov(Z(x),Z(u)) = 0-2R(x,u) 10 où o-2 représente la variance de Z et R(x,u) représente la fonction de corrélation. Cette fonction de corrélation, dans le cadre du prototypage, est une fonction de corrélation exponentielle. Sur la base de ces hypothèses, le principe de l'algorithme de 15 krigeage est de déterminer un ensemble de modèles de krigeage K valides (i.e. respectant un critère d'erreur prédéfini), puis d'identifier au sein de cet ensemble le modèle de krigeage donnant le meilleur résultat, et enfin, d'optimiser les paramètres d'entrée pour minimiser l'erreur. Ce processus est répété jusqu'à atteindre un critère de convergence prédéfini (seuil sur 20 l'erreur, nombre d'itérations maximal). Plus précisément, l'algorithme comprend les étapes suivantes : 1. Sélection au sein de l'espace des solutions d'un ensemble r d'échantillons des paramètres recherchés (x1, ..., xr) où xi E D c Pour chacun de ces échantillons, exécution du modèle MTAE et 25 génération d'un vecteur (yi, , yr) d'erreurs par rapport au modèle MTAE. 2. Création et calibrage d'un ensemble de modèles de krigeage Ki, i = 1 r, chaque modèle Ki, i = 1 r étant calibré à l'aide de l'ensemble des échantillons à l'exception de xi. 30 3. Vérification de la validité des modèles de krigeage à l'aide d'une validation croisée : chaque modèle est testé sur la base de l'échantillon xi et le critère de validité se base sur l'erreur constatée. 3062723 16 o Dans le cas où l'un (ou plusieurs) des modèles de krigeage est (sont) invalide(s), on ajuste l'échantillon (xl, xr) et on exécute de nouveau le modèle MTAE sur les éléments de l'échantillon modifiés. On réitère ensuite l'étape 3 et ce 5 jusqu'à l'obtention d'un échantillon valide. 4. Recherche au sein de l'échantillon du modèle de krigeage le plus performant et sur la base de ce dernier, recherche des valeurs des paramètres xmin minimisant l'erreur (algorithme du gradient). o Dans le cas où le critère de convergence (critère basé sur 10 l'erreur) n'est pas atteint, on étend l'échantillon et on relance l'algorithme. La figure 4 représente un exemple de zone 11 à laquelle le procédé selon l'invention a été appliqué pour en déterminer la valeur des 15 grandeurs thermiques. Il s'agit de déterminer les caractéristiques thermiques d'une pièce (référencée D2.5) d'un bâtiment scolaire 10 situé à Cardiff (Royaume-Uni). Cette zone 11 est de forme rectangulaire, possède une paroi 51 donnant sur l'extérieur (orientation nord, nord-ouest) équipée d'un ouvrant, et trois parois opaques 52, 53, 54 donnant sur des zones 20 adjacentes. L'ensemble des relevés géométriques (hauteur et largeur des parois et ouvrants) ont été effectués et utilisés pour la configuration initiale du modèle MTAE. Les données de températures : température extérieure, température d'ambiance de la zone ciblée et des zones adjacentes ont été 25 mesurées pendant deux semaines (quinze jours). La zone n'était pas occupée mais chauffée, l'apport de puissance issu du chauffage étant connu. Le rayonnement a été négligé en raison de l'orientation (nord, nord-ouest) et de la saison. Les mesures ont ensuite été utilisées pour identifier les 30 paramètres du modèle MTAE. Les paramètres visés étaient : la résistance thermique R (K.m2/VV) et la capacité thermique C (kJ/m2/K). Deux configurations du modèle MTAE ont été étudiées. La première, dans laquelle toutes les parois 51, 52, 53, 54 ont été fusionnées, visait à déterminer la capacité de la zone thermique 11 (Czone). La 35 deuxième, dans laquelle l'ensemble des parois internes 52, 53, 54 ont été 3062723 17 fusionnées mais la paroi externe 51 est conservée indépendante, visait à déterminer la résistance de la paroi externe 51 (Rext). Les résultats obtenus ont été les suivants : Czone = 27.9 kJ/m2/K (configuration 1) 5 Rext = 1.7 K.m2/VV (configuration 2). Afin de valider ces résultats, la résistance thermique de la paroi externe 51 a été calculée à l'aide de la méthode ISO 9869-1 (« Isolation thermique -- éléments de construction -- mesurage in-situ de la résistance 10 thermique et du coefficient de transmission thermique -- Partie 1 : méthode du fluxmètre », aussi disponible en version anglaise : « Thermal insulation -Building elements -- In-situ measurement of thermal resistance and thermal transmittance -- Part 1: Heat flow meter method », accessible depuis le site internet : http://www.iso.org/iso/catalogue_detail.htm?csnumber=59697). La 15 capacité de la zone 11 a été calculée à l'aide de la méthode QUB/e (disponible dans la publication Comptes Rendus Physique, 13(4), pp. 383390, « Quick measurements of energy efficiency of buildings » de Mangematin, E., Pandraud, G. & Roux, D., 2012). Les résultats confirment la précision des calculs réalisés à l'aide du procédé selon l'invention : 20 Czone = 27.6 +/- 0.8 kJ/m2/K (QUB/e) Rext = 1.7 +/- 0.15 K.m2/VV (ISO 9869-1) Comme mentionné précédemment, plusieurs applications de l'invention sont envisageables. Elle peut être mise en oeuvre pour la 25 qualification des performances thermiques de l'enveloppe ou une partie d'un bâtiment ancien (par exemple en vue de la rénovation totale ou partielle de ce bâtiment ancien), dont les caractéristiques sont non ou imparfaitement connues. L'invention peut aussi être mise en oeuvre pour la qualification des performances thermiques de l'enveloppe ou une partie d'un bâtiment 30 nouvellement construit, en vue de vérifier que ces performances sont conformes aux spécifications données en conception et, potentiellement, d'identifier des défauts de construction (garantie de performances). Dans ce cadre, l'invention peut donc servir de base à un processus de certification des performances thermiques du bâtiment en entier ou d'une zone du 35 bâtiment. 3062723 18 L'invention permet de calculer la résistance thermique des parois par opposition aux approches de l'art antérieur qui visent uniquement à fiabiliser les prédictions réalisées par le modèle calibré. Elle permet de réaliser une détermination de la résistance thermique des parois locales par 5 opposition aux approches de l'art antérieur qui ne permettent que de caractériser les grandeurs globales. De plus, le procédé selon l'invention ne s'appuie que sur les mesures de températures d'ambiance (température à l'intérieur) de la zone considérée et des zones limitrophes, et de la température extérieure, par opposition aux approches de l'art antérieur qui 10 nécessitent des données additionnelles telles que le taux d'humidité, la consommation par exemple et a fortiori aux approches qui nécessitent une instrumentation avancée : température de surface et/ ou flux de chaleur. En outre, le procédé selon l'invention s'appuie sur des relevés de caractéristiques géométriques simples à réaliser dans la zone considérée 15 (hauteur, largeur des murs et ouvrants) par opposition aux approches de l'art antérieur nécessitant une connaissance de la géométrie complète du bâtiment. L'invention concerne aussi un dispositif de détermination d'une 20 valeur d'au moins une grandeur thermique d'une zone d'un bâtiment soumise à au moins une donnée thermique environnementale et ayant un comportement thermique réel, le dispositif comprenant des moyens pour mettre en oeuvre les étapes du procédé selon l'invention décrit précédemment. 25 L'invention concerne également un produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes du procédé de détermination selon l'invention lorsque ledit programme est exécuté sur un ordinateur. 30 L'homme du métier considèrera que la présente invention peut s'implémenter à partir d'éléments matériel et/ou logiciel et opérer sur un ordinateur, un tel ordinateur comprenant des moyens permettant d'opérer les étapes du procédé. L'ordinateur peut être couplé à un système de 35 numérisation, un scanner, un appareil photo. L'invention peut aussi opérer 3062723 19 sur une tablette équipée d'une caméra pour l'acquisition des caractéristiques géométriques de la zone du bâtiment. L'invention peut être disponible en tant que produit programme d'ordinateur sur un support lisible par un ordinateur. Le support peut être électronique, magnétique, optique, électromagnétique 5 ou être un support de diffusion de type infrarouge. De tels supports sont par exemple, des mémoires à semi-conducteur (Random Access Memory RAM, Read-Only Memory ROM), des bandes, des disquettes ou disques magnétiques ou optiques (Compact Disk - Read Only Memory (CD-ROM), Compact Disk - Read/VVrite (CD-R/VV) et DVD). 10 REVENDICATIONS1. Procédé, mis en oeuvre par ordinateur, de détermination d'une valeur d'au moins une grandeur thermique (Rth1, Rth2, C) d'une zone (11, 12) d'un bâtiment (10) soumise à au moins une donnée thermique environnementale et ayant un comportement thermique réel, caractérisé en ce qu'il comprend : - Une étape (201) de relevé des caractéristiques géométriques de la zone (11, 12) du bâtiment (10), - Une étape (203) de construction d'un modèle thermique initial de la zone ayant un comportement thermique théorique, le modèle thermique initial utilisant la valeur de la au moins une grandeur thermique (Rth1, Rth2, C), la valeur de la au moins une grandeur thermique étant initialisée à une valeur initiale (Rth10, Rth20, Co), l'étape (203) de construction du modèle thermique initial étant basée sur le relevé des caractéristiques géométriques de la zone (11, 12) du bâtiment (10), - Une étape d'entrée dans le modèle thermique initial de la zone d'au moins une donnée thermique environnementale issue d'une mesure, - Une étape (205) d'exécution d'un algorithme d'identification au modèle thermique initial modifiant la valeur initiale (Rth10, Rth20, Co) de la au moins une grandeur thermique jusqu'à obtenir un modèle thermique calibré ayant un comportement thermique sensiblement égal au comportement thermique réel, - Une étape (207) d'extraction de la valeur de la au moins une grandeur thermique (Rth1, Rth2, C) du modèle thermique calibré. 2. Procédé de détermination selon la revendication 1, caractérisé en ce qu'il comprend après l'étape (207) d'extraction de la valeur de la au moins une grandeur thermique du modèle thermique calibré une étape d'affichage de la valeur de la au moins une grandeur thermique. 3. Procédé de détermination selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que l'étape (203) de construction du modèle thermique initial est réalisée par analogie électrique. 3062723 21 4. Procédé de détermination selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend, préalablement à l'étape d'entrée de la au moins une donnée thermique environnementale issue d'une mesure, une étape (209) de mesure de la au moins une donnée 5 thermique environnementale sur une période prédéterminée pour obtenir la au moins une donnée thermique environnementale issue d'une mesure. 5. Procédé de détermination selon la revendication 4, caractérisé en ce que l'étape (209) de mesure de la au moins une donnée thermique 10 environnementale sur une période prédéterminée comprend la mesure d'au moins une température intérieure de la zone (11, 12) et une température extérieure à la zone (11, 12). 6. Procédé de détermination selon l'une quelconque des 15 revendications 4 ou 5, caractérisé en ce que l'étape (209) de mesure de la au moins une donnée thermique environnementale sur une période prédéterminée comprend la mesure du rayonnement global du bâtiment (10). 7. Procédé de détermination selon l'une quelconque des 20 revendications précédentes, caractérisé en ce que l'étape (205) d'exécution d'un algorithme consiste à exécuter un algorithme de recherche d'une solution optimale avec une approche de régression. 8. Procédé de détermination selon l'une quelconque des 25 revendications précédentes, caractérisé en ce que l'algorithme d'identification est le krigeage ou un algorithme Particle Swarm Optimization (PSO) ou un algorithme Support Vector Machines (SVM) ou une combinaison d'un algorithme Particle Swarm Optimization (PSO) et Support Vector Machines (SVM). 30 9. Procédé de détermination selon l'une quelconque des revendications 1 à 8, dans lequel la au moins une grandeur thermique est la résistance thermique (Rth1, Rth2) d'une paroi de la zone (11, 12) et/ou la capacité thermique (C) de la zone (11, 12). 3062723 22 10. Dispositif de détermination d'une valeur d'au moins une grandeur thermique (Rth1, Rth2, C) d'une zone (11, 12) d'un bâtiment (10) soumise à au moins une donnée thermique environnementale et ayant un comportement thermique réel, le dispositif comprenant des moyens pour 5 mettre en oeuvre les étapes du procédé selon l'une quelconque des revendications 1 à 9. 11. Produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes du 10 procédé de détermination selon l'une quelconque des revendications 1 à 9, lorsque ledit programme est exécuté sur un ordinateur.