Abstract:
A method determining heat loss coefficient K of a premises, includes: in unoccupied premises, performing a campaign of measurements of at least one temperature inside the premises T ik  at closely-spaced time intervals over at least two successive time periods D k  corresponding to distinct heating powers P totk  of the premises; determining temperature of outside air T ek  at the same closely-spaced times; for each time period D k , on the basis of evolution T ik (t) of a quantity T ik  as a function of time, selecting a time interval Δt k  for which the evolution T ik (t) is substantially linear, then determining the slope α k  of the tangent to the evolution T ik (t) over this time interval Δt k , and deducing the value of the heat loss coefficient K of the premises on the basis of the slopes α k .

Description:
BACKGROUND 
     The present invention relates to a method and a device for determining the heat loss coefficient of a premises. Within the meaning of the invention, a premises is an individual house or a building, in particular for residential or tertiary use, or else a part of such a building, for example an apartment in a building with multiple floors. 
     The heat loss coefficient of a premises, denoted K, is equal to the heat loss power of the premises (in Watts) per degree (Kelvin or Celsius) of difference between the temperature of the air inside the premises and the temperature of the outside air. This coefficient K is representative of the energy performance of the envelope of the premises. 
     The heat loss coefficient K of a premises is influenced, on the one hand, by the heat losses by transmission through the walls of the premises and, on the other hand, by infiltrations of air. The heat losses by transmission are represented by a factor H T =UA T , where U is the heat transfer coefficient of the envelope of the premises, also called the specific heat transfer coefficient of the premises, and A T  is the total area of the walls of the premises. The infiltrations of air into the premises are represented by a factor m′.Cp, where m′ is the air renewal flow rate and Cp is the heat capacity of the air. Consequently, the heat loss coefficient K is given by the relation:
 
 K=H   T   +m′.Cp=UA   T   +m′.Cp.  
 
     The coefficient U is used, within the framework of thermal regulations such as RT 2005 in France or the regulation EnEV in Germany, to estimate the overall energy consumption of the premises. Its determination is useful in order to perform a diagnosis of the thermal insulation of a premises, in particular after its construction, to verify that the constructor has adhered to the standards in force in terms of thermal insulation both in terms of choice of materials and their implementation, or when a refurbishment of the premises is considered, with a view to evaluating the measures that should be taken to improve the thermal performance. 
     It is known to determine the heat loss coefficient of a premises by means of computational softwares which involve modeling the envelope of the premises. Such softwares are relatively tricky to implement and have the drawback of providing only a theoretical result, which does not assess such real parameters as the effective implementation of the insulation materials, the construction technology, etc. 
     Moreover, it is known to determine the heat loss coefficient of a premises by carrying out in situ measurements in the premises over long periods, which generally extend over several weeks or months, and then by undertaking a statistical analysis of these measurements. The statistical analysis is made necessary by the multiplicity of parameters which influence the thermal behavior of the premises over the measurement period, in particular the meteorological conditions and the occupancy conditions of the premises. These known in situ measurement methods have the drawback of being lengthy and of involving significant and expensive equipment. 
     BRIEF SUMMARY 
     It is these drawbacks that the invention is intended more particularly to remedy by proposing a method and a device making it possible to determine in a fast manner the heat loss coefficient of a premises, with moderate cost and reasonable accuracy. 
     For this purpose, one subject of the invention is a method for determining the heat loss coefficient K of a premises, characterized in that it comprises steps in which:
         in the unoccupied premises, a campaign of measurements of at least one temperature inside the premises T ik  at closely-spaced time intervals is performed over at least two successive time periods D k  corresponding to distinct heating powers P totk  of the premises;   the temperature of the outside air T ek  at said same closely-spaced times is determined;   for each time period D k , on the basis of the evolution T ik (t) of the quantity T ik  as a function of time, this evolution is quantitatively analyzed through a simple mathematical model:
           either, if there exists a time interval Δt k  for which the evolution T ik (t) is substantially linear, the slope α k  of the tangent to the curve T ik (t) is determined over this time interval Δt k  and the value of the heat loss coefficient K of the premises is deduced on the basis of the slopes α k ;   or, if there does not exist any time interval for which the evolution T ik (t) is substantially linear, a time interval Δt k ′ is selected over which the evolution T ik (t) is substantially exponential of type exp(−t/τ), with τ the thermal time constant of the premises, and the value of the heat loss coefficient K of the premises is deduced, which is the value such that the curve   
               

             Ln   ⁡     [       (         θ   k     ⁡     (   t   )       -       P     tot   ⁢           ⁢   k       K       )     /     (         θ   k     ⁡     (   0   )       -       P     tot   ⁢           ⁢   k       K       )       ]           
is a straight line, where θ k (t)=T ik (t)−T ekm ′ with T ekm ′ the average of the temperature measurements of the outside air T ek  over the time interval Δt k ′.
 
     The principle underlying the invention is to use the transient variations of the inside temperature of the premises subjected to controlled internal impulses and in a measured external environment. Quantitative analysis of the variation of the inside temperature of the premises makes it possible to quantitatively determine the energy efficiency of the premises over a short period, extending over a few hours, while limiting the number of parameters liable to influence the thermal behavior of the premises. In particular, the brevity of the measurements makes it possible to circumvent the influence of the conditions of use of the premises and of the variations of the outside climatic conditions. 
     Of course, the method according to the invention does not necessarily require the setting up of a graphical representation of the evolution T ik (t). 
     In particular, in the case where there exists a time interval Δt k  for which the evolution T ik (t) is substantially linear, the slope α k  of the tangent to the curve T ik (t) over the time interval Δt k  is equal to the derivative of the evolution T ik (t) over this interval Δt k . Hence, the step of determining the slope α k  of the tangent to the curve T ik (t) over the time interval Δt k  may be carried out, within the framework of the invention, by computing the derivative of the evolution T ik (t) over the time interval Δt k , without resorting to a graphical representation of the evolution T ik (t). 
     In the same manner, in the case where there does not exist any time interval for which the evolution T ik (t) is substantially linear, the determination of the value of the heat loss coefficient K of the premises such that the curve 
             Ln   ⁡     [       (         θ   k     ⁡     (   t   )       -       P     tot   ⁢           ⁢   k       K       )     /     (         θ   k     ⁡     (   0   )       -       P     tot   ⁢           ⁢   k       K       )       ]           
is a straight line does not necessarily involve a graphical representation of this curve, and can be done via a mathematical linearization of the evolution
 
     
       
         
           
             
               Ln 
               ⁡ 
               
                 [ 
                 
                   
                     ( 
                     
                       
                         
                           θ 
                           k 
                         
                         ⁡ 
                         
                           ( 
                           t 
                           ) 
                         
                       
                       - 
                       
                         
                           P 
                           
                             tot 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             k 
                           
                         
                         K 
                       
                     
                     ) 
                   
                   / 
                   
                     ( 
                     
                       
                         
                           θ 
                           k 
                         
                         ⁡ 
                         
                           ( 
                           0 
                           ) 
                         
                       
                       - 
                       
                         
                           P 
                           
                             tot 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             k 
                           
                         
                         K 
                       
                     
                     ) 
                   
                 
                 ] 
               
             
             . 
           
         
       
     
     The computation steps of the method, in particular for the determination of the slopes α k  or for the linearization of the evolution 
               Ln   ⁡     [       (         θ   k     ⁡     (   t   )       -       P     tot   ⁢           ⁢   k       K       )     /     (         θ   k     ⁡     (   0   )       -       P     tot   ⁢           ⁢   k       K       )       ]       ,         
