# Optimal Hierarchical Decision-Making for Heat Source Selection of District Heating Systems

Authors: Fang Fang; Nan Wang
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

With the rapid development of China’s urbanization, the proportion between the heating consumption and the energy consumption of the whole society keeps rising in recent years. For a district heating system, the selection of the heat source makes significant impact on the energy efficiency and the pollutant emissions. By integrating the methods of the Analytic Hierarchy Process (AHP) and the Preference Ranking Organization Method for Enrichment Evaluation (PROMETHEE), a multiple-attribute decision-making scheme for the heat source selection of district heating systems is proposed in this paper. As a core part of this scheme, a comprehensive benefit index with hierarchical parallel structure is constructed. The economic benefit, environment benefit, and technical benefit can be reflected with a certain percentage in the comprehensive benefit index. To test the efficiency of the proposed scheme, a case study for a large-scale district heating system in Beijing is carried out, where five kinds of heat sources (water source heat pump, ground source heat pump, gas-fired boiler, coal-fired boiler, and oil-fired boiler) are taken into account. The analysis and instructions for the final sorting result are also demonstrated.


---

## Body

## 1. Introduction

The district heating is a system for the distributing heat that is generated in a centralized location for residential and commercial heating requirements. Compared with the localized heating system, the district heating system can reach higher efficiency and provide better pollution control. With the rapid development of China’s urbanization, the district heating system has become an indispensable infrastructure. Till the end of 2013, the district heating area in China has reached 63.6 billion square meters and this number is six times as many as in 2000 [1]. But the higher energy consumption, the lower equipment efficiency, and the unreasonable load match between the capacity of heating equipment and the users’ load at the end of the heating network still restrict the healthy development of the district heating systems.As we know, there are three essential elements of the district heating systems: the heat sources, the heating pipe networks, and the heat users. Among them, the heat source is not only the heat provider for the whole heating system, but also the direct consumer of the primary energy. And the overall safety, efficiency, and emissions of the heating system have close relationship with heat sources. Furthermore, the operation mechanism, the energy efficiency, and the primary energy consumption are of significant difference for various types of heat sources. All of the above reasons make the selection of heat sources become the top priority in the planning of district heating systems.In the face of such a typical complex networked system, the advanced real-time control [2] and the online optimization [3] will be instructive for the daily operation of the district heating system, while performance evaluation and the energy-loss analysis can give important guidance for the planning and construction of heat sources. Hu et al. [4] gave comparative analysis of the characteristics and the energy consumption for five kinds of common heat sources. Chen et al. [5] put forward suggestions and measures for improving the efficiency of heat sources by analyzing the consumptions of fuel, electrical power, water, and heat. Zheng [6] discussed several key problems about energy saving from heat sources, heating networks, heating stations, and secondary-level heating networks.Referring to the advanced experience and considering the needs of social and economic development, multisource and complementarity of heat sources have become the developing trend of the configuration for district heating systems [7]. This situation makes the selection of heat sources a multiobjective and multiconstraint decision-making problem, where the technical performance, the social benefit, and environmental benefit should be considered comprehensively [8]. Wang et al. [9] studied the environmental assessment of heat sources in urban area using fuzzy multiobjective optimization method. Zheng et al. [10] introduced AHP in the evaluation of district heating systems and made a comprehensive evaluation based on the total economic cost. Ghafghazi et al. [11] proposed a multicriteria approach, which can take into account the interests of investors, operators, users, and other stakeholders, to evaluate the district heating system options.All of the above results suggest the following.(i)
The district heating system has complicated network characteristics that determine the existence of strong technical correlation between the heat source and other elements.(ii)
The heat source, located in a special geographical and social environment, also has significant social and economic attribute.Therefore, before making the decision of heat source selection, the multilevel and multiattribute index system should be established first. And the advanced evaluation and decision-making methods may help to obtain the final result, which can satisfy the global demand under most conditions.Following the previous researches and referring to advanced decision-making methods applied in complicated systems, in this paper, the Analytic Hierarchy Process (AHP) and the Preference Ranking Organization Method for Enrichment Evaluation (PROMETHEE) are integrated to produce an optimal hierarchical decision-making scheme for heat source selection of district heating systems. The AHP is used for the structure analysis of the candidate heat sources and for the determination of the relative weights of evaluation indices; the PROMETHEE is introduced to sort the candidate heat sources and conclude the optimal results. The combination of the AHP and the PROMETHEE can coordinate the qualitative indices and quantitative indices effectively and can make complicated decision-making problems coherent and structured.The rest of this paper is organized as follows: a scheme of multiple-attribute decision-making is discussed in Section2; in Section 3, a comprehensive evaluation index including economic benefit, environment benefit, and technical benefit is proposed for the selection of heat sources; a large-scale district heating system in Beijing is taken as an example for a case study based on the proposed hierarchical decision-making scheme in Section 4; finally, the results and discussion of this work are addressed in Section 5.


## 2. Multiple-Attribute Evaluation Scheme

### 2.1. Determine the Preference of Alternative by the PROMETHEE

