Synthesis and Optimal Operation of Smart Microgrids Serving a Cluster of Buildings on a Campus with Centralized and Distributed Hybrid Renewable Energy Units

Micro-district heating networks based on cogeneration plants and renewable energy technologies are considered efficient, viable and environmentally-friendly solutions to realizing smart multi-energy microgrids. Nonetheless, the energy production from renewable sources is intermittent and stochastic, and cogeneration units are characterized by fixed power-to-heat ratios, which are incompatible with fluctuating thermal and electric demands. These drawbacks can be partially overcome by smart operational controls that are capable of maximizing the energy system performance. Moreover, electrically driven heat pumps may add flexibility to the system, by shifting thermal loads into electric loads. In this paper, a novel configuration for smart multi-energy microgrids, which combines centralized and distributed energy units is proposed. A centralized cogeneration system, consisting of an internal combustion engine is connected to a micro-district heating network. Distributed electric heat pumps assist the thermal production at the building level, giving operational flexibility to the system and supporting the integration of renewable energy technologies, i.e., wind turbines, photovoltaic panels, and solar thermal collectors. The proposed configuration was tested in a hypothetical case study, namely, a University Campus located in Trieste, Italy. The system operation is based on a cost-optimal control strategy and the effect of the size of the cogeneration unit and heat pumps was investigated. A comparison with a conventional configuration, without distributed heat pumps, was also performed. The results show that the proposed configuration outperformed the conventional one, leading to a total-cost saving of around 8%, a carbon emission reduction of 11%, and a primary energy saving of 8%.


Introduction
Cogeneration of useful heat and electrical power for an urban district or a cluster of buildings is a technically mature, environmentally-friendly and cost-effective solution, supported by the European Union Directive [1] on energy efficiency, together with the use of renewable energy sources (RES).Indeed, the Directive [2] on energy performance of buildings indicates four high-efficiency technologies, whose feasibility should be evaluated prior to construction of any new building: (a) decentralized energy supply systems based on RES; (b) cogeneration; (c) district or block heating or cooling; and (d) heat pumps.
In view of this, several different configurations of distributed energy systems (DES) have been investigated in recent years, mostly focusing on cogeneration units and RES technologies.Pagliarini and Ranieri [3] studied the effectiveness of thermal storage coupled with a cogeneration engine to satisfy the energy requirements of a university campus, and stressed the importance of the sizing of the storage.Bracco et al. [4] dealt with the topic of distributed generation, by presenting the University of Genoa polygeneration microgrid, which is based on RES and cogeneration units.A tool for the optimal integrated design and operation of a trigeneration system serving a cluster of buildings was proposed by Piacentino et al. [5,6].The optimal design and operation of a hybrid renewable energy system based on an internal combustion engine and photovoltaic panels was investigated by Destro et al. [7].
In addition to the synthesis and design problems, the optimal operational strategy of energy microgrids has also received considerable attention.Indeed, the adoption of smart control techniques can significantly improve the economic and environmental performances of those systems [8].For example, Roldán-Blay et al. [9] developed an algorithm for the optimal management of a RES-based electric microgrid.Similarly, Phan et al. [10] investigated schedule strategies to minimize the operating cost of a building energy system with photovoltaic panels and a wind micro-turbine.Asaleye et al. [11] proposed a decision-support tool that identifies the optimal operation of renewable energy microgrids by considering forecast of climate variables.
All the above-mentioned works show the energy, environmental and economic effectiveness of cogeneration systems and renewable energy technologies in smart energy grids.Nevertheless, some unresolved issues remain.Indeed, the intermittent and stochastic nature of RES limits their use, and cogeneration units are characterized by a fixed power-to-heat ratio, thus, they fail to match both fluctuating thermal and electric demands.In this context, electrically driven heat pumps may represent an interesting solution due to their ability to shift thermal loads into electric loads.Moreover, heat pumps are a mature and efficient technology, and they are especially suited to the implementation of smart control strategies [12].
For those reasons, the present work discusses a novel configuration for smart multi-energy microgrids, which consists of distributed energy units and a centralized cogeneration unit feeding a micro-district heating network.Specifically, we investigate the benefits of integrating reversible heat pumps for heating and cooling purposes at the building level.The heat pumps represent an interconnection between the electricity and heating networks, therefore, they can be used to increase the operational flexibility of the microgrid and support the integration of renewable energy technologies, i.e., wind turbine, photovoltaic panels, and solar thermal collectors.
The paper is structured as follows.Section 2 presents the design and modeling methodology of the smart multi-energy microgrid.Sections 3 and 4 present the case study and the optimization problem and methodology, which are used in Section 5 to compare the proposed configuration using distributed heat pumps, with a more conventional solution that employs a centralized CHP (Combined Heat and Power) system and natural-gas boilers.Finally, Section 6 presents the concluding remarks.

