Integration of Micro-Cogeneration Units and Electric Storages into a Micro-Scale Residential Solar District Heating System Operating with a Seasonal Thermal Storage

: A micro-scale district heating network based on the operation of solar thermal collectors coupled to a long-term borehole thermal storage is modeled, simulated and investigated over a period of ﬁve years. The plant is devoted to covering the domestic hot water and space heating demands of a district composed of six typical residential buildings located in Naples (southern Italy). Three alternative natural gas-fueled back-up auxiliary systems (condensing boiler and two di ﬀ erent technologies of micro-cogeneration) aiming at balancing the solar energy intermittency are investigated. The utilization of electric storages in combination with the cogeneration systems is also considered with the aim of improving the self-consumption of cogenerated electric energy; heat recovery from the distribution circuit is also evaluated to pre-heat the mains water for domestic hot water production. The performances of the proposed plant schemes are contrasted with those of a typical Italian decentralized heating plant (based on the utilization of natural gas-fueled non-condensing boilers). The comparison highlighted that the proposed conﬁgurations can decrease the primary energy consumption (up to 11.3%), the equivalent emissions of carbon dioxide (up to 11.3%), and the operation costs (up to 14.3%), together with an acceptable simple pay-back period (about 4.4 years).


Introduction
The building sector is responsible of a significant portion of the world's greenhouse gas emissions because it consumes around 40% of the overall energy production; a large amount of this energy consumption is used to cover the heating demand [1].
District heating (DH) networks, integrating renewable energy sources, are characterized by several benefits with respect to decentralized heating plants [2,3]. In particular, during the last years solar district heating systems achieved a great attention in the scientific community, with a significant number of worldwide applications [4][5][6]; this is because solar energy is easily and economically accessible during all year in most part of the world and its exploitation is supported by several economic incentives. However, one of the hurdles to solar energy utilization is related to the mismatch between the availability of solar energy and thermal demand; this is one of the reasons why the concept of thermal energy storage is considered as one of the best options to address this temporal misalignment [1,7]. A seasonal thermal energy storage allows to store thermal energy over long periods (weeks or months); according to the review of Rad and Fung [8], borehole thermal energy storage (BTES) is characterized by several advantageous conditions for energy storage over long periods [8][9][10][11][12] with respect to other technologies. BTESs consist of closed-loops where heat is charged or discharged by vertical or horizontal Borehole Heat Exchangers (BHEs) which are installed into boreholes below the ground surface. After drilling, a "U" pipe is inserted into the borehole; the borehole is then filled with a high thermal conductivity grouting material. BHEs can be single or double U-pipes and they can be hydraulically connected in series or in parallel. The underground is used as storage material in a BTES.
However, intermittent energy sources with stochastic behavior have to be coupled and compensated by other generation devices. A number of technologies could be used in solar DH networks as back-up units with the aim of addressing the intermittency of solar energy. Micro-cogeneration (MCHP) units are systems simultaneously producing thermal and electric power with an electric generation below 50 kW el ; they are considered as one of the most efficient solutions to reduce primary energy consumption as well as greenhouse gas emissions [13][14][15][16]. Micro-cogeneration devices applied in residential units generally operate under a heat-led control logic in order to maintain a desired temperature level [17] inside a thermal energy storage. In this case, a large fraction of the cogenerated electricity is usually produced at times when the electric demand of buildings is lower than MCHP electric output (due to typical mismatch between thermal and electric demands [14]). Therefore, substantial amounts of electricity have to be exported to the grid with reduced revenues in the case of heat-led control logic. As a consequence, the utilization of distributed electrochemical battery energy storage systems could be a way to increase the own use of cogenerated electricity, enhancing the profitability of MCHP devices. Batteries currently represent mature devices with high energy densities [15,18] and they could act as a key technology for load-leveling and peak-shaving, damping energy oscillations as well as improving power reliability and quality [19][20][21]. Gomes et al. [22] highlighted that Li-ion batteries are nowadays the major technology applied in electric vehicles and represent a very good candidate to stationary applications for a number of reasons, such as high energy and power per unit as well as they are lighter and smaller than other types. They can be used with a depth-of-discharge up to 80% of their total capacity, with a great number of cycles varying between 2000 and 5000 according to some manufacturers and greater efficiency with great loads as their power inverters have almost no Peukert's loss. High battery costs have been the greatest barrier to the massive adoption of battery storage systems in the last few years. However, according to [22], since 2008, Li-ion battery costs have been reduced by a factor of four and their costs are projected to decrease to around 100 $/kWh el by 2030 benefiting from important economies of scale driven by the growth in the electro-mobility sector. Several scientific studies expect the economic viability of residential hybrid Li-ion battery storages/solar photovoltaic system in the near future [23]. Naumann et al. [24] used a multi-parameter economic model for estimating the profitability of 'home storage' battery usage for buffering the surplus of photovoltaic (PV) generation upon varying both technical and economical parameters; they found one viable scenario for 2018, with a strong decrease in battery system prices below 450 €/kWh el and an increase in retail electricity prices. Linssen et al. [25] performed an optimization analysis from a techno-economic perspective of photovoltaic battery systems; according to their simulation results, a photovoltaic battery system can be economically operated with a battery cost of 1200 €/kWh el if subsidies are taken into consideration. In contrast, Moshövel et al. [26] developed a tool which is able to calculate a proper battery size for a wide range of different combinations of households and PV panels; they highlighted a possibility of economical operation at battery prices below 500 €/kWh el , and only partly economical operation at prices between 500 ÷ 700 €/kWh el . Parra and Patel [27] investigated the performance of lead-acid and Li-ion battery systems in combination with PV generation for a single home in Switzerland using a time-dependent analysis; they found that a storage medium cost equal to 375 CHF/kWh el , a durability of 5000 equivalent full cycles or a 6.5% retail price increase per year (based on a retail electricity price equal to 0.22 CHF/kWh el in 2015) are necessary to achieve positive internal rate of return in Geneva. Pawel [28] developed a new framework for the calculation of levelized cost of stored energy for a PV and storage combined power plant; he concluded that economic viability is given at battery prices lower than 300 €/kWh el . Vonsien and Madlener [29] evaluated the economic efficiency of a Li-ion battery storage unit attached to a photovoltaic system in a private household. They found that, regarding the battery prices, a cost of 700 €/kWh el could be identified to be economically efficient for a battery lifetime of 21.6 years and, according to their estimations, a price of 700 €/kWh el will be reached in 2023 by assuming an average decrease of the costs for stationary Li-ion batteries. Dufo-López and Bernal-Agustín [30] presented a methodology to evaluate the technical and economic performance of a grid-connected system with storage under a time-of-use electricity tariff; the methodology was applied for lead-acid or Li-ion batteries. The authors focused on the storage system's profitability for the electricity consumer, analyzing the total net present cost of a system with storage and comparing it with a system without storage; the results showed that the storage system is economically profitable in the case of cost of battery storage (battery bank + bidirectional inverter + control) is 0.075 €/kWh el in the case of efficient Li-ion batteries (with bi-directional converter of 300 €/kW el ). Two deployment cases of Li-ion batteries (for frequency regulation in the eastern United States and electricity bill reduction for commercial or industrial customers by limiting their peak load in California) have been defined and examined in cost-benefit and sensitivity analyses by Günter and Marinopoulos [21]; the results showed that energy storage is cost-efficient in these cases even if frequency regulation market prices and subsidies drop below 2016 s level. According to the findings of the conducted analyses, the authors also concluded that electric energy storage offers valuable alternatives to other grid resources and can integrate fluctuating generation and demand, if market structures that recognize the value of energy storage and the services it can provide are in place. Martinez-Bolanos [31] assessed the economic feasibility of replacing conventional peak power plants, such as diesel generators, by using distributed electric energy storages, to implement energy time shift during peak hours for commercial consumers; the economic analysis was performed by calculating the break-even point for four different battery technologies using the HOMER Energy software in the simulations. The results showed that a reduction, with respect to the levels of 2018, of about 30% in the costs of battery makes the installation of batteries viable; in addition, they estimated that, in a 5-to 6-year horizon, these technologies would become economically attractive as a result of the strong decrease of the expected costs for the years to come.
One more concern is related to the fact that it is predicted that lithium required for the production of Li-ion electric storages will face a severe shortage in the foreseeable future [32]; in addition, it should be considered that the emissions generated in the procurement and production stage of supply chain represent one of the main contributors through electric vehicle manufacturing [32]. Vandepaer et al. [33] also showed how Li-ion batteries cause significantly more impacts than lithium metal polymer units in terms of global warming and ozone depletion, mainly during the battery manufacturing stage; they also found that centralized battery system configurations bring smaller environmental effects than distributed systems with more and smaller storage units. However, according to the scientific literature [34,35], the feasibility of electric storages coupled with electric generation devices is still an open question and needs to be further investigated taking into account several factors with high impact on the economics and the environment. Gomes et al. [22] indicated that the economic assessment of DC microgrids connected with DC generators and DC loads (stationary batteries) compared to AC microgrids is lacking and additional studies have to be performed. In particular, the feasibility of batteries depends on several factors, such as electric demand profile, electric generation profile, self-consumption of generated electric energy, unit cost of electric energy purchased from the central grid, unit price of electric energy sold to the central grid, capacity and technology of electric storages and electric generation units, economic incentives, etc.. Some scientific works analyzed the operation of cogeneration devices while coupled with batteries [36][37][38]. Darcovich et al. [36] assessed the energy consumption of a Canadian house served by a Stirling engine-based micro-cogeneration device operating in combination with a 2 kW/6 kWh battery. Gimelli et al. [37] suggested a methodology to select the optimal configuration of a modular combined heat and power plant equipped with an electric storage device. Darcovich et al. [38] also analyzed the operation of a plant consisting of an internal combustion engine-based MCHP system, photovoltaic panels (PV) and electric battery; the authors highlighted that the utilization of PV systems and battery could generate economic advantages, even if to a smaller degree in comparison to the MCHP device.
To the knowledge of the authors, only five works [39][40][41][42][43] studied the performances of Italian district heating networks operating with seasonal thermal storages combined with solar thermal collectors. In 1998 Oliveti et al. [40] analyzed the annual performance of a plant located in Reggio Calabria (southern Italy) substantially composed of a 500 m 3 long-term storage coupled with 91.2 m 2 of solar thermal collectors dedicated to covering the heating demands (111 GJ) of an office (floor area of 1750 m 2 ); they found a value of solar fraction equal to about 28%. Buoro et al. [41] investigated the combination of distributed energy supply units, solar systems and seasonal thermal storage, integrated into a district heating plant; the proposed scheme was designed to cover the cooling, thermal, and electric requirements of 9 industrial buildings situated in the north of Italy. The data demonstrated that the utilization of solar systems is potentially able to reduce by 5% the annual costs as well as by 15% the consumption of primary energy in comparison to a conventional scheme; moreover, they found that the solar field should be sized in order to satisfy about 55-60% of the annual thermal requirements. In [42] the authors analyzed, through the software TRNSYS [44], the performance of a DH network (including 200 m 2 solar thermal collectors combined with a BTES) devoted to satisfying the energy demand of a school situated in Palermo (south of Italy). The simulation data highlighted a significant energy saving in contrast with conventional plants, confirming the suitability of BTES while applied under Italian climatic conditions. In 1982 [43], a DH network based on the operation of a 2250 m 3 seasonal borehole storage coupled with 180 m 2 of solar thermal collectors was analyzed form an experimental point of view while devoted to satisfy the demand for domestic hot water (DHW) production and space heating of 280 MWh/year (north of Italy); the analysis of system performance highlighted a significant solar fraction equal to about 80%. In 1985 [39,43], the operation of a solar DH plant was experimentally analyzed while serving residential buildings (situated in northern Italy) characterized by a total floor area of about 9200 m 2 with an energy demand for DHW and space heating equal to about 3456 GJ. The plant was based on the operation of a long-term storage with a volume of 43,000 m 3 combined with 2727 m 2 of solar thermal collectors; the authors reported a relevant solar fraction equal to about 72%.
The analysis of scientific literature indicates that the studies carried out with reference to Italian applications are limited and generally dated, with just two works [40,42] performed under climatic conditions associated to southern Italy; this represents a relevant research gap taking into account that it is very well known that climatic conditions significantly influence both energy requirements as well as solar energy recovery and, consequently, the performance of plants. Moreover, most of these studies are based on an energy analysis only, generally neglecting the economic as well as the environmental aspects. In addition, it should be underlined that the size of the analyzed districts is significantly bigger than that one investigated in this work. Finally, the published papers are not focused on systems integrating cogeneration units and/or electric energy storages. In conclusion, the literature review highlights a significant lack of studies focused on micro-scale solar hybrid district heating networks including borehole thermal energy storages, micro-cogeneration devices and batteries; this confirms the necessity to carry out new analyses for Italian applications aiming to evaluate their potential benefits and drawbacks and, therefore, the novelty of the present study.
In this work, a DH network based on the operation of solar thermal collectors connected to a long-term borehole thermal storage is modeled, simulated and investigated for a period of 5 years. The plant is used to covering the domestic hot water and space heating requirements of a micro-scale residential district composed of six typical homes located in Naples (southern Italy). Three alternative natural gas-fueled auxiliary back-up systems ((i) a 26.1 kW condensing boiler, (ii) two parallel-connected 12.5 kW th Internal Combustion Engine (ICE)-based micro-cogeneration units and (iii) a 26.0 kW th Stirling Engine (SE)-based MCHP device) are considered for balancing the solar energy intermittency; the utilization of batteries in combination with the cogeneration systems is also considered with the aim of improving the self-consumption of cogenerated electric energy; heat recovery from the distribution circuit is adopted to pre-heat the mains water for DHW production by using decentralized tanks. The performances of the proposed schemes are compared with those of a typical Italian decentralized heating plant (based on the utilization of natural gas-fueled boilers) in terms of consumption of primary energy, equivalent emissions of carbon dioxide, operation and capital costs.
This paper mainly aims to assess: • the feasibility of micro-scale solar district heating networks in comparison to conventional heating systems in the case of Italian scenario; • the performance of long-term thermal energy storages in the case of micro-scale solar DH networks; • the most suitable back-up technology to be adopted into micro-scale solar DH networks; • the feasibility of micro-cogeneration units in satisfying the requirement of electricity associated to small Italian district composed of residential buildings only; • the capacity of batteries in improving the self-consumption of cogenerated electric energy and, consequently, decreasing the amount of power bought from the central electric grid; • the benefits deriving from the possibility to pre-heat the mains water for DHW production by recovering heat from the distribution circuit through decentralized storages.
In this paper, the residential district served by the proposed plants is described in the Section 2 in terms of characteristics of buildings' envelopes as well as thermal and electric loads. The Section 3 depicts the investigated configurations of the solar district heating network, the main characteristics of both plant components as well as selected dynamic simulation models, together with the strategies controlling the operation of the systems. Details about the typical Italian decentralized heating plant compared with the proposed solar district heating networks are reported in the Section 4. The Section 5 summarizes the methods used to perform the energy, environmental and economic comparisons, while the Section 6 shows and comments the results of the study.