By introducing a preference function, the PROMETHEE method can establish a preference relationship for each alternative. According to the calculation of incoming flow and outgoing flow, an incomplete rank for each alternative can be obtained [12].Assuming there arek plans in a certain alternative set X=x1,x2,…,xk and h indices are included in evaluation index set F=f1,f2,…,fh.For arbitrary two alternativesxa∈X and xb∈X, define the preference structure under the index fc∈F:
(1)aPbfcifffcxa>fcxb,aIbfcifffcxa=fcxb,
where aPbfc represents that the preference degree of xa is larger than that of xb for decision makers under the index fc, aIbfc represents that the preference degrees of xa and xb are equal, and fcxa and fcxb are attribute values of xa and xb under fc.There are two main types of methods to get the attribute valuefcxi: quantitative method and qualitative methods. The former works with the practical engineering calculation or experiment results; the latter is mainly based on rating scale method to decide the importance of alternatives for decision-making.In order to describe the precedence relationships of the alternatives, the preference functionPcxa,xb for index fc is defined to transfer these relationships to numerical values located in [0,1]. The preference function usually has six common types [11]. According to the need of this paper, two types of preference functions are selected for decision-making as shown in Table 1.Table 1
Two types of preference functions.
Type
Function name
Form of preference function
I
Normal preference function
P
c
(
x
a
,
x
b
)
=
{
1,
if
f
c
(
x
a
)
-
f
c
(
x
b
)
>
0,
0,
if
f
c
(
x
a
)
-
f
c
(
x
b
)
≤
0,
II
Quasi preference function
P
c
(
x
a
,
x
b
)
=
{
1,
if
f
c
(
x
a
)
-
f
c
(
x
b
)
>
e,
0,
if
f
c
(
x
a
)
-
f
c
(
x
b
)
≤
eIn order to give a quantitative description for the preference degree between any two alternatives clearly, the preference coefficientπxa,xb is defined as follows:
(2)πxa,xb=∑c=1hωcpcxa,xb∑c=1hωc,
where ωc is the relative weight of fc. It can be calculated by AHP, which will be discussed in Section 2.2.Furthermore, the outgoing flowϕ+xa, incoming flow ϕ-xa, and net flow ϕxa are defined for any alternative xa as follows:
(3)ϕ+xa=∑i=1kπ(xa,xi),ϕ-xa=∑i=1kπ(xi,xa),ϕxa=ϕ+xa-ϕ-xa.Outgoing flowϕ+xa represents the superior degree of alternative xa to other alternatives. The bigger the value of outgoing flow is, the better the alternative xa is. Incoming flow ϕ-xa represents the worse degree of alternative xa than other alternatives. The bigger the value of incoming flow is, the worse the alternative xa is. Net flow ϕxa represents the total preference relationships and is used to determine the final ranking result of alternatives.


### 2.2. Determine the Relative Weight of Evaluation Index by the AHP

The AHP is usually applied in the analysis for decision-making of complex system with multiattribute, multifactors, multilevels, and multicandidates [13]. Here we use it to calculate the relative weights of evaluation indices.In order to simplify the calculation, the relative weights of evaluation indices are divided into two levels: the objective layer and the attribute layer. The objective layer includes the relative weight final values of evaluation indices; the attribute layer includes the relative importance of evaluation indices obtained by experts’ scoring.Here we also take the alternative setX=x1,x2,…,xk and evaluation index set F=f1,f2,…,fh as examples. The evaluation index set represents the objective layer. The judgment matrix Dh, established based on h evaluation indices as shown in Table 2, makes up the attribute layer.Table 2
The judgment matrixDh.
f
1
f
2
⋯
f
h
f
1
d
11
d
12
⋯
d
1
h
f
2
d
21
d
22
⋯
d
2
h
⋮
⋮
⋮
⋮
⋮
f
h
d
h
1
d
h
2
⋯
d
h
hIn Table2, dpq is the relative importance of index fp∈F comparing with index fq∈F, and there is dpq=1/dqp. When fp=fq, dpq=1. The rating scale method is often used to get the relative importance according to experts’ scoring. Table 3 shows a typical nine-class scale.Table 3
The definition table of nine-class scale.
Rating scale
Definition (comparing indexfp with index fq)
1
Indicesfp and fq are important equally
3
Indexfp is slightly more important than fq
5
Indexfp is more important than fq
7
Indexfp is more important than fq
9
Indexfp is absolutely more important than fq
2, 4, 6, 8
Middle important value between above scalesThe AHP determines the relative weight for each evaluation index by calculating the weight coefficient of the attribute layer to the objective layer. The geometric mean algorithm can be used for deciding the weight with the advantage of judgment matrix.(1)
Calculate thehth root of the product of h elements inpth row of the judgment matrix:
(4)ωp′=∏q=1hdpqh.(2)
For indexfp, the relative weight can be calculated by
(5)ωp=ωp′∑p=1hωp′.Through building judgment matrix, the relative weight vector ofh evaluation indices can be obtained as ω=ω1,ω2,…,ωh.


## 2.1. Determine the Preference of Alternative by the PROMETHEE

By introducing a preference function, the PROMETHEE method can establish a preference relationship for each alternative. According to the calculation of incoming flow and outgoing flow, an incomplete rank for each alternative can be obtained [12].Assuming there arek plans in a certain alternative set X=x1,x2,…,xk and h indices are included in evaluation index set F=f1,f2,…,fh.For arbitrary two alternativesxa∈X and xb∈X, define the preference structure under the index fc∈F:
(1)aPbfcifffcxa>fcxb,aIbfcifffcxa=fcxb,
where aPbfc represents that the preference degree of xa is larger than that of xb for decision makers under the index fc, aIbfc represents that the preference degrees of xa and xb are equal, and fcxa and fcxb are attribute values of xa and xb under fc.There are two main types of methods to get the attribute valuefcxi: quantitative method and qualitative methods. The former works with the practical engineering calculation or experiment results; the latter is mainly based on rating scale method to decide the importance of alternatives for decision-making.In order to describe the precedence relationships of the alternatives, the preference functionPcxa,xb for index fc is defined to transfer these relationships to numerical values located in [0,1]. The preference function usually has six common types [11]. According to the need of this paper, two types of preference functions are selected for decision-making as shown in Table 1.Table 1
Two types of preference functions.
Type
Function name
Form of preference function
I
Normal preference function
P
c
(
x
a
,
x
b
)
=
{
1,
if
f
c
(
x
a
)
-
f
c
(
x
b
)
>
0,
0,
if
f
c
(
x
a
)
-
f
c
(
x
b
)
≤
0,
II
Quasi preference function
P
c
(
x
a
,
x
b
)
=
{
1,
if
f
c
(
x
a
)
-
f
c
(
x
b
)
>
e,
0,
if
f
c
(
x
a
)
-
f
c
(
x
b
)
≤
eIn order to give a quantitative description for the preference degree between any two alternatives clearly, the preference coefficientπxa,xb is defined as follows:
(2)πxa,xb=∑c=1hωcpcxa,xb∑c=1hωc,
where ωc is the relative weight of fc. It can be calculated by AHP, which will be discussed in Section 2.2.Furthermore, the outgoing flowϕ+xa, incoming flow ϕ-xa, and net flow ϕxa are defined for any alternative xa as follows:
(3)ϕ+xa=∑i=1kπ(xa,xi),ϕ-xa=∑i=1kπ(xi,xa),ϕxa=ϕ+xa-ϕ-xa.Outgoing flowϕ+xa represents the superior degree of alternative xa to other alternatives. The bigger the value of outgoing flow is, the better the alternative xa is. Incoming flow ϕ-xa represents the worse degree of alternative xa than other alternatives. The bigger the value of incoming flow is, the worse the alternative xa is. Net flow ϕxa represents the total preference relationships and is used to determine the final ranking result of alternatives.


## 2.2. Determine the Relative Weight of Evaluation Index by the AHP

The AHP is usually applied in the analysis for decision-making of complex system with multiattribute, multifactors, multilevels, and multicandidates [13]. Here we use it to calculate the relative weights of evaluation indices.In order to simplify the calculation, the relative weights of evaluation indices are divided into two levels: the objective layer and the attribute layer. The objective layer includes the relative weight final values of evaluation indices; the attribute layer includes the relative importance of evaluation indices obtained by experts’ scoring.Here we also take the alternative setX=x1,x2,…,xk and evaluation index set F=f1,f2,…,fh as examples. The evaluation index set represents the objective layer. The judgment matrix Dh, established based on h evaluation indices as shown in Table 2, makes up the attribute layer.Table 2
The judgment matrixDh.
f
1
f
2
⋯
f
h
f
1
d
11
d
12
⋯
d
1
h
f
2
d
21
d
22
⋯
d
2
h
⋮
⋮
⋮
⋮
⋮
f
h
d
h
1
d
h
2
⋯
d
h
hIn Table2, dpq is the relative importance of index fp∈F comparing with index fq∈F, and there is dpq=1/dqp. When fp=fq, dpq=1. The rating scale method is often used to get the relative importance according to experts’ scoring. Table 3 shows a typical nine-class scale.Table 3
The definition table of nine-class scale.
Rating scale
Definition (comparing indexfp with index fq)
1
Indicesfp and fq are important equally
3
Indexfp is slightly more important than fq
5
Indexfp is more important than fq
7
Indexfp is more important than fq
9
Indexfp is absolutely more important than fq
2, 4, 6, 8
Middle important value between above scalesThe AHP determines the relative weight for each evaluation index by calculating the weight coefficient of the attribute layer to the objective layer. The geometric mean algorithm can be used for deciding the weight with the advantage of judgment matrix.(1)
Calculate thehth root of the product of h elements inpth row of the judgment matrix:
(4)ωp′=∏q=1hdpqh.(2)
For indexfp, the relative weight can be calculated by
(5)ωp=ωp′∑p=1hωp′.Through building judgment matrix, the relative weight vector ofh evaluation indices can be obtained as ω=ω1,ω2,…,ωh.


## 3. Comprehensive Evaluation Index for the Heat Source Selection

Combined with the actual needs of the current social and economic development in China, a comprehensive evaluation model for the heat source selection of the district heating system is presented in this paper. This model, giving consideration to social, environmental, and technical interests, consists of the following evaluation indices.(1)
Economic evaluation index: it is used to reflect the influence of heat demand from users on the economic benefit of the heating system.(2)
Environment evaluation index: it is used to reflect the impact on the natural and social environments during the construction and operation of the heating system.(3)
Technical evaluation index: it is used to reflect the influence of energy utilization efficiency on the energy saving operation of the heating system.Assuming that there arek alternatives X=x1,x2,…,xk for the selection of heat sources; the life period of the heating equipment is m years, and the number of installed equipment for each alternative is n. Accordingly, some variables and symbols are defined in Nomenclature attached at the end of this paper.


### 3.1. Economic Benefit Index

Economic benefit index includes the initial construction cost and operation cost. In order to give a quantitative description of economic benefit for each alternative, the dynamic annual annuity method [14] is used here for converting the initial investment cost and the operation cost to the equivalent annual annuity rcl(xj) oflth equipment injth alternative in its life cycle. Consider
(6)rclxj=Cinil(xj)-Lml(xj)11+Im+mCopelxj11+Im1+Im1+Im-1,Copel(xj)=∑l=1ncfuell(xj)+cdepl(xj)+cwell(xj).The value of economic benefit forjth alternative is
(7)RCxj=∑l=1nrclxj=fCinilxj,Lmlxj,Copelxj,I.HereRc(xj), reflecting the equivalent annual operation cost of an alternative in the entire operation period, is the function of initial investment, operation, and management cost, the equipment salvage, and the yield rate of the alternative.The economic benefit index ofjth alternative is the normalizing value of economic benefit. Consider
(8)PCxj=min⁡RmxjRmxj.


### 3.2. Environment Benefit Index

Due to the burning of fossil fuels, the boiler-type heat sources will produce gaseous, liquid, and solid pollutants inevitably. Moreover, although they do not consume primary energy, the electric power consumption of the heat-pump-type heat sources can also cause the pollutant emission during the energy conversion process. In this paper, the annual emissions of SO2, NOx, CO, CO2, and other pollutants are considered in the environment benefit evaluation, and the PROMETHEE is selected to calculate the quantitative result of the environment benefit index.Referring to Section2.1, the environment benefit equations ofjth alternative can be deduced from (3). Consider
(9)ϕE+xj=∑i=1kπxj,xi,(10)ϕE-xj=∑i=1kπxi,xj,(11)ϕExj=ϕE+xj-ϕE-xj,
where π(xj,xi) is relevant to the preference function and the weight coefficient. According to the characteristics of real data, the type II preference function is selected to describe NOx and SO2; the type I preference function is used for other pollutants. The weight vector set of five kinds of gas pollutants is denoted as ωE={ω1,ω2,ω3,ω4,ω5}. Referring to Table 2, a five-order judgment matrix D5 can be built. The value of every element in D5 is rated by decision makers and experts based on Table 3.According to the judgment matrixD5, the relative weight for each pollutant attribute can be obtained referring to (4) and (5). Furthermore, the economic benefit index ofjth alternative can be normalized as
(12)PExj=ϕExjmax⁡ϕExj.ϕ
E
(
x
j
)
<
0 means that the environment benefit ofjth alternative has no advantage compared with other alternatives. This alternative should be excluded directly from the alternative set.


### 3.3. Technical Benefit Index

Technical benefit index contains five aspects: the technical sophistication, the heating quality, safety and suitability, local resource condition, and energy utilization efficiency. Same as the calculation of environment benefit index, the PROMETHEE is adopted. Consider(13)ϕT+xj=∑i=1kπxj,xi,ϕT-xj=∑i=1kπxi,xj,ϕTxj=ϕT+xj-ϕT-xj.Since the real values of the five technical attributes are difficult to be measured in practice, the rating scale method, which draws support from expert scoring or questionnaire, is the common choice.The normalized technical benefit index ofjth alternative can be got by
(14)PTxj=ϕTxjmax⁡⁡ϕTxj.


### 3.4. Comprehensive Benefit Index

Based on the above analysis, a comprehensive benefit index for the heat source selection can be established:(15)Pxj=αcPCxj+αePExj+αtPTxj,
where αc+αe+αt=1 to ensure P(xj)∈0,1. The objective of the heat source selection is to find the alternative with maximal P(xj).The comprehensive benefit index evaluates the performance of the heat sources from three aspects: economy, environment, and technology. A typical hierarchical parallel structure can be seen between the three aspects and the comprehensive benefit index as Figure1.Figure 1
The hierarchical parallel structure of the comprehensive benefit index.


## 3.1. Economic Benefit Index

Economic benefit index includes the initial construction cost and operation cost. In order to give a quantitative description of economic benefit for each alternative, the dynamic annual annuity method [14] is used here for converting the initial investment cost and the operation cost to the equivalent annual annuity rcl(xj) oflth equipment injth alternative in its life cycle. Consider
(6)rclxj=Cinil(xj)-Lml(xj)11+Im+mCopelxj11+Im1+Im1+Im-1,Copel(xj)=∑l=1ncfuell(xj)+cdepl(xj)+cwell(xj).The value of economic benefit forjth alternative is
(7)RCxj=∑l=1nrclxj=fCinilxj,Lmlxj,Copelxj,I.HereRc(xj), reflecting the equivalent annual operation cost of an alternative in the entire operation period, is the function of initial investment, operation, and management cost, the equipment salvage, and the yield rate of the alternative.The economic benefit index ofjth alternative is the normalizing value of economic benefit. Consider
(8)PCxj=min⁡RmxjRmxj.


## 3.2. Environment Benefit Index

Due to the burning of fossil fuels, the boiler-type heat sources will produce gaseous, liquid, and solid pollutants inevitably. Moreover, although they do not consume primary energy, the electric power consumption of the heat-pump-type heat sources can also cause the pollutant emission during the energy conversion process. In this paper, the annual emissions of SO2, NOx, CO, CO2, and other pollutants are considered in the environment benefit evaluation, and the PROMETHEE is selected to calculate the quantitative result of the environment benefit index.Referring to Section2.1, the environment benefit equations ofjth alternative can be deduced from (3). Consider
(9)ϕE+xj=∑i=1kπxj,xi,(10)ϕE-xj=∑i=1kπxi,xj,(11)ϕExj=ϕE+xj-ϕE-xj,
where π(xj,xi) is relevant to the preference function and the weight coefficient. According to the characteristics of real data, the type II preference function is selected to describe NOx and SO2; the type I preference function is used for other pollutants. The weight vector set of five kinds of gas pollutants is denoted as ωE={ω1,ω2,ω3,ω4,ω5}. Referring to Table 2, a five-order judgment matrix D5 can be built. The value of every element in D5 is rated by decision makers and experts based on Table 3.According to the judgment matrixD5, the relative weight for each pollutant attribute can be obtained referring to (4) and (5). Furthermore, the economic benefit index ofjth alternative can be normalized as
(12)PExj=ϕExjmax⁡ϕExj.ϕ
E
(
x
j
)
<
0 means that the environment benefit ofjth alternative has no advantage compared with other alternatives. This alternative should be excluded directly from the alternative set.


## 3.3. Technical Benefit Index

Technical benefit index contains five aspects: the technical sophistication, the heating quality, safety and suitability, local resource condition, and energy utilization efficiency. Same as the calculation of environment benefit index, the PROMETHEE is adopted. Consider(13)ϕT+xj=∑i=1kπxj,xi,ϕT-xj=∑i=1kπxi,xj,ϕTxj=ϕT+xj-ϕT-xj.Since the real values of the five technical attributes are difficult to be measured in practice, the rating scale method, which draws support from expert scoring or questionnaire, is the common choice.The normalized technical benefit index ofjth alternative can be got by
(14)PTxj=ϕTxjmax⁡⁡ϕTxj.


## 3.4. Comprehensive Benefit Index

Based on the above analysis, a comprehensive benefit index for the heat source selection can be established:(15)Pxj=αcPCxj+αePExj+αtPTxj,
where αc+αe+αt=1 to ensure P(xj)∈0,1. The objective of the heat source selection is to find the alternative with maximal P(xj).The comprehensive benefit index evaluates the performance of the heat sources from three aspects: economy, environment, and technology. A typical hierarchical parallel structure can be seen between the three aspects and the comprehensive benefit index as Figure1.Figure 1
The hierarchical parallel structure of the comprehensive benefit index.


## 4. Case Study

A large-scale district heating system for a community with 330000 m2 heat supply area in Beijing is taken as an example for analysis. This heating system mainly undertakes the supply of daily living hot water and heating supply in winter. The heating supply period of Beijing in winter is 129 days, and the lowest indoor average temperature is 24°C. Referring to the algorithm given in [11], the peak heat load of this community can be estimated as 16 MW and the maximum total heat demand in winter is about 49536 MWh.There are two types, five alternative heat sources for this district heating system. The first type is boiler-type, which includes oil-fired boiler, coal-fired boiler, and gas-fired boiler. The second type is heat-pump-type, which includes ground source heat pump and water source heat pump. The significant differences which exist between the two types are as follows: boiler-type heat sources usually burn fossil fuels, and heat-pump-type heat sources do not consume the fossil fuels during the operation; the rated load and return water temperature difference of boiler-type heat sources are higher than those of heat-pump-type heat sources; in the users side, radiation heat exchangers are commonly used for boiler-type heat sources, and coil pipe exchangers are usually adopted for heat-pump-type heat sources. For case studies, operation parameters of the five heat sources are listed in Table4.Table 4
The operation parameters of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
Rated heat capacity (MW)
4.2
4.2
5.6
2.8
2.8
Temperature difference between supply water and return water (°C)
20
20
20
10
10
Energy source
Oil/electric power
Gas/electric power
Coal/electric power
Geothermy/electric power
Water/electric power
Fuel consumption (Kg/h)
393
480 (m3/h)
1120
—
—
Rated power consumption (Kw)
18.5
18.5
12.5
757
660
Power consumption of accessory equipment (Kw)
37
37
65
85
90


### 4.1. Economic Benefit Evaluation

Initial construction cost and operation cost are two key factors affecting the economic benefit of each alternative directly. The sample data for the calculation of the economic benefit index comes from the completed engineering projects and the related literatures [15–17].For boiler-type heat sources, construction cost of main body equipment, such as boiler and turbine, accounts for a large percentage of initial investment. For heat-pump-type heat sources, the initial investment mainly concentrates in the hot well, the underground heat exchanger, and so forth. The initial construction cost of each alternative is shown in Table5.Table 5
The initial construction investments of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
Construction cost for terminal equipment (104 RMB)
1158.03
1158.03
1158.03
2427.46
2427.46
Construction cost for pipeline (104 RMB)
150
150
150
230
230
Construction cost for heat station and accessory equipment (104 RMB)
330
310
400
2618
1880
Total investment cost for initial construction (104 RMB)
1638.03
1618.03
1708.03
5275.46
4537.46Annual operation and management cost includes annual fuel cost, annual electric power cost, labor cost, maintenance cost, and annual depreciation cost. For boiler-type heat sources, fuel cost accounts for a large percentage of the annual operation and management cost. Although the heat-pump-type heat sources have no energy consumption during the operation, the huge amount of electric power consumption has to be considered.For the 16 MW peak heat load, the heat source number of each alternative can be determined according to the rated heating capacity shown in Table4. Combining with the system scale, heating period, energy consumption in rated condition, and fuel prices, the cost of maximum primary energy consumption and electric power consumption can be got for alternatives. However, referring to the real heating load in winter, heating supply of the heat source must be changed to meet the user’s needs and the rated output cannot be held throughout day and night. In order to describe the actual energy consumption reasonably, an equivalent energy consumption coefficient (μh reflects the average load rate of main equipment in a heating period) and equivalent electric power consumption coefficient (με reflects the average power consumption rate in a heating period) are introduced here. The real energy consumption of alternatives can be calculated by multiplying the rated maximum energy consumption and the equivalent energy consumption coefficient. By the same way, the electric power consumption can be got. For the heat-pump-type heat sources, there is no consumption of primary energy and the overall energy conversion efficiency is higher; an attenuation coefficient (η=0.75) should be considered when calculating the real electric power consumption [18].A boiler-type 16 MW district heating system usually needs more than 30 employees. Due to the high level of automation, the heat-pump-type system with the same capacity only requires the half number of employees. The wage for one employee is assumed as 24000 RMB per winter and the welfare is 11% of the wage. The maintenance cost in one running period is accounted as 2.5% of the total initial investment. The depreciation of equipment is calculated as 5% of the initial investment. The detailed operation costs of alternatives are displayed in Table6.Table 6
The annual operation costs of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
The number of installed heat sources
4
4
5
6
6
Primary energy consumption (104 kg/year)
267.68
326.92 (m3)
572.13
—
—
Fuel cost (104 RMB/year)
1338.40
745.40
514.90
—
—
Power consumption for main equipment (104 KWh/year)
16.04
16.04
8.13
788.27
643.65
Power cost (104 RMB/year)
9.62
9.62
4.88
442.96
386.19
Power consumption for accessory equipment (104 KWh/year)
32.07
32.07
42.26
106.62
117.01
Power cost of accessory equipment (104 RMB/year)
28.86
28.86
38.034
95.96
105.31
Wages and welfare (104 RMB/year)
7.92
7.92
7.92
3.96
3.96
Maintenance cost (104 RMB/year)
40.95
40.45
42.70
68.81
59.59
Depreciation of equipment (104 RMB/year)
81.9
80.9
85.4
137.63
119.18Equivalent annual annuity method is used to evaluate the economic benefit. Accounting the salvage value after service life of one alternative is 10% of the initial investment, the service life ism=15 years and the yield rate isI = 11% (including the fixed rate I0=7% and the risk rate Ir=4%). Based on (6), (7), the economic benefits of five alternatives are calculated and shown in Table 7.Table 7
The economic benefits of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
R
m (104 RMB/year)
880.34
618.42
536.06
685.83
602.86
P
c
j
0.609
0.887
1
0.782
0.889The sorting result of economic benefits is as follows:coal-fired boiler > water source heat pump > gas-fired boiler > ground source heat pump > oil-fired boiler.Energy consumption and price are the two important factors influencing the economic benefit. The amount of primary energy consumption of coal-fired boiler is the largest in all alternatives. However, the fact of lower coal price in China helps reduce the operation cost in service cycle. Although heat-pump-type systems use new energy resources and avoid the consumption of fossil fuel energy, the much more electric power consumption in operation cycle and the initial investment make their economic benefits lower than traditional boiler-type systems.The prices of oil, coal, natural gas, and electric power are fluctuating with time and location. So the different decision results may be obtained for different engineering projects. At the same time, the selection of heat-pump-type heat source should also depend on the local natural condition.


### 4.2. Environment Benefit Evaluation

During the winter heating period, environmental pollution caused by the district heating system depends on the fuel type and consumption. Here, gaseous pollutants, which influence atmospheric environment directly, are selected as indicator for environment benefit evaluation.The empirical method is usually used for estimating the gaseous pollutions during one heating period. By analyzing and accounting the real data of running process, the pollutant emission coefficients are defined for various fuel combustion processes. Thus, the quantitative relationship between pollutant emission and fuel consumption can be established [19]. In this work, NOx, SO2, CO, CO2, and other gaseous pollutants are taken into account and the detailed emission data is listed in Table 8.Table 8
The annual pollutant emissions of alternatives.
Items
Oil-fired boiler(g/m2)
Gas-fired boiler(g/m2)
Coal-firedboiler (g/m2)
Ground sourceheat pump (g/m2)
Water sourceheat pump (g/m2)
N
O
x
11.97
7.85
40.54
2.18
1.854
SO2
13.10
7.914
12.34
3.84
3.264
CO
0.361
0.016
1.974
0.0236
0.0206
CO2
4159.38
2695.65
2810.13
570.31
484.76
Others
1.87
2.2886
5.097
0.529
0.449Following the method stated in Section2.2, the relative weights of gaseous pollutions are calculated to establish the judgment matrix as Table 9.Table 9
The judgment matrix for environment evaluation index.
Items
N
O
x
SO2
CO
CO2
Others
N
O
x
1
1/3
7
2
8
SO2
3
1
8
3
7
CO
1/7
1/8
1
1/5
4
CO2
1/2
1/3
5
1
6
Others
1/8
1/7
1/4
1/6
1Based on (11) and (12), the relative weight set of environment benefit is acquired for five alternatives:
(16)ωE=0.267,0.448,0.087,0.178,0.02.Furthermore, referencing to the PROMETHEE, the type II preference function is adopted forNOx and SO2, and the type I preference function is used for other gaseous pollutions. According to (2) and (9)–(12), the environment benefits of the five alternatives are calculated as shown in Table 10. Here, the negative value of the environment benefit means that both of the coal-fired boiler and the oil-fired boiler have no advantage in environmental benefit compared to other alternatives.Table 10
The environment benefits of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
ϕ
E
+
(
j
)
1.02
2.2445
0
2.689
2.3405
ϕ
E
-
(
j
)
2.8465
1.398
4
0.9625
0.087
ϕ
E
(
j
)
−1.8265
0.8465
−4
1.7265
2.2535
P
E
(
j
)
—
0.3756
—
0.7661
1The sorting result of environment benefits is as follows:water source heat pump > ground source heat pump > gas-fired boiler > oil-fired boiler = coal-fired boiler.


### 4.3. Technical Benefit Evaluation

According to Section3.3, quantitative technical benefits of alternatives are obtained as Table 11.Table 11
Attribute table of technical benefit*.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
Technical sophistication
4
4
5
3
3
Heat quality
4
4
5
3
3
Safety and suitability
2
3
4
5
5
Local source condition
2
3
2
5
5
Energy utilizing efficiency
3
3
2
5
5
*Five-class scale is used here, and the larger is better.The relative weights of technical attributes are calculated to establish the judgment matrix as Table12.Table 12
The judgment matrix for technical attributes.
Items
Technicalsophistication
Heat quality
Safety andsuitability
Local sourcecondition
Energy utilizingefficiency
Technical sophistication
1
1/2
1/2
5
1/2
Heat quality
2
1
3
4
1/3
Safety and suitability
2
1/3
1
5
1/2
Local source condition
1/5
1/4
1/5
1
1/3
Energy utilizing efficiency
2
3
2
3
1Referring to the calculation process of environment benefit, the relative weight set of technical benefit for five alternatives is got. Consider(17)ωT=0.154,0.257,0.188,0.054,0.347.As discussed in Sections2.1 and 3.3, type I preference function is used here. According to (2), (13)–(14), the values of incoming flow, outgoing flow, and net flow are calculated and listed in Table 13.Table 13
The technical benefits of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
ϕ
T
+
(
j
)
1.169
1.653
1.928
2.114
1.767
ϕ
T
-
(
j
)
2.019
1.589
2.114
1.281
1.628
ϕ
T
(
j
)
−0.85
0.064
−0.186
0.833
0.139
P
T
(
j
)
—
0.076
—
1
0.166The sorting result of technical benefits is as follows:ground source heat pump > water source heat pump > gas-fired boiler > oil-fired boiler = coal-fired boiler.


### 4.4. Comprehensive Evaluation

Based on the above analysis and calculation about economic benefit, environment benefit, and technical benefit, the comprehensive benefit of alternatives can be calculated by AHP with (15). Consider
(18)Pxj=0.277,0.583,0.455,0.797,0.898,
where the weight factor of each benefit is selected as αc,αe,αt=0.455,0.461,0.084. So the final sorting result of five alternatives is as follows:water source heat pump > ground source heat pump > gas-fired boiler > coal-fired boiler > oil-fired boiler.


## 4.1. Economic Benefit Evaluation

Initial construction cost and operation cost are two key factors affecting the economic benefit of each alternative directly. The sample data for the calculation of the economic benefit index comes from the completed engineering projects and the related literatures [15–17].For boiler-type heat sources, construction cost of main body equipment, such as boiler and turbine, accounts for a large percentage of initial investment. For heat-pump-type heat sources, the initial investment mainly concentrates in the hot well, the underground heat exchanger, and so forth. The initial construction cost of each alternative is shown in Table5.Table 5
The initial construction investments of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
Construction cost for terminal equipment (104 RMB)
1158.03
1158.03
1158.03
2427.46
2427.46
Construction cost for pipeline (104 RMB)
150
150
150
230
230
Construction cost for heat station and accessory equipment (104 RMB)
330
310
400
2618
1880
Total investment cost for initial construction (104 RMB)
1638.03
1618.03
1708.03
5275.46
4537.46Annual operation and management cost includes annual fuel cost, annual electric power cost, labor cost, maintenance cost, and annual depreciation cost. For boiler-type heat sources, fuel cost accounts for a large percentage of the annual operation and management cost. Although the heat-pump-type heat sources have no energy consumption during the operation, the huge amount of electric power consumption has to be considered.For the 16 MW peak heat load, the heat source number of each alternative can be determined according to the rated heating capacity shown in Table4. Combining with the system scale, heating period, energy consumption in rated condition, and fuel prices, the cost of maximum primary energy consumption and electric power consumption can be got for alternatives. However, referring to the real heating load in winter, heating supply of the heat source must be changed to meet the user’s needs and the rated output cannot be held throughout day and night. In order to describe the actual energy consumption reasonably, an equivalent energy consumption coefficient (μh reflects the average load rate of main equipment in a heating period) and equivalent electric power consumption coefficient (με reflects the average power consumption rate in a heating period) are introduced here. The real energy consumption of alternatives can be calculated by multiplying the rated maximum energy consumption and the equivalent energy consumption coefficient. By the same way, the electric power consumption can be got. For the heat-pump-type heat sources, there is no consumption of primary energy and the overall energy conversion efficiency is higher; an attenuation coefficient (η=0.75) should be considered when calculating the real electric power consumption [18].A boiler-type 16 MW district heating system usually needs more than 30 employees. Due to the high level of automation, the heat-pump-type system with the same capacity only requires the half number of employees. The wage for one employee is assumed as 24000 RMB per winter and the welfare is 11% of the wage. The maintenance cost in one running period is accounted as 2.5% of the total initial investment. The depreciation of equipment is calculated as 5% of the initial investment. The detailed operation costs of alternatives are displayed in Table6.Table 6
The annual operation costs of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
The number of installed heat sources
4
4
5
6
6
Primary energy consumption (104 kg/year)
267.68
326.92 (m3)
572.13
—
—
Fuel cost (104 RMB/year)
1338.40
745.40
514.90
—
—
Power consumption for main equipment (104 KWh/year)
16.04
16.04
8.13
788.27
643.65
Power cost (104 RMB/year)
9.62
9.62
4.88
442.96
386.19
Power consumption for accessory equipment (104 KWh/year)
32.07
32.07
42.26
106.62
117.01
Power cost of accessory equipment (104 RMB/year)
28.86
28.86
38.034
95.96
105.31
Wages and welfare (104 RMB/year)
7.92
7.92
7.92
3.96
3.96
Maintenance cost (104 RMB/year)
40.95
40.45
42.70
68.81
59.59
Depreciation of equipment (104 RMB/year)
81.9
80.9
85.4
137.63
119.18Equivalent annual annuity method is used to evaluate the economic benefit. Accounting the salvage value after service life of one alternative is 10% of the initial investment, the service life ism=15 years and the yield rate isI = 11% (including the fixed rate I0=7% and the risk rate Ir=4%). Based on (6), (7), the economic benefits of five alternatives are calculated and shown in Table 7.Table 7
The economic benefits of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
R
m (104 RMB/year)
880.34
618.42
536.06
685.83
602.86
P
c
j
0.609
0.887
1
0.782
0.889The sorting result of economic benefits is as follows:coal-fired boiler > water source heat pump > gas-fired boiler > ground source heat pump > oil-fired boiler.Energy consumption and price are the two important factors influencing the economic benefit. The amount of primary energy consumption of coal-fired boiler is the largest in all alternatives. However, the fact of lower coal price in China helps reduce the operation cost in service cycle. Although heat-pump-type systems use new energy resources and avoid the consumption of fossil fuel energy, the much more electric power consumption in operation cycle and the initial investment make their economic benefits lower than traditional boiler-type systems.The prices of oil, coal, natural gas, and electric power are fluctuating with time and location. So the different decision results may be obtained for different engineering projects. At the same time, the selection of heat-pump-type heat source should also depend on the local natural condition.


## 4.2. Environment Benefit Evaluation

During the winter heating period, environmental pollution caused by the district heating system depends on the fuel type and consumption. Here, gaseous pollutants, which influence atmospheric environment directly, are selected as indicator for environment benefit evaluation.The empirical method is usually used for estimating the gaseous pollutions during one heating period. By analyzing and accounting the real data of running process, the pollutant emission coefficients are defined for various fuel combustion processes. Thus, the quantitative relationship between pollutant emission and fuel consumption can be established [19]. In this work, NOx, SO2, CO, CO2, and other gaseous pollutants are taken into account and the detailed emission data is listed in Table 8.Table 8
The annual pollutant emissions of alternatives.
Items
Oil-fired boiler(g/m2)
Gas-fired boiler(g/m2)
Coal-firedboiler (g/m2)
Ground sourceheat pump (g/m2)
Water sourceheat pump (g/m2)
N
O
x
11.97
7.85
40.54
2.18
1.854
SO2
13.10
7.914
12.34
3.84
3.264
CO
0.361
0.016
1.974
0.0236
0.0206
CO2
4159.38
2695.65
2810.13
570.31
484.76
Others
1.87
2.2886
5.097
0.529
0.449Following the method stated in Section2.2, the relative weights of gaseous pollutions are calculated to establish the judgment matrix as Table 9.Table 9
The judgment matrix for environment evaluation index.
Items
N
O
x
SO2
CO
CO2
Others
N
O
x
1
1/3
7
2
8
SO2
3
1
8
3
7
CO
1/7
1/8
1
1/5
4
CO2
1/2
1/3
5
1
6
Others
1/8
1/7
1/4
1/6
1Based on (11) and (12), the relative weight set of environment benefit is acquired for five alternatives:
(16)ωE=0.267,0.448,0.087,0.178,0.02.Furthermore, referencing to the PROMETHEE, the type II preference function is adopted forNOx and SO2, and the type I preference function is used for other gaseous pollutions. According to (2) and (9)–(12), the environment benefits of the five alternatives are calculated as shown in Table 10. Here, the negative value of the environment benefit means that both of the coal-fired boiler and the oil-fired boiler have no advantage in environmental benefit compared to other alternatives.Table 10
The environment benefits of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
ϕ
E
+
(
j
)
1.02
2.2445
0
2.689
2.3405
ϕ
E
-
(
j
)
2.8465
1.398
4
0.9625
0.087
ϕ
E
(
j
)
−1.8265
0.8465
−4
1.7265
2.2535
P
E
(
j
)
—
0.3756
—
0.7661
1The sorting result of environment benefits is as follows:water source heat pump > ground source heat pump > gas-fired boiler > oil-fired boiler = coal-fired boiler.


## 4.3. Technical Benefit Evaluation

According to Section3.3, quantitative technical benefits of alternatives are obtained as Table 11.Table 11
Attribute table of technical benefit*.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
Technical sophistication
4
4
5
3
3
Heat quality
4
4
5
3
3
Safety and suitability
2
3
4
5
5
Local source condition
2
3
2
5
5
Energy utilizing efficiency
3
3
2
5
5
*Five-class scale is used here, and the larger is better.The relative weights of technical attributes are calculated to establish the judgment matrix as Table12.Table 12
The judgment matrix for technical attributes.
Items
Technicalsophistication
Heat quality
Safety andsuitability
Local sourcecondition
Energy utilizingefficiency
Technical sophistication
1
1/2
1/2
5
1/2
Heat quality
2
1
3
4
1/3
Safety and suitability
2
1/3
1
5
1/2
Local source condition
1/5
1/4
1/5
1
1/3
Energy utilizing efficiency
2
3
2
3
1Referring to the calculation process of environment benefit, the relative weight set of technical benefit for five alternatives is got. Consider(17)ωT=0.154,0.257,0.188,0.054,0.347.As discussed in Sections2.1 and 3.3, type I preference function is used here. According to (2), (13)–(14), the values of incoming flow, outgoing flow, and net flow are calculated and listed in Table 13.Table 13
The technical benefits of alternatives.
Items
Oil-firedboiler
Gas-firedboiler
Coal-firedboiler
Ground sourceheat pump
Water sourceheat pump
ϕ
T
+
(
j
)
1.169
1.653
1.928
2.114
1.767
ϕ
T
-
(
j
)
2.019
1.589
2.114
1.281
1.628
ϕ
T
(
j
)
−0.85
0.064
−0.186
0.833
0.139
P
T
(
j
)
—
0.076
—
1
0.166The sorting result of technical benefits is as follows:ground source heat pump > water source heat pump > gas-fired boiler > oil-fired boiler = coal-fired boiler.


## 4.4. Comprehensive Evaluation

Based on the above analysis and calculation about economic benefit, environment benefit, and technical benefit, the comprehensive benefit of alternatives can be calculated by AHP with (15). Consider
(18)Pxj=0.277,0.583,0.455,0.797,0.898,
where the weight factor of each benefit is selected as αc,αe,αt=0.455,0.461,0.084. So the final sorting result of five alternatives is as follows:water source heat pump > ground source heat pump > gas-fired boiler > coal-fired boiler > oil-fired boiler.


## 5. Results and Discussion

Considering three benefits (the economic benefit, the environment benefit, and the technical benefit), the water source heat pump is selected by the proposed hierarchical decision-making scheme as the best heat source of the district heating system for a community in Beijing. Compared to the ground source heat pump that can only be applied at geothermal location, the water source heat pump has a variety of energy sources (such as the groundwater, the surface water, the industrial waste water, or the living waste water) and therefore is more likely to be accepted.In boiler-type alternatives, gas-fired boiler has the best comprehensive benefit. In the geothermal and water scarcity area, gas-fired boiler will be the preference selection for the district heating system. Although the coal-fired boiler has the best economic benefit, it still cannot be selected as an ideal solution for its worse environment benefit.This work shows that the integration of the PROMETHEE and the AHP makes the qualitative and quantitative attributions easy to be treated together in the decision-making process. The case study further indicates that, in broader areas of industrial and civil engineering, the proposed hierarchical decision-making scheme has promising developing trend.



---

*Source: 594862-2014-11-16.xml*