Energy System Overview and Modeling
In this work, we refer to multi-energy microgrids of medium dimensions with different buildings and loads, using an integrated thermal and electrical energy production system fed by traditional and renewable sources to concurrently satisfy various services (heating, cooling, electrical energy and domestic hot water).The considered generators are: (i) a CHP consisting in an internal combustion engine (ICE); (ii) natural gas boilers; (iii) heat pumps and chillers; (iv) solar thermal collectors; (v) wind turbines; and (vi) photovoltaic modules.Thermal storage is also considered. Figure 1 shows a simplified classification scheme of the reference energy system.
As is well-known, the traditional design approach based on a separate analysis of each component represents a suboptimal design method for multi-energy systems [13].The so-called simulation-based optimization methods are the most recognized procedures to investigate the best synthesis, sizing and control of integrated system through the simulation of the operative performances.Therefore, in the following sub-section we present the operative dynamic model of each block listed in Figure 1.As is well-known, the traditional design approach based on a separate analysis of each component represents a suboptimal design method for multi-energy systems [13].The so-called simulation-based optimization methods are the most recognized procedures to investigate the best synthesis, sizing and control of integrated system through the simulation of the operative performances.Therefore, in the following sub-section we present the operative dynamic model of each block listed in Figure 1.

Modeling of the System Components
The components models must be based on a proper trade-off between the accuracy of the results and computational effort.The latter feature is essential to allow their employment within the optimization procedure to identify the most efficient design and operation of the smart multi-energy microgrid.
The ICE is modeled through performance curves taken from [14], which provide the thermal and electric efficiency, and , respectively, as a function of the engine load factor .The electric ( ) and thermal ( ) output are evaluated based on the load factor, , and the ICE nominal electric power capacity ( ), as follows: The boilers can be modeled with a constant efficiency ( ) over their whole operating range.Photovoltaic panels are simulated through the model provided by [15], which considers the PV performance as a function of the solar irradiance, PV characteristics and cell array temperature:

Modeling of the System Components
The components models must be based on a proper trade-off between the accuracy of the results and computational effort.The latter feature is essential to allow their employment within the optimization procedure to identify the most efficient design and operation of the smart multi-energy microgrid.
The ICE is modeled through performance curves taken from [14], which provide the thermal and electric efficiency, η el and η th , respectively, as a function of the engine load factor L ICE .The electric (E ICE ) and thermal (Q ICE ) output are evaluated based on the load factor, L ICE , and the ICE nominal electric power capacity (E nom ICE ), as follows: The boilers can be modeled with a constant efficiency (η B ) over their whole operating range.Photovoltaic panels are simulated through the model provided by [15], which considers the PV performance as a function of the solar irradiance, PV characteristics and cell array temperature: The thermal performances of the ST are evaluated through the classical model illustrated in [16], based on the characteristics of the panel in terms of transmittance and absorptance factors for normal irradiance (< τα > n ), removal factor (F r ) and frontal losses (U l ).The equations read: Energies 2019, 12, 745 4 of 17 The electrical power generated by the wind turbine varies as the cube of the wind speed, between a cut-in speed and a nominal speed; the latter corresponds to the nominal electrical power generated by the wind turbine.The nominal electrical power is generated between the nominal wind speed and a cut-out speed.The equations, in accordance with [17], read: where k is a coefficient that depends on the characteristic curve of the generator.
The internal energy variation in the thermal storage is calculated considering the thermal fluxes provided to the water volume from all the connected generators and the heat delivered to the load: where the heat losses of the storage tank are evaluated as Reversible heat pumps and chiller performance are evaluated by means of the so-called second-law efficiency [18].The method reads: (11) where COP id and EER id are the coefficients of performance of a reversed Carnot cycle operating between the source and sink temperatures.According to manufacturers, both η I I H and η I I C can be assumed as constant.
The generators are connected to the thermal storages and the buildings through a district heating network (DHN), whose heat losses can be modeled as follows: Finally, the heating/cooling loads of the buildings can be evaluated through a model that correlates the sol-air temperature [19] with the energy load of the building, based on the standard EN 15306 [20].This model is further improved by considering the effect of the building thermo-physical properties in shifting the influence of the external climate on the heating/cooling load.
Energies 2019, 12, 745 5 of 17 Further details on Equations ( 13)-( 16) can be found in [21].This model represents a good trade-off between simplified models (e.g., the energy signature method [22]), which simply correlates external temperature and heating/cooling load, and dynamic models (e.g., TRNSYS or EnergyPlus), which include the building inertia characteristics, solar radiation, and internal loads, providing more accurate results, but requiring a detailed knowledge of the building envelope and heat gain profiles.

Case Study
In this work, we refer to an integrated energy system serving a hypothetical campus, located in Trieste, Italy.This city has a favorable climate, where RES (solar thermal, photovoltaic modules, wind turbines) can provide a significant amount of energy.The Italian Thermotechnical Committee (CTI) provides hourly profiles of external temperature, global solar irradiance on the horizontal plane and wind speed [23].The monthly-average values of the external temperature and irradiance on horizontal plane are reported in Figure 2.
Further details on Equations (13-16) can be found in [21].This model represents a good tradeoff between simplified models (e.g., the energy signature method [22]), which simply correlates external temperature and heating/cooling load, and dynamic models (e.g., TRNSYS or EnergyPlus), which include the building inertia characteristics, solar radiation, and internal loads, providing more accurate results, but requiring a detailed knowledge of the building envelope and heat gain profiles.

Case Study
In this work, we refer to an integrated energy system serving a hypothetical campus, located in Trieste, Italy.This city has a favorable climate, where RES (solar thermal, photovoltaic modules, wind turbines) can provide a significant amount of energy.The Italian Thermotechnical Committee (CTI) provides hourly profiles of external temperature, global solar irradiance on the horizontal plane and wind speed [23].The monthly-average values of the external temperature and irradiance on horizontal plane are reported in Figure 2. The campus is located far enough from the city to avoid airflow obstructions and shading, which would reduce the renewable energy share.As a whole, the 1000-student campus occupies a surface of 1 × 0.5 km 2 and includes five dormitories, a dining hall, a gym, a students' center with classrooms and administrative offices.A schematic representation of the campus is shown in Figure 3.The campus is located far enough from the city to avoid airflow obstructions and shading, which would reduce the renewable energy share.As a whole, the 1000-student campus occupies a surface of 1 × 0.5 km 2 and includes five dormitories, a dining hall, a gym, a students' center with classrooms and administrative offices.A schematic representation of the campus is shown in Figure 3.
Further details on Equations (13-16) can be found in [21].This model represents a good tradeoff between simplified models (e.g., the energy signature method [22]), which simply correlates external temperature and heating/cooling load, and dynamic models (e.g., TRNSYS or EnergyPlus), which include the building inertia characteristics, solar radiation, and internal loads, providing more accurate results, but requiring a detailed knowledge of the building envelope and heat gain profiles.