Description of the Residential District
The proposed DH plants are analyzed while serving six typical Italian homes situated in Naples (latitude = 40 • 51 46" 80 North; longitude = 14 • 16 36" 12 East; Heating Degree Days = 1034). Three different types of residences are analyzed (building type A, building type B, building type C). In more detail, the community consists of 2 residences for each type; each type varies from the others with respect to the floor area, volume, area of windows, maximum number of occupants (as indicated in Table 1). To accurately evaluate the performances associated to the proposed schemes, the occupant-driven loads are taken into account.
For each type of building, a specific annual stochastic profile (consisting of 365 different daily profiles) characterized by a 1-min time resolution is determined by means of the models suggested in [45] with the aim of modeling the presence of occupants.
A specific annual stochastic profile (consisting of 365 different daily profiles) characterized by a 1-min time resolution is considered for each type of building in order to represent the electric demand (due to domestic appliances and lighting systems); these stochastic profiles are defined on the basis of the models of Richardson and Thomson [45].
The annual profiles representing the internal gains for each type of building are obtained by assuming the thermal gains associated to the lighting appliances equal to 70% of their nominal power consumption, the thermal gains of occupants equal to 70 W/person and the thermal gains due to the operation of domestic appliances according to the data suggested by the manufacturers.
Thermal transmittances of buildings are assumed equal to the values (0.38 W/m 2 K for external vertical walls, 0.40 W/m 2 K for floors, 0.36 W/m 2 K for roofs, 2.40 W/m 2 K for windows) imposed by the Italian legislation requirements with reference to the year 2018 [46], whatever the building type is.
The calculation of the air exchange rate is based on the European Standard EN 12831 [47], resulting in 0.24 volumes/hour for each house.
In this paper, stochastic DHW demand profiles in the time scale of 1 min are considered. In particular, for both the building types A and B, a DHW demand profile with an average load equal to 100 l/day is selected, while for the building type C a DHW demand profile with an average load equal to 200 l/day is adopted. The selected DHW profiles are derived from the sets of yearly profiles identified within the Task 26 of IEA-SHC [48], where the time of occurrence and the values of every flow rate is defined by statistical means. Table 2 reports the values of (i) annual space heating energy demand, (ii) annual energy demand for DHW production and (iii) annual electricity demand (due to domestic appliances, lighting systems, fan-coils, and individual pumps only) associated to each building as well as all end-users.
The schematics of the proposed configurations are reported in Figure 2a-e, where Figure 2a corresponds to the plant configuration named SCHEME 1, Figure 2b corresponds to both the plant configurations named SCHEMES 2 and 3, Figure 2c corresponds to the plant configuration named SCHEME 4, Figure 2d corresponds to both the plant configurations named SCHEMES 5 and 6, and Figure 2e corresponds to both the plant configurations named SCHEMES 7 and 8. In these figures, the following major components are indicated: end-users (6 houses), solar thermal collectors (SC), short-term vertical cylindrical thermal energy storage (STTES), seasonal single U-pipe borehole thermal energy storage (BTES), heat dissipator (HD), back-up system (BS), two counter-flow plate heat exchangers (HE1 and HE2), 6 local/decentralized individual non-condensing boilers (B) for DHW production (one per each building), 6 short-term DHW vertical cylindrical tanks (DHWT) integrated with an internal heat exchanger (IHE) to pre-heat the mains water (one per each building), fan-coils (FC), inverter/charge controller, electric batteries, pumps (P), 3-way valves and pipes.  Proposed CSHPSS configurations: schematic of the SCHEME 1 (a), schematic of the SCHEME 2 and the SCHEME 3 (b), schematic of the SCHEME 4 (c), schematic of the SCHEME 5 and the SCHEME 6 (d), schematic of the SCHEME 7 and the SCHEME 8 (e). The 8 proposed schemes exhibit some identical characteristics. In particular, whatever the CSHPSS configuration is: • the heat carrier fluid used in this paper is a mixture of ethylene glycol and water (40%/60% by volume) to achieve a boiling temperature of about 105 • C at 1 bar; • solar energy recovered by the solar collectors is firstly moved, thanks to the HE1, into the STTES; • excess solar energy is dissipated by blowing air through the HD (consisting of a finned coil heat exchanger) in the case of the temperature of the heat carrier fluid at the outlet of the solar collectors becomes larger than 95 • C; • in the case of space heating requirements, solar energy is moved from the STTES, by means of the HE2, into the fan-coils installed into the buildings via the distribution circuit to obtain the desired indoor temperature level during the heating period (15 November-31 March, according to [49]); • solar energy stored into the STTES can be moved into the BTES during the entire year ("BTES charging mode") in the case of it is not instantaneously requested for space heating. Solar energy from the BTES can go back to the STTES ("BTES discharging mode") only during the heating period with the aim of integrating the level of temperature inside the STTES. The direction of the heat carrier fluid is from the BTES center to the storage boundaries during the charging phase, with the aim of obtaining larger temperatures in the BTES center and smaller ones at the storage boundaries; the direction of the heat carrier fluid is inverted during the discharging phase; • when the solar energy stored into both the STTES and BTES is not able to fully cover the energy demands, a back-up auxiliary system is operated with the aim of supplementing the thermal energy needs in order to achieve the desired supply temperature level (55.0 • C); • the DHW is produced by means of 6 local individual boilers (B), one per house.
The proposed 8 configurations only differ in terms of: (a) back-up auxiliary system, and/or; (b) production of DHW, and/or; (c) utilization of electric energy storage.
The differences between the proposed schemes are specified and summarized in Table 3.  The target of supply temperature for the distribution network is set to 55.0 • C; therefore, the auxiliary unit (MB or SE-MCHP or ICE-MCHP) is operated only when integrating the temperature level at the outlet of the heat exchanger HE2 is required to get the desired value. Therefore, the MCHP units are operated under a thermal-load following control logic according to the temperature level at the outlet of the HE2, providing their maximum thermal and electric outputs while activated. The electric energy provided by the MCHP systems is used to meet the electric demand; any excess of the MCHP electric output is charged into the batteries (if used) or sold to the power line; the electric storages (if used) and the central grid are also utilized to satisfy all peaks of electric demand. Figure 2a-e and Table 3, with reference to the DHW production, it can be noticed that:

As indicated in
(a) the DHW is produced via 6 natural gas-fired non-condensing individual boilers only (one per house) for the SCHEMES 1-3, 7 and 8; in particular, the mains water enters the boilers and it is heated up to 45.0 • C; (b) for the remaining configurations (SCHEMES 4-6) the DHW is produced thanks to the operation of 6 natural gas-fired individual non-condensing boilers in combination with 6 local vertical cylindrical DHW tanks (one per building) equipped with one inlet, one outlet as well as one internal heat exchanger (IHE). In particular, during the heating period (15 November-31 March), mains water enters the heat exchanger immersed into the DHW tank, while the heat carrier fluid exiting the fan-coils is transferred into the DHWT with the aim of pre-heating the mains water. The individual boiler is activated with the aim of achieving the desired temperature of 45.0 • C only when the temperature of mains water exiting the heat exchanger is lower than the given target. During the remaining part of the year, the DHW at 45.0 • C is obtained through the 6 individual/decentralized boilers only, without using the DHWTs (as in the configurations SCHEMES 1-3).
As highlighted in Figure 2a-e and Table 3, three-series connected electric energy storages are used for the SCHEMES 5-8 (all including a micro-cogeneration system as back-up unit), while the batteries are not used for the remaining four configurations (SCHEMES 1-4). Whatever the season is (during both heating and cooling periods), the cogenerated electricity is firstly used to cover the electric load for the SCHEMES 5-8. In the case of the MCHP electric output is bigger than the electric demand and state of charge of the batteries is 100%, then the cogenerated electricity is sold to the central grid.
In the case of the charge level of the electric storages is lower than the maximum (100%), then the eventual surplus of electric energy generation with respect the electric demand is charged into the batteries. In the case of the cogenerated electric output is not enough to fully satisfy the electric needs, the electric storages are discharged only when their state of charge is larger than 10%; the electric storages discharging is interrupted in the case their state of charge becomes smaller than 10%. The grid is used to cover all peaks of electric demand in the case of the cogenerated electricity is not enough to fully cover the electric demand and the state of charge of the batteries is lower than 10%.
In the case of the configurations with the main boiler as back-up system (SCHEMES 1 and 4), the electric demands are satisfied by means of the electricity provided by the electric grid only. Table 4 describes and summarizes the principal characteristics of the leading components included in the proposed configurations reported in Figure 2a-e.

Simulation Models
Accurate numerical models are adopted to simulate the plant components in order to take into account (i) the occupant-driven loads and thermal behavior of buildings, (ii) the partial load operation of all components, (iii) the coupling between heating loads and simulation outputs of the components, and (iv) the logics controlling the operation of the plant.
The TRaNsient SYStems [44] environment is selected for modeling and analyzing the proposed plants during 5 years by setting a time step equal to 1 min. The length of the period of simulation (5 years) is selected with the aim of taking into account that it requires a significant time to totally charge the BTES; the first day of simulations is set to 1 January.
In the TRNSYS software, individual models (named "Types") are used to model each sub-system. In this study, the Types are chosen from the TRNSYS library and calibrated based on information provided by the manufactures or data available in the literature. In the following the TRNSYS Types selected in the project for running the simulations of the proposed systems are briefly described.

