Next Article in Journal
Online Knowledge Learning Model Based on Gravitational Field Theory
Next Article in Special Issue
Opportunities for Integrating Underground Railways into Low Carbon Urban Energy Networks: A Review
Previous Article in Journal
An EKF-Based Method and Experimental Study for Small Leakage Detection and Location in Natural Gas Pipelines
Previous Article in Special Issue
Contracts, Business Models and Barriers to Investing in Low Temperature District Heating Projects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integration of Reversible Heat Pumps in Trigeneration Systems for Low-Temperature Renewable District Heating and Cooling Microgrids

1
DESTEC Department of Energy, Systems, Territory and Constructions Engineering, University of Pisa, Largo L. Lazzarino, 56122 Pisa, Italy
2
Department of Mechanical Engineering, Universitat Rovira i Virgili, Avda. Països Catalans 26, 43007 Tarragona, Spain
*
Author to whom correspondence should be addressed.
Submission received: 21 June 2019 / Revised: 15 July 2019 / Accepted: 31 July 2019 / Published: 5 August 2019
(This article belongs to the Special Issue Urban District Heating and Cooling Technologies)

Abstract

:
District heating and cooling networks based on trigeneration systems and renewable energy technologies are widely acknowledged as an energy efficient and environmentally benign solution. These energy systems generally include back-up units, namely fossil-fuel boilers and electric chillers, to enhance system flexibility and cover peak energy demand. On the other hand, 4th generation district heating networks are characterized by low-temperature water distribution to improve energy and exergy efficiencies. Moreover, reversible heat pumps are a versatile technology, capable of providing both heating and cooling, alternately. In this paper, the integration of reversible heat pumps as single back-up units in hybrid renewable trigeneration systems serving low-energy micro-district heating and cooling networks is investigated. A detailed modeling of the system is provided, considering part-load and ambient condition effects on the performance of the units. Size and annual operation of the proposed system are optimized in a case study, namely a large office building located in Pisa (Italy), by means of a genetic algorithm-based procedure. A comparison with the conventional trigeneration system is performed in terms of economic and environmental perspectives. Results show that the integration of reversible heat pumps is an economically viable solution capable of reducing by 7% the equivalent annual cost, increasing the installed power of renewables up to 23%, and lowering by 11% carbon dioxide emissions, compared to the energy system with conventional back-up units.

1. Introduction

Energy consumption worldwide grew by 2.3% in 2018, due in large part to higher heating and cooling needs. As a result, global energy-related CO2 emissions increased by 1.7%, up to 33.1 Gt/yr. Energy efficiency and renewable energy sources (RES) are considered the main weapons to reduce carbon dioxide emissions in the energy sector and improve primary energy intensity [1].
In this context, Smart Energy Systems (SES) are acknowledged as a viable and effective solution to shift to a more sustainable energy paradigm [2]. Smart energy systems aim to optimally integrate different energy sectors (i.e., electricity, heating, cooling, industry, transportation), in order to improve the overall flexibility and efficiency. In this way, both the economic viability and environmental sustainability of energy systems can be enhanced.
Within this framework, Distributed Energy Systems (DES), which bring the production of energy close to end-users, are a clear opportunity to realize the smart integration of different energy technologies and sources [3]. Likewise, District Heating and Cooling (DHC) networks are key elements in achieving higher security of supply, lower costs, and lower carbon dioxide emissions [4]. Indeed, smart thermal grids can exploit hybrid substations (e.g., distributed heat pumps) and usually implement optimal control strategies, thus fostering the coupling between the electrical and thermal grids [5].
In recent years, the concept of 4th generation district heating (4GDH) has emerged, which outlines the main features that modern thermal grids must have to meet the challenges of a sustainable energy system [6]. The first and foremost feature is the low temperature of the water distribution, so that grid losses are lowered, heat from low-temperature and renewable sources can be exploited, and high-efficiency technology, such as heat pumps, can be integrated within the network. In this respect, several works can be found in recent literature. Sameti and Haghighat [7] proposed a mathematical programming model to optimize design and planning of a 4GDH, from an economic perspective. Cai et al. [8] demonstrated that lower supply temperatures can improve the energy efficiency of thermal networks, by analyzing the operation of electric heat boosters in low-energy districts. Likewise, a case study was investigated by Baldvinsson and Nakata [9], which proved that lower temperatures provide higher exergy efficiency in district heating systems.
At the same time, the issue of possible layouts of DES and plants feeding DHC networks has been thoroughly investigated in recent years. Indeed, the optimization problem of SES must deal with the synthesis of the plant, i.e., the existence and size of each component/technology considered in the superstructure of the energy system [10]. The possibility of integrating RES in district cooling systems has been exhaustively reviewed in [11]. Soltero et al. [12] analyzed biomass-based district heating systems, combining forestry biomass and agriculture residues as fuel. Testi et al. [13] proposed the adoption of distributed heat pumps to assist the thermal production at building level in a hybrid microgrid. The integration of a trigeneration system with an Organic Rankine Cycle and a ground source heat pump was studied by Kang et al. [14]. Possible configurations and operational modes of heat pumps in district heating systems have been proposed and analyzed by Ommen et al. [15] and by Sayegh et al. [16]. Similarly, Ostergaard et al. [17] investigated the effect of the introduction of heat pumps in small-scale district heating systems.
Basic trigeneration systems feeding DHC networks are traditionally composed of a prime mover (e.g., an internal combustion engine), a heat recovery system, and a thermally driven cooling technology (usually an absorption chiller). Moreover, back-up units are included in the system, to ensure that the demand can always be met and to improve the flexibility of the system. Generally, those back-up units are electrically driven chillers and fossil-fuel boilers, for the cooling and heating loads, respectively [18]. Nevertheless, fossil-fuel boilers are exergetically inefficient, especially when the temperature of the delivered heat is low, as in 4GDH.
On the other hand, electrically driven reversible heat pumps are a versatile technology, capable of providing both heating and cooling, alternately [19,20,21]. Furthermore, among other factors, the performance of heat pumps is highly dependent on the supply temperature; therefore, this technology is especially suitable for low-temperature heating applications. Moreover, unlike fossil-fuel boilers, electric heat pumps represent a hub between the electric and thermal networks, thus enhancing system flexibility and integration of RES.
In view of all this, this paper aims to investigate the integration of reversible electric heat pumps in trigeneration systems, as single back-up units to replace traditional fossil-fuel boilers and non-reversible electric chillers. An energy supply system for a district heating and cooling microgrid integrating RES is considered, consisting of an internal combustion engine, a single-effect absorption chiller, a reversible electric heat pump, photovoltaic panels, and a wind turbine. The thermal network serves office buildings equipped with fan coil terminals, thus requiring low-temperature heat. Energy demand data from a case study located in Pisa, Italy, are considered. The integrated optimal sizing and hourly operation of the proposed system are evaluated from the economic point of view. A two-level optimization algorithm is developed: the optimal operation problem is solved by exhaustive search, while a genetic algorithm is applied to the sizing problem. Moreover, a comparison with the conventional trigeneration system (i.e., with fossil-fuel boiler and non-reversible electric chiller) is also performed in terms of economic and environmental perspectives.

