FR 2962216 A1 20120106 FR 1055419 A 20100705 PROCEDE DE DETERMINATION DE LA CONSOMMATION EN ENERGIE THERMIQUE D'UNE CELLULE La présente invention a pour objet un procédé de détermination de la consommation en énergie thermique d'une cellule appartenant à un ensemble de cellules chauffées chacune à l'aide d'au moins un aérotherme associé à un fluide caloporteur. L'invention trouve tout particulièrement son application dans la gestion thermique des entrepôts. Dans un entrepôt dont les cellules sont occupées par différents locataires, il est important de disposer d'une mesure de la consommation en énergie thermique des cellules reflétant aussi près que possible la consommation énergétique réelle. La consommation en énergie thermique de chaque cellule est généralement déterminée en fonction de la surface de la cellule. Or, selon les habitudes des occupants des cellules, cette répartition peut ne pas être représentative de la réalité. Il est également connu de déterminer cette consommation en énergie thermique à l'aide de calorimètres placés dans les cellules. Toutefois, les calorimètres sont chers, leur coût d'installation est élevé et ils nécessitent un entretien régulier pour rester fiables. La présente invention vise à remédier à ces inconvénients. L'invention propose en particulier un procédé permettant de déterminer de façon fiable et économique la consommation en énergie 25 thermique de chaque cellule. L'invention a ainsi pour objet un procédé de détermination de la consommation en énergie thermique d'une cellule (i) appartenant à un ensemble de cellules chauffées chacune à l'aide d'au moins un aérotherme associé à un fluide caloporteur issu d'un dispositif de chauffage de fluide. 30 Le procédé selon l'invention comprend les étapes consistant à : a) déterminer l'énergie thermique (Eboucle_i) transmise par le fluide caloporteur à ladite cellule (i), b) déterminer l'énergie thermique (Eboucle) transmise par le fluide caloporteur à l'ensemble des cellules, c) déterminer l'énergie thermique (Eaéros) prélevée au fluide caloporteur par l'ensemble des aérothermes des cellules, par la relation : Eaéros = Ecalo - Eboucle dans laquelle Ecalo désigne l'énergie thermique produite par le dispositif de chauffage de fluide et Eboucle désigne l'énergie thermique transmise par le fluide caloporteur à l'ensemble des cellules, d) déterminer l'énergie thermique (Eaéros_i) transmise à ladite cellule (i) par le ou les aérothermes présent(s) dans ladite cellule (i), en fonction de l'énergie thermique (Eaéros) prélevée au fluide caloporteur par l'ensemble des aérothermes des cellules, déterminée à l'étape c), et e) déterminer la consommation en énergie thermique de ladite cellule (i), comme étant égale à la somme de l'énergie thermique (Eboucle_i) transmise par le fluide caloporteur dans ladite cellule (i), déterminée à l'étape a), et de l'énergie thermique (Eaéros_i) transmise à ladite cellule (i) par le ou les aérothermes présent(s) dans ladite cellule (i), déterminée à l'étape d). Le fluide caloporteur pouvant être mis en circulation par une pompe, 20 l'énergie thermique (Ebouclei) transmise par le fluide caloporteur dans ladite cellule (i) peut être donnée par la relation : Eboucle_i =f Statut pompe* Kglobal * Li * ((Tdépart+Tretour)/ 2-Tambiante_i ) - dt dans laquelle: - Statut pompe vaut O si la pompe est à l'arrêt et vaut 1 si la pompe est en 25 marche, - Kglobal désigne un coefficient d'échange thermique calibrable, - le fluide caloporteur circulant dans un circuit de fluide, Li désigne le rapport de la longueur du circuit de fluide dans la cellule (i) sur la longueur totale du circuit de fluide, 30 - Tdépart désigne la température du fluide caloporteur partant du dispositif de chauffage de fluide, - Tretour désigne la température du fluide caloporteur revenant au dispositif de chauffage de fluide, - Tambiante_i désigne la température dans la cellule (i). L'énergie thermique (Eboucle) transmise par le fluide caloporteur dans l'ensemble des cellules est la somme des énergies thermiques (Eboucle_i) transmises par le fluide caloporteur à chaque cellule (i). L'énergie 5 thermique (Eboucle) peut ainsi être donnée par la relation : Eboucle -JStatut pompe* Kglobal * Li*((Tdépart+Tretour)/2-Tambiante_i )-dt Le coefficient d'échange thermique Kglobal peut être donné par la relation : Kglobal - Eboucle 10 J(Statut pompe*li Li*(Tmoy-boucle - Tambiante_i))-dt dans laquelle: - Eboucle désigne l'énergie thermique transmise par le fluide caloporteur dans l'ensemble des cellules, - Statut pompe vaut O si la pompe est à l'arrêt et vaut 1 si la pompe est en 15 marche, - le fluide caloporteur circulant dans un circuit de fluide, Li désigne le rapport de la longueur du circuit de fluide dans la cellule (i) sur la longueur totale du circuit de fluide, - Tmoy-boucle est égal à (Tdépart+Tretour)/2, Tdépart désignant la température 20 du fluide caloporteur partant du dispositif de chauffage de fluide, et Tretour désignant la température du fluide caloporteur revenant au dispositif de chauffage de fluide, - Tambiante_i désigne la température dans la cellule (i). L'énergie thermique (Eaéros_i) transmise à ladite cellule (i) par le ou les 25 aérothermes présent(s) dans ladite cellule (i) peut être donnée par la relation : Eaéros_i - Statut pompe *Caéros_i *(Tmoy-boucle - Tambiante_1) * Eaéros Caéros_i *(Tmoy-boucle - Tambiante_1) 30 dans laquelle : - Statut pompe vaut O si la pompe est à l'arrêt et vaut 1 si la pompe est en marche, - Caéros désigne la consommation électrique du ou des aérothermes de la cellule (i), - Tmoy-boucle est égal à (Tdépart+Tretour)/2, Tdépart désignant la température du fluide caloporteur partant du dispositif de chauffage de fluide, et Tretour désignant la température du fluide caloporteur revenant au dispositif de chauffage de fluide, - Tambiante_i désigne la température dans la cellule i, Eaéros désigne l'énergie thermique prélevée au fluide caloporteur par l'ensemble des aérothermes des cellules. L'énergie thermique (Eaéros_i) transmise à ladite cellule (i) par le ou les aérothermes présent(s) dans ladite cellule (i) peut être déterminée sans utiliser de compteur de l'énergie électrique consommée par le ou les aérothermes. L'énergie thermique (Eaéros_i) transmise à ladite cellule (i) par le ou les aérothermes présent(s) dans ladite cellule peut ainsi être donnée par la relation : Eaéros_i = Statut pompe * Statut aéros_i *(Tmoy-boucle - Tambiante_i) * Eaéros Statut aéros_i *(Tmoy-boucle - Tambiante_i) dans laquelle : - Statut pompe vaut O si la pompe est à l'arrêt et vaut 1 si la pompe est en 20 marche, - 5tatut aéros_i vaut 0 si le ou les aérothermes présent(s) dans la cellule (i) sont à l'arrêt et vaut 1 si le ou les aérothermes présent(s) dans la cellule (i) sont en marche, - Tmoy-boucle est égal à (Tdépart+Tretour)/2, Tdépart désignant la température 25 du fluide caloporteur partant du dispositif de chauffage de fluide, et Tretour désignant la température du fluide caloporteur revenant au dispositif de chauffage de fluide, - Tambiante_i désigne la température dans la cellule i, - Eaéros désigne l'énergie thermique prélevée au fluide caloporteur par 30 l'ensemble des aérothermes des cellules. L'énergie thermique (Ecalo) produite par le dispositif de chauffage de fluide peut être donnée par un calorimètre mesurant la chaleur produite par le dispositif de chauffage. L'énergie thermique (Ecalo) produite par le dispositif de chauffage de fluide peut également être déterminée sans utiliser de calorimètre mesurant la chaleur produite par le dispositif de chauffage. Elle peut ainsi être donnée par la relation Ecalo = f 3600*d*c*p*(Tdépart-Tretour)-dt dans laquelle : - d désigne le débit du fluide caloporteur en m3/h, - c désigne la chaleur massique du fluide caloporteur en kJ/kg.K, - p désigne la densité du fluide caloporteur en kg/m3, - Tdépart désigne la température du fluide caloporteur partant du dispositif de chauffage de fluide, - Tretour désigne la température du fluide caloporteur revenant au dispositif de chauffage de fluide. Le procédé peut comprendre en outre une étape f) de détermination 15 du ratio (rconso_i) de consommation en énergie thermique de ladite cellule (i) à l'aide de la relation : rconso_i = (Eaéro_i + Eboucle_i) (Eaéroi + Eboucle_i ) D'autres caractéristiques et avantages de la présente invention apparaîtront plus clairement à la lecture de la description suivante donnée à 20 titre d'exemple illustratif et non limitatif et faite en référence à la figure 1 annexée qui illustre un dispositif soumis au procédé selon l'invention. Tel qu'illustré à la figure 1, un entrepôt 1 comprend par exemple six cellules 2 occupées par différents locataires. Une chaudière 3 produit de la chaleur véhiculée par un fluide caloporteur, typiquement de l'eau, qui circule 25 dans une boucle d'eau chaude 4 (ou circuit d'eau chaude 4) au moyen d'une pompe d'un circulateur, non représentée, par exemple une pompe de la chaudière 3. La boucle 4 comprend une boucle aller partant de la chaudière 3 et une boucle retour retournant vers la chaudière 3. Dans chaque cellule 2, un ou 30 plusieurs aérothermes 5 prélèvent de la chaleur à la fois à la boucle aller et à la boucle retour. Un aérotherme est un générateur soufflant de l'air chaud à l'aide par exemple d'un ventilateur afin d'élever la température d'un local. Les aérothermes 5 prélèvent de la chaleur sur la boucle d'eau chaude 4 et la transmettent aux cellules 2. Le dispositif comprend en outre un capteur de température 6 disposé dans chaque cellule 2, ainsi qu'un capteur de température 7 placé au départ de la boucle d'eau chaude 4 et un capteur de température 8 placé au retour de la boucle d'eau chaude 4. Le dispositif peut également comprendre un compteur d'énergie électrique 9 relié aux aérothermes 5 dans chaque cellule 2, ainsi qu'un calorimètre 10 mesurant la chaleur produite par la chaudière 3. Le dispositif enregistre en permanence les paramètres mesurés. L'hypothèse de départ est que l'énergie thermique produite par la chaudière 3 de l'entrepôt 1, notée Ecalo, est communiquée à l'entrepôt 1 sous forme d'énergie dissipée par déperditions thermiques au niveau de la boucle d'eau chaude 4 parcourant les cellules 2 (non calorifugées en général), notée Eboucle, et sous forme d'énergie prélevée par les aérothermes 5 des cellules 2 sur la boucle d'eau chaude 4 de la chaudière 3, notée Eaéros. Toute l'énergie produite par la chaudière 3 est transmise via la boucle d'eau chaude 4 aux cellules 2 de l'entrepôt 1 sous l'une ou l'autre des deux formes décrites ci-dessus. On en déduit la relation : Ecalo = Eaéros + Eboucle. L'énergie thermique Ecalo est mesurée au niveau de la chaudière 3 par le calorimètre 10. Le procédé selon l'invention va permettre de calculer pour 25 chaque cellule (i), Eboucle_i et d'y ajouter Eaéros _i- Le calcul de l'énergie dissipée via la boucle d'eau chaude 4 repose sur trois hypothèses. La première hypothèse et une hypothèse d'uniformité des caractéristiques physiques de la boucle d'eau chaude 4. Afin de calculer 30 l'énergie Eboucle dissipée par les déperditions de la boucle d'eau chaude 4, on fait l'hypothèse que la boucle 4 a les mêmes caractéristiques physiques sur toute sa longueur, du départ au retour vers la chaudière 3. Les déperditions dans chaque cellule 2 sont caractérisées de façon uniforme par un coefficient Kglobal mais dépendent uniquement de la différence entre la température de l'eau de la boucle 4 et la température ambiante de la cellule 2. Tout au long de la boucle 4, la température de l'eau chaude diminue linéairement à mesure que la distance entre la portion de boucle 4 étudiée et la chaudière 3 augmente. La boucle 4 part de la chaudière 3, traverse chacune des cellules 2 puis effectue le chemin inverse. On peut donc considérer dans une deuxième hypothèse que la température de la boucle 4, dans toute sa longueur, est égale à la moyenne des températures de départ et de retour de la boucle 4 : Tmoy-boucle = (Tdepboucle + Tretboucle )/2 La troisième hypothèse est une hypothèse de répartition linéique de la boucle d'eau chaude 4 dans les cellules 2. La longueur de boucle 4 traversant chaque cellule 2 est également un paramètre pouvant influer sur le calcul de Eboucle. La répartition linéique de boucle 4 dans les cellules 2 peut ne pas être homogène. Dans ce cas, on définit pour chaque cellule (i) un coefficient Li de répartition linéique de la boucle 4. Ce coefficient implique une variation de la dissipation de chaleur en fonction de la quantité de boucle 4 présente dans chacune des cellules. Il est défini comme suit : Li = Lcanalisation i /Lcanalisation total Il vaut 1 si la canalisation est répartie dans l'entrepôt 1 de façon homogène. Il est différent de 1 dans le cas contraire. Finalement, l'énergie dissipée via la boucle d'eau chaude 4 par la cellule (i) est : Eboucle_i = J (Statut pompe* Li *Kglobal * (Tmoy-boucle - Tamb_i)) - dt Au niveau de l'entrepôt 1, l'énergie dissipée via la boucle d'eau chaude 4 se calcule donc comme la somme des déperditions via la boucle 4, 30 de toutes les cellules 2 : Eboucle = JStatut pompe* Kglobal*li Li*(Tmoy-boucle - Tamb_i))-dt Le coefficient Kglobal est calculé durant une étape de calibrage, qui sera décrite plus loin, et le coefficient Li est lui évalué à partir des plans et des données techniques des équipements. Ayant calculé l'énergie Eboucle, on déduit du calcul l'énergie prélevée par 5 les aérothermes 5 des cellules 2 sur la boucle d'eau chaude 4. Cette énergie s'exprime donc comme suit : Eaéros = Ecalo - Eboucle. Eaéros est transmise globalement aux cellules 2 de l'entrepôt 1. Elle se répartit en outre, selon les cellules 2, en fonction de l'activité des 10 aérothermes 5 et de la différence entre la température de la boucle d'eau chaude 4 et la température ambiante de la cellule 2. Or, l'activité des aérothermes 5 d'une cellule 2 est décrite par leur consommation électrique Caéros_i- L'énergie diffusée dans chacune des cellules i via les aérothermes 5 est 15 donc calculée comme suit : Eaéros_i = Statut pompe *Caéros_i *(Tmoy-boucle - Tambiante_i) * Eaéros Caéros_i *(Tmoy-boucle - Tambiante 20 Pour calculer le coefficient Kglobal, il faut auparavant isoler les dates auxquelles : - aucun aérotherme 5 ne fonctionne (Caéros_i =0 donc Eaéros =0), - la pompe de la chaudière 3 est en marche (Statut pompe=l). A ces dates précises, on obtient les relations suivantes : 25 Ecalo = Eboucle D'où : Kglobal = Efluide f (Statut pompe*£i Li*(Tmoy-boucle - Tambiante_i))-dt 30 Obtenant un Kglobal pour chacune des dates de mesure, on réalise la moyenne sur l'ensemble des calculs effectués. Pour éviter les distorsions liées aux effets de bord, le calibrage peut être réalisé sur une période de temps supérieure à 24 heures et les valeurs relevées pendant le régime transitoire sont alors écartées du calcul. Les données relevées nécessaires au calcul de répartition des consommations locatives sont : - les températures ambiantes de chacune des cellules (Tambi), - la température de départ de la boucle d'eau chaude 4 et la température de retour de la boucle 4, desquelles on déduit la température moyenne de la boucle (Tdepboucle, Tretboucle), - la consommation électrique de l'ensemble des aérothermes 5 de chaque cellule (Caéros_i). - l'énergie thermique produite par la chaudière (Ecalo en kWh), de laquelle on déduit également l'état de marche de la pompe (Statut pompe). Les données sont traitées mois par mois de façon à pouvoir calculer une répartition des consommations mensuelles. A partir du Kglobal déterminé par le calibrage, on calcule l'énergie dissipée via la boucle d'eau chaude 4 dans la cellule (i). Pour chaque cellule et pour toute date de mesure, on calcule : Eboucle_i =f Statut pompe* Kglobal * Li * ((Tdépart+Tretour)/ 2-Tambiante_i ) - dt On réalise la somme sur le mois, et sur les cellules 2, des valeurs trouvées. On obtient Eboucle. On accède alors à la quantité d'énergie dissipée via la boucle d'eau 20 chaude 4 au sein de l'entrepôt 1 et sur le mois étudié. Pour déterminer l'énergie diffusée via les aérothermes 5, la première étape est de calculer Eaéros à partir de Ecalo et de Eboucle. On effectue donc d'abord : Eaéros = Ecalo - Eboucle 25 Ensuite, pour chaque cellule, on calcule l'énergie thermique diffusée via les aérothermes 5 pour les dates où : - la pompe de la chaudière 3 est en fonctionnement (Statut pompe =1) - les aérothermes 5 de la cellule fonctionnent (Caéros_i >0). Pour chaque date isolée, on applique la formule suivante : 30 Eaéros_i = Statut pompe * Statut aéros_i * (Tmoy-boucle - Tambiante_1) * Eaéros Statut aéros_i *(Tmoy-boucle - Tambiante_1) On réalise la somme sur le mois de ces valeurs pour obtenir l'énergie thermique diffusée dans la cellule via les aérothermes 5 sur le mois. La consommation d'énergie thermique totale de chaque cellule (i) s'exprime comme la somme de l'énergie dissipée via la boucle d'eau chaude Eboucle_i et de l'énergie diffusée via les aérothermes Eaéro_i- On évalue finalement le ratio de consommation global pour chaque cellule (i) : rconso = (Eaéro_i + Eboucle/ (Eaéroi + Eboucle_i ) Pour chaque mois étudié, on calcule le rendement global des aérothermes 5 de chaque cellule 2. Il est défini par Eaéro i/Caérosi. Ce paramètre permet d'évaluer I« Les données nécessaires au calcul de répartition des consommations de chaque cellule sont cette fois-ci : - les températures ambiantes de chacune des cellules (Tambi), - la température de départ de la boucle d'eau chaude 4 et la température de retour de la boucle 4, desquelles on déduit la température moyenne de la boucle (Tdepboucle, Tretboucle), - l'activité de l'ensemble des aérothermes 5 de chaque cellule (Statut aéros_i) Les données sont traitées mois par mois de façon à pouvoir calculer 10 une répartition des consommations mensuelle. A partir des températures de départ et de retour de la boucle d'eau chaude 4 ainsi que de la donnée initiale du débit du circulateur, on calcule, pour tous les incréments de mesures, l'énergie produite par la chaudière 3 avec la formule décrite ci-dessus : 15 Ecalo = J(3600*d*c*p*(Tdepboucle - Tretboucle))-dt L'énergie dissipée via la boucle d'eau chaude 4 est calculée comme décrit précédemment. La première étape reste de calculer Eaéros à partir de Ecalo et de Eboucle. On effectue donc : 20 Eaéros = Ecalo - Eboucle Ensuite, pour chaque cellule, on calcule l'énergie thermique diffusée via les aérothermes 5 pour les dates où : - la pompe de la chaudière 3 est en fonctionnement (Statut pompe =1) - les aérothermes 5 de la cellule fonctionnent (Statut aéros_i =1) 25 à l'aide de la formule suivante : Eaéros_i = Statut pompe * Statut aéros_i * (Tmoy-boucle - Tambiante_1) * Eaéros Statut aéros_i *(Tmoy-boucle - Tambiante_1) On réalise la somme sur le mois des résultats exprimés en kWh pour 30 obtenir l'énergie thermique diffusée dans la cellule via les aérothermes 5. La consommation d'énergie thermique totale de chaque cellule (i) s'exprime comme la somme de l'énergie thermique diffusée via les 13 aérothermes Eaéro_i et de l'énergie thermique dissipée via la boucle d'eau chaude Eboucle l- On détermine ensuite le ratio de consommation global pour chaque cellule (i) : rconso = (Eaéro_i + Eboucle/ (Eaéroi + Eboucle_i ) Le procédé selon l'invention permet ainsi de déterminer de manière précise la répartition des consommations de chaque cellule. Ce procédé repose, non plus sur la surface de chacune des cellules, mais sur le comportement thermique des équipements au sein des cellules. Le procédé permet d'obtenir une répartition fiable des consommations de chauffage, avec ou sans comptage des consommations électriques des aérothermes et avec ou sans calorimètre dans la chaudière. REVENDICATIONS1. Procédé de détermination de la consommation en énergie thermique d'une cellule (i) appartenant à un ensemble de cellules (2) chauffées chacune à l'aide d'au moins un aérotherme (5) associé à un fluide caloporteur issu d'un dispositif (3) de chauffage de fluide, caractérisé en ce que le procédé comprend les étapes consistant à : a) déterminer l'énergie thermique (Eboucle_i) transmise par le fluide caloporteur à ladite cellule (i), b) déterminer l'énergie thermique (Eboucle) transmise par le fluide caloporteur à l'ensemble des cellules (2), c) déterminer l'énergie thermique (Eaéros) prélevée au fluide caloporteur par l'ensemble des aérothermes (5) des cellules (2), par la relation : Eaéros = Ecalo - Eboucle dans laquelle Ecalo désigne l'énergie thermique produite par le dispositif (3) de chauffage de fluide et Eboucle désigne l'énergie thermique transmise par le fluide caloporteur à l'ensemble des cellules (2), d) déterminer l'énergie thermique (Eaéros_i) transmise à ladite cellule (i) par le ou les aérothermes (5) présent(s) dans ladite cellule (i), en fonction de l'énergie thermique (Eaéros) prélevée au fluide caloporteur par l'ensemble des aérothermes (5) des cellules (2), déterminée à l'étape c), et e) déterminer la consommation en énergie thermique de ladite cellule (i), comme étant égale à la somme de l'énergie thermique (Eboucle_l) transmise par le fluide caloporteur dans ladite cellule (i), déterminée à l'étape a), et de l'énergie thermique (Eaéros_i) transmise à ladite cellule (i) par le ou les aérothermes (5) présent(s) dans ladite cellule (i), déterminée à l'étape d). 2. Procédé selon la revendication 1, caractérisé en ce que, le fluide caloporteur étant mis en circulation par une pompe, l'énergiethermique (Ebouclel) transmise par le fluide caloporteur dans ladite cellule (i) est donnée par la relation : Eboucle_i = JStatut pompe* Kglobal * Li * ( (Tdépart+Tretour)/2-Tambiante_1 ) - dt dans laquelle: - Statut pompe vaut 0 si la pompe est à l'arrêt et vaut 1 si la pompe est en marche, - Kglobal désigne un coefficient d'échange thermique calibrable, - le fluide caloporteur circulant dans un circuit de fluide (4), Li désigne le rapport de la longueur du circuit de fluide (4) dans la cellule (i) sur la longueur totale du circuit de fluide (4), - Tdépart désigne la température du fluide caloporteur partant du dispositif (3) de chauffage de fluide, - Tretour désigne la température du fluide caloporteur revenant au dispositif (3) de chauffage de fluide, - Tambiante_i désigne la température dans la cellule (i). 3. Procédé selon la revendication 2, caractérisé en ce que l'énergie thermique (Eboucle) transmise par le fluide caloporteur dans l'ensemble des cellules (2) est donnée par la relation : Eboucle - JStatut pompe* Kglobal *El Li * ( (Tdépart+Tretour)/2-Tambiante_1 ) - dt 4. Procédé selon la revendication 2 ou 3, caractérisé en ce que le coefficient d'échange thermique Kglobal est donné par la relation : Kglobal = Eboucle J(Statut pompe*£i Li*(Tmoy-boucle - Tambiante_i))-dt dans laquelle: - Eboucle désigne l'énergie thermique transmise par le fluide caloporteur dans l'ensemble des cellules (2), - Statut pompe vaut 0 si la pompe est à l'arrêt et vaut 1 si la pompe est en marche, - le fluide caloporteur circulant dans un circuit de fluide (4), Li désigne le rapport de la longueur du circuit de fluide (4) dans la cellule (i) sur la longueur totale du circuit de fluide (4), - Tmoy-boucle est égal à (Tdépart+Tretour)/2, Tdépart désignant la température du fluide caloporteur partant du dispositif (3) de chauffage de fluide, et Tretour désignant la température du fluide caloporteur revenant au dispositif (3) de chauffage de fluide, - Tambiante_i désigne la température dans la cellule (i). 5. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que l'énergie thermique (Eaéros_i) transmise à ladite cellule (i) par le ou les aérothermes (5) présent(s) dans ladite cellule (i) est donnée par la relation : Eaéros = Statut pompe *Caéros_i *(Tmoy-boucle - Tambiante_i) * Eaéros Caéros_i *(Tmoy-boucle - Tambiante dans laquelle : - Statut pompe vaut 0 si la pompe est à l'arrêt et vaut 1 si la pompe est en marche, - Caéros_i désigne la consommation électrique du ou des aérothermes (5) de la cellule (i), - Tmoy-boucle est égal à (Tdépart+Tretour)/2, Tdépart désignant la température du fluide caloporteur partant du dispositif (3) de chauffage de fluide, et Tretour désignant la température du fluide caloporteur revenant au dispositif (3) de chauffage de fluide, - Tambiante_i désigne la température dans la cellule (i),- Eaéros désigne l'énergie thermique prélevée au fluide caloporteur par l'ensemble des aérothermes (5) des cellules (2). 6. Procédé selon l'une des revendications 1 à 5, caractérisé en ce que l'énergie thermique (Eaéros_i) transmise à ladite cellule (i) par le ou les aérothermes (5) présent(s) dans ladite cellule (i) est déterminée sans utiliser de compteur (9) de l'énergie électrique consommée par le ou les aérothermes (5). 7. Procédé selon l'une des revendications 1 à 4 ou 6, caractérisé en ce que l'énergie thermique (Eaéros_i) transmise à ladite cellule (i) par le ou les aérothermes (5) présent(s) dans ladite cellule (i) est donnée par la relation : Eaéros_i = Statut pompe * Statut aéros_i *(Tmoy-boucle - Tambiante_i) * Eaéros Statut aéros_i *(Tmoy-boucle - Tambiante_i) dans laquelle : - Statut pompe vaut 0 si la pompe est à l'arrêt et vaut 1 si la pompe est en marche, - Statutaéros_i vaut 0 si le ou les aérothermes (5) présent(s) dans la cellule (i) sont à l'arrêt et vaut 1 si le ou les aérothermes (5) présent(s) dans la cellule (i) sont en marche, - Tmoy-boucle est égal à (Tdépart+Tretour)/2, Tdépart désignant la température du fluide caloporteur partant du dispositif (3) de chauffage de fluide, et Tretour désignant la température du fluide caloporteur revenant au dispositif (3) de chauffage de fluide, - Tambiante_i désigne la température dans la cellule (i), - Eaéros désigne l'énergie thermique prélevée au fluide caloporteur par l'ensemble des aérothermes (5) des cellules (2). 8. Procédé selon l'une des revendications 1 à 7, caractérisé en ce que l'énergie thermique (Ecalo) produite par le dispositif (3) de chauffage de fluide est déterminée sans utiliser de calorimètre (10) mesurant la chaleur produite par le dispositif (3) de chauffage. 9. Procédé selon l'une des revendications 1 à 8, caractérisé en ce que l'énergie thermique (Ecalo) produite par le dispositif (3) de chauffage de fluide est donnée par la relation : Ecalo = J3600*d*c*p*(Tdépart-Tretour)-dt dans laquelle : - d désigne le débit du fluide caloporteur en m3/h, - c désigne la chaleur massique du fluide caloporteur en kJ/kg.K, - p désigne la densité du fluide caloporteur en kg/m3, - Tdépart désigne la température du fluide caloporteur partant du dispositif (3) de chauffage de fluide, - Tretour désigne la température du fluide caloporteur revenant au dispositif (3) de chauffage de fluide. 10.Procédé selon l'une des revendications 1 à 9, caractérisé en ce qu'il comprend en outre une étape f) de détermination du ratio (rconsoi) de consommation en énergie thermique de ladite cellule (i) à l'aide de la relation : rconso_i = (Eaéro_i + Eboucle_i) (Eaéroi + Eboucle_i )