Next Article in Journal
Performance Assessment of a Building-Integrated Photovoltaic Thermal System in a Mediterranean Climate—An Experimental Analysis Approach
Next Article in Special Issue
A Multi-Objective Approach toward Optimal Design of Sustainable Integrated Biodiesel/Diesel Supply Chain Based on First- and Second-Generation Feedstock with Solid Waste Use
Previous Article in Journal
SPEAR (Solar Pyrolysis Energy Access Reactor): Theoretical Design and Evaluation of a Small-Scale Low-Cost Pyrolysis Unit for Implementation in Rural Communities
Previous Article in Special Issue
Numerical Investigation of Gravity-Driven Granular Flow around the Vertical Plate: Effect of Pin-Fin and Oscillation on the Heat Transfer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Framework for Design and Operation Optimization for Utilizing Low-Grade Industrial Waste Heat in District Heating and Cooling †

1
State Key Laboratory of Heavy Oil Processing, China University of Petroleum, Beijing 102249, China
2
School of Chemical Engineering & Technology, Xi’an Jiaotong University, Xi’an 710049, China
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in Chemical Engineering Transactions for The 23rd Conference on Process Integration, Modelling and Optimisation for Energy Saving and Pollution Reduction (PRES’20), Xi’an, China, 17–21 August 2020; Volume 81, pp. 511–516.
Energies 2021, 14(8), 2190; https://doi.org/10.3390/en14082190
Submission received: 11 March 2021 / Revised: 3 April 2021 / Accepted: 12 April 2021 / Published: 14 April 2021

Abstract

:
In the process industry, a large amount of low-grade waste heat is discharged into the environment. Furthermore, district heating and cooling systems require considerable low-grade energy. The integration of the two systems has great significance for energy saving. Because the energy demand of consumers varies in periods, the design and operation of an industrial waste heat recovery system need to match with the fluctuations of district energy demand. However, the impact of the periodic changes on the integration schemes are not considered enough in existing research. In this study, a framework method for solving above problem is proposed. Industrial waste heat was integrated with a district heating and cooling system through a heat recovery loop. A three-step mathematical programming method was used in design and operation optimization for multiperiod integration. A case study was conducted, and the results show that the multiperiod optimization method can bring significant benefits to the system. By solving the mixed integer nonlinear programming model, the optimal operation plans of the integration in different periods can be obtained.

1. Introduction

In the process industry, a large amount of waste heat usually exists. According to the different temperature levels, industrial waste heat can be divided into high-grade waste heat, medium-grade waste heat, and low-grade waste heat. Although there are no determinate boundaries between these three levels, it is acceptable that 650 °C and 230 °C are the dividing points [1]. As an effective method to reduce energy consumption, Heat Integration (HI) has been studied in depth in the past 40 years [2]. Such methods led to significant energy saving in the process industry around the world. The very first work involving HI was about heat exchanger network (HEN) synthesis. There are two mainstream methods for the design of HENs: pinch analysis (PA) and mathematical programming (MP). The former is based on thermodynamic principles and a graphical method [3], and the core of the latter is superstructure formulation, whereby the generalized superstructure [4] and the stage-wise superstructure are widely used [5]. In the process industry, high-grade and medium-grade waste heat can be effectively utilized by HENs. In addition, these two grades of waste heat can also be recovered for steam production by a waste heat boiler [6]. However, since heat sinks that need low-grade heat are rare, most low-grade waste has to discharge into the environment, causing serious energy waste [7]. In Europe, the amount of low-grade heat consumed by the chemical industry is up to 1142 PJ every year, while 20–50% of it eventually becomes waste heat [8]; According to the estimation method of Fang et al. [9], 760 GJ (a coal equivalent of nearly 260 million tons) of low-grade waste heat is lost in northern China every year.
Due to the above-mentioned reasons, effective utilization of low-grade waste heat has become a hotspot of research. Huang et al. [10] comprehensively summarized the technologies of low-grade waste heat utilization and classified these technologies into two categories: “heat upgrading” and “heat converting”. The heat pump is the most popular “heat upgrading” technology. It takes high-grade energy as the driving force to convert low-grade heat into higher-grade heat. The “heat converting” technology can transform waste heat into mechanical work via thermodynamic cycles, and then electricity can be generated. There are three kinds of thermodynamic cycles that are usually used: organic Rankine cycle [11], Kalina cycle [12], and trilateral cycle [13]. Although these technologies have been researched intensively, the efficiency of low-grade waste heat utilization is still relatively low due to strict thermodynamic constraints.
In recent years, the integration of low-grade industrial waste heat into district heating and cooling (DHC) systems has become an attractive option. The DHC system delivers steam, hot water, or cold water from a central energy station to consumers via large-scale underground pipeline networks. DHC systems have developed in three generations. The media adopted for district heating (DH) are steam, pressurized hot water, and atmospheric hot water. The method of district cooling (DC) has changed from transporting refrigerant to delivering cold water. Technologies such as natural cooling from lakes and absorption refrigeration have gradually been applied [14]. Therefore, more non-fossil-fuel energy sources can be integrated with the DHC systems in the future [15]. In the fourth-generation DHC system, utilization of low-grade industrial waste is more promising. The advantages of integrating low-grade industrial waste heat into the DHC system include improving energy efficiency, reducing fossil fuel combustion and greenhouse gas emissions, and saving a large amount of cooling water.
A number of studies were reported on utilizing industrial waste heat in DH systems. Tveit et al. [16] proposed a system framework for analyzing the integration of large-scale energy systems. They designed a flexible heat exchanger network with bypasses for the integration of four waste heat streams in a paper plant into a DH system. Their work was constructive in the early research on the integration of low-grade waste heat with DH systems, but they did not utilize waste heat for district cooling. Kapil et al. [17] proposed a method to identify the quantity and quality of waste heat from the total site that can be utilized in a DH system. Their method is based on economic benefit evaluation. They emphasized the influence of distance factors, operating conditions, and economic parameters on the economic benefits of the integration, but a concrete scheme for the integration including HEN was not designed. Fang et al. [9] combined the low-grade waste heat of two plants (a copper smeltery and a cement plant) for district heating. A total of 1896 m3/h of hot water at about 80 °C can be obtained using their program of waste heat recovery. They discussed the three key issues related to using low-grade industrial waste heat for district heating in their subsequent work [18]: (1) “tangency technology” (a method similar to PA) was applied for the collection of multiple-grade waste heat sources; (2) technologies such as heat pump were suggested for long-distance delivery of waste heat; (3) a combined heat and power (CHP) plant was recommended for peak shaving. Wang et al. [19] improved the “tangency technology” using the exergy analysis method. Fitó et al. [20] discussed two design directions for the design of a waste heat recovery system for DH: energy-based and exergy-based optimal design. The methods of the above two studies can provide a clear energy target for integration, but the cost target and economic-related factors were not discussed enough. Recently, Khosravi et al. [21] suggested that the waste heat from data centers and fifth generation (5G) smart poles can be recovered for DH networks, and they confirmed the feasibility of utilizing low-grade waste heat in future urban smart thermal grids.
Another effective method to utilize low-grade industrial waste heat is absorption refrigeration [12]. Through this technology, low-grade waste heat can be used in DC systems. A comprehensive review of this technology can be found in the work of Srikhirin et al. [22]. According to the total site Heat Integration (TSHI), Liew et al. [23] used low-grade waste heat in the site to drive the absorption chiller. Combined with the electric compression chiller, a lower cost was achieved in the case study.
Although many studies were presented on the utilization of low-grade industrial waste heat in DHC systems, almost all of them separately discussed waste heat utilization in DH systems or DC systems. In practice, both the heating demand and the cooling demand of consumers vary in different periods; thus, DH and DC systems should be considered simultaneously. Therefore, different heat recovery schemes for industries and different energy supply schemes in DHC systems are required in different periods. As the important methods of HI, PA and exergy analysis have been widely applied, but it is difficult for them to deal with complex problems involving multiple parameters and multiple dimensions. In cases with multiperiod HI, the compatibility of operation schemes in different periods is important. Therefore, the above methods are not suitable. To solve the above problem, a three-step method based on mathematical programming is proposed in this paper. This method is an extension of our previous work in a conference paper [24]. In previous work, a scheme for DHC systems was not involved. By applying this method, the design and operation optimization of the heat recovery system and DHC can be obtained, and the operating parameters in each period can be determined.

