1. Introduction
With the continuous increase of energy consumption in industrial processes, the energy crisis has further intensified on account of the limited fossil energy reserves. Energy prices are also on the rise at the same time, which drives the pursuit of energy-saving technologies and methods. Heat exchanger networks (HENs) are an inevitable part in processing enterprises for heat recovery. Hot process streams that need to be cooled down and cold process streams that need to be heated up widely exist in process enterprises. Synthesis of HENs is able to obtain HEN structures with reasonable matches between these hot and cold process streams and recover waste heat to the greatest extent. Thus the additional utility consumption will be greatly reduced. On the other hand, utility systems based on steam power cycles are one of the major sources which can export multi-level steam and power to industrial enterprises from a single primary energy source. Energy consumption in utility systems accounts for a large proportion of total industrial energy consumption, which makes it meaningful to improve energy efficiency of utility system. These two parts, the HEN and the utility system, are involved the whole process of energy supply, energy consumption and energy saving. While in implementation, the utility system is usually designed independently with its connections with the processes ignored, so the simultaneous synthesis and optimization of HENs and utility systems should be improved.
Determining the use of utilities is always one important content in HEN synthesis, either within a sequential synthesis method or a simultaneous synthesis method. In pinch technology [
1] and the trans-shipment model-based method [
2], the consumption of cold and hot utilities is targeted by monitoring of pinch point locations, before designing the network structure in terms of paring hot and cold streams; in the superstructure-based simultaneous method [
3], utility consumption is optimized by the trade-off with the capital costs of heat exchangers. Although the study of HEN synthesis has been developed for decades, there are still deficiencies. For example, most previous works on HEN synthesis only used one type of utility at the stream end and assumed that the temperatures of these utilities are able to satisfy all the hot and cold demands, while the fact is that the utility can be in many forms, such as flue gas, steam with different temperatures and hot water. Thus, studies have been launched to investigate the reasonable use of these utilities in HEN synthesis. Costa and Queiroz [
4] introduced an extension of the problem table algorithm to optimize and analyze multiple utilities selection and utilization instead of using a grand composite curve. Although the advantages of multi-level utility selection are not analyzed, it is an effective way to optimize multi-level utility utilization with a table algorithm. Salama [
5] developed a simple and direct numerical geometry-based technique to target optimal assignment of multiple utilities, but only the heat load was optimized without considering the equipment investment. Shenoy et al. [
6] presented a cost-optimal targeting methodology considering the trade-off between energy consumption and equipment investment simultaneously. Optimal selections and loads for multiple utilities were determined based on pinch analysis and the cheapest utility principle. In essence it was still a sequential method, which makes it impossible to achieve a real trade-off between energy and equipment costs. Isafiade and Fraser [
7] studied an interval-based MINLP superstructure where the intervals were defined according to the supply and target temperatures of hot or cold process streams. The superstructure model was then applied to HEN synthesis with multi-level available utilities, trading off operating and capital costs for each type of utility. Many nonlinear terms were ignored by mixing split streams at equal temperature, which made the obtained results improvable. Ponce-Ortega et al. [
8] developed a stage-wise superstructure that allowed existence of intermediate placement of multiple utilities within each stage. Load and placement of multiple utilities were treated as optimization variables by employing disjunctive programming formulation, rather than regarding utilities as process streams and setting their capacity flowrates as optimization variables. Na et al. [
9] proposed a modified superstructure that contained utility sub-stages between adjacent stages. Series utility locations were fixed to facilitate convergence. These studies greatly expanded the design space of the HEN structure and utilization of multiple utilities, but the selection of multiple utilities was mainly based on their price, and the interactions with other systems was not investigated in depth. Zhang et al. [
10] explored a new HEN superstructure presentation named the stage-wise chessboard model for management of feasible research regions. A random walk algorithm was employed to lower the calculation load. The problem can be solved faster, but many possible matches are not included in the proposed superstructure because it was equivalent to a stage-wise superstructure with only one stage. Pavão et al. [
11] considered the solving difficulty of complex mathematical models resulting from an enhanced stage-wise superstructure which includes the use of multiple utilities at single stream branches. An enhanced meta-heuristic solution method was presented to handle the complex mathematical model, and the superstructure had no difference from the superstructure in Ponce et al. [
8] mentioned above. Ma et al. [
12] launched multi-objective optimization of interplant HENs operated for multi-periods, which used steam as the heat transfer medium. This research analyzed the conflict between environmental impacts and exchanger investment according to the trade-off between utility consumption and exchanger areas in a case study, but the influence of multi-level steam distribution on the environmental and economic objectives was ignored. Besides, power generation in the utility system was not included in the objective function either, thus the HEN and utility system were not globally integrated.
HENs have been further integrated with utility systems based on their close interactions. Klemeš et al. [
13] studied the simultaneous synthesis of a production process and a utility system based on pinch analysis methodology in a total site integration, but the optimization was carried out after integration within a single plant was completed in advance. Thus the heat recovery and utility system were not optimized simultaneously. Liew et al. [
14] optimized the design and operation of a centralized utility system to adapt to shutdowns or process upsets, making a trade-off between operational adaptability and operating costs. Detailed synthesis of the HEN was not included either. Chen and Lin [
15] proposed a MINLP model to design a steam network and heat recovery network simultaneously. Hot and cold utilities were only placed at the stream ends, while steam can be generated within an inner-stage to promote energy synergy among plants. Different from previous total site integrations, the heat demand of the HEN is unknown before the design stage, but this research did not consider the utilization of multiple types of utilities within the inner-stage. Hipólito-Valencia et al. [
16] combined an organic Rankine cycle with HEN and converted waste heat into mechanical energy, which provided more inspiration for the combination of HEN and other processes. Goh et al. [
17] synthesized a HEN and utility system simultaneously, but the minimum operating cost, hot and cold utility demands were determined in advance through multiple cascade automated targeting, without realizing the overall optimization. Luo et al. [
18] integrated a HEN with a utility system, using the sensible heat of steam condensate to heat cold process streams. Besides, waste heat of the hot process streams was recovered to preheat boiler feedwater, but the combination was allocated at stream ends, which limited the optimization space compared with the superstructure considering inner-stage utility utilization. Martelli et al. [
19] developed a two-stage sequential synthesis algorithm to solve the nonconvex MINLP problem derived from the simultaneous synthesis of HEN and utility systems. An isothermal mixing assumption was made and the utility systems were set as process streams, which will increase the problem-solving complexity. Elsido et al. [
20] have worked on the simultaneous synthesis of utility systems, Rankine cycles and heat exchanger networks. However, their mathematical model was established based on a
p-h superstructure but not a stage-wise superstructure. The main purpose was proposing an ad hoc bilevel decomposition method to improve solution efficiency but not optimizing the utility utilization through improving the superstructure. The simultaneous optimization of multi-plant heat integration using steam as intermediate fluid described in Chang et al. [
21] reflected the utilization of energy within the inner-stage, but the steam used within the inner-stage was generated from process streams, and the utility system was not considered completely here. Huang et al. [
22] extended their research to the simultaneous optimization of a heat exchanger network, steam Rankine cycle and organic Rankine cycle, but the steam utilization and ORC evaporator were placed at stream ends, and inner- and inter-stage improvements were not considered in the superstructure.
Research about HEN synthesis and utility system is summarized here. On one hand, improvements have been made to provide more matching possibilities in HENs and obtain HEN configurations closer to an optimal solution. On the other hand, utility systems and HENs were considered simultaneously to strengthen the heat integration. In total site heat integration using the sequential method, the utility system was usually designed after the determination of heating demand, making it unable to obtain the optimal trade-off between equipment investment and operating cost. As for the mathematical programming method which can perform the design of HENs while optimizing the operation of the utility system, either the heaters were placed at stream ends, or the operating parameters and power generation of utility system were not included, even though the utilization of utilities was considered within the inner-stage. Thus, in order to achieve a better network management, an improved superstructure considering multiple utilities utilization within the inner- and inter-stage is presented and integrated with the utility system based on a Rankine cycle. Multiple utilities (steam in different pressure levels) are produced in cascade in the utility system with concurrently power generation. The structure of the HEN and the operation of the utility system are simultaneously influenced by the placement and distribution of multiple utilities, leading to a trade-off among capital costs, fuel costs and power generation profits. It should be noted that steam is not regarded as process stream but rather as an additional heating source in this study, so the established mathematical model corresponding to the special HEN superstructure including inner- and inter-stage heaters is different from previous formulations, by which the steam in any alternative level can be selected as long as the temperature difference demand is met. Finally, cases are illustrated to show the goals of this paper, demonstrating that multi-level steam selection and utilization within the inner- and inter-stage will provide a larger optimization space for steam distribution, power generation and fuel consumption of the utility system, achieving better economic performance of the whole system.
2. Problem Statement
A synthesis of HEN involving the optimization of steam supply and utilization is desired in this study, and the entire problem can be described as follows: a set of hot and cold process streams which must be cooled and heated are given with their heat capacity flowrates, supply temperatures, target temperatures and heat transfer coefficients. To compensate the heating gap, a utility system is optimized to produce appropriate steam, in terms of multiple paralleled mains/levels with different pressures and as well temperatures. These steam branches are usually classified as high pressure steam, medium pressure steam, and low pressure steam (HPS, MPS, LPS), etc. In this study, a set of steam alternatives with certain pressures and temperatures are given for selection. The utility system follows the basic process of a Rankine cycle, but the properties and quantity of generated steam must be optimized towards the best benefit of the whole system. Additionally, all the cost-related parameters are also given in order to assess cold utility cost, fuel cost, equipment investment (for heat exchangers, boilers and turbines) and the profit from selling power. The integration problem aims at achieving the most cost-efficient configuration of steam supply HEN, by making trade-offs amongst the contribution of the mentioned economic sectors in the total annualized cost. For this purpose, the interaction between the two sub-systems must be comprehensively explored and the network details, such as fuel consumption, steam generation and allocation, generated power, stream matches, heat loads, operating temperatures, and area of heat exchangers must be optimized. In order to alleviate the solving difficulties caused by a complicated mathematical presentation, it is assumed that: (1) the heat capacity and heat transfer coefficient are constant throughout the whole process; (2) the heat exchange task is completed in a countercurrent heat exchanger; (3) all operations are adiabatic, ignoring the heat loss and mass loss; (4) transportation problems are not taken into account; (5) only the primary units in the utility system are considered, ignoring the auxiliary devices.
4. Superstructure
Figure 1 presents the superstructure of the proposed integration problem. All potential configuration alternatives of the HEN are displayed by considering the utilization of steam. For the HEN side, two hot process streams and two cold streams are used to exhibit the features of the stage-wise superstructure. Since the heaters involving steam could be before, after and in parallel with stream-stream heat exchangers, new stages are introduced in the superstructure to fully explore the utilization of steam. It should be noticed that the stages are divided into two types—the inner-stage and the inter-stage—which are alternately arranged. Within each inner-stage, a cold process stream can be split into multiple branches to perform the heat exchange with hot process streams and steam simultaneously. After that, the branches are mixed non-isothermally into a main stream and go into the inter-stage, where only steam is allowed to heat the cold streams. Inter-stage heaters are allowed at both the start and end of cold streams, so the inter-stage number is one more than that of inner-stages. This arrangement could provide much more matching options between cold process streams and steam with limited stages, enriching the configurations of the superstructure. Hot streams are operated in the traditional manner, cooled down by cold process streams within inner-stages and cooling water at stream ends.
The operation of utility system is in typical Rankine cycle: the highest pressure superheated steam is generated in a boiler through combusting fuel, and then sent into a turbine to generate multi-level superheated steam with turbine condensed steam extracted and condensed at the end, producing power at the same time. The superheated steam needs to be desuperheated by mixing with water from the deaerator before being sent to the HEN as saturated steam. After heating the process streams, the condensed water returns to the deaerator, and splits into two parts after the operation, part is sent to the boiler and evaporated into superheated steam, and the residue is mixed with multi-level superheated steam to obtain saturated steam, finishing a cycle. Although the flowsheet of the utility system is quite certain, the possibilities of steam generation and utilization in the HEN need to be determined within the design.
Based on the presented superstructure, a mathematical model can be formulated to determine the optimal network structures and network parameters. In this work, considering the deaerator is an energy-consuming unit, another system without deaerator is also investigated, e.g., using a chemical deaeration technique. The resultant variation in structure is that part of the condensed water is mixed with the corresponding superheated steam, while the rest flows back to the boiler directly.
5. Mathematical Formulation
The synthesis of the overall network is performed with a mathematical model in MINLP. The mass balances, heat balances, feasibility constraints, and formulas for device related and cost related calculations are all given to describe the characteristics contained in the superstructure. The model is decomposed into three parts for understanding: model for HEN synthesis, constraints for the utility system, and the objective function. Implementation of the model is also introduced.
As mentioned in problem statement, the given parameters in the model contain:
- (1)
Stream parameters: supply and target temperatures, flowrates, heat transfer coefficients;
- (2)
Cost-related parameters: unit price of cold utility, fuel and power; cost parameters for heat exchanger, boiler and turbine;
- (3)
Parameters of the steam options: enthalpy of saturate steam, superheated steam, condensate, and water from deaerator.
The notation for these parameters is listed in the Nomenclature section. The rest of the notations in the model indicate the variables that should be determined within model implementation. The binary variables (, , , ) indicate the existences of heat exchangers, determining the structure of the network, and the continuous variables (including inlet and outlet temperatures, stream flowrate, heat load, heat exchange area of each exchanger, quantity of generated steam and power, and cost results) refer to the operation and assessment aspects of the obtained network.
5.1. Model for HEN Synthesis
The model of the HEN is formulated to determine the existences of heat exchangers, operating parameters of heat exchangers, required heat transfer area, capital cost and utility cost. The model formulation presented below consists mainly of mass and energy balance equations. Temperature feasibility constraints are also included. Several sets are defined to assist the expression of involved parameters and variables: i∈I, j∈J denote hot and cold process streams, k∈K denotes temperature interval of inner- and inter-stage, m∈M is the order of temperature location which is corresponding to stage k, n∈N denotes the steam levels in different pressures.
Heat balances of hot and cold process streams are defined by Equations (1) and (2), wherein
FHi and
FCj denote the total heat capacity flowrates of a hot process stream
i and a cold process stream
j;
thini,
thouti,
tcinj,
tcoutj denote supply and target temperatures of streams;
qi,j,k, represents the heat load of a heat exchanger between hot and cold streams;
qhuj,k,n and
qhusj,k,n are the heat satisfied by steam
n within inter-stage and inter-stage, respectively.
qhuj,0,n is the heat satisfied by steam
n at the end of cold stream
j. Equation (1) is used to ensure each hot stream is cooled from its supply temperature to target by performing adequate heat exchanges with cold streams and cold utility. Analogously, Equation (2) is involved to guarantee that the required heat of each cold stream is satisfied by inner-stage stream-stream exchange, inter-stage stream-steam exchange (including the matching at stream end) and inner-stage stream-steam exchange:
Equations (3)–(6) present the energy balances of hot and cold process streams within inner- and inter-stage.
tii,m and
tjj,m are temperature locations of process streams corresponding to stage
k. Equation (3) shows that, within inner-stage, hot streams only exchange heat with cold streams. Equation (4) defines that cold process streams can be heated up by hot process stream and steam
n in parallel form within inner-stage
k, and more than this, the inter-stage heating from steam is also available, as indicated in Equation (5). The relationship of temperature location
m and stage
k is worth mentioning here. In view of the existence of inner- and inter-stage, the presentation of temperature location becomes complicated. In order to enhance the universality of the model, constraints are employed to determine the relationship of
m and
k. For inner-stage
k, the left temperature location
m of process stream corresponds to 2
k−1, right temperature location
m of process stream corresponds to 2
k. While for inter-stage
k, its left temperature location is the same as right temperature of inner-stage
k, so its left temperature location is 2
k, and the right temperature location corresponds to 2
k+1:
Equations (7)–(9) express the mass balances of hot and cold process streams within the inner- and inter-stages. A hot stream must be split into branches to perform the potential heat exchanges with cold streams within each inner-stage (Equation (7)), but for a cold stream, there are two kinds of matches within inner-stage
k, with hot process streams and with steam (Equation (8)). For inter-stages, split branches can only match with steam (Equation (9)):
The characteristics of an exchanger include its location, inlet and outlet temperatures, hot and cold stream (branch) flowrates, heat load and heat transfer area. Due to the variety of the exchangers in the presented superstructure, the models for exchangers are presented in five sets: inner-stage stream-stream exchangers, inter-stage stream-steam exchangers, inner-stage stream-steam exchangers, heaters at end of cold process streams, and coolers at the end of hot streams. All these exchangers are expressed with similar heat balances and area equations in the modeling.
Equations (10)–(16) are used to calculate the heat load
qi,j,k and heat transfer area
Ai,j,k of heat exchangers between hot and cold streams within inner-stage
k. Heat balance Equations (10) and (11) are used to constraint the inlet and outlet temperatures of the unit that, the released heat from hot stream
i equals to that obtained by the pairing cold stream
j.
dtij1
i,j,k and
dtij2
i,j,k in Equations (12) and (13) denote the temperature difference at the two sides of the heat exchanger. In order to circumvent the presence of logarithm and the resultant difficulty in model solving, the logarithmic mean temperature difference used for transfer area calculation (Equation (15)) is approximated by Equation (14) [
23]. In the model, binary variable
zi,j,k is employed to indicate the existence of heat exchangers, that is,
zi,j,k = 1 means the heat exchanger exists, otherwise the exchanger is not involved. Heat loads are highly related to the locations of exchangers, so Equation (16) is used by defining sufficiently large constant
Qmax to assist with the constraint:
In similar way, Equations (17)–(23) are included to describe the steam heaters within each inner-stage. It is stipulated that steam, at most one type, is allowed for each possible stream-steam match, so Equation (22) is introduced into the model, by using
zhusj,k,n to denote the existences of the heaters:
Equations (24)–(30) are formulated to determine the steam heater within each inter-stage.
zhuj,k,n denotes the existence of the heater having steam
n involved. Equation (29) is employed to ensure that a cold process stream can only be heated up by one steam at most within each inter-stage:
Equations (31)–(37) are employed to determine the size of the heaters at the end of cold process streams. As mentioned in the superstructure section, the number of inter stages is one more than that of inner stages, so the exchange for this stage is expressed with the subscript ‘0’:
Equations (38)–(44) present the calculation of heat load
qci and heat transfer area
Acui of the cooler at
ith stream end.
tcuin and
tcuout are the specified inlet and outlet temperatures of cold utility.
zcui denotes the existence of a cooler:
5.2. Constraints for the Utility System
The utility system consists of a boiler, turbine and deaerator. Steam at different pressure levels is extracted from the turbine as a utility. Equation (44) presents the required heat load of steam
n,
qsteamn, equaling to that consumed in inter-stage steam-steam exchanges and inner-stage steam-steam exchanges. The mass flow is deduced from Equation (45), wherein
lheatn is the latent heat of saturated steam
n:
Superheated steam extracted from the turbine is desuperheated by boiler feedwater from the deaerator. Equations (46)–(48) give the mass and heat balances for this process, according to the model of Luo et al. [
18]. Saturated steam is obtained after mixing superheated steam and the boiler feed water which comes from the deaerator. Equations (46) and (47) describe the mass and heat balances of this desuperheating process, wherein
enth denotes the enthalpy of each involved stream/steam. As all condensate water of saturated steam is mixed with exhaust steam in the deaerator at the specified operating temperature of the device, the flow rate of the exhaust steam must be constrained by the balance of the deaerator as indicated in Equation (48), wherein
mext denotes the mass flow of exhaust steam from turbine.
Equations (49)–(52) are used to calculate economy-related factors. Equation (49) gives the power generation model of the steam turbine, wherein
msteamz is the total mass flow through subsection
z of the turbine, equaling the sum of the steam at lower pressure level. As shown, the quantity of generated power is determined by the enthalpy difference of the inlet and outlet superheated steam as well as the steam flow rate through each section. The costs of fuel, turbine and boiler, which are highly-related to steam distribution and power generation, are obtained through Equations (50)–(52) [
15], respectively. Fuel cost is calculated in Equation (50), wherein the required heat load for generating steam from condensate to superheated state is determined by steam flow rate,
, and unit enthalpy difference of the material in and out the boiler,
. Accordingly, the consumption of fuel can be calculated based on the heat content of fuel,
, and boiler efficiency,
. As shown in Equations (51)–(52), the capital cost of the turbine is determined by the quantity of generated power and the capital cost of the boiler depends on the amount of steam generated:
5.3. Objective Function
The aim the study is to obtain an integrated system with the minimum total annual cost (TAC). As presented in the objective function of Equation (53), the TAC equals the total economic cost minus the profit from selling power. Economic cost is composed by the capital cost for exchangers, turbine, boiler, and the operating cost for fuel and cold utility consumption. The exchangers are assessed in general form
CF +
CE ×
Aβ, in which
CF,
CE and
β denote the fixed charge of heat exchanger, coefficient of area cost and exponent for area cost, respectively.
lheatcondenser is latent heat of exhaust steam from turbine end. The costs of turbine, boiler and fuel are obtained according to Equations (50)–(52). The quantities of consumed cold utility and generated power are determined by Equations (38) and (49), then the corresponding cost and profit are calculated by multiplying the unit prices
Ccu and
Cpow:
5.4. Model Implementation
The formulated model is in MINLP presentation, containing a large number of continuous variables, binary variables and non-linear terms. The former kind refers to operational aspects, such as designed load, flowrate, temperature, fuel consumption, steam quantity, generated power and device size, while the later kind account for structure aspect, related to the existence of heat exchange matches and the utilization of steam. Non-isothermal mixing, heat exchange area calculation and cost calculation contribute to the non-linear terms. The model is programmed in the general modeling software GAMS and solved by the solver BARON [
24], within the Windows 10 computing environment using an Intel(R) Core(TM) i7-8,700k CPU @ 3.70 GHz, 32.0 GB RAM. Here it should be noted that the global optimum cannot be always guaranteed if handling big size problems, but appropriate initial values of the variables will facilitate the search to improve the quality of solutions.