2. Methodology

2.1. The Energy System

As mentioned, the main purpose of this work is the comparison of two different configurations of trigeneration system integrating renewable energy sources. The schematics of the energy systems are shown in Figure 1. Both configurations comprise an internal combustion engine (ICE) as the cogeneration unit, an absorption chiller (AC), photovoltaic panels (PV), and a wind turbine (WT). The main difference between the configurations lies in the back-up units. In the traditional trigeneration system, a natural-gas boiler (B) and a non-reversible electric chiller (EC) assist the heating and cooling production, respectively, while, in the proposed configuration, a reversible heat pump (RHP) is included as single back-up unit for both thermal demands.

2.2. Subsystem Models

In this subsection, models and features adopted for the components of the simulated energy system are described.
The cogeneration unit is a spark ignition reciprocating internal combustion engine fueled by natural gas. Part-load efficiencies have been considered [22] and shown in Figure 2, where the load factor is defined as L I C E = E I C E / E I C E , n o m . The heat recovery system recovers heat from the exhaust gas, jacket cooling water system and lube oil of the ICE in the form of hot water at around 90 °C. The power capacity lower bound has been set to 50% of the nominal power capacity.
The reversible heat pump and the non-reversible electric chiller are air-to-water units. The single-effect absorption chiller, fed by the hot water recovered from the ICE, has a cooling tower as heat sink. The heat pump and the chillers have been modelled considering both ambient temperature and part-load condition effects on the coefficient of performance. The so-called “locally constant exergy efficiency” (or Carnot efficiency) approach allows to consider how the temperature of the thermal sources affects the performance of the heat pump/chiller. Moreover, part-load behavior can be considered by a correction factor. The coefficient of performance (COP) of the units is expressed as the product of the ideal COP (i.e., the COP of the Reversed Carnot Cycle), the exergy efficiency (defined as the ratio between the actual COP of the unit and the ideal COP at those thermal source temperatures), and a part-load correction factor [18]:
C O P = C O P i d · η e x · C F
The ideal COP of the electrically driven heat pumps under the heating and cooling modes are, respectively:
C O P R H P H , i d = T h o t + 273.15 T h o t T c o l d
C O P R H P C , i d / E C , i d = ( T c o l d + 273.15 ) T h o t T c o l d
The ideal COP of the absorption chiller is:
C O P A C , i d = 1 ( T h o t + 273.15 ) T g e n + 273.15 1 ( T c o l d + 273.15 ) T h o t + 273.15
Moreover, the effect of the ambient temperature on the capacity of the heat pump must be considered, by means of temperature coefficients ( κ ):
Q R H P H = Q n o m , R H P H 1 + κ R H P H T a m b T n o m , R H P H
C R H P C / E C = C n o m , R H P C / E C 1 + κ R H P C / E C T n o m , R H P C / E C T a m b
where T n o m , R H P H = 7   ° C , T n o m , R H P C / E C = 35   ° C , κ R H P H = 0.04   kW / ° C and κ R H P C / E C = 0.015 kW / ° C based on a survey on technical specification sheets by [23].
Based on technical specification sheets [24,25,26,27], the following values for the exergy efficiencies have been adopted: η e x , R H P H = 35 % , η e x , R H P   C / E C = 29 % , η E X , A C = 34 % . The part-load correction factors [28] are shown in Figure 3.
As mentioned, a cooling tower is included in the system to dissipate the heat from the condenser/absorber of the absorption chiller. The inlet temperature of the condenser water (outlet temperature of the cooling tower) changes as the outdoor weather conditions vary (i.e., as the ambient wet bulb temperature varies). The so called “fixed approach control” method has been considered as the control method of the cooling tower. Such a method consists in varying the cooling tower air flow rate to maintain a constant temperature difference between the cooling tower outlet water temperature and the ambient air wet bulb temperature [18]. This temperature difference is commonly known as “approach”. Therefore, the inlet temperature of the cooling water circuit of the absorption chiller can be evaluated as the sum of the ambient air wet bulb temperature and the approach temperature, set at 6 °C. The electric power consumption due to cooling tower fans is evaluated as 0.026 kW per kW of heat rejection capacity [29].
The natural-gas boiler has a nominal efficiency of 90% and it has been modelled considering part-load efficiency, as follows [30]:
η B = 0.9 0.0951 + 1.525 L B 0.6249 L B 2
The PV generator has been modeled by considering the ambient temperature and solar irradiance; the output power is evaluated as follows [31].
E P V = E P V , S T C G G S T C 1 γ T c e l l 25
where P P V , S T C is the rated power of the PV generator at standard test conditions (STC), G is the solar irradiance (kW/m2), γ is the power temperature coefficient, and T c e l l is the cell temperature. The STC are T c e l l , S T C = 25   ° C , G S T C = 1   kW / m 2 , and wind speed of 1 m/s. The cell temperature T c e l l is obtained as follows.
T c e l l = T a m b + G 0.8 N O C T 20
where T a m b is the ambient temperature and N O C T is the nominal operating cell temperature. In this paper, the PV generator is composed of polycrystalline modules that have all been considered to be at the same temperature. The temperature coefficient γ is 0.043 % / ° C , and the N O C T is 45.5   ° C [31].
The electrical power output of the horizontal-axis wind turbine E W T has been evaluated as follows [32].
E W T = 0   E W T = E n o m v k v c u t i n k / ( v n o m k v c u t i n k ) E W T = E n o m   E W T = 0     v v c u t i n v c u t i n v v n o m v n o m v v c u t o f f ( v > v c u t o f f )
where v is the wind speed and k is the Weibull shape parameter of the wind speed distribution [32]. The wind speed data, usually recorded at 10 m, must be corrected to consider the turbine hub height, as follows:
v = v r e f e r e n c e   h e i g h t log z t u r b i n e   h u b / z 0 log z r e f e r e n c e / z 0
where z 0 is the roughness length, which depends on the terrain characteristics.