2. Problem Statement

In the production area of a plant, there are several hot streams that need to be cooled. The low-grade waste heat of hot streams is originally discharged into the environment by cold utility. In this study, waste heat was integrated with the DHC system through a heat recovery loop (HRL), as shown in Figure 1. Hot water was chosen as the intermediate fluid. In order to establish this indirect HI, a small station for placing absorption chillers was built on the empty land near the plant. There was a long distance between the production area and the station. A pump was used to transport hot water in the HRL. This is called the “waste heat recovery (WHR) system” in this paper. After recovering the waste heat, hot water with a higher temperature is received, which can be directly used in district heating or to drive the operation of the absorption chiller. This is referred to as the “district energy supply (DES) system” in this paper.
A number of consumers with heating or cooling demands exist in the region around the plant. Normally, the amount of waste heat from industries cannot satisfy the demands of all consumers; thus, the heating and cooling duties of the plant were determined in the first step of this study. In the WHR system, not all the waste heat is economic for recovery considering costs related to pumping, piping, and heat exchange. Only when the benefit is greater than the cost should the waste heat of the hot stream be recovered. Therefore, it is necessary to consider both cost and benefit in the objective function, as shown in Equation (1). By formulating and solving the mixed integer nonlinear programming (MINLP) model, a two-period HEN structure and the heating and cooling potential can be obtained. The mathematical model is introduced in Section 3.
min {TAC = Ccu + Cexchanger + Cstation + Cpipe,loop + CpumpIDHC},
where TAC is the total annual cost in the first step, Ccu is the cost of cold utility that is still needed, Cexchanger is the capital cost of heat exchangers installed to recover waste heat, Cstation is the capital cost of the station, including the cost of the absorption chiller, Cpipe,loop and Cpump are the capital cost of pipeline and the cost of pump in the HRL, And IDHC is the income of the plant by district heating and cooling. In all cases, capital costs are annualized.
Since each consumer has a different heating/cooling demand and distance to the plant, it is necessary to determine the suitable consumer for providing energy. Therefore, the second step of this study was to determine the match between industry and consumers. In a district energy system, the peak periods of heating and cooling are winter and summer. Some types of buildings (such as hospitals and data centers) have cooling demand throughout the year. Considering the above situation, the low-grade waste heat resource of the plant was used for district heating only in winter and district cooling throughout the year. It should be noted that the cooling demand is the highest in summer and lowest in winter, while the cooling demands in spring and autumn were assumed to be the same in this work.
There are two major costs in a DES system: the capital cost of pipelines between the energy station and consumers, and the capital cost of the absorption chiller. Both are directly related to the capacity of energy supply [25] (the cost of the absorption chiller is included in the cost of the station Cstation). In order to determine the optimal scheme of energy supply, the objective function of the second step was to maximize the total profit of DES system, as shown in Equation (2). To avoid complexity in calculation, in this step, the operation in industry was not optimized.
min {TAP = IDHCCpipe,DHCCstation},
where TAP is the total annual profit of the DES system, and Cpipe,DHC is the capital cost of pipelines in the DHC system.
After the second step, the total heating or cooling capacity of the DES system in each period is known, but this does not mean that the operating scheme of the WHR system is fixed. Chang et al. [26] emphasized that there is a tradeoff in the HRL between the heat exchanger cost and the distance-related cost when the capacity of heat recovery is determined. Therefore, the third step of this study was to optimize the multiperiod operation schemes of the WHR system under a certain energy supply task. In this step, the energy supply constraint in each period was added into the model of the first step. The items Cstation and IDHC that were determined in the second step were deleted in the objective function, which is expressed by Equation (3).
min {TAC’ = Ccu + Cexchanger + Cpipe,loop + Cpump},
where TAC’ is the total annual cost in the third step.
In summary, this work was carried out in three steps: (1) determination of the heating and cooling potential of a plant; (2) determination of the regional energy supply objects in each period; (3) design and operation optimization of the multiperiod HRL under the definite energy supply task.

3. Mathematical Formulation

In this paper, hot streams of the plant are defined by index i, II, I = {H1, H2, …, Hm}, periods are defined by index s, sS, S = {spring, summer, autumn, winter}, and consumers are defined by index j, jJ, J = {N1, N2, …, Nn}.

3.1. The First Step

In the DHC system, peak periods of heating and cooling are winter and summer. Therefore, for the convenience of expression, s = winter represents the period with the maximum heating potential, and s = summer represents the period with the maximum cooling potential.

3.1.1. Model of Absorption Chiller

Although dozens of compounds have been suggested as absorbents or refrigerants, the commonly used working fluid pairs of absorption refrigeration are still H2O/NH3 and LiBr/H2O [22]. Compared with LiBr/H2O, H2O/NH3 requires a higher heat source temperature, while the double-effect and triple-effect absorption chillers require a higher heat source temperature than single-effect ones [27]. In order to utilize low-grade waste heat as much as possible, a single-effect LiBr-H2O absorption chiller was selected in this study.
The COP (coefficient of performance) of the chiller is an important variable in absorption refrigeration. It is acceptable that the electricity consumed by the chiller can be ignored, whereby the cooling output of a single-effect LiBr/H2O refrigerator is calculated using Equation (4) [28].
Q s c o l d = C O P s · Q s r e c o v e r ,
where Q s c o l d is the cooling output of the chiller, and Q s r e c o v e r is the capacity of the waste heat recovered. As mentioned, absorption refrigeration works in three periods: spring, summer, and autumn (only s = summer is involved in Section 3.1); thus, the subscript of the periodic element was added to the related symbols in this section.
Wang et al. [28] explored the variation in COP of a single-effect LiBr/H2O refrigerator at three evaporator temperatures (5 °C, 8 °C, and 11 °C). Pressurized hot water at 100–150 °C was used as the heat source in their study. In this paper, the temperature of the evaporator was selected as 5 °C. In this condition, there is a nonlinear relation between COPs and the temperature of hot water entering the chiller T s w , c h i l l e r , i n [28]; accordingly, piecewise linearization should be used to describe this relationship, as expressed as Equation (5).
In mathematical programming, Equation (5) is replaced by Equations (6)–(12).
In the simulation conditions of Wang et al. [28], there is a one-to-one correspondence between T s w , c h i l l e r , i n and the temperature of hot water leaving the chiller T s w , c h i l l e r , o u t . According to their data, the relationship between T s w , c h i l l e r , i n and T s w , c h i l l e r , o u t is shown in Figure 2. It can be seen that there is an obvious linear relationship between T s w , c h i l l e r , i n and T s w , c h i l l e r , o u t . Equation (13) can be obtained by linear fitting.
C O P s = 0.0437 · T s w , c h i l l e r , i n 4.217   100     T s w , c h i l l e r , i n     110 , 0.0082 · T s w , c h i l l e r , i n 0.312   110     T s w , c h i l l e r , i n     120 , 0.0043 · T s w , c h i l l e r , i n 0.156   120     T s w , c h i l l e r , i n     130 , 0.0018 · T s w , c h i l l e r , i n 0.488   130     T s w , c h i l l e r , i n     150 ,
y a s · 100 T s w , c h i l l e r , i n , a y a s · 110 ,
y b s · 110 T s w , c h i l l e r , i n , b y b s · 120 ,
y c s · 120 T s w , c h i l l e r , i n , c y c s · 130 ,
y d s · 130 T s w , c h i l l e r , i n , d y d s · 150 ,
yas + ybs + ycs + yds = 1,
T s w , c h i l l e r , i n = T s w , c h i l l e r , i n , a + T s w , c h i l l e r , i n , b + T s w , c h i l l e r , i n , c + T s w , c h i l l e r , i n , d ,
C O P s = ( 0.0437 · T s w , c h i l l e r , i n , a y a s · 4.217 ) + ( 0.0082 · T s w , c h i l l e r , i n , b y b s · 0.312 ) + ( 0.0043 · T s w , c h i l l e r , i n , c y c s · 0.156 ) + ( 0.0018 · T s w , c h i l l e r , i n , d y d s · 0.488 ) ,
where yas, ybs, ycs, and yds are binary variables (0 or 1), while T s w , c h i l l e r , i n , a , T s w , c h i l l e r , i n , b , T s w , c h i l l e r , i n , c , and T s w , c h i l l e r , i n , d are the piecewise elements of T s w , c h i l l e r , i n .
T s w , c h i l l e r , o u t = 0.426 · T s w , c h i l l e r , i n + 52.8 .
The investment of the station including the capital cost of the absorption chiller is calculated using Equation (14).
C s t a t i o n = A f · ( b 1 + b 2 · Q s c o l d ) ( s = summer ) ,
where Af is the annual factor, and b1 and b2 are the parameters of the cost of the small station.

