A Comparison of Di ﬀ erent District Integration for a Distributed Generation System for Heating and Cooling in an Urban Area

: The paper proposes a comparison of di ﬀ erent district integration options for a distributed generation system for heating and cooling in an urban area. The system considered includes several production units located close to the users, a central unit and the district heating and cooling network which can connect all the users to each other and to a central unit, where a cogeneration system and a solar plant can be placed. Thus, each user can be regarded as isolated from the others, satisfying its energy needs by means of an autonomous production unit. Alternatively, it can be connected to the others through the district heating and cooling network. When a district heating and cooling network is included in the design option the synthesis-design and operation problems cannot be solved separately, because the energy to be produced by each production site is not known in advance, as the ﬂows through the district heating and cooling network are not deﬁned. This paper uses a mixed integer linear programming (MILP) methodology for the multi-objective optimization of the distributed generation energy system, considering the total annual cost for owning, operating and maintaining the whole system as the economic objective function, while the total annual CO 2 emissions as the environmental objective function. The energy system is optimized for di ﬀ erent district integration option, in order to understand how they a ﬀ ect the optimal solutions compared with both the environmental and economic objects.


Introduction
The reduction of pollutant emissions is one of the current main targets fixed by international authorities. A lower energy need in the residential and tertiary sectors can help to achieve this goal, as it represents one of the dominant energy consuming sectors in industrialized societies. However, the adoption of a defined energy system still depends on technical and financial evaluations, while environmental aspects are not generally regarded as design goals.
A review of the open literature on these topics shows that the current research works can be grouped into three major groups: • research focusing on the optimization of the operation of energy systems, ranging from the optimization of a single component, to the operation of the overall DG system; • research dealing with the optimization of the system synthesis; and • research focusing on synthesis, design and operation optimization.
Each group can be further subdivided, considering single and multi-objective optimization targets. Over the last decade an increasing number of papers dealing with energy system optimization have been produced [1]. One of the first optimization models was developed by Henning in 1992 [2], and consists of a linear programming model to minimize the operating cost of an energy supply system for local Swedish utilities. In 1997 he presented a linear programming model called MODEST [3] for the minimization of capital and operation costs of energy supply and demand side management. Curti et al. [4] proposed an optimization model for aiding the design of a mixed energy production system, including heat pump based district heating, conventional boilers and decentralized heat pump. Yokoyama et al. [5] in 2002 proposed a method for optimal structural design, to determine the structures of energy supply systems in consideration of their multi-period operation. Karlsson [6] has recently presented the MIND method, a decision support for optimization of industrial energy systems.
For a general overview of models, methods and applications of multi-utility energy systems, the authors refer readers to extended reviews presented recently [1,[7][8][9][10]. Among a large quantity of research presented in this field, some common points can be outlined: • almost all models rely on linear programming or mixed integer linear programming (MILP). However, some approaches based on meta-heuristics (simulated annealing, genetic algorithms, etc.) have been proposed, but they present some difficulties concerning the determination of search parameters and the judgment about optimality [11][12][13]. • the research normally focus only on specific targets: operation or synthesis optimization, economic and/or environmental optimization, unit or district heating network (DHN) optimization, etc.
To deal with optimization of DG energy systems, including DHCN and thermal storage, and focusing on different objectives (economic rather than environmental), it is necessary to consider all aspects at the same time, and not in successive steps. This is because the operation optimization strongly affects the optimal synthesis of the system and, in addition, the economic optimum does not correspond to the environmental one.
Some recent papers seem to go in this direction, performing a single objective optimization, generally economic: Chinese proposed a MILP model for the optimization of a DHCN in a DG context [14], Soderman and Petterson [15] presented a structural and operational optimization of a DG energy system. Pavicevic et al. [16] performed the optimization of sizing and operation of a DH system, focusing on the technological options for the TS. Pérez-Mora et al. [17] optimized DHC systems, focusing on the alternative between absorption and compression chillers. Ameri and Besharati [18] defined a model for determining the optimal capacity and operation of seven combined cooling, heating and power (CCHP) systems in the heating and cooling network of a residential complex located east of Theran.
Finally, Ren et al. [19] proposed a multi-objective optimization model to analyse the optimal operating strategy of a distributed energy system, while combining the minimization of energy costs with the minimization of environmental impact, which is assessed in terms of CO 2 emissions. Carvalho [20] presented a model for the synthesis and operation optimization of residential units, considering environmental and economic aspects.
In this paper, the optimization MILP model is applied to a real case study, made up of nine tertiary sector users located in a small town city centre situated in the northeast of Italy. A preliminary energy audit allowed the determination of the users' energy needs, so that the latter is regarded as the input of the optimization procedure.
This paper proposes a comparison of different district integration option for a distributed generation (DG) system for heating and cooling in an urban area. The distributed generation energy system considered in the paper includes several production units located close to the users, a central unit and the district heating and cooling network (DHCN) which can connect all the users to each other and to a central unit, where a cogeneration system and a solar plant can be placed. Thus, each user can be regarded as isolated from the others, satisfying its energy needs by means of an autonomous production unit. Alternatively, it can be connected to the others through the DHCN. In this case, it

Decision Variables and Constraints
The optimization variables to be considered in order to determine the optimal configuration and operation of the DG energy system can be divided into two main groups: • binary variables: they represent the existence/absence of each component and the operation status (on/off) of each component in each time interval. There are other additional binary variables which do not represent any physical quantity, added to linearize some relations; and • continuous variables: they represent the size of components, the size of pipelines, the load of components in each time interval, the energy content of the storages and the connection flows.
The performance characteristics of equipment and energy balance relationships are the fundamental constraints of the optimization problem. Other constraints, such as relationships between maximum contract demands and consumption of energy purchased and operational restrictions, have also to be considered. In the MILP model, equality constraints express fix relations and balances, while inequality constraints, express limits and feasibility conditions. The constraints can be generally grouped into four categories which describe: • components; • district heating and cooling network; • thermal storage; and • energy balances.

Components
All cogenerators, heat pumps and absorption chillers included in the superstructure are of fixed size, while boilers, compression chillers, solar and photovoltaic panels are of variable size. Up to j components of the same fixed size can be adopted in parallel by the same user. The components of

Decision Variables and Constraints
The optimization variables to be considered in order to determine the optimal configuration and operation of the DG energy system can be divided into two main groups: • binary variables: they represent the existence/absence of each component and the operation status (on/off) of each component in each time interval. There are other additional binary variables which do not represent any physical quantity, added to linearize some relations; and • continuous variables: they represent the size of components, the size of pipelines, the load of components in each time interval, the energy content of the storages and the connection flows.
The performance characteristics of equipment and energy balance relationships are the fundamental constraints of the optimization problem. Other constraints, such as relationships between maximum contract demands and consumption of energy purchased and operational restrictions, have also to be considered. In the MILP model, equality constraints express fix relations and balances, while inequality constraints, express limits and feasibility conditions. The constraints can be generally grouped into four categories which describe:

Components
All cogenerators, heat pumps and absorption chillers included in the superstructure are of fixed size, while boilers, compression chillers, solar and photovoltaic panels are of variable size. Up to j components of the same fixed size can be adopted in parallel by the same user. The components of the central unit are all of variable size. The following constraints (Equations (1)- (5)) describe the linear model of ICE which can be installed in each production site.
The first set of constraints ensures that a consistent set of binary variables (X ice , O ice ) is taken into account in each time interval and throughout the year: the component j can be installed only if the component j − 1 has been already adopted (Equation (1)), and the component j can never be in operation if it has not been adopted (Equation (2)): The second group describes the partial load performance of energy conversion devices by means of linear relations: The coefficients Kh ice and Kf ice can be obtained through a linear regression of the load curves. The constraints which describe the MGT can be easily inferred by changing in each variable or coefficient the subscript "ice" with the subscript "mgt". A variable size ICE can be installed in the central unit. The constraints which describe this component (Equations (6)- (12)) are different from the previous constraints, because both size and load are decision variables. Therefore, it is necessary to introduce additional constraints and decision variables in order to maintain the linearity of the problem.
The first set of constraints limits the size of the ICE which can be adopted and put in relation the operation with the existence of the component: The second set of constraints expresses the relations among the main product (E ice,c ), the sub-product (H ice,c ) and the fuel flows (F ice,c ): The last set of equations is required to constrain the additional variables ξ ice,c (m,d,h), which allow us to introduce a linear relation with two independent variables-size S ice,c and load E ice,c (m, d, h)-avoiding inconsistent results when the engine is off: Appl. Sci. 2019, 9, 3521 6 of 40 The boiler which can be installed in the central unit is described by the following equations (Equations (13)-(16)), where a minimum load limit (H boi,lim,c = 0.1) has been taken into account and auxiliary variables ψ boi,c (m, d, h) play a role analogous to variables ξ ice,c (m, d, h), previously introduced for the ICE: The boilers and the compression chiller installed in the production sites have been modelled by means of the efficiency and the COP alone, without a minimum load limit.
Absorption chillers are modelled like local ICE and GT, as components of fixed size which can be installed as multiple units operating in parallel. Each ABS operation is allowed only if the heat produced by ICE and MGT, in the same site, is greater than the heat required by the former and the ABS existence is constrained by the existence of ICE and MGT, as expressed by Equations (17) and (18): The heat pump is modelled as a black box that can produce either heat-O hp,h (m, d, h, j, u) = 1-or cold-O hp,c (m,d,h,j,u) = 1, but not simultaneously: Equations (23) and (24) refer to the characteristic curve of the heat pump in the heating and cooling mode, respectively. The solar thermal collectors and the photovoltaic panels are modelled considering their production proportional to the size of the plants. The unitary production-K stp (m, d, h, u) and K pv (m, d, h, u)-is evaluated a priori considering inclination, orientation angle of installation and hourly solar radiation of the site of the plant.

District Heating and Cooling Network
The modelling of the DHCN is important for the optimization of the DG energy system, because it strongly affects the optimal solution [25][26][27]. The aim of the DG energy system optimization, is to define the lay-out of the DHCN and the dimension of each single pipeline, taking into account the operation of the whole system. The heat which can be transferred by a DHCN pipeline can be expressed by: Appl. Sci. 2019, 9, 3521 7 of 40 Assuming that the velocity v p is fixed, the transferred heat . Q p reported in Equation (28) depends on two variables, which are the cross section of the pipeline A p and the temperature difference between inlet and outlet pipelines ∆t [28]. Assuming a fixed temperature of the network and a fixed temperature difference between the inlet and outlet, the modelling introduces a constant ratio between the size of the pipeline and the maximum flow which can be transferred, and considers the size and the layout of the network as decision variables, constrained by the pipeline super-structure and the flow rate limits of each pipe. The ∆t adopted normally ranges between 15-25 • C depending on the application, while the medium velocity v p ranges between 1.5-2.5 m/s. The thermal losses are considered proportional to the length of each pipeline through the coefficient δ t : Equations (30)-(32) describe the existence conditions of thermal pipelines, the energy flow through each pipeline during operation and the maximum flow characterizing the size of each pipeline: The same modelling can be obtained for the district cooling network, by changing properly the subscripts.

Thermal Storage
A long-term thermal storage might be able to extend the operation time of the combined heat and power (CHP) unit. In fact, during summer time, for instance, the heat demand can be so low that the CHP plant must be shut down and a boiler, which is often expensive in terms of operational cost, must be brought into operation. By using a thermal storage, power generation can be increased and the use of fossil fuels can be reduced [29]. The problem of intermittent energy sources is especially severe for solar energy, because thermal energy is usually needed most when solar availability is lowest, namely, in winter. Small TS can cover periods of inadequate sunshine, while large TS, operating over long period, can partially cover the lower winter solar thermal production [30]. In any case, the design of the integrated system is very difficult and normally several rules of thumb are used to define the system. An optimization procedure can be very helpful to reach a full exploitation of the potential benefit of CHP and solar systems.
The thermal and cooling storages which can be installed in each production unit or in the central unit can be modelled in the same way, accepting the approximations of perfect stratification of medium (water) inside the thermal storage. This approximation correspond to the hypothesis that, if the storage is not completely discharged, the residual energy is stored at the same temperature required by the DHN. The energy stored in the thermal energy storage, can be evaluated through: As well as for the DHCN, the temperature difference ∆t between inlet and outlet temperature is considered constant. Therefore, the thermal energy stored in the thermal storage is proportional to the volume of the medium inside the storage and it is considered as a decision variable. A set of equations is required to describe the energy balance of the thermal storage. In order to allow for the seasonal charging/discharging cycle, the thermal storage has to be modelled throughout the whole year, without any time decomposition as can be done for the other components, where a set of similar days can be represented by only one typical day. Applying the typical day approach to all operation variables, except for the variables which represent the energy stored in the thermal storages, allows us, in any case, to reduce the overall number of variables, keeping the option of properly representing the charge and discharge phases of the thermal storages during all the year.
In the proposed model, the year is decomposed into 24 typical days of 24 h, one typical working day and one typical non-working day each month. Therefore, each single month is made up of four similar weeks, in turn made up of five similar working days and two similar non-working days. In this way, the optimal operation of the system is similar in each working or non-working day of the month and each month of the year is made up of 28 days. This approximation can be accepted selecting particular typical days which describe the whole year producing the same total consumption. Specific procedures for the identification of the proper typical days can be found in [31,32].
The energy balance of the thermal storage is approximated considering that the energy contained in the storage in a general time interval t is equal to the energy stored in the time interval t-1 multiplied by a thermal loss coefficient plus the input energy in the time interval t: Additional constraints have to be added to connect the storage's condition at the end of a time period with that at the beginning of the following one. For example, Equation (35) connects two days of the same kind (some adjustment has to be introduced to take into account the transition from a working day, to a non-working day, or between different weeks, months, etc.): Finally, the heat stored inside the storage Q ts , has to be lower than storage size S ts :

Energy Balances
Energy balances are a set of constraints which ensure that in each node and in each time interval the input energy is equal to the output energy. With reference to the superstructure represented in Figure 1, there are nodes three for each site (electric, thermal and cooling) plus the two nodes of the central site (electric and thermal). The energy balance of the network is included in the node energy balance. In the following the electric balance for a site u is described in detail, while the thermal and cooling balances are reported with no description.
The electricity and the thermal balances of a site are expressed by: The thermal production of the central unit can be shared among the users, which receive heat from the DHN through node u. Therefore, the thermal balance related to the central unit site reads: All continuous variables have to be greater than zero, except the variables related to the thermal storage input/output heat flow (H ts,c , H ts , C ts ) which are free. Positive values represent input flows, while negative values mean an energy extraction from the thermal storage. The temperatures of the thermal flows are not taken into account because it would have compromised the linearity of the problem. However, this is not a strong approximation considering that the thermal energy required by the users is normally supplied at a temperature of 50-55 • C and all components are able to produce the thermal energy at higher temperatures. Some restrictions due to the coupling of components related to the operating temperatures (ICE, MGT together with ABS) have been considered through a particular conformation of the superstructure and with additional constraints which consider this matter (e.g., Equations (17) and (18)).

Objective Functions
The economic objective function to be minimized is the total annual cost for owning, operating and maintaining the whole system: The annual cost for the investment (c inv ) is the sum of the investment cost of the sites, of the central unit, and of the network. The investment cost of a site can be evaluated through: c inv,u (u) = Σ j [f mgt ·X mgt (j,u)·c mgt (j,u) + f ice ·X ice (j,u)·c ice (j,u) + f hp ·X hp (j,u)·c hp (j,u) + f abs ·X abs (j,u)·c abs (j,u)] + f boi ·S boi (u)·c boi + f cc ·S cc (u)·c cc + f pvp ·S pvp (u)·c pvp + f stp ·S stp (u)·c stp + f ts ·S ts (u)·c ts + f ts ·S cs (u)·c ts (45) The investment cost of the central unit is: The investment cost of the network is: Maintenance costs related to the components are considered proportional to the products, while the operation costs include the costs for fuel, electricity bought from the grid and the eventual income from the sale of electricity. For each site u the annual operation cost is evaluated as: The other objective function to be minimized is the annual operation greenhouse emissions (CO 2 emissions). The total annual emissions are related to the net electric energy received from the grid (bought minus sold) and to the fuel consumption by the DG energy system (boilers and/or CHP). Therefore, the total annual emissions can be evaluated through: The greenhouse emissions due to the consumption of fuel by the DG system are mainly related to the kind of fuel itself [33], while the electricity carbon intensity largely depends on the technology mix of each electricity supplier.
These two objectives are conflicting ones because the adoption of environmental efficient energy systems are costly. Likewise, the solution which allows the minimum annual cost does not permit us to obtain minimum total annual operation emissions.
The ε-constrained method [34] has been adopted to obtain the Pareto front solutions. First, the economic and environmental optimal solutions have been obtained. Secondly, the difference between the environmental objective in the two optimal cases has been calculated. Thirdly, a set of predefined intermediate emission levels has been identified and each level has been introduced as an additional constraint, in a further economic optimization. In this way, a set of dominant solutions, on the Pareto front of the multi-objective optimization (economic and environmental), can be obtained.

Case Study
A systematic approach for the selection of an appropriate DG energy supply system requires a detailed knowledge of heat, cooling, and electricity user demands. The detail level affects the model complexity and one of the factors which plays against the model compactness is the number of time intervals considered. These time periods are defined by the number of different energy demands that have to be covered, and the periodicity considered in the model (hourly, weekly, monthly). Long time periods, such as weeks or months, can be considered for industrial applications, characterized by quite constant energy demands that are fairly independent of environmental conditions. Hourly energy demand data is very important when tertiary sector/residential energy systems are analysed, where the influence of environmental conditions is quite important. In this last case, one solution to contain the model complexity is to represent the whole year through some typical days [32]. In the current optimization case study, the whole year has been subdivided into 24 typical days made up of 24 h each: 12 typical days refer to working days, while the remaining 12 refer to non-working days, so that each month is represented by one working and one non-working day. All values related to each single time interval are weighted through the parameter wgt, which consider their weights in the overall year. The grouping through typical days can be done for all variables, with the exception of the variables related to the thermal storages (either thermal or cooling) for which the whole year has to be considered (detailed explanation can be found in Section 2).
This case study is made up of nine tertiary sector users located close to each other in the centre of a small city (60,000 inhabitants) in the northeast of Italy. The users considered are all owned by the public service and this gave us the possibility to access the energy demand data. The users considered are the town hall, a theatre, a library, a primary school, a retirement home, the archive, a hospital, a secondary school and a swimming pool ( Figure 2). The heterogeneous choice of the buildings under consideration, characterized by different kinds of energy demands, allows us to consider the achieved results not affected by a specific user profile. Furthermore, a similar mix of users is expected to be easily recognized in a lot of other small and medium-sized towns in Europe. Figure 2 shows the location of the nine buildings involved in the study and of the central unit together with solar field. In addition, the possible path of the DHCN is outlined in red.
The determination of the possible paths is the result of a preliminary study which considers: • layout of the roads which connect the buildings; • position of the underground utilities (waterworks, sanitation, gas network, etc.); and • location of the boiler rooms of the buildings.
The users are close to each other and the maximum distance, between user 1 and the central unit, is about 2.5 km.
The users considered require thermal energy for space heating and for sanitary hot water, which is conventionally produced by boilers. The thermal energy for space heating is required at a temperature of about 65-70 • C. The cooling energy is required only for space cooling typically during the summer season and it is produced by compression chillers. In the optimizations, the design values of the temperatures for the operation of the DHN and DCN have been fixed at 70 • C and 12 • C, respectively. At the moment, the electric energy is bought from the grid and covers both the electricity demands and the electricity required to power the mechanical chillers. Table 1 reports the annual energy consumptions and peak power of the nine users. The hospital requires about 50% of the electric energy consumption, the second energy consumer is the swimming pool, while the other users require less than 7% each. Similar situations can be found for thermal and cooling demand, with the exception that cooling energy is not required by the schools, as in summer there are no students and the cooling plants have not been provided. The last row of Table 1 shows the sum of all power peaks, and is, as can be expected, greater than the total power peak by about 10%. The electric demand shown in the table does not account for the energy required by the compression chillers for the production of cooling energy.  Figure 3 shows the trend of electric, thermal and energy demands for all buildings. It shows clearly that the trends are characteristic of tertiary sector users and also of a typical European continental climate, where during winter the thermal demand is higher than in summer period because of the need for space heating. During summer the thermal demand is only due to hot sanitary water demand, while there is a requirement of cold for space cooling. The electric energy is slightly higher during winter and lower in summer. This trend is related to the daylight hours, as during the summer there is a lower need of lighting.
because of the need for space heating. During summer the thermal demand is only due to hot sanitary water demand, while there is a requirement of cold for space cooling. The electric energy is slightly higher during winter and lower in summer. This trend is related to the daylight hours, as during the summer there is a lower need of lighting. The energy demand patterns of the hospital in winter and summer are shown in Figure 4, for two representative working days. The winter trend shows greater electric energy demand during the daylight hours due to a higher occupancy factor, a heating demand higher in the morning and in the evening and a very low cooling demand, which is required by the air conditioning system. In summer, the cooling demand is higher compared with the winter season and reaches a maximum at about 3 p.m., the electric demand is similar to that of winter, while the thermal demand employed as sanitary water is higher during daylight hours and lower during the night. Each building has The energy demand patterns of the hospital in winter and summer are shown in Figure 4, for two representative working days. The winter trend shows greater electric energy demand during the daylight hours due to a higher occupancy factor, a heating demand higher in the morning and in the evening and a very low cooling demand, which is required by the air conditioning system. In summer, the cooling demand is higher compared with the winter season and reaches a maximum at about 3 p.m., the electric demand is similar to that of winter, while the thermal demand employed as sanitary water is higher during daylight hours and lower during the night. Each building has different energy patterns which depend on the occupancy factor, thermal insulation, night lighting, etc.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 13 of 38 different energy patterns which depend on the occupancy factor, thermal insulation, night lighting, etc.
The type and the size of the equipment which is part of the system superstructure must be appropriate to allow their integration, and proportionate to the user energy demands. All components considered in the optimization are commercially available and the prices have been obtained through a market survey. Two different kinds of components have been considered: fixed size components and variable size components. The optimal size/configuration of the energy system is obtained by defining the number of fixed size equipment installed and the size of variable size components (boilers, compression chillers, thermal storages). Table 2 shows the sizes of components which can be installed in the different units. These sizes have been chosen consistently with the peak demand expected for each units (Table 1) taking into account that up to six elements of the same kind can be adopted in parallel. The technical characteristics of all pieces of equipment are described in File S1, jointly with their investment costs. The latter have been obtained by a market survey and also include transportation, installation, connection, engineering costs, etc.
The amortization factors which multiply the investment costs are a function of the interest rate and of the life span of each component. The interest rate is assumed equal to 6% and is a sum of the real economic interest rate (4%) and a risk rate, assumed in this case equal to 2%. The life span of DHCN has been fixed equal to 30 years, for PV and ST panels and TS equal to 20 years, for ICE, MGT, ABS and HP equal to 15 years, and for BOI and CC equal to 10 years.
Maintenance costs are reported in File S1 and they are proportional to the energy produced by each component.  The type and the size of the equipment which is part of the system superstructure must be appropriate to allow their integration, and proportionate to the user energy demands. All components considered in the optimization are commercially available and the prices have been obtained through a market survey. Two different kinds of components have been considered: fixed size components and variable size components. The optimal size/configuration of the energy system is obtained by defining the number of fixed size equipment installed and the size of variable size components (boilers, compression chillers, thermal storages). Table 2 shows the sizes of components which can be installed in the different units. These sizes have been chosen consistently with the peak demand expected for each units (Table 1) taking into account that up to six elements of the same kind can be adopted in parallel. The technical characteristics of all pieces of equipment are described in File S1, jointly with their investment costs. The latter have been obtained by a market survey and also include transportation, installation, connection, engineering costs, etc. The amortization factors which multiply the investment costs are a function of the interest rate and of the life span of each component. The interest rate is assumed equal to 6% and is a sum of the real economic interest rate (4%) and a risk rate, assumed in this case equal to 2%. The life span of DHCN has been fixed equal to 30 years, for PV and ST panels and TS equal to 20 years, for ICE, MGT, ABS and HP equal to 15 years, and for BOI and CC equal to 10 years.
Maintenance costs are reported in File S1 and they are proportional to the energy produced by each component.
Operation costs are related to the costs of fuels and electricity. The Italian gas and electricity market has been liberalized since 2007, after a process which began in the 1999 with the "Bersani Decree" [35] and lasted eight years. Since 2007 all consumers can freely choose a supplier and leave the regulated-rate system or can remain connected to the old regulated market. Herein, the regulated-rate system has been considered as reference for the electricity and gas prices. The natural gas cost has been considered constant all year long, as well as the price of the bought electricity, while the price of the sold electricity has been assumed variable in each time interval, based on the hourly market prices. In particular, the bought electricity price has been considered equal to 0.17 €/kWh, while electricity sold price variable between 0.05-0.12 €/kWh, and natural gas price equal to 0.06 €/kWh. The prices considered include the pure cost of energy (about 40%) and taxes (about 60%).
The CO 2 emissions related to the consumption of electricity and natural gas have been assumed from literature [33]. The natural gas CO 2 emissions depend on the chemical composition of the gas and then on its provenience. However, the slight difference can be ignored considering the same value for the natural gas CO 2 emissions. The same approximation cannot be made for electricity. In fact, electricity carbon intensity depends heavily on the national electricity system. The reference case study has been optimized assuming the average electricity carbon intensity of the European Union in 2007-2009 (0.356 kgCO 2 /kWh), while a second set of optimizations has been performed assuming the average electricity carbon intensity of the OECD Americas (Canada, United States, Mexico, Chile) in 2007-2009 (0.485 kgCO 2 /kWh). The natural gas carbon intensive has been assumed equal to 0.202 kgCO 2 /kWh.
To take into account the incentives that many European countries have adopted to promote the development of new energy technologies, two different incentives have been considered for the reference case study: • natural gas detaxation for cogeneration use; and • renewable energy production incentives.
Cogeneration systems complying with restrictions defined by the Dlgs 20/07 [36] are classified in Italy as "high efficiency cogeneration systems" and can operate with natural gas detaxation. Therefore, the cost of the natural gas for cogeneration use in the case study has been considered equal to 0.045 €/kWh (25% less compared with the conventional price of natural gas).
The second incentive considered in the case study is applied to the electricity produced by renewable energy sources. The rates considered in the case study for electricity produced by PV panels are 0.199 €/kWh if put in the grid and 0.111 €/kWh if directly used by the user, in line with Italian "Quinto Conto Energia" [37].

Results of the Optimizations
The aim of the optimization is to determine the optimal configuration of a complex DG energy system together with the optimal operational strategy on an hourly basis throughout one year. The objective functions of the optimizations are the total annual cost for owing, operating and maintaining the whole system and the total annual operation CO 2 emissions. Using the -constrained method the Pareto fronts have been obtained for different plant configurations.
The model has been optimized to obtain the optimal solution in 5 cases, corresponding to different district integration options: distributed generation solution without central unit and district cooling network; 4.
distributed generation solution with central unit but without cooling network; and 5.
complete distributed generation solution.
As can be noted, the last configuration includes all the other configurations, while the conventional solution is a subset of all the other configurations. In this way it is possible to understand which is the influence of the different configurations and how they contribute to the achievement of the minimization of the objective functions.
The MILP model has been implemented in the X-press ® Optimization Suite. X-press ® is a commercial software produced by FICO ® (Fair Isaac Corporation, San Jose, CA, USA) for solving large optimization problems by means of the application of resident algorithms. The mathematical model has been implemented through Mosel, a modelling and programming language that allows users to formulate problems, to solve them by using the solver engines, and to analyse the solutions.
The optimization toolbox uses evolutionary algorithms, cut generations and heuristic algorithms, together with the branch and bound technique and revised simplex techniques. The branch and bound method starts with the optimization of the relaxed MILP problem, and explores the solution tree step by step, by fix one discrete decision variable at a time. In this way, the absolute optimal value of the objective function is not generally obtainable without exploring the complete solution tree, and a near optimal solution has to be accepted, stopping the optimization when a determined gap between the objective function of the relaxed problem and that of the current problem is reached. The case study has to be considered a very large problem as there are about 600,000 decision variables and 950,000 constraints. For this reason, the optimization procedure is stopped when a gap lower than 1% is reached.
The optimizations have been performed with a PC equipped with an Intel ® processor Core TM i7CPU 920@2.67 GHz, 6.00 GB RAM and a 64 bit operating system. An optimization of the overall problem, accepting a gap of 1%, takes about 100 h.

Conventional Solution
The Conventional Solution is considered as the reference solution, assuming that all thermal energy is produced by BOI, all cooling energy is produced by electric CC and all electricity is bought from the grid. Furthermore, the thermal and cooling energy can be stored in separated energy storages. The optimization is performed adopting a reduced superstructure where only BOI and CC can be installed ( Figure 5).     Table 3 shows the optimal configuration of the conventional solution, obtained by minimizing the total annual cost. The table reports the energy peaks of each single user. It can be noted that the boilers installed in each unit are of smaller sizes compared with the thermal peaks, as all units are provided with a proper sized thermal storage. On the other hand, the sizes of the compression chillers correspond to the cooling peaks as the cooling storages are not included in the optimal solution.  Table 4 shows the economic and environmental results of the optimization performed for the conventional solution. The total annual cost of the conventional solution is 2622 k€ per year, and it is made up of about 93% of operation costs (costs for thermal energy and electricity), 6% of investment costs and 1% of maintenance costs. The results show also that the hospital contributes to the total annual cost by about 50%, as it is the largest energy consumer. The table shows also the environmental results and that they depend basically on the energy consumption of each user. The energy balances are reported in Table 5. The thermal energy produced by BOI is slightly higher than the thermal energy required by the users, because of the heat losses in the TS. It should also be noticed that, as can be expected, the thermal and cooling energies produced by the components are all used by the users without any waste. The results presented in this section, and obtained by minimizing the economic objective function of the conventional solution, will be used as a reference for the evaluation of the forthcoming optimizations.

Isolated Solution
Two different superstructures have been considered for the optimization of the isolated solution: one including all components considered in Figure 6 and another including all those components, with the exception of the local thermal storages (either heating or cooling). These two different optimizations allow us to understand what the influence of the thermal storages in the optimal solution is. Pure economic and environmental optimizations have been conducted for the isolated solutions.
total annual emissions.
The environmental optimization shows an increase of the operation costs of about 60% compared with the economic optimizations, while it allows a reduction of about 25% of the total annual emissions with respect to the conventional solution. It can be also noted that the amount of emissions due to electricity usage significantly increases, while the emissions saved due to the electricity sold to the grid, and that in consequence of natural gas usage, significantly decrease.  The optimal configurations of the isolated solutions are compared to each other and to the conventional solution in Table 6. For matters of clarity, only the total power installed, for each kind of component, has been reported. Focusing on the economic optimizations, the ICEs have been adopted as cogenerators, while the MGTs are never installed. The ST panels are never adopted as well, because PV panels are more convenient and the space intended for the installation of both kind of panels is limited. The introduction of the TS allows a measurable reduction of the sizes of BOI and HP. The CS are not adopted in the optimal solution but, due to the arrangement of the optimal operation caused by the availability of the TS, the ABSs size decreases while the CCs size has to increase, in order to satisfy the cooling peak demand. Focusing on the environmental optimization, it should be observed that the adoption of large components does not imply any penalty for the environmental objective function. However, it can be outlined that the ST panels replace the PV panels adopted in the economic optimization. Table 7 shows the results of the optimizations performed for the isolated solutions. Comparing the economic optimizations to the conventional solution, the cost for natural gas used by BOIs significantly decreases, as well as that of the bought electricity. The cost for natural gas used by CHP did not exist for the conventional solution. However, the operating cost in the optimal isolated solutions is halved compared with the conventional one. The optimal isolated solutions, without and with the TS and CS, obtained minimizing the economic objective functions allow us to reduce the total annual cost by 36.7% and 38.8%, and the total annual emissions by 15.9% and 16.5%, respectively. Therefore, the adoption of TSs allows us to reduce either the total annual cost or the total annual emissions. The environmental optimization shows an increase of the operation costs of about 60% compared with the economic optimizations, while it allows a reduction of about 25% of the total annual emissions with respect to the conventional solution. It can be also noted that the amount of emissions due to electricity usage significantly increases, while the emissions saved due to the electricity sold to the grid, and that in consequence of natural gas usage, significantly decrease.
The energy balances of the optimizations performed for the isolated solutions are reported in Table 8. It can be observed that the bought electricity is negligible in the economic optimizations, while a significant amount of electricity is sold to the grid. The latter is higher when the TSs are adopted, as they allow us to decouple the thermal demand from the electric one, and to operate the ICE when it is more convenient. Almost all thermal demand is satisfied by the ICEs, while the cooling demand is covered by all three kinds of components which can be adopted (CC, ABS, HP). In the environmental optimizations, the quantity of electricity produced by CHPs and sold to the grid decreases noticeably, while the quantity of electricity bought from the grid increases. The heat produced by the HPs increases significantly, while the cooling energy required by the users is produced by CCs and HPs.

Distributed Generation Solution
The distributed generation solution is obtained by adding the DHN to the isolated solution ( Figure 7). The thermal energy produced in a production unit can be used directly either by the user of the site, or by other users, exchanging the heat through the DHN. This solution gives the opportunity to create a central production unit where all energy is produced and sent to the users through the DHN. The electric energy, if not used directly to the user has to be sent to the electric network. However, the isolated solution would still be adopted if it were more convenient, based on economic, or environmental evaluations. Table 9 shows the optimal configurations obtained for the distributed cogeneration solution. As previously introduced, the Pareto front of the bi-objective optimization has been identified by means of the ε-constraint method, by fixing the pollution levels at the 30%, 60% and 90% of the difference between the two environmental objectives of the economic and environmental optimal solutions. of the CCs is quite constant, while ABSs and HPs have opposite trends, moving from economic towards environmental optimum. The HPs result more suitable in the environmental optimal solution. In all optimal solutions, the thermal and cooling wasted energy is negligible.   From a general overview of Table 9, it can be noted that the DHN is convenient either from an economic or environmental point of view, while the MGTs, the CSs and the BOIs are actually never adopted, with the exception of the Environmental optimum.
Moving from the economic towards the environmental optimal solution, the number of the DHN pipelines increase, as well as the size of the HPs and the ST panels installed. The size of the CCs remains quite constant while the size of ICEs, of ABSs, of PV panels and of TSs adopted decreases. The sizes of the components related to the environmental optimal solution do not follow this trend, because the size of the components is free of penalty. Figure 8 shows the Pareto front of the distributed generation solutions, and is compared with the isolated solutions and conventional solution. It can be clearly seen that, with the exception of the environmental optimal solutions, all the other optimal solutions dominate the conventional solution and allow both the objective functions to be improve. However, the distributed generation solution does not lead to a significant improvement of the economic objective function, compared with the Isolated Solutions (less than 1%).  If two (or more) objectives are considered, the choice of the proper configuration has to be based on further technical evaluation of the designer and on an economic evaluation of the stakeholders. If the only aim is to achieve the lowest total annual cost, the economic optimal solution will be adopted, but in a different context, where the sensitivity to environmental problems is more important, the 90% environmental solution could be chosen. It allows the total annual emissions to be significantly reduced, controlling, at the same time, the Total annual cost.
In this specific case study, the 60% environmental optimal solution has been identified as the proper compromise. Figure 9 shows the optimal configuration and the lay out of the network in this case. This solution provides a subdivision of the DHN into two different subsystems: the first one is made up of sites 1-6, the second one is made of the sites 7-9. In the first subsystem the site 2 can be identified as a central node, where the greatest part of the thermal energy is produced and sent to the other sites through the DHN. In the second subsystem the sites are integrated with each other and a main site cannot be identified.  Table 10 shows the economic and environmental results of the optimizations. Moving from the economic to the environmental optimum the operation cost increase by about 80%, resulting from a significant increase of the bought electricity cost and a reduction of the natural gas cost and sold electricity income. The maintenance cost decreases noticeably, while the investment cost increases significantly. The environmental optimal solutions are characterized by a noticeable reduction of the total annual emissions, but a raise in total annual cost.  Table 11 shows the total optimal annual energy magnitudes for the different solutions. Moving from the economic towards the environmental optimal solution, it can be noted that the electricity produced by the ICEs decreases significantly. On the other hand, the electricity bought from the grid increases to feed electric components, such as CCs and HPs. In the economic optimal solution, the thermal energy produced by ICEs satisfies almost completely the demand, while it decreases significantly in the environmental optimal solution and is partially replaced by the thermal production of HPs and ST panels. As far as the cooling energy balance is concerned, the production of the CCs is quite constant, while ABSs and HPs have opposite trends, moving from economic towards environmental optimum. The HPs result more suitable in the environmental optimal solution. In all optimal solutions, the thermal and cooling wasted energy is negligible. If two (or more) objectives are considered, the choice of the proper configuration has to be based on further technical evaluation of the designer and on an economic evaluation of the stakeholders. If the only aim is to achieve the lowest total annual cost, the economic optimal solution will be adopted, but in a different context, where the sensitivity to environmental problems is more important, the 90% environmental solution could be chosen. It allows the total annual emissions to be significantly reduced, controlling, at the same time, the Total annual cost.
In this specific case study, the 60% environmental optimal solution has been identified as the proper compromise. Figure 9 shows the optimal configuration and the lay out of the network in this case. This solution provides a subdivision of the DHN into two different subsystems: the first one is made up of sites 1-6, the second one is made of the sites 7-9. In the first subsystem the site 2 can be identified as a central node, where the greatest part of the thermal energy is produced and sent to the other sites through the DHN. In the second subsystem the sites are integrated with each other and a main site cannot be identified.
Distributed generation solutions allow a slight improvement of the economic objective function (less than 1%) to be achieved, jointly with an improvement of the environmental one of about 2%, compared with the isolated solutions.

Distributed Generation Solution Integrated with the Central Solar System
The following results have been obtained by considering a superstructure which embeds all components included in the previous one, but a central system has been added, made up by an ICE, a BOI, the central storage and the solar field ( Figure 10).
Four optimizations have been performed: one economic optimization, one environmental optimization and two intermediate optimizations, obtained constraining the environmental objective function. Looking at the optimal configurations (Table 12) it can be noted that in all optimizations the solar field is adopted together with the central thermal storage, while the central BOI is not adopted because it is not economic to produce energy in the central unit and then transfer it to the sites, losing thermal energy through the DHN. Neither is the ICE is adopted, because the marginal cost of the heat results to be much higher when the electricity is sold to the grid (and therefore it has a low value) compared with the situation when the electricity is directly used by the users. Therefore, as the electric energy produced by the central ICE cannot be sent directly to the users, but can only be sold to the grid, the adoption of the central ICE is not economic. seasonal charging/discharging cycle: from April to August, part of the heat produced by the solar field is stored in the thermal storage, and is used from September to November. In the other months the storage is operated with weekly charging/discharging cycles (see Figure 13).
Larger thermal storage and solar field would probably allow a real seasonal operation of the storage: the heat stored from April to August would be used from September to March. However, this solution would have led to higher investment costs which would not have been paid back by a sufficient reduction of the operation costs.   A difference to be highlighted concerning the previous section is the presence of the boilers to cover the thermal peaks, the size of which increases as we move towards the environmental solution.
A separation in the DHN is still present in the economic and 30% environmental optimizations, forming two separated subsystems: sites 1-6, and sites 7-9, with the second group only connected to the central unit. Starting from the 70% environmental optimization, all users are connected to each other, whilst the storages installed in the local production units disappear, replaced by a larger central thermal storage, with lower heat loss. In parallel, the solar thermal field increases threefold compared with the economic optimum. The size of ABSs increases, too, while the ST panels at the user site are never adopted. Table 13 shows the economic and environmental results of the optimizations. The trends of the costs and of the emissions are similar to the ones observed in the previous optimizations. Comparing Table 10 with Table 13, it can be observed that the operating cost of the solutions integrated with the central solar system is lower (−6%), while they are characterized by higher investment costs (+30%). It is worth noting that the Pareto front obtained with integrating the central solar system (Figure 11) dominates the other Pareto fronts, as all the solutions obtained allow lower annual costs together with lower annual emissions to be achieved.  Table 14 shows the total annual energy magnitudes obtained for the optimizations of the distributed generation solution integrated with the central solar system. The trends of each item are similar to the ones observed in previous Table 11. The most important differences are the noticeable increase of the thermal production from the ST panels, and the greater amount of wasted thermal energy. They both are due to the presence of a large ST field, the production of which allows a great fraction of the thermal demand to be satisfied by renewable energy, but cannot be usefully exploited during all the year, even adopting the optimal operation strategy. In the pure and 30% economic optimizations, the solar field covers about 50% of the thermal demand and a small central thermal storage is adopted. It is operated with a daily/weekly charging/discharging cycles.  During summer, when the solar field is producing a lot of thermal energy, the thermal storage is full, but the users do not require thermal energy, therefore, the heat produced by the solar field is wasted. If a larger thermal storage were adopted, it would imply an investment cost noticeably greater, which would not be paid back by the related savings. Hence, it is cheaper to waste the exceeding heat (instead of storing it) and produce it with boilers and/or cogenerators, when it is necessary. Figure 12 shows the trends of the thermal energy demand of sites 7, 8 and 9, the central storage level and the in/out storage thermal flow in a typical week, operated following the economic optimal strategy. The other sites are not considered because the central storage is not connected to them. It can be noted that, during the weekend, the thermal demand is lower than during the working days, consequently the heat produced by the solar field is saved in the thermal storage. This heat is then used in the first hours of the following working days, when the operation of the ICE is not economic.   Figure 13 shows the optimal operation of the storage, with charging/discharging cycles, for the 70% environmental optimization of the distributed generation solution integrated with the central solar system. The figure shows clearly that each month is characterized by weeks with the same operation, consistently with the time discretization adopted in the optimization model. In fact, the In the pure environmental and 70% optimizations, the solar field and the thermal storage are larger, compared with the other optimizations. In this case, the thermal storage is operated with a seasonal charging/discharging cycle: from April to August, part of the heat produced by the solar field is stored in the thermal storage, and is used from September to November. In the other months the storage is operated with weekly charging/discharging cycles (see Figure 13).
Larger thermal storage and solar field would probably allow a real seasonal operation of the storage: the heat stored from April to August would be used from September to March. However, this solution would have led to higher investment costs which would not have been paid back by a sufficient reduction of the operation costs. Figure 13 shows the optimal operation of the storage, with charging/discharging cycles, for the 70% environmental optimization of the distributed generation solution integrated with the central solar system. The figure shows clearly that each month is characterized by weeks with the same operation, consistently with the time discretization adopted in the optimization model. In fact, the year is made up of 48 weeks (instead of 54) grouped by four, representing the 12 months. Each single week is made up of five working days plus two non-working days. Therefore, all weeks of the same month are characterized by the same operation. Additionally, in this case, a compromise between the best economic result and the most environmental friendly solution have been looked for on the Pareto front, and it has been identified with the 70% environmental optimal solution. It allows a total cost reduction by about 32% and a total annual emissions reduction by about 41%, compared with the conventional solution. Figure 14 shows the optimal configuration and the lay-out of the network. It can be noted that the DHN connects the users to each other differently from the optimal solution shown in Figure 9, where the users were subdivided into two sub networks. In this case the thermal energy produced in the central unit can be sent to all the users through the DHN. The dimension of the pipes decreases moving towards sites 1 and 7. The introduction of the central solar field and of a large thermal storage, allows the solution shown in Figure 14 to obtain a noticeable reduction of the CO2 emissions (−22%), compared to that shown in Figure 9.  Additionally, in this case, a compromise between the best economic result and the most environmental friendly solution have been looked for on the Pareto front, and it has been identified with the 70% environmental optimal solution. It allows a total cost reduction by about 32% and a total annual emissions reduction by about 41%, compared with the conventional solution. Figure 14 shows the optimal configuration and the lay-out of the network. It can be noted that the DHN connects the users to each other differently from the optimal solution shown in Figure 9, where the users were subdivided into two sub networks. In this case the thermal energy produced in the central unit can be sent to all the users through the DHN. The dimension of the pipes decreases moving towards sites 1 and 7. The introduction of the central solar field and of a large thermal storage, allows the solution shown in Figure 14 to obtain a noticeable reduction of the CO 2 emissions (−22%), compared to that shown in Figure 9.

Complete Distributed Generation Solution
The complete distributed generation solutions include also the DCN into the superstructure and corresponds to the most general superstructure presented in Section 2 (see Figure 1). Table 15 shows the optimal configuration of the four optimizations performed: economic,

Complete Distributed Generation Solution
The complete distributed generation solutions include also the DCN into the superstructure and corresponds to the most general superstructure presented in Section 2 (see Figure 1). Table 15 shows the optimal configuration of the four optimizations performed: economic, environmental and the two intermediate optimizations, obtained by constraining the environmental objective function. The optimal configurations and the trends are very similar to the ones obtained in the previous paragraph, where the DCN was not included in the superstructure. Also in this case, the first two optimizations (economic and 30% optimizations) subdivide the whole system in two sub-DHN: sites 1-6, and sites 7-9 connected to the central unit.
The last two optimizations (pure and 70%environmental optimizations) provide a DHN which connects all sites to each other. The DCN is adopted and always connects site 2 to sites 5 and 6, and site 9 to 7. Starting from the 70% environmental optimization, also site 1 is connected to site 6. Compared with the previous optimizations (see Table 12) the size of the CC decreases noticeable, because of the presence of the DCN.
The introduction of DCN does not affect significantly the economic and environmental performance of the system. In fact, an improvement of only a few percentage points is obtained for the annual cost and for the annual emissions (see Table 16), in comparison with the case considered in the previous paragraph. Focusing on the economic optimization, the operation costs decrease by about 60 k€. The annual investment cost increases by 10 k€, due to the adoption of the three pipes of the DCN. Thus, the total annual cost decreases by about 50 k€. The environmental optimization leads to a slight improvement of the total annual emissions (10 tons).
From Figure 15, it can be easily inferred that the complete distributed generation solutions dominates all the other solutions analysed for this case study. The economic optimal solution allows us to obtain a reduction by about 43% of the annual cost and by about 25% of the annual emissions (see Table 16), compared with the conventional solution. The environmental optimal solution allows a reduction of about 50% to be obtained, compare with the conventional solution, but the total annual cost increases. The corresponding performance of the 70% environmental solution, which has been identified as the best compromise among the complete distributed solutions, are a reduction of about 32% for the annual cost and of about 41% for the annual emissions.  A comparison among the annual energy magnitudes of the economic optimizations for the complete distributed generation and for the distributed generation integrated with the central solar system (Tables 14 and 17) shows some differences, which highlight the effect of the DCN in the economic optimal solutions. The electricity used by the HPs and the cooling energy produced by the ABSs are higher by about 67% and 8%, respectively, when the DCN is adopted. At the same time, the sold electricity and the electrical consumption of the CC are lower by about 17% and 68%, respectively. Nevertheless, the additional operation options allowed by the adoption of the DCN do not lead us to a strong improvement of the environmental, or of the economic performance of the system, but only small improvements by about 1% and 3%, respectively, are obtained.
The evaluation of Pareto fronts of the multi-objective optimizations performed allows the best compromise between the minimum annual cost and the minimum annual emissions to be identified for each different heat/cooling integration options. Table 18 summarizes these compromise solutions and compares them with the economic optimum obtained for the conventional and isolated solutions.
As mentioned before, the 70% environmental optimization has been identified as the best  A comparison among the annual energy magnitudes of the economic optimizations for the complete distributed generation and for the distributed generation integrated with the central solar system (Tables 14 and 17) shows some differences, which highlight the effect of the DCN in the economic optimal solutions. The electricity used by the HPs and the cooling energy produced by the ABSs are higher by about 67% and 8%, respectively, when the DCN is adopted. At the same time, the sold electricity and the electrical consumption of the CC are lower by about 17% and 68%, respectively. Nevertheless, the additional operation options allowed by the adoption of the DCN do not lead us to a As mentioned before, the 70% environmental optimization has been identified as the best compromise for the Complete distributed solution. Figure 16 shows its optimal configuration: the layout of the DHCN is shown together with the size of the components installed in each site. The DCN connects sites 2, 5, 6 and 1, and a single pipe connects site 7 to 9. The layout of the DHN is very similar to the one shown in Figure 14 for the distributed generation solution integrated with the central solar system, while the component sizes installed in each production unit are slightly changed. These differences are consequence of the additional operation options allowed by the adoption of the DCN.

Conclusions
The MILP model presented in this paper has been used to optimize the configuration and operation for a case study of a distributed generation system for heating and cooling in an urban area. The model allows us to obtain the optimal configuration of the system, satisfying the energy requirements of the users, that minimizes the total annual cost and the total CO2 emissions during operation. The aim of the work is to compare on a common basis, different district heating integration options. In more detail, five cases have been considered: • conventional solution; Nevertheless, in the evaluation of the optimization results, we have to keep in mind that the MILP problems have been optimized with a gap equal to 1%. Therefore, two feasible solutions with a very low difference in the objective functions (in principle, lower than 1%) may be regarded by the optimization algorithm as equivalent, even if their decision variable vectors are different. When the DHCN is included in the superstructure, the number of similar solutions increases exponentially, and solutions which have a different location, or size, of some components may result to be equivalent for the MILP algorithm, because of their similar optimal values of the objective functions.

Conclusions
The MILP model presented in this paper has been used to optimize the configuration and operation for a case study of a distributed generation system for heating and cooling in an urban area. The model allows us to obtain the optimal configuration of the system, satisfying the energy requirements of the users, that minimizes the total annual cost and the total CO 2 emissions during operation. The aim of the work is to compare on a common basis, different district heating integration options. In more detail, five cases have been considered: Starting from the conventional solution, the further cases consider, at each step, additional components and a more thorough district heating/cooling integration. The best solution, both from the economic or environmental points of view can be achieved with the complete distributed generation solution, which includes various kinds of energy components, a DHCN, a solar field and a seasonal thermal storage. The best economic solution allows a 43% reduction of the total annual cost compared with the conventional solution, while the best environmental solution allows about a 50% reduction of the total annual emissions (see Table 16). This result can be achieved only with the optimal size and/or multiplicity of the components and with the optimal operation of the energy system.
The economic optimization of the Isolated Solution, which comprises only the distributed cogenerators and absorption machines, without any district integration, permits a consistent reduction of the total cost compared with the conventional solution (37%), together with a reduction of the CO 2 emissions (16%). In this case, the adoption of local thermal storages is suggested from both the economic and environmental points of view). This performance is similar to the one obtained for the economic optimization of the complete distributed generation; the isolated solution also implies a much lower investment cost (Table 7). On the other hand, the complete distributed generation solution is cheaper in the long term view (over 20 years) and allows also a greater reduction of the annual emissions. In fact, important reductions of the annual emissions can be obtained only with the adoption of the solar field, of the seasonal storage and of the district heating network.
The distributed generation solution, with heating district integration, but without the solar field and the seasonal storage, permits us to slightly improve both the economic and the environmental objective functions (less than 2%) compared with the isolated solution. The further integrations of the DG system allow incremental benefit of similar magnitude for the economic optimal solutions (1-3% each), reaching the higher reduction with the complete distribution generation solution, as highlighted before.
From the analysis of the Pareto fronts obtained, it can be inferred that a very high environmental benefit (greater than 40%) can be obtained by complete distribution generation solution, without giving up a large long-term cost reduction of about 32% (see Table 18). At the same time, it is evident that the advantages obtained by introducing the cooling district network are very small, at least when the optimal synthesis and operation of the system is considered.
When the central solar field with central TS can be chosen, the optimal compromise solution adopts a central TS of about 172,000 kWh and a solar field of about 22,000 m 2 . The optimal operations identified in this paper show that the central thermal storage is operated with seasonal charging/discharging cycles only when the environmental objective function is considered in the optimizations (Figure 13). The heat produced by the solar field during warmer months is used during the first colder months.
Meanwhile, pure economic optimizations provide a weekly operation of the central thermal storage: the heat produced by the solar field during week-end, when the energy demand is lower, is used during the following working days.