Analysis of the methodology to obtain several Key Indicators Performance (KIP), by energy retrofitting of the actual building to the district heating fuelled by biomass, focusing on nZEB Goal: Case of study

energies Article Analysis of the Methodology to Obtain Several Key Indicators Performance (KIP), by Energy Retrofitting of the Actual Building to the District Heating Fuelled by Biomass, Focusing on nZEB Goal: Case of Study Rosaura Castrillón Mendoza 1,2, Javier M. Rey Hernández 2,3,*, Eloy Velasco Gómez 2, Julio F. San JoséAlonso 2and Francisco J. Rey Martínez 2 1Department of Energy and Mechanics, University Autónoma de Occidente Cali (UAO), Cali 760030, Colombia; [email protected] 2Department of Energy and Fluid Mechanics, School of Engineering (EII), University of Valladolid (UVa), 47002 Valladolid, Spain; [email protected] (E.V.G.); [email protected] (J.F.S.J.A.); [email protected] (F.J.R.M.) 3Higher Polytechnic College, European University Miguel de Cervantes (UEMC), 47012 Valladolid, Spain *Correspondence: javier[email protected]; Tel.: +34-983-423-366 Received: 20 November 2018; Accepted: 21 December 2018; Published: 28 December 2018   Abstract: In order to achieve the objectives of the European 20/20/20 strategy, and to obtain a greater energy efficiency, integration of renewable energies and the reduction of carbon emissions, a District Heating (DH) system has been designed by the University of Valladolid (UVa), Spain, one of the most important DH fed by biomass fuel in Spain, supplying heating and domestic hot water (DHW) to 31 buildings in Valladolid, the majority of them, educational buildings on the University Campus. The aims of this paper were to study the change from an energy system fueled by natural gas to District Heating by biomass in a building on the campus of the University of Valladolid—the School of Engineering (EII)—studying its consumption from its connection to the District Heating system. An energy management methodology such as ISO 50001 is carried out, applied to efficiency systems in buildings, thus establishing new criteria of sustainability and economic value. In this paper, energy management will also be analyzed in accordance with the proposed tools of an Energy Management System (EMS) applied to the EII building, through the measurement of energy parameters, calculation of thermal consumption, thermal energy savings as a result of the change from system to District Heating by biomass, economic savings, reduction of environmental impact and indicators of thermal efficiency I 100 and CUSUM indicator. Finally, the primary renewable and non-renewable energy efficiency indicators for the new District Heating system will be determined. The concept of the near Zero Energy Buildings is defined in the European Union (EU) in order to analyze an approach to an nZEB which results from replacing the natural gas heating system by a biomass District Heating system. Keywords: district heating; biomass; energy management in renovated building; nZEB 1. Introduction The spread of university campuses in recent years has led to a significant increase in energy consumption. The total amount of energy consumption includes lighting, DHW, heating or air conditioning systems in all campus buildings. Due to the large amount of energy consumption in these buildings, a policy of promoting energy savings, energy management, supply of energy useful to the final energy consumption, perfectly planned at a local scale, can maintain a standard of energy consumed on campuses in a controlled Energies 2019,12, 93; doi:10.3390/en12010093 www.mdpi.com/journal/energies Energies 2019,12, 93 2 of 20 environment. The University as an institution is one that can best stimulate energy saving among other institutions by setting an example of its actions, pertaining to the consumption that it requires (Figure 1). Energies 2017, 10, x FOR PEER REVIEW 2 of 20 controlled environment. The University as an institution is one that can best stimulate energy saving among other institutions by setting an example of its actions, pertaining to the consumption that it requires (Figure 1). Figure 1. Plan view of connected buildings. Among the different actions carried out by the University of Valladolid to reduce its carbon footprint, is the implementation of a DH system, making it to also the most important biomass DH in Spain, supplying heating and DHW to 31 buildings. The thermal biomass generation plant has a power of 14 MW with a total distribution network of 12 km in length, to supply an annual energy consumption of 22 million kWh. The expected economic saving is 30% and the reduction in CO2 emissions of 6,800 tonCO2 year . This represents a reduction in CO2 emissions of around 30% of all our energy consumption (electricity and thermal) and an increase of 40% of energy production through renewables in our facilities. The main purpose of this project is to achieve energy efficiency objectives and reduction of carbon emissions, through the implementation of the DH. The DH uses forest biomass as fuel, wood residues from near forests, a renewable energy source with low greenhouse gas emissions, which is cheaper than conventional fossil fuels, and allows us to achieve the objectives of the European strategy 20/20/20. The use of renewable energy is an important commitment by the University of Valladolid to help control its carbon footprint through a significant reduction in CO2 emissions. It also means reducing energy dependence on fossil fuels. The current natural gas boilers will be maintained to ensure the supply of heating in the event of a disconnection of District Heating, thus ensuring thermal comfort in buildings. The biofuel to be used are wood chips, with a grain size range from G50 to G100 and humidity between 20% and 40%. Consumption is 7,886 tons/year, of which the University will consume 6,140 tonCO2 year (77.87%), 183.74 tons per year (2.33%) by buildings of the City Council of Valladolid, and 1,561.43 tonCO2 year (19.80%) by buildings of the Government of Castilla y León. The annual energy generation in 2016–2017 was 22,069 MWh per year to the whole of the network, of which 17,188 MWh correspond to the consumption of the UVa (77.87%), 515,180 kWh correspond to buildings of the City Council of Valladolid (2.33%) and 4,366,685 kWh belong to buildings of the Junta de Castilla y León (19.80%). Some current research focuses on the study and application of the Thermal District concept. Mazhar et al. [1] conducted a review of district heating systems, which highlights the importance of modernizing heating systems to reduce greenhouse gas emissions and contribute to sustainability, for which district heating (DH) systems are considered the future trend. In the view of Gao et al. [2], Figure 1. Plan view of connected buildings. Among the different actions carried out by the University of Valladolid to reduce its carbon footprint, is the implementation of a DH system, making it to also the most important biomass DH in Spain, supplying heating and DHW to 31 buildings. The thermal biomass generation plant has a power of 14 MW with a total distribution network of 12 km in length, to supply an annual energy consumption of 22 million kWh. The expected economic saving is 30% and the reduction in CO 2 emissions of 6800 tonCO2 year . This represents a reduction in CO 2 emissions of around 30% of all our energy consumption (electricity and thermal) and an increase of 40% of energy production through renewables in our facilities. The main purpose of this project is to achieve energy efficiency objectives and reduction of carbon emissions, through the implementation of the DH. The DH uses forest biomass as fuel, wood residues from near forests, a renewable energy source with low greenhouse gas emissions, which is cheaper than conventional fossil fuels, and allows us to achieve the objectives of the European strategy 20/20/20. The use of renewable energy is an important commitment by the University of Valladolid to help control its carbon footprint through a significant reduction in CO 2 emissions. It also means reducing energy dependence on fossil fuels. The current natural gas boilers will be maintained to ensure the supply of heating in the event of a disconnection of District Heating, thus ensuring thermal comfort in buildings. The biofuel to be used are wood chips, with a grain size range from G50 to G100 and humidity between 20% and 40%. Consumption is 7886 tons/year, of which the University will consume 6140 tonCO2 year (77.87%), 183.74 tons per year (2.33%) by buildings of the City Council of Valladolid, and 1561.43 tonCO2 year (19.80%) by buildings of the Government of Castilla y León. The annual energy generation in 2016–2017 was 22,069 MWh per year to the whole of the network, of which 17,188 MWh correspond to the consumption of the UVa (77.87%), 515,180 kWh correspond to buildings of the City Council of Valladolid (2.33%) and 4,366,685 kWh belong to buildings of the Junta de Castilla y León (19.80%). Some current research focuses on the study and application of the Thermal District concept. Mazhar et al. [ 1 ] conducted a review of district heating systems, which highlights the importance of modernizing heating systems to reduce greenhouse gas emissions and contribute to sustainability, for which district heating (DH) systems are considered the future trend. In the view of Gao et al. [ 2 ], District Heating offers greater advantages than traditional heating systems, including energy savings, Energies 2019,12, 93 3 of 20 regulation and control, which represents a great potential for development and a broad market outlook for the future. Regarding DH development, Lund et al. [ 3 ] presented a vision of the future of district heating systems and technologies and their role in a sustainable energy future. Wener [ 4 ] presented the background and current position of district heating and cooling in Sweden. The review structure considers market, technical, supply, environmental, institutional and future contexts. Lygnerud and Werner. [ 5 ] concluded in their study of 107 heat recovery systems in Swedish industry that the recovery of excess industrial heat using DH systems can be characterized by high political interest, high potential, low utilization and often high profitability. A paper submitted by Schmidt [ 6 ] exposes the cooperative work carried out together with the International Energy Agency (IEA) on thermal districts, determined that low-temperature district heating is a key technology that allows to increase the integration of renewable and residual energy for heating and cooling. Turski-Michal and Sekret-Robert [ 7 ] analyzed 63 heating stations in Poland in order to determine the energy effect of using buildings and the district heating network as thermal energy storage to compensate for the reduced heat output of the district heating system. In the investigation some parameters are often underestimated, such as the effect of the incidence of external temperature and the duration of outdoor temperature on the production of heat from DH systems. In Vienna, Fallahnejad et al. [ 8 ] studied the impacts of key DH parameters (ceilings, network costs, DH participation and extension), in which they determined that increasing DH market share in HD areas under a given network cost limit significantly reduces average network costs. Furthermore, that expanding the DH system without increasing market share in the DH areas does not effectively increase the DH share of the total heat demand and leads to higher average network costs. Restoration of existing buildings is considered a way to reduce energy use and carbon emissions in building stock. In [ 9 ] the impact of renovating a single-family home with different district heating systems under restoration packages chosen to represent typical but innovative ways to improve the energy efficiency and indoor climate of a single-family home was studied. In Bavaria, the performance of six small-scale district heating systems (DH) was monitored for 12 months in order to identify typical optimization potentials and develop standards for performance monitoring, analyze extensive operational data and evaluate Key Performance Indicators (KPIs). The KPI showed strong fluctuation and variation between different DH indicating that the main potential was found in the control of DH [10]. One of the challenges for near Zero Energy Consumption Buildings (nZEB) is thermal conditioning systems, Wu et al. [ 11 ] compared the energy performance of HVAC systems for a zero-grid energy residential building (nZEB) in different climate zones in the United States, where air conditioning represented 23.8% to 72.9% of the total energy of the building, depending on the air conditioning option and climate zone. In [ 12 ] whether nZEB standards applied to Italian school buildings guarantee good thermal conditions inside and which building configuration is the most favored was studied. Nielsen and Möller [ 13 ] examined the excess heat production of nZEBs in DH, determined that the excess heat of nZEBs can benefit DH systems by decreasing the production of fuel-using production units. In DH areas where the heat demand in the summer months is already covered by renewable energy, adding seasonal heat storage is essential to achieve nZEB status in terms of efficiency indicators as an alternative to track changes in efficiency and savings obtained according to the improvements implemented. Sekki et al. [ 14 ] showed that energy efficiency can be measured using alternative indicators and confirmed that different indicators have a different impact on the results showing efficiency. In the cases studied, energy savings can be achieved by investing in technical measures or operating the building automation system based on actual occupancy. On the other hand, Abu et al. [ 15 ] reviewed the Energy Efficiency Index (EEI) as an indicator to measure the performance of energy consumption in a building, with the forecast of energy consumption in a building being an important strategy to achieve the goal of reducing energy demand as well as improving energy efficiency. A review of how to quantify the Environmental Building Performance (EBP) is conducted by Maslesa et al. [ 16 ] through a systematic review of the literature where the information obtained is valuable to decision makers and facility managers in the process of implementing an environmental strategy and focusing on improving the EBP. It also concludes that Energies 2019,12, 93 4 of 20 buildings as products are complex structures with a long service life compared to most other products and induce considerable environmental impacts throughout their life cycle. As for the importance of integrating SGEn in the industrial and service sector, ISO 50001 [ 17 ] provides benefits for the industrial and service sectors. It is estimated that the standard, applied to different economic sectors, could influence savings of up to 60% in the world’s energy consumption according to studies by Van der Hoeven [ 18 ]. In terms of methodologies to implement SGEn, Castrillón et al. [ 19 ] reported their results of the implementation of SGIE in a wet cement production facility, which showed an increase in energy efficiency associated with a reduction in electricity consumption of 4.6%, achieved without investment. This means that only through the innovation of processes through applied management technologies, in addition to the adoption of a culture of efficient energy management and continuous improvement. Benedetti et al. [ 20 ] also presented a new methodology for managing energy performance through the development, analysis and maintenance of energy performance indicators in manufacturing plants, taking into account the requirements of ISO 50001:2011 and ISO 50006:2014. The proposed methodology allows for the immediate identification of energy performance deviations from the manufacturing plant through the monitoring of Energy Performance Indicators over time and the identification of possible causes and responsibilities for such deviations. Saadi et al. [ 21 ] described the results of a study that has been carried out to reduce the energy consumption of a library building in the warm Omani climate. The project follows a systematic approach of collecting data from the maintenance department and projects, performing an energy audit and generating a building simulation model. To learn the real scenario of building energy management facilities, Afroz et al. [ 22 ] conducted a study in an institutional building at Murdoch University (Australia), which has incorporated leading-edge technologies over the past two decades. Through this case study analysis, relevant information is revealed that will bring benefits to energy management staff as well as researchers in this area, showing that an efficient energy management system in commercial or institutional buildings can reduce energy consumption and operating costs and provide a comfortable and healthy indoor environment. In terms of management tools and techniques as well as statistics, Castrillón and González [ 23 ] set out the procedures for validating the energy indicators and baselines needed in the energy planning of any institution and its subsequent implementation of an energy planning management system. This represents a reduction of 75% of emissions associated with the thermal supply of heating, allowing us to reduce our carbon footprint. This will also allow us to improve the energy certification of campus buildings. A building that is connected to a network improves its energy rating, an improvement that can be even greater when the network of biomass heaters is used, as in this case. 2. Case Study The Engineering School (EII), is an educational building, mainly dedicated to teaching and research. It has a base area of 20,397 m 2 (16,691.4 m 2 useful), distributed over four floors and the orientation of the main façade is southwest. (Figure 2). Energies 2017, 10, x FOR PEER REVIEW 4 of 20 managers in the process of implementing an environmental strategy and focusing on improving the EBP. It also concludes that buildings as products are complex structures with a long service life compared to most other products and induce considerable environmental impacts throughout their life cycle. As for the importance of integrating SGEn in the industrial and service sector, ISO 50001 [17] provides benefits for the industrial and service sectors. It is estimated that the standard, applied to different economic sectors, could influence savings of up to 60% in the world’s energy consumption according to studies by Van der Hoeven [18]. In terms of methodologies to implement SGEn, Castrillón et al. [19] reported their results of the implementation of SGIE in a wet cement production facility, which showed an increase in energy efficiency associated with a reduction in electricity consumption of 4.6%, achieved without investment. This means that only through the innovation of processes through applied management technologies, in addition to the adoption of a culture of efficient energy management and continuous improvement. Benedetti et al. [20] also presented a new methodology for managing energy performance through the development, analysis and maintenance of energy performance indicators in manufacturing plants, taking into account the requirements of ISO 50001:2011 and ISO 50006:2014. The proposed methodology allows for the immediate identification of energy performance deviations from the manufacturing plant through the monitoring of Energy Performance Indicators over time and the identification of possible causes and responsibilities for such deviations. Saadi et al. [21] described the results of a study that has been carried out to reduce the energy consumption of a library building in the warm Omani climate. The project follows a systematic approach of collecting data from the maintenance department and projects, performing an energy audit and generating a building simulation model. To learn the real scenario of building energy management facilities, Afroz et al. [22] conducted a study in an institutional building at Murdoch University (Australia), which has incorporated leading-edge technologies over the past two decades. Through this case study analysis, relevant information is revealed that will bring benefits to energy management staff as well as researchers in this area, showing that an efficient energy management system in commercial or institutional buildings can reduce energy consumption and operating costs and provide a comfortable and healthy indoor environment. In terms of management tools and techniques as well as statistics, Castrillón and González [23] set out the procedures for validating the energy indicators and baselines needed in the energy planning of any institution and its subsequent implementation of an energy planning management system. This represents a reduction of 75% of emissions associated with the thermal supply of heating, allowing us to reduce our carbon footprint. This will also allow us to improve the energy certification of campus buildings. A building that is connected to a network improves its energy rating, an improvement that can be even greater when the network of biomass heaters is used, as in this case. 2.Case Study The Engineering School (EII), is an educational building, mainly dedicated to teaching and research. It has a base area of 20,397 m2 (16,691.4 m2 useful), distributed over four floors and the orientation of the main façade is southwest. (Figure 2). Figure 2. School of Industrial Engineering (EII) building. Figure 2. School of Industrial Engineering (EII) building. Energies 2019,12, 93 5 of 20 Table 1lists the features and facilities of the EII building. Table 1. Features of the building. System/Fuel Lighting: Heating: HVAC Other Equipment Working Time and Set Point Gas Boilers (NG) LED and fluorescence with electromagnetic ballast T5 Natural Gas boilers of 540 kW “REMEHA” model “Gas-3B/13-Duo” η=0.78 Individual Split-type equipment. 19 Air Handling Units (AHU) Elevators, computers and laptops, printers, laboratory equipment, etc. Weekdays (8 a.m. to 10 p.m.). Saturdays (9 a.m. to 2 p.m.). Set point (21 ◦C) DH (Biomass) LED and fluorescence with electromagnetic ballast T5 119,000 L of backup tank DH substation with flat plate heat exchanger Biomass Boilers 19 MW η=0.85 Individual Split-type equipment 19 Air Handling Units (AHU) Elevators, computers and laptops, printers, laboratory equipment, etc. Weekdays (8 a.m. to 10 p.m.) Saturdays (9 a.m. to 2 p.m.) Set point (21 ◦C) The heat exchange system between the DH and the EII is done through a secondary installation in the facilities of the building, formed by a flat plate exchanger, elements regulating and controlling the heating operation, and a discharge pump, which distributes the fluid to the different areas of the building (Figure 3). Energies 2017, 10, x FOR PEER REVIEW 5 of 20 Table 1 lists the features and facilities of the EII building. Table 1. Features of the building. System/Fuel Lighting: Heating: HVAC Other Equipment Working Time and Set Point Gas Boilers (NG) LED and fluorescence with electromagnetic ballast T5 Natural Gas boilers of 540kW “REMEHA” model “Gas-3B/13-Duo” 𝜂 =0.78 Individual Split-type equipment. 19 Air Handling Units (AHU) Elevators, computers and laptops, printers, laboratory equipment, etc. Weekdays (8 a.m. to 10 p.m.). Saturdays (9 a.m. to 2 p.m.). Set point (21 °C) DH (Biomass) LED and fluorescence with electromagnetic ballast T5 119,000 L of backup tank DH substation with flat plate heat exchanger Biomass Boilers 19 MW 𝜂 =0.85 Individual Split-type equipment 19 Air Handling Units (AHU) Elevators, computers and laptops, printers, laboratory equipment, etc. Weekdays (8 a.m. to 10 p.m.) Saturdays (9 a.m. to 2 p.m.) Set point (21 °C) The heat exchange system between the DH and the EII is done through a secondary installation in the facilities of the building, formed by a flat plate exchanger, elements regulating and controlling the heating operation, and a discharge pump, which distributes the fluid to the different areas of the building (Figure 3). Figure 3. Heat exchanger substation in the facilities at the EII Building. The system consists of two circuits, the primary circuit, from the Central Heat Power (CHP) to the building facilities where the heat exchanger is located, and the secondary circuit, from the heat exchanger to the heat supply to the building. The connection between the DH and the substation of the building is from underground pipes that lead to the installations of the building. Once the supply and return pipes reach the facilities at EII building, the pressure and temperature from the DH is regulated with the heat exchanger to the necessary conditions for the heating demand of the building. Both the supply and return temperatures of the system ought to be regulated in such a way that the temperature that reaches the heat exchanger from the central heating plant is 90 ℃ and the temperature at the outlet of the heat exchanger in the return pipe towards the central heating plant is close to 70 ℃, with a thermal gap of approximately 20 °C. In order to visualize, control and manage at any time the operation of the installation, a dynamic monitoring is developed through sensors and automatons that register data, which are exported and saved in a Supervisory Control And Data Acquisition (SCADA). In this way, any possible failure or warning can be read quickly and the exact location of the problem can be visualized (Figure 4). Figure 3. Heat exchanger substation in the facilities at the EII Building. The system consists of two circuits, the primary circuit, from the Central Heat Power (CHP) to the building facilities where the heat exchanger is located, and the secondary circuit, from the heat exchanger to the heat supply to the building. The connection between the DH and the substation of the building is from underground pipes that lead to the installations of the building. Once the supply and return pipes reach the facilities at EII building, the pressure and temperature from the DH is regulated with the heat exchanger to the necessary conditions for the heating demand of the building. Both the supply and return temperatures of the system ought to be regulated in such a way that the temperature that reaches the heat exchanger from the central heating plant is 90 °C and the temperature at the outlet of the heat exchanger in the return pipe towards the central heating plant is close to 70 °C , with a thermal gap of approximately 20 ◦ C. In order to visualize, control and manage at any time the operation of the installation, a dynamic monitoring is developed through sensors and automatons that register data, which are exported and saved in a Supervisory Control And Data Acquisition (SCADA). In this way, any possible failure or warning can be read quickly and the exact location of the problem can be visualized (Figure 4). Energies 2019,12, 93 6 of 20 Energies 2017, 10, x FOR PEER REVIEW 6 of 20 Figure 4. Screenshot of the EII Heat Exchanger substation SCADA. The data provided by the implemented monitoring system is accessed via SCADA, which manages the DH generation and the heat exchange substations connected to DH. With the SCADA, all the energy parameters being measured by sensors can be visualized, represented and studied in depth. The consumption by the generation and energy demand are displayed in real time, in addition to being able to observe an instant control of the heat exchangers of each building. The SCADA implemented in the system facilitates the analysis and studies of historical consumption and can generate reports on the development of each system together or separately, comparing values between the different buildings connected to the DH. In addition all this data can be exported to Excel files and simulation files, enabling their study through other software programs. In addition, there are different sensors spread throughout the building that record the different energy and comfort parameters in real-time, delivered through a Modbus network to a dynamic monitoring tool, where are stored for subsequent study, using them as an energy management tool for the building. The dynamic monitoring of energy parameters through SCADA, provides information on the use of energy consumed as an energy management tool. In this way, it has been possible to obtain the thermal and electrical energy consumption, where this consumption occurs and the conditions of the spaces where this energy is being supplied. Figure 5 shows the dynamic monitoring of the thermal parameters in the heating system, such as the power generated at the EII building during a period of nine standard winter days. Figure 5. Screenshot of EII SCADA. Thermal power consumed in a winter period. Figure 4. Screenshot of the EII Heat Exchanger substation SCADA. The data provided by the implemented monitoring system is accessed via SCADA, which manages the DH generation and the heat exchange substations connected to DH. With the SCADA, all the energy parameters being measured by sensors can be visualized, represented and studied in depth. The consumption by the generation and energy demand are displayed in real time, in addition to being able to observe an instant control of the heat exchangers of each building. The SCADA implemented in the system facilitates the analysis and studies of historical consumption and can generate reports on the development of each system together or separately, comparing values between the different buildings connected to the DH. In addition all this data can be exported to Excel files and simulation files, enabling their study through other software programs. In addition, there are different sensors spread throughout the building that record the different energy and comfort parameters in real-time, delivered through a Modbus network to a dynamic monitoring tool, where are stored for subsequent study, using them as an energy management tool for the building. The dynamic monitoring of energy parameters through SCADA, provides information on the use of energy consumed as an energy management tool. In this way, it has been possible to obtain the thermal and electrical energy consumption, where this consumption occurs and the conditions of the spaces where this energy is being supplied. Figure 5shows the dynamic monitoring of the thermal parameters in the heating system, such as the power generated at the EII building during a period of nine standard winter days. Energies 2017, 10, x FOR PEER REVIEW 6 of 20 Figure 4. Screenshot of the EII Heat Exchanger substation SCADA. The data provided by the implemented monitoring system is accessed via SCADA, which manages the DH generation and the heat exchange substations connected to DH. With the SCADA, all the energy parameters being measured by sensors can be visualized, represented and studied in depth. The consumption by the generation and energy demand are displayed in real time, in addition to being able to observe an instant control of the heat exchangers of each building. The SCADA implemented in the system facilitates the analysis and studies of historical consumption and can generate reports on the development of each system together or separately, comparing values between the different buildings connected to the DH. In addition all this data can be exported to Excel files and simulation files, enabling their study through other software programs. In addition, there are different sensors spread throughout the building that record the different energy and comfort parameters in real-time, delivered through a Modbus network to a dynamic monitoring tool, where are stored for subsequent study, using them as an energy management tool for the building. The dynamic monitoring of energy parameters through SCADA, provides information on the use of energy consumed as an energy management tool. In this way, it has been possible to obtain the thermal and electrical energy consumption, where this consumption occurs and the conditions of the spaces where this energy is being supplied. Figure 5 shows the dynamic monitoring of the thermal parameters in the heating system, such as the power generated at the EII building during a period of nine standard winter days. Figure 5. Screenshot of EII SCADA. Thermal power consumed in a winter period. Figure 5. Screenshot of EII SCADA. Thermal power consumed in a winter period. Energies 2019,12, 93 7 of 20 2.1. Thermal Consumption The thermal primary energy consumption per year in kWh and the thermal efficiency indicator kWh m2 of Natural Gas in the EII Building at the UVa, between 2006 and 2017, are shown in Table 2, where it can be observed that in 2010 there was a maximum due to the extreme temperatures that winter, and the necessity of higher heating consumption. There is also a reduction in the thermal consumption of Natural Gas from 2015 to 2017, when the DH fueled by biomass begins to work, thus the NG boilers just work as backup. Table 2. Natural Gas Consumption from 2006 to 2017. Year Natural Gas KWh KWh/m2 2006 1,264,550 75.76 2007 1,415,661 84.81 2008 1,632,630 97.81 2009 1,585,339 94.98 2010 2,569,220 153.92 2011 1,639,719 98.24 2012 1,805,694 108.18 2013 1,510,476 90.49 2014 1,229,459 73.66 2015 1,131,687 67.8 2016 87,716 5.26 2017 14,912 0.88 In Table 3, the natural gas thermal consumption of in kWh and the thermal efficiency indicator kWh m2are shown monthly during the last years, from 2014 to 2017. Table 3. Natural Gas consumption and indicators from 2014 to 2017. Month 2014 2015 2016 2017 kWh kWh/m2kWh kWh/m2kWh kWh/m2kWh kWh/m2 January 258,821 15.51 377,448 22.61 33,617 2.01 2834 0.17 February 276,638 16.57 360,341 21.59 2294 0.14 0 0 March 208,852 12.51 243,979 14.62 12,292 0.74 851 0.05 April 116,325 6.97 95,293 5.71 12,395 0.74 5184 0.31 May 18,881 1.13 31,520 1.89 2912 0.17 4587 0.27 June 0 0 0 0 0 0 0 0 July 0 0.02 0 0 0 0 0 0 August 0 0 0 0 0 0 0 0 September 0 0 0 0 574 0.03 0 0 October 993 0.06 5905 0.35 11,345 0.68 672 0.04 November 173,818 10.41 13,498 0.81 8908 0.53 464 0.03 December 174,875 10.48 3703 0.22 3379 0.2 320 0.02 In the variation of the monthly NG thermal efficiency indicator kWh m2 , of the EII in the years 2014 to 2017, it can be noted that the months of June, July, August and in some years the month of September, due to the fact that the external temperature is high and no heating is needed, there is no thermal consumption. However, there is a cost due to the payment of the fee for contracting this NG service. From September 2015, the DH was incorporated, this year still experimentally, which reduced the thermal consumption of NG as well as the cost of heating. As for the cost, the annual consumption of NG for heating shows that the NG tariff was stable, with service contracting rates distributed as follows: 135.13 € kWh in 2014, 137.42 € kWh in 2015, 136.54 € kWh Energies 2019,12, 93 8 of 20 in 2016 and 108.66 € kWh in 2017. Another important consideration is that in September 2015, both the energy and cost decreased thanks to the supply of the DH network. The biomass District Heating (DH) system began to provide heating service to the entire campus of the UVa and therefore the EII building in 2015/2016, with the year 2016/2017 being the second year of operation of the DH network. Table 4shows the variation in thermal consumption of DH fed with biomass in kWh and the thermal efficiency indicator of biomass in the EII during the different months of the year. Table 4. Primary energy consumption and energy efficiency indicator for DH fueled by biomass from 2015 to 2017. Month 2015 2016 2017 KWh KWh/m2KWh KWh/m2KWh KWh/m2 January 0 0 167,500 10.04 283,300 16.97 February 0 0 224,400 13.44 205,600 12.32 March 0 0 164,700 9.87 169,800 10.17 April 0 0 126,700 7.59 53,800 3.22 May 0 0 18,900 1.13 13,800 0.83 June 0 0 0 0 0 0 July 0 0 0 0 0 0 August 0 0 0 0 0 0 September 0 0 0 0 0 0 October 53,000 3.18 28,700 1.72 8100 0.49 November 173,200 10.38 186,100 11.15 137,100 8.21 December 126,100 7.55 177,400 10.63 219,700 13.16 The evolution of thermal consumption of biomass, as well as the thermal efficiency indicator, shows the impact that the incorporation of viomass, through the DH system, has had for the EII as a substitute fuel for NG in recent years. Figure 6shows the evolution of the thermal efficiency indicator for biomass as a function of the months in the annual period 2015–2017. Energies 2017, 10, x FOR PEER REVIEW 8 of 20 136.54 € kWh in 2016 and 108.66 € kWh in 2017. Another important consideration is that in September 2015, both the energy and cost decreased thanks to the supply of the DH network. The biomass District Heating (DH) system began to provide heating service to the entire campus of the UVa and therefore the EII building in 2015/2016, with the year 2016/2017 being the second year of operation of the DH network. Table 4 shows the variation in thermal consumption of DH fed with biomass in kWh and the thermal efficiency indicator of biomass in the EII during the different months of the year. Table 4. Primary energy consumption and energy efficiency indicator for DH fueled by biomass from 2015 to 2017. Month 2015 2016 2017 KWh KWh/m² KWh KWh/m² KWh KWh/m² January 0 0 167,500 10.04 283,300 16.97 February 0 0 224,400 13.44 205,600 12.32 March 0 0 164,700 9.87 169,800 10.17 April 0 0 126,700 7.59 53,800 3.22 May 0 0 18,900 1.13 13,800 0.83 June 0 0 0 0 0 0 July 0 0 0 0 0 0 August 0 0 0 0 0 0 September 0 0 0 0 0 0 October 53,000 3.18 28,700 1.72 8100 0.49 November 173,200 10.38 186,100 11.15 137,100 8.21 December 126,100 7.55 177,400 10.63 219,700 13.16 The evolution of thermal consumption of biomass, as well as the thermal efficiency indicator, shows the impact that the incorporation of viomass, through the DH system, has had for the EII as a substitute fuel for NG in recent years. Figure 6 shows the evolution of the thermal efficiency indicator for biomass as a function of the months in the annual period 2015–2017. Figure 6. Evolution of the thermal efficiency indicator for biomass use 2015–2017. With regards to the composition of the EII’s energy matrix, in 2014 the NG heating system was installed, having only high thermal consumption in the months when the heating started. At that time the distribution of annual energy consumption corresponded to 56% thermal consumption of NG and the remaining 44% is a consumption of electricity. In 2015, the implementation of the DH began, which began to heat the EII in September, however due to climate issues, heating consumption in this month was not necessary. Once the DH began to operate, there was a decrease in NG consumption, resulting in an annual distribution of thermal consumption for that year of 14% biomass, and 44% NG. Due to the fact that the DH was 0 2 4 6 8 10 12 14 16 18 kWh/m2 2015 2016 2017 Figure 6. Evolution of the thermal efficiency indicator for biomass use 2015–2017. With regards to the composition of the EII’s energy matrix, in 2014 the NG heating system was installed, having only high thermal consumption in the months when the heating started. At that time the distribution of annual energy consumption corresponded to 56% thermal consumption of NG and the remaining 44% is a consumption of electricity. In 2015, the implementation of the DH began, which began to heat the EII in September, however due to climate issues, heating consumption in this month was not necessary. Once the DH began to operate, there was a decrease in NG consumption, resulting in an annual distribution of thermal consumption for that year of 14% biomass, and 44% NG. Due to the fact that the DH was operating in experimental mode with frequent stoppages for verification and balancing, the EII was still using Energies 2019,12, 93 9 of 20 the old NG system to allow for proper heating. The remaining 42% corresponds to the electricity consumption of the EII. In 2016, the DH fueled by biomass started to operate year round, with the support of the NG boiler system of the EII. The distribution of the energy matrix throughout this year corresponded to a thermal consumption of 51% of biomass consumption, and 4% of NG, and the remainder is electricity consumption. In this year and subsequent ones the percentage of NG consumption is very low compared to the biomass consumption. In 2017, the second year of operation of the DH network, the installation was already at an optimum point of operation due to its balance. In this way the distribution of energy consumption percentage of the EII, corresponds to a thermal consumption of 72% of biomass consumption and 1% of NG. As a result, NG consumption is negligible compared to biomass consumption. The annual analysis is shown in Table 5, where it is possible to decrease both the kWh/m2indicator and the economic expenditure €/m2. Table 5. Indicator of thermal efficiency kWh m2 and energy cost Indicator of biomass from 2015 to 2017. Biomass €€/m2KWh KWh/m2 2015 28,359.93 1.7 352,300 21.11 2016 88,060.89 5.28 1,094,400 65.57 2017 87,803.40 5.26 1,091,200 65.43 Figure 7shows the thermal consumption of the EII building from 2014 to 2017. Energies 2017, 10, x FOR PEER REVIEW 9 of 20 operating in experimental mode with frequent stoppages for verification and balancing, the EII was still using the old NG system to allow for proper heating. The remaining 42% corresponds to the electricity consumption of the EII. In 2016, the DH fueled by biomass started to operate year round, with the support of the NG boiler system of the EII. The distribution of the energy matrix throughout this year corresponded to a thermal consumption of 51% of biomass consumption, and 4% of NG, and the remainder is electricity consumption. In this year and subsequent ones the percentage of NG consumption is very low compared to the biomass consumption. In 2017, the second year of operation of the DH network, the installation was already at an optimum point of operation due to its balance. In this way the distribution of energy consumption percentage of the EII, corresponds to a thermal consumption of 72% of biomass consumption and 1% of NG. As a result, NG consumption is negligible compared to biomass consumption. The annual analysis is shown in Table 5, where it is possible to decrease both the kWh/m2 indicator and the economic expenditure €/m2. Table 5. Indicator of thermal efficiency ( kWh m2 ) and energy cost Indicator of biomass from 2015 to 2017. Biomass € €/m² KWh KWh/m² 2015 28,359.93 1.7 352,300 21.11 2016 88,060.89 5.28 1,094,400 65.57 2017 87,803.40 5.26 1,091,200 65.43 Figure 7 shows the thermal consumption of the EII building from 2014 to 2017. (a) (b) (c) (d) Figure 7. Thermal consumption of the EII building: (a) 2014; (b) 2015; (c) 2016; (d) 2017. Figure 8 shows how from 2015 onwards, the energy cost of natural gas decreases and is replaced by the energy cost of biomass itself; each year the cost is lower as the operation of DH stabilizes. 0 50 100 150 200 250 300 MWh 0 50 100 150 200 250 300 350 400 MWh 0 50 100 150 200 250 MWh 0 50 100 150 200 250 300 MWh Natural Gas Biomass b) c) d) Figure 7. Thermal consumption of the EII building: (a) 2014; (b) 2015; (c) 2016; (d) 2017. Figure 8shows how from 2015 onwards, the energy cost of natural gas decreases and is replaced by the energy cost of biomass itself; each year the cost is lower as the operation of DH stabilizes. Energies 2019,12, 93 16 of 20 If there is no change in consumption patterns, the CUSUM curve will simply fluctuate around a horizontal line. However, any change that produces savings or overconsumption will cause the curve to change direction: a downward trend in the curve means savings, while upward trend means losses. 2.1.3. The Primary Energy Indicator All supplied and exported energies are shown in a unique indicator. The Primary Energy Indicator [ 24 ]. It is estimated from the energy supplied and exported using the national primary energy factor provided by Equation (20): Ep,nrem =∑ i (Edel,i fdel,nren,i)−∑ iEexp,i fdel,exp,i(20) EPp=Ep,nrem Anet (21) The near zero performance is determined by a Member State in terms of the achievable use of primary energy, the amount of primary energy supplied by renewable energies, the funding sources available for renewable energies or measures of energy efficiency, and the funding requirements as well as the successful rate of the definition. The nZEBs such as a national use of energy which is both technologically and economically feasible of > 0 kWh m2·y , but no higher than the existing Member State threshold for primary non-renewable energy, generated by a mix of best-practice eco-efficiency measures plus renewable energy systems that might or might not be cost-optimal. Spain has not yet laid down this limited value; however a draft document already exists for the Spanish Building Standard, in which bases for the nZEB definition are argued. December 2016, Spain launched a report to improve the evaluation methodology in high-efficiency buildings. This report suggests several predefined thresholds based on those defined in July 2016, in Commission Recommendation (EU) 2016/1318. In Valladolid, which is marked by a continental climate, such upper limits are 85 to 100 kWh m2·y , of primary energy consumption in total on these buildings per year. At this point, the upper limits for non-renewable sources are between 45 and 55 kWh m2·y , and the contribution of renewable energies to primary energy must be greater than 45 kWh m2·y. In order to assess compliance with upper limits for primary energy, primary energy intensity at the buildings is established. The methodology for the determination applies the approach defined in the general EPBD regulation. Primary energy is obtained as the biomass heat being multiplied by its respective Primary Energy Factor (PEF). Moreover, in Spain, these factors are established in its own normative, the Spanish Thermal Systems Standard. [ 25 ]. The Table 12 shows the primary energy factors by each energy sources. Table 12. Spanish Building Standards. PEF established. Energy Source Use fp,ren fp,nren Electricity Grid Input 0.414 1.954 Biomass On-site Input 1.003 0.034 Natural Gas On-site Input 0.005 1.190 Non-renewable conversion factor fi applied to the electricity, in Spain, is currently 1.954. Biomass conversion factor is 1. At the EII Building there is still no export energy, so the primary energy indicator is: EPp=118.15 kWh m2·y(22) Energies 2019,12, 93 17 of 20 This is caused by the main factors established by Spain, which has a higher percentage of renewable energy (17%) from the main electricity grid. by calculating the on-site RES share; it is possible to obtain a high value of 91.21 kWh m2·y of contribution of renewable primary energies in situ at EII Building. A primary energy comparison between the EII Building versus the assumptions of the European Union Directive is shown in Figure 13. Energies 2017, 10, x FOR PEER REVIEW 17 of 20 This is caused by the main factors established by Spain, which has a higher percentage of renewable energy (17%) from the main electricity grid. by calculating the on-site RES share; it is possible to obtain a high value of 91.21 kWh m2·y of contribution of renewable primary energies in situ at EII Building. A primary energy comparison between the EII Building versus the assumptions of the European Union Directive is shown in Figure 13. Figure. 13. Primary energy comparison between EII Building and EU Directive. The detailed boundary values of the system for determining the energy supplied and exported are added based on the reference limits of EN15603. As started in the EPBD recast, the positive influence of renewable energy produced on site is taken into account, so that it reduces the amount of energy supplied needed and can be exported when the demand of buildings is exceeded [26]. 2.1.4. Renewable Energy Ratio To calculate the amount of renewable energies used through the Renewable Energy Ratio (RER), all sources of renewable energy must be taken into account [24]. Solar thermal energy, photovoltaic energy, wind and hydro power, renewable energy from environmental heat sources such as heat pumps and free cooling, renewable biofuels and off-site renewable energy are included. The environmental heat sources used in heat pumps and free-cooling should be considered within the limit of the use of renewable energy system, because in the calculation of the RER, free cooling and heat pumps are not only accounted for in the calculation of the energy supplied on the basis of the COP, it is also taken into account in the energy obtained from the environmental heat sources. The proportion of renewable energies are obtained in relation to the energy use at the building, as total primary energy. It keeps in mind that the exported energy compensates for the energy delivered. By default, the export energies are deemed to compensate for the grid mix or for thermal energy, the grid mix of district cooling or heating. For near and on-site renewable energies, the overall primary energy factor is 1.0, while the primary non-renewable energy factor is 0. The total primary energy-based RER is given by the following Equation (23): 𝑅𝐸𝑅p=∑𝐸ren,ii +∑((𝑓del,tot,i −𝑓del,nrem,i)𝐸del,i) i ∑𝐸ren,ii +∑(𝐸del,i𝑓del,tot,i) i−∑(𝐸exp,i𝑓exp,tot,i) i (23) In order to get the 𝑅𝐸𝑅p for the EII building, 𝑓i= 1.954 is taken for electricity and 𝑓i= 1.003 for the biomass and 𝑓i= 0.005 for natural gas: 𝑅𝐸𝑅p= 0.424 (24) Thus, the RER, the third indicator, in the EII Building is 0.424, a figure that represents a pretty strong result for a nZEB. 3. Conclusions 0 50 100 150 EU EII building Energy non-renewable Energy renewable Figure 13. Primary energy comparison between EII Building and EU Directive. The detailed boundary values of the system for determining the energy supplied and exported are added based on the reference limits of EN15603. As started in the EPBD recast, the positive influence of renewable energy produced on site is taken into account, so that it reduces the amount of energy supplied needed and can be exported when the demand of buildings is exceeded [26]. 2.1.4. Renewable Energy Ratio To calculate the amount of renewable energies used through the Renewable Energy Ratio (RER), all sources of renewable energy must be taken into account [ 24 ]. Solar thermal energy, photovoltaic energy, wind and hydro power, renewable energy from environmental heat sources such as heat pumps and free cooling, renewable biofuels and off-site renewable energy are included. The environmental heat sources used in heat pumps and free-cooling should be considered within the limit of the use of renewable energy system, because in the calculation of the RER, free cooling and heat pumps are not only accounted for in the calculation of the energy supplied on the basis of the COP, it is also taken into account in the energy obtained from the environmental heat sources. The proportion of renewable energies are obtained in relation to the energy use at the building, as total primary energy. It keeps in mind that the exported energy compensates for the energy delivered. By default, the export energies are deemed to compensate for the grid mix or for thermal energy, the grid mix of district cooling or heating. For near and on-site renewable energies, the overall primary energy factor is 1.0, while the primary non-renewable energy factor is 0. The total primary energy-based RER is given by the following Equation (23): RERp=∑iEren,i +∑i(( fdel,tot,i −fdel,nrem,i)Edel,i) ∑iEren,i +∑i(Edel,i fdel,tot,i)−∑i(Eexp,i fexp,tot,i)(23) In order to get the RERp for the EII building, fi= 1.954 is taken for electricity and fi= 1.003 for the biomass and fi=0.005 for natural gas: RERp=0.424 (24) Thus, the RER, the third indicator, in the EII Building is 0.424, a figure that represents a pretty strong result for a nZEB. Energies 2019,12, 93 18 of 20 3. Conclusions A new analysis methodology is developed, based on the ISO 50001 standard, for the implementation of new thermal systems in building renovations using DH. Following the objectives proposed by the EU 20-20-20, in this research paper several indicators have been studied to analyze how the change from a NG heating system to a new heating system by a DH fueled by biomass on an educational building, Engineering School (EII) at the University of Valladolid in Spain. The proposed indicators as well as the methodology carried out fulfill the established methodology in an SGEn ISO 50001, in order to analyze the reduction of thermal energy and the heating efficiency of the building with the new DH as well as the economic costs and environmental impact. Finally, three new indicators of non-renewable primary energy, renewable energy and RER are determined in accordance with the methodology proposed by EU EN15603 and EPBD 2018 on nZEB buildings. With them it is possible to determine the level of impact in the change towards the DH, without altering the energy demand of the building, looking for this existing building to look like a building nZEB. The results obtained show that the new DH system fueled by biomass in the period 2016–2017 achieved a percentage of thermal savings of 19.9 %, an economic saving of 59.2%, a carbon emissions saving of 93.2 %, and an increase in renewable energy of 77.5 %. The I 100 indicator has also been calculated and it can be seen that the DH system works in an energy efficiency zone, as well as the Cusum indicator where its graph shows that the DH over the entire time interval has a tendency to improve energy efficiency. The I 100 indicator has also been calculated and it can be seen that the DH system works in an energy efficiency zone, as well as the Cusum Indicator where its graph shows that the DH over the entire time interval has a tendency to improve energy efficiency. The non-renewable and renewable primary energy indicators of the new DH system fueled by biomass have been calculated at 118.15 kWh m2·y , 91.21 kWh m2·y , respectively. They have been compared by the values that the EU proposes for a city like Valladolid and it is observed that the non-renewable primary energy is of the order of the double, this is due that although the DH reduces in a 40% with respect to the old system of NG, to the lack of any improvement in the energetic demand of the building EII does that this value is still high. On the other hand, however, the comparison of primary renewable energy from DH reaches a value of double, so that the integration of renewable energy by the change from NG to biomass gives a high RER of 0.4. Use of these indicators proposed by the EU enables us to conclude that the change in the EII building of the new DH system fueled by biomass achieves an approximation towards the nZEB building status. Author Contributions: Conceptualization F.J.R.M., J.S.J.A. and E.V.G.; Methodology, F.J.R.M., R.C.M. and J.M.R.H.; Validation, F.J.R.M., J.S.J.A. and E.V.G.; Formal Analysis, R.C.M. and J.M.R.H.; Investigation, R.C.M. and J.M.R.H.; Resources, F.J.R.M., J.S.J.A. and E.V.G.; Data Curation, R.C.M.; Writing—Original Draft Preparation, R.C.M.; Writing—Review & Editing, J.M.R.H.; Visualization, F.J.R.M., J.S.J.A. and E.V.G.; Supervision, F.J.R.M. Funding: This research received no external funding. Acknowledgments: The support of La Dirección de Investigaciones y Desarrollo Tecnológico (DIDT) de la Universidad Autónoma de Occidente (UAO) in Cali, Colombia, and the RETO GIRTER Project (New Intelligent Manager for Thermal Networks), project funded by the European Regional Development Fund through the Ministry of Economy, Industry and Competitiveness of the Government of Spain, have made this research paper feasible. Conflicts of Interest: The authors declare no conflict of interest. Energies 2019,12, 93 19 of 20 Abbreviations DH: District Heating. NG: Natural Gas. DHW: Domestic Hot Water. EII: Engineering School of Valladolid. UVa: University of Valladolid. EMS: Energy Management System. nZEB: near Zero Energy Building. IEA: International Energy Agency. KPI: Key Performance Indicators. HVAC: Heat Ventilation Air Conditioning. EEI: Energy Efficiency Index. EBP: Environmental Building Performance. LED: Light Emitting Diode. ISO: International Organization for Standardization. SCADA: Supervisory Control And Data Acquisition. EPp:Primary energy indicator (kWh m2·y). Ep,nrem:Primary Non-renewable energy (kWh y). Edel,i:Energy supplied on-site or around the site (kWh y). fdel,exp,i: Primary Non-renewable energy of energy supplied compensated by export-energy, which is by default the same as the factor of the energy supplied, if not defined at national level otherwise. fdel,nren,i:Primary Non-renewable energy for the energy supplied. Eexp,i:Export-Energy on site or around the site (kWh y). Anet:Usable space (m2). RERp:Renewable Energy Ratio from primary energy. Eren,i:Renewable energy generated on-site or around the site, (kWh y). fdel,tot,i:Total primary energy factor for the delivered energy carrier i. fexp.tot,i:Total primary energy factor of the delivered energy compensated by the exported energy for energy carrier i. References 1. Mazhar, A.R.; Liu, S.; Shukla, A. A state of art review on the district heating systems. Renew. Sustain. Energy Rev. 2018,96, 420–439. [CrossRef] 2. Gao, L.; Cui, X.; Ni, J.; Lei, W.; Huang, T.; Bai, C.; Yang, J. Technologies in Smart District Heating System. Energy Procedia 2017,142, 1829–1834. [CrossRef] 3. Lund, H.; Duic, N.; Østergaard, P.A.; Mathiesen, B.V. Future district heating systems and technologies: On the role of smart energy systems and 4th generation district heating. Energy 2018,165, 614–619. [CrossRef] 4. Werner, S. District heating and cooling in Sweden. Energy 2017,126, 419–429. [CrossRef] 5. Lygnerud, K.; Werner, S. Risk assessment of industrial excess heat recovery in district heating systems. Energy 2018,151, 430–441. [CrossRef] 6. Schmidt, D. Low Temperature District Heating for Future Energy Systems. Energy Procedia 2018 ,149, 595–604. [CrossRef] 7. Turski, M.; Sekret, R. Buildings and a district heating network as thermal energy storages in the district heating system. Energy Build. 2018,179, 49–56. [CrossRef] 8. Fallahnejad, M.; Hartner, M.; Kranzl, L.; Fritz, S. Impact of distribution and transmission investment costs of district heating systems on district heating potential. Energy Procedia 2018,149, 141–150. [CrossRef] 9. Lidberg, T.; Gustafsson, M.; Myhren, J.A.; Olofsson, T.; Ödlund (former Trygg), L. Environmental impact of energy refurbishment of buildings within different district heating systems. Appl. Energy 2018 ,227, 231–238. [CrossRef] 10. Bücker, D.; Jell, P.; Botsch, R. Performance monitoring of rural district heating systems. Energy Procedia 2018 , 149, 5–14. [CrossRef] Energies 2019,12, 93 20 of 20 11. Wu, W.; Skye, H.M. Net-zero nation: HVAC and PV systems for residential net-zero energy buildings across the United States. Energy Convers. Manag. 2018,177, 605–628. [CrossRef] 12. Zinzi, M.; Pagliaro, F.; Agnoli, S.; Bisegna, F.; Iatauro, D. Assessing the overheating risks in Italian existing school buildings renovated with nZEB targets. Energy Procedia 2017,142, 2517–2524. [CrossRef] 13. Nielsen, S.; Möller, B. Excess heat production of future net zero energy buildings within district heating areas in Denmark. Energy 2012,48, 23–31. [CrossRef] 14. Sekki, T.; Airaksinen, M.; Saari, A. Effect of energy measures on the values of energy efficiency indicators in Finnish daycare and school buildings. Energy Build. 2017,139, 124–132. [CrossRef] 15. Abu Bakar, N.N.; Hassan, M.Y.; Abdullah, H.; Rahman, H.A.; Abdullah, M.P.; Hussin, F.; Bandi, M. Energy efficiency index as an indicator for measuring building energy performance: A review. Renew. Sustain. Energy Rev. 2015,44, 1–11. [CrossRef] 16. Maslesa, E.; Jensen, P.A.; Birkved, M. Indicators for quantifying environmental building performance: A systematic literature review. J. Build. Eng. 2018,19, 552–560. [CrossRef] 17. International Organization for Standardization (ISO). Available online: https://www.iso.org/ (accessed on 14 September 2018). 18. Van der Hoeven, M. Energy Efficiency Indicators: Essentials for Policy Making. Available online: https://webstore. iea.org/energy-efficiency-indicators-essentials-for-policy-making (accessed on 5 November 2018). 19. González, A.J.; Castrillón, R.; Quispe, E.C. Energy Efficiency Improvement in the Cement Industry through Energy Management. In Proceedings of the 2012 IEEE-IAS/PCA 54th Cement Industry Technical Conference, San Antonio, TX, USA, 14–17 May 2012. [CrossRef] 20. Benedetti, M.; Cesarotti, V.; Introna, V. From energy targets setting to energy-aware operations control and back: An advanced methodology for energy efficient manufacturing. J. Clean. Prod. 2017 ,167, 1518–1533. [CrossRef] 21. Al-Saadi, S.N.J.; Ramaswamy, M.; Al-Rashdi, H.; Al-Mamari, M.; Al-Abri, M. Energy Management Strategies for a Governmental Building in Oman. Energy Procedia 2017,141, 206–210. [CrossRef] 22. Afroz, Z.; Higgins, G.; Urmee, T.; Shafiullah, G. Technological Advancement of Energy Management Facility of Institutional Buildings: A Case Study. Energy Procedia 2017,142, 3088–3095. [CrossRef] 23. Castrillon, R.; González, A. Metodología Para la Planificación Energética a Partir de la Norma ISO 50001; Editorial Universidad Autónoma de Occidente: Cali, Colombia, 2018; ISBN 978-958-8994-59-8. 24. Rey-Hernández, J.M.; Velasco-Gómez, E.; San José-Alonso, J.F.; Tejero-González, A.; Rey-Martínez, F.J. Energy analysis at a near zero energy building. A case-study in Spain. Energies 2018,11, 857. [CrossRef] 25. CTE (Spanish Technical Building Code). Available online: http://www.codigotecnico.org (accessed on 15 September 2018). 26. European Directives. Available online: http://eur-lex.europa.eu/ (accessed on 14 September 2018). © 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).