 ORIGINAL ARTICLE
 Open Access
 Published:
Upscaling a district heating system based on biogas cogeneration and heat pumps
Energy, Sustainability and Society volumeÂ 5, ArticleÂ number:Â 16 (2015)
Abstract
Background
The energy supply of the Meppel district Nieuwveense landen is based on biogas cogeneration, district heating, and ground source heat pumps. A centrally located combined heat and power engine (CHP) converts biogas from the municipal wastewater treatment facility into electricity for heat pumps and heat for district heating purposes. Development of the urban district is influenced by the current economic and building decline. For the district heating energy concept, a migration strategy for the required infrastructure is required. The migration spans the districtâ€™s smallscale starting phase involving 40 houses up to a scale of 176 houses.
Methods
An optimization model which maximizes profitability is developed which includes data from district heating and cooling demand patterns.
Results
With the optimization model, optimal CHP size, boiler size, and operational hours are determined for various scenarios.
Conclusions
From the scenario analysis, a migration strategy is developed which starts with a simple system concept supported by boilers to a larger system which includes a CHP. Sustainability in terms of CO _{2} emission savings of the energy concept is compared with other possible energy concepts.
Background
In the Netherlands, at least 40% of the national carbon dioxide emissions are related to the heating demand of buildings. Due to the natural gas reserves, conversion by natural gas boilers dominates the heating supply for space heating and hot tap water. But natural gas reserves are gradually depleting and largescale carbon dioxide emissions cause climate change. Hence, Dutch policies are aimed to gradually diminish the consumption of fossil fuels in general and increase the share of renewable energy. Recently, local governments take up their responsibility and create favorable policies to stimulate reduction of energy consumption and increase the share of renewable energy on a local scale. One of these initiatives is the smart grid demonstration project MeppelEnergie, which is funded by the Dutch program Switch2SmartGrids. For the projectâ€™s energy system, we develop a smart grid control. The Meppel energy concept consists of a biogas combined heat and power engine (CHP), backup boilers, high temperature (HT) water storage, heat pumps, and aquifer underground thermal storage. The CHP generates electric and thermal energy. The thermal energy of the CHP is used for district heating, the electricity is used to supply heat pumps placed at houses with no connection to the district heating, or is sold through the external grid.
Cooling energy for the houses is provided by an underground aquifer consisting of a warm and cold well. During the heating season, the warm well provides low temperature (LT) heat for the heat pumps. Cooling energy of the houses also provides part of the required regeneration energy to maintain temperature balance within the underground thermal storage. Another part of the regeneration is provided either by a dry cooler or an effluent stream from the municipal wastewater treatment plant (MTP).
Due to the recent economic recession in the Dutch housing market, the planning of the districtâ€™s expansion is shifted from ambitious to a more realistic planning. The district is now built in smaller steps, i.e., phase 1 for the period up to 2015 contains 40 houses. After that, phase 2 contains another 136 houses and is built up to 2017. Even after phase 2, the scale is not yet reached to justify all investments in the pictured energy concept. Hence, the energy concept is to be treated as a horizon relating to 300 to 400 houses in the district, which is not expected before 2020. We investigate in this paper a possible scenario for the evolution of the district heating system from a starting system with less investments for the first building phase as shown in Figure 1 up to the final system which includes aquifer thermal storage. The starting energy concept consists of heat generators which run on natural gas. It has to be determined which type of generators (boilers or a CHP with supporting boilers) are optimal in this case. Without biogas supply during the first years, integration of renewable energy is realized by purchasing certified green energy (green gas) from the grid. The starting concept also includes a refrigeration chiller to provide cooling for the houses.
The first goal of this paper is to determine the required heating and cooling generator capacity rates for a district heating and cooling system. The second goal is to determine the most economic heat converter (boiler or CHP) in relation to the number of houses. The third goal is to compare the sustainability performance of the starting energy concept with other options. Contributions of this paper are twofold, first we develop a simplified methodology to estimate heating and cooling demand data for a district heating network and a method to determine the operating time of a heating or cooling device based on such data. Second, we demonstrate the application of such methods for an integrative study on financial and sustainability aspects of a complex energy system.
The paper is structured as follows: the â€˜Backgroundâ€™ section introduces the energy concept and problem statement. The â€˜Related workâ€™ subsection outlines work related to our investigation. The â€˜Methodsâ€™ section defines the energy model, equations and case related parameters. The â€˜Results and discussionâ€™ section includes results, and finally, in the â€˜Conclusionsâ€™ section, the conclusions are drawn.
Related work
Experience with district heating in the Netherlands is traditionally related to largescale steam power plants (e.g., Amsterdam and Almere city heating network and many smaller city projects). More recently, decentral projects for new urban districts apply either biomass (wood) thermal conversion in boilers (e.g., muziekwijk Zwolle) or biogas cogeneration (e.g., Apeldoorn, Zeewolde, and Leeuwarden) [1]. An overview of Dutch district heating projects and profitability investigation is presented in [2].
Vallios et al. [3] develop a modeling approach for sizing biomassfueled district heating systems. The district heating (DH) system only includes biomass boilers. Optimal size of a residential micro CHP and storage in Japan is investigated by Ren et al. [4] who apply a mixed integer nonlinear programming model. They determine optimal CHP size in relation to varying electricity prices. In a similar study, Cho et al. [5] determine with a linear programming model, for a micro CHP an optimal control scheme to reach economic optimal operation. Chinese et al. [6] develop two optimization models to investigate the design of biomass based district heating networks in Italy in relation to profit maximization and greenhouse gas minimization. The optimization is applied to the design of a new industrial district heating network. Curti et al. [7] investigate configurations for district heating heat pumps which draw low temperature heat from Swiss lakes, supported by cogeneration. They use genetic algorithms to find optimal system sizes. Hlebnikov et al. [8] investigate different scenarios to renovate an existing Estonian district heating network and they apply EnergyPLAN and RETScreen softwares. Lund et al. [9] investigate the role of district heating for the future energy system of Denmark, developing a model of the Danish energy system in EnergyPLAN. Pirouti [10] develops a nonlinear programming model to design district heating networks which are minimized on total energy and exergy consumption and operational costs. Finally, Ã˜stergaard [11] reports on the regional energy system of the Danish city Frederikshavn in which an EnergyPLAN model of the regionâ€™s energy production and consumption is used to predict the outcome of geothermal energy application.
Most authors develop specific models for optimization of case specific energy systems, while others apply modeling capabilities of specific software like EnergyPLAN or RETScreen. Both approaches require specific input generation, e.g., home heating and cooling demand, energy price, and tax schemes. Besides that, constraints of the energy supply system are required, e.g., amount of available biogas, supplying order between the CHP and boilers, and system efficiencies. We consider using EnergyPLAN at a later stage of our research. Due to the complexity involved, we develop our own optimization model within this paper which at a later stage serves as validation tool for results generated by EnergyPLAN.
An important difference between the optimization carried out in, e.g., the papers of Vallios et al., Ren et al., Chinese et al., and Curti et al. and our optimization is that we have limited the amount of variables by developing an additional constraint which relates size of each generator (CHP, boiler) to its operational time. Another difference is that within the present paper, our focus is not limited to technical size but also to compare energy concepts on profit and sustainability performance.
Methods
District heating system energy model
The starting concept of the DH system is characterized by converter efficiencies, energy sources, energy demands, and energy flows, as indicated in Figure 1. In the following, we develop an energy balance model for which the nomenclature is given in Table 1.
Efficiency of the thermal storage is determined by the size, the wall insulation, and the average temperature difference between water within the storage and outside air. For the model, we assume that the size of the storages is sufficient to store and supply thermal energy within a daily cycle. As a rule of thumb, 125 L/household storage capacity is estimated for a DH system with maximum storage temperature of 90Â°C. With this, we determine an average storage efficiency. Hence, capacity or size of the storage is not a variable within the model.
The DH system is conceived as a utility which buys energy (biogas, natural gas) from the upstream market, converts it to electricity and heat, and sells heat to the downstream market, i.e., the households. The electricity is used by circulation pumps and by heat pumps. Surplus electricity produced by the CHP and not used by the local energy system is sold to the grid. Our first emphasis is on the primary choice between either a DH system with only natural gas boilers or with a biogas CHP and some supporting natural gas boilers.
Besides differences in capital and operational costs of the generators, we have to take cost differences of the downstream system (i.e., storage and network assets) into account caused by the type of converter being applied. The network assets are not influenced but we expect a larger HT storage in case of a CHP. When only boilers are applied, it is easy to adjust the generated thermal power instantly within a large range due to cascading and power modulation. With a CHP the thermal modulation range is not as large and partload conditions lead to a loss of electrical generation efficiency. Besides that, for best efficiency and lifetime performance, a CHP should run more or less continuous when it runs and should not constantly follow the actual demand of the network. Therefore, increased HT storage capital costs are calculated for the CHP case.
Based on average climatic data, the model determines the base generator sizes for heating (CHP, boilers). How this is done is explained in the â€˜Relation between CHP size and operational timeâ€™ subsection. In practice, an additional peak boiler is added to provide sufficient heating capacity during extremely cold days which occur infrequently. The peak boiler capacity is calculated from the coldest day during the last 5 years. For cooling, we determine the generator size directly from the warmest day during the last 5 years.
We assume energy tariffs (natural gas, electricity) are only related to total yearly amounts, although on the electricity spot market and day ahead market, prices vary hourly. But these variations are small and often fixed electricity prices are contracted for longer periods. In this way, the complexity and the number of required iterations for the optimization problem to find the number of houses for break even profit is significantly reduced.
The yearly heat balance for the converters is:
When the CHP and the boilers together generate the required heat, the CHP has priority and the boilers have a supportive function. Hence, an additional expression is required which relates the operational hours of the CHP to the size of the CHP, which is developed in the â€˜Relation between CHP size and operational timeâ€™ subsection.
The total yearly heat demand of the DH system (Q _{dh,demand,h }) is related to the average household heat demand. For the average household demand, suitable expressions are developed in the â€˜Thermal demand specificationâ€™ subsection. The following relation applies for n houses:
Equally, for m houses heated with heat pumps the following relation applies:
The source heat demand of the heat pumps is related to the household demand by:
CHP heat generation rate is related to CHP electric power rate as follows:
The biogas and natural gas flows are related to the CHP electric and boiler thermal power rate, respectively, as follows:
The total yearly electric energy demand including generation of cooling energy is as follows:
f _{netw,p } is an estimated average percentage which discounts the required electric pump energy for network pumps as a percentage of the thermal energy transmitted through the piping network. Real pump energies are to be estimated from network fluid friction calculations and frequency controlled pump characteristics. A value of 5% is derived from practical calculations.
Optimization model equations
Financial profit of the starting energy concept shown in Figure 1 is analyzed with an equation in which profit terms are considered related to selling energy to houses (heating and cooling) and to the network (surplus CHP electricity). Cost terms that we take into account are as follows: the input energy costs (fuel and bought electricity) and the true capital and operational costs of specific energy generators (boilers, CHP, heat pumps, and related storages). In Table 2, the nomenclature of the profit function is given. We excluded infrastructural costs of the downstream network because in practice, all infrastructural investments are payed for by homeowners through connection fees, so these are not carried by the utility. However, there are also financial reservations to be made for reinvestments in the infrastructure and in that case, there is a difference between the DH system and heat pumps as the infrastructure for the heat pumps is relatively cheaper. This difference is accounted for in the fixed part of the heat tariff and is translated into additional income in the case of heat pumps. We explain this in the â€˜Profit rate for selling heatâ€™ subsection.
The equation which expresses the net yearly profit as objective function for maximization is defined as follows:
Electricity generated by the CHP is partly consumed by the electric equipment of the DH network (pumps, heat pumps, chiller). One of the objectives of the smart grid control is to balance power generation and consumption as much as possible. The objective is to match the daily CHP operating times with the operating times of electric and thermal consumers. The smaller the CHP, the longer it will run daily and the higher the chance of matching its generation power with consumption. Therefore, we have implemented a simple algorithm which calculates the fraction of CHP operating time matching operational times of electrical demand by the heat pumps, chiller, and network pumps. This is combined with an algorithm that determines the electricity costs and profit for buying energy from the grid and feeding into the grid.
Investment capital rate I and investments (INV) (unit: â‚¬) for the CHP, boiler, heat pumps, and storage are related to size or capacity and are described with power functions. Suitable power function coefficients are obtained from cost engineering handbooks, field experts, and company quotations.
Capital costs are expressed as net present value and are determined by initial investments, the interest rate r which we assume the same for all equipment, and equipment life time (LT) are as follows:
Operational costs for each equipment are assumed as a fixed percentage of the initial investment.
Thermal demand specification
Approach and assumptions
Thermal demand includes demand for domestic hot water and space heating. As the district does not yet exist, we develop a straightforward design calculation method in subsequent sections to generate an average demand profile from average daily ambient temperature data for the district.
Domestic hot water demand
For a district, the demand is determined by the amount of households and the number of persons per household. The following equation based on [12] expresses the variation of daily heat demand for domestic hot water due to holidays throughout the year:
In which n the number of households, p the fraction of households not on holiday during the holiday season, and day_shift the number of days that the minimum of the cosine function (peak day of the holiday season with the least domestic hot water demand) is shifted forward from the first of July. In reality, demand profiles are not smooth functions, hence a random daily variation on the profile is generated between 80% and 120%. This is shown in the case application in the subsequent subsection.
Equations for space heating demand
The energy demand for space heating can be predicted by taking a number of assumptions into account:

Heat loss to the environment. In general, the method we apply to construct a daily demand profile is to relate space heating demand to the temperature difference between a constant base temperature (T _{base}) and the average daily ambient temperature (\(\bar {T}_{a,i}\)). The base temperature is defined as the average daily ambient temperature for which no space heating is required. Due to solar and internal thermal gains within the interior of a house, the base temperature is lower than the heating setpoint temperature and can be found from experience (i.e., observed number of heating days) or dynamic simulation. The following constraints signifying the space heating demand limit apply:
$$ \ \Delta T_{\text{heat},i}= \left\{ \begin{array}{ll} T_{\text{base}}\bar{T}_{a,i} & ~~\text{if}~~ (T_{\text{base}}\bar{T}_{a,i}>0\\ 0 & ~~\text{if}~~ (T_{\text{base}}\bar{T}_{a,i}\leq 0 \end{array} \right. $$((11)) 
For the daily sum of the space heating demand, Q _{sh,i } (unit MJ/day), a simple prediction profile for a single household is constructed from the following relation:
$$ Q_{\text{sh},i}=U \ \Delta T_{\text{heat},i} $$((12))In which parameter U (unit MJ/K) is a constant total heat transfer coefficient, signifying all heat loss to the ambient due to envelope conduction, air ventilation and infiltration, and gains by solar radiation, residents, and appliances. U is determined from a known yearly sum of space heating demand as follows:
$$ U=\frac{\sum_{i} Q_{\text{sh},i}}{\sum_{i} \Delta T_{\text{heat},i}} $$((13)) 
The yearly sum Q _{sh} (unit MJ/y) indicated in Equation 13 is determined either by simulation, from literature or energy performance coefficient (EPC) calculation tools which are developed specifically for Dutch houses. We used the EPC verification tool EPG & Kosten [13] which is distributed by the Dutch government free of charge. Within this tool, various house types (AgentschapNL reference houses, [14]) with reference dimensions based on comparative studies are defined. Building and installation details can be selected and the program calculates the relevant yearly energy demands, based on EPC reference calculations, involving the reference Dutch climate year. The EPC value is the most important parameter for determining the yearly space heating demand. According to current building regulations, an EPC value of 0.4 will be effective from 1 January 2015 for all newly built houses within the Netherlands. Houses within a district heating system which is based on renewable energy receive an EPC reduction, according to [15]. An equivalent generation efficiency of the DH system is calculated with a calculation tool [16], also distributed by the Dutch government.
Case related thermal demand specification
Based on the methodology developed in the previous subsections, the following specific parameters are defined for the Meppel case:

Domestic hot water demand. Blokker and Poorteman [17] investigate hot water patterns of an urban district built in 2002 which consists of 3,000 houses. The population composition (age, type of households) compares with the target population of the Meppel case. On average, 2.6 persons per household are calculated. According to [17], [18], and [19] the daily hot water demand per person is 60 L/day of 45Â°C. With an average cold water supply temperature of 12Â°C, this amounts to 8.3 MJ/day Â· person and for an average household of 2.6 persons, this yields Q _{dhw,hh}=21.6 MJ/day and 7.88 GJ/year.

Space heating and cooling demand. Ambient temperature data for Meppel is taken from weather station in Hoogeveen which is close to Meppel. We determine the average daily temperature for the preceding 5 years. By comparing with simulations, we calculate T _{base}=14Â°C. The number and type of houses is estimated from the current progress of the Meppel house building project. The target EPC value is estimated at 0.4. A mixed DH system based on biogas CHP and natural gas boilers yields an equivalent generation efficiency of 1.5 and results in an EPC reduction of 0.15. Hence, we estimate that the real heat demand of houses is designed for an EPC value of 0.55. With this, the heating and cooling demand values presented in Table 3 are calculated. The resulting district heating demand profile per average household is shown in Figure 2. The profile is constructed by applying Equations 12 and 10 which are added together. The cooling demand profile is shown in Figure 3.
Relation between CHP size and operational time
When a CHP is operated, it is supported by boilers. A suitable size relation between CHP and boiler size is determined from the coldest day generation balance:
On the coldest day, both heat generators will run continuously for 24 h to generate the required thermal energy. Hence, the following relation between the heat demand and required peak thermal power applies for the coldest day:
If the daily heat demand data is sorted to a data array Q _{demand} from days with the largest to the smallest heat demand, Equation 15 is equal to the first value in this array. In Figure 4, the sorted array is drawn. We have also drawn the shaded heat demand coverage of a certain CHP which has just enough capacity to generate the heat demand indicated with a dot. The remaining area above the dotted line is covered by natural gas boilers.
For the sorted array which we convert from GJ/day to kWh/day, let \(\dot {Q}_{i}\) be a variable CHP generation capacity such that on day i, 24h operating time with this capacity results in exactly the heat demand of that day, i.e., Q _{demand,i }. For the shaded area of Figure 4, it can be demonstrated that the following equations determine the yearly operating time t _{op} for CHP generation capacity \(\dot {Q}_{i}\):
In this relation, \(\dot {Q}_{i}\) can be any CHP thermal capacity between 0 and \(\dot {Q}_{\text {max}}\). Equation 16 then yields the operating time of that CHP. Figure 5 shows the resulting relation between CHP size and its full power total operation time for the whole year, t _{op,chp}, calculated with Equation 16.
Case related energy cost parameters
Profit and cost rate of electricity
In the Meppel case, the utility company Meppel Energie buys electric energy from the grid and in case of a CHP it will also feed in surplus electricity into the grid. Contractual buying and selling prices of energy are usually based on the expected amount of energy drawn from and fed into the grid within a certain time period. Following Dutch CHP practice, we assume the buying and selling prices are equal for equal amounts of energy transported to and from the grid over the period of 1 year.
In general, the price of energy drops when the amount of energy being bought increases. For very large amounts (i.e., above 500.000 kWh/year) electricity prices are close to the spot market or day ahead market price of electricity, which is currently around 5.5 â‚¬ct/kWh. In case of CHP generation with various electricity consumers, if the amount of sold (surplus) energy exceeds the amount of bought energy, network utilities are often reluctant to pay any compensation, or spot market prices at most. Specific energy price equations are developed for the Dutch situation and implemented into the model.
Profit rate for selling heat
As we explain in the â€˜Optimization model equationsâ€™ subsection, infrastructural investments are payed off by the residents with a onetime connection fee. This fee is usually part of the building price of a house. Besides the connection fee, residents are charged for the consumed heat according to the following equation, i.e., the Dutch heat law, effective from January 2014 which defines the maximum possible tariff: P _{ Q,h }=209,92+19.86Â·Q _{demand,h } (excluding VAT).
Financial reservations have to be made for reinvestments in the required infrastructure. For this, there is a difference between houses connected to the DH system and houses with heat pumps:

In the case of houses connected to the DH system, we estimate the yearly costs of reservations for infrastructural reinvestments at approximately: â‚¬146.00 /(year Â· household). This leaves as income for the profit model, refer to Equation 8: P _{ Q,dh,h }=63.00+19.86Â·Q _{demand,dh,h }

In case of houses with a heat pump, we estimate lower reservations for infrastructural reinvestments due to the use of noninsulated pipes and a generation building is not required for heat pumps, i.e., â‚¬65.00 /(year Â·household). This leaves as income: P _{ Q,hp,h }=145.00+19.86Â·Q _{demand,hp,h }.
Profit rate of cooling
Originally, the Meppel project was planned for a larger starting scale, including the use of the aquifer for cooling energy and source energy for the heat pumps. Cooling of the houses provides part of the required energy to balance temperature of the aquifers and is therefore offered for free to the residents. Meppel energy wants to keep its promise of free cooling to the residents and therefore the profit rate for cooling is 0, P _{ Q,c }=0.
Cost rate of biogas for CHP and natural gas for boilers
Cost relations for natural gas from the national grid, specific relations for biogas, and applicable energy taxes are developed and implemented as algorithms into the optimization model.
Results and discussion
Generator and storage capacities
For the two building phases (phase 1: 40 DH houses; phase 2: 160 DH houses, 16 HP houses) maximum required generation capacity results are shown in Table 4. The peak boiler capacity mentioned is only required for infrequent extremely cold days and also serves as a heat generation backup in case of CHP or boiler failure or during maintenance periods.
Profit generation, optimal CHP size, and operating times
The relation between profit and optimal size of the CHP in relation to boiler size is investigated by maximizing the net profit equation (refer to Equation 8) with \(\dot {Q}_{\text {chp}}\) as variable. For the two building phases with (n,m) houses, the optimal sizes in Table 5 are calculated.
Table 5 shows a negative profit for both building phases. Care should be taken in the interpretation of the profit values as they are influenced by assumptions explained in the â€˜Profit rate for selling heatâ€™ subsection. We may have exaggerated the reinvestment costs, lacking detailed insight in the Meppel Energie business case. On the other hand, practical experience indicates problematic profitability of smallscale DH systems, which support our results. In the Meppel case, there is more pressure on profitability because the houses have a lower heating demand (and thus lower revenue per household) than usual, and there are generation costs involved for cooling while on the other hand, cooling is offered free of charge. If we calculate a fee for cooling based on break even operation, this yields â‚¬64/year and â‚¬70/year per household for phases 1 and 2, respectively. To place this in perspective of what households would normally pay if they generate the required cooling with homebased air conditioning equipment, the average cooling demand of 6 GJ/year requires 555 kWh per household electric energy which would costs each household â‚¬127/year only for electricity. Hence, the calculated fees are reasonable propositions to the residents.
What is the influence of houses with heat pumps on the profit? By maximizing the net profit Equation 8 with \(\dot {Q}_{\text {chp}}\) and m as variables, we find that for various numbers of n houses, m=0. So in all cases, houses with heat pumps have a negative influence on profit. There are two reasons for this. First, in the starting energy concept, the source heat for the heat pumps is generated by the CHP and boilers, with fuel costs as a consequence. Second, individual heat pumps for each house are relatively expensive in comparison with DH delivery sets as HPâ€™s require more investments and maintenance and also have a shorter life time. Clearly, the compensation offered by a higher tariff (â€˜Profit rate for selling heatâ€™ subsection) is not sufficient for an equally profitable operation of heat pumps compared to houses connected to the DH system.
If the number of houses connected to the DH system expands, when will operation become profitable? This depends mostly on (a) the applied interest rate r and (b) the household heating demand. If r=4% then operation breaks even at n=490 houses. The corresponding optimal generation sizes and operating times are as follows: \(\dot {E}_{\text {chp}}=426\) kW, t _{op,chp}=5,293 h/year, \(\dot {Q}_{b}=448\) kW, and t _{op,b }=1,340 h/year. For r=2% the break even point is n=277 houses.
CHP migration steps
When the number of houses in the district increases, ideally the CHP and boiler sizes also have to increase according to the optimal profitability found in the previous subsection. It is relatively easy to install more boiler capacity, but replacing or adding a single CHP requires much higher investments and more accompanying equipment. CHP expansion is possible in two ways, either by leasing CHP units for a period of time or by investing in a starting CHP unit and investing into more CHP units when this is required. The latter option has the advantage that CHP capacity is gradually expanded and eventually the generation capacity exists of a number of smaller CHP units, which enables smoother power regulation and operation of each CHP under full load conditions, which is beneficial for fuel efficiency and technical reliability. However, this option requires more investments, although the Dutch tax scheme offers substantial tax reductions for CHP investments, so this is selected as the most interesting approach. Therefore, we investigate the relation between CHP size and profit in four steps from building phase 1 to phase 2, to estimate the proper CHP sizes. The result is shown in Figure 6.
For each line in the Figure, the value for (n,m) is indicated. Around the optimum of each line, the relation between profit and CHP size is rather flat. This gives flexibility to choose the CHP size according to logical migration steps for the district. An attractive migration route is to start with a 30kWe CHP for a small district with 40 houses and to install a second, larger CHP, e.g., 140 kWe when 160 houses are connected.
Comparison on sustainability
Sustainability of the Meppel energy concept shown in Figure 1 is investigated by comparing fossil fuel requirements for building phase 2 (n=160, m=16 houses), for which we define four cases:

1.
Case 1: starting energy concept with optimal sized CHP on biogas supported by a boiler on natural gas. The Sankey or energy flow diagram for this case is shown in Figure 7. Total natural gas input: 1,648 GJ/year. Total biogas input: 6,177 GJ/year. If the CHP runs on natural gas, the total natural gas input is: 7,825 GJ/year.

2.
Case 2: starting energy concept without a CHP. In this case, the boiler delivers 4,860 GJ/year with an equal natural gas input. The required electricity for heat pumps, cooling and network pumps, and energy delivered to the region total 2,222 GJ/year. This requires 4,938 GJ/year of natural gas input at a gridbased electrical power plant, taking Dutch national fuel efficiency of 45% into account. Hence, total natural gas input: 9,798 GJ/year.

3.
Case 3: Dutch conventional house heating by individual homefitted natural gas boilers. In this case for heating, 4,318 GJ/year natural gas is required. The required electricity for cooling and comparable electricity supplied to the region total 1,816 GJ/year. If a gridbased electrical power plant is used, 4,036 GJ/year natural gas input is required. In that case, the natural gas input totals: 8,354 GJ/year.

4.
Case 4: as an alternative for case 3, the 1,816 GJ/year electricity is generated by home installed solar PV. For 176 houses, this is 2,866 kWh/year per household. Under Dutch circumstances, this is feasible. Total natural gas input: 4,318 GJ/year only for heating is required.
The resulting natural gas input is shown in Figure 8. Case 1 shows the best sustainability results; however, if natural gas is used for the CHP, the concept is only 6% more sustainable than case 3, i.e., conventional house heating with natural gas boilers. Case 2 shows that the DH system without a CHP actually has the worst performance on sustainability.
Case 4 with home PV installations is interesting as improved sustainability is combined with less infrastructural investments than case 1 and is hence easier to implement in an unsure house building market. Case 4 also leads to another objective, i.e., to invest into more energy efficient houses, in order to further reduce natural gas consumption. We recalculated case 4 with passive house (PH) space heating requirements, and this results in the fifth column in Figure 8. This case can be taken further, replacing natural gas boilers with heat pumps. The district then contains passive houses, each with a solar PV roof and a heat pump connected to a ground source or aquifer. Natural gas input is then eliminated totally. The drawback of such a concept however is the imbalance between periods of major PV electricity generation and periods with the highest heating demand. Further comparisons of these concepts including financial aspects and loads on the electricity network are left for future work.
The use of biogas instead of natural gas reduces carbon dioxide (CO _{2}) emissions, while biogas is formed with CO _{2} from the earthâ€™s atmosphere and the CO _{2} is released when the biogas is burned. For natural gas, a conversion factor of 52 kg CO _{2}/GJ is derived. Case 1 with biogas results in 86 tons CO _{2}/year and case 2 in 509 tons CO _{2}/year. Hence, case 1 results in a CO _{2} reduction of 423 tons/year.
The last question of interest for the Meppel case is how many houses and what size of the CHP are appropriate if all the available biogas of the municipality (i.e., 9,200 GJ/year) is applied? As optimum, we find n=232 houses and \(\dot {E}_{\text {chp}}=159\) kW. The natural gas input required for boilers is then 1,775 GJ/year, resulting in net CO _{2} emissions of 92 tons/year. The net profit (or loss in this case) is calculated at âˆ’â‚¬5,711/year. When a cooling fee of â‚¬50/year is asked from each household, the net profit is positive, i.e., â‚¬5,899/year.
Conclusions
In this paper, optimal heat and power generation sizes are determined for the new urban district Nieuwveenslanden situated in Meppel for two building phases starting at 40 houses and expanding towards 176 houses. Besides optimal CHP and boiler sizes, we determine corresponding profits and sustainability of the network on fossil fuel requirements. Profits appear to be negative, however more than 50% improvement is possible if a CHP is used for heat and power generation compared to operation with only boilers. With an interest rate of 4% on invested capital, profits are positive when more than 490 houses are connected, with 2% this is reduced to 277 houses. If households are charged a modest fee for cooling, the energy concept with a CHP is profitable from the start. Based on maximizing profit, a migration scenario for the CHP is developed starting with a single 30kWe CHP for 40 houses connected to the DH system and adding a 140kWe CHP when 160 houses are connected.
Sustainability of the energy concept is compared with the Dutch reference, i.e., houses heated by individual natural gas boilers. If a CHP on biogas is applied, the Meppel energy concept reduces 351 tons CO _{2}/year or 80% reduction. However, if the CHP runs on natural gas, the CO _{2} reduction is only 6%. It has to be considered that a scenario based on the Dutch reference combined with home solar PV is another possible route towards improved sustainability but involves less infrastructural investments. However, the Meppel DH system concept with a CHP on biogas still performs better on sustainability than this alternative route.
For the Meppel DH system, houses with a heat pump have a negative influence on the profitability. This is because individual household heat pumps are more expensive to operate than the DH system. It is interesting to compare this also for the final energy concept in which source heat and cooling are generated by an underground aquifer. In that case, we expect to reach better profitability. Investigating the final energy concept is left for future work. Also a more comprehensive comparison between the Meppel energy concept and more individual energy generation solutions requiring less infrastructure is left for future work.
References
AgentschapNL (2014) Biogaswarmtekrachtkoppeling (WKK) voor Een Wijk. http://www.rvo.nl/onderwerpen/duurzaamondernemen/gebouwen/energieprestatienieuwbouwepn/ontwerpen. referentiewoningennieuwbouw. visited January 2014.
Schepers BL, Valkengoed MPJv (2009) Warmtenetten in nederland  overzicht van grootschalige en kleinschalige warmtenetten in nederland. Technical Report 09.3031.45, CEDelft http://www.ce.nl/publicatie/warmtenetten_in_nederland/976. visited May 2015.
Vallios I, Tsoutsos T, Papadakis G (2009) Design of biomass district heating systems. Biomass and bioenergy 33(4): 659â€“678.
Ren H, Gao W, Ruan Y (2008) Optimal sizing for residential CHP system. Appl Therm Eng 28(5): 514â€“523.
Cho H, Luck R, Eksioglu SD, Chamra LM (2009) Costoptimized realtime operation of CHP systems. Energy and Buildings 41(4): 445â€“451.
Chinese D, Meneghetti A (2005) Optimisation models for decision support in the development of biomassbased industrial districtheating networks in Italy. Appl energy 82(3): 228â€“254.
Curti V, Favrat D, von Spakovsky MR (2000) An environomic approach for the modeling and optimization of a district heating network based on centralized and decentralized heat pumps, cogeneration and/or gas furnace. Part II: application. Int J Ther Sci 39(7): 731â€“741.
Hlebnikov A, Dementjeva N, Siirde A (2009) Optimization of Narva district heating network and analysis of competitiveness of oil shale CHP building in Narva. Oil Shale 26: 269â€“282.
Lund H, MÃ¶ller B, Mathiesen BV, Dyrelund A (2010) The role of district heating in future renewable energy systems. Energy 35(3): 1381â€“1390.
Pirouti M (2013) Modelling and analysis of a district heating network.
Ã˜stergaard PA, Lund H (2011) A renewable energy system in Frederikshavn using lowtemperature geothermal energy for district heating. Appl Energy 88(2): 479â€“487.
Jordan U, Vajen K (2001) Realistic domestic hotwater profiles in different time scales. Report for IEASHC Task 26. http://sel.me.wisc.edu/trnsys/trnlib/ieashctask26/ieashctask26loadprofilesdescriptionjordan.pdf, visited May 2015.
RVO (2014) EPG Rekentool. https://www.rvo.nl/onderwerpen/duurzaamondernemen/gebouwen/energieprestatienieuwbouwepn/epcberekenen/hulpmiddelenepcberekening?wss%l=1, visited January 2014.
AgentschapNL (2013) Referentiewoningen Nieuwbouw 2013. Agentschap NL. http://www.rvo.nl/onderwerpen/duurzaamondernemen/gebouwen/energieprestatienieuwbouwepn/ontwerpen/referentiewoningennieuwbouw visited January 2014.
AgentschapNL (2014). http://www.rvo.nl/sites/default/files/Gebiedsmaatregelen\%20voor\%20het\%20eerst\%20gewaardeerd\%20in\%20de. %20EPC_bepaling_0.pdf, visited January 2014.
RVO (2014) EMG Rekentool. https://www.nen.nl/NENShop/Norm/Updaterekentoolenergieprestatiegebiedbeschikbaar.htm, visited January 2014.
Blokker EJM, Poortema K (2007) Effecten levering warm tapwater door derden. Technical report, sKIWA.
Geudens PJJG (2012) Drinkwaterstatistieken 2012. VEWIN 41(4): 445â€“451.
Blokker EJM (2006) Modelleren van afnamepatronen; beschrijving en validatie van het simulatiemodel simdeum. Technical report, KIWA. http://www.kwrwater.nl/WorkArea/DownloadAsset.aspx?id=8787. visited May 2015.
Acknowledgements
The authors would like to thank the Dutch national program TKISwitch2SmartGrids for supporting the project MeppelEnergie and the STW organization for supporting the projects iCare and Dream that resulted in this work.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authorsâ€™ contributions
RVL carried out the research for this paper, developed the models, performed all simulations, and drafted the manuscript. JF participated in the model development and critically revised the paper. JDW defined alternative energy systems to compare on sustainability and supplied background case related data to compare with. GS critically revised the paper and gave final approval for the paper to be published. All authors read and approved the final manuscript.
Authorsâ€™ information
JF is affiliated to the STW project iCare and Dream.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
van Leeuwen, R.P., Fink, J., de Wit, J.B. et al. Upscaling a district heating system based on biogas cogeneration and heat pumps. Energ Sustain Soc 5, 16 (2015). https://doi.org/10.1186/s137050150044x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s137050150044x
Keywords
 Cogeneration
 Biogas
 District heating
 Heat pumps
 Renewable energy
 Lowenergy building
 System migration