Case Study
In this work, we refer to an integrated energy system serving a hypothetical campus, located in Trieste, Italy.This city has a favorable climate, where RES (solar thermal, photovoltaic modules, wind turbines) can provide a significant amount of energy.The Italian Thermotechnical Committee (CTI) provides hourly profiles of external temperature, global solar irradiance on the horizontal plane and wind speed [23].The monthly-average values of the external temperature and irradiance on horizontal plane are reported in Figure 2. The campus is located far enough from the city to avoid airflow obstructions and shading, which would reduce the renewable energy share.As a whole, the 1000-student campus occupies a surface of 1 × 0.5 km 2 and includes five dormitories, a dining hall, a gym, a students' center with classrooms and administrative offices.A schematic representation of the campus is shown in Figure 3.All the buildings have a similar structure in terms of thermal transmittance of the walls, roofs, floors, and windows (0.29 W/m 2 •K, 0.19 W/m 2 •K, 0.26 W/m 2 •K, 1.80 W/m 2 •K, respectively).Specific profiles of internal gains and electrical energy requirements were chosen for each building type, according to their use, periods of presence, and appliances.The internal gains and electrical energy requirements have lower values during weekends and holidays.The buildings have different terminal units that require two supply temperature levels (i.e., low and medium).The power peak and total energy for the heating, cooling, domestic hot water (DHW) and electrical energy services are reported in Table 1.In this work, we compare a "centralized" configuration (see Figure 4a), in which the DHN is fed by a cogeneration unit and a centralized gas boiler, with a "distributed" one (see Figure 4b), in which reversible heat pumps are installed in the buildings.In both configurations, the sizing of the solar and wind generators remains the same, while an optimization analysis is performed for the CHP and heat pumps, together with the optimal control strategy (see Section 4).All the buildings have a similar structure in terms of thermal transmittance of the walls, roofs, floors, and windows (0.29 W/m 2 •K, 0.19 W/m 2 •K, 0.26 W/m 2 •K, 1.80 W/m 2 •K, respectively).Specific profiles of internal gains and electrical energy requirements were chosen for each building type, according to their use, periods of presence, and appliances.The internal gains and electrical energy requirements have lower values during weekends and holidays.The buildings have different terminal units that require two supply temperature levels (i.e., low and medium).The power peak and total energy for the heating, cooling, domestic hot water (DHW) and electrical energy services are reported in Table 1.In this work, we compare a "centralized" configuration (see Figure 4a), in which the DHN is fed by a cogeneration unit and a centralized gas boiler, with a "distributed" one (see Figure 4b), in which reversible heat pumps are installed in the buildings.In both configurations, the sizing of the solar and wind generators remains the same, while an optimization analysis is performed for the CHP and heat pumps, together with the optimal control strategy (see Section 4).The generators concur to satisfy the heating/cooling/DHW/electrical energy requirements of the campus with the following strategy:

•
The electrical demand consists of the energy input of the electrical appliances, DHN circulation pumps, and heat pumps or chillers.Electrical energy is purchased from the grid (E P ) when the power generated by the photovoltaic modules (E PV ), the wind turbine (E WT ) and ICE (E ICE ) is not sufficient to satisfy the requirements.On the contrary, if the electrical production is higher than loads, the overproduction is sold to the grid (E S ).

•
The DHW service (Q HT,D ) is required only for dormitories and the gym.The thermal storages are heated by the solar thermal panels (Q ST ), and by the DHN or back-up boilers (Q B ) when the temperature drops below the setpoint value (i.e., 55 • C).In the centralized configuration, the DHN represents the back-up generator, while in the distributed configuration both DHN or local boilers can be used (see Figure 4).The thermal storages are only used for the DHW service.

•
The heating service is required at two different levels of temperature (low and medium, Q LT,D and Q MT,D ).In the centralized configuration, it is delivered by the DHN, in the distributed layout it is provided by the DHN or by the local heat pumps (Q HP ).

•
Also, the cooling service is required at two different levels of temperature (C LT,D cooling requirement through radiant panels, C MT,D cooling requirement through fan-coils).In the centralized configuration, air-to-water electrically-driven chillers are used whereas in the distributed configuration, reversible heat pumps are used.

•
In the distributed configuration, the reversible heat pumps are in the buildings and provide both heating and cooling service at all temperature levels.
Table 2 reports all the parameters of the analyzed smart multi-energy microgrid.

Optimization Problem and Methodology
In this work, we aim to compare the centralized and the distributed configurations defined in Section 3, from an economical point of view.The selected performance index is the annual total cost, TC, defined as: where I NV/t li f e , O&M, and OC are the yearly capital, maintenance, and energy operational costs, respectively, and t li f e is the considered lifetime of the microgrid (i.e., 20 years).Since we are making a comparative analysis, TC only includes the costs that differ in the two configurations, i.e., the purchased equipment cost (PEC) for the ICE unit, reversible HPs, chillers, and boilers, the associated operations and maintenance costs (O&M), and the net cost of the energy purchased from the gas and power grids.
Energies 2019, 12, 745 8 of 17 We do not consider the sizing of the RES technologies, which are assumed to have the same design and energy production in both the distributed and centralized configurations; therefore, they are not included in the economic analysis and optimization process.The terms in Equation ( 17) read: The timestep length adopted for the energy system simulation is one hour.The investment and maintenance cost functions are presented in Table 3.The costs of HPs, chillers and boilers were obtained through a linear regression of actual manufacturers' data.Manufacturers data and [27,28] According to the operating strategy described in Section 3, both the sizing and control optimization can be written as a function of the ICE nominal electrical capacity, E nom ICE , and load factor profile, L i ICE . 1.

2.
For each n-th ICE size, an exhaustive search is performed to find the ICE load factor, L i ICE,j , which minimizes the operational energy cost at any i-th timestep, OC i , (i.e., control optimization).Further details on this point are provided in Section 4.1.

3.
The optimal L i ICE,n sequence determines the thermal and electrical output profiles of the ICE.Subsequently, also the other generators output can be easily found as the residual load.The maximum output times a precautionary factor of 1.1 gives the nominal capacity of each generator.

4.
Capital and annual maintenance costs are evaluated for any n-th sizing and corresponding optimal control sequence.Then, the total annual cost is evaluated.

5.
Finally, we selected the best sizing and corresponding control strategy as the one with the minimum total annual cost, evaluated at point 4.

Operational Optimization Problem
The operational optimization problem consists in identifying the scheduling of the generators that meets the energy demands at minimum cost (i.e., cost for purchasing electricity from the grid, income for selling electricity to the grid, cost of natural gas).The operational optimization problem is therefore defined as the minimization of the total annual energy cost.
where i = 1, . . ., 8760 timesteps, c F is the fuel price (0.04 €/kWh), and c el,P and c el,S are the prices of purchased and sold electricity (0.18 and 0.04 €/kWh), respectively.Consequently, the following decision variables are considered: and demand constraints and balance equations and inequalities are defined as follows To solve the optimal operation problem, an ad-hoc dispatch strategy algorithm has been developed, based on the following considerations: • the problem can be considered "static"; • three orders of priority must be considered for the DHN dispatch.
Indeed, in the energy system under investigation, the overall optimum coincides with the sum of optimums of every single timestep, since the behavior of the TSs (Thermal Storages) linked to the solar thermal is independent from the operational control.For this reason, as already shown in [29], the overall operational problem can be split into 8760 subproblems, one for each timestep, and the problem can be considered "static".Moreover, the DHN must always satisfy with higher priority the medium-temperature heat demand Q i MT,D , as opposed to the low-temperature heat demand Q i LT,D , since the HPs operate with higher COP at lower supply temperatures.Therefore, only the three following combinations must be evaluated: (1) HT-MT-LT (2) MT-HT-LT (3) MT-LT-HT For each possible dispatch priority, once the L ICE is set and the amount of electricity and heat produced by the ICE is defined, the thermal losses and the net amount of heat available at the DHN are known from Equations ( 23) and (24).Then, Equations ( 25)-( 27) state that the boiler and the HP production must meet the remaining heat demand, if any.Furthermore, Equations ( 27) and ( 28) require that the electrical chillers or reversible heat pumps meet the chilled water demand, and Equation ( 29) defines the electrical energy exchange with the grid.Therefore, as mentioned, the nine decision variables are bound to each other and the problem is conveniently reduced to finding the dispatch priority order and the optimal ICE load factor L ICE that minimize the cost of energy at each timestep.Eleven discrete values of L ICE have been considered, and an exhaustive search algorithm was adopted to identify the optimal solution, among all the possible combinations.This allows the development of a low computational-cost algorithm, compatible with the need for a quick response for real-time implementation [30] and further advanced analyses (e.g., optimal design and uncertainty analysis, as in [29]).