3.1.2. Model of HEN

In the WHR system, the stage-wise superstructure [5] is adopted to formulate the heat exchanger network, which is shown in Figure 3. The parameters T i h s , i n , T h s , o u t I, and F i h s are the initial temperature, target temperature, and flow capacity of hot stream i, respectively. The variables T s w , H E N , i n , T s w , H E N , o u t , and m s w , t o t a l are the inlet temperature, outlet temperature, and mass flow rate of hot water in period s. m i , k , s w , s p l i t and qi,k,s represent the mass flow rate of hot water and the heat load of exchanger for hot stream i at stage k in period s. T i , k , s h s and T k , s w , H E N represent the temperature of hot stream i and hot water at hot end of stage k. q i , s c u is the heat load of cooler between hot stream i and cold utility.
The inlet and outlet temperatures of hot streams and hot water in HEN are defined by Equations (15)–(17).
T i , 1 , s h s = T i h s , i n .
T 1 , s w , H E N = T s w , H E N , o u t .
T NT + 1 , s w , H E N = T s w , H E N , i n .
In current DH systems, the temperature of atmospheric hot water sent to consumers ( T s w , H E N , o u t ) is generally required to be higher than 70 °C [14]. The temperature of water sent back to the station ( T s w , H E N , i n ) is always not lower than 40 °C [15]. In order to reduce the solution difficulty, the latter temperature was fixed to be 40 °C rather than being a variable in this work. Therefore, the constraints for the temperature of hot water in the WHR system in winter are expressed by Equations (18) and (19).
70 T s w , H E N , o u t 100   ( s = winter ) .
T s w , H E N , i n = 40   ( s = winter ) .
In summer, the inlet and outlet temperature of hot water in the HEN are equal to the temperature leaving and entering the absorption chiller; thus, Equation (20) emerges on the basis of Equation (13).
T s w , H E N , i n = 0.426 · T s w , H E N , o u t + 52.8   ( s = summer ) .
The total capacity of waste heat recovered is calculated by Equation (21). cpw represents the specific heat capacity of water.
Q s r e c o v e r = c p w · m s w , t o t a l · ( T s w , H E N , o u t T s w , H E N , i n ) .
Energy and mass balance constraints for hot water and hot streams are expressed by Equations (22)–(27).
q i , k , s = F i h s · ( T i , k , s h s T i , k + 1 , s h s ) .
q i , k , s = c p w · m i , k , s w , s p l i t · ( T k , s w , H E N T k + 1 , s w , H E N ) .
m s w , t o t a l = Σ i   m i , k , s w , s p l i t .
Σ i   q i , k , s = c p w · m s w , t o t a l · ( T k , s w , H E N T k + 1 , s w , H E N ) .
q i , s c u = F i h s · ( T N T + 1 , s h s T i h s , o u t ) .
F i h s · ( T i h s , i n T i h s , o u t ) = Σ k   q i , k , s + q i , s c u .
The Big-M constraints for heat loads and temperature approaches are used to determine the existence of heat exchangers, which are expressed by Equations (28)–(30).
q i , k , s z i , k · q i h s ,
d t i , k , s T i , k , s h s T k , s w , H E N + ( 1 z i , k ) · Γ ,
d t i , k + 1 , s T i , k + 1 , s h s T k + 1 , s w , H E N + ( 1 z i , k ) · Γ ,
where zi,k is the binary variable to indicate whether the heat exchanger between hot water and hot stream i at stage k exists. q i h s is the total heat load of hot stream i. dti,k,s is the temperature approach between hot stream i and hot water at hot end of stage k in period s. The parameter Γ is the upper bound for temperature difference, which is 250 °C in this paper. If the heat exchanger exists, the binary variable zi,k is equal to 1; if the heat exchanger does not exist, zi,k can be 0 or 1 in Equations (28)–(30). In order to minimize the cost of heat exchangers, zi,k is forced to be equal to 0 by objective function.
The temperature approaches must be greater than the minimum temperature difference ΔTmin in order to meet the requirement for heat transfer. This relationship is shown by Equation (31).
dti,k,s ≥ ΔTmin.
Equation (32) is employed to calculate the heat transfer area needed for hot water to recover the heat of hot stream i at stage k in period s ( A i , k , s p e r i o d ).
A i , k , s p e r i o d = q i , k , s · ( 1 / h w + 1 / h i h s ) / L M T D i , k , s ,
where LMTDi,k,s is the logarithmic mean temperature difference of this heat exchanger. hw and h i h s are the film heat transfer coefficients of hot water and hot stream i.
In order to reduce the solution difficulty, LMTDi,k,s is calculated using the Chen approximation [29], which is expressed by Equation (33).
LMTDi,k,s = [dti,k,s ·dti,k + 1,s · (dti,k,s + dti,k+1,s)/2]1/3.
In multiperiod HENs, the design area of exchanger ( A i , k d e s i g n ) is the maximum area needed in all periods [30], which is expressed by Equation (34).
A i , k d e s i g n A i , k , s p e r i o d .
Equations (35) and (36) are used to calculate the cost of cold utility and the capital cost of heat exchangers, respectively.
C c u = Σ i Σ s u c u · ( t s p e r i o d / t t o t a l ) · q i , s c u ,
C e x c h a n g e r = A f · [ α · Σ i Σ k   z i , k + β · Σ i Σ k   A i , k d e s i g n ] γ ,
where ucu is the unit cost of cold utility. t s p e r i o d and ttotal are the time of period s and total year. α, β, and γ are the cost parameters of the heat exchanger.

3.1.3. Model of Pipeline

There are two types of pipelines in the HRL: a long pipeline between production area and station and short pipelines in the production area. This work only considers the long pipeline, while short pipelines are ignored. Peters et al. [31] provided an equation for designing the inner diameter of the pipeline using an appropriate velocity. However, the flow rate of hot water varies in different periods, resulting in a variable velocity. In this study, the inner diameter of the pipeline takes the maximum value in each period. It is indicated that hot water flows at the velocity recommended by Peters et al. only in the period with the largest flow rate [31]. In other periods, hot water flows more slowly. Because of the high temperature, desalinated water should be used as the intermediate fluid; hence, the situation of fouling is ignored in this study. Therefore, the inner diameter (Din) of the pipeline between production area and station is expressed by Equation (37) according to the method of Peters et al. [31].
D i n = 0.363 · ( m s w , t o t a l / ρ s ) 0.45 · ρ s 0.13 ,
where ρs is the density of hot water. It should be declared that the densities of atmospheric water and pressurized water are different.
Because the pressurized hot water is used for absorption refrigeration, the schedule 80 steel pipe is employed in HRL. The cost of the pipeline per unit length in HRL (Pculloop) was calculated using the method of Stijepovic et al. [32], which is shown by Equations (38)–(40).
Dout = 1.101 · Din + 0.006349.
Wt = 1330 · Din2 + 75.18 · Din + 0.9268.
Pculloop = a1 · Wt + a2 · Dout0.48 + a3 + a4 · Dout.
The total cost of the pipeline in the HRL is calculated using Equation (41). Lloop is the distance between station and production area.
Cpipe,loop = Af · 2 ·Lloop · Pculloop

3.1.4. Model of Pump