Battery [56]
Number of series-connected batteries 3 Single battery capacity (kWh) 13.5 Efficiency round-trip (%)/Depth of discharge (%) 90/100 Single battery power (kW) 5 (continuous)/7 (peak) Inverter and charge controller [44] Efficiency of regulator (%)/Inverter (%) 78 The Type 56 is selected for modeling the thermal behavior of the buildings of the district. The BTES is modeled by means of the Type 557a [57], i.e., the Duct ground heat STorage (DST) model, representing the state-of-the-art in modeling borehole heat exchangers (BHEs) [13]. In this model, the storage is characterized by a cylindrical shape with vertical symmetry axis. The layout of boreholes is fixed in a hexagonal shape and uniformly within the BTES field and the ground is assumed to be homogeneous. The top of the storage is covered with insulating material. The BTES numerical model is described in greater detail in [57]. Pärisch et al. [58] analyzed the performance of 3 double U-tube borehole heat exchangers with a depth of about 70 m by means of 2 h of experimental tests keeping constant the inlet temperature; the measured data have been compared to three different simulations models, including the DST model used in this paper. The results of the comparison between the measured outlet temperature and the outlet temperature predicted by the DST model highlighted how the predicted trend (i) follows the curve of measured inlet temperature, (ii) underestimates the thermal power injected in the borehole heat exchanger only during the first 0.25 h, (iii) is able to better approach the measured data in the remaining part of the simulation time. The authors selected the BTES volume (435.8 m 3 ), the number (8), depth (12.43 m) and connection scheme (in series) of the boreholes thanks to the results of a large parametric analysis carried out in a previous paper [9]. In particular, concerning the configuration of the BTES, the authors investigated [9] the impacts of (i) the volume of the BTES, (ii) the depth of boreholes, (iii) the borehole heat exchangers (BHEs) number, as well as (iv) the BHEs type of connection on the performances of the solar DH plant while satisfying the energy demands of the same residential buildings; the values allowing to maximize the primary energy savings have been adopted. The thermal conductivities of soil and grout are defined based on information available in the current literature; in particular, Casasso and Sethi [59] indicated that thermal conductivity of soil could range between 0.5 and 3.0 W/mK, while thermal conductivity of grouting material can vary in the range 0.5 ÷ 5.0 W/mK.
The Type 534, selected for modeling the STTES, allows to divide the storage into 10 fully mixed equal sub-volumes with 10 isothermal temperature layers (top layer = node 1 and bottom layer = node 10) with the aim of taking into account the stratification in the storage. The tank model is detailed according to the data suggested by the manufacturer [51]. The authors selected the STTES volume (6.0 m 3 ) thanks to the results of a large parametric analysis carried out in a previous paper [9].
The performances of DHWTs are also analyzed by using the Type 534; 10 isothermal temperature layers (where the bottom layer is 10 and the top layer is 1) are considered to take into account the stratification in the decentralized storages. Manufacturer specifications are taken into account to calibrate the numerical model [52].
The flat-plate solar collectors are modeled by means of the Type 1b. In more detail, the thermal efficiency of collectors is defined by the following equation: where G is the solar irradiation incident on collectors surface in kJ/hm 2 , T amb is the ambient temperature, T in,SC is the fluid inlet temperature. The values of 0.7484 and 16.17 hm 2 K/kJ are, respectively, the intercept efficiency and the efficiency slope. In the model, the intercept efficiency is adjusted for off-normal solar incidence by using the following incidence angle modifier K: where S = (1/cosθ − 1) with θ being the incident angle for beam radiation, b 0 = 0.103 and b 1 = 0. The coefficients adopted in Equations (1) and (2) are defined according to manufacturer specifications [50]. The orientation and tilted angle of solar collector are selected for optimizing the performance with reference to the city of Naples. The values of the thermal efficiency of a single solar thermal collector measured by the manufacturer [50] according to the Standard [60] have been compared by the authors against the simulation values obtained by using the model (Type 1b) used in this work under the same operating conditions. The comparison indicated that the percentage differences between the manufacturer data and the outputs predicted by the model are very small (between a minimum of −1.52% up to a maximum of 1.14%) and, therefore, confirmed the suitability of the adopted model. The authors selected the area of the solar field (55.44 m 2 ) thanks to the results of a large parametric analysis carried out in a previous paper [9]. The Type 511, representing a system for cooling a liquid by means of external air, is selected for modeling the heat dissipator.
HE1 and HE2 are simulated as counter-flow plate heat exchangers by using the Type 5b; an overall heat transfer coefficient, respectively, of 54.3 W/K per m 2 of collector area for the HE1 and 25 W/K per GJ of heat demand for the HE2, are set based on the criteria indicated in [61].
As mentioned in the previous section, 3 alternative different auxiliary systems are analyzed in this work: First of all, the capacity required to the auxiliary back-up system in order to maintain the level of temperature at the outlet of the heat exchanger HE2 up to 55 • C (i.e., desired target of supply temperature) has been estimated by performing dedicated simulations by means of the software TRNSYS (with a time step of 1 min). The simulation results highlighted that the thermal power to be supplied by the back-up system to obtain 55 • C as supply temperature is always lower than 53.4 kW, and, in particular, is lower than 25.0 kW during 98% of simulation time. Therefore, the capacities of back-up units have been selected by considering devices commercially available on the Italian market and taking into account this simulation result.
The performance of the condensing main boiler (used in the configurations named SCHEMES 1 and 4)) is simulated in this paper by means of the Type 700 by considering its thermal efficiency (i.e., the ratio between thermal output and primary input of the unit) depending on the load factor (i.e., the ratio between the current and the maximum thermal outputs of the system). In particular, the following equation is implemented in the Type 700 for estimating the thermal efficiency of the condensing boiler: where P CSHPSS th,out,MB is the thermal power supplied by the condensing boiler to meet the district needs and P CSHPSS th,nom,MB is the nominal thermal power of the condensing boiler (26.1 kW th ). The experimental thermal efficiency of the selected condensing boilers is measured by Di Perna et al. [62] and compared by the authors against the predictions obtained by using Equation (3) upon varying the load factor. The comparison highlighted that the percentage differences between the measured performance and the simulation outputs are in the range −2.37% ÷ 4.8%, confirming the suitability of the model (Type 700) adopted in this study. According to the information reported by Di Perna et al. [62], a minimum turn-down ratio (defined as the ratio between the minimum and maximum capacities of the system) equal to 0.21 is assumed.
The decentralized non-condensing boilers used for producing DHW in all CSHPSS configurations have been simulated through the Type 700; a constant thermal efficiency of 90% is assumed according to the average performance of devices available in the current Italian market.
The operation of the MCHP devices (used in the SCHEMES 2, 3, 5 and 6) has been analyzed through the Type 907. Thomas [63] measured both the electric and overall efficiencies of the ICE-MCHP unit (model Dachs HKA G 5.5) and the SE-MCHP device (model SOLO 161) selected in this paper as a function of the supply temperature (with a constant return temperature). The electric efficiency η el of a MCHP unit is defined as the ratio between the electric output and the primary input, while the overall efficiency η tot is the ratio between the total (electric and thermal) output and the primary input; the difference between η tot and η el is the thermal efficiency, i.e., the ratio between the thermal output and the primary power input. The following regression equations are derived by the authors by interpolating the experimental overall efficiencies measured by Thomas [63] for both the SE-MCHP and ICE-MCHP units: supply − 0.0443 · T 2 supply + 102.74 (4) η ICE−MCHP tot = 0.00006 · T 3 supply − 0.0134 · T 2 supply + 0.8581 · T supply + 77.106 (5) where T supply is the temperature at the outlet of the MCHP devices.
The outputs of the above-mentioned equations have been compared with the experimental values measured by Thomas [63]. The comparison between the experimental data [63] and the predicted ones highlighted maximum deviations equal to about 1.90% and about −0.44% for the ICE-MCHP and the SE-MCHP devices, respectively. Therefore, Equations (4) and (5) have been implemented in the Type 907 and used in this project for predicting the performance of the MCHP devices upon varying the supply temperature.
The operation of the batteries is evaluated through the Type 47a, while the inverter/charge controller is simulated by means of the Type 48b [44]. The Type 47a allows to specify how the charge level of batteries changes as a function of the time, while the inverter converts the DC power to AC and supplies it to the load and/or feeds it back to cover the electric demand. The overall capacity of the 3 series-connected batteries (equal to 40.5 kWh el ) is selected based on the procedure indicated by the manufacturer [56] as a function of the daily average overall electric requirement (equal to about 44.1 kWh el ). Francis [64] investigated the performance of a solar-powered electric vehicle station including a Tesla Powerwall battery [56] with a capacity of 6.4 kWh el . Francis [64] measured the charge/discharge profiles of the 6.4 kWh el Tesla Powerwall along with the corresponding state of charge during the experimental validation in the Electric Sustainable Power Lab at TU Delft; the authors compared the curve representing the state of charge predicted by the Type 47a, used in this paper for modeling the batteries, against that one measured by Francis [64] in the case of same charging/discharging profiles; the comparison indicated that the percentage difference between predicted and measured states of charge is relatively small, ranging from 0% up to 6.7%, and, therefore, demonstrating the suitability of the adopted model.
The Type 742 is adopted for simulating the operation of the district heating network pump; the remaining pumps are analyzed through the Type 656.
The operation of the fan-coils is modeled by means of the Type 753d; this numerical model is detailed according to data provided by the manufacturer [65]. In particular, 5, 6, 8 fan-coils are considered inside the building types A, B and C, respectively. Each fan-coil is characterized by a nominal heating capacity of 0.65 kW as well as a nominal power consumption of 8 W (minimum fan velocity) based on manufacturer specifications [65].
The network of solar thermal collectors as well as the distribution circuit are modeled by considering a single pair of supply and return pipes. The pipes are modeled through the Type 31, where the heat losses are calculated by considering the heat loss coefficient equal to 0.05 kJ/hm 2 K.
An EnergyPlus weather data file [66] is used for estimating the external climatic conditions associated to the city of Naples; the weather data are the same each year due to the fact that they are one year long.

Control Strategies
The heating period covers the period from 15 November up to 31 March, according to the Italian legislation related to the activation of heating systems in the city of Naples [49].
The mass flow rate of the heat carrier fluidṁ FC flows into the fan-coils only in the case of there is a call for thermal energy triggered by a thermostat installed inside the houses. The target of indoor air temperature T room,set-point is set to 20.0 • C (with a deadband of 0.5 • C) only when at least one person being inside the houses during the heating period.
The recovery of solar power is performed depending on the temperature difference DTS = (T SC,out − T 10,STTES ) between the outlet temperature of the solar collectors T SC,out and the temperature T 10,STTES at the node 10 of the STTES. In particular, the flow rate of the solar pump could vary with the time while the solar pump is on; the solar pump is activated and the mass flow rate achieves its lowest value (19.1 kg/h per m 2 of gross area of solar collectors) when DTS is not lower than 10 • C; then, the flow of the solar circuit is modulated between the above-mentioned minimum and the maximum of 38.1 kg/h per m 2 of gross area of solar collectors according to the value of DTS: when DTS is larger than 10.0 • C, the flow is automatically increased every time-step in order to guarantee DTS = 10.0 • C; the flow is automatically decreased every time-step if DTS is lower than 10.0 • C (but greater than 2.0 • C). The minimum and maximum values of the mass flow rate managed by the solar pump are defined according to the manufactured data of the solar collectors [50].
Whatever the plant configuration is, the flow into the HE1 pump is equal to that of one of the solar pumps.
The DH network pump is never deactivated, and it is always operated with a minimum flow rate in order to avoid a significant temperature drop/increase, even if there is no space heating demand. The flow rate of this pump is not constant, and it ranges between a minimum of 497.7 kg/h up to a maximum of 3782.7 kg/h. The maximum flow rate is set according to the nominal flow rates suggested by the fan-coils manufacturer [65] and it is achieved in the case of all fan-coils of all buildings are activated; the minimum flow rate is set according to the nominal flow rates suggested by the fan-coils manufacturer [65] by considering the case when all fan-coils of a single building (A1) only are activated.
Whatever the plant configuration is, the flow rate of the HE2 pump, when it is activated, is equal to that one corresponding to the DH network pump. In particular, the HE2 pump operates according to the temperature levels of hot and cold sides of the heat exchanger HE2 (T in,HE2,hot , T in,HE2,cold ).
The target for the supply temperature of the distribution circuit is set to 55.0 • C; therefore, the back-up system (main condensing boiler of MCHP units) is activated/deactivated depending on the inlet and outlet temperature levels (T in,MB /T in,MCHP and T out,MB /T out,MCHP ).
The set-point temperature to be achieved by the DHW mass flow rateṁ DHW is set to 45.0 • C; the decentralized boilers operate in order to produce DHW at 45.0 • C according to the inlet/outlet levels of temperature (T in,B /T out,B ).
The BTES discharging/charging is operated depending on the temperature levels (T 1,STTES and T 10,STTES , respectively) at nodes 1 and 10 of the STTES, the temperature T BTES,center in the center of the BTES and the given set-point inside the houses T room,set-point . A constant charging flow rate is set to the half of the minimum flow rate of the circuit of solar collectors (equal to 19.1 kg/h per square meter of gross area of solar collectors according to manufacturer data [50]); the discharging flow rate is not constant and it is assumed equal to that one flowing into the distribution network. The selection of the charging/discharging flow rates is based on the criteria suggested by Pahud [61]. Charging and discharging modes of BTES field are alternative; the BTES cannot be charged and discharged simultaneously.
Several studies demonstrated that operating the micro-cogeneration units under an electric-load following control strategy is not efficient and convenient. Therefore, the MCHP units are operated under a thermal-load following control logic; this means that the MCHP units are activated to achieve the given target temperature (55 • C) at the outlet of the heat exchanger HE2, providing their maximum thermal and electric outputs in the case of activation.
The electricity provided by the MCHP units is firstly transferred to the inverter/charge controller and then used to satisfy the electrical requirements, whatever the season is. The cogenerated electricity is sold to the central grid when the electric output of MCHP units becomes larger than the electric demand as well as the state of charge of the batteries is 100%. In the case of the charge level of the electric storages is lower than the maximum, then the possible excess of electric energy generation with respect the electric demand is injected into the batteries; when the electric production is not sufficient to completely satisfy the electric load, then the electric storages are discharged only when their state of charge is bigger than 10%; the discharging of the batteries is interrupted in the case their level of charge becomes lower than 10%. The grid is used to cover eventual peaks of electric demand in the case of the cogenerated electricity is not enough to fully cover the electric demand and the state of charge of the electric storage is lower than 10%. Table 5 reports the main control logics adopted in this paper. Table 5. Control logics of the investigated CSHPSS.

Fan-coils blower & Individual pumps
Heating season AND T room ≤ 19.

Description of the Heating System Assumed as Reference
A Conventional heating System (CS), representative of the current Italian scenario, is modeled and simulated (in TRNSYS environment) while satisfying the energy demands for space heating and DHW of the same district. In the heating system assumed as reference, each house is equipped only with a 26.6 kW th non-condensing boiler fueled with natural gas (modeled through the TRNSYS Type 700 by considering a constant thermal efficiency equal to 90.0%). The heat carrier fluid flows through radiators (modeled through the TRNSYS Type 1231) installed inside the houses. The system is operated with the same control logic used for the CSHPSS schemes. The target temperature of the boilers is set to 55.0 • C for space heating (with a deadband of 0.5 • C). The same profiles of DHW demand, occupancy and electric requirements adopted for the CSHPSS configurations are also used for the conventional system. The target for DHW production is set to 45.0 • C (similarly to the CSHPSS schemes). The electric requirements are fully covered through the electricity provided by the central grid.

Methods of Analysis
In this paper, the simulation results associated to the proposed plant are contrasted with the performance of the conventional heating system. The comparison is carried out by contrasting the consumption of primary energy, the equivalent emissions of carbon dioxide, the operation and capital costs.
In particular, the energy analysis is carried out by calculating the following parameter named Primary Energy Saving (PES): where E CS p is the primary energy consumption of the conventional system and E CSHPSS p is the primary energy consumption of the CSHPSS schemes. The following formulas are used for calculating the PES: E CS p =E th,B /η B + E el,individual pumps + E el,lighting +E el,domestic appliances /η PP (8) where E th,MB is the thermal energy supplied by the main condensing boiler, E p,MCHP is the primary energy consumption of the MCHP devices, E th,DHW is the thermal energy provided by the individual boilers for DHW production for the CSHPSS schemes, E el,import is the electricity imported from the central grid for the CSHPSS schemes, E th,B is the thermal energy provided by the individual boilers for both space heating and domestic hot water production in the case of the reference plant, E el,individual pumps is the electric energy consumption of the individual pumps (P A1 , P A2 , P B1 , P B2 , P C1 , P C2 ) of the reference plants, E el,lighting is the electricity consumptionof the lighting systems of the reference plant, E el,domestic appliances is the electricity consumption of the domestic appliances of the reference plant, η MB is the thermal efficiency of the main condensing boiler [53], η B is the thermal efficiency of the decentralized boilers [53], η PP is the power plant average efficiency in Italy (the value of η PP is considered equal to 0.42 according to the data suggested in [9,67] taking into account also the transmission losses). The values of E el,import are calculated, time-step by time-step, as the difference between the overall electricity demand E el,demand (due to the pumps, heat dissipator, fan-coils, the lighting systems as well as domestic appliances) and the electricity cogenerated by the MCHP units E el,MCHP and/or provided by the batteries (E el,batteries ) in the case of E el,demand > (E el,MCHP + E el,batteries ). The environmental impact associated to both the CSHPSS configurations and the conventional plant is assessed by means of the energy output-based emission factor approach indicated by Chicco and Mancarella [68]. In particular, the equivalent emissions of carbon dioxide are calculated and compared through the following parameter ∆CO 2 : =β · E th,MB /η MB + E p,MCHP + E th,DHW /η B + α · E el,import (10) m CS CO 2 = β · E th,B /η B + α · E el,individual pumps + E el,lighting +E el,domestic appliances (11) where β is the CO 2 equivalent emission factor related to the consumption of natural gas and α represents the CO 2 equivalent emission factor associated to the electricity generation. In particular, α and β are set to 573 g CO2 /kWh el and 207 g CO2 /kWh p , respectively, according to the data indicated in [67,69] according to the current Italian scenario. The comparison from an economic point of view is carried out firstly by contrasting the operation costs of the CSHPSS schemes with those associated to the conventional plant. The comparison is done by calculating the following parameter ∆OC: where OC CS represents the operating costs of the reference plant and OC CSHPSS represents the operating costs of the CSHPSS schemes. The values of OC CS and OC CSHPSS are calculated as indicated below: OC CSHPSS = UC ng · E th,MB / LHV ng · ρ ng · η MB + UC ng,MCHP · E p,MCHP / LHV ng · ρ ng + +UC ng · E th,DHW / LHV ng · ρ ng · η B + UC el · E el,import − UC el,sold · E el,MCHP,exp (13) OC CS = UC ng · E th,B / LHV ng · ρ ng · η B + +UC el · E el,individual pumps + E el,lighting + E el,domestic appliances (14) where UC ng is the unit cost of natural gas consumed by the decentralized boilers  [70] in Italy; this is a consequence of the fact that economic incentives are adopted by the Italian Government to support the diffusion of micro-cogeneration units [72].
In Italy, electricity tariffs that change as a function of energy time of use, in which the electricity costs are divided into more periods throughout the day based on the profile of the demand; this pricing structure provides price signals for the final users, aiming at shifting energy consumption to periods outside the utility expensive peak hours. In particular, the values of UC el depend on (i) the day of the week, (ii) the energy time of use during the day, (iii) the cumulated annual electric energy imported from the grid (in particular, four different intervals of consumption are specified: 0 ÷ 1800 kWh el , 1801 ÷ 2640 kWh el , 2641 ÷ 4440 kWh el , more than 4440 kWh el ) and (iv) the region where the energy is consumed. Figure 3 indicates the weekly profile (from Monday to Sunday) of the unit cost of the electric energy purchased from the central grid (including the excise tax and the VAT) for the city of Naples according to [70]. According to the data reported in this figure, it can be underlined that the unit cost of the electric energy purchased from the central grid: • during the weekdays is larger or equal with respect to the values related to the weekends, whatever the cumulated level of annual electricity imported from the grid is; • increases at increasing the cumulated annual electric energy imported from the grid, whatever the hour and day are; • ranges between a minimum of 12.1 c€/ kWh el up to a maximum equal to 30.1 c€/ kWh el .
In this study, the unit price of the electric energy sold to the central grid UC el,sold is calculated month by month by means of the following formula according to the Italian scenario: UC el,sold = minimum PUN avg ; PZ avg + UC Sf (15) where PUN avg is the monthly average national single price, PZ avg is the monthly average zonal price and UC Sf is the rate devoted to support electric energy generation from renewable sources and cogeneration devices. The values of PUN vary as a function of the hour, ranging between a minimum of 10 €/MWh el up to a maximum of 170 €/MWh el . The values of PZ depend on the hour of the day as well as the month of the year, varying between 32 €/MWh el and 78 €/MWh el . The value of UC Sf is 47.51 €/MWh el in the case of the annual amount of electric energy sold to the grid is not larger than 1800 kWh el , while it is equal to 79.08 €/MWh el for an annual amount of electric energy sold to the grid bigger than 1800 kWh el . The economic analysis is carried out also in terms of capital costs by calculating the so-called Simple Pay-Back (SPB) period, representing the time required to recover the extra initial investment cost. The following formula is used: (16) where CC CSHPSS is the capital cost of the CSHPSS schemes, CC CS is the capital cost of the reference system, EI CSHPSS represents the economic incentives guaranteed by the Italian Government every year (for a maximum of 5 years), OC CSHPSS represents the operation costs of the CSHPSS (Equation (13)) and OC CS represents the operation costs of the conventional heating plant (Equation (14)). The maintenance costs are neglected. The calculation of the difference (CC CSHPSS − CC CS ), representing the extra investment cost of the proposed plant with respect to the conventional system, is carried out by means of the following formula: where CC CSHPSS SC is the investment cost of the solar field collectors (equal to 169 €/m 2 in agreement with [73]), CC CSHPSS SCpumps is the investment cost of the solar field pumps (equal to 1488 € in agreement with the manufacturer data [74]), CC CSHPSS STTES is the investment cost of the short-term thermal energy storage (equal to 5651 € in agreement with the values suggested by Pahud [61]), CC CSHPSS BTES is the investment cost of the BTES equal to 22,917 € (calculated in agreement with the formula reported in [61] as a function of number and depth of boreholes, depth of the soil layer and thickness of the insulation on the top of storage, BTES area and borehole spacing), CC CSHPSS BTESpump is the investment cost of the BTES pump (equal to 2826 € in agreement with the manufacturer data [74]), CC CSHPSS MB is the investment cost of the main condensing boiler (equal to 1436 € in agreement with the data suggested in [75]), CC CSHPSS ICE−MCHP is the investment cost of the two internal combustion engine-based MCHP units (assumed equal to 44,000 € in agreement with the values suggested by Angrisani et al. [76]), CC CSHPSS SE−MCHP is the investment cost of the Stirling engine-based MCHP (equal to 38,950 € in agreement with the data reported in [77]), CC CSHPSS DHnetwork is the investment cost of both pipes and pumps included into the distribution network (equal to 14,500 € in agreement with the price list on public works of the Campania Region [78]) and CC CSHPSS DHWT is the investment cost of the 6 DHWT (equal to 2202 € based on the values suggested by Pahud [61]). The investment cost of the three electric batteries CC CSHPSS Batteries is equal to 18,984 € in agreement with the manufacturer data [56], i.e. about 468.7 €/kWh el ; according to the scientific literature, this specific cost of selected batteries could make economically feasible the adoption of the electric storages.
With respect to the economic incentives, it should be underlined that the Italian Government has dedicated 55 billion € to address the economic consequences of Covid-19 lockdown. These measures include a jump supporting energy-efficient measures associated with such renovation projects, from 50% to 110% of capital costs up to a maximum of 30,000 € for each residential end-user [79]. Homeowners will have three ways in which to secure the eco-bonus: via the transfer of the tax-deductible allowance for installers and product suppliers, through the tax deduction for a period of maximum 5 years, and through invoiced discounts. In this paper, it is assumed EI CSHPSS = 1.10·(CC CSHPSS − CC CS )/5 every year (up to a maximum of 30,000 € for each residential end-user [79]) by considering a tax deduction over a period of maximum 5 years. The decree is already in force and the economic incentive will be used to all expenses incurred between 1 July 2020 and 31 December 2021.

Results and Discussion
The CSHPSS schemes detailed in Figure 2a-d and Table 3 are modeled, simulated and evaluated with reference to a period of 5 years; the performances of the proposed configurations are contrasted with those of the reference plant (described in the Section 3) by using the parameters defined by Equations (6), (9), (12) and (15). Table 6 reports the annual values of primary energy consumption (Equation (8)), equivalent emissions of carbon dioxide (Equation (11)) and operation costs (Equation (14)) for space heating, producing DHW as well as covering the electric demands of the conventional system. The weather data used in the simulation are one year long and, therefore, the simulation results are the same each year. The values reported in this table represent the baseline performance information of the conventional plant to be contrasted with those related to the proposed CSHPSS schemes.  Figure 4 indicates the values of PES (Equation (6)) upon varying the simulation year as a function of the plant configuration. This figure underlines that, whatever the plant configuration is, the PES greatly rises up from the first to the second simulation year; then, it becomes almost constant. Similar trends are derived for the parameters ∆CO 2 (Equation (9)) and ∆OC (Equation (12)); therefore, the performance associated to the fifth simulation year can be considered as representative of a steady-state operation (whatever the scheme of the plant is). The trends indicated in this figure are substantially related to the average temperature of the seasonal storage. Figure 5a indicates the average temperature of the ground in the storage volume T avg,BTES as a function of the time during the first 2 years of operation, while Figure 5b reports the monthly amount of solar energy charged into the BTES and discharged from the BTES during the first two years of simulation. The results indicated in this figure correspond to the SCHEME 5 (that can be considered as representative of all plant configurations).  • the average temperature of the BTES field is initially equal to 10 • C on 1 January and, after that, it rises up to about 73 • C during the summer of the first simulation year; then, it decreases (due to both heat losses and BTES discharging) during the heating season between the first and the second operation years, achieving a value of about 48 • C; finally, it increases again up to a maximum of about 77 • C during the summer of the second simulation year (Figure 5a); • heat discharged from the BTES field (only during the winter of the first two years of the simulation) reaches a maximum of about 1.67 MWh; • the solar energy injection into the seasonal storage occurs every month, except December; it increases up to reaching a maximum during May-July and then reduces; the largest value of solar energy charged into the boreholes is about 3.8 MWh during the first simulation year (in the second year it is lower due to the higher average temperature of the BTES field). Figure 6 shows typical daily operation for the proposed 8 plant schemes during the same day (1 February of the fifth simulation year) in terms of thermal energy flows. In particular, Figure 6a-h corresponds to the SCHEMES 1-8, respectively. These figures indicate the values of net solar power recovered by the solar collectors, thermal power charged into the BTES, thermal power discharged from the BTES and thermal power supplied by the auxiliary system (boiler or MCHP units) as a function of the time. These figures highlight that: (i) the back-up system is often activated in order to compensate the lack of solar energy and maintain the desired supply temperature; (ii) BTES discharging generally occurs during the first part of the day to compensate the temperature levels in the STTES; (iii) BTES charging generally occurs during afternoons/evenings thanks to the solar energy previously stored into the STTES; (iv) solar energy recovered by the solar collectors increases with the time from about 9 a.m. up to reaching a maximum at about 12 a.m. and then reduces.  year) of the 8 configurations: SCHEME 1 (a), SCHEME 2 (b), SCHEME 3 (c), SCHEME 4 (d), SCHEME 5 (e), SCHEME 6 (f), SCHEME 7 (g) and SCHEME 8 (h). These figures shows the daily profiles of the electric power demand, the electric power generated by the micro-cogeneration devices, the electric power charged into the batteries, the electric power discharged from the batteries, the state of charge of the electric storages, the unit cost of electric energy purchased from the grid and the unit price of electric energy sold to the grid as a function of the time. These figures denote that: • the electric power generation occurs according to the activation of the MCHP units (that are operated under a heat-led control logic in order to maintain the desired temperature level at the exit of the STTES) and corresponds to the nominal electric output of the cogeneration devices (11.0 kW el for the SCHEMES 5 and 7 and 9.5 kW el for the SCHEMES 6 and 8); • as detailed in the previous "Section 2. Description of the residential district", the electric power demand is stochastic according to the models suggested by [45]; • for the selected day, the electric power generation is always larger than the electric power demand when the MCHP devices are activated. Whatever the day is, the cogenerated electricity is firstly used to cover the electric load; in the case of the MCHP electric output is bigger than the electric demand and state of charge of the batteries is 100%, then the cogenerated electricity is sold to the central grid; the grid is used to cover eventual peaks of electric demand in the case of the cogenerated electricity is not enough to fully cover the electric demand and the state of charge of batteries is lower than 10%; • the charging/discharging of the batteries can occur during the entire year (both heating and cooling periods) depending on the levels of electric power cogenerated, electric power demanded as well as state of charge of the electric storages (as detailed in the previous "Section 3.2. Control strategies"). In particular, during the selected day the state of charge of the batteries is in the range 53.0% ÷ 84.5% for the SCHEME 5, in the range 40.0% ÷ 74.2% for the SCHEME 6, in the range 34.0% ÷ 79.2% for the SCHEME 7 and in the range 20.0% ÷ 71.5% for the SCHEME 8. The state of charge of the electric storages increases when the electric power production becomes larger than the electric power demanded and vice versa; • electric energy is not sold to the central grid during the selected day thanks to the fact that the state of charge of the electric storages is always lower than 100% and, therefore, there is always room to store the surplus of cogenerated electricity into the batteries; • during the selected day the electric energy demanded is not enough to fully discharge the batteries and, therefore, the state of charge of the electric storages is always larger than the given minimum (10%); • electric energy is not imported from the central grid during the selected day thanks to the fact that the state of charge of the electric storages is always larger than 10%. As a consequence, the required electricity can be discharged from the batteries in case of need; • the unit cost of electric energy purchased from the grid is larger than the unit price of electric energy sold to the grid, with a percentage difference ranging between 15.8% and 21.4% for the selected day during which the cumulated level of annual electricity imported from the grid is in the range 0 ÷ 1800 kWh el . With respect to this point, it should be considered that this percentage difference becomes more and more relevant at increasing the cumulated annual electric energy consumption (as detailed in the "Section 5. Methods of analysis" according to the Italian scenario [70]): (i) in the case of the annual electricity imported from the grid is in the range 1801 ÷ 2640 kWh el , then the above-mentioned percentage difference ranges between 43.4% and 47.2%; (ii) in the case of the annual electricity imported from the grid is in the range 2641 ÷ 4440 kWh el , then the above-mentioned percentage difference ranges between 59.0% and 61.8%; (iii) in the case of the annual electricity imported from the grid is larger than 4440 kWh el , then the above-mentioned percentage difference ranges between 65.5% and 67.8%. The annual electric energy consumptions of the single households served by the plant proposed in this study are larger than 1800 kWh el ; as a consequence, the percentage of energy lost during the batteries charge/discharge process (due to the limited efficiency of the regulator (78%) and the inverter (96.0%)) is more relevant than the percentage difference between the unit costs of electric energy purchased/sold only during the first period of the year (up to the time where the annual electricity imported from the grid is in the range 0 ÷ 1800 kWh el ); this is not true during the remaining part of the year, i.e., when the annual electricity imported from the grid becomes larger than 1800 kWh el . year) of the SCHEME 5 (a), SCHEME 6 (b), SCHEME 7 (c), SCHEME 8 (d).
Similar trends can be derived for the other days of the year, so that it can be stated that (i) it is worth to include an electric storage and (ii) the selected capacity of batteries (determined according to the manufacturer criteria [56]) is feasible. Figure 8a reports the state of charge of the electric batteries (obtained every time-step as one of the outputs of the simulations) as a function of the time during the fifth simulation year for the SCHEME 5. Figure 8b shows the mean daily profile of the state of charge of the electric batteries, where each value is calculated as the arithmetic average (performed every time-step) of the state of charge of the batteries at the same time of different days by considering all days of the heating season only (when the MCHP units, and related electric generation, could be switched on) for the SCHEME 5. These figures demonstrate that: • the state of charge of the batteries reaches the maximum value of 100% during several time-steps in the period between January and Mach ( Figure 8a); • the mean state of charge of the batteries ranges between a minimum of about 11.9% and a maximum equal to about 34.3%; the maximum is obtained at 8:30 a.m., while the minimum is achieved at 3:20 p.m. (Figure 8b). The values of the mean state of charge of the batteries are affected by the limited daily operation time of the MCHP units during November and December because the BTES has been fully charged during the summer and, therefore, it has been able to provide the most part of thermal demand required to achieve the desired supply temperature; • the mean state of charge at first hours of the day (about 25%) is almost the same that at last hours of the day (Figure 8b).  the SCHEME 5 is that one allowing to obtain the best values of PES 5th-year , ∆CO 2 5th-year and ∆OC 5th-year ; in particular, with respect to the reference plant, this scheme allows to decrease the consumption of primary energy as well as the equivalent emissions of carbon dioxide by about 11.3%, while the operating costs are lowered by about 14.3% during the fifth year; • the SCHEME 2 is characterized by the lowest values of PES 5th-year , ∆CO 2 5th-year and ∆OC 5th year ; in particular, with respect to the conventional heating plant, in this case the consumption of primary energy, the equivalent emissions of carbon dioxide and the operation costs are reduced during the fifth year by about 4.1%, 3.4% and 6.4%, respectively; • the simulation data related to the SCHEMES 1-3 highlight that the best auxiliary unit is the natural gas-fired condensing boiler (when the DHW production is performed with the individual boilers only and the electric energy storage is not used); in particular, the values of PES 5th-year , ∆CO 2 5th-year and ∆OC 5th-year are, respectively, about 2.3, 2.2 and 1.4 times larger in the case of the MB is adopted with respect to the cases when the MCHP units are considered; this is mainly due to the temporal mismatch between power generation of the MCHP units (operated under thermal load-following control logic) and electric requirement; • the comparison between the SCHEMES 2 and 3 shows that using the SE-based MCHP device allows to achieve slightly better performance with respect to the case when the ICE-based MCHP unit is adopted (in the case of DHW is produced with the individual boilers and the battery is not used) thanks to slightly higher overall efficiency of the SE-based system compared to that one of the ICE-based co-generator. However, the SCHEME 5 with the ICE-based MCHP unit is more advantageous than that one integrated with the SE-based MCHP device (SCHEME 6) in the case of both the electric energy storage is adopted and the DHW production is obtained by pre-heating the mains water through the DHWTs; this is because the electricity supplied by the ICE-based MCHP system and charged into the batteries (SCHEME 5) is about 1.3 times larger than that one related to the SE-based MCHP coupled with the electric storages (SCHEME 6); • the SCHEMES 7 and 8 are characterized by significantly larger values of PES 5th-year , ∆CO 2 5th-year and ∆OC 5th year in comparison to SCHEMES 2 and 3, respectively. In particular, in comparison to the SCHEME 2, the SCHEME 7 allows to improve the values of PES 5th-year , ∆CO 2 5th-year and ∆OC 5th year by about 2.5, 2.9 and 2.1 times, respectively; similarly, the SCHEME 8 enhances the values of PES 5th-year , ∆CO 2 5th-year and ∆OC 5th year by about 2.1, 2.4 and 1.7 times, respectively, with respect to the SCHEME 3. These results demonstrate how adding the electric storages to the micro-cogeneration devices (without any other modification) greatly enhances the energy, environmental and economic performance of the plant mainly thanks to the self-consumption of MCHP electric output and, therefore, the reduced import of electricity from the central grid; • the comparison between the SCHEME 1 and the SCHEME 4 (differing only in terms of configuration for DHW production) shows that the domestic hot water production by pre-heating the mains water through the DHWTs increases the values of PES 5th-year , ∆CO 2 5th-year and ∆OC 5th-year by about 1.10, 1.11 and 1.01 times, respectively; similar conclusions can be derived by contrasting the SCHEME 5 with the SCHEME 7 (the values of PES 5th-year , ∆CO 2 5th-year and ∆OC 5th-year increase by about 1.1 times) as well as the SCHEME 6 with the SCHEME 8 (the values of PES 5th-year , ∆CO 2 5th-year and ∆OC 5th-year improve by about 1.1, 1.1 and 1.0 times, respectively).  Figure 10 shows the main thermal energy flows of the CSHPSS schemes during the fifth operation year upon varying the configuration plant; in particular, the total solar energy falling on the tilted surface of thermal collectors (equal to 96.8 MWh, whatever the year is, due to the fact that the weather data are the same each year), the net solar energy recovered from the solar collectors, the thermal energy charged into the BTES field, the thermal energy discharged from the BTES field, the thermal energy supplied by the auxiliary unit (main boiler or MCHP units), the thermal energy supplied by the individual boilers for producing DHW are indicated.  Figure 11 highlights the main electric flows, i.e., the electricity cogenerated by the MCHP devices, the electric energy sold to the power line, the electricity charged/discharged into/from the electric storages, the electricity supplied by the MCHP units to directly cover the load and the electricity imported from the central grid upon varying the scheme of the plant during the fifth simulation year. In the same figure, the overall electric energy demand (16.2 MWh) associated to the operation of solar pump, heat dissipator, HE1 pump, HE2 pump, BTES charging/discharging pump, DH network pump, lighting and domestic appliances, individual pumps and fan-coils is also reported.  • the configuration allowing to recover the largest amount of solar energy is the SCHEME 4, while the least significant solar energy recovery is that one associated to the SCHEME 3; • the annual average thermal efficiency of solar collectors, i.e., the ratio between the net solar energy recovered and the solar energy incident on the surface of collectors, ranges between a minimum of 20.5% (corresponding to the SCHEME 3) up to a maximum of 21.9% (corresponding to the SCHEME 4); • the efficiency of the BTES (defined as the ratio between the energy discharged from the BTES and the energy charged into the BTES) ranges between 23.1% and 28.6%, assuming the largest value for the SCHEME 4 and the lowest value in the case of the SCHEME 3. The values of BTES efficiency obtained in this paper are quite consistent with those found in the study carried out by Zhu and Chen [80]; in particular, they performed a huge sensitivity analysis on BTES performance under the climatic conditions of Tianjin (East Coast of China), characterized by a similar latitude with respect to Naples, upon varying the thermal conductivity of soil, the boreholes spacing and depth; they found a BTES efficiency (i) between 20% and 30% in about 17% of simulation cases, and (ii) lower than 50% in about 54% of investigated configurations; • the thermal energy supplied by the auxiliary device is the minimum one when the SCHEME 1 is considered, while it assumes the largest value for the SCHEME 6; • the thermal energy provided by the decentralized boilers for producing DHW is about 0.7 times lower for the SCHEMES 4-6 with respect to the SCHEMES 1-3, respectively; this is thanks to the fact that for the SCHEMES 4-6 the energy demand for DHW production is partially covered by pre-heating the mains water through the decentralized DHWTs; • for the SCHEMES 2 and 3 the electricity imported from the central grid is about 0.95 times lower than that one associated to the case SCHEME 1 thanks to the cogenerated electricity partially covering the overall electric demand; • in comparison to the SCHEMES 2 and 3 (without the electric energy storages), respectively, the SCHEMES 7 and 8 (integrated with the batteries) allow to significantly decrease the import of electric energy from the central grid by about 10.2% and 7.6%, respectively. • the SCHEMES 5-8 (including micro-cogeneration units as back-up system together with electric storages) allow to avoid exporting the cogenerated electricity into the central grid with low revenues (so that the produced electric energy is directly used to cover the load or charged into the batteries).
With respect to the last point, it should be considered that, according to the Italian scenario, the values of UC ng,MCHP range from a minimum of 0.43754425 €/Nm 3 up to a maximum of 0.73602551 €/Nm 3 depending on the cumulated level of natural gas annual consumption. Therefore, assuming both a lower heating value of natural gas equal to 49,599 kJ/kg (representative of the scenario associated to the city of Naples) as well as a density of natural gas equal to 0.72 kg/m 3 , the maximum value of UC ng,MCHP in the case of the proposed application is around 0.022 € for producing 1 kWh el of electric energy by means of the micro-cogeneration devices investigated in this paper (according to the nominal electric and thermal efficiencies suggested by the manufacturers [53,54]). This means that the unit cost of electric energy cogenerated is lower than the unit cost of electric energy sold to the grid (about 0.097 €/kWh el ) and, therefore, the losses of batteries charge/discharge process could be compensated, and the proposed application could be economically viable.
In addition, according to the Italian scenario [67,69], the proposed application also appears environmentally feasible taking into account that: • the equivalent emissions of carbon dioxide associated to the import of 1 kWh el of electric energy from the central grid are equal to 573 g CO2 ; • the equivalent emissions of carbon dioxide associated to the cogeneration of 1 kWh el of electric energy through the MCHP devices selected in this study are equal to about 62 g CO2 .
The feasibility of batteries depends on several factors, such as electric demand profile as well as electric generation profile. Figure 12 reports the difference between the electric power produced by the microgeneration units and the demanded electric power as a function of the time with reference to the entire fifth year of operation of the SCHEME 5 (characterized by the best performances). The positive values of this difference represent the excess of power generation, i.e., the electric power that could be potentially charged into the batteries and/or sold to the central grid; the negative values in the graph are related to the electric power to be discharged from the electric storages and/or purchased form the central grid in order to cover the electric demand. there is an excess of electric energy generation (cogenerated electricity larger than demanded electricity) only during the heating period because of the fact that, as indicated in Table 5, the micro-cogeneration units (with associated electric output) can be activated only during this interval of time (under heat-led control logic); as a consequence, electric energy can be charged into the batteries and/or sold to the central grid only during the heating season; • during the cooling season the MCHP units are switched off and, therefore, the electric demand is almost entirely covered by the central grid.
The application of electric storages is also affected by cost unit cost of the electric energy purchased from the central grid UC el as well as the unit price of the electric energy sold to the central grid UC el,sold . In this study, the values of UC el vary between a minimum of 0.121 €/kWh el up to a maximum equal to 0.301 €/kWh el (according to the contents of Figure 3), while the average value of UC el,sold is about 0.097 €/kWh el . This means that the values of UC el are significantly larger that the values of UC el,sold , justifying the inclusion of electric storages into the plant schemes. In addition, it could be underlined that the Tesla Powerwall battery solution (adopted in this work) can be considered as a viable investment in regions where the values of UC el are over 0.25 €/kWh el according to [35].
The solar fraction SF 5th-year of the CSHPSS is calculated as the ratio between the thermal energy demand that is covered thanks to the solar source and the overall thermal input to the CSHPSS schemes during the fifth operation year: where E CSHPSS th,HE2 is the solar energy moved into the DH circuit by means of the heat exchanger HE2, E CSHPSS th,back−up is the thermal energy supplied by the auxiliary unit (main boiler or MCHP devices) and E CSHPSS th,DHW is the thermal energy provided by the decentralized boilers for producing DHW during the fifth simulation year. Figure 13 indicates the values of SF 5th-year as a function of the configuration plant. This figure shows that: • the largest value of SF 5th-year , equal to 31.5%, is associated to the SCHEME 4 (using the boiler as auxiliary unit); this means that 31.5% of the overall thermal energy demand is satisfied by means of solar source in the above-mentioned configuration; • the lowest value of SF 5th year is that one obtained for both the SCHEMES 3 and 8; it is equal to 26.1%.
The values of solar fraction found in this study are consistent with those associated to similar studies carried out under Italian climatic conditions [40,42,43], highlighting values of SF in the range between 28.2% and 80.0%. Table 7 displays the values of the Simple-Pay Back (Equation (16)), i.e., the time required to recover the extra initial investment cost related to the CSHPSS schemes in comparison to the conventional scheme, as a function of the plant configurations. • the variation of SPB period is almost negligible upon varying the plant configurations, ranging from a minimum of 4.39 years (corresponding to the SCHEME 1) and a maximum of 4.49 years (corresponding to the SCHEME 2); • the proposed configurations are all feasible from an economic point of view thanks to the short time required to recover the investment.
With respect to the values reported in Table 7, it should be highlighted that they are greatly affected by the relevant economic incentives put in place by the Italian Government; without this support, the SPB periods would be much larger and the proposed configurations would not be economically feasible. In particular, the capital costs of the MCHP units are significantly consistent (about 44,000 € for the ICE-MCHP unit and about 38,950 € for the SE-MCHP device) and significantly contribute to the overall investment cost associated to the CSHPSS; the capital cost related to the BTES is also relevant, accounting for about 22,917 €.

Conclusions
During the last years, solar district heating systems achieved a great attention in the scientific community, with a significant number of worldwide applications, thanks to several benefits they can obtain with respect to decentralized plants. The review of scientific literature indicates that the studies carried out with reference to Italian applications are really limited and generally dated, with just two works [40,42] performed under climatic conditions associated to the southern Italy. Moreover, most of these studies are based on an energy analysis only, generally neglecting the economic as well as the environmental aspects. In addition, the size of the analyzed districts is significantly bigger than that one investigated in this work. Finally, the papers available in the literature are not focused on systems integrating cogeneration units and/or electric energy storages.
In this paper, the performances of a solar hybrid DH network including a seasonal borehole thermal energy storage are investigated during five years upon varying the auxiliary back-up system (boiler or micro-cogeneration units), with and without the utilization of electric energy storages in combination of the MCHP units; heat recovery from the distribution circuit for DHW production is also analyzed. The performances of the proposed plants have been contrasted with those of typical Italian decentralized heating plants.
The main scopes of the study can be summarized as follows: • the feasibility of micro-scale solar district heating networks in comparison to conventional heating systems in the case of Italian scenario; • the performance of long-term thermal storages in the case of micro-scale solar DH networks; • the most suitable back-up technology to be adopted into micro-scale solar district heating networks; • the capability of micro-cogeneration units in satisfying the requirement of electricity associated to small Italian district composed of residential buildings only; • the feasibility of batteries in improving the self-consumption of cogenerated electric energy and, consequently, decreasing the amount of power bought from the central electric grid; • the benefits deriving from the possibility to pre-heat the mains water for DHW production by recovering heat from the distribution circuit through decentralized tanks.