may be implemented with the aid of any appropriate computation means. This may entail in particular an electronic computation unit comprising at the same time means for acquiring the measurements of temperature inside the premises T ik  and means for computing the heat loss coefficient K of the premises on the basis of these measurements.
 
     According to the invention, the expression “heating power of the premises” is intended to mean any operative condition generating a variation of the inside temperature of the premises for given conditions of outside temperature. It is understood that the heating power P totk  may be positive, zero or negative. A positive heating power corresponds to a supply of heat into the premises, whilst a negative heating power corresponds to a supply of cold into the premises. In the case of a zero heating power, the variation of the inside temperature of the premises may result from a difference between the inside temperature of the premises and the outside temperature, or else from a variation of the outside temperature. According to the invention, provision is made for at least one of the heating powers P totk  to be non-zero. 
     In an advantageous manner, for each time period D k , the heating power P totk  of the premises comprises an imposed heating power P impk  imposed by means of a controlled power source. If no power source other than that used to apply the imposed heating power P impk  is active in the premises in the course of the time period D k , the total heating power P totk  of the premises is equal to the imposed heating power P impk . If on the other hand there exists, in the course of the period D k , an additional power P supk  in the premises supplementing the power P impk , the total heating power of the premises is equal to P impk +P supk . In particular, in the case where the solar radiation over the time period D k  is significant, the contribution of the solar radiation to the heating of the premises forms part of the additional power P supk . 
     In practice, the conditions of implementation of the method are adapted so as to limit the supplies of additional power P supk  other than the imposed heating power P impk , in particular by taking care that the premises is unoccupied. 
     In an advantageous manner, the method is implemented over time periods D k  for which the solar radiation is weak, preferably zero. In a preferred manner, the method is implemented over time periods D k  chosen during the night, or optionally in the daytime in the morning or the evening. It is thus possible to reduce the contribution of the solar radiation and to limit the fluctuations of the temperature of the outside air. 
     The time periods D k  may be either disjoint, or follow one another in immediate succession. In the latter case, it may be considered that the method is carried out in its entirety over a continuous time period, formed by the succession of time periods D k . Preferably, with a view to limiting the implementation time of the method while reducing the contribution of the solar radiation, the method is carried out in its entirety continuously over a single nocturnal period. 
     Preferably, over each time period D k , any fixed ventilation system fitted to the premises is deactivated and all the ventilation inlets are closed or plugged, so as to limit the exchanges of air with the outside. 
     As a variant, the fixed ventilation systems of the premises can operate in the course of the method over each time period D k . However, this introduces an additional term relating to air renewal into the expression for the heat loss coefficient K:
 
 K=H   T   +m′   1   .Cp+m′   2   .Cp,  
 
where m′ 1  the renewal flow rate of air by infiltration and m′ 2  the renewal flow rate of air due to the fixed ventilation systems are correlated, the value of one depending on the value of the other.
 
     The method according to the invention rests upon a modeling of the premises in the form of an isothermal box characterized, on the one hand, by its heat loss coefficient K, and, on the other hand, by its inertia or its effective heat capacity C. The effective heat capacity C of a premises, which corresponds to the heat capacity of the materials situated in the insulating envelope of the premises, is defined as the energy necessary to increase the ambient temperature of the premises by 1 K at constant outside temperature. 
     In a general manner, for each time period D k , the energy balance of the premises may be written: 
                       K   ⁡     (       T     i   ⁢           ⁢   k       -     T   ek       )       +     C   ⁢       ⅆ     (       T   ik     -     T   ek       )         ⅆ   t           =     P   totk             (   1   )               
with
         K the heat loss coefficient of the premises,   T ik (t) the inside temperature of the premises,   T ek (t) the temperature of the outside air,   C the effective heat capacity of the premises,   P totk  the total heating power of the premises.       

     Equation (1) admits as solution: 
     
       
         
           
             
               
                 
                   T 
                   ik 
                 
                 ⁡ 
                 
                   ( 
                   t 
                   ) 
                 
               
               - 
               
                 
                   T 
                   ek 
                 
                 ⁡ 
                 
                   ( 
                   t 
                   ) 
                 
               
             
             = 
             
               
                 
                   P 
                   
                     tot 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     k 
                   
                 
                 K 
               
               + 
               
                 
                   [ 
                   
                     
                       ( 
                       
                         
                           
                             T 
                             ik 
                           
                           ⁡ 
                           
                             ( 
                             0 
                             ) 
                           
                         
                         - 
                         
                           
                             T 
                             ek 
                           
                           ⁡ 
                           
                             ( 
                             0 
                             ) 
                           
                         
                       
                       ) 
                     
                     - 
                     
                       
                         P 
                         
                           tot 
                           ⁢ 
                           
                               
                           
                           ⁢ 
                           k 
                         
                       
                       K 
                     
                   
                   ] 
                 
                 ⁢ 
                 
                   
                     ⅇ 
                     
                       
                         
                           - 
                           K 
                         
                         · 
                         t 
                       
                       C 
                     
                   
                   . 
                 
               
             
           
         
       
     
     This solution may be linearized to short times: 
     
       
         
           
             
               
                 
                   T 
                   ik 
                 
                 ⁡ 
                 
                   ( 
                   t 
                   ) 
                 
               
               - 
               
                 
                   T 
                   ek 
                 
                 ⁡ 
                 
                   ( 
                   t 
                   ) 
                 
               
             
             = 
             
               
                 
                   P 
                   
                     tot 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     k 
                   
                 
                 K 
               
               + 
               
                 
                   
                     [ 
                     
                       
                         ( 
                         
                           
                             
                               T 
                               ik 
                             
                             ⁡ 
                             
                               ( 
                               0 
                               ) 
                             
                           
                           - 
                           
                             
                               T 
                               ek 
                             
                             ⁡ 
                             
                               ( 
                               0 
                               ) 
                             
                           
                         
                         ) 
                       
                       - 
                       
                         
                           P 
                           
                             tot 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             k 
                           
                         
                         K 
                       
                     
                     ] 
                   
                   ⁡ 
                   
                     [ 
                     
                       
                         
                           - 
                           K 
                         
                         · 
                         t 
                       
                       C 
                     
                     ] 
                   
                 
                 . 
               
             
           
         
       
     
     The slope α k  of the tangent to the curve representative of the evolution of the quantity T ik −T ek  as a function of time is then: 
     
       
         
           
             
               α 
               k 
             
             = 
             
               
                 
                   P 
                   
                     tot 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     k 
                   
                 
                 C 
               
               - 
               
                 
                   ( 
                   
                     
                       
                         T 
                         ik 
                       
                       ⁡ 
                       
                         ( 
                         0 
                         ) 
                       
                     
                     - 
                     
                       
                         T 
                         ek 
                       
                       ⁡ 
                       
                         ( 
                         0 
                         ) 
                       
                     
                   
                   ) 
                 
                 ⁢ 
                 
                   
                     K 
                     C 
                   
                   . 
                 
               
             
           
         
       
     
     In practice, within the framework of the method of the invention, within the time period D k , a time interval Δt k  is sought for which the evolution T ik (t) is substantially linear. Over this time interval Δt k , it may be considered that the temperature of the outside air T ek  is substantially constant and equal to the mean temperature over the time interval Δt k , denoted T ekm . Furthermore, as a time interval Δt k  is selected in the heating period D k , the exact positioning of the time t=0 with respect to the heating period is arbitrary and it is preferable to consider a mean value T ikm  of T ik (t) over the time interval Δt k  in the expression for the slope α k . Hence, the slope α k  of the curve θ k (t)=T ik (t)−T ekm  over the time interval Δt k , which is equal to the slope of the curve T ik (t) over the time interval Δt k , is, in this linear approximation: 
     
       
         
           
             
               
                 α 
                 k 
               
               = 
               
                 
                   
                     P 
                     
                       tot 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       k 
                     
                   
                   C 
                 
                 - 
                 
                   
                     θ 
                     km 
                   
                   ⁢ 
                   
                     K 
                     C 
                   
                 
               
             
             , 
             
               
                 with 
                 ⁢ 
                 
                     
                 
                 ⁢ 
                 
                   θ 
                   km 
                 
               
               = 
               
                 
                   T 
                   ikm 
                 
                 - 
                 
                   
                     T 
                     ekm 
                   
                   . 
                 
               
             
           
         
       
     
     If the value of the effective heat capacity of the premises C is known, the value of the heat loss coefficient K of the premises is thus derived directly. 
     If the value of the effective heat capacity of the premises C is not known, the value of the heat loss coefficient K of the premises can be derived by applying over two successive time periods, two heating powers P imp1  and P imp2  of different values, and by measuring the evolution of at least one temperature inside the premises T i1 (t) or T i2 (t) over each of these two time periods. The value of the coefficient K can then be obtained by selecting a time interval Δt 1  or Δt 2  for which the evolution T i1 (t) or T i2 (t) is substantially linear, and by determining, over this time interval Δt 1  or Δt 2 , the slope α 1  or α 2  of the tangent to the curve (T ik (t)) k=1 or 2 . The ratio of the slopes 
               α   1       α   2           