In the HRL, the pressure drop of hot water should be overcome by the pump. The total pressure drop of hot water Δ p s t o t a l is expressed in Equation (42).
Δ p s t o t a l = Δ p s d i s t a n c e + Δ p s H E N ,
where Δ p s d i s t a n c e and Δ p s H E N are the pressure drop of hot water in the long pipeline and HEN in period s.
Δ p s d i s t a n c e is calculated using the Fanning formula. The detailed process of calculation is not introduced here; a calculation case can be found in [26].
The pressure drop of hot water in the tube side of each heat exchanger (Δ p i , k , s t u b e ) is calculated using the method of Soltani et al. [33], which is shown by Equations (43) and (44).
K t i , k , s = ( Φ 4.5 · d i 0.5 · μ s 11 / 6 ) · ( d i / d e ) / [ 0.023 2.5 · m i , k , s w , s p l i t · ρ s · k 7 / 3 · ( c p w ) 7 / 6 ] .
Δ p i , k , s t u b e = K t i , k , s · A i , k , s p e r i o d · ( h w ) 3.5 .
In the HEN based on the stage-wise superstructure, the pipes of hot water are in a series–parallel structure, whereas stages are connected in series. Therefore, the total pressure drop of hot water in the HEN is equal to the sum of pressure drop in each stage (Δ p k , s s t a g e ); this relationship is shown in Equation (45).
Δ p s H E N = Σ k   Δ p k , s s t a g e .
According to the principle of fluid flow, the pressure drops of fluid on branches are equal in a stable condition. Therefore, the pressure drop of hot water at a stage is equal to the maximum pressure drop of each branch, as shown in Equation (46).
Δ p k , s s t a g e Σ k   Δ p i , k , s t u b e .
The power of pump in each period ( P s p e r i o d ) and the rated power of pump (Prate) are expressed by Equations (47) and (48), respectively.
P s p e r i o d = ( m s w , t o t a l · Δ p s t o t a l ) / ρ s .
P r a t e P s p e r i o d
The cost of pump Cpump is equal to the sum of capital cost and operation cost, which are calculated using Equations (50) and (51), respectively.
Cpump = Cpump,cap + Cpump,op,
Cpump,cap = Af · [c1 + c2 · (Prate)c3],
C p u m p , c a p = u e · Σ s ( P s p e r i o d · t s p e r i o d ) / ( 1000 · η ) ,
where c1, c2, and c3 are the parameters of capital cost of the pump, ue is the unit price of electricity, and η is the efficiency of pump.

3.1.5. Calculation of Income

The maximum heating and cooling potentials (QDH,max and QDC,max) are represented by Equations (52) and (53). The total income of the plant in these two periods is calculated by Equation (54). uDC and uDH are the unit prices of heating and cooling respectively.
Q D H , m a x = Q s = w i n t e r r e c o v e r .
Q D C , m a x = Q s = s u m m e r r e c o v e r .
I D H C = ( u D C · Q D C , m a x · t s = summer p e r i o d ) + ( u D H · Q D H , m a x · t s = winter p e r i o d ) .

3.2. The Second Step

Based on the heating and cooling potential of the plant obtained in the first step, the classic “knapsack problem” model in integer programming was used to determine the actual energy supply task in each period. The “knapsack problem” can be expressed as follows: a traveler has to choose among items with different weights and values, and then put the selected items into a knapsack with a certain capacity, where the maximum total value of selected items is expected. In this study, the plant is the “knapsack”, and consumers are the “items” to be selected.
In the second step, all the period-related elements were included. In a DHC system, the loss of heat and cold are often more than 10% in the distribution network [15]. The difference between the actual capacity of heat/cold supplied to consumer j from plant ( E j , s s u p p l y ) and the demand of consumer j ( E j , s n e e d ) cannot be ignored. The loss of energy per unit distance in DHC system can be taken as 1%/km [17]. Equation (55) is used to express the relationship between E j , s s u p p l y and E j , s n e e d . L j c o n s u m e r is the distance between the plant and consumer j.
E j , s s u p p l y = E j , s n e e d / [ 0.99 ( L j c o n s u m e r / 1000 ) ] .
According to the “knapsack problem” model, the binary variable y j , s s u p p l y is used for the selection of consumers, which is shown in Equation (56).
Σ j   y j , s s u p p l y · E j , s s u p p l y Q D C , m a x   s = spring , summer , autumn , Q D C , m a x s = winter .
In a DHC system, the temperature difference between supply and return of water in the DC system is much smaller than that in the DH system, which means that the diameters of pipelines in the DC system is often significantly larger than that in the DH system [15]. Therefore, the pipelines cannot be shared by these two systems. However, Best et al. [34] proposed that the flow rate of water transported by the pipeline with a certain diameter can be within a range in a DHC system; that is, the minimum flow rate should be greater than 0.475 times the maximum flow rate. This recommendation makes it possible to combine the pipeline for district cooling in spring/autumn and summer. In this study, the supply and return temperature of cold water in DC system were set to 7 °C and 12 °C. Therefore, the flow rate of cold water supplied to each consumer is proportional to its demand of cooling. When the demands of a consumer in periods are given, it can be known whether it is feasible that the DC pipeline can be shared in spring/autumn and summer. This possibility is expressed by the binary parameter y j c o m b i n e , which is shown in Equation (57).
y j c o m b i n e = 0   ( E j , s = s p r i n g s u p p l y     0.475   ·   E j , s = s p r i n g s u p p l y = spring ) , 1   ( E j , s = s p r i n g s u p p l y     0.475   ·   E j , s = s p r i n g s u p p l y = spring ) .
When consumer j is cooled by the plant in spring/autumn, if it is also cooled in summer, and the condition y j c o m b i n e = 1 is satisfied, the pipeline of consumer j can be shared in spring, summer, and autumn. Therefore, on the condition that the combination of DC pipeline is considered, the total capital cost of pipelines in the DHC system Cpipe, DHC can be calculated using Equations (58)–(61).
Cpipe, DHC = Cpipe, DHC, wiNTer + Cpipe, DHC, summer + Cpipe, DHC, spring,
C p i p e , D H C , w i n t e r = Σ j   y j , s = winter s u p p l y · ( L j c o n s u m e r · P c u l j , s = winter D H C ) ,
C p i p e , D H C , s u m m e r = Σ j   y j , s = summer s u p p l y · ( L j c o n s u m e r · P c u l j , s = summer D H C ) ,
C p i p e , D H C , s p r i n g = Σ j   y j , s = spring s u p p l y · ( 1 y j c o m b i n e · y j , s s u p p l y ) · ( L j c o n s u m e r · P c u l j , s = spring D H C ) ,
where P c u l j , s D H C is the cost of the pipeline per unit length from the plant to consumer j in period s.
The calculation methods for the cost of pipeline in the DHC system and in the HRL are different. Best et al. [34] provided a range of flow rates of water corresponding to a series of inner pipeline diameters in a DHC system. The inner diameter is the unique variable of the cost of pipeline per unit length ( P c u l j , s D H C ), while the relationship between P c u l j , s D H C and E j , s s u p p l y can be obtained indirectly according to the data of Best et al. [34]. It should be explained that, because the supply and return temperatures of cold water in a DC system are fixed, the mass flow rate of cold water supplied to a certain consumer is constant. The maximum flow rate in the range can be adopted to determine the relationship between P c u l j , s D H C and E j , s s u p p l y . However, the supply temperature of hot water in a DH system (namely, T s w , H E N , o u t in this paper) is a variable to be optimized that ranges from 70 °C to 100 °C. Because the return temperature was set to 40 °C, the difference between the supply and return temperatures of hot water ranges from 30 °C to 60 °C. It is indicated that the maximum flow rate of hot water can reach twice the minimum value in conditions of a certain heating capacity. Therefore, the average flow rate of the range should be adopted to determine the relationship between P c u l j , s D H C and E j , s s u p p l y for the DH system. Since the data provided by Best et al. [34] are in a large range, the fitting range should be narrowed when E j , s s u p p l y is known in order to reduce the error of the relationship as much as possible. In Section 5, an equation of the relationship between P c u l j , s D H C and E j , s s u p p l y is provided in a case study.
In the second step, the optimal scheme for the DES system in each period can be obtained. The total capacity of energy supplied ( E s t o t a l ) is expressed by Equation (62).
E s t o t a l = Σ j   y j , s s u p p l y · E j , s s u p p l y .
The total income of the plant in this step is calculated using Equation (63).
I D H C = u D C · ( Σ j y j , s = spring s u p p l y · E j , s = spring n e e d ) · t s = spring p e r i o d + u D C · ( Σ j y j , s = summer s u p p l y · E j , s = summer n e e d ) · t s = summer p e r i o d + u D C · ( Σ j y j , s = autumn s u p p l y · E j , s = autumn n e e d ) · t s = autumn p e r i o d + u D H · ( Σ j y j , s = winter s u p p l y · E j , s = winter n e e d ) · t s = winter p e r i o d .

3.3. The Third Step

The set of periods in the third step is as follows: S = {spring, summer, autumn, winter}. The constraint of the energy supply task of plant in each period obtained in the second step is added into the model, which is shown by Equation (64); then, the optimization of the operation scheme for the WHR system can be carried out. The rest of the model is the same as the first step.
E s t o t a l = C O P s · Q s r e c o v e r   s =   spring , summer , autumn , Q s r e c o v e r   s =   winter .

4. Solution Strategy

The MINLP model in this paper is very complex, and many nonlinear terms are involved in the model; therefore, it is very difficult to obtain a good solution. In order to overcome such a problem, two methods are adopted in this paper.
First, the linearization of the cost of pipeline in the HRL is employed. It can be found that the inner diameter of pipeline Din is the unique variable for Pculloop. Hence, the linear fitting of Pculloop can be performed taking Din as the variable, and the nonlinear Equations (38)–(40) can be replaced by a linear equation to reduce the complexity of calculation.
Second, because the cooling capacity in spring and autumn is less than that in summer, the solution process of the third step can be divided into two sub-steps. Only summer and winter are involved in the first sub-step. In the second sub-step, according to the results obtained in the first sub-step, the upper bound of the heat exchanger areas and the upper bound of the pump rated power are given. The optimization of schemes in spring and autumn can be conducted individually. In the objective function of the second sub-step, the capital costs of heat exchangers and pump are not involved.
The solution of the MINLP model can be made much easier by these two methods.

5. Case Study

5.1. Data

In this case, there are 10 hot streams in the plant; the data are shown in Table 1.
There are eight consumers in the region; the data are shown in Table 2.
The time of each period and the physical properties of water in each period are presented in Table 3.
Other physical parameters and economic parameters are listed in Table 4 and Table 5.
According to the demand information of consumers in this case, Equation (65) can be obtained by fitting the data provided by Best et al. [34].
P c u l j , s D H C =   0.4213   ·   E j , s s u p p l y 1000 2 + 8.752   ·   E j , s s u p p l y 1000 + 164.7   s =   spring , summer , autumn , 0.1356   ·   E j , s s u p p l y 1000 2 + 3.519   ·   E j , s s u p p l y 1000   + 156.4     s = winter .
According to the first solution strategy, the approximate linear relationship between Pculloop and Din can be obtained, as shown in Equation (66). Equations (37)–(40) can be replaced by Equation (66).
Pculloop = 1410.8 · Din − 86.31.

5.2. Results

In this study, the MINLP model of each step was solved using the DICOPT solver of GAMS software [36] on a 3.00 GHz personal computer (PC) with 16 GB random-access memory (RAM). The computational time of the three steps was 11.2 s, 0.4 s, and 4.6 s, respectively. The results of the first step are shown in Figure 4 and Figure 5 and Table 6. In the HEN, the waste heat of all hot streams is recovered. There are 13 heat exchangers; 11 exchangers (blue) operate for the maximum cooling potential, and 12 exchangers (red) operate for the maximum heating potential. The district heating potential of the plant is QDH,max = 27,016 kW, while the district cooling potential is QDC,max = 8641 kW.
The scheme obtained in the second step for the DES system is shown in Table 7, in which “-” means that the consumer has no cooling or heating demand, while “0” represents that the consumer has cooling or heating demand but it is not supplied by the plant. In this scheme, the actual total heating capacity of the plant in winter is 26,229.1 kW, and the actual total cooling capacity in summer is 8552.3 kW. The actual capacity of energy supply reaches more than 95% of the potential. In spring and autumn, the total demand of cooling in this region is relatively small and can be fully supplied by the plant (5264.1 kW).
The operation scheme for the WHR in the year is provided by the third step, which is shown in Figure 6 and Figure 7 and Table 8. There are 11 heat exchangers; six exchangers (green) operate in spring and autumn, nine exchangers (blue) operate in summer, and 10 exchangers (red) operate in winter.
Economic data of the second and third steps of this case are shown in Table 9.

5.3. Discussion

In this case, it can be calculated that, when the low-grade waste heat is not recovered, the annual cost of cold utility is 4.204 × 105 USD. By integrating waste heat with the DHC system, this number drops to 1.854 × 105 USD with a decrease of 55.9%. Especially in winter when the capacity of waste heat utilization is the largest, the cost of cold utility drops from 1.401 × 105 USD to 8.985 × 103 USD, with a decrease of 93.6%. The annual income of the plant by district heating and cooling is up to 9.953 × 106 USD. The above two parts of revenue indicate that the annual benefit of the integration amounts to 1.019 × 107 USD, while the total annual cost is 4.366 × 106 USD; therefore, the annual profit of the plant is 5.822 × 106 USD. It is proven that there is a strong motivation to integrate low-grade waste heat into the DHC system.
The district cooling potential of the plant QDC,max is only about one-third of the district heating potential QDH,max. There are two reasons for this difference. First, the hot water used for absorption refrigeration requires a higher temperature (>100 °C); therefore, less heat can be recovered. In this case, the temperature of hot streams H9 and H10 cannot meet the requirement of refrigeration; their heat is not recovered in spring, summer, and autumn. Second, the COP of a single-effect LiBr/H2O absorption chiller is lower than 1.0; therefore, the cold capacity is less than the heat recovered. Both reasons make QDC,max lower than QDH,max. If it is desired to obtain a large cooling capacity, the temperature of hot water entering the chiller T s w , c h i l l e r , i n (namely, the outlet temperature of hot water in HEN T s w , H E N , o u t ) should be in a moderate range. It can be seen from Equation (5) that, when T s w , c h i l l e r , i n increases, COPs increases rapidly at first, and then increases slowly. Although the efficiency of the chiller increases simultaneously, the corresponding temperature of hot water leaving the chiller T s w , c h i l l e r , o u t (i.e., the inlet temperature of hot water in HEN T s w , H E N , i n ) rises correspondingly, leading to a reduction in recovered heat. Therefore, there is an optimal temperature of hot water entering the chiller to ensure the largest cooling capacity.
It can be calculated that the energy loss in the DHC system is about 9%. By considering energy loss and piping cost, consumers with low energy demand and long distance from the plant, such as consumer N8 in winter, are not recommended to be supplied by the plant. A consumer where the cooling demand in spring/autumn is close to the demand in summer, such as consumers N5 and N6 in this case, is suggested to be preferentially supplied by the plant. This can be explained as follows: the DC pipeline can be shared in three periods, whereby more benefit can be acquired under the same investment.
In this case, there is a wide range of temperature of heat streams. If there is only one stage in HEN, the cascade recovery of waste heat cannot be realized. When the number of stages grows, the hot water flows in multistage series and the pressure drop increases, whereby a greater cost of pump is needed. Combining the above two factors, it is reasonable that two stages exist in the HEN. Two heat exchangers exist in hot stream H3, since it has the highest temperature and cascade use of heat is needed. In winter, the heat of hot streams H9 and H10 is used to increase the capacity of heat recovery in the low-grade stage. At this time, the heat exchanger E4 in the low-grade stage for heat stream H3 is shut down, while the heat of H3 is fully recovered by heat exchanger E3 in high-grade stage. On the other hand, in spring, summer, and autumn, the temperatures of hot streams H9 and H10 are too low to meet the requirement of absorption refrigeration; thus, the heat exchangers E10 and E11 are shut down. At this time, the heat exchanger E4 is turned on so that the total amount of heat recovery in the low-grade stage increases. The above operation keeps the ratio of the heat recovery capacity in the two stages basically unchanged to make the total area smaller.
In the WHR system, the outlet temperature of hot water T s w , H E N , o u t is an important variable to be optimized, which has a direct impact on the results. As mentioned, a medium T s w , H E N , o u t should be adopted in summer when the cooling capacity is large. On the other hand, in spring and autumn, the cooling capacity is relatively small. Because the waste heat resource is sufficient in the periods to achieve the required cooling capacity, a lower T s w , H E N , o u t can be adopted. Although the COP of the chiller is relatively low, waste heat can be recovered as much as possible. By doing this, the cost of cold utility can be reduced. In spring and autumn, the difference between the inlet and outlet temperature of hot water in HEN is relatively small, and the task of waste heat recovery is completed in one stage to reduce pressure drop. When the system is operated in winter, a tradeoff between Cexchanger and distance-related cost (including Cpipe,loop and Cpump) under the same capacity of heat recovery should be made. When Tw,HEN,in is fixed at 40 °C, the heat transfer temperature difference decreases with the rise in Tw,HEN,out, and the total heat transfer area needed for exchangers increases, so that the Cexchanger increases. Meanwhile, Cpipe,loop and Cpump can be reduced due to the reduction in mw,total. In this study, the scheme of the WHR system in winter should cooperate with that in the other three periods, and the temperature of hot water gained in winter is moderate (78.9 °C). It can be found that the total mass flow rate of hot water in each season is almost equal, indicating that the velocity of hot water is appropriate in each period.
In the third step, if the operation scheme of WHR system is not optimized, the waste heat of hot streams is recovered by hot water in a purely parallel structure. In spring, summer, and autumn, the inlet and outlet temperatures of the hot water are the same as that in the first step. In winter, the outlet temperature of hot water is fixed at 85 °C as in [20]. It can be calculated that the TAC’ would change to 381,546 USD, which is greater than the optimal solution. Therefore, it is of significance to optimize the multiperiod operation scheme of the WHR system. According to the previous discussion, T s = s u m m e r w , H E N , o u t can be regarded as the most impactful variable. A sensitivity analysis of T s = s u m m e r w , H E N , o u t and TAC’ was performed after the third step. The result is shown in Figure 8. The range of T s = s u m m e r w , H E N , o u t is relatively narrow. In this case, if the value of T s = s u m m e r w , H E N , o u t exceeds the range in Figure 8, the district cooling task in summer cannot be accomplished. On the left and right sides of 124 °C in the figure, a higher cost of pipeline and a higher cost of cold utility respectively led to a worse solution.

6. Conclusions

A systematic framework method for utilizing industrial low-grade heat in DHC with a multiperiod condition was proposed. Both design and operation variables in different periods could be determined. A case study was carried out, and the results showed that the heat exchanger network in the WHR system is not a simple parallel structure. A two-stage series/parallel network of intermediate fluids should be established. In different periods, the operation plan of heat exchangers should be flexible, while the total flow rate of hot water should be basically the same. The temperatures of hot water obtained in WHR ( T s w , H E N , o u t ) should be different in periods when the cooling capacities are significantly different. In seasons with a larger cooling capacity, T s w , H E N , o u t is suggested to be close to 121 °C, while it can take a smaller value in seasons with a smaller cooling capacity. In this case, the total annual cost of WHR system (TAC’) will increase by more than 6000 USD without the above multiperiod operation optimization. The cost saving proves the worth of the third step. In a DES system, priority should be given to the consumers where (1) the cooling demand in spring/autumn is close to that in summer, and (2) the distance from the plant is relatively short while the heating demand is large. The economic data of the case study show that the total annual benefit of the integration is up to 1.019 × 107 USD, which is greater than the total annual cost 4.366 × 106 USD. Strong motivation for an integration of low-grade industrial waste heat and a DHC system is indicated.

Author Contributions

Conceptualization, X.F. and Y.W.; methodology, Y.W.; software, L.Z.; validation, L.Z. and Y.W.; formal analysis, Y.W.; investigation, L.Z.; resources, X.F.; data curation, L.Z.; writing—original draft preparation, L.Z.; writing—review and editing, Y.W.; visualization, L.Z.; supervision, Y.W.; project administration, Y.W.; funding acquisition, Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

Financial support from the National Natural Science Foundation of China under Grant No. 22022816 and 22078358 is gratefully acknowledged.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Abbreviations
DESDistrict energy supply
DCDistrict cooling
DHDistrict heating
DHCDistrict heating and cooling
HENHeat exchanger network
HIHeat Integration
HRLHeat recovery loop
MPMathematical programming
MINLPMixed integer nonlinear programming
WHRWaste heat recovery
Subscripts
iHot stream index
jConsumer index
sPeriod index
Parameters and Variables
A i , k , s p e r i o d Area of heat exchanger in period (m2)
A i , k d e s i g n Design area of heat exchanger (m2)
AfAnnual factor
CcuCost of cold utility (USD)
CexchangerCapital cost of heat exchangers (USD)
Cpipe,loopCapital cost of pipeline in HRL (USD)
Cpipe,DHCCapital cost of pipeline in DHC system (USD)
CpumpTotal cost of pump (USD)
Cpump,capCapital cost of pump (USD)
Cpump,opOperation cost of pump (USD)
CstationCapital cost of the station (USD)
COPsCoefficient of performance of the absorption chiller
cpwSpecific heat capacity of water (kJ·kg−1·°C−1)
dti,k,sTemperature approach of heat exchanger (°C)
E j , s n e e d Energy demand of consumer (kW)
E j , s s u p p l y Actual energy capacity should be supplied to consumer (kW)
E s t o t a l Total capacity of energy supplied (kW)
F i h s Flow capacity of hot stream (kW·°C−1)
hwFilm heat transfer coefficient of hot water (kW·m−2·°C−1)
h i h s Film heat transfer coefficient of hot stream (kW·m−2·°C−1)
IDHCIncome of the plant by district heating and cooling (USD)
LloopDistance between station and production area (m)
L j c o n s u m e r Distance between the plant and consumer (m)
LMTDi,k,sLogarithmic mean temperature difference of heat exchanger (°C)
m s w , t o t a l Total mass flow rate of hot water (kg·s−1)
m i , k , s w , s p l i t Mass flow rate of hot water in heat exchangers (kg·s−1)
P s p e r i o d Power of pump in each period (kW)
PrateRated power of pump (kW)
PculloopCapital cost of the pipeline per unit length in HRL (USD·m−1)
P c u l j , s D H C Capital cost of the pipeline per unit length in DHC system (USD·m−1)
qi,k,sHeat load of heat exchanger (kW)
q i , s c u Heat load of cooler (kW)
q i h s Total heat load of hot stream (kW)
Q s c o l d Cooling output of chiller (kW)
Q s r e c o v e r Capacity of the waste heat recovered (kW)
QDC,maxDistrict cooling potential of the plant (kW)
QDH,maxDistrict heating potential of the plant (kW)
t s p e r i o d Time of each period (h)
ttotalTime of total year (h)
T i h s , i n Initial temperature of hot stream (°C)
T i h s , o u t Target temperature of hot stream (°C)
T i , k , s h s Temperature node of hot stream in HEN (°C)
T s w , c h i l l e r , i n Temperature of hot water entering the chiller (°C)
T s w , c h i l l e r , o u t Temperature of hot water leaving the chiller (°C)
T s w , H E N , i n Inlet temperature of hot water in HEN (°C)
T s w , H E N , o u t Outlet temperature of hot water in HEN (°C)
T k , s w , H E N Temperature node of hot water in HEN (°C)
TACTotal annual cost in the first step (USD)
TAC’Total annual cost in the third step (USD)
TAPTotal annual profit of the DES system (USD)
ucuUnit price of cold utility (USD·kW−1·year−1)
uDCUnit price of district cooling (USD·MWh−1)
uDHUnit price of district heating (USD·MWh −1)
ueUnit price of electricity (USD·kWh −1)
y j c o m b i n e Binary parameter for pipeline combination
y j , s s u p p l y Binary variable for energy supply selection
zi,kBinary variable for heat exchanger
Δ p s d i s t a n c e Pressure drop of hot water in long pipeline (Pa)
Δ p s H E N Pressure drop of hot water in HEN (Pa)
Δ p k , s s t a g e Pressure drop of hot water in each stage (Pa)
Δ p s t o t a l Total pressure drop of hot water (Pa)
Δ p i , k , s t u b e Pressure drop of hot water in tube side of heat exchanger (Pa)
ΔTminMinimum temperature difference (°C)
ηEfficiency of pump
μsViscosity of hot water (mPa·s−1)
ρsDensity of hot water (kg·m−3)
ΓUpper bound of temperature difference (°C)

References

  1. Ji, C.; Qin, Z.; Dubey, S.; Choo, F.H.; Duan, F. Three-dimensional transient numerical study on latent heat thermal storage for waste heat recovery from a low temperature gas flow. Appl. Energy 2017, 205, 1–12. [Google Scholar] [CrossRef]
  2. Klemeš, J.J.; Kravanja, Z. Forty years of Heat Integration: Pinch Analysis (PA) and Mathematical Programming (MP). Curr. Opin. Chem. Eng. 2013, 2, 461–474. [Google Scholar] [CrossRef]
  3. Linnhoff, B.; Hindmarsh, E. The pinch design method for heat exchanger networks. Chem. Eng. Sci. 1983, 38, 745–763. [Google Scholar] [CrossRef]
  4. Floudas, C.A.; Ciric, A.R.; Grossmann, I.E. Automatic synthesis of optimum heat exchanger network configurations. AIChE J. 1986, 32, 276–290. [Google Scholar] [CrossRef]
  5. Yee, T.F.; Grossmann, I.E. Simultaneous optimization models for heat integration II: Heat exchanger network synthesis. Comput. Chem. Eng. 1990, 14, 1165–1184. [Google Scholar] [CrossRef]
  6. Shang, J.; Ji, Z.; Qiu, M.; Ma, L. Multi-objective optimization of high-sulfur natural gas purification plant. Pet. Sci. 2019, 16, 1430–1441. [Google Scholar] [CrossRef] [Green Version]
  7. Forman, C.; Muritala, I.K.; Pardemann, R.; Meyer, B. Estimating the global waste heat potential. Renew. Sustain. Energy Rev. 2016, 57, 1568–1579. [Google Scholar] [CrossRef]
  8. Van de Bor, D.M.; Infante Ferreira, C.A.; Kiss, A.A. Low grade waste heat recovery using heat pumps and power cycles. Energy 2015, 89, 864–873. [Google Scholar] [CrossRef]
  9. Fang, H.; Xia, J.; Zhu, K.; Su, Y.; Jiang, Y. Industrial waste heat utilization for low temperature district heating. Energy Policy 2013, 62, 236–246. [Google Scholar] [CrossRef]
  10. Huang, F.; Zheng, J.; Baleynaud, J.M.; Lu, J. Heat recovery potentials and technologies in industrial zones. J. Energy Inst. 2017, 90, 951–961. [Google Scholar] [CrossRef]
  11. Gao, X.; Yin, X.; Yang, S.; Yang, D. Improved organic rankine cycle system coupled with mechanical vapor recompression distillation for separation of benzene-toluene mixture. Process Integr. Optim. Sustain. 2019, 3, 189–198. [Google Scholar] [CrossRef]
  12. Xia, L.; Liu, R.; Zeng, Y.; Zhou, P.; Liu, J.; Cao, X.; Xia, S. A review of low-temperature heat recovery technologies for industry processes. Chin. J. Chem. Eng. 2019, 27, 2227–2237. [Google Scholar] [CrossRef]
  13. Garcia, S.I.; Garcia, R.F.; Carril, J.C.; Garcia, D.I. A review of thermodynamic cycles used in low temperature recovery systems over the last two years. Renew. Sustain. Energy Rev. 2017, 81, 760–767. [Google Scholar] [CrossRef]
  14. Lund, H.; Werner, S.; Wiltshire, R.; Svendsen, S.; Thorsen, J.E.; Hvelplund, F.; Mathiesen, B.V. 4th generation district heating (4gdh): Integrating smart thermal grids into future sustainable energy systems. Energy 2014, 68, 1–11. [Google Scholar] [CrossRef]
  15. Li, Y.; Rezgui, Y.; Zhu, H. District heating and cooling optimization and enhancement—Towards integration of renewables, storage and smart grid. Renew. Sustain. Energy Rev. 2017, 72, 281–294. [Google Scholar] [CrossRef]
  16. Tveit, T.M.; Aaltola, J.; Laukkanen, T.; Laihanen, M.; Fogelholm, C.J. A framework for local and regional energy system integration between industry and municipalities—Case study UPM-Kymmene Kaukas. Energy 2006, 31, 2162–2175. [Google Scholar] [CrossRef]
  17. Kapil, A.; Bulatov, I.; Smith, R.; Kim, J.K. Process integration of low grade heat in process industry with district heating networks. Energy 2012, 44, 11–19. [Google Scholar] [CrossRef]
  18. Fang, H.; Xia, J.; Jiang, Y. Key issues and solutions in a district heating system using low-grade industrial waste heat. Energy 2015, 86, 589–602. [Google Scholar] [CrossRef]
  19. Wang, J.; Wang, Z.; Zhou, D.; Sun, K. Key issues and novel optimization approaches of industrial waste heat recovery in district heating systems. Energy 2019, 188, 116005. [Google Scholar] [CrossRef]
  20. Fitó, J.; Hodencq, S.; Ramousse, J.; Wurtz, F.; Stutz, B.; Debray, F.; Vincent, B. Energy- and exergy-based optimal designs of a low-temperature industrial waste heat recovery system in district heating. Energy Convers. Manag. 2020, 211, 112753. [Google Scholar] [CrossRef]
  21. Khosravi, A.; Laukkanen, T.; Vuorinen, V.; Syri, S. Waste heat recovery from a data centre and 5G smart poles for low-temperature district heating network. Energy 2021, 218, 119468. [Google Scholar] [CrossRef]
  22. Srikhirin, P.; Aphornratana, S.; Chungpaibulpatana, S. A review of absorption refrigeration technologies. Renew. Sustain. Energy Rev. 2011, 5, 343–372. [Google Scholar] [CrossRef]
  23. Liew, P.Y.; Walmsley, T.G.; Wan Alwi, S.R.; Manan, Z.A.; Klemeš, J.J.; Varbanov, P.S. Integrating district cooling systems in Locally Integrated Energy Sectors through Total Site Heat Integration. Appl. Energy 2016, 184, 1350–1363. [Google Scholar] [CrossRef]
  24. Zhang, L.; Wang, Y.; Feng, X. Design and Operation Optimization of Industrial Waste Heat Recovery for District Heating and Cooling. In Proceedings of the 23rd Conference on Process Integration, Modelling and Optimisation for Energy Saving and Pollution Reduction (PRES’20), Xi’an, China, 17–21 August 2020; p. 265. [Google Scholar]
  25. Söderman, J. Optimisation of structure and operation of district cooling networks in urban regions. Appl. Therm. Eng. 2007, 27, 2665–2676. [Google Scholar] [CrossRef]
  26. Chang, C.; Chen, X.; Wang, Y.; Feng, X. An efficient optimization algorithm for waste Heat Integration using a heat recovery loop between two plants. Appl. Therm. Eng. 2016, 105, 799–806. [Google Scholar] [CrossRef]
  27. Arabkoohsar, A.; Andresen, G.B. Supporting district heating and cooling networks with a bifunctional solar assisted absorption chiller. Energy Convers. Manag. 2017, 148, 184–196. [Google Scholar] [CrossRef]
  28. Wang, Y.; Wang, C.; Feng, X. Optimal match between heat source and absorption refrigeration. Comput. Chem. Eng. 2017, 102, 268–277. [Google Scholar] [CrossRef]
  29. Chen, J.J.J. Comments on improvements on a replacement for the logarithmic mean. Chem. Eng. Sci. 1987, 42, 2488–2489. [Google Scholar] [CrossRef]
  30. Verheyen, W.; Zhang, N. Design of flexible heat exchanger network for multi-period operation. Chem. Eng. Sci. 2006, 61, 7730–7753. [Google Scholar]
  31. Peters, M.S.; Timmerhaus, K.D.; West, R.E. Plant Design and Economics for Chemical Engineers, 5th ed.; McGraw-Hill: Boston, MA, USA, 2003. [Google Scholar]
  32. Stijepovic, M.Z.; Linke, P. Optimal waste heat recovery and reuse in industrial zones. Energy 2011, 36, 4019–4031. [Google Scholar]
  33. Soltani, H.; Shafiei, S. Heat exchanger networks retrofit with considering pressure drop by coupling genetic algorithm with LP (linear programming) and ILP (integer linear programming) methods. Energy 2011, 36, 2381–2391. [Google Scholar] [CrossRef]
  34. Best, R.E.; Rezazadeh Kalehbasti, P.; Lepech, M.D. A novel approach to district heating and cooling network design based on life cycle cost optimization. Energy 2020, 194, 1–16. [Google Scholar] [CrossRef]
  35. Jabbari, B.; Tahouni, N.; Ataei, A.; Panjeshahi, M.H. Design and optimization of CCHP system incorporated into kraft process, using Pinch Analysis with pressure drop consideration. App. Therm. Eng. 2013, 61, 88–97. [Google Scholar] [CrossRef]
  36. Rosenthal, R.E. GAMS—A User’s Guide; GAMS Development Corporation: Washington, DC, USA, 2008. [Google Scholar]
Figure 1. The indirect integration between waste heat and the district heating and cooling (DHC) system.
Figure 1. The indirect integration between waste heat and the district heating and cooling (DHC) system.
Energies 14 02190 g001
Figure 2. The relationship between T s w , c h i l l e r , i n and T s w , c h i l l e r , o u t [28].
Figure 2. The relationship between T s w , c h i l l e r , i n and T s w , c h i l l e r , o u t [28].
Energies 14 02190 g002
Figure 3. The stage-wise superstructure.
Figure 3. The stage-wise superstructure.
Energies 14 02190 g003
Figure 4. The heat exchanger network (HEN) of the first step.
Figure 4. The heat exchanger network (HEN) of the first step.
Energies 14 02190 g004
Figure 5. Heat load of heat exchangers in the first step.
Figure 5. Heat load of heat exchangers in the first step.
Energies 14 02190 g005
Figure 6. The HEN of the third step.
Figure 6. The HEN of the third step.
Energies 14 02190 g006
Figure 7. Heat load of heat exchangers in the third step.
Figure 7. Heat load of heat exchangers in the third step.
Energies 14 02190 g007
Figure 8. Sensitivity analysis of T s = s u m m e r w , H E N , o u t .
Figure 8. Sensitivity analysis of T s = s u m m e r w , H E N , o u t .
Energies 14 02190 g008
Table 1. The data of hot streams.
Table 1. The data of hot streams.
i T i h s , i n (°C) T i h s , o u t (°C) F i h s (kW/°C) q i h s (kW) h i h s (kW/m2·°K)
H1140752415602
H2186683440122
H32001053836102
H4144403031202
H5175933024602
H6160424047202
H7153563231042
H8130503528002
H9100805410802
H1095352615602
Table 2. The data of consumers.
Table 2. The data of consumers.
j L j u s e r (m) E j , s n e e d (kW) y j c o m b i n e
s = Springs = Summers = Autumns = Winter
N1520000040000
N2630010002500100020000
N3700001800030000
N4790000070000
N5950020004000200050001
N610,40018002000180060001
N711,00000080000
N811,80003500050000
Table 3. The period-related data.
Table 3. The period-related data.
s t s p e r i o d (h)ρs (kg/m3)μs (mPa·s)
spring15009450.242
summer28809450.242
autumn15009450.242
winter28809800.430
Table 4. Physical parameters of the case.
Table 4. Physical parameters of the case.
ItemsRemarksData
ΔTminMinimum approach temperature difference10 °C
hwFilm heat transfer coefficient of hot water1.5 kW/m2 °C
cpwSpecific heat capacity of hot water4.2 kJ/kg °C
ΦViscosity correction factor1.05
kConductivity0.6 W/m °C
diTube internal diameter of heat exchanger15.4 mm
deTube external diameter of heat exchanger19.1 mm
LloopDistance between station and production area400 m
ηPump efficiency0.7
Table 5. Economic parameters of the case.
Table 5. Economic parameters of the case.
ItemsRemarksData
AfAnnualized factor [26]0.264
αCost parameters of heat exchanger [32]11,000 $
β150 $/m2
γ1
b1Cost parameters of cooling station [25]400,000 $
b2400 $/kW
a1Cost parameters of pipe [32]0.82 $/kg
a2185 $/m1.48
a36.8 $/m
a4295 $/m2
c1Cost parameters of pump [35]8600 $
c27310 $/W0.2
c30.2
ucuUnit price of cold utility [26]15 $/kW·y
ueUnit price of electricity [26]0.1 $/kW·h
uDCUnit price of cooling [15]60 $/MW·h
uDHUnit price of heating [15]100 $/MW·h
Table 6. Operation parameters of the first step.
Table 6. Operation parameters of the first step.
s m s w , t o t a l
(kg/s)
T s w , H E N , i n
(°C)
T s w , H E N , o u t
(°C)
P s p e r i o d
(kW)
Q s r e c o v e r
(kW)
Q s D H
(kW)
Q s D C
(kW)
COPs
summer177.440.075.74.20427,01627,016--
winter180.2103.9121.05.34812,744-86410.678
Table 7. The result data of the second step.
Table 7. The result data of the second step.
j y j , s s u p p l y · E j , s s u p p l y (kW)
s = Springs = Summers = Autumns = Winter
N1---4214.6
N21065.401065.40
N3-1931.2-0
N4---7578.4
N52200.44400.72200.45500.9
N61998.32220.41998.30
N7---8935.2
N8-0-0
Table 8. Operation parameters of the third step.
Table 8. Operation parameters of the third step.
s m s w , t o t a l
(kg/s)
T s w , H E N , i n
(°C)
T s w , H E N , o u t
(°C)
P s p e r i o d
(kW)
Q s r e c o v e r
(kW)
COPs
spring and autumn156.6101.5114.33.6598420.00.625
summer156.3105.2124.04.39112,341.00.693
winter160.440.078.95.09726,229.10
Table 9. The economic data.
Table 9. The economic data.
IDHC (USD)Cpipe,DHC (USD)Cstation (USD)TAP (USD)
9,953,2803,168,5911,007,9495,776,740
Ccu (USD)Cexchanger (USD)Cpipe,loop (USD)Cpump (USD)TAC’ (USD)
185,43975,44796,06118,312375,259
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, L.; Wang, Y.; Feng, X. A Framework for Design and Operation Optimization for Utilizing Low-Grade Industrial Waste Heat in District Heating and Cooling. Energies 2021, 14, 2190. https://doi.org/10.3390/en14082190

AMA Style

Zhang L, Wang Y, Feng X. A Framework for Design and Operation Optimization for Utilizing Low-Grade Industrial Waste Heat in District Heating and Cooling. Energies. 2021; 14(8):2190. https://doi.org/10.3390/en14082190

Chicago/Turabian Style

Zhang, Lingwei, Yufei Wang, and Xiao Feng. 2021. "A Framework for Design and Operation Optimization for Utilizing Low-Grade Industrial Waste Heat in District Heating and Cooling" Energies 14, no. 8: 2190. https://doi.org/10.3390/en14082190

APA Style

Zhang, L., Wang, Y., & Feng, X. (2021). A Framework for Design and Operation Optimization for Utilizing Low-Grade Industrial Waste Heat in District Heating and Cooling. Energies, 14(8), 2190. https://doi.org/10.3390/en14082190

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

Article Metrics

Back to TopTop