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

: 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 signiﬁcance 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 ﬂuctuations 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 signiﬁcant beneﬁts to the system. By solving the mixed integer nonlinear programming model, the optimal operation plans of the integration in different periods can be obtained. study show that the total annual beneﬁt of the integration is up to 1.019 × 10 7 USD, which is greater than the total annual cost 4.366 × 10 6 USD. Strong motivation for an integration of low-grade industrial waste heat and a DHC system is indicated.


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 m 3 /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.

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 = C cu + C exchanger + C station + C pipe,loop + C pump − I DHC }, (1) where TAC is the total annual cost in the first step, C cu is the cost of cold utility that is still needed, C exchanger is the capital cost of heat exchangers installed to recover waste heat, C station is the capital cost of the station, including the cost of the absorption chiller, C pipe,loop and C pump are the capital cost of pipeline and the cost of pump in the HRL, And I DHC 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 C station ). 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.
where TAP is the total annual profit of the DES system, and C pipe,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 C station and I DHC that were determined in the second step were deleted in the objective function, which is expressed by Equation (3).
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.

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.

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 H 2 O/NH 3 and LiBr/H 2 O [22]. Compared with LiBr/H 2 O, H 2 O/NH 3 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-H 2 O 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/H 2 O refrigerator is calculated using Equation (4) [28].
where Q cold s is the cooling output of the chiller, and Q recover s 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/H 2 O 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 COP s and the temperature of hot water entering the chiller T w,chiller,in s [28]; accordingly, piecewise linearization should be used to describe this relationship, as expressed as Equation (5).
In the simulation conditions of Wang et al. [28], there is a one-to-one correspondence between T w,chiller,in s and the temperature of hot water leaving the chiller T w,chiller,out s . According to their data, the relationship between T w,chiller,in s and T w,chiller,out s is shown in Figure 2. It can be seen that there is an obvious linear relationship between T w,chiller,in s and T w,chiller,out s . Equation (13) can be obtained by linear fitting.
ya s + yb s + yc s + yd s = 1, where ya s , yb s , yc s , and yd s are binary variables (0 or 1), while T w The investment of the station including the capital cost of the absorption chiller is calculated using Equation (14).
where Af is the annual factor, and b 1 and b 2 are the parameters of the cost of the small station.

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 hs,in i , T hs,out I, and In current DH systems, the temperature of atmospheric hot water sent to consumers (T w,HEN,out s ) is generally required to be higher than 70 • C [14]. The temperature of water sent back to the station (T w,HEN,in s ) 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).
The total capacity of waste heat recovered is calculated by Equation (21). cp w represents the specific heat capacity of water.
Energy and mass balance constraints for hot water and hot streams are expressed by Equations (22)-(27).
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).
where z i,k is the binary variable to indicate whether the heat exchanger between hot water and hot stream i at stage k exists. q hs i is the total heat load of hot stream i. dt i,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 z i,k is equal to 1; if the heat exchanger does not exist, z i,k can be 0 or 1 in Equations (28)- (30). In order to minimize the cost of heat exchangers, z i,k is forced to be equal to 0 by objective function.
The temperature approaches must be greater than the minimum temperature difference ∆T min in order to meet the requirement for heat transfer. This relationship is shown by Equation (31).
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 period i,k,s ).
where LMTD i,k,s is the logarithmic mean temperature difference of this heat exchanger. h w and h hs i are the film heat transfer coefficients of hot water and hot stream i.
In order to reduce the solution difficulty, LMTD i,k,s is calculated using the Chen approximation [29], which is expressed by Equation (33).
In multiperiod HENs, the design area of exchanger (A design i,k ) is the maximum area needed in all periods [30], which is expressed by Equation (34).
Equations (35) and (36) are used to calculate the cost of cold utility and the capital cost of heat exchangers, respectively.
where u cu is the unit cost of cold utility. t period s and t total are the time of period s and total year. α, β, and γ are the cost parameters of the heat exchanger.

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 (D in ) of the pipeline between production area and station is expressed by Equation (37) according to the method of Peters et al. [31].
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 (Pcul loop ) was calculated using the method of Stijepovic et al. [32], which is shown by Equations (38)-(40).
The total cost of the pipeline in the HRL is calculated using Equation (41). L loop is the distance between station and production area.

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 total s is expressed in Equation (42).
where ∆p distance s and ∆p HEN s are the pressure drop of hot water in the long pipeline and HEN in period s. ∆p distance s 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 tube i,k,s ) is calculated using the method of Soltani et al. [33], which is shown by Equations (43) 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 stage k,s ); this relationship is shown in Equation (45).
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).
The power of pump in each period (P period s ) and the rated power of pump (P rate ) are expressed by Equations (47) The cost of pump C pump is equal to the sum of capital cost and operation cost, which are calculated using Equations (50) and (51), respectively.
where c 1 , c 2 , and c 3 are the parameters of capital cost of the pump, u e is the unit price of electricity, and η is the efficiency of pump.

Calculation of Income
The maximum heating and cooling potentials (Q DH,max and Q DC,max ) are represented by Equations (52) and (53). The total income of the plant in these two periods is calculated by Equation (54). u DC and u DH are the unit prices of heating and cooling respectively.

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 supply j,s ) and the demand of consumer j (E need j,s ) cannot be ignored. The loss of energy per unit distance in DHC system can be taken as 1%/km [17]. Equation (55) 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 combine j , which is shown in Equation (57).
When consumer j is cooled by the plant in spring/autumn, if it is also cooled in summer, and the condition y combine j = 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 C pipe, DHC can be calculated using Equations (58)-(61). C pipe, DHC = C pipe, DHC, wiNTer + C pipe, DHC, summer + C pipe, DHC, spring , C pipe,DHC,winter = Σ j y supply j,s=winter ·(L consumer j ·Pcul DHC j,s=winter ), C pipe,DHC,summer = Σ j y supply j,s=summer ·(L consumer j ·Pcul DHC j,s=summer ), C pipe,DHC,spring = Σ j y supply j,s=spring ·(1 − y combine j ·y supply j,s )·(L consumer j ·Pcul DHC j,s=spring ), where Pcul DHC j,s 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 (Pcul DHC j,s ), while the relationship between Pcul DHC j,s and E supply j,s 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 Pcul DHC j,s and E supply j,s . However, the supply temperature of hot water in a DH system (namely, T w,HEN,out s 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 Pcul DHC j,s and E supply j,s 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 supply j,s is known in order to reduce the error of the relationship as much as possible.
In Section 5, an equation of the relationship between Pcul DHC j,s and E supply j,s 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 total The total income of the plant in this step is calculated using Equation (63).

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. (64)

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 Pcul loop . Hence, the linear fitting of Pcul loop 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.

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 Tables 4 and 5.  Distance between station and production area 400 m η Pump efficiency 0.7 Table 5. Economic parameters of the case.

Items Remarks Data
Af Unit price of electricity [26] 0.1 $/kW·h u DC Unit price of cooling [15] 60 $/MW·h u DH Unit price of heating [15] 100 $/MW·h 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]. (66)

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 Figures 4 and 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 Q DH,max = 27,016 kW, while the district cooling potential is Q DC,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 Figures 6 and 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.

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 × 10 5 USD. By integrating waste heat with the DHC system, this number drops to 1.854 × 10 5 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 × 10 5 USD to 8.985 × 10 3 USD, with a decrease of 93.6%. The annual income of the plant by district heating and cooling is up to 9.953 × 10 6 USD. The above two parts of revenue indicate that the annual benefit of the integration amounts to 1.019 × 10 7 USD, while the total annual cost is 4.366 × 10 6 USD; therefore, the annual profit of the plant is 5.822 × 10 6 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 Q DC,max is only about one-third of the district heating potential Q DH,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/H 2 O absorption chiller is lower than 1.0; therefore, the cold capacity is less than the heat recovered. Both reasons make Q DC,max lower than Q DH,max . If it is desired to obtain a large cooling capacity, the temperature of hot water entering the chiller T w,chiller,in s (namely, the outlet temperature of hot water in HEN T w,HEN,out s ) should be in a moderate range. It can be seen from Equation (5) that, when T w,chiller,in s increases, COP s 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 w,chiller,out s (i.e., the inlet temperature of hot water in HEN T w,HEN,in s ) 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 w,HEN,out s is an important variable to be optimized, which has a direct impact on the results. As mentioned, a medium T w,HEN,out s 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 w,HEN,out s 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 C exchanger and distance-related cost (including C pipe,loop and C pump ) under the same capacity of heat recovery should be made. When T w,HEN,in is fixed at 40 • C, the heat transfer temperature difference decreases with the rise in T w,HEN,out , and the total heat transfer area needed for exchangers increases, so that the C exchanger increases. Meanwhile, C pipe,loop and C pump can be reduced due to the reduction in m w,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 w,HEN,out s=summer can be regarded as the most impactful variable. A sensitivity analysis of T w,HEN,out s=summer and TAC' was performed after the third step. The result is shown in Figure 8. The range of T w,HEN,out s=summer is relatively narrow. In this case, if the value of T w,HEN,out s=summer 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.

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 w,HEN,out s ) should be different in periods when the cooling capacities are significantly different. In seasons with a larger cooling capacity, T w,HEN,out s 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 × 10 7 USD, which is greater than the total annual cost 4.366 × 10 6 USD. Strong motivation for an integration of low-grade industrial waste heat and a DHC system is indicated.