leads to the value of the heat loss coefficient K of the premises by removing the dependency with regard to the effective heat capacity of the premises C:
 
     
       
         
           
             
               
                 
                   
                     
                       α 
                       1 
                     
                     
                       α 
                       2 
                     
                   
                   = 
                   
                     
                       
                         
                           P 
                           
                             tot 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             1 
                           
                         
                         - 
                         
                           
                             θ 
                             
                               1 
                               ⁢ 
                               m 
                             
                           
                           · 
                           K 
                         
                       
                       
                         
                           P 
                           
                             tot 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             2 
                           
                         
                         - 
                         
                           
                             θ 
                             
                               2 
                               ⁢ 
                               m 
                             
                           
                           · 
                           K 
                         
                       
                     
                     . 
                   
                 
               
               
                 
                   ( 
                   2 
                   ) 
                 
               
             
           
         
       
     
     More precisely, the method for determining the heat loss coefficient K of the premises then comprises steps in which:
         in the unoccupied premises and over two successive time periods D 1  and D 2  are performed:       

     i. over the first time period D 1 , the application of a first imposed heating power P imp1  of the premises by means of a controlled power source, and a campaign of measurements of at least one temperature inside the premises T i1  at closely-spaced time intervals, as well as the determination of the temperature of the outside air T e1  at said same closely-spaced times, and then 
     ii. over the second time period D 2 , the application of a second imposed heating power P imp2  of the premises by means of a controlled power source, where the second imposed heating power P imp2  is different from the first power P imp1 , and a campaign of measurements of at least one temperature inside the premises T i2  at closely-spaced time intervals, as well as the determination of the temperature of the outside air T e2  at said same closely-spaced times;
         for each of the first and second time periods D 1  and D 2 , a time interval Δt 1  or Δt 2  is selected for which the evolution T i1 (t) or T i2 (t) is substantially linear and the slope α 1  or α 2  of the tangent to the curve (T ik (t)) k=1 or 2  over this time interval Δt 1  or Δt 2  is determined;   the value of the heat loss coefficient K of the premises is deduced on the basis of the ratio of the slopes       

     
       
         
           
             
               
                 α 
                 1 
               
               
                 α 
                 2 
               
             
             . 
           
         
       
     
     According to a feature of the invention, the difference between the imposed heating powers P imp1  and P imp2  is maximized. In an advantageous manner, one from among the powers P imp1  and P imp2  is a zero power, whilst the other power is a strictly positive power making it possible to obtain a variation of the inside temperature T ik  of at least 1° C. over a time interval Δt k  for which the evolution T ik (t) is substantially linear. Over the time period corresponding to the application of the non-zero imposed heating power, the rise in temperature of the premises as a function of time is measured. Over the time period corresponding to the application of the zero imposed heating power, this corresponding to an absence of heating in the premises, the fall in temperature of the premises as a function of time is measured. The heat loss coefficient K of the premises is then determined on the basis of the ratio of the slopes of the temperature rise and fall curves of the premises. 
     Preferably, for each time period D k , the campaign of measurements of the inside temperature of the premises T ik  is carried out over a time period sufficient to obtain a variation of the inside temperature T ik  of at least 1° C., preferably of between 1° C. and 10° C. 
     By way of example, it is possible to perform successively the application of the zero first imposed heating power P imp1 , corresponding to an absence of heating in the premises, and then the application of the non-zero second imposed heating power P imp2 . 
     As a variant, the non-zero imposed heating power, which generates a temperature rise, may be applied prior to the zero imposed heating power, that is to say the stopping of the heating in the premises, which generates a temperature fall. 
     According to an advantageous feature, the controlled power source for the heating of the premises may be a fixed item of equipment of the premises, that is to say a heating means installed in the premises independently of the implementation of the method, provided that this heating means is sufficiently powerful to ensure fast heating of the premises and that the power delivered can be measured accurately. 
     This may in particular entail a heat pump whose coefficient of performance (COP) is known. The COP, which is the ratio of the thermal power produced to the electrical power consumed, is customarily of the order of 3 to 5. The COP varies with the temperatures of the cold and hot sources. Now, whereas the hot source is regulated and exhibits a substantially constant temperature, the cold source is for its part generally taken outside and is therefore not controllable. The determination of the temperature of the outside cold source over a time period encompassing the period of implementation of the method of the invention makes it possible to adjust the value of the COP more accurately. 
     As a variant, the controlled power source for the heating of the premises may be a source brought into the premises specifically for the implementation of the method. 
     The heating elements for the premises may be of convective, conductive or radiative type, or may combine several of these technologies. Preferably, the heating elements are electrical appliances, thereby making it possible to determine the heating power directly and accurately. Examples of electrical heating appliances comprise in particular appliances of convective type involving the blowing of air heated by means of electrical resistances; heating blankets or heating films, in particular deposited against the walls of the premises; parasol radiant heaters. As a variant, the heating elements may be appliances that operate using gas or fuel oil, provided that the efficiencies of the burners and the fuel flow rates can be estimated sufficiently accurately to derive the heating power. 
     In an advantageous manner, when the heating elements for the premises are of convective type, it is possible to combine these heating elements, which ensure heat dissipation, with at least one fan, which ensures good spatial distribution of the heat. Such a heating assembly combining heating elements of convective type and fans makes it possible to guarantee homogeneous heating of the premises. 
     Preferably, the source for the heating of the premises is chosen to allow heating which ensures that the mean temperature of the walls inside the premises is substantially equal to the temperature of the ambient air inside the premises. In practice, this condition is fulfilled when the heating elements, fixed or brought into the premises specifically for the method, heat the mass of the premises directly, that is to say the heat capacity of the premises, and not only the air inside the premises. Direct heating of the mass of the premises, which may be obtained for example by means of a fixed floor-based heating system of the premises, or else by means of heating films added onto the floor of the premises, is therefore generally preferable to heating of convective type which heats the air inside the premises by priority. 
       FIGS. 1 and 2  annexed illustrate the advantages, for obtaining homogeneous heating of a premises, of the use of a heating system which heats the mass of the premises directly, rather than a heating system of convective type. These figures show the curves of air temperature rise in various rooms of a 5-room two-storey house, respectively in the case of heating by means of electric convectors for  FIG. 1  and in the case of heating by means of electrical heating films for  FIG. 2 . 
     In the case of  FIG. 1 , an electric convector is placed at the center of each measured room of the house, with a temperature setting of 31° C. In the case of  FIG. 2 , the floor of each measured room of the house is covered with heating films of 150 W/m 2  power, marketed by the company DOMOTECK in the “Aluminium Mat” range, with a temperature setting of 31° C. The results of  FIGS. 1 and 2  correspond to one and the same total heating power dissipated in the house, respectively by the electric convectors for  FIG. 1  and by the heating films for  FIG. 2 . For  FIGS. 1 and 2 , the temperature measurements are carried out, in each measured room of the house, with the aid of a thermocouple positioned in the ambient air at the center of the room at a height of 160 cm. 
     As emerges from a comparison of  FIGS. 1 and 2 , heating using heating films makes it possible to achieve better temperature homogeneity in all the rooms of the house than convective heating. In particular, with heating using heating films, fewer disparities are observed between the temperatures of the first-floor rooms, on the one hand, and the temperatures of the second-floor rooms, on the other hand. It is also noted that the slopes of the temperature rise curves are more homogeneous from one room to another with the heating using heating films, this being particularly sought-after within the framework of the method according to the invention. These two advantages, namely better temperature homogeneity and better homogeneity in the temperature rise slopes, originate from the fact that the mass of the house is heated directly in the case of heating using heating films. 
     According to an advantageous feature, each campaign of measurements of the temperature inside the premises comprises measurements of the ambient temperature inside the premises, measurements of the temperature of walls of the premises and/or measurements of the radiant mean temperature inside the premises. In practice, any known measurement method may be used to derive these temperatures, in particular the measurement methods described in the standard NF EN ISO 7726. By way of example, the measurements of the ambient temperature inside the premises and of the temperature of the walls of the premises may be carried out with the aid of thermocouples or, preferably, Pt100 probes. For the measurements of the radiant mean temperature inside the premises, a black globe thermometer may advantageously be used. 
     Preferably, at least the temperature of a wall inside the premises is measured. When the heating of the premises ensures that the ambient temperature is sufficiently close to the temperature of the walls inside the premises, which may be obtained in particular with a heating which heats the mass of the premises by priority, a measurement of the ambient temperature inside the premises may be substituted for the measurement of the wall temperature inside the premises. 
     If the heating of the premises is very homogeneous, so that the inside temperature is the same throughout the premises, or in all the rooms of the premises if the latter comprises internal partitions, then the measurements of the temperature inside the premises may be limited to measurements inside a single room of the premises. 
     If the method of the invention is implemented in a premises for which the heating is less homogeneous, it may be envisaged to measure the temperature in several rooms of the premises and to consider that the temperature inside the premises at each time t is the average of the temperature measurements obtained at the time t in the various rooms of the premises, on condition that they are not too different, which would indicate a lack of ventilation of the premises. Provision may also be made for several different temperature measurements in each room of the premises. Thus, provision may be made to carry out in each room at the same time a measurement of the ambient temperature and/or a measurement of the temperature of a wall of the envelope of the premises and/or a measurement of the radiant mean temperature. 
     According to an advantageous feature, the determination of the temperature of the outside air T ek  takes place, within the framework of the method of the invention, by way of a campaign of measurements which are simultaneous with the measurements of the inside temperature of the premises T ik , that is to say at the same closely-spaced times. 
     As a variant, the determination of the temperature of the outside air T ek  at these closely-spaced times may be obtained by interpolation of meteorological data local to the premises. 
     Preferably, the method according to the invention is implemented over a time period for which the temperature of the outside air T ek  is stable. 
     When the variation of the temperature of the outside air T ek (t) is significant during the implementation of the method, this variation can be taken into account by regarding it as a linear evolution, thus resulting in a non-linearity of the evolution of the inside temperature of the premises T ik (t). The curve θ k (t)=T ik (t)−T ek (t) can then be modeled as a second-degree polynomial whose coefficients involve the characteristics K and C of the premises. In this case, it is possible to determine K and C by applying a single imposed heating power P impk  of the premises, with a single campaign of measurements of the temperature inside the premises T ik , since, by expanding the solution of equation (1) to second order, equality is obtained between two second-degree polynomials, with equality for the term in t and equality for the term in t 2 . Such an expedient is however less preferred, because of its greater complexity. 
     As explained previously, the heat loss coefficient K determined in accordance with the invention integrates the contributions of the heat losses by transmission and infiltrations of air, that is to say:
 
 K=H   T   +m′.Cp=UA   T   +m′.Cp.  
 
     If it is desired to derive the heat transfer coefficient U of the premises, it is possible to decouple the contribution of the heat losses by transmission, on the one hand, and that of the air infiltrations, on the other hand, by evaluating the renewal flow rate of the air m′ in the premises. 
     When no fixed ventilation system of the premises is active in the course of the method for determining the coefficient K, the flow rate m′ is equal to the renewal flow rate of the air by infiltration. This flow rate m′ may be determined by any appropriate scheme, in particular by a detection scheme based on tracer gases or by a blower door infiltrometry test. 
     In an advantageous manner, the detection scheme based on tracer gases provides an instantaneous value of the flow rate m′. This scheme based on tracer gases also makes it possible to take accurate account of the contribution of the ventilation in the energy balance of the premises, this being advantageous in cases where the method is implemented whilst a ventilation system fitted to the premises is active, in particular in cases where it is not possible to deactivate, prior to the implementation of the method, the fixed ventilation systems fitted to the premises. 
     If one opts for a blower door test, it is possible to compute a mean value of the flow rate m′ on the basis of the measured value of the leakage flow rate. 
     By way of example, if a blower door test is carried out in which the indicator n50 is measured such as defined in the standard NF EN 13829, that is to say the leakage flow rate under 50 Pa divided by the heated volume of the premises, it is known to deduce the value of the air renewal flow rate m′ on the basis of the measured value of the n50, by using the rule established by Drubul: 
                       m   ′     =       n   ⁢           ⁢     50   ·   V       20       ,           (   3   )               