Results and Discussion
Figure 5 shows the total annual costs depending on the nominal electrical capacity of the CHP unit.We note that the optimal CHP size is practically the same for both of the configurations (75 kW el ), Energies 2019, 12, 745 10 of 17 with a total cost reduction of about 8% for the distributed solution.Economic details are presented in Table 4.The size of all generators are shown in Table 5.   Figure 6 shows some examples of the heat and power profiles resulting from the optimization procedure in three weeks of the year for the distributed configuration.In addition to the areas and lines explained in the chart legend, we specify that the white area under the blue curve in the thermal plots represents the thermal overproduction by the CHP unit; the white area under the green, orange, or yellow lines in the electricity plots quantifies the electrical energy sold to the grid.Figure 6 shows some examples of the heat and power profiles resulting from the optimization procedure in three weeks of the year for the distributed configuration.In addition to the areas and lines explained in the chart legend, we specify that the white area under the blue curve in the thermal plots represents the thermal overproduction by the CHP unit; the white area under the green, orange, or yellow lines in the electricity plots quantifies the electrical energy sold to the grid.
demand of the appliances.The HPs are rarely used as it is more convenient to use the boilers to meet the few thermal peaks with respect to purchasing more electrical energy from the grid.In summer, the CHP still operates at high load factors to meet the residual electrical load from the variable solar and wind production.The thermal overproduction mainly occurs during the mid and summer seasons because it is more convenient to maximize the electrical production of the CHP unit, with respect to purchasing electricity from the grid.The annual values reported in Figures 7 and 8 highlight the energy and economic advantages of the distributed configuration.The distributed configuration reduces both the DHN thermal losses and the CHP overproduction (-27% and -9%, respectively), thus, reducing the overall thermal energy production (-6%).The heat pumps deliver about 24 % of the heating load and the boilers reduce their output to about 68%.In the centralized configuration, the CHP unit delivers 77% of the thermal load, The CHP delivers the base load during the winter season and the heat pumps meet the heating demand.The boilers are rarely used as the electricity produced by the renewables sources and CHP is enough to run the heat pumps.During the mid-seasons, the CHP mainly meets the high electricity demand of the appliances.The HPs are rarely used as it is more convenient to use the boilers to meet the few thermal peaks with respect to purchasing more electrical energy from the grid.In summer, the CHP still operates at high load factors to meet the residual electrical load from the variable solar and wind production.The thermal overproduction mainly occurs during the mid and summer seasons because it is more convenient to maximize the electrical production of the CHP unit, with respect to purchasing electricity from the grid.
The annual values reported in Figures 7 and 8 highlight the energy and economic advantages of the distributed configuration.The distributed configuration reduces both the DHN thermal losses and the CHP overproduction (−27% and −9%, respectively), thus, reducing the overall thermal energy production (−6%).The heat pumps deliver about 24 % of the heating load and the boilers reduce their output to about 68%.In the centralized configuration, the CHP unit delivers 77% of the thermal load, being more cost-effective than the gas boiler.However, the corresponding total electricity production (CHP, PV, and wind generator) is significantly greater than the electrical load, resulting in an unprofitable overproduction (about 28% of the total electricity production is sold to the grid at an uneconomical price).On the contrary, the introduction of the heat pumps shifts part of the thermal load to the electrical one, increasing the self-consumption of electricity and reducing both the CHP energy production (−20% of thermal and −21% of electrical energy, respectively) and the power sold to the grid (−51%).In total, the presence of the HPs reduces the power exchange with the electrical grid from 259 to 159 MWh (−40%) considering both sold and purchased quantities.
Although the objective function of the optimization process refers to an economic index, the energy and the environmental benefits of the distributed configuration are shown in Table 6.The net no-RES primary energy consumption was reduced by about 8%, and the equivalent CO 2 emissions were reduced by about 11%.These values were evaluated considering the primary energy factors and the specific CO 2 emissions of the Italian energy systems [31,32].According to a grid perspective, the energy system can be thought of as an electricity generator with a specific emission of 232 g/kWh in the distributed configuration and 238 g/kWh in the centralized one.We note that the average value of specific CO 2 emissions for Italian power production is about 313 g/kWh.Overall, the results show that the proposed hybrid centralized-distributed configuration outperforms the more conventional centralized configuration from an economic, environmental, and efficiency perspectives.Indeed, the introduction of heat pumps at the building level enhances the operational flexibility of the system by enabling the interconnection between the thermal and electric networks.In this way, RES-based energy production can be used mainly on-site-instead of being sold to the regional grid-and the use of inefficient technologies, such as natural-gas boilers, can be drastically reduced.

Conclusions
An innovative configuration for smart multi-energy microgrids serving clusters of buildings has been presented.The energy system combines both centralized and distributed generation units, optimally integrating cogeneration-based micro-district heating, RES technologies, and reversible heat pumps.
The proposed system was tested in a hypothetical case study, namely, a University Campus located in Trieste (Italy).A detailed modeling of the building load demands, district heating network, and all energy units has been provided in order to simulate the energy system in a reference-year scenario.Moreover, an operational optimization algorithm was specifically developed to identify the generator scheduling that meets the energy demand while minimizing the operational cost.The optimal size of the cogeneration unit and reversible heat pumps has also been found.
The proposed configuration was compared to a more conventional layout based completely on centralized heat production.The results show how the introduction of distributed heat pumps to assist the thermal production at the building level enhances the flexibility and cost-effectiveness of the energy system.Indeed, an 8% total-cost saving, 11% carbon emission reduction, and 8% primary energy saving were achieved compared to the centralized reference case.Moreover, the proposed configuration significantly reduced the electric energy exchange with the regional grid (around 40% less).
Future work will address the current limitations of the work: the optimal sizing of the whole system will be investigated, the effect of uncertainty in weather conditions and economic parameters will be analyzed, and the effectiveness of energy storage managed by predictive control will be evaluated.

Figure 1 .
Figure 1.Schematic of the energy system.

Figure 2 .
Figure 2. Average monthly temperature and daily irradiance on horizontal plane.

Figure 2 .
Figure 2. Average monthly temperature and daily irradiance on horizontal plane.

Figure 2 .
Figure 2. Average monthly temperature and daily irradiance on horizontal plane.

Figure 3 .
Figure 3. Scale representation of the campus: buildings, district heating, and generation systems.

Energies 2019 , 18 Figure 3 .
Figure 3. Scale representation of the campus: buildings, district heating, and generation systems.

Figure 5 .
Figure 5. Installation, energy, and maintenance annual costs using the optimal control strategy.
(a)1st week of January (b)1st week of January (c)1st week of April (d)1st week of April

Figure 6 .
Figure 6.Heat and power profiles in 3 weeks of the year (areas are for consumption, lines for production).

Figure 6 .
Figure 6.Heat and power profiles in 3 weeks of the year (areas are for consumption, lines for production).

Table 1 .
Peak values and energy needs for the four services of the campus and for the five types of buildings.

Table 1 .
Peak values and energy needs for the four services of the campus and for the five types of buildings.

Table 2 .
Parameters used in the analysis.

Table 4 .
Annual costs of the optimal distributed and centralized configurations, k€/yr.

Table 5 .
Optimal sizes for the distributed and centralized configurations, kW.

Table 4 .
Installation, energy, and maintenance annual costs using the optimal control strategy.Annual costs of the optimal distributed and centralized configurations, k€/yr.

Table 5 .
Optimal sizes for the distributed and centralized configurations, kW.

Table 6 .
Energy and environmental performance indexes.An innovative configuration for smart multi-energy microgrids serving clusters of buildings has been presented.The energy system combines both centralized and distributed generation units, Electrical energy balance (MWh/yr).

Table 6 .
Energy and environmental performance indexes.