2.3. Optimization Problem and Method

The objective function to be minimized is the Equivalent Annual Cost (EAC) of the system, which is calculated over the period of a year and composed of the annualized investment cost for the technologies, I n v e s t , and the total annual operating cost, O p .
E A C : f = I n v e s t + O p
I n v e s t = i a i C A P i b i C R F i
where C A P i is the capacity of the i   th technology to be installed (expressed in kW), a i and b i are the correlation parameters of the equipment cost as a function of the capacity, and C R F is the capital-recovery factor.
C R F = r r + 1 l i f e t i m e r + 1 l i f e t i m e 1
The total annual operating cost comprises the cost for purchasing electricity and natural gas and the revenue from selling electricity to the grid. An hourly timestep has been considered.
O p = t = 1 8760 c F F I C E , t + F B , t + c P E E P , t c S E E S , t
Demand constraints must be satisfied in each t   th time step.
E I C E , t + E P V , t + E W T , t + E P , t E S , t E E C / R H P C ,   t E R H P H , t E C T , t = E D , t
Q I C E , t + Q B / R H P H , t Q D , t
C A C , t + C E C / R H P C , t 0
where t = 1 , 2 , ,   8760 , and the electric consumption due to the cooling-tower fans depends on the amount of heat to be rejected from the absorption chiller:
E C T , t = C A C , t 1 + C O P A C , t C O P A C , t ε C T
Governing equations and capacity constraints must also be considered:
E I C E , t E I C E , n o m δ I C E , t
E I C E , t E I C E , m i n δ I C E , t
F I C E , t = E I C E , t / η E , I C E , t
Q I C E , t = E I C E , t η H , I C E , t η E , I C E , t
C A C , t = Q I C E , t / C O P A C , t
C A C , n o m = E I C E , n o m η H , I C E , n o m C O P A C , n o m
Q B , t Q B , n o m
F B , t Q B , t / η B , t
E E C , t E E C , n o m
E E C / R H P C ,   t = C E C / R H P C ,   t / C O P E C / R H P C ,   t
E R H P H , t = Q R H P H , t / C O P R H P H , t
Q R H P H , t C R H P C , t = 0
E R H P , t E R H P , n o m
C R H P C , n o m / C O P R H P C , n o m = Q R H P H , n o m / C O P R H P H , n o m
where t = 1 , 2 , ,   8760 .
Ultimately, the overall problem consists in the minimization of the Equivalent Annual Cost:
m i n i m i z e f = E A C
which results in a mixed integer non-linear optimization problem.
To solve this optimization problem, a two-stage algorithm has been specifically developed, which is schematically summarized in Figure 4. The optimization variables can be distinguished in two main groups: sizing variables (i.e., E I C E , n o m ,   Q B , n o m , C E C , n o m , Q R H P H , n o m , E P V , E W T ) and operation variables (i.e., E I C E , t , C E C / R H P C , t , Q B / R H P H , t , E S , t , E P , t , C A C , t with t = 1 , , 8760 ). Therefore, the problem can be decomposed into two levels: the upper level, where the optimal capacity of each energy unit is identified, and the lower level, where the optimal operating conditions of the system in each timestep are found. Obviously, the two problems must be solved simultaneously.
At the upper level, the sizing problem is addressed by means of a Genetic Algorithm, which is a widely used optimization tool for the design of DESs [33]. Clearly, it is also possible that in the optimal solution a technology is not included in the system, i.e., its resulting optimal size is null. At the lower level, since the operation problem under analysis is static, it can be solved by means of a “greedy approach” [18,34]. That means the overall operation problem can be split into 8760 sub-problems, one for each timestep, to be solved independently. Then, for each timestep, the problem is reduced to find the optimal load factor of the ICE. Indeed, once the load factor of the ICE is set, the ICE heating production and the cooling production by the absorption chiller are known. Therefore, the back-up units (either the boiler and the non-reversible electric chiller, or the reversible heat pump) just need to meet the remaining thermal demand, if any. The electric balance is also trivial, once the load factor of the ICE is set, since the production by RES in each timestep is also known. Consequently, the optimal operation problem consists in finding the optimal load factor of the ICE, in each timestep. This task is fulfilled with an exhaustive-search, which guarantees the global optimality of the solution.
All the simulations and optimization are performed using scripts written in MATLAB environment. The MATLAB Genetic Algorithm Solver [35] is used, whose settings and stopping criteria are shown in Table 1. A simple in-house code is adopted to implement the greedy approach, and a 0.05 step in the load factor of the ICE is considered in the exhaustive search.