where V is the heated volume of the premises.
 
     As a variant, the mean value of the air renewal flow rate m′ may be obtained on the basis of the measurement of indicators other than the indicator n50, or else with the aid of regulatory computation schemes other than that involving relation (3) hereinabove. In particular, an alternative approach to Drubul&#39;s rule for determining the value of the air renewal flow rate m′ on the basis of the measured value of the n50 consists in estimating the air infiltration rate based on the empirical model proposed in the standard NF EN ISO 13790, annex G. 
     When a fixed ventilation system of the premises is active in the course of the method for determining the coefficient K, the flow rate m′ incorporates the contribution of air renewal due to this ventilation system, that ought to be taken into account in addition to the air renewal by infiltration. The air renewal flow rate due to the fixed ventilation system may be determined by measuring the air flow rate at each ventilation vent. 
     When the trials are carried out by preferentially heating the mass rather than the interior air of the premises, or when the heat loss coefficient is computed with the aid of wall temperatures rather than temperatures of the ambient air, the computed loss coefficient K calc  is not perfectly representative of the actual losses, since it takes no account of the heat exchange between the wall and the ambient air. This may be partially corrected by determining the convecto-radiative exchange coefficient h i . The coefficient h i  may be estimated through the knowledge of the ambient and black globe inside temperatures of the premises, as well as of the air speed, on the basis of the following relation, defined in RT 2005: h i =2.5+4.σ.ε i .T mi   3 , where T mi  is the radiant mean temperature which, in accordance with the standard NF EN ISO 7726, is related, for a standard globe 85 mm in diameter and with emissivity 0.95, to the globe temperature T g , the ambient temperature T α  and the speed of the air v α :
 
 T   mi =[( T   g +273) 4 +2.5×10 8   ×v   α   0.6 ( T   g   −T   α )] 1/4 −273.
 
     As a variant, if the radiant mean temperature T mi  is not available, the convective exchange coefficient h i  may be assumed known. In particular, the convective exchange coefficient h i  may be estimated as being of the order of 8 W/m 2 .K, which is the order of magnitude given by RT 2005 in the Th-U rules. 
     The corrected loss coefficient K corr  then equals: 
     
       
         
           
             
               
                 
                   
                     1 
                     
                       K 
                       corr 
                     
                   
                   = 
                   
                     
                       1 
                       
                         K 
                         calc 
                       
                     
                     + 
                     
                       
                         1 
                         
                           
                             h 
                             i 
                           
                           ⁢ 
                           
                             A 
                             T 
                           
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   4 
                   ) 
                 
               
             
           
         
       
     
     One subject of the invention is also a device for the implementation of a method such as described hereinabove, comprising at least one temperature sensor which measures a temperature inside the premises T ik  and a heating device for homogeneous heating of the premises comprising a controlled power source. 
     According to an advantageous feature of such a device, the heating device heats the mass of the premises directly, that is to say the heat capacity of the premises and not only the air inside the premises, and the temperature sensor measures the temperature in the air inside the premises. As mentioned previously, examples of heating devices which heat the mass of the premises directly comprise heating films added on the floor of the premises or else fixed floor-based heating systems of the premises. 
     The selection of a device which combines direct heating of the mass of the premises and measurement of the temperature in the air inside the premises is particularly advantageous for the implementation of the method according to the invention. Indeed, if it is chosen to carry out measurements of wall temperature, it is necessary, in order to obtain a good estimation of the mean temperature in a room, to carry out measurements of temperatures on several walls of the room, and then to determine an average of these temperatures of walls, and to do so whatever the homogeneity of the heating. Conversely, if it is chosen to carry out measurements of air temperature, a single measurement in the air volume substantially at the center of the room suffices to obtain a value representative of the mean temperature in the room, on condition that the heating is sufficiently homogeneous. Measurement of temperature in the air therefore makes it possible to reduce the number of measurements to be performed within the framework of the method according to the invention, provided that it is carried out in a room heated in a globally homogeneous manner. Now, as explained previously, homogeneous heating is obtained more easily and more rapidly with a heating device which heats the mass of the premises directly, rather than with a convective heating device which heats the air inside the premises by priority. Hence, by combining direct heating of the mass of the premises and measurement of the temperature in the air inside the premises, the implementation of the method according to the invention is simplified and its duration is limited. 
     According to another advantageous feature, the device furthermore comprises an electronic central unit comprising means for acquiring the temperature measurements inside the premises, means for computing the heat loss coefficient K of the premises on the basis of the acquired temperature measurements and means of automatic control of the power source as a function of the acquired temperature measurements. 
     An automatic control device such as this uses the temperature signals acquired, not only to compute the heat loss coefficient K of the premises, but also to automatically drive the power source ensuring the heating of the premises. In an advantageous manner, the processing of the temperature information is ensured by the electronic central unit, which is associated with control software internal to the device which is parametrized in such a way that any temperature evolution corresponding to predefined criteria of the method according to the invention, in particular criteria describing the linearity of the evolution T ik (t) over a time interval, the stability of the temperature of the outside air T ek , etc., makes it possible to control the intensity of the heating power in the premises. 
     By way of example, a test that might be run in an autonomous manner by an automatic control device such as this comprises the following sequence of steps:
         running of the procedure;   turning on of the heating device;   when the evolution T i1 (t) of the temperature measured inside the premises in the course of the heating of the premises is substantially linear and the temperature of the outside air T e1  is stable with regard to the predefined criteria, computation and storage of the value of the slope α 1  of the tangent to the heating curve T i1 (t) and turning off of the heating device;   when the evolution T i2 (t) of the temperature measured inside the premises in the course of the cooling of the premises is substantially linear and the temperature of the outside air T e2  is stable with regard to the predefined criteria, computation and storage of the value of the slope α 2  of the tangent to the cooling curve T i2 (t);   on the basis of the stored values of the slopes α 1  and α 2 , computation of the heat loss coefficient K of the premises.       

     The heating device of the automatic control device may be a heating device tied to the premises tested or a heating device added specifically for carrying out the trials. Likewise, the temperature measurement sensors of the automatic control device may be tied to the premises or add-ons. As explained previously, the temperature measurement sensors are preferably sensors for measuring the temperature in the air inside the premises, rather than sensors for measuring the temperature of walls inside the premises. 
     According to an advantageous feature, the control software internal to the device is designed in such a way that each test run in an autonomous manner by the automatic control device is run preferably at night, and, if previous measurements have already been carried out, the automatic control device optimizes the cycle of each test so as to minimize its duration and to maximize the accuracy of characterization. 
     Preferably, the criteria of the method that are used by the automatic control device take account of the accuracy of the temperature measurement, that is to say of the accuracy in determining the slope of the temperature evolution. The lower the accuracy of the measurement, the longer the measurement time must be to ensure correct slope determination. During the cooling phases, the criterion must be defined in such a way that the absolute value of the slope has a maximum value, any quasi-zero slope having to be proscribed. 
     In an advantageous manner, in the case of a device using heating systems and sensors internal to the premises, it may be envisaged to optimize the driving of the heating by the automatic control device as a function of the occupancy of the premises and of its thermal performance K and C. 
     One subject of the invention is also a medium for recording information, comprising instructions for the implementation, within the framework of a method such as described hereinabove:
         of the steps of computing the heat loss coefficient K of the premises on the basis of acquired temperature measurements, and   of the steps of automatic control of the power source as a function of the acquired temperature measurements,
 
