Brevet FR3049045A1 — publié le 20170922

PROCEDE ET DISPOSITIF DE DETERMINATION DU TAUX DE RENOUVELLEMENT D'AIR D'UN LOCAL La présente invention a trait à un procédé et un dispositif pour déterminer le taux de renouvellement d'air d'un local. Au sens de l'invention, on entend par "local" tout espace délimité par une enveloppe, en particulier un espace d'habitat ou de stockage. Il peut s'agir d'un espace d'habitat ou de stockage fixe, tel qu'une maison individuelle ou un bâtiment, notamment à usage d'habitation ou tertiaire, une partie d'un tel bâtiment, par exemple un appartement dans un immeuble à plusieurs étages, ou encore une cuve de stockage. Il peut également s'agir d'un espace d'habitat ou de stockage transportable, tel qu'une cabine ou une cuve de stockage dans un véhicule, notamment dans une voiture, un camion, un train, un navire, un sous-marin, un avion, un vaisseau spatial. Le taux de renouvellement d'air d'un local, noté ACH, est le volume d'air intérieur du local qui est renouvelé par heure, exprimé en m3.h-1. Il peut s'agir d'un taux de renouvellement d'air par ventilation naturelle, ou d'un taux de renouvellement d'air par ventilation forcée (c'est-à-dire à l'aide d'une ventilation mécanique contrôlée). Pour un bâtiment ancien à rénover, la détermination du taux de renouvellement d'air par ventilation naturelle, sans ventilation mécanique, est utile pour évaluer si le débit d'air est suffisant en regard de la qualité de l'air intérieur et pour connaître l'impact des infiltrations naturelles sur le bilan énergétique du bâtiment. Pour un bâtiment neuf comprenant un système de ventilation mécanique à récupération de chaleur, on peut déterminer le taux de renouvellement d'air par ventilation naturelle, en configuration désactivée de la ventilation mécanique, ce qui permet d'évaluer l'étanchéité du bâtiment ; on peut également déterminer le taux de renouvellement d'air par ventilation forcée, en configuration active de la ventilation mécanique, ce qui permet d'évaluer le débit du système de ventilation mécanique et de vérifier qu'il fonctionne de manière optimale. Il est connu de déterminer le taux de renouvellement d'air dans un bâtiment par un test de type "porte soufflante" ("blower door" en anglais), 3049045 2 dans lequel on obtient une pressurisation ou une dépressurisation mécanique du bâtiment et on mesure les débits d'air qui en résultent dans une plage de différences de pression statique entre l'intérieur et l'extérieur du bâtiment. Un coefficient classiquement mesuré est le débit de fuite sous 50 5 Pa divisé par le volume chauffé, noté n50, qui correspond au volume d'air renouvelé par heure dans le bâtiment soumis à une différence de pression de 50 Pa entre l'intérieur et l'extérieur du bâtiment. Dans cette méthode de "porte soufflante", on installe dans une ouverture du bâtiment une fausse porte formée par une bâche étanche dans 10 laquelle est fixé un ventilateur, puis on met en fonctionnement le ventilateur afin de créer une différence de pression entre l'intérieur et l'extérieur du bâtiment, et on procède à des mesurages du débit d'air et de la différence de pression sur une plage de différences de pression allant généralement de 20 Pa à 70 Pa. Ces différences de pression importantes ne permettent pas de 15 rendre compte de l'effet des conditions extérieures sur le bâtiment. Il est nécessaire d'utiliser des modèles physiques pour extrapoler les résultats à de plus faibles différences de pression et remonter à des paramètres physiques réalistes, un exemple étant le modèle de Persily-Kronvall. Toutefois, il est difficile d'obtenir un résultat satisfaisant de manière simple, 20 les modèles physiques simplifiés n'étant pas adaptés à tous les types de bâtiments et les modèles physiques plus élaborés étant complexes à mettre en oeuvre. De plus, avec cette méthode, il est nécessaire de réaliser la mesure en condition de faible vent, ce qui n'est pas toujours possible. Par ailleurs, en physique du bâtiment, il est connu de déterminer le 25 taux de renouvellement d'air par la méthode du gaz traceur, qui a l'avantage de prendre en compte les conditions naturelles de mesure. Dans la variante la plus classique de cette méthode, on injecte une grande quantité d'un gaz spécifique tel que SF6 ou CO2 à l'intérieur du bâtiment, et on mesure la décroissance de la concentration de ce gaz dans l'air pendant plusieurs 30 heures. Le taux de renouvellement d'air est alors estimé par le rapport de la variation temporelle de la concentration du gaz à l'intérieur du bâtiment et de la différence de quantité absolue de celui-ci entre l'intérieur et l'extérieur. 3049045 3 Toutefois, cette méthode n'est pas très utilisée sur les chantiers après construction, principalement en raison de son coût et des longues durées de mesure. C'est à ces inconvénients qu'entend plus particulièrement remédier 5 l'invention en proposant un procédé et un dispositif permettant de déterminer le taux de renouvellement d'air d'un local de manière simple et rapide, notamment sur une période de quelques heures, avec un coût modéré et une précision raisonnable, pour tout type de local. A cet effet, l'invention a pour objet un procédé de détermination du 10 taux de renouvellement d'air ACH d'un local, caractérisé en ce qu'il comprend des étapes dans lesquelles : - sur au moins deux périodes de temps Di, successives correspondant à des débits d'un gaz donné distincts appliqués dans le local, on procède à une campagne de mesures permettant de déterminer la concentration du gaz 15 à l'intérieur du local ca, à intervalles de temps rapprochés, ainsi qu'à la détermination de la concentration du gaz à l'extérieur du local ceic à intervalles de temps rapprochés ; - on détermine la valeur du taux de renouvellement d'air ACH du local en faisant converger : d'une part, un modèle diffusif exprimant la variation 20 temporelle de la concentration du gaz à l'intérieur du local ca, en fonction de la concentration du gaz à l'extérieur du local ceic et de paramètres physiques du local à partir desquels le taux de renouvellement d'air ACH du local est calculable ; et, d'autre part, l'évolution mesurée cac(t) de la concentration du gaz à l'intérieur du local en fonction du temps. 25 Au sens de l'invention, la "concentration du gaz" peut être la teneur en gaz dans l'air en grammes de gaz par mètre cube d'air (g/m3), ce qui correspond à l'humidité absolue dans le cas du gaz H20, ou la teneur en gaz dans l'air en grammes de gaz par gramme d'air (g/g d'air), ce qui correspond à l'humidité spécifique dans le cas du gaz H20. Le taux de renouvellement 30 d'air ACH du local peut alors être donné soit en m3/h, soit en g d'air/h. 3049045 4 De plus, dans le cadre de l'invention, "l'intérieur du local" et "l'extérieur du local" désignent les deux masses d'air qui sont séparées par l'enveloppe du local et situées respectivement à l'intérieur et à l'extérieur de l'enveloppe du local. 5 Le gaz utilisé dans le cadre de l'invention, dont on applique des débits distincts dans le local et dont on suit l'évolution de la concentration, peut être choisi, notamment, parmi H20, CO2, He, SF6, H2, N2, ou d'autres gaz traceurs réfrigérants, cette liste n'étant pas limitative. Plusieurs critères importants sont à considérer pour le choix du gaz, en particulier le gaz doit 10 être facilement identifiable et mesurable dans l'atmosphère ; le gaz ne doit pas être dangereux lors d'une inhalation longue ou fréquente, même si le procédé ne requiert pas la présence de personnes tout au long des mesures ; de préférence, le gaz n'est pas nocif pour l'environnement. De plus, l'utilisation d'un gaz qui est présent dans l'atmosphère, tel que H20 ou 15 CO2, impose la prise en compte de la concentration de celui-ci dans l'atmosphère et l'utilisation de rectifications lors du calcul du taux de renouvellement d'air. L'utilisation de la vapeur d'eau H20 dans le cadre de l'invention, en tant que gaz dont on applique des débits distincts dans le local et dont on suit l'évolution de la concentration, est avantageuse du fait 20 de son faible coût, de son innocuité et de la facilité à générer un débit de vapeur d'eau. L'invention permet une détermination in situ du taux de renouvellement d'air d'un local. Le principe à la base de l'invention est d'utiliser les variations transitoires de la concentration d'un gaz donné à 25 l'intérieur du local lorsque le local est soumis à des sollicitations internes contrôlées et dans un environnement externe mesuré. L'analyse quantitative de la variation de la concentration du gaz à l'intérieur du local permet de déterminer quantitativement le renouvellement d'air entre l'intérieur et l'extérieur du local sur une période courte, de l'ordre de quelques heures, en 30 limitant le nombre de paramètres susceptibles d'influencer le comportement du local. En particulier, la brièveté des mesures permet de s'affranchir de l'influence des conditions d'utilisation du local, des variations des conditions 3049045 5 extérieures et des phénomènes d'absorption ou de désorption du gaz par les matériaux de construction et/ou le mobilier du local. Dans le cadre de l'invention, on entend par "débit du gaz appliqué dans le local" toute condition opératoire générant une variation de la 5 concentration du gaz à l'intérieur du local, pour des conditions données de concentration du gaz à l'extérieur du local. Il est entendu que le débit du gaz appliqué dans le local peut être positif, nul ou négatif. Un débit du gaz positif correspond à une injection du gaz dans le local, alors qu'un débit du gaz négatif correspond à une extraction du gaz hors du local. Dans le cas d'un 10 débit du gaz nul ou sensiblement nul appliqué dans le local, la variation de la concentration du gaz à l'intérieur du local peut résulter d'un déséquilibre entre les concentrations initiales du gaz à l'intérieur et à l'extérieur du local. Selon l'invention, il est prévu qu'au moins l'un des débits dudit gaz donné est non nul. 15 Le modèle diffusif utilisé pour déterminer le taux de renouvellement d'air du local peut être de tout type connu de l'homme du métier. Il peut s'agir, notamment, d'un modèle R-C avec un nombre adapté de résistances et de capacités. De préférence, le modèle diffusif utilisé pour déterminer le taux de renouvellement d'air du local est un modèle R-C simple avec une 20 résistance et une capacité. En variante, le modèle diffusif peut être, notamment, un modèle R-C plus complexe que le modèle R-C simple tel qu'un modèle dit "3R2C" avec trois résistances et deux capacités. Au sens de l'invention, le fait de faire converger le modèle diffusif et l'évolution mesurée cac(t) signifie que l'on ajuste la valeur de paramètres 25 physiques du local utilisés dans le modèle diffusif de manière à minimiser l'écart, au moins sur un intervalle de temps compris dans chaque période de temps Dk, entre l'évolution temporelle de la concentration du gaz à l'intérieur du local calculée à partir du modèle diffusif et l'évolution temporelle de la concentration du gaz à l'intérieur du local effectivement mesurée cac(t). 30 L'ajustement peut ainsi se faire sur toute l'étendue de chaque période de 3049045 6 temps Dk, ou sur un ou plusieurs intervalles de temps compris dans chaque période de temps Dk. A titre d'exemple, dans le cas où le modèle diffusif est un modèle R-C simple avec une résistance et une capacité et où, pour chaque période de 5 temps Dk, il existe un intervalle de temps Atk pour lequel l'évolution mesurée cik(t) de la concentration du gaz à l'intérieur du local en fonction du temps est sensiblement linéaire, on peut faire converger le modèle R-C simple et l'évolution mesurée cik(t) sur les intervalles de temps Atk de la manière suivante : pour chaque période de temps Dk on détermine sur l'intervalle de 10 temps Atk la pente ak de la tangente à l'évolution cik(t) et on déduit la valeur du taux de renouvellement d'air ACH du local à partir des pentes ak, des débits du gaz appliqués dans le local et de la différence de concentration du gaz entre l'intérieur et l'extérieur du local. Selon un autre exemple, dans le cas où le modèle diffusif est un 15 modèle R-C plus complexe, tel qu'un modèle "3R2C", on fait converger le modèle R-C plus complexe et l'évolution mesurée cik(t) en ajustant la valeur des paramètres physiques du local utilisés dans le modèle de manière à minimiser l'écart (lifting"), sur l'ensemble des périodes de temps Dk, entre l'évolution temporelle de la concentration du gaz à l'intérieur du local calculée 20 à partir du modèle R-C plus complexe et l'évolution temporelle de la concentration du gaz à l'intérieur du local effectivement mesurée cik(t). En pratique, des données d'entrée sont injectées dans le modèle diffusif, comprenant le volume du local, la concentration du gaz à l'extérieur du local cek déterminée sur chaque période de temps Dk. 25 Des exemples de paramètres physiques du local susceptibles d'intervenir dans le modèle diffusif et d'être ajustés de manière à faire converger le modèle diffusif et l'évolution mesurée cik(t) comprennent, notamment, la surface totale d'enveloppe du local, le nombre et la taille des menuiseries, permettant par exemple de déterminer une surface équivalente 30 de trous. 3049045 7 Selon un aspect de l'invention, le procédé comprend des étapes dans lesquelles : - sur au moins deux périodes de temps Dk successives correspondant à des débits du gaz distincts appliqués dans le local, on procède à une 5 campagne de mesures permettant de déterminer la concentration du gaz à l'intérieur du local ca, à intervalles de temps rapprochés, ainsi qu'à la détermination de la concentration du gaz à l'extérieur du local cek à intervalles de temps rapprochés ; - pour chaque période de temps Dk, à partir de l'évolution mesurée cik(t) 10 de la concentration du gaz à l'intérieur du local ca, en fonction du temps : o soit, s'il existe un intervalle de temps At!, pour lequel l'évolution cik(t) est sensiblement linéaire, on détermine sur cet intervalle de temps At!, la pente ak de la tangente à l'évolution cik(t) et on déduit la valeur du taux de renouvellement d'air ACH du local à partir des pentes ak; 15 o soit, s'il n'existe pas d'intervalle de temps pour lequel l'évolution cik(t) est sensiblement linéaire, on sélectionne un intervalle de temps Atk' sur lequel l'évolution cik(t) est sensiblement exponentielle de type exp(-t/ -c), avec 'r le temps au bout duquel le volume d'air intérieur du local a été renouvelé, et on déduit la valeur du taux de renouvellement d'air ACH du 20 local, qui est la valeur telle que l'évolution Lnkt9k(t) - ACH)/(ek(0) - AeicH)1 est une droite, avec Ok(t) = cik(t) - cekm: où Cekm' est la moyenne de la concentration du gaz à l'extérieur du local cek sur l'intervalle de temps Atk'. Bien entendu, le procédé selon l'invention ne requiert pas nécessairement la mise en place d'une représentation graphique de 25 l'évolution cik(t). En particulier, sur chaque intervalle de temps Atk, la pente ak de la tangente à l'évolution cik(t) est égale à la dérivée de l'évolution cik(t) sur l'intervalle Atk. Dès lors, l'étape de détermination de la pente ak de la tangente à l'évolution cik(t) sur l'intervalle Atk peut être réalisée, dans le 30 cadre de l'invention, en calculant la dérivée de l'évolution cik(t) sur 3049045 8 l'intervalle Atk, sans recourir à une représentation graphique de l'évolution cik (t). Les étapes de calcul du procédé, en particulier pour la détermination des pentes ak, peuvent être mises en oeuvre à l'aide de tout moyen de calcul 5 approprié. Il peut s'agir notamment d'un terminal comportant un système d'acquisition pour acquérir les mesures requises par le procédé et des moyens de calcul pour exécuter tout ou partie des étapes du procédé à partir des mesures acquises. Dans le cadre de l'invention, les périodes de temps Dk peuvent être 10 soit disjointes, soit immédiatement successives les unes aux autres. Dans ce dernier cas, on peut considérer que le procédé est réalisé dans sa globalité sur une période de temps continue, formée par la succession des périodes de temps Dk. De préférence, le procédé est mis en oeuvre avec deux périodes de 15 temps D1 et D2 successives correspondant à deux consignes de débits du gaz et distinctes appliquées dans le local. Selon un aspect de l'invention, le premier débit du gaz appliqué sur la première période de temps D1 et le deuxième débit du gaz appliqué sur la deuxième période de temps D2 ont des valeurs éloignées l'une de l'autre. 20 Dans un mode de réalisation avantageux, le premier débit du gaz appliqué sur la première période de temps D1 est strictement positif ou strictement négatif, alors que le deuxième débit du gaz appliqué sur la deuxième période de temps D2 est nul ou sensiblement nul. De manière avantageuse, pour chaque période de temps Dk, le débit 25 du gaz appliqué dans le local comprend un débit Rimpk imposé au moyen d'au moins un appareil à débit contrôlé. On entend ici par "appareil à débit contrôlé" un appareil tel que la quantité du gaz injectée ou extraite par l'appareil peut être déterminée précisément. De préférence, la commande du ou de chaque appareil à débit contrôlé est automatisée. 30 Lorsque le gaz utilisé est la vapeur d'eau, le ou chaque appareil à débit contrôlé peut être un humidificateur, par exemple un humidificateur à 3049045 9 ultrasons. Une méthode pour déterminer précisément le débit de vapeur d'eau injecté sur chaque période de temps Di, est alors de peser l'humidificateur au début et à la fin de la période Di, et, connaissant la durée d'allumage de l'humidificateur, de déduire le débit moyen de vapeur d'eau 5 appliqué dans le local sur la période Dk. Dans le cadre de l'invention, les débits du gaz distincts appliqués dans le local sur les différentes périodes de temps Di, peuvent être différentes consignes de débit du gaz, c'est-à-dire que pour chaque période de temps Di, le débit du gaz est constant sur toute la période Dk. En 10 variante, le débit du gaz peut ne pas être constant sur une ou plusieurs périodes Di, et varier autour d'une valeur de débit du gaz moyen km, pour autant que le débit du gaz moyen km sur la période Di, soit distinct des débits du gaz (constants ou moyens) appliqués sur les périodes de temps qui entourent la période Dk. Dans ce cas, le débit du gaz considéré sur la 15 période de temps Di, est le débit du gaz moyen km. La variation du débit du gaz sur chaque période de temps Di, doit être faible par rapport à la différence entre les débits du gaz appliqués sur deux périodes de temps consécutives, de préférence la variation du débit du gaz sur chaque période de temps Di, est inférieure à 30%, encore de préférence inférieure à 20%, 20 encore de préférence inférieure à 10%, de la différence entre les débits du gaz appliqués sur deux périodes de temps consécutives. Si aucune source du gaz donné autre que les appareils utilisés pour appliquer le débit imposé n'est active dans le local au cours de la période de temps Dk, le débit du gaz appliqué dans le local est égal au 25 débit imposé En revanche s'il existe, au cours de la période de temps Dk, un débit du gaz supplémentaire Rsupk appliqué dans le local en plus du débit le débit du gaz appliqué dans le local est égal à + gsupk Il convient donc, soit d'être capable de mesurer le débit supplémentaire, soit de couper toutes les sources autres que les appareils utilisés pour appliquer 30 le débit imposé. 3049045 10 En particulier, lorsque le gaz utilisé est la vapeur d'eau, un débit de vapeur d'eau supplémentaire dans le local peut provenir de sources telles que la respiration des occupants, la cuisson, la présence de linge qui sèche, l'utilisation d'une douche, l'utilisation d'un sèche-linge. De préférence, le 5 procédé est mis en oeuvre alors que le local est inoccupé. De plus, les éventuels appareils électroménagers sont préférentiellement éteints. Le procédé selon l'invention repose sur le fait que la variation de la concentration du gaz à l'intérieur du local est proportionnelle à la quantité d'air renouvelée par les infiltrations. Si la concentration du gaz à l'intérieur du 10 local n'est pas homogène, il existe un risque qu'une quantité plus importante d'air que du gaz soit échangée, ou vice-versa, ce qui est susceptible de rendre les mesures imprécises. Le même risque peut résulter de stratifications du gaz, causées par des écarts trop importants de température à l'intérieur du local. Afin d'améliorer l'homogénéité de la concentration du 15 gaz et de la température à l'intérieur du local, il est avantageux d'ouvrir les portes de communication à l'intérieur du local et de mettre en place un système de brassage de l'air à l'intérieur du local. Lorsque le local présente un grand volume intérieur, plusieurs appareils à débit contrôlé sont avantageusement répartis dans le local, en étant associés à un système de 20 brassage de l'air. Selon un aspect de l'invention, les mesures permettant de déterminer la concentration du gaz à l'intérieur du local ca, sont effectuées à l'aide d'un ou de plusieurs capteurs dudit gaz placés dans le volume d'air à l'intérieur du local. Dans le cas où le gaz utilisé est la vapeur d'eau, chaque capteur dudit 25 gaz est un capteur d'humidité configuré pour mesurer la concentration de vapeur d'eau, également appelée humidité absolue, à l'intérieur du local ci/c. De préférence, chaque capteur dudit gaz est associé à un capteur de température apte à mesurer la température de l'air à l'intérieur du local. La concentration du gaz à l'intérieur du local cik, et la température à l'intérieur 30 du local Tac, peuvent alors être considérées chacune comme étant la moyenne des mesures des différents capteurs distribués dans le volume d'air à l'intérieur du local. Le nombre de capteurs doit être d'autant plus important 3049045 11 que l'homogénéité de la concentration du gaz et de la température à l'intérieur du local est imparfaite. L'utilisation de plusieurs capteurs dans différents endroits peut également permettre de localiser les fuites d'air, par exemple en repérant qu'un capteur positionné près d'une menuiserie peu 5 étanche a un comportement différent d'un capteur positionné au centre du volume d'air à l'intérieur du local. De préférence, sur chaque période de temps Dk, la température de l'air à l'intérieur du local Ta, est stable. Des variations de la température de l'air à l'intérieur du local peuvent en effet influer sur la concentration du gaz 10 dans l'air à l'intérieur du local. Dès lors, il est préférable, sur chaque période de temps Dk, de couper ou de maintenir à une consigne constante les sources de chaleur ou de froid dans le local, notamment les dispositifs de chauffage ou de climatisation. En particulier, lorsque le gaz utilisé est la vapeur d'eau, une augmentation de la température à l'intérieur du local peut 15 entraîner un séchage des matériaux dans le local, d'où une évaporation de vapeur d'eau dans l'air, tandis qu'une diminution de la température à l'intérieur du local peut entraîner une condensation de vapeur d'eau présente dans l'air à la surface des matériaux. Il est également préférable, sur chaque période de temps Dk, d'avoir 20 un rayonnement solaire faible, de préférence nul, afin d'éviter une augmentation de la température à l'intérieur du local sous l'effet du rayonnement solaire. Les volets ou autres éléments occultants du local peuvent être fermés au cours de chaque période de temps Dk, notamment du côté du local orienté vers le sud. De préférence, le procédé selon 25 l'invention est mis en oeuvre de nuit. Selon un aspect avantageux, en vue de limiter le temps de mise en oeuvre du procédé tout en réduisant la contribution du rayonnement solaire, le procédé est réalisé dans sa globalité en continu sur une seule période nocturne. 30 Pour que l'évolution de l'humidité absolue à l'intérieur du local lors de la mise en oeuvre du procédé soit uniquement liée au renouvellement d'air, il faut se placer dans des conditions telles que les échanges de vapeur d'eau 3049045 12 contenue dans l'air et dans les matériaux de construction et/ou le mobilier sont négligeables. Or, à titre d'exemple, une surface de 50 m2 de chêne non vernis peut stocker environ 2 kg d'eau en 3 heures lorsque l'humidité à l'intérieur du local varie de 60% à 80%. Il est donc favorable d'appliquer les 5 débits de vapeur d'eau distincts et d'effectuer les mesures conformément à l'invention sur des périodes de temps Di, aussi courtes que possible, afin de limiter ces échanges. De préférence, chaque période de temps Di, a une durée inférieure ou égale à 2 heures, plus préférentiellement inférieure ou égale à 1 heure. 10 Si l'on souhaite mesurer le taux de renouvellement d'air par infiltration naturelle d'un local équipé d'un système de ventilation mécanique, la ventilation mécanique doit être désactivée sur chaque période de temps Dk, afin de ne pas introduire un débit d'air additionnel. En variante, il est possible de maintenir la ventilation mécanique en configuration active à un débit 15 souhaité sur les périodes de temps Dk, le procédé selon l'invention permettant alors de mesurer le débit d'air engendré par la ventilation mécanique et ainsi de vérifier son fonctionnement. De préférence, sur chaque période de temps Dk, la concentration du gaz à l'extérieur du local ceic est stable. 20 La concentration du gaz à l'extérieur du local ceic peut être déterminée par des mesures à l'aide de capteurs dudit gaz placés dans l'air à l'extérieur du local, notamment des capteurs d'humidité et de température lorsque le gaz utilisé est la vapeur d'eau. Ces capteurs d'humidité et de température sont placés dans le volume d'air extérieur en évitant d'être à proximité de 25 sources de chaleur ou de production de vapeur d'eau, afin d'avoir des mesures représentatives de la masse d'air globale au voisinage de la surface extérieure de l'enveloppe du local. En variante, lorsque le gaz utilisé est la vapeur d'eau et le local donne directement à l'extérieur, la concentration de vapeur d'eau ou humidité 30 absolue à l'extérieur du local ceic peut être déterminée par interpolation de données météorologiques du lieu du local. 3049045 13 En général, l'air extérieur contient de la vapeur d'eau, surtout en période humide. Le procédé selon l'invention est d'autant plus fiable que les débits de vapeur d'eau appliqués dans le local conformément au procédé de l'invention sont adaptés aux conditions extérieures et que la différence 5 d'humidité absolue entre l'intérieur et l'extérieur du local est grande. Dans le cadre de l'invention, on peut utiliser un modèle R-C simple pour décrire un local, avec deux noeuds de concentration du gaz homogène (en g.m-3), l'un à l'intérieur du local et l'autre à l'extérieur du local, qui sont séparés par une résistance représentant l'inverse du taux de renouvellement 10 d'air ACH du local (en m3.1-1-1). Le noeud de concentration du gaz à l'intérieur du local est relié à un condensateur qui représente le volume effectif du local (en m3). Le débit du gaz appliqué dans le local (en g.1-1-1) est compensé par l'échange du gaz à travers l'enveloppe du local et la quantité du gaz emmagasinée dans la structure de l'enveloppe, ce qui est décrit par 15 l'équation : ak = ACH(cik-Cek) + V dt où cik est débit du gaz total appliqué dans le local, cik et cek sont respectivement la concentration du gaz à l'intérieur du local et la concentration du gaz à l'extérieur du local, ACH est le taux de renouvellement d'air du local et V est le volume effectif du local. 20 On suppose que la réponse du local est une simple exponentielle décroissante et que sa constante de temps est le produit du taux de renouvellement d'air ACH et du volume effectif V du local. En réalité, la réponse du local est plus complexe et est la superposition d'un grand nombre d'exponentielles décroissantes, mais il a été validé 25 expérimentalement qu'en adaptant les conditions du test, en particulier la durée du test et la valeur de débit de gaz appliqué dans le local, seule la constante de temps la plus grande joue un rôle et le modèle décrit précédemment est valide. CICik 3049045 14 En appliquant deux débits du gaz et de valeurs différentes dans le local sur deux périodes de temps D1 et D2, il est alors possible de déterminer le taux de renouvellement d'air ACH du local selon l'équation : a1 X R2 - a2 X V = (2)a1 X Ac2,, - a2 X Acim 5 où (Rk)k-i ou 2 est le débit du gaz appliqué sur la période de temps Dk, (a1c)k=1 ou 2 est la pente sur l'intervalle de temps Atk de la tangente à l'évolution cik(t) de la concentration du gaz à l'intérieur du local en fonction du temps, et (ACkm)k=1 ou 2 est la différence entre la concentration du gaz moyenne à l'intérieur du local et la concentration du gaz moyenne à 10 l'extérieur du local sur l'intervalle de temps Atk. Selon un aspect de l'invention, suite à la détermination de la valeur calculée Vcaic du volume effectif V du local conformément au procédé, on vérifie que cette valeur Vcaic correspond bien au volume réel du local, notamment à 20% près. Si ce n'est pas le cas, il est nécessaire de modifier 15 le traitement des mesures pour obtenir une valeur corrigée du taux de renouvellement d'air ACH du local. En effet, il a été constaté que le procédé selon l'invention peut conduire à une valeur calculée Vcaic du volume effectif du local différente du volume réel du local. En particulier, lorsque le gaz utilisé est de la vapeur 20 d'eau injectée dans le local, il a été observé que, quand la durée de chaque période de temps Dk est trop courte, la valeur calculée Vcaic est plus faible que le volume réel, alors que quand la durée de chaque période de temps Dk est trop longue, la valeur calculée Vcaic est plus importante que le volume réel. Cela peut s'expliquer par le fait que la valeur calculée Vcaic du volume 25 effectif V du local ne correspond pas au volume réel du local, mais plutôt au volume effectivement sollicité par le gaz utilisé dans le cadre du procédé. Une valeur calculée Vcaic sous-estimée par rapport au volume réel du local peut indiquer que le gaz n'a pas eu le temps de sonder tout le volume du ACH = a1 X Ac2,, - a2 X Acim et le volume effectif V du local selon l'équation : X Ac27, - a2 X ACim (1) 3049045 15 local, tandis qu'une valeur calculée Vcaic surestimée par rapport au volume réel du local peut indiquer un effet lié à l'absorption du gaz par les matériaux présents dans le local. Il convient alors, pour le traitement des mesures obtenues selon le 5 procédé de l'invention, de se placer à des durées suffisamment longues pour le que gaz ait sollicité tout le volume du local, mais suffisamment courtes pour éviter tout effet lié à l'absorption du gaz par les matériaux de construction et/ou le mobilier du local, et ainsi ne mesurer que le renouvellement d'air responsable de la variation de la concentration du gaz 10 dans l'air. Dans le cas où le gaz utilisé est de la vapeur d'eau injectée dans le local, il a été constaté qu'une durée favorable pour le traitement des mesures sur chaque période de temps Di, est de l'ordre de 10 minutes à 1 heure. De manière très avantageuse, la comparaison de la valeur calculée 15 Vcaic du volume effectif du local par rapport au volume réel permet ainsi de valider si le taux de renouvellement d'air ACH du local a été déterminé correctement, d'aider à définir des durée favorables pour le traitement des mesures sur chaque période de temps Di, en fonction des matériaux présents dans le local, et même de corriger a posteriori la valeur calculée 20 ACHcaic du taux de renouvellement d'air du local dans le cas où une durée trop longue a été appliquée pour le traitement des mesures. En particulier, on peut procéder à une telle correction a posteriori de la valeur ACHccac en déplaçant les intervalles de temps Atic dans chaque période de temps Dk, de manière à réduire la durée de traitement des mesures, jusqu'à obtenir une 25 valeur du volume effectif Valu sensiblement égale au volume intérieur réel du local. Selon un mode de réalisation, le procédé comprend des étapes dans lesquelles : - on procède, sur deux périodes de temps D1 et D2 successives : 30 i. sur la première période de temps D1, à l'application dans le local d'un premier débit du gaz et à une campagne de mesures permettant de 3049045 16 déterminer la concentration du gaz à l'intérieur du local cil à intervalles de temps rapprochés, ainsi qu'à la détermination de la concentration du gaz à l'extérieur du local cet à intervalles de temps rapprochés, le premier débit du gaz ACHre f .ACi(0) étant tel que le paramètre a = 1 est inférieur ou égal à 6/1 5 0,8, avec Ac1(0) = c11(0) - cern, où t = 0 est le point de départ de la première période de temps D1, Cern est la concentration moyenne du gaz à l'extérieur du local sur l'ensemble des périodes de temps D1 et D2, et ACH'f est une valeur de référence du taux de renouvellement d'air du local, puis ii. sur la deuxième période de temps D2, à l'application dans le local d'un 10 deuxième débit du gaz sensiblement nul, et à une campagne de mesures permettant de déterminer la concentration du gaz à l'intérieur du local cil à intervalles de temps rapprochés, ainsi qu'à la détermination de la concentration du gaz à l'extérieur du local cet à intervalles de temps rapprochés ; 15 - on détermine la valeur du taux de renouvellement d'air ACH du local en faisant converger : d'une part, un modèle diffusif exprimant la variation temporelle de la concentration du gaz à l'intérieur du local ca, en fonction de la concentration du gaz à l'extérieur du local ceic et de paramètres physiques du local à partir desquels le taux de renouvellement d'air ACH du local est 20 calculable ; et, d'autre part, l'évolution mesurée cac(t) de la concentration du gaz à l'intérieur du local ca, en fonction du temps. Dans ce mode de réalisation, on sélectionne une sollicitation spécifique du local, qui permet d'accéder au taux de renouvellement d'air ACH du local avec une bonne précision et sur un temps réduit, cette 25 sollicitation spécifique étant l'application d'un premier débit du gaz strictement positif ou strictement négatif propre à générer une évolution forcée de la concentration du gaz à l'intérieur du local ca, suivie de l'application d'un deuxième débit du gaz sensiblement nul donnant lieu à une évolution libre de la concentration du gaz à l'intérieur du local cil. 3049045 17 De préférence, le premier débit du gaz est tel que le paramètre ACHre f (CO a = 1 est supérieur ou égal à 0,3, de préférence supérieur ou égal à 0,4. En effet, pour des milieux présentant peu d'infiltrations, lorsque le paramètre a est inférieur à 0,3 ou 0,4, la sensibilité des capteurs de mesure 5 classiques ne permet pas d'obtenir des données satisfaisantes concernant l'évolution de la concentration du gaz à l'intérieur du local c,' sur la première période de temps D1, d'où une augmentation de l'incertitude sur la valeur du taux de renouvellement d'air ACH du local déterminé selon l'invention. La détermination de la valeur du premier débit du gaz à appliquer 10 sur la première période de temps D1 pour satisfaire aux critères sur le paramètre a impose de connaître une valeur de référence ACH'f du taux de renouvellement d'air ACH du local. Une méthode pour accéder à une valeur de référence ACH'f du taux de renouvellement d'air ACH du local est l'utilisation d'une grandeur issue 15 d'un test de type "porte soufflante" du local. D'autres méthodes pour accéder à une valeur de référence ACH'f sont également envisageables, en particulier la valeur de référence peut être la valeur de renouvellement d'air indiquée dans l'étude thermique déposée avec le permis de construire. Dans un mode de réalisation avantageux, pour chaque période de 20 temps D1 et D2, il existe un intervalle de temps At, ou At2 pour lequel l'évolution mesurée (c,k(t))k=1 ou 2 de la concentration du gaz à l'intérieur du local en fonction du temps est sensiblement linéaire, et on fait converger un modèle diffusif R-C et l'évolution mesurée (C 'kW)Ic=1 ou 2 de la manière suivante : pour chaque période de temps D k on détermine sur l'intervalle de 25 temps Atk la pente ak de la tangente à l'évolution cik(t), et on détermine la valeur du taux de renouvellement d'air ACH du local à partir des pente ak et des débits du gaz appliqués dans le local. De préférence, les intervalles de temps At, et At2 sur lesquels on effectue le traitement des mesures présentent une "symétrie", où on entend 30 par "symétrie" des intervalles de temps At, et At2 le fait que les deux 3049045 18 intervalles ont, d'une part, une même durée et, d'autre part, un point de départ situé à une même distance temporelle du début de la période D1 ou D2 (notamment, x minutes après le début de chaque période D1 ou D2). Telle que décrite précédemment, l'invention propose d'imposer des 5 débits du gaz distincts dans le local sur au moins deux périodes de temps Di, successives et de mesurer pour chaque période de temps Di, l'évolution temporelle de la concentration du gaz à l'intérieur du local cac(t). En variante, il est également possible d'imposer des concentrations du gaz ca, distinctes à l'intérieur du local sur au moins deux périodes de temps 10 Di, successives et de mesurer pour chaque période de temps Di, l'évolution temporelle du débit du gaz à l'intérieur du local ic(t). Selon cette variante, l'invention a pour objet un procédé de détermination du taux de renouvellement d'air ACH d'un local, caractérisé en ce qu'il comprend des étapes dans lesquelles : 15 - sur au moins deux périodes de temps Di, successives correspondant à des concentrations d'un gaz donné ca, distinctes appliquées dans le local, on procède à une campagne de mesures du débit du gaz à l'intérieur du local à intervalles de temps rapprochés, ainsi qu'à la détermination de la concentration du gaz à l'extérieur du local ceic à intervalles de temps 20 rapprochés ; - on détermine la valeur du taux de renouvellement d'air ACH du local en faisant converger : d'une part, un modèle diffusif exprimant la variation temporelle du débit du gaz à l'intérieur du local en fonction de la concentration du gaz à 25 l'extérieur du local ceic et de paramètres physiques du local à partir desquels le taux de renouvellement d'air ACH du local est calculable, et d'autre part, l'évolution mesurée k(t) du débit du gaz à l'intérieur du local en fonction du temps. 3049045 19 De préférence, le procédé est mis en oeuvre avec deux périodes de temps D1 et D2 successives correspondant à deux consignes de concentration du gaz cil et cil distinctes appliquées dans le local. Dans un mode de réalisation, au moins une partie des étapes du 5 procédé de détermination du taux de renouvellement d'air ACH du local sont déterminées par des instructions de programmes d'ordinateurs. En conséquence, l'invention a également pour objet un programme d'ordinateur sur un support d'enregistrement, ce programme étant susceptible d'être mis en oeuvre dans un terminal, ou plus généralement 10 dans un ordinateur, ce programme comportant des instructions adaptées à la mise en oeuvre de tout ou partie des étapes d'un procédé tel que décrit ci-dessus. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire 15 entre code source et code objet, tel que dans une forme partiellement compilée. L'invention a aussi pour objet un support d'enregistrement lisible par un ordinateur, et comportant des instructions d'un programme d'ordinateur tel que mentionné ci-dessus. 20 Le support d'enregistrement peut être n'importe quelle entité ou dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une mémoire morte, une mémoire non volatile réinscriptible, par exemple une clé USB, une carte SD, une EEPROM, ou encore un moyen d'enregistrement magnétique, par exemple 25 un disque dur. Le support d'enregistrement peut aussi être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé. Le support d'enregistrement peut être un support transmissible tel 30 qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet. 3049045 20 Un autre objet de l'invention est un dispositif de détermination du taux de renouvellement d'air ACH d'un local, ce dispositif comprenant : - au moins un appareil configuré pour appliquer, sur au moins deux périodes de temps Di, successives, des débits distincts d'un gaz donné 5 dans le local ; - au moins un capteur configuré pour mesurer une concentration du gaz à l'intérieur du local ca, à intervalles de temps rapprochés ; - un terminal comportant un module de traitement configuré pour faire converger, d'une part, un modèle diffusif exprimant la variation temporelle de 10 la concentration du gaz à l'intérieur du local ca, en fonction de la concentration du gaz à l'extérieur du local ceic et de paramètres physiques du local à partir desquels le taux de renouvellement d'air ACH du local est calculable et, d'autre part, l'évolution mesurée cac(t) de la concentration du gaz à l'intérieur du local ca, en fonction du temps, de manière à obtenir la 15 valeur du taux de renouvellement d'air du local. Chaque débit du gaz appliqué dans le local sur une période de temps Di, à l'aide du ou des appareils d'application d'un débit du gaz peut être un débit positif, correspondant à une injection du gaz dans le local, un débit négatif, correspondant à une extraction du gaz hors du local, ou un 20 débit nul. Selon un aspect de l'invention, lorsque le gaz utilisé est la vapeur d'eau, le dispositif peut comprendre, pour appliquer chaque débit de vapeur d'eau dans le local, un ou plusieurs humidificateurs, notamment de type humidificateur à ultrasons, et pour mesurer la concentration de vapeur d'eau 25 dans l'air à l'intérieur du local, un ou plusieurs capteurs d'humidité destinés à être positionnés dans le volume d'air à l'intérieur du local. Selon une caractéristique, le dispositif comprend également au moins un capteur configuré pour mesurer une concentration du gaz à l'extérieur du local ceic à intervalles de temps rapprochés. 3049045 21 Selon une autre caractéristique, le dispositif comprend en outre au moins un capteur configuré pour mesurer la température de l'air à l'intérieur du local Tac. L'invention a également pour objet un terminal comportant un module 5 de traitement configuré pour faire converger, d'une part, un modèle diffusif exprimant la variation temporelle de la concentration d'un gaz donné à l'intérieur d'un local ca, en fonction de la concentration du gaz à l'extérieur du local ceic et de paramètres physiques du local à partir desquels le taux de renouvellement d'air ACH du local est calculable et, d'autre part, une 10 évolution mesurée cac(t) de la concentration du gaz à l'intérieur du local ca, en fonction du temps, de manière à obtenir la valeur du taux de renouvellement d'air du local. Selon un aspect, le module de traitement du terminal selon l'invention comporte un programme d'ordinateur tel que mentionné ci-dessus, ce 15 programme étant enregistré sur un support d'enregistrement conforme à l'invention et constitué par une mémoire non volatile réinscriptible du terminal, les instructions du programme étant interprétables par un processeur du terminal. Le terminal, le programme d'ordinateur et le support d'enregistrement 20 présentent, selon l'invention, les mêmes caractéristiques que le procédé selon l'invention. L'invention peut être mise en oeuvre avec tout type de terminal, par exemple un ordinateur portable ou non, une tablette, un téléphone intelligent ("smartphone" en anglais). Dans le cadre de l'invention, le ou chaque capteur de mesure d'une 25 concentration du gaz à l'intérieur du local ca, peut être un capteur indépendant du terminal. En variante, le ou chaque capteur de mesure d'une concentration du gaz à l'intérieur du local ca, peut être un capteur intégré dans le terminal. En particulier, il est possible selon l'invention d'effectuer les mesures à l'aide de capteurs de gaz et éventuellement de capteurs de 30 température qui sont intégrés dans une tablette ou un smartphone, et d'avoir 3049045 22 une application installée sur la tablette ou le smartphone pour l'acquisition des mesures et la mise en oeuvre des étapes du procédé. Dans un mode de réalisation, le dispositif selon l'invention comprend des moyens de liaison, notamment sans fil tels que Bluetooth, 5 WiFi, etc., entre le ou chaque capteur de mesure d'une concentration du gaz et le terminal. De manière avantageuse, le terminal comporte des moyens de commande du ou de chaque appareil d'application d'un débit dudit gaz dans le local. 10 Les caractéristiques et avantages de l'invention apparaîtront dans la description qui va suivre de deux exemples de mise en oeuvre d'un procédé et d'un dispositif selon l'invention, donnée uniquement à titre d'exemple et faite en se référant aux figures annexées dans lesquelles : - la figure 1 est une vue schématique d'un local dont on souhaite 15 déterminer le taux de renouvellement d'air ACH conformément au procédé de l'invention, en utilisant un dispositif conforme à l'invention comprenant au moins un humidificateur à débit contrôlé pour injecter de la vapeur d'eau dans le local, au moins un capteur d'humidité et un terminal ; - la figure 2 est un graphe montrant l'évolution de l'humidité absolue 20 ca, à l'intérieur d'un bungalow en fonction du temps, au cours de la mise en oeuvre du procédé de l'invention comprenant une première période de temps D1 sur laquelle un premier débit de vapeur d'eau positif est appliqué dans le bungalow, suivie d'une deuxième période de temps D2 sur laquelle un deuxième débit de vapeur d'eau sensiblement nul est appliqué dans le 25 bungalow de manière à laisser évoluer librement l'humidité absolue ca, à l'intérieur du bungalow, l'évolution de l'humidité absolue ceic à l'extérieur du bungalow en fonction du temps étant également montrée sur cette figure ; - la figure 3 est un graphe montrant l'évolution de l'humidité absolue ca, à l'intérieur d'un appartement en fonction du temps, au cours de la mise 30 en oeuvre du procédé de l'invention comprenant une première période de temps D1 sur laquelle un premier débit de vapeur d'eau positif est appliqué 3049045 23 dans l'appartement, suivie d'une deuxième période de temps D2 sur laquelle un deuxième débit de vapeur d'eau sensiblement nul est appliqué dans l'appartement de manière à laisser évoluer librement l'humidité absolue ca, à l'intérieur de l'appartement, l'évolution de l'humidité absolue ceic à l'extérieur 5 de l'appartement en fonction du temps étant également montrée sur cette figure ; - la figure 4 est un schéma d'un modèle dit "3R2C" d'un local, avec trois résistances et deux capacités. 10 PROTOCOLE EXPERIMENTAL Le procédé selon l'invention est mis en oeuvre pour la détermination du taux de renouvellement d'air ACH d'un bungalow (exemple 1) et d'un appartement (exemple 2), qui constituent chacun un local au sens de l'invention, référencé 1 sur la figure 1. Le protocole expérimental est similaire 15 pour les deux exemples, comme décrit ci-dessous. Pour chaque exemple, le procédé comprend une première phase d'humidification sur une première période de temps D1, ayant une durée d'une heure pour le bungalow, de minuit à 1:00 du matin, et de deux heures pour l'appartement, de minuit à 2:00 du matin, au cours de laquelle un 20 premier débit de vapeur d'eau strictement positif est appliqué dans le local 1, ce qui correspond à une injection de vapeur d'eau dans le local. La première phase est suivie d'une deuxième phase de "déshumidification naturelle" sur une deuxième période de temps D2, ayant une durée d'une heure pour le bungalow, de 1:00 à 2:00 du matin, et de deux heures pour 25 l'appartement, de 2:00 à 4:00 du matin, au cours de laquelle un deuxième débit de vapeur d'eau nul est appliqué dans le local 1, donnant lieu à une diminution libre de la concentration de vapeur d'eau à l'intérieur du local. Dans chaque exemple, la phase de déshumidification naturelle a la même durée que la phase d'humidification. 30 L'application des débits de vapeur d'eau et dans le local 1 est assurée par le biais d'un ou plusieurs humidificateurs 20 de type "GOTA", 3049045 24 commercialisés par la société Air Naturel, placés dans le local. Chaque humidificateur 20 est un humidificateur à ultrasons avec un réservoir de 3 L, pouvant humidifier jusqu'à un taux de 300 g.h-1. Chaque humidificateur 20 est branché sur une prise programmable afin de pouvoir être commuté en 5 ON/OFF de manière automatique. La programmation des prises est effectuée de manière à démarrer le ou chaque humidificateur 20 au début de la première période de temps D1, et à couper le ou chaque humidificateur 20 à la fin de la première période de temps D1. Afin de connaître le débit de vapeur d'eau appliqué dans le local au cours de la première période de 10 temps D1, chaque humidificateur 20 est pesé au début et à la fin de la période D1. Connaissant la durée d'allumage de chaque humidificateur 20, il est possible de déterminer le débit moyen de vapeur d'eau appliqué au cours de la période D1. Au cours des deux périodes de temps D1 et D2, l'évolution de 15 l'humidité et de la température à l'intérieur et à l'extérieur du local 1 est mesurée à l'aide de capteurs 40 d'humidité et de température de type SHT15, commercialisés par la société Sensirion. Pour le bungalow (exemple 1), le nombre de capteurs 40 utilisés est un capteur 40 placé dans l'air à l'intérieur du bungalow et un capteur 40 placé dans l'air à l'extérieur du 20 bungalow. Pour l'appartement (exemple 2), le nombre de capteurs 40 utilisés est cinq capteurs 40 répartis en étant placés dans l'air à l'intérieur de l'appartement et un capteur 40 placé dans l'air à l'extérieur de l'appartement. Le positionnement des capteurs 40 dans le local 1 est ajusté de manière à obtenir une mesure représentative de l'humidité moyenne et de la 25 température moyenne de l'air à l'intérieur du local. Pour chaque exemple, le procédé est mis en oeuvre alors que le local 1 est inoccupé et toutes les sources de production de vapeur d'eau autres que les humidificateurs 20 sont coupées dans le local. Aucune ventilation mécanique ni source de chaleur ou de froid n'est active dans le local 1 au 30 cours du procédé. De plus, dans le cas de l'appartement (exemple 2), toutes les portes de communication internes sont ouvertes. Dans chaque exemple, 3049045 25 des ventilateurs domestiques sont installés dans le local 1, afin d'assurer un léger brassage de l'air à l'intérieur du local. Les données de mesures brutes obtenues dans le local 1 sont ensuite acquises par un système d'acquisition d'un terminal 10, qui les traite en 5 exécutant des instructions d'un programme d'ordinateur PG conforme à l'invention installé dans le terminal 10, de manière à déterminer le taux de renouvellement d'air du local 1. Comme montré de façon schématique sur la figure 1, le terminal 10 est dans cet exemple constitué par une tablette ou un téléphone intelligent 10 ("smartphone"), offrant à l'utilisateur, au-delà des fonctions de communication, l'accès à diverses applications dès lors que celles-ci ont été installées dans le terminal. D'un point de vue matériel, le terminal 10 comporte notamment un processeur 11, une mémoire morte de type ROM 12 dans laquelle sont 15 enregistrées des fonctions systèmes, en particulier des pilotes logiciel ("drivers" en anglais) et le système d'exploitation du terminal, un écran 15, un ou plusieurs modules de communication (3G, 4G, Bluetooth, WiFi, ...) 17 et une mémoire non volatile réinscriptible 18 comportant des applications APP et des données utilisateurs non représentées sur la figure, ces éléments 20 étant reliés entre eux par un système de bus. De façon connue, l'écran 15 constitue une interface homme-machine tactile sur laquelle sont représentées des icônes 11, 12, IT, correspondant aux applications système et aux différentes applications APP installées par l'utilisateur du terminal. 25 Parmi ces icônes, une icône IT permet au terminal d'accéder à distance, via un réseau de télécommunications, à un portail d'applications téléchargeables compatibles avec le système d'exploitation du terminal et d'installer, éventuellement moyennant paiement et/ou authentification, de nouvelles applications APP dans la mémoire non volatile réinscriptible 18. 30 Dans le mode de réalisation décrit ici, le programme d'ordinateur PG conforme à l'invention peut être téléchargé à partir de ce portail d'applications, et une icône associée présentée sur l'interface tactile 15. 3049045 26 En particulier, le terminal 10 sur lequel est installé le programme d'ordinateur PG comporte un module de traitement qui est configuré pour mettre en oeuvre les étapes suivantes. - Calcul de l'humidité absolue cik à l'intérieur du local 5 L'évolution de l'humidité absolue (ou concentration de vapeur d'eau) Ca, à l'intérieur du local 1 est déterminée à partir des humidités relatives et températures mesurées dans le volume intérieur du local. Une moyenne des mesures des différents capteurs 40 répartis dans le local 1 est calculée, pondérée par le volume d'air représentatif de chaque point de mesure. 10 L'humidité absolue ca, à l'intérieur du local 1 est alors déterminée en calculant tout d'abord la pression de vapeur saturante, notamment par l'équation : Pws T, 1n- = - (C1i9 + C291.5 + C3193 + C4193.5 + C5194 + C5194 + C6197.5) (3) T avec 19 = 1 - -T, où T est la température en K, Pws la pression de vapeur Tc saturante en hPa, Tc. = 647,096 K la température critique, /3, = 220640 hPa la 15 pression critique, et C. des coefficients sans unité donnés dans le Tableau 1 ci-dessous. Cette relation donne une bonne précision pour des températures allant de 0°C à 373°C. Bien entendu, d'autres équations que l'équation (3) ci-dessus sont également utilisables pour le calcul de l'humidité absolue. C1 -7,85951783 C2 1,84408259 C3 -11,7866497 C4 22,6807411 C5 -15,9618719 C6 1,80122502 20 Tableau 1 3049045 27 La pression de vapeur est alors calculée comme étant Pw = Pws.RH, où RH est l'humidité relative de l'air, et l'humidité absolue comme étant Cik = C.141 où C = 2,16679 g. K/J est une constante. - Calcul du débit de vapeur d'eau 91 au cours de la période D1 5 Le débit de vapeur d'eau appliqué dans le local 1 au cours de la première période de temps D1 est calculé par différence entre la masse totale des humidificateurs 20 au début et à la fin de la période D1 et en connaissant la durée d'injection de vapeur d'eau : Mtot,fin Mtot,début tin] OÙ Mtot,fin et Mtot,début sont respectivement les masses totales des 10 humidificateurs dans le local au début et à la fin de la période D1 en grammes et tin] est la durée d'injection de vapeur d'eau en heures. - Calcul de l'humidité absolue cek à l'extérieur du local L'évolution de l'humidité absolue (ou concentration de vapeur d'eau) Cek à l'extérieur du local 1 est déterminée à partir des humidités relatives et 15 températures mesurées à l'extérieur du local, en utilisant la même équation (3) que pour le calcul de l'humidité absolue cik à l'intérieur du local. - Calcul du taux de renouvellement d'air ACH et du volume effectif V Le taux de renouvellement d'air ACH et le volume effectif V du local 1 sont calculés en utilisant, respectivement, l'équation (1) et l'équation (2) 20 précédentes. Dans la description qui suit, on détaille les étapes de calcul du taux de renouvellement d'air ACH',, et du volume effectif V',, de chaque local 1 en faisant référence au graphe de la figure 2 ou au graphe de la figure 3, qui montrent l'évolution de l'humidité absolue cik à l'intérieur du local en fonction 25 du temps. Ces graphes illustrent les grandeurs intervenant dans les équations (1) et (2), à savoir les pentes a1 et a2 de l'humidité absolue à l'intérieur du local cil et ci2 sur chaque intervalle de temps At, ou At2 sélectionné sur la période D1 ou D2, et les différences d'humidité absolue = 3049045 28 &lm et Ac27, entre l'intérieur et l'extérieur du local sur chaque intervalle de temps At, ou Ott. Même si les étapes de calcul détaillées ci-après peuvent être effectuées manuellement, elles sont de préférence effectuées de manière 5 automatisée par un terminal tel que le terminal 10 décrit précédemment qui peut, en particulier, ne pas recourir à une représentation graphique de l'évolution cik(t). Le terminal requiert auprès de l'utilisateur des données d'entrée, qui peuvent être le volume du local, les débits du gaz appliqués sur les 10 périodes D1, D2 (ou, en variante, la valeur du paramètre a et une valeur de référence du taux de renouvellement d'air ACH'f). Le terminal fournit en sortie les valeurs calculées du taux de renouvellement d'air ACHcaIc et du volume effectif Vcaic du local, ainsi que d'éventuels autres paramètres tels que ceux indiqués dans les tableaux 2 et 15 3 ci-dessous. Dans le cadre de l'invention, il peut être intéressant de stocker les données de mesures obtenues dans le local. Ces données peuvent par exemple être utilisées comme un historique dans le cadre d'une rénovation. De plus, ces données peuvent être réutilisées pour optimiser les valeurs 20 calculées du taux de renouvellement d'air ACHcaIc et du volume effectif Vcaic du local, en modifiant des paramètres du traitement pour un même set de mesures, comme illustré par les étapes de correction dans les exemples ci-dessous. Le stockage des données de mesures peut être réalisé dans la mémoire d'un terminal utilisé pour la mise en oeuvre de l'invention ou dans 25 une mémoire externe de tout type adapté à cette fonction de stockage. EXEMPLE 1 : BUNGALOW La figure 2 montre les résultats obtenus en mettant en oeuvre le procédé dans un bungalow ayant une surface au sol de 13 m2, un volume de 30 33 m3 et une surface totale d'enveloppe de 70 m2. La paroi externe du bungalow est constituée de panneaux sandwich isolants comprenant une 3049045 29 couche de polyuréthane d'épaisseur 4 cm insérée entre deux plaques de métal, d'une porte et de deux fenêtres triple vitrage. Une isolation supplémentaire a été ajoutée à l'enveloppe, qui comprend les matériaux suivants : 5 - des membranes de type STOPVAP, commercialisées par la société Saint-Gobain Isover, - des panneaux VIP (Vacuum Insulated Panels) de type va-Q-vip F, commercialisés par la société va-Q-tec, qui recouvrent les murs ; - 3 cm de polystyrène expansé pour le sol et le plafond, le sol étant 10 aussi recouvert d'un panneau de bois à lamelles minces orientées (OSB). L'enveloppe du bungalow présente un coefficient de transmission thermique UBAT d'environ 0,6 W/m2K. Le coefficient n50 du bungalow, obtenu par un test de porte soufflante, 15 est de 8,6 h-1. En utilisant le modèle de Persily-Kronvall, cela correspond à une valeur de référence du taux de renouvellement d'air ACH'f. = 15 m3.11-1. Dans la phase d'humidification du bungalow, sur la première période de temps D1 de minuit à 1:00 du matin, on applique à l'aide des humidificateurs 20 présents dans le bungalow le premier débit de vapeur 20 d'eau strictement positif, qui est choisi tel que le paramètre a = 1 - ACHref3Ci(0) est égal à 0,46. Dans cet exemple, la valeur de référence ACH'f est de 15 m3.11-1 et la différence d'humidité absolue initiale Ac1(0) est de 6,4 g.m-3, ce qui correspond à une valeur du premier débit de vapeur d'eau de l'ordre de 180 g.h-1. 25 La courbe représentative de l'évolution de l'humidité absolue c,' à l'intérieur du bungalow en fonction du temps pendant la première période de temps D1 est montrée sur la figure 2. Comme visible sur cette figure, la courbe d'augmentation de l'humidité absolue à l'intérieur du bungalow présente une partie sensiblement linéaire sur l'intervalle de temps At, de 30 0:45 à 1:00 du matin. La mise en équation de cette partie linéaire de la courbe donne : c11 = 17,8 g. M-3 + 0,00630 t, avec t en secondes. 3049045 30 La figure 2 fait également apparaître l'évolution de l'humidité absolue cet à l'extérieur du bungalow en fonction du temps pendant la première période de temps D1. L'humidité absolue ce' à l'extérieur du bungalow sur l'intervalle de temps At, est suffisamment stable pour qu'on puisse la 5 considérer sensiblement constante et égale à l'humidité absolue moyenne sur l'intervalle de temps At', à savoir dans cet exemple celin = 11,0 g.m-3. Dans la phase de "déshumidification naturelle" du bungalow, sur la deuxième période de temps D2 de 1:00 à 2:00 du matin, les humidificateurs 20 présents dans le bungalow sont éteints de manière à appliquer le 10 deuxième débit de vapeur d'eau nul. La figure 2 montre la courbe représentative de l'évolution de l'humidité absolue cil à l'intérieur du bungalow en fonction du temps pendant la deuxième période de temps D2. Comme visible sur cette figure, la courbe de diminution de l'humidité absolue à l'intérieur du bungalow présente une partie 15 sensiblement linéaire sur l'intervalle de temps At2 de 1:45 à 2:00 du matin. L'intervalle de temps At2 de traitement des mesures sur la deuxième période D2 est choisi de manière à présenter une "symétrie" avec l'intervalle de temps At, de traitement des mesures sur la première période D1, c'est-à-dire de sorte que les deux intervalles de temps At, et At2 ont, d'une part, une 20 même durée de 15 minutes et, d'autre part, un point de départ situé, pour chaque intervalle Atic, 45 minutes après le début de la période Di,. La mise en équation de la partie linéaire de la courbe sur l'intervalle de temps At2 donne : ci2 = 19 g.m-3- 0,00311 t, avec t en secondes. L'évolution de l'humidité absolue ce2 à l'extérieur du bungalow en 25 fonction du temps pendant la deuxième période de temps D2 est également montrée sur la figure 2. Comme dans la première étape, l'humidité absolue Ce2 à l'extérieur du bungalow sur l'intervalle de temps Ott est suffisamment stable pour qu'on puisse la considérer sensiblement constante et égale à l'humidité absolue moyenne sur l'intervalle de temps At2, à savoir dans cet 30 exemple ce27, = 11,0 g.m-3. 3049045 31 aixii2-a2xCh Comme ACH = d'après l'équation (1) précédente, en a1xAc2m-a2xAcim prenant a1 = 22,7 g.m-3.h-1, a2 = - 11,2 g.m-3.h-1, Aclm = 7,7 g.m-3, Oc2m= 7,0 g.rn = 180 g.h-1, = 0 g.h-1, on obtient la valeur du taux de renouvellement d'air ACH du bungalow : 5 ACHcaic = 8,2 m3.11-1. Cette valeur ACHcaIc du taux renouvellement d'air obtenue selon l'invention présente un écart significatif par rapport à l'estimation issue du test de porte soufflante ACH'f. = 15 m3.11-1, mais l'ordre de grandeur reste cohérent. En effet, pendant la mesure de type "porte soufflante", le 10 renouvellement d'air se produit grâce à un forçage à différence de pression très importante, qui n'est pas rencontré de manière naturelle, et il faut utiliser un modèle d'extrapolation à basse pression, ici le modèle de PersilyKronvall, pour remonter au paramètre physique réaliste. Comme V = X'n'C2M el2 X'n'Clin d'après l'équation (2) précédente, on a1xAc2m-a2xAcim 15 obtient également la valeur du volume effectif du bungalow : Vcaic = 122,9 m3. On peut observer que le volume effectif déterminé Vau,. = 122,9 m3 est bien supérieur au volume intérieur réel du bungalow Vrée, = 33 m3. Comme expliqué précédemment, cela peut être lié à l'absorption de la vapeur d'eau 20 par les matériaux présents dans le bungalow et au fait qu'une durée trop longue a été appliquée pour le traitement des mesures. Lorsqu'on obtient une telle valeur trop élevée du volume effectif Vcaic, l'invention propose de procéder à une correction a posteriori de la valeur calculée ACHcaIc du taux de renouvellement d'air, en réduisant la durée 25 appliquée pour le traitement des mesures. En pratique, on procède à cette correction en déplaçant les intervalles de temps At, et At2 dans chaque période de temps D1 et D2, tout en conservant la "symétrie" des intervalles de temps At, et At2 telle que définie précédemment, c'est-à-dire en choisissant pour les deux intervalles de temps At, et At2, d'une part, une même durée et, 30 d'autre part, un point de départ situé à une même distance temporelle du 3049045 32 début de la période D1 ou D2 (notamment, x minutes après le début de chaque période D1 ou D2), jusqu'à obtenir une valeur du volume effectif Vcaic sensiblement égale au volume intérieur réel. Avec une telle démarche, on obtient pour le bungalow des intervalles 5 de temps "corrigés" At,' de minuit à 0:15 du matin et At2' de 1:00 à 1:15 du matin tels que montrés sur la figure 2, qui correspondent à une durée réduite pour le traitement des mesures, ce qui donne les résultats corrigés suivants : ACHccac' = 9,6 m3.11-1 ; Vcaic' = 56,6 m3. 10 Dans le cas du bungalow, le volume effectif corrigé Valu' = 56,6 m3 est toujours supérieur au volume intérieur réel du bungalow Vréei = 33 m3. Cela est dû au fait que les revêtements intérieurs du bungalow ne sont pas finalisés, ce qui conduit à une absorption très rapide de la vapeur d'eau par les matériaux de construction du bungalow. Pour un tel bungalow, il serait 15 préférable de mettre en oeuvre l'invention avec un autre gaz que la vapeur d'eau pour gagner en précision sur la valeur du volume effectif et du taux de renouvellement d'air. Il ressort des résultats corrigés précédents que le renouvellement d'air dans le bungalow est d'environ 0,29 volume par heure, ce qui est insuffisant 20 pour assurer une bonne qualité de l'air intérieur. La perte thermique du bungalow liée à l'infiltration est telle que = p.Cp.V.ACH, où p est la masse volumique de l'air sec à 20°C et Cp est la chaleur spécifique de l'air sec à 20°C. A partir de la valeur du taux de renouvellement d'air ACliccac' = 9,6 m3.11-1, on obtient donc une valeur de la 25 perte thermique liée à l'infiltration Kin,. = 3,3 W.K-1. Une évaluation du coefficient de déperdition thermique totale du bungalow a aussi été réalisée, qui fournit une valeur Kt0t = 29 W.K-1. Ainsi, la perte thermique du bungalow liée à l'infiltration représente environ 11% de la perte thermique totale du bungalow, pour un bilan annuel de l'ordre de 13 30 kWh. m-2. an-1. 3049045 33 Les données d'infiltration et de déperdition thermique pour le bungalow sont résumées dans le Tableau 2 ci-dessous. Renouvellement d'air - Procédé de l'invention Renouvellement d'air - Porte soufflante Déperdition thermique V ACHcau. V cale (m3) (m3. h-1) (m3) 33 15 ACH K (h-1) (W. K-1) 0,43 4,9 ssol (kWh. m-2. an-1) 13 20 56,6 0,29 29 118 Tableau 2 (bungalow) 5 EXEMPLE 2 : APPARTEMENT La figure 3 montre les résultats obtenus en mettant en oeuvre le 10 procédé dans un appartement situé dans une copropriété ancienne (année de construction : 1879) situé à Levallois-Perret, France. L'appartement n'est pas isolé, il est équipé de doubles vitrages récents mais présente une mauvaise étanchéité à l'air. La surface au sol est de 54 m2, le volume intérieur de 151 m3, et l'appartement possède deux faces déperditives 15 d'environ 47 m2. L'enveloppe de l'appartement présente un coefficient de transmission thermique UBAT d'environ 1,9 W/m2K. Le coefficient n50 de l'appartement, obtenu par un test de porte soufflante, est de 7,3 h-1. En utilisant le modèle de Persily-Kronvall, cela 20 correspond à une valeur de référence du taux de renouvellement d'air ACH'f. = 55 m3.11-1. Dans la phase d'humidification de l'appartement, sur la première période de temps D1 de minuit à 2:00 du matin, on applique à l'aide des 3049045 34 humidificateurs 20 répartis dans l'appartement le premier débit de vapeur d'eau strictement positif, qui est choisi tel que le paramètre a = 1 - ACH re f .ACi (0) est égal à 0,78. Dans cet exemple, la valeur de référence ACH'f est de 55 m3.11-1 et la différence d'humidité absolue initiale Ac1(0) est de 2,2 5 g.m-3, ce qui correspond à une valeur du premier débit de vapeur d'eau de l'ordre de 550 g.h-1. La courbe représentative de l'évolution de l'humidité absolue c,' à l'intérieur de l'appartement en fonction du temps pendant la première période de temps D1 est montrée sur la figure 3. Comme visible sur cette figure, la 10 courbe d'augmentation de l'humidité absolue à l'intérieur de l'appartement présente une partie sensiblement linéaire sur l'intervalle de temps At, de 1:30 à 2:00 du matin. La mise en équation de cette partie linéaire de la courbe donne : c11 = 7,4 g. M-3 + 0,00578 t, avec t en secondes. La figure 3 fait également apparaître l'évolution de l'humidité absolue 15 cet à l'extérieur de l'appartement en fonction du temps pendant la première période de temps D1. L'humidité absolue ce, à l'extérieur de l'appartement sur l'intervalle de temps At, est suffisamment stable pour qu'on puisse la considérer sensiblement constante et égale à l'humidité absolue moyenne sur l'intervalle de temps At', à savoir dans cet exemple celin = 4,3 g.m-3. 20 Dans la phase de "déshumidification naturelle" de l'appartement, sur la deuxième période de temps D2 de 2:00 à 4:00 du matin, les humidificateurs 20 présents dans l'appartement sont éteints de manière à appliquer le deuxième débit de vapeur d'eau nul. La figure 3 montre la courbe représentative de l'évolution de l'humidité 25 absolue cil à l'intérieur de l'appartement en fonction du temps pendant la deuxième période de temps D2. Comme visible sur cette figure, la courbe de diminution de l'humidité absolue à l'intérieur de l'appartement présente une partie sensiblement linéaire sur l'intervalle de temps At2 de 3:30 à 4:00 du matin. Comme pour l'exemple précédent, l'intervalle de temps At2 de 30 traitement des mesures sur la deuxième période D2 est choisi de manière à présenter une "symétrie" avec l'intervalle de temps At, de traitement des 3049045 mesures sur la première période D1, c'est-à-dire de sorte que les deux intervalles de temps At, et Ott ont, d'une part, une même durée de 30 minutes et, d'autre part, un point de départ situé, pour chaque intervalle Atic, 1h30 après le début de la période Di,. La mise en équation de la partie 5 linéaire de la courbe sur l'intervalle de temps Ott donne : ci2 = 9,0 g.m-3 - 0,00334 t, avec t en secondes. L'évolution de l'humidité absolue cet à l'extérieur de l'appartement en fonction du temps pendant la deuxième période de temps D2 est également montrée sur la figure 3. Comme dans la première étape, l'humidité absolue 10 Cet à l'extérieur de l'appartement sur l'intervalle de temps Ott est suffisamment stable pour qu'on puisse la considérer sensiblement constante et égale à l'humidité absolue moyenne sur l'intervalle de temps Ott, à savoir dans cet exemple ce27, = 4,2 g.m-3. aixii2-a2xCh Comme ACH = d'après l'équation (1) précédente, en a1xAc2m-a2xAcim 15 prenant a1 = 20,8 g.m-3.h-1, a2 = - 12,0 g.m-3.h-1, Aclm = 4,8 g.m-3, AC2m= 2,9 g.m-3, = 550 g.1-1-1, = 0 g.1-1-1, on obtient la valeur du taux de renouvellement d'air ACH de l'appartement : ACHcaic = 56,1 m3.1-1-1. Cette valeur ACHcaIc du taux renouvellement d'air obtenue selon 20 l'invention est en accord avec l'estimation issue du test de porte soufflante ACHref = 55 m3.11-1. Comme V = eilx°c2mei2x°clin d'après l'équation (2) précédente, on a1xAc2m-a2xAcim obtient également la valeur du volume effectif de l'appartement : Vccac = 321,7 m3. 25 On peut observer que le volume effectif déterminé Veau. = 321,7 m3 est bien supérieur au volume intérieur réel de l'appartement Vréei = 151 m3, ce qui peut être lié à l'absorption de la vapeur d'eau par les matériaux présents dans l'appartement et au fait qu'une durée trop longue a été appliquée pour le traitement des mesures. 3049045 36 Comme dans l'exemple du bungalow, on peut alors procéder à une correction a posteriori de la valeur calculée ACHcaIc du taux de renouvellement d'air de l'appartement, en déplaçant les intervalles de temps At, et At2 dans chaque période de temps D1 et D2, tout en conservant la 5 "symétrie" de ces intervalles de temps At, et At2, de manière à réduire la durée de traitement des mesures, jusqu'à obtenir une valeur du volume effectif Vcaic sensiblement égale au volume intérieur réel de l'appartement. Avec une telle démarche, on obtient pour l'appartement des intervalles de temps "corrigés" At,' de minuit à 0:30 du matin et At2' de 2:00 à 2:30 du 10 matin tels que montrés sur la figure 3, qui correspondent à une durée réduite pour le traitement des mesures, ce qui donne les résultats corrigés suivants : ACHccac' = 73,3 m3.1-1-1 ; Vcaic' = 159,1 m3. La valeur corrigée du volume effectif Valu' = 159,1 m3 est en bon 15 accord avec le volume intérieur réel de l'appartement Vréei = 151 m3, ce qui tend à montrer que la valeur corrigée du taux de renouvellement d'air ACliccac' = 73,3 m3.11-1 est une meilleure estimation que la première valeur ACHccac = 56,1 m3.1-1-1. Il ressort des résultats ci-dessus que le renouvellement d'air dans 20 l'appartement est d'environ 0,48 volume par heure, ce qui est presque suffisant pour assurer une bonne qualité de l'air intérieur. La perte thermique de l'appartement liée à l'infiltration est telle que Kin,. = p. Cp.V.ACH, où p est la masse volumique de l'air et Cp est la chaleur spécifique de l'air sec à 20°C. A partir de la valeur du taux de 25 renouvellement d'air ACHccac = 73,3 m3.1-1-1, on obtient donc une valeur de la perte thermique liée à l'infiltration Kin,. = 24,4 W.K-1. Une évaluation du coefficient de déperdition thermique totale de l'appartement a aussi été réalisée, qui fournit une valeur Kt0t = 81 W.K-1. Ainsi, la perte thermique de l'appartement liée à l'infiltration représente 30 environ 30% de la perte thermique totale de l'appartement, pour un bilan annuel de l'ordre de 25 kWh.m-2.an-1. 3049045 37 Les données d'infiltration et de déperdition thermique pour l'appartement sont résumées dans le Tableau 3 ci-dessous. Renouvellement V ACH'u. V cale ACH K ssol d'air - Procédé de (m3) (m3. h-1) (m3) (h-1) (W. K-1) (kWh. m-2. an-1) l'invention 151 73,3 159,1 0,48 24,4 25 55 0,36 19 19 81 83 Renouvellement d'air - Porte soufflante Déperdition thermique Tableau 3 (appartement) 5 De manière avantageuse, dans les exemples précédents, les étapes de sélection des intervalles de temps Atic pour le traitement des données, de 10 linéarisation, et de calcul du taux de renouvellement d'air ACH et du volume effectif V du local, sont réalisées au moyen d'un terminal tel que le terminal 10 décrit précédemment comportant le programme PG. Comme illustré dans les exemples précédents, le procédé de l'invention propose de réaliser des sollicitations du local en au moins deux 15 phases ayant des débits différents du gaz (qui dans les exemples est la vapeur d'eau). Cette expérience dynamique en au moins deux phases permet de réduire le temps de mesure, tout en conservant une bonne précision sur le résultat. Des applications du procédé et du dispositif selon l'invention 20 comprennent, notamment : - la prescription lors de la rénovation de logements anciens, en permettant d'évaluer la quantité d'énergie perdue à cause des 3049045 38 infiltrations naturelles et ainsi de prescrire une solution adaptée pour améliorer l'enveloppe du local ; - la validation de la perméabilité de l'enveloppe à réception d'un bâtiment neuf ; 5 - l'évaluation du taux de renouvellement d'air créé par une ventilation mécanique, afin de vérifier que le débit est suffisant pour assurer une bonne qualité de l'air intérieur. L'invention n'est pas limitée aux exemples décrits et représentés. En particulier, le procédé selon l'invention peut être mis en oeuvre 10 avec tout gaz adapté autre que la vapeur d'eau, notamment CO2, He, SF6, H2, N2, ou d'autres gaz traceurs réfrigérants. De plus, sur au moins une période de temps Dk, le débit du gaz appliqué dans le local peut être négatif au lieu d'être positif, correspondant à une extraction du gaz hors du local au lieu d'une injection du gaz dans le 15 local. Par ailleurs, dans les exemples précédents, la méthode de traitement des données correspond au cas où le modèle diffusif utilisé est un modèle R-C simple avec une résistance et une capacité et où on procède à une linéarisation sur des intervalles de temps Atic sélectionnés dans chaque 20 période de temps Dk. En variante, les données de l'évolution de la concentration du gaz à l'intérieur du local ca, en fonction du temps peuvent être traitées différemment, par exemple en utilisant un modèle récursif sous forme discrète en temps avec une période d'échantillonnage donnée décrivant le comportement en régime transitoire du local, ou encore avec un 25 autre modèle qu'un modèle R-C simple, par exemple un modèle "3R2C" du local avec trois résistances et deux capacités. Le modèle "3R2C" est alors ajusté ("fitting") sur l'évolution de la concentration du gaz à l'intérieur du local en fonction du temps. La figure 4 montre un schéma d'un modèle "3R2C" d'un local avec 30 trois résistances et deux capacités. Dans ce modèle "3R2C", on considère un réseau à trois noeuds, qui sont l'air à l'intérieur du local (I), l'air à l'extérieur du local (E) et les parois du local (W). L'ambiance extérieure est 3049045 39 considérée à une concentration du gaz constante imposée cE. Deux noeuds cw et cI représentent schématiquement la concentration du gaz dans les parois et dans l'air intérieur, qui ont chacun une inertie Cw, CI associée représentant la capacité de stockage du gaz dans les parois et dans l'air 5 intérieur. La résistance REW, placée entre l'ambiance extérieure et le noeud des parois, et la résistance R1w, placée entre l'ambiance intérieure et le noeud des parois, représentent la résistance à la diffusion du gaz à travers les parois. La troisième résistance RIE, placée entre l'ambiance intérieure et l'ambiance extérieure, représente la résistance à la diffusion du gaz par 10 infiltrations. Le taux de renouvellement d'air ACH du local est alors l'inverse de la résistance RIE. La résistance totale RT est telle que : RIE(RIW + REW) RT = RIE+Riw + REW En pratique, la résistance RIE est très inférieure à la somme des résistances REW et R1w car le transfert du gaz via les infiltrations est très 15 rapide, alors que la diffusion à travers des parois opaques est un processus lent, ce qui explique que le modèle R-C simple est une bonne approximation. Selon une variante, il est également possible, dans le cadre de l'invention, d'effectuer les mesures à l'aide de capteurs de gaz (notamment d'humidité) et de capteurs de température qui sont intégrés dans le terminal 20 utilisé pour l'acquisition des mesures et la mise en oeuvre des étapes de traitement des données. De manière avantageuse, de tels capteurs d'humidité et de température sont couramment intégrés dans les tablettes et smartphones. Selon une autre variante, le terminal peut également comporter des 25 moyens de commande du ou des appareils d'application d'un débit du gaz dans le local, et communiquer avec ceux-ci par des moyens de liaison sans fil tels que Bluetooth ou WiFi. REVENDICATIONS1. Procédé de détermination du taux de renouvellement d'air ACH d'un local, caractérisé en ce qu'il comprend des étapes dans lesquelles : - sur au moins deux périodes de temps Di, successives correspondant à des débits d'un gaz donné distincts appliqués dans le local, on procède à une campagne de mesures permettant de déterminer la concentration du gaz à l'intérieur du local ca, à intervalles de temps rapprochés, ainsi qu'à la détermination de la concentration du gaz à l'extérieur du local ceic à intervalles de temps rapprochés ; - on détermine la valeur du taux de renouvellement d'air ACH du local en faisant converger : d'une part, un modèle diffusif exprimant la variation temporelle de la concentration du gaz à l'intérieur du local ca, en fonction de la concentration du gaz à l'extérieur du local ceic et de paramètres physiques du local à partir desquels le taux de renouvellement d'air ACH du local est calculable, et d'autre part, l'évolution mesurée cac(t) de la concentration du gaz à l'intérieur du local ca, en fonction du temps. 2. Procédé selon la revendication 1, caractérisé en ce qu'il comprend des étapes dans lesquelles : - sur au moins deux périodes de temps Di, successives correspondant à des débits du gaz distincts appliqués dans le local, on procède à une campagne de mesures permettant de déterminer la concentration du gaz à l'intérieur du local ca, à intervalles de temps rapprochés, ainsi qu'à la détermination de la concentration du gaz à l'extérieur du local ceic à intervalles de temps rapprochés ; - pour chaque période de temps Dk, à partir de l'évolution mesurée cac(t) de la concentration du gaz à l'intérieur du local ca, en fonction du temps : o soit, s'il existe un intervalle de temps Atic pour lequel l'évolution cik(t) est sensiblement linéaire, on détermine sur cet intervalle de temps Atic 3049045 41 la pente ak de la tangente à l'évolution cik(t) et on déduit la valeur du taux de renouvellement d'air ACH du local à partir des pentes ak ; o soit, s'il n'existe pas d'intervalle de temps pour lequel l'évolution cik(t) est sensiblement linéaire, on sélectionne un intervalle de temps Atk' 5 sur lequel l'évolution cik(t) est sensiblement exponentielle de type exp(-t/ -c), avec 'r le temps au bout duquel le volume d'air intérieur du local a été renouvelé, et on déduit la valeur du taux de renouvellement d'air ACH du local, qui est la valeur telle que l'évolution Lnkt9k(t) - AcH)/(61k(0) - AeicH)1 est une droite, avec °k(t) = cik(t) - cekin: où Cekm' est la moyenne de la 10 concentration du gaz à l'extérieur du local cek sur l'intervalle de temps Atk'. 3. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il est mis en oeuvre avec deux périodes de temps D1 et D2 successives correspondant à deux consignes de débits du gaz et distinctes appliquées dans le local. 15 4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le gaz est H2O, CO2, He, SF6, H2, N2, ou un gaz réfrigérant. 5. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que, pour chaque période de temps Dk, le débit du gaz 20 appliqué dans le local comprend un débit Rimpk imposé au moyen d'au moins un appareil (20) à débit contrôlé. 6. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que les mesures permettant de déterminer la concentration du gaz à l'intérieur du local cik sont effectuées à l'aide d'un ou de plusieurs 25 capteurs (40) dudit gaz placés dans le volume intérieur du local. 7. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que, sur chaque période de temps Dk, la température à l'intérieur du local Tik est stable. 8. Procédé selon l'une quelconque des revendications précédentes, 30 caractérisé en ce que, sur chaque période de temps Dk, la concentration du gaz à l'extérieur du local cek est stable. 3049045 42 9. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que, sur chaque période de temps Dk, le rayonnement solaire est faible, de préférence nul. 10. Procédé selon l'une quelconque des revendications précédentes, 5 caractérisé en ce que le procédé est mis en oeuvre alors que le local est inoccupé. 11. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on vérifie que la valeur V',, du volume effectif du local calculée à partir du modèle diffusif et de l'évolution mesurée cac(t) 10 correspond au volume réel du local, notamment à 20% près. 12. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend des étapes dans lesquelles : - on procède, sur deux périodes de temps D1 et D2 successives : i. sur la première période de temps D1, à l'application dans le local d'un 15 premier débit du gaz et à une campagne de mesures permettant de déterminer la concentration du gaz à l'intérieur du local cil à intervalles de temps rapprochés, ainsi qu'à la détermination de la concentration du gaz à l'extérieur du local cet à intervalles de temps rapprochés, le premier débit du gaz étant tel que le paramètre ACHre f.ACi (0) 20 a = 1 est inférieur ou égal à 0,8, avec Ac1(0) = c11(0) - Q1 Cern, où t = 0 est le point de départ de la première période de temps D1, Cern est la concentration moyenne du gaz à l'extérieur du local sur l'ensemble des périodes de temps D1 et D2, et ACH'f est une valeur de référence du taux de renouvellement d'air du local, puis 25 ii. sur la deuxième période de temps D2, à l'application dans le local d'un deuxième débit du gaz sensiblement nul, et à une campagne de mesures permettant de déterminer la concentration du gaz à l'intérieur du local cil à intervalles de temps rapprochés, ainsi qu'à la détermination de la concentration du gaz à l'extérieur du local cet à 30 intervalles de temps rapprochés ; 3049045 43 - on détermine la valeur du taux de renouvellement d'air ACH du local en faisant converger : d'une part, un modèle diffusif exprimant la variation temporelle de la concentration du gaz à l'intérieur du local ca, en fonction de la concentration 5 du gaz à l'extérieur du local ceic et de paramètres physiques du local à partir desquels le taux de renouvellement d'air ACH du local est calculable, et d'autre part, l'évolution mesurée cac(t) de la concentration du gaz à l'intérieur du local ca, en fonction du temps. 10 13. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le modèle diffusif est un modèle R-C. 14. Programme d'ordinateur comportant des instructions pour l'exécution de tout ou partie des étapes d'un procédé selon l'une quelconque des revendications 1 à 13 lorsque ledit programme est exécuté par un 15 ordinateur. 15. Support d'enregistrement lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur comprenant des instructions pour l'exécution de tout ou partie des étapes d'un procédé selon l'une quelconque des revendications 1 à 13. 20 16. Dispositif de détermination du taux de renouvellement d'air ACH d'un local, ce dispositif étant caractérisé en ce qu'il comprend : - au moins un appareil (20) configuré pour appliquer, sur au moins deux périodes de temps Di, successives, des débits d'un gaz donné distincts dans le local ; 25 - au moins un capteur (40) configuré pour mesurer une concentration du gaz à l'intérieur du local ca, à intervalles de temps rapprochés ; - un terminal (10) comportant un module de traitement configuré pour faire converger, d'une part, un modèle diffusif exprimant la variation temporelle de la concentration du gaz à l'intérieur du local ca, en fonction de 30 la concentration du gaz à l'extérieur du local ceic et de paramètres physiques du local à partir desquels le taux de renouvellement d'air ACH du local est 3049045 44 calculable et, d'autre part, l'évolution mesurée cac(t) de la concentration du gaz à l'intérieur du local ca, en fonction du temps, de manière à obtenir la valeur du taux de renouvellement d'air du local. 17. Dispositif selon la revendication 16, caractérisé en ce qu'il 5 comprend en outre au moins un capteur (40) configuré pour mesurer une concentration du gaz à l'extérieur du local ceic à intervalles de temps rapprochés. 18. Dispositif selon l'une quelconque des revendication 16 ou 17, caractérisé en ce qu'il comprend en outre au moins un capteur de 10 température configuré pour mesurer la température à l'intérieur du local Tac. 19. Dispositif selon l'une quelconque des revendications 16 à 18, caractérisé en ce qu'il comprend des moyens de liaison, notamment sans fil, entre le ou chaque capteur (40) de mesure d'une concentration du gaz et le terminal (10). 15 20. Dispositif selon l'une quelconque des revendications 16 à 19, caractérisé en ce que le terminal (10) comporte des moyens de commande de l'appareil (20) configuré pour appliquer un débit du gaz. 21. Terminal (10) comportant un module de traitement configuré pour faire converger, d'une part, un modèle diffusif exprimant la variation 20 temporelle de la concentration d'un gaz donné à l'intérieur d'un local ca, en fonction de la concentration du gaz à l'extérieur du local ceic et de paramètres physiques du local à partir desquels le taux de renouvellement d'air ACH du local est calculable et, d'autre part, une évolution mesurée cac(t) de la concentration du gaz à l'intérieur du local ca, en fonction du temps, de 25 manière à obtenir la valeur du taux de renouvellement d'air du local. 22. Terminal (10) selon la revendication 21, caractérisé en ce que le module de traitement comporte un programme d'ordinateur selon la revendication 14, ledit programme étant enregistré sur un support d'enregistrement selon la revendication 15 constitué par une mémoire non 30 volatile réinscriptible (18) du terminal, les instructions dudit programme étant interprétables par un processeur (11) du terminal.