1. Introduction
District heating and cooling (DHC) is considered more efficient than individual, distributed systems for heating and cooling, especially because DHC solutions can benefit from locally available, low-cost energy sources, like environmental heat and cool, industrial waste heat and solid waste incineration [
1], but DHC, as a heat/cool demand aggregator with thermal storage capacity, can offer also flexibility in managing energy demand. In the future power grid dominated by renewable non-programmable electricity, this feature would contribute to make DHC solutions based on power-to-heat/cool technologies of strategic importance [
2]. Existing DHC systems using large size, electricity driven heat pumps are examples in this field. These systems can provide social and environmental benefits when used to replace heating technologies relying on combustion in densely populated urban areas, where air pollution is an impelling problem. Concerning heat pump technologies, different low-grade heat sources have been used in existing district heating (DH) plants, including industrial excess heat, ambient water and sewage water. Sewage water is widely adopted in Sweden, where DH plants utilizing sewage water represents about 50% of the heat generated by power-to-heat DH plants [
3]. Regarding power-to-cool conversion, examples of large DHC systems have been proposed in Japan [
4], in which river water constituted the most efficient heat sink.
Wastewater heat recovery applications based on heat pumps are becoming widespread in energy saving applications for both heating and cooling. Heat recovery can be performed inside the buildings (domestic), from sewerage lines (urban) and from wastewater treatment plants (municipal). In the review conducted in [
5], COP values in the range of 1.77–10.63 for heating and 2.23–5.35 for cooling are reported based on the experimental and simulated values. Moreover, the sewage heat exchanger is one of the key components in wastewater-source heat pumps. Different types of sewage heat exchangers have been used, including shell & tube, plate, spiral, gravity-film and channel type [
6]. In domestic utilization, a consistent amount of heat energy can be recovered from washers, dishwashers and showers. Gravity film and spiral heat exchangers can be used directly in the drainage system. In wastewater heat recovery at urban scale, both channel type and external (e.g. shall & tube, plate) heat exchangers are used. External heat exchangers are more effective although they require wastewater screening and extra pumping and piping system. The heat recovery at municipal scale is technically easier, although water treatment plants are seldom close to the consumer.
In densely populated areas, heat recovery from sewage at urban scale has a large potential, as shown in [
7] where, through a GIS-based analysis that matched availability of sewage and heat demand, high utilization factors of the heat theoretically recoverable at the final treatment plant are found for different sewerage lines in Tokyo. However, untreated urban sewage is not widely used due to the problem of filth. Auto-avoiding-clogging equipment can be used to continuously capture suspended solids in the sewage as in [
8], where an untreated sewage source heat pump (USSHP) system is experimentally investigated showing COP of the heat pump unit and of the system of 4.3 and 3.6, respectively. The operational experience with another type of filth block device is reported in [
9], where an urban sewage source heat pump system composed of a filth block device, a wastewater heat exchanger, and a heat pump is demonstrated. The results indicated that in typical conditions the heating COP is about 4.3 and the cooling COP is about 3.5.
Based on these premises, the energy and economic feasibility of a DHC system based on a large size heat pump and complemented by sewage heat/cool recovery is investigated. The system is conceived to alternatively supply heat and cool to a commercial district located in northern Italy, a densely populated area where air quality is a main concern and commercial buildings, characterized by heating and cooling loads of comparable magnitude (respectively about 100 kWh/m
2 year and 75 kWh/m
2 year [
10,
11]), highly contribute to the local thermal energy demand and the associated emissions of air pollutants. Communal sewage provides the heat source in heating and the heat sink in cooling. With respect to sewage potential and availability, it is worth noticing that Italy is characterized by a large per capita water consumption, 175 liter/day/person [
12]. Such a system would offer the following main advantages: (1) ability to exploit the increasing share of renewable electricity in the electricity mix, which implies lower CO
2 emissions as compared to gas boilers for heating; (2) ability to purchase electricity on the wholesale market at competitive tariffs, as compared to the average electricity consumer; (3) zero local emissions of air pollutants (e.g., PM, NOx), as compared to gas boilers and DH conversion systems relying on combustion (e.g., CHP); (4) possibility to exploit PV electricity generated on-site by a large capacity PV plant that can benefit from economies of scale. The analysis aims to estimate the overall plant efficiency, including heat losses and parasitic energy consumption, and the economic competitiveness and the environmental benefit as compared to individual electrical heat pumps, which represent an alternative of comparable social and environmental impact for the location considered. To answer these questions, a detailed mathematical model of the system is built and simulated in Trnsys [
13].
2. Plant Configuration and Operation
The general scheme of a district heating and cooling system using sewage heat recovery is displayed in
Figure 1. Communal sewage is accessible at relatively close distance from the thermal energy consumer. After screening, sewage exchanges heat with a vapor compression heat pump, located in the machinery room, that supplies thermal energy to the consumer through a heating and cooling network.
In the plant configuration considered within this work, the DHC network supplies heat, in winter, and cool, in summer, to a commercial district comprising distinct department stores. The seasonal switch from heating to cooling occurs typically at end of April and the cooling season lasts until mid-October. Due to the simultaneity of the thermal demand across buildings of similar characteristics, a two-pipe network is used. At the end-user, hot water is supplied at 50 °C with a peak load temperature drop of 10 °C and chilled water at 6 °C with peak rise of 7 °C. PV panels are installed on a portion of the overall flat surface available on buildings roofs (see
Figure 2).
The machinery room (see
Figure 3), located underground, comprises sewage basin with screening unit, elevation pump, sewage heat exchanger, heat pump and hot/cold water storage. The heat pump operation mode can change from heating to cooling through inversion on the external water loops (not shown in
Figure 3).
The electricity of the PV system is used to drive the system pumps and the heat pump, in conjunction with the electricity purchased from the grid. If PV electricity is generated in excess, the surplus is supplied to the electricity grid. The remuneration of the electricity supplied to the grid is determined according to the net metering regulation currently in force [
14], which allows selling at a tariff constituted by the sum of the wholesale electricity price and a contribution related to the system charges and general transmission and distribution costs.
The control strategy is defined with the objective to limit parasitic energy consumption. Therefore, variable speed pumps are used. Pump P2 operates to keep the water storage at the desired set point temperature. The heat pump compressor frequency is modulated to: (i) deliver water at the desired temperature to the water storage when heat pump capacity is overabundant (partial load); (ii) limit above 0 °C the temperature of the water sent to the SHX in heating operation to prevent the evaporator from freezing (antifreeze protection); (iii) limit below 50 °C the temperature of the water sent to the SHX in cooling operation to prevent excessive heating of sewage (overheating protection). Pumps EP and P1 modulate their speed based on the compressor speed. Pump P3 is controlled according to the network return temperature, as further explained in the following section.
3. Mathematical Model
In the following, the mathematical models of the main system components are presented. The models are calibrated and validated against the manufacturer’s data of a reference unit. Their outputs, once scaled with respect to the capacity of the reference unit, are assumed representative for units slightly different in size.
3.1. Heat Pump
The heat pump model predicts condenser heat duty (
), evaporator heat duty (
) and compressor power input (
) at different values of chilled water inlet temperature (
), chilled water mass flow rate (
), hot water intel temperature (
), hot water mass flow rate (
) and compressor frequency (
) by means of thermodynamic modelling. The reference unit is a screw chiller of 788 kW heating capacity using R134a as refrigerant [
15]. Based on the refrigerant property [
16], the cycle state points are determined by imposing constant temperature differences for the following temperature pairs: refrigerant condensation and hot water outlet (0 K), hot water inlet and refrigerant at condenser outlet (1 K), chilled water outlet and refrigerant evaporation (1 K), chilled water inlet and superheated vapor at evaporator outlet (1 K). Isentropic efficiency as function of compression ratio (
) and frequency ratio (
) is used for calibrating the model output with the performance data derived from the manufacturer’s datasheet. The resulting functional dependency is shown in
Figure 4.
Volumetric flow rate at compressor inlet is assumed proportional to frequency ratio and its maximum value (
) is calibrated based on the heat pump nominal heating capacity (
). The resulting scale factor
is 0.331 m
3/s MW. The comparison between model output and manufacturer’s data is presented in
Table A1, showing very good accuracy in different operating conditions.
3.2. Sewage Heat Exchanger
The sewage heat exchanger (SHX) comprises screening unit, elevation pump, and shell and tube heat exchanger. Sewage circulates in the tubes and clean technical water flows inside the shell, in a counterflow arrangement (see
Figure 5). Concerning the parasitic energy consumption of the SHX, sewage elevation is the main cost. As the partial load value in the simulations is generally very high, the exponential dependence of the parasitic energy consumption on the sewage water flow rate is linearized. The proportionality factor (
), estimated based on experimental data, is provided in
Table A2.
One of the main issues with the SHX is sludge accumulation at the inlet sections of the tube banks, which progressively decreases flow passage area (
), sewage mass flow rate (
), and ultimately heat transfer rate. Cleaning is thus necessary to periodically restore the nominal performance. To calculate the influence of performance degradation, a model for the prediction of sewage flow rate is developed starting from the balance equation for the mass of sludge:
where
is sludge deposition,
is sludge removal and
is the mass of sludge accumulated at the entrance region of the tubes bank. According to the fouling model of Kern & Seaton [
17],
and
where
is the shear stress. With the hypothesis of (1) negligible removal, (2) constant pressure drop across inlet and outlet sections, (3) constant friction factor and (4) constant hydraulic diameter (i.e. assuming occlusion of one pipe at a time), Expression (1) can be recast as:
where
is a constant,
is the pressure drop across inlet and outlet sections, and
is the flow passage area which is equal to
after cleaning. The reduction of
is not directly proportional to sludge accumulation because initially sludge is likely to accumulate in stagnation zones which do not contribute to flow passage area. For the derivation of an approximated empirical law, a critical sewage mass (
), after which
becomes proportional to
through a constant proportionality factor (
), is assumed:
In conclusion, a differential equation for
is obtained in the form:
Since under the current simplifying assumptions
, the following expressions for
can be derived:
where
and
are identified experimentally (see
Figure 6) and their values are provided in
Table A2.
The result can be generalized to the case of variable pressure drop, assuming constancy of pressure drop during a timestep interval [
;
]:
where
. Expressions (7) and (8) are easily implemented in transient simulations and allow considering variable flow control (
) and periods during which the SHX is not in operation (
). Lastly, the UA value of the SHX is determined from the well-known expression of
valid for counterflow heat exchangers, with
the fraction of the heat capacities of the media considered smaller than 1:
Effectiveness () and the mass flow rates of each stream are measured in clean conditions and the corresponding UA is calculated. With the accumulation of sludge at the entrance sections, the degradation of the UA is idealized as the consequence of the decrease in effective heat transfer area resulting from the progressive blocking of the internal tubes. Therefore, it is assumed that UA is directly proportional to .
3.3. District Heating and Cooling Network
The pipeline diameter is sized according to the design volumetric flow rate and pressure drop. Since the simulation timestep (one hour) is about four/five times as large as the water transition time through the pipeline, a one-node lumped capacity model with heat losses is used for both forward and return pipes. The control law of the mass flow rate (
) is based on the return temperature and assumes that load (
) is initially modulated by lowering the mass flow rate in the attempt to keep constant the temperature difference across supply and return (
), and by lowering the
when mass flow rate is towards its minimum value. A suitable, arbitrary mathematical formula that reproduces the control law of
with
is Expression (10), from which Expression (11) can be derived:
The parasitic consumption of the pump is estimated based on mass flow rate, hydraulic characteristic of the pipeline, and pump efficiency. The hydraulic characteristic can be estimated based on the network design parameters, whose values are presented in
Section 4.3.
3.4. Heating and Cooling Loads
The district heating plant must be sized according to both heating and cooling loads, therefore suitable heating and cooling hourly profiles must be generated using a building energy model.
The two-node capacitive building model is shown in
Figure 7, whose states are room air temperature (
) and inner building envelope temperature (
) and
refers to the ambient temperature. Such a model is selected because as compared to more sophisticated models, it requires a minimum set of input data that can be easily tuned based on floor area, indoor volume, building envelope insulation characteristics, and target specific heating and cooling demands:
UA between room air and outdoor ambient air (, W/K)
UA between room air and the inner building envelope (, W/K)
Thermal capacity associated to the room air (, J/K)
Thermal capacity associated to the inner building envelope (, J/K)
Occupation density (p/m2)
Internal gain, e.g. lights (W/m2)
Ventilation (vol/h)
Infiltration (vol/h)
Fraction of solar radiation transmitted to room air (-)
Daily comfort hours (from-to h)
Heating and cooling daily schedule (from-to h)
Cooling setpoint temperature (°C)
Heating setpoint temperature (°C)
Based on outdoor ambient temperature and radiation on the horizontal plane, the model calculates solar input (
), internal gain (
), heat gain of the occupants (
), sensible and latent heat gain of ventilation and infiltration (
and
), and predicts heating
and cooling (
loads applying the following energy balances:
3.5. Key Performance Figures
The following energy and economic performance figures are introduced to evaluate the performance of the system:
- (1)
Heat pump seasonal COP:
where
is the cumulated heat delivered by the condenser in heating mode operation and
is the associated electrical consumption of the compressor.
- (2)
Heat pump seasonal EER:
where
is the cumulated cool delivered by the evaporator in heating mode operation and
is the associated electrical consumption of the compressor.
- (3)
System level seasonal COP:
where
is the cumulated heat delivered to the users net of heat losses through storage and network and the denominator comprises the associated electricity consumption of the heat pump compressor
and all pumps (
,
,
,
).
- (4)
System level seasonal EER:
where
is the cumulated cool delivered to the users net of cool losses through storage and network and the denominator comprises the associated electricity consumption of the heat pump compressor and all pumps.
- (5)
Full load equivalent hours:
where
is the nominal power of the condenser and
is the nominal power of the evaporator.
- (6)
Levelized cost of heating and cooling:
where
is calculated for each plant component by dividing investment cost times annuity factor, function of discount rate and useful life.
includes maintenance costs, evaluated as a percentage of the investment costs on machineries, and cost of cleaning for the sewage heat recovery system. Moreover, annual purchases of electricity contribute to total
. Annual benefits (
) include the sales at wholesale price (
) of PV electricity generated in excess and the contribution related to avoided cost of transmission and distribution (
), calculated according to the following expression [
14]:
where the valorization of the electricity sold to the grid (
) is capped by the electricity purchased (
) and valorized at the national average selling price (PUN).
4. Case Study
The study focuses on a commercial district located in Milan, Italy and comprising distinct department stores. The Milan area is characterized by a large per capita water consumption, about 175 liter per day, and a large sewage network that collects wastewaters and supplies three wastewater treatment plants located in the external ring of the city. With an estimated total average wastewater flow rate of about 19,000 m3/h, the potential sewage heat recovery amounts to about 110 MWt (assuming a temperature difference of 5 K). In the following, the external conditions and the input figures of the study are presented.
4.1. External Conditions
The external conditions for the selected location consist of the following sets of hourly data:
Outdoor ambient temperature
Global solar radiation on the horizontal plane
Diffuse solar radiation on the horizontal plane
Wholesale electricity price
Sewage water temperature
The hourly profiles have been collected for the same reference year, 2017. Such profiles are needed to calculate heating and cooling loads, predict PV electricity generation, and estimate sales of PV electricity generated in excess. A specific year (2017) has been selected instead of a meteorological standard year because actual electricity prices are considered and the deviations in heating and cooling degree days of 2017 were minor compared to the typical meteorological year.
Monthly solar irradiation and outdoor temperatures for Milan [
18] are shown in
Figure 8. The weather conditions in Milan are characterized by an average daily irradiation on the horizontal plane equal to 3.8 kWh/m
2, mild cold temperatures during most of winter months and hot temperatures in the central summer months.
Energy prices on the electricity wholesale market [
19] are shown in
Figure 9, including the national average selling price and minimum, maximum, mean zonal prices for northern Italy. The maximum selling price deviates significantly from the average selling price, thus showing the potential economic margin deriving from PV electricity sold to the grid. The overall tariff related to the excess electricity sold to or purchased from the grid includes the contribution related to system charges and transmission and distribution costs, which is estimated equal to 55 €/MWh.
Concerning sewage water, its temperature is subject to seasonal variation due to the influence of weather. Based on direct measurements of treated water temperatures in Milano [
12], a sinusoidal fluctuation from 13 °C in winter to 23 °C in summer can be assumed. This assumption is also supported by the measurements of the sewage temperature at the DHC in Budapest corrected by the shift between the average ambient temperatures of Budapest and Milano (2.5 K).
4.2. Heating and Cooling Loads
The total floor area of the commercial district is estimated at 12,000 m
2. Based on the building model input figures (see
Table A3), heating and cooling hourly loads are calculated. The corresponding monthly values are shown in
Figure 10 and
Figure 11. The specific annual heating and cooling demands are 104 kWh/m
2 and 76 kWh/m
2, while the respective specific peak thermal loads are 87 W/m
2 and 83 W/m
2. These values are in line with the typical heating and cooling demands in the commercial sector in northern Italy. The overall heating and cooling annual demands are 1251 MWh and 912 MWh, and the peak loads are respectively 1050 kW
t and 1000 kW
t.
4.3. Plant Sizing
The DHC network supplies heat at 50 °C/40 °C and cool at 6 °C/13 °C. Based on the temperature differences and thermal loads, mass flow rates are derived. The preliminary sizing of the heat pump is performed according to the peak loads. As a first design value, a heat pump of
equal to 1100 kW is selected, and the heat duty of the sewage heat exchanger is determined as the maximum between evaporator heat duty in heating mode and condenser heat duty in cooling mode. Concerning the DHC network, the main design parameters are diameter and length of the pipeline. Pipeline diameter influences pressure drops and, ultimately, the parasitic energy consumption associated to the circulation pump of the DHC network. To limit this consumption, the design pressure drop at the circulation pump is fixed at 350 kPa, of which 50 kPa are associated to heat exchangers. Thus, the pipeline dimeter is selected to meet the net design pressure drop of 300 kPa at the required design volumetric flow rate. The DHC network design parameters are shown in
Table A4. The PV system area is an optimization parameter, since the remuneration of PV electricity sold the grid is penalized if the electricity in excess is larger than the electricity purchased. The optimum is achieved when total electricity exported to the grid is equal to the total electricity imported from grid on a yearly basis.
4.4. Cost Parameters
The estimated specific investment costs [
20,
21,
22] are reported in
Table A5, along with the associated useful life and the resulting annuity factor, evaluated at an interest rate of 2.5%. Concerning maintenance, yearly maintenance costs are estimated equal to 1.5% of the investment cost of machineries. Lastly, the cost of each cleaning operation of the sewage heat exchanger is estimated at 500 €, assuming that the cleaning operation requires the work of one specialized technician for ten hours.
4.5. Reference System
A reference system based on multiple, independent, reversible air-to-water heat pumps of small capacity is chosen as term for comparison with the system of interest (see
Figure 12).
The individual heat pumps are assumed to deliver hot water at 50 °C for heating and chilled water at 7 °C for cooling. The specific investment cost of each unit is estimated at 340 €/kWh [
23]. Hydraulics and electrical wiring costs are 30% of investment, and yearly maintenance is set to 1.5% of the investment. Since air-to-water units are installed outdoor, useful life is set to 10 years. At 2.5% interest rate, the corresponding annuity factor amounts to 8.97. In the reference scenario, the average cost of electricity for commercial users, equal to 150 €/MWh [
24], is used. The overall installed capacity is determined by the peak heating load, since the heating capacity of air-to-water heat pumps is greatly reduced at the heating design condition. The installed heating capacity amounts to 1500 kW
t (at air temperature equal to 7 °C), which corresponds to 1150 kW
t at heating design condition (air temperature equal to −5 °C). The seasonal energy performance in heating and cooling is calculated dividing the heating and cooling hourly loads respectively by the instantaneous values of COP and EER, as provided by the manufacturer of a typical air-to-water heat pump [
23]. The calculated seasonal COP and EER achieve respectively 2.41 and 3.27. The costs and the energy performance of the reference system are shown in
Table A6,
Table A7 and
Table A8, while the calculation of its LCOHC, equal to 95.6 €/MWh, is reported in
Table A9.
6. Conclusions
According to the results of the present research, which are based on the detailed modelling of partial load operation and year-round simulations, a power to heat/cool system comprising sewage heat recovery, heat pump, storage, distribution network with substations, can be competitive with respect to individual, distributed, electricity driven heat pumps. With the addition of large-scale PV, the system can be even more competitive thanks to the remuneration of electricity sold to the grid and the large share of self-consumed PV electricity associated to the heating and cooling.
In the case of interest of this study, a commercial district located in northern Italy, the proposed DHC system can be an attractive solution. The overall LCOHC is 79.1 €/MWh against 95.6 €/MWh for the reference system considered (individual heat pumps). The reasons for this are found in: the superior energy performance of the DHC system, with system COP and EER of 3.10 and 3.64, higher than the respective values for the reference system (respectively 2.41 and 3.27); the longer useful life of the DHC system; and the largely positive contribution of the PV system that lowers LCOHC by 10.8 €/MWh. It shall be noticed that the superior energy performance of the DHC system can justify the investment even in absence of PV, with a LCOHC of 89.9 €/MWh. In this case, the electricity saving and the associated avoided CO2 emission amount to 211 MWh/y and 78 tCO2/y, respectively.
Although the results are specific for the location of interest and the regulation on distributed renewable electricity generation currently in place, they are promising. Commercial districts are interesting targets for DHC since heating and cooling intensity is expected to remain high. With new energy efficiency measures in buildings, cooling is expected to gain more importance than heating. Moreover, sewage is available with continuity of flow in large cities and represents an interesting thermal medium for heat source (in heating) and rejection (in cooling). The coexistence of heating and cooling loads in commercial buildings with sewage availability should be the driver for future researches on the potential of heating and cooling technologies complimented by sewage heat recovery. However, the analysis has also shown that the sewage heat recovery system is a rather expensive technology, and margin should exist to lower the associated investment costs and foster the replication of this solution. Moreover, frequent cleaning operations are needed. Cleaning frequency is tightly linked to sizing and an optimum size of sewage heat recovery system, heat pump and heat storage that allows maintaining cleaning at a manageable, cost-effective frequency must be pursued at design stage.