when these instructions are executed by an electronic computation unit.
       

    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         FIGS. 1 and 2  graphically illustrate the advantages, for obtaining homogeneous heating of a premises, of the use of a heating system which heats the mass of the premises directly; 
       The features and advantages of the invention will be apparent in the following description of several embodiments of a method and of a device according to the invention, given solely by way of example and made with reference to the appended  FIGS. 1 to 8  in which: 
         FIG. 3  is a schematic view of an individual house whose heat loss coefficient K it is desired to determine in accordance with the invention, this house being equipped with a heat pump as controlled power source which feeds a floor-based heating system; 
         FIG. 4  is a curve representative of the evolution of the temperature T ik  inside the house of  FIG. 3  as a function of time, in the course of the implementation of the method in accordance with the invention; 
         FIG. 5  is a curve representative of the evolution of the temperature T i1  inside the house of  FIG. 3  as a function of time, for a first time period D 1  in the course of which a zero imposed heating power P imp1  is applied in the house, that is to say this first time period corresponds to an absence of heating in the house; 
         FIG. 6  is a curve representative of the evolution of the temperature T i2  inside the house of  FIG. 3  as a function of time, for a second time period D 2  in the course of which a non-zero imposed heating power P imp2  is applied in the house; 
         FIG. 7  is a schematic view of a bungalow whose heat loss coefficient K it is desired to determine in accordance with the invention, the heating of the bungalow being ensured by electric convectors; and 
         FIG. 8  is a curve representative of the evolution of the temperature T ik  inside the bungalow of  FIG. 7  as a function of time in the course of the implementation of the method in accordance with the invention, showing a first time period D 1  in the course of which a zero imposed heating power P imp1  is applied in the bungalow, that is to say this first time period corresponds to an absence of heating in the bungalow, followed by a second time period D 2  in the course of which a non-zero imposed heating power P imp2  is applied in the bungalow. 
     
    
    
     DETAILED DESCRIPTION 
     Example 1 
     With reference to  FIG. 3 , the method according to the invention is implemented for the determination of the heat loss coefficient K of an individual house  1  of recent construction, equipped with a heat pump  2 . The heat pump  2  feeds a floor-based heating system  3  allowing homogeneous heating of the house. In particular, the floor-based heating system  3  ensures heating of the house  1  such that the temperature of the walls inside the house is substantially equal to the temperature of the ambient air inside the house. The heat loss coefficient K of the house  1  may be sought within the framework of a diagnosis of the energy performance of the house, for example to verify that the house  1  satisfies certain quality labels in terms of thermal insulation, such as the BBC label or the Passivhaus label. 
     The heating power provided by the heat pump  2  for the heating of the house is readily determinable, in particular on the basis of the COP of the heat pump as explained previously. The heat pump  2  therefore constitutes a controlled power source that is well adapted for generating the heating pulses for the house, that is to say the imposed heating powers P impk , which are required by the method of the invention. In this example, the value of the COP of the heat pump  2  is equal to 4.23. 
     The method for the determination of the heat loss coefficient K of the house  1  is implemented whilst the house is unoccupied. Moreover, the house  1  is equipped with a fixed ventilation system comprising two hygro-adjustable single flow controlled mechanical ventilations (CMVs), which are not deactivated during the implementation of the method. However, as the house  1  is unoccupied and the measurement periods of the method are relatively short, it may be considered that these CMVs do not work in the course of the method. 
     As mentioned previously, the method according to the invention is preferably implemented continuously in its entirety over a single time period. This is the case for the example represented in  FIGS. 4 to 6  since, as clearly visible in  FIG. 4 , the second time period D 2  in the course of which a non-zero imposed heating power P imp2  is applied in the house  1  immediately follows the first time period D 1  in the course of which a zero imposed heating power P imp1  is applied in the house. 
     In the example of  FIGS. 4 to 6 , the method takes place over a continuous time period which lasts about 1500 minutes, and which begins at nightfall and extends into the following daytime. Over this time period of 1500 minutes, the recorded solar radiation is low since the weather was cloudy during the implementation of the method. Furthermore, all the shutters of the house were closed during the implementation of the method. Under these conditions, the contribution of the solar radiation to the heating of the house  1  may be neglected. 
     Moreover, no power source other than that used for the application of the heating pulses is active in the house  1  in the course of the method. Thus, for each time period D 1  and D 2 , the only power to be considered in the energy balance is the imposed heating power P imp1  or P imp2 . 
     In a first step of the method, which corresponds to the first time period D 1 , the zero first imposed heating power P imp1  is applied in the house  1 , commencing from a starting temperature T i1d =17° C., that is to say the heating system  3  does not operate during this period D 1 . The ambient temperature inside the house T i1  is then measured every ten minutes in two different rooms of the house, namely a main room and a bedroom. For this purpose, a temperature sensor is installed in each of these two rooms, in the ambient air at a height of 180 cm. 
     In this example, the measured evolution of the ambient temperature inside the main room and the measured evolution of the ambient temperature inside the bedroom are almost identical. Indeed, the heating of the house  1  is particularly homogeneous, so that the ambient temperature is the same in all the rooms of the house. In  FIGS. 4 to 6 , only the evolution of the ambient temperature inside the main room is represented, it being understood that the evolution of the ambient temperature inside the bedroom has a similar profile. 
     The curve representative of the evolution of the inside temperature of the house T i1  as a function of time during the time period D 1  is shown in  FIG. 5 . As visible in this figure, the temperature fall curve for the house  1  exhibits a substantially linear part over the time interval Δt 1 . Fitting an equation to this linear part of the curve gives: T i1 =16.5° C.−0.003 (t−800), with t in minutes. 
       FIG. 5  also depicts the evolution of the temperature of the outside air T e1  during the time period D 1 . The temperature of the outside air T e1  over the time interval Δt 1  is sufficiently stable for it to be possible to consider it substantially constant and equal to the mean temperature over the time interval Δt 1 , namely in this example T e1m =5° C. 
     In a second step of the method, which corresponds to the second time period D 2 , the heating in the house  1  is restored by applying the second imposed heating power P imp2  in the house  1 , equal to the maximum power delivered by the heat pump  2 , i.e. P imp2 =COP×5000 W=4.23×5000 W, commencing from a starting temperature T i2d =15° C. As in the first step, the ambient temperature inside the house T i2  is then measured every ten minutes, by means of two temperature sensors installed respectively in the main room and in the bedroom of the house, each time in the ambient air at a height of 180 cm. Here again, the measurements show that the evolution of the ambient temperature inside the main room and the evolution of the ambient temperature inside the bedroom are almost identical. 
       FIG. 6  shows the curve representative of the evolution of the inside temperature of the house T i2  as a function of time during the time period D 2 . As visible in this figure, the temperature rise curve for the house  1  exhibits a substantially linear part over the time interval Δt 2 . Fitting an equation to this linear part of the curve gives: T i2 =16.2° C.+0.014 (t−1730), with t in minutes. 
     The evolution of the temperature of the outside air T e2  during the same time period D 2  is also shown in  FIG. 6 . As in the first step, the temperature of the outside air T e2  over the time interval Δt 2  is sufficiently stable for it to be possible to consider it substantially constant and equal to the mean temperature over the time interval Δt 2 , namely in this example T e2m =14° C. 
     The ratio 
               α   1       α   2           