3. Case Study: Energy Demand and Economic Parameters

Both energy system configurations under analysis have been tested in a case study. Energy demand data of a large office building located in Pisa (Italy) have been considered. The office is located in an 8-story building, with a gross volume about 73,000 m3 and a floor area about 24,000 m2, having 25-cm concrete walls and double-glazed windows (the window-to-wall ratio is 0.16). DHW demand is met by small local electrically driven heaters, so the plant must satisfy only space heating and cooling requirements. The complete low-temperature district heating and cooling network is not modelled. Just to test the viability of the energy supply plant using reversible heat pumps, the energy demand is simplified to cover the demand of an office building that could be representative of the energy profiles expected to be covered in urban DHC networks. The proposed methodology and system configuration could be easily applied to a bigger-scale system. Figure 5 shows the hourly energy demand data for the case study. Three energy services are needed: electricity, cooling (chilled water from 7 to 12 °C) and heating (hot water from 45 to 35 °C). Hourly data of electricity consumption were available from measurements, while thermal load demands were evaluated by a dynamic simulation in TRNSYS 17. The hourly profiles of the external temperature, relative humidity, solar irradiance, and wind speed of the typical meteorological year [36] of Pisa have been adopted.
Table 2 shows the values of all the economic parameters adopted in the simulations. The considered lifetime for all technologies is 20 years [13].

4. Results