•
The comparison between the proposed plant configurations and the conventional heating system highlighted that: • all schemes are able to reduce the consumption of primary energy, the equivalent emissions of CO 2 and the operation costs with respect to the reference system during the fifth operation year; • the largest reduction in terms of primary energy consumption, equivalent emissions of carbon dioxide and operation costs are equal to 11.3%, 11.3% and 14.3%, respectively; • the efficiency of the seasonal storage is significant, ranging between 23.1% and 28.6%; • the largest solar fraction achieved during the fifth simulation year is about 31.5%; • in the case of the batteries are not used, the best auxiliary technology is the natural gas-fired condensing boiler; • the configuration using the ICE-based MCHP device as back-up system coupled with the batteries and pre-heating the mains water through the heat recovery from the district network (SCHEME 5) is that one allowing to obtain the largest savings; • the thermal energy provided by the decentralized boilers for producing DHW is about 0.7 times lower for the 3 configurations where the pre-heating of mains water is performed through the heat recovery from the district network in comparison to the cases where the domestic hot water production is based on the operation of individual boilers only; • with respect to the schemes using the MCHP units without the batteries, the addition of the electric energy storages only allows to significantly reduce the import of electricity from the central grid, the equivalent emissions of carbon dioxide and operation costs up to 1.1, 2.9, and 2.1 times, respectively; • the proposed schemes are economically feasible (SPB periods of about 4.5 years) mainly thanks to the significant economic incentives put in place by the Italian Government.
In the future, the authors would extend this work to consider the integration with eco-friendly chillers that could be driven by solar energy (such as adsorption units) with the aim of covering the cooling demand of the district and assess its energy, environmental and economic feasibility with respect to conventional heating and cooling systems. In addition, one of the next research steps would consider the utilization of excess of electrical energy generation (produced energy less demanded energy) to power vapor-compression electric heat pumps instead of adopting electric storages. Moreover, future researches will be also focused on the potential application of the proposed schemes under different operating scenarios upon varying the buildings loads and/or the climate conditions. Finally, the assessment of environmental impact will be performed by means of a life-cycle analysis (LCA) in order to take into account all the stages of the life-cycle associated to the proposed components.

Conflicts of Interest:
The authors declare no conflict of interest.