of, on the one hand, the slope α 1  of the straight line representative of the evolution of the quantity T i1 (t) to, on the other hand, the slope α 2  of the straight line representative of the evolution of the quantity T i2 (t) therefore equals
 
     
       
         
           
             
               
                 α 
                 1 
               
               
                 α 
                 2 
               
             
             = 
             
               - 
               
                 
                   0.003 
                   0.014 
                 
                 . 
               
             
           
         
       
     
     As the ratio of the slopes 
               α   1       α   2           
is equal to
 
                 P     tot   ⁢           ⁢   1       -       θ     1   ⁢           ⁢   m       ·   K           P     tot   ⁢           ⁢   2       -       θ     2   ⁢   m       ·   K             
according to the previous equation (2), by taking θ 1m =10.9° C., θ 2m =3.9° C., P imp1 =0 W, P imp2 =4.23×5000 W, the value of the heat loss coefficient K of the house  1  is obtained:
 
K calc =386 W/K.
 
     In practice, the steps of acquiring the inside temperature T ik  and outside temperature T ek  and of computing the heat loss coefficient K on the basis of the acquired temperatures may be carried out by means of a microprocessor or any other appropriate electronic computation unit. 
     The heating having been supplied via the mass of the building, relation (4) is applicable and, with h i =8 W/m 2 .K and A T =736 m 2 , we obtain:
 
K corr =362 W/K.
 
     The value of the heat loss coefficient K corr  of the house  1  determined by the method of the invention may be compared with a mean static value K s  of the heat loss coefficient. This mean static value K s  is estimated on the basis of the energy consumption of the heat pump  2 , which is measured once a week for thirteen winter weeks, in the occupied configuration of the house and for an imposed static temperature inside the house T is  of 20° C. By deducting the share of consumption due to the sanitary hot water, the mean static value K s  of the heat loss coefficient is estimated at about 430 W/K. This mean static value K s  is indeed of the same order of magnitude as the value of the heat loss coefficient K corr  determined by the method of the invention. 
     The difference between the two values K s  and K corr  may originate from the different occupancy conditions of the house  1  for, on the one hand, the implementation of the method in accordance with the invention and, on the other hand, the evaluation of the mean static value K s . Indeed, the method of the invention has been implemented in the unoccupied house, while the mean static value K s  has been evaluated with the house occupied by four people, thus implying additional energy contributions, in particular in terms of air renewal flow rate due to the hygro-adjustable CMVs, and in terms of power due to the solar radiation or to the operation of domestic electrical appliances. 
     A blower door test has also been carried out on the house  1 . This test leads to a value of the indicator n50 of 7.35 for a heated volume of the house  1  of 688 m 3 . With the aid of relation (3), it is then possible to deduce the mean value of the air renewal flow rate m′ in the house  1 :
 
m′=250 m 3 /h.
 
     It follows from this that the contribution of the air infiltrations represents:
 
m′.Cp=85 W/K.
 
     By assuming that the two hygro-adjustable CMVs of the house  10  did not work in the course of the time periods D 1  and D 2 , it may be deduced that the contribution of the heat losses by transmission is:
 
 H   T   =UA   T   =K   corr   −m′.Cp= 277 W/K.
 
     Furthermore, the method of the invention makes it possible to derive the value of the inertia or effective heat capacity C of the house  1  as defined previously, namely the energy necessary to increase the ambient temperature of the house by 1 K at constant outside temperature:
 
C=80 MJ/K.
 
     Example 2 
     With reference to  FIGS. 7 and 8 , the method according to the invention is implemented for the determination of the heat loss coefficient K of a bungalow  10  which has an interior volume V i  of 5.8 m×2.1 m×2.6 m and which is fitted with two triple-glazed windows. The envelope of the bungalow  10  consists of insulating sandwich panels assembled via a metallic structure. An additional insulation has been added to the envelope, in the form of a 40 mm thickness of glass wool and of a plasterboard as interior lining of the envelope. The permeability to air, measured by tracer gas, is 0.43 vol/h. The method is implemented whilst the bungalow is unoccupied. 
     The heating of the bungalow  10  is ensured by electric convectors  20  having a measured actual power of 1880 W. The convectors make it possible to heat the air in the bungalow and, on account of the limited volume of the bungalow, the heating of the bungalow is homogeneous. The convectors  20  constitute a controlled power source that is well adapted for generating the heating pulses for the bungalow, that is to say the imposed heating powers P impk , which are required by the method of the invention. During heating, the setting is 32° C., measured by a black globe placed at the center of the air volume. 
     The method according to the invention is implemented continuously in its entirety over a single nocturnal time period, so as to circumvent the contribution of the solar radiation to the heating of the bungalow  10 . Commencing from an inside temperature of the bungalow T i1d =32° C., a cooling of the bungalow is undertaken over a first time period D 1  from 11 pm to 3 am, this corresponding to a zero imposed heating power P imp1 , and then a heating of the bungalow is undertaken over a second time period D 2  from 3 am to 8 am with an imposed heating power P imp2  of 1880 W. The second time period D 2  thus immediately follows the first time period D 1 . 
     No power source other than that used for the application of the heating pulses is active in the bungalow  10  in the course of the method. Thus, for each time period D 1  and D 2 , the only power to be considered in the energy balance is the imposed heating power P imp1  or P imp2 . 
     In a first step of the method, which corresponds to the first time period D 1 , the zero first imposed heating power P imp1  is applied in the bungalow  10 , commencing from the starting temperature T i1d =32° C., that is to say the convectors  2  do not operate during this period D 1 . The ambient temperature inside the bungalow T i1  is then measured every ten seconds. For this purpose, a black globe is placed at the center of the bungalow. 
     The curve representative of the evolution of the inside temperature of the bungalow T i1  as a function of time during the time period D 1  is shown in  FIG. 8 . As visible in this figure, the temperature fall curve of the bungalow  10  exhibits a substantially linear part over the time interval Δt 1 . The fitting of an equation to this linear part of the curve gives: T i1 =26.5° C.−0.00053 t, with t in seconds. 
       FIG. 8  also depicts the evolution of the temperature of the outside air T e1  during the time period D 1 . The temperature of the outside air T e1  over the time interval Δt 1  is sufficiently stable for it to be possible to consider it substantially constant and equal to the mean temperature over the time interval Δt 1 , namely in this example T e1m =13.4° C. 
     In a second step of the method, which corresponds to the second time period D 2 , the heating in the bungalow  10  is restored by applying the second imposed heating power P imp2  in the bungalow, equal to 1880 W, commencing from a starting temperature T i2d =22° C. As in the first step, the ambient temperature inside the house T i2  is then measured every ten seconds, by means of a black globe placed at the center of the bungalow. 
       FIG. 8  shows the curve representative of the evolution of the inside temperature of the bungalow T i2  as a function of time during the time period D 2 . As visible in this figure, the temperature rise curve of the bungalow  10  exhibits a substantially linear part over the time interval Δt 2 . The fitting of an equation to this linear part of the curve gives: T i2 =28.9° C.+0.00179 t, with t in seconds. 
     The evolution of the temperature of the outside air T e2  during the same time period D 2  is also shown in  FIG. 8 . As in the first step, the temperature of the outside air T e2  over the time interval Δt 2  is sufficiently stable for it to be possible to consider it substantially constant and equal to the mean temperature over the time interval Δt 2 , namely in this example T e2m =11.8° C. 
     The ratio 
               α   1       α   2           
of, on the one hand, the slope α 1  of the straight line representative of the evolution of the quantity T i1 (t) to, on the other hand, the slope α 2  of the straight line representative of the evolution of the quantity T i2 (t) therefore equals
 
     
       
         
           
             
               
                 α 
                 1 
               
               
                 α 
                 2 
               
             
             = 
             
               - 
               
                 
                   0.00053 
                   0.00179 
                 
                 . 
               
             
           
         
       
     
     As the ratio of the slopes 
               α   1       α   2           
is equal to
 
                 P     tot   ⁢           ⁢   1       -       θ     1   ⁢   m       ·   K           P     tot   ⁢           ⁢   2       -       θ     2   ⁢   m       ·   K             
according to the previous equation (2), by taking θ 1m =11.7° C., θ 2m =18.1° C., P imp1 =0 W, P imp2 =1880 W, the value of the heat loss coefficient K calc  of the bungalow  10  is obtained:
 
K calc =32.6 W/K.
 
     In practice, the steps of acquiring the inside temperature T ik  and outside temperature T ek  and of computing the heat loss coefficient K on the basis of the acquired temperatures may be carried out by means of a microprocessor or any other appropriate electronic computation unit. 
     In this example, the interior air of the bungalow  10  is heated and, in order to compute the heat loss coefficient, use is made of black globe inside temperatures which are close to the temperatures of the ambient air inside the bungalow. Hence, it is not necessary to correct the computed value of the heat loss coefficient K calc  with relation (4). 
     The value of the heat loss coefficient K calc  of the bungalow  10  determined by the method of the invention may be compared with a mean static value K s  of the heat loss coefficient. This mean static value K s  is estimated under steady conditions. In practice, it is considered that steady conditions are reached when the temperature of the walls of the bungalow is stable, the temperature of the air then being likewise. The estimation of K s  was done in the course of a long trial during which the mean power P m  necessary to maintain the ambient temperature inside the bungalow T i  stable was measured for a period of 8 hours from 11 pm to 7 am, the temperature of the outside air T e  preferably also being stable during this period. The mean static value K s  of the heat loss coefficient is thereafter estimated on the basis of the relation:
 
 P   m   =K   s   A   T ( T   i   −T   e ).
 
     The mean static value K s  thus estimated of the heat loss coefficient is of the order of 32.7 W/K, this being close to the value of the heat loss coefficient K calc  as determined by the method of the invention. This validates the model used. 
     The mean value of the air renewal flow rate m′ in the bungalow is:
 
m′=0.43 V i =13.62 m 3 /h.
 
     It follows from this that the contribution of the air infiltrations represents:
 
m′.Cp=4.6 W/K.
 
     It may be deduced that the contribution of the heat losses by transmission is:
 
 H   T   =UA   T   =K   calc   −m′.Cp= 28.0 W/K.
 
     Furthermore, the method of the invention makes it possible to derive the value of the inertia or effective heat capacity C of the bungalow  10  as defined previously, namely the energy necessary to increase the ambient temperature of the bungalow by 1 K at constant outside temperature:
 
C=720 kJ/K.
 
     By way of comparative example, the method according to the invention was implemented again for the bungalow  10 , but this time the bungalow  10  is fitted with two plasterboards as interior lining of the envelope of the bungalow, instead of a single plasterboard as previously. All the other parameters are identical, in particular the permeability to air measured by tracer gas is still 0.43 vol/h and the method is implemented continuously in its entirety over a single nocturnal time period during which the bungalow is unoccupied. The heating-cooling cycles are the same as previously. 
     The value of the heat loss coefficient K calc  of the bungalow  10  then obtained is:
 
K calc =30.1 W/K.
 
     The contribution of the air infiltrations still being m′.Cp=4.6 W/K, it may be deduced that the contribution of the heat losses by transmission is:
 
 H   T   =UA   T   =K   calc   −m′.Cp= 25.5 W/K.
 
     The value of the inertia or effective heat capacity C of the bungalow fitted with two plasterboards as interior lining, namely the energy necessary to increase the ambient temperature of the bungalow by 1 K at constant outside temperature, is then:
 
C=1071 kJ/K,
 
this corresponding to an increase of the order of 350 kJ/K in the effective heat capacity of the bungalow with respect to the same bungalow fitted with a single plasterboard as interior lining. An estimation of the plaster area added in order to go from one to two plasterboards as interior lining of the bungalow indicates an addition of inertia of 400 kJ/K. Thus, it emerges that the method according to the invention is capable of discriminating K and C.
 
     Example 3 
     In order to verify the validity of the assumptions upon which the method of the invention rests, virtual trials have been carried out with the TRNSYS software on a fictitious house having an inhabited part of 12.10 m×9.90 m×2.50 m and a total loss area S=350 m 2 . More precisely, two series of computations have been performed:
         a first series corresponding to steady conditions obtained with a non-realistic weather file, with no sun and by fixing the outside temperature at 10° C. and the power at 30 kW during the time required for the temperatures to stabilize (400 hours at the maximum);   a second series corresponding to transient conditions obtained with a realistic weather file, namely 3 days in March for a climate matched to the town of Chambéry, by supplying power distributed homogeneously, either via the floor at the surface (not depth-wise), or via the air, of 30 kW from 7 pm to midnight, and of 3 kW during the following 36 hours, after regulation to 19° C. for a day.       

     The two series of trials were carried out without occupancy or internal heat supply. Each time, the infiltrations were firstly not taken into account, and then they were. 
     The steady case makes it possible to obtain the theoretical heat loss coefficient K of the house. By considering heating via the air and no infiltration, we obtain K calc =167 W/K. Under the same conditions in the transient case, by applying the method in accordance with the invention, we obtain K calc =164 W/K, i.e. less than 2% difference with the value obtained in the steady case. This validates that the method according to the invention furnishes a good estimation of the heat loss coefficient K. 
     By repeating the two series of computations under the same conditions, but with supply of heat via the floor instead of the air, we obtain respectively K calc =177 W/K in the steady case, and K calc =181 W/K in the transient case by applying the method in accordance with the invention. By considering that the coefficient of convective exchange h between the walls and the ambient air is 8 W/m 2 .K and by using relation (4), we obtain K corr =166 W/K in the steady case and K corr =170 W/K in the transient case. This validates the fact that the direct heating of the mass of the premises may be used in the absence of infiltrations. 
     Finally, by considering a more realistic case, with supply of heat via the floor and infiltrations of 0.4 vol/h, we obtain K calc =220 W/K in the transient case by applying the method in accordance with the invention. By making the assumption that it is possible to apply the correction of relation (4) described previously even in the presence of infiltrations, to take account of the heat exchange between the walls and the ambient air with the convective exchange coefficient h taken at 8 W/m 2 .K, we obtain K corr =204 W/K at the air level. By taking account of the value of the infiltrations of 0.4 vol/h=120 m 3 /h, the loss by infiltrations is 120 m 3 /h/3600 s×1.2 kg/m 3 ×1 kJ/kg/K=0.04 kJ/s.K=40 W/K. The heat loss coefficient of the envelope alone is then 204 W/K−40 W/K=164 W/K, a value very close to that estimated under steady conditions. This validates the fact that the direct heating of the mass of the premises may be used even in the presence of infiltrations, and that it is possible by measuring these infiltrations to separate the heat losses by transmission and the losses by infiltrations. 
     The invention is not limited to the examples described hereinabove. In particular, as already mentioned, the method according to the invention may be implemented equally with heating means with which the premises are equipped in a fixed manner or with heating means which are brought into the premises specifically for the implementation of the method, as long as the power provided by these heating means for the pulses required by the method can be accurately determined. Thus, in the example of the house  1 , the method according to the invention could have been implemented by deactivating the heating system using the heat pump  2  of the house and by applying the heating pulses by means of a heating device brought into the house, such as a device comprising electrical heating films or else a device combining electrical heating appliances of convective type and fans. 
     As illustrated previously, when the value of the effective heat capacity of the premises C is not known, the method of the invention involves at least two different imposed heating powers P impk  of the premises. Preferably, the difference between the two imposed powers P impk  is maximized. Thus, in the previous examples, a maximum power and a zero power have been selected. As a variant, it would have been possible to choose two non-zero imposed powers P impk , in particular a relatively low heating power and a maximum heating power. 
     Moreover, in the case of a premises of large size, such as an apartment building with multiple floors, the method of the invention may be used either for the determination of the heat loss coefficient K of the premises as a whole, in which case the heating and the intermingling of the air must be ensured throughout the premises, or for the determination of the heat loss coefficient K of just a part of the premises. Thus, in the case of an apartment building, it is possible to test just one apartment in the building. Accordingly, it is necessary:
         either to take into account the heat losses via computation, on condition however that the parts adjoining the measured apartment are in a thermal state representative of their normal occupancy state, in particular that the normally inhabited adjoining parts are at an ambient temperature of the order of 20° C.;   or to minimize the heat losses as far as possible, for example by over-insulating the adjoining walls by means of added insulation, or else by conditioning the adjoining parts in the same manner as the measured apartment so as to ensure a temperature difference on either side of the adjoining wall that is as close to zero as possible.       

     Finally, as emerges from the previous examples, the method according to the invention is very suitable for the determination of the heat loss coefficient K of premises exhibiting good thermal insulation. In this case, it is indeed easy to dissipate a thermal power such that the curve T ik (t) may be considered to be a straight line. For other configurations of premises, in particular old and less well thermally insulated premises, the response time to the heating pulses P impk  may be too short for the evolution of the quantity T ik (t) to be considered linear. The non-linear evolution of T ik (t) can then be modeled by an exponential of the type 
               (         θ   k     ⁡     (   0   )       -       P   impk     K       )     ⁢     exp   ⁡     (       -   t     /   τ     )             
over at least one time interval Δt k ′, where
 
             τ   =     K   C           
is the thermal time constant of the premises. It is possible to determine the coefficient K by applying a single heating pulse for the premises, with a necessarily non-zero power P impk , the scheme consisting thereafter in determining the value K* of the coefficient K such that the curve
 
             Ln   ⁡     [           θ   k     ⁡     (   t   )       -       P   impk       K   *               θ   k     ⁡     (   0   )       -       P   impk       K   *           ]           
is a straight line, where θ k (t)=T ik (t)−T ekm ′ and T ekm ′ is the average of the temperature measurements of the outside air T ek  over the time interval Δt k ′.