In this section, the results of the integrated optimal sizing and operation in terms of Equivalent Annual Cost are shown and discussed. As mentioned, the optimization has been performed for both the energy system configurations shown in Figure 1, i.e., the proposed configuration with a reversible heat pump as single thermal back-up unit, and the traditional configuration with a boiler and a non-reversible electric chiller as thermal back-up units. Table 3 summarizes the main results obtained from the optimizations. In the evaluation of the carbon dioxide emissions, the values of the carbon intensity of the electricity from the grid and the natural gas are 325.2 and 202.0 gCO2eq/kWh, respectively [42].
First and foremost, the proposed configuration with reversible heat pump achieves the lowest EAC, with a 7% saving compared to the traditional configuration. The discounted payback period, considering the initial extra investment and annual cost savings of the proposed configuration compared to the traditional one, is around 2.5 years. Moreover, the installed capacity of RES is 23% higher in the proposed configuration. This is due to the installation of the reversible heat pump; indeed, the larger share of cooling covered by the RHP compared to the non-reversible electric chiller allows a higher self-consumption of the electricity produced by RES. Both the optimal design solutions do not comprise the wind turbine, because of the low wind speed at the site. For this reason, all the installed RES capacity is represented exclusively by the PV system. Another significant effect of the integration of the reversible heat pump is the low value of the annual carbon dioxide emissions, which is around 11% lower than that of the traditional configuration. Indeed, the use of natural gas is reduced by 18% with the integration of the reversible heat pump, as a result of the elimination of the boiler and the 34% downsizing of the internal combustion engine.
Moreover, according to the European Commission Decision establishing the guidelines for Member States on calculating renewable energy from heat pumps [43], part of the energy produced by the heat pump may be considered as renewable energy if the seasonal coefficient of performance (SCOP) is higher than 2.5. In the case study, the SCOP of the reversible heat pump is equal to 3.3. The amount of renewable energy supplied by the reversible heat pump, Θ R E S , R H P , shall be calculated with the following formula:
Θ R E S , R H P = Θ R H P 1 1 S C O P R H P
where Θ R H P is the total thermal energy delivered by the heat pump, both heat and cooling. As a result, around 70% of the energy produced by the reversible heat pump can be considered renewable, for a total of 570 MWh. Instead, in the traditional configuration, the electric chiller can be considered to produce only 190 MWh of renewable energy.
Figure 6 and Figure 7 summarizes the annual share of production and consumption of the different energy services (i.e., electricity, heating, and cooling) for the two configurations under analysis. The share of electricity produced by the PV panels increases from the 31.8% to 36.6% when the reversible heat pump is included. Moreover, the amount of electricity sold to the grid remains approximately the same, despite the increase of RES installed capacity. Overall, the reversible heat pump plays a key role in covering the heating and cooling demands, with an increase of around 18% compared to the share covered by the boiler and the non-reversible electric chiller in the traditional configuration.
Finally, Figure 8, Figure 9, Figure 10 and Figure 11 show how the energy demands (electricity, heating, and cooling) are met, in typical winter and summer weeks, with the novel configuration of the energy system. The internal combustion engine is generally continuously on to cover the electric base-load, except when the PV production is particularly high. To meet peak-demands, electricity is bought from the grid, while the sale of electricity to the grid only occurs with high PV production. The heating and cooling base-loads are covered by the ICE and the absorption chiller, respectively, while peak-demands are covered by the reversible heat pump during both winter and summer.
Obviously, the performance of the proposed configuration is strongly dependent on weather conditions and demand profiles of the case study. More specifically, a sufficiently high RES availability is needed, so that investment on PV panels or wind turbine is justified. Moreover, both peak and average heating and cooling demands must be comparable, so that the capacity of the reversible heat pump is suitable for both energy services. In summary, the proposed configuration is especially suitable for climates where the heating and cooling seasons are comparable in terms of duration and intensity and the RES electricity production is significant.

5. Conclusions

The integration of reversible heat pumps in hybrid renewable trigeneration systems serving low-energy district heating and cooling microgrids was proposed and investigated in this paper. Traditional thermal back-up units, i.e., natural gas boiler and non-reversible electric chiller, have been replaced by an electrically driven reversible heat pump as single thermal back-up unit. The energy system also includes an internal combustion engine, an absorption chiller, photovoltaic panels, and a wind turbine.
The proposed configuration was compared to the benchmark configuration of hybrid trigeneration system with conventional back-up units. A case study was considered, and the integrated optimal sizing and operation were identified by means of a genetic algorithm-based procedure. The proposed configuration achieved the best economic performance, with a 7% saving in the equivalent annual cost. Moreover, with the integration of the reversible heat pump, a 23% higher capacity of renewables was installed in the optimal design of the system. As a result, annual carbon dioxide emissions were lowered by 11% compared to the traditional configuration.
Overall, the integration of electrically driven reversible heat pumps in trigeneration systems was found to be an economically viable solution capable of increasing the share of renewables, reducing the use of fossil fuels, and lowering greenhouse gas emissions.

Author Contributions

Conceptualization: L.U., D.T. and J.C.B.; Methodology: L.U. and D.T.; Software: L.U.; Data curation: L.U.; Writing-original draft: L.U.; Writing-review & editing: L.U., D.T. and J.C.B.; Supervision: D.T. and J.C.B.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

Acronyms
4GDH4th Generation District Heating
DHCDistrict Heating and Cooling
GAGenetic Algorithm
NOCTNominal Operating Cell Temperature
RESRenewable Energy Source
STCStandard Test Conditions
Parameters
aCoefficient of the equipment cost curve, €/kW
bExponent of the equipment cost curve, dimensionless
cCost per unit of energy, €/kWh
CAPCapacity, kW
CFCorrection factor, dimensionless
COPCoefficient of performance, dimensionless
CRFCapital recovery factor, dimensionless
EACEquivalent annual cost, €
InvestInvestment cost, €
kWeibull shape parameter, dimensionless
OpAnnual operating cost, €
rInterest rate, dimensionless
SCOPSeasonal coefficient of performance, dimensionless
tTimestep, hour
zHeight, m
γ PV power temperature coefficient, °C−1
ε Electric energy consumption per unit of rejected heat, dimensionless
η Efficiency, dimensionless
κ Heat pump capacity temperature coefficient, kW/°C
Continuous variables
CCooling power, kW
EElectric power, kW
GSolar irradiance, kW/m2
LLoad factor (defined as the ratio between the power output of the unit and its capacity), dimensionless
vWind speed, m/s
QHeating power, kW
TTemperature, °C
Θ Thermal energy, kWh
Binary variables
δ On-off state
Subscripts
ACAbsorption chiller
BBoiler
CCooling mode
CTCooling Tower
DDemand
EElectric energy
ECElectric chiller
exExergy
FFuel
genGenerator
HHeating mode
ICEInternal Combustion Engine
nomNominal
PPurchased from the grid
PVPhotoVoltaics
RHPReversible Heat Pump
SSold to the grid
WTWind Turbine

References

  1. International Energy Agency (IEA). Global Energy and CO2 Status Report 2018. Available online: https://www.iea.org/geco (accessed on 5 August 2019).
  2. Lund, H.; Østergaard, P.A.; Connolly, D.; Mathiesen, B.V. Smart energy and smart energy systems. Energy 2017, 137, 556–565. [Google Scholar] [CrossRef]
  3. Karmellos, M.; Georgiou, P.N.; Mavrotas, G. A comparison of methods for the optimal design of Dsitrbuted Energy Systems under uncertainty. Energy 2019, 178, 318–333. [Google Scholar] [CrossRef]
  4. Werner, S. International review of district heating and cooling. Energy 2017, 137, 617–631. [Google Scholar] [CrossRef]
  5. Buffa, S.; Cozzini, M.; D’Antoni, M.; Baratieri, M.; Fedrizzi, R. 5th generation district heating and cooling systems: A review of existing cases in Europe. Renew. Sustain. Energy Rev. 2019, 104, 504–522. [Google Scholar] [CrossRef]
  6. Lund, H.; Werner, S.; Wiltshire, R.; Svendsen, S.; Thorsen, J.E.; Hvelpund, F.; Vad Mathiesen, B. 4th Generation District Heating (4GDH) Integrating smart thermal grids into future sustainable energy systems. Energy 2014, 68, 1–11. [Google Scholar] [CrossRef]
  7. Sameti, M.; Haghighat, F. Optimization of 4th generation distributed district heating system: Design and planning of combined heat and power. Renew. Energy 2019, 130, 371–387. [Google Scholar] [CrossRef]
  8. Cai, H.; You, S.; Wang, J.; Bindner, H.W.; Klyapovskiy, S. Technical assessment of electric heat boosters in low-temperature district heating based on combined heat and power analysis. Energy 2018, 150, 938–949. [Google Scholar] [CrossRef] [Green Version]
  9. Baldvinsson, I.; Nakata, T. A feasibility and performance assessment of a low temperature district heating system–a North Japanese case study. Energy 2016, 95, 155–174. [Google Scholar] [CrossRef]
  10. Sameti, M.; Haghighat, F. Optimization approaches in district heating and cooling thermal network. Energy Build. 2017, 140, 121–130. [Google Scholar] [CrossRef]
  11. Inayat, A.; Raza, M. District cooling system via renewable energy sources: A review. Renew. Sustain. Energy Rev. 2019, 107, 360–373. [Google Scholar] [CrossRef]
  12. Soltero, V.M.; Chacartegui, R.; Ortiz, C.; Lizana, J.; Quirosa, G. Biomass District Heating Systems Based on Agriculture Residues. Appl. Sci. 2018, 8, 476. [Google Scholar] [CrossRef]
  13. Testi, D.; Conti, P.; Schito, E.; Urbanucci, L.; D’Ettorre, F. Synthesis and Optimal Operation of Smart Microgrids Serving a Cluster of Buildings on a Campus with Centralized and Distributed Hybrid Renewable Energy Units. Energies 2019, 12, 745. [Google Scholar] [CrossRef]
  14. Kang, L.G.; Yang, J.H.; An, Q.S.; Deng, S.; Zhao, J.; Li, Z.L.; Wang, Y.Z. Complementary configuration and performance comparison of CCHP-ORC system with a ground source heat pump under different operation strategies. Energy Convers. Manag. 2017, 135, 244–255. [Google Scholar] [CrossRef]
  15. Ommen, T.S.; Markussen, W.B.; Elmegaard, B. Heat pumps in combined heat and power systems. Energy 2014, 76, 989–1000. [Google Scholar] [CrossRef]
  16. Sayegh, M.; Jadwiszczak, P.; Axcell, B.; Niemierka, E.; Brys, K.; Jouhara, H. Heat pump placement, connection and operational modes in European district heating. Energy Build. 2018, 166, 122–144. [Google Scholar] [CrossRef]
  17. Ostergaard, P.A.; Jantzen, J.; Marczinkowski, H.M.; Kristensen, M. Business and socioeconomic assessment of introducing heat pumps with heat storage in small-scale district heating systems. Renew. Energy 2019, 139, 904–914. [Google Scholar] [CrossRef]
  18. Urbanucci, L.; Testi, D.; Bruno, J. An operational optimization method for a complex polygeneration plant based on real-time measurements. Energy Convers. Manag. 2018, 170, 50–61. [Google Scholar] [CrossRef]
  19. Kinab, E.; Marchio, D.; Rivière, P.; Zoughaib, A. Reversible heat pump model for seasonal performance optimization. Energy Build. 2010, 42, 2269–2280. [Google Scholar] [CrossRef]
  20. Madonna, F.; Bazzocchi, F. Annual performances of reversible air-to-water heat pumps in small residential buildings. Energy Build. 2013, 65, 299–309. [Google Scholar] [CrossRef]
  21. Dongellini, M.; Naldi, C.; Morini, G.L. Sizing effects on the energy performance of reversible air-source heat pumps for office buildings. Appl. Therm. Eng. 2017, 114, 1073–1081. [Google Scholar] [CrossRef]
  22. Ebrahimi, M.; Keshavarz, A. Combined Cooling, Heating and Power: Decision-Making, Design and Optimization; Elsevier: Amsterdam, The Netherlands, 2015. [Google Scholar]
  23. Cinquini, F. Dimensionamento e Gestione Ottimizzati di Sistemi a Pompa di Calore Aria-Acqua per la Climatizzazione Annuale di Edifici Terziari. Master’s Thesis, University of Pisa, Pisa, Italy, 2018. [Google Scholar]
  24. Aermec. Heat Pumps and Chillers: Catalogue and Technical Documentation. Available online: https://support.aermec.com/default.asp (accessed on 5 August 2019).
  25. Galletti. Heat Pumps and Chillers: Catalogue and Technical Documentation. Available online: http://www.galletti.com/users/login-utente (accessed on 5 August 2019).
  26. Rossato. Heat Pumps and Chillers: Catalogue and Technical Documentation. Available online: https://www.rossatogroup.com/prodotti/pompe-di-calore/pompe-di-calore-aria-acqua.html (accessed on 5 August 2019).
  27. YORK. Heat Pumps and Chillers: Catalogue and Technical Documentation. Available online: https://www.johnsoncontrols.com/-/media/jci/global-capabilities/be/files/be_york_industrial_commercial_hvac_2018.pdf (accessed on 5 August 2019).
  28. European Committee for Standardization. Heating Systems in Buildings–Method for Calculation of System Energy Requirements and System Efficiencies—Part 4-2: Space Heating Generation Systems, Heat Pump Systems; EN 15316-4-2; CEN: Brussels, Belgium, 2017. [Google Scholar]
  29. DiPippo, R. Geothermal Power Generation: Developments and Innovation; Woodhead Publishing: Sawston/Cambrige, UK, 2016. [Google Scholar]
  30. Zhou, Z.; Liu, P.; Li, Z.; Pistikopoulos, E.N.; Georgiadis, M.C. Impacts of equipment off-design characteristics on the optimal design and operation of combined cooling, heating and power systems. Comput. Chem. Eng. 2013, 48, 40–47. [Google Scholar] [CrossRef]
  31. Riffonneau, Y.; Bacha, S.; Barruel, F.; Ploix, S. Optimal Power Flow Management for Grid Connected PV Systems with Batteries. IEEE Trans. Sustain. Energy 2011, 2, 309–320. [Google Scholar] [CrossRef]
  32. Borowy, B.S.; Salameh, Z.M. Optimum photovoltaic array size for a hybrid wind/PV system. IEEE Trans. Energy Convers. 1994, 9, 482–488. [Google Scholar] [CrossRef]
  33. Abdmouleh, Z.; Gastli, A.; Ben-Brahim, L.; Haouari, M.; Al-Emadi, N.A. Review of optimization techniques applied for the integration of distributed generation from renewable energy sources. Renew. Energy 2017, 113, 266–280. [Google Scholar] [CrossRef]
  34. Urbanucci, L.; Testi, D. Optimal integrated sizing and operation of a CHP system with Monte Carlo risk analysis for long-term uncertainty in energy demands. Energy Convers. Manag. 2018, 157, 307–316. [Google Scholar] [CrossRef]
  35. MATLAB. Genetic Algorithm Solver. Available online: https://it.mathworks.com/discovery/genetic-algorithm.html (accessed on 5 August 2019).
  36. Italian Thermotechnical Committee (CTI). National Typical Meteorological Years; CTI: Milan, Italy, 2016. [Google Scholar]
  37. Carvalho, M.; Lozano, M.A.; Ramos, J.; Serra, L.M. Synthesis of Trigeneration Systems: Sensitivity Analyses and Resilience. Sci. World J. 2013, 2013, 1–16. [Google Scholar] [CrossRef] [Green Version]
  38. Urbanucci, L.; D’Ettorre, F.; Testi, D. A Comprehensive Methodology for the Integrated Optimal Sizing and Operation of Cogeneration Systems with Thermal Energy Storage. Energies 2019, 12, 875. [Google Scholar] [CrossRef]
  39. Renewable Energy Report, Politecnico di Milano. 2018. Available online: http://www.energystrategy.it/report/renewable-energy-report.html (accessed on 5 August 2019).
  40. Urbanucci, L.; Bruno, J.C.; Testi, D. Thermodynamic and economic analysis of the integration of high-temperature heat pumps in trigeneration systems. Appl. Energy 2019, 238, 516–533. [Google Scholar] [CrossRef]
  41. Eurostat, European Commission. Energy Statistics–Prices of Natural Gas and Electricity. Available online: https://ec.europa.eu/eurostat/web/energy/data/database (accessed on 5 August 2019).
  42. Energy efficiency report 2018, Politecnico di Milano. Available online: http://www.energystrategy.it/assets/files/EER_18_web.pdf (accessed on 5 August 2019).
  43. Commission Decision of 1 March 2013 establishing the guidelines for Member States on calculating renewable energy from heat pumps from different heat pump technologies pursuant to Article 5 of Directive 2009/28/EC of the European Parliament and of the Council. Official Journal of the European Union 2013, 62, 27, Decision 2013/114/EU. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32013D0114&from=EN (accessed on 5 August 2019).
Figure 1. Schematics of the compared energy systems.
Figure 1. Schematics of the compared energy systems.
Applsci 09 03194 g001
Figure 2. Part-load efficiencies of the ICE.
Figure 2. Part-load efficiencies of the ICE.
Applsci 09 03194 g002
Figure 3. Part-load correction factors for the absorption chiller, reversible heat pump, and electric chiller.
Figure 3. Part-load correction factors for the absorption chiller, reversible heat pump, and electric chiller.
Applsci 09 03194 g003
Figure 4. Optimization procedure.
Figure 4. Optimization procedure.
Applsci 09 03194 g004
Figure 5. Hourly energy demand data of the case study.
Figure 5. Hourly energy demand data of the case study.
Applsci 09 03194 g005
Figure 6. Sankey diagrams of the annual electricity balance for (a) traditional and (b) novel configurations.
Figure 6. Sankey diagrams of the annual electricity balance for (a) traditional and (b) novel configurations.
Applsci 09 03194 g006
Figure 7. Sankey diagrams of the annual thermal energy balances for (a) traditional and (b) novel configurations.
Figure 7. Sankey diagrams of the annual thermal energy balances for (a) traditional and (b) novel configurations.
Applsci 09 03194 g007
Figure 8. Optimal operation: electricity flows in a typical winter week.
Figure 8. Optimal operation: electricity flows in a typical winter week.
Applsci 09 03194 g008
Figure 9. Optimal operation: heating demand and load shares in a typical week.
Figure 9. Optimal operation: heating demand and load shares in a typical week.
Applsci 09 03194 g009
Figure 10. Optimal operation: electricity flows in a typical summer week.
Figure 10. Optimal operation: electricity flows in a typical summer week.
Applsci 09 03194 g010
Figure 11. Optimal operation: cooling demand and load shares in a typical week.
Figure 11. Optimal operation: cooling demand and load shares in a typical week.
Applsci 09 03194 g011
Table 1. Settings and parameters adopted for the Genetic Algorithm.
Table 1. Settings and parameters adopted for the Genetic Algorithm.
P o p u l a t i o n S i z e   E l i t e C o u n t     T o l F u n     C r o s s o v e r F r a c t i o n     M i g r a t i o n D i r e c t i o n  
  250     10     1 e 6     0.8     f o r w a r d  
  S t a l l G e n L i m i t     G e n e r a t i o n s     T o l C o n     M i g r a t i o n F r a c t i o n     M i g r a t i o n I n t e r v a l  
  50     400     1 e 3     0.2     20  
Table 2. Price of energy carriers and investment cost parameters (used in Equations (13)–(15)).
Table 2. Price of energy carriers and investment cost parameters (used in Equations (13)–(15)).
ParameterValueParameterValueParameterValueParameterValue
  a I C E     7789 [13]   b C T     1 [37]   a B     345.9 [38]   b W T     1 [39]
  b I C E     0.6 [13]   a R H P     1036 [23]   b B     0.79 [38]   r     7 % [13]
  a A C     375 [40]   b R H P     0.74 [23]   a P V     1500 [39]   c F     0.052   / kWh [41]
  b A C     0.65 [40]   a E C     609 [23]   b P V     1 [39]   c P E     0.166   / kWh [41]
  a C T     28.75 [40]   b E C     0.79 [23]   a W T     1570 [39]   c S E     0.040   / kWh [13]
Table 3. Optimal Equivalent Annual Cost: sizing and operation—results.
Table 3. Optimal Equivalent Annual Cost: sizing and operation—results.
Energy System ConfigurationNovel Configuration with Reversible Heat PumpTraditional Configuration with Boiler and Non-Reversible Electric Chiller
Internal Combustion Engine (Nominal Electric Power)165 kW250 kW
Absorption Chiller
(Cooling Capacity)
190 kW290 kW
Reversible Heat Pump (Heating/Cooling Capacity)795 kW/660 kW\
Boiler (Heating Capacity)\615 kW
Electric Chiller
(Cooling Capacity)
\575 kW
Photovoltaic Panels
(Power at STC)
505 kW410 kW
Wind Turbine
(Nominal Power)
0 kW0 kW
Annual Operating Cost175.6 k€203.3 k€
Investment Cost1199.0 k€1138.0 k€
Equivalent Annual Cost288.8 k€310.8 k€
Annual Carbon
Dioxide Emissions
580.7 ton694.9 ton

Share and Cite

MDPI and ACS Style

Urbanucci, L.; Testi, D.; Bruno, J.C. Integration of Reversible Heat Pumps in Trigeneration Systems for Low-Temperature Renewable District Heating and Cooling Microgrids. Appl. Sci. 2019, 9, 3194. https://0-doi-org.brum.beds.ac.uk/10.3390/app9153194

AMA Style

Urbanucci L, Testi D, Bruno JC. Integration of Reversible Heat Pumps in Trigeneration Systems for Low-Temperature Renewable District Heating and Cooling Microgrids. Applied Sciences. 2019; 9(15):3194. https://0-doi-org.brum.beds.ac.uk/10.3390/app9153194

Chicago/Turabian Style

Urbanucci, Luca, Daniele Testi, and Joan Carles Bruno. 2019. "Integration of Reversible Heat Pumps in Trigeneration Systems for Low-Temperature Renewable District Heating and Cooling Microgrids" Applied Sciences 9, no. 15: 3194. https://0-doi-org.brum.beds.ac.uk/10.3390/app9153194

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop