Owned by: China Association for Science and Technology
Sponsored by: China Coal Society
Published by: Springer Nature
About issue

The International Journal of Coal Science & Technology is a peer-reviewed open access journal. It focuses on key topics of coal scientific research and mining development, serving as a forum for scientists to present research findings and discuss challenging issues.


Coverage includes original research articles, new developments, case studies and critical reviews in all aspects of scientific and engineering research on coal, coal utilizations and coal mining. Among the broad topics receiving attention are coal geology, geochemistry, geophysics, mineralogy, and petrology; coal mining theory, technology and engineering; coal processing, utilization and conversion; coal mining environment and reclamation and related aspects.


The International Journal of Coal Science & Technology is published with China Coal Society, who also cover the publication costs so authors do not need to pay an article-processing charge.


The journal operates a single-blind peer-review system, where the reviewers are aware of the names and affiliations of the authors, but the reviewer reports provided to authors are anonymous.


  • A forum for new research findings, case studies and discussion of important challenges in coal science and mining development

  • Offers an international perspective on coal geology, coal mining, technology and engineering, coal processing, utilization and conversion, coal mining environment and reclamation and more

  • Published with the China Coal Society

Show More
Editors-in-Chief
Suping Peng, Shimin Liu
Managing Editor
Wanjie Wang
Associate Editors
Bo Hyun Kim, Dongjie Xue, Pedram Roghanchi, Wu Xiao, Zhiqiang Wu
Publishing model
Open Access. Learn about publishing OA with us
Home > Volumes and issues > Volume 8, issue 4

Flexibility analysis and design of heat exchanger network for syngas-to-methanol process

Research Article

Open Access

Published: 16 April 2021

0 Accesses

International Journal of Coal Science & Technology Volume 8, 1468-1478, (2021)

Abstract

The heat exchanger network (HEN) in a syngas-to-methanol process was designed and optimized based on pinch technology under stable operating conditions to balance the energy consumption and economic gain. In actual industrial processes, fluctuations in production inevitably affect the stable operation of HENs. A flexibility analysis of the HEN was carried out to minimize such disturbances using the downstream paths method. The results show that two-third of the downstream paths cannot meet flexibility requirements, indicating that the HEN does not have enough flexibility to accommodate the disturbances in actual production. A flexible HEN was then designed with the method of dividing and subsequent merging of streams, which led to 13.89% and 20.82% reductions in energy consumption and total cost, respectively. Owing to the sufficient area margin and additional alternative heat exchangers, the flexible HEN was able to resist interference and maintain production stability and safety, with the total cost increasing by just 4.08%.

1.Introduction

Methanol is one of the most important raw materials in the chemical industry for producing formaldehyde, acetic acid, methyl formate, etc. It has also been used as fuel for methanol automobiles (Li et al. 2010; Riaz et al. 2013). Along with the rapid development of the coal and chemical industry in recent years, methanol has played an important role in comprehensive coal utilization processes, including coal-to-olefins, coal-to-dimethyl ether, and other coal-to-chemical conversions (Galadima and Muraza 2015; Gao et al. 2018). Generally, syngas from coal gasification technology is used to synthesize methanol via catalytic processes. Coal-based syngas-to-methanol technology has been the main route for methanol production, especially in China, because of its unique situation in terms of energy resources, that is, it is rich in coal but poor in oil and natural gas (Li et al. 2010).

Several studies have focused on the energy utilization efficiency of the coal-based syngas-to-methanol process because of its high energy consumption (Riaz et al. 2013; Bessa et al. 2012; Cui et al. 2017; Rashid et al. 2011; Sun et al. 2012). For example, the entire coal-to-olefins process consists of six subprocesses, including coal mining, transportation, coal-to-methanol (CTM), methanol-to-olefins, product delivery, and carbon capture and storage (CCS) (Gao et al. 2018). The CTM subprocess consumes the most energy, accounting for 71.04% of the total energy consumption, which is five times greater than that of the CCS subprocess (13.81%) (Gao et al. 2018). Extensive efforts have been made to develop energy-saving strategies for the CTM process; among these, crude methanol purification technology has attracted much attention. A five-column heat integrated methanol distillation scheme was proposed by adding a medium-pressure column, which significantly reduced the total energy loss by 21.5% compared with the previous four-column purification scheme (Sun et al. 2012; Chien et al. 2005). A hybrid methanol purification process was developed to improve energy efficiency, which combined heat pump distillation and double-effect thermal integration by designing an intermediate heater to shunt the heat load of the reboiler (Bessa et al. 2012; Douglas and Hoadley 2006). Although these strategies can reduce energy consumption of certain subprocesses in the syngas-to-methanol process, other energy-saving strategies that can improve the energy efficiency of the overall process should also be considered.

The study of heat exchanger networks (HENs) could lead to a comprehensive heat exchange strategy for the entire syngas-to-methanol process by considering energy utilization, energy efficiency, operation cost, equipment cost, etc. Energy consumption can be reduced by using an optimized heat exchange strategy according to the results of global HEN analysis (Kang and Liu 2019). Currently, the optimal HEN can be calculated via pinch technology, which was developed for analyzing potential energy savings and economic gains (Rashid et al. 2011). Previous studies (Kang and Liu 2019; Payet et al. 2018) have reported that the energy consumption and operating cost of the syngas-to-methanol process can be reduced by restructuring the HEN according to the results of pinch analysis. However, the optimization of the HEN has always been carried out under stable heat transfer operating conditions, while ignoring some inevitable uncertainties arising from various external and internal factors (e.g., feed status, product output, heat transfer coefficients, fouling, etc.). To make the optimal HEN from pinch analysis more precise and practical, a flexibility analysis was introduced by considering these uncertain factors under real operating conditions (Payet et al. 2018; Zhu et al. 1996). However, few studies have focused on the flexibility analysis of the HEN for methanol production.

In this work, the syngas-to-methanol process was first simulated with an optimized HEN by pinch technology under the assumption of stable operating conditions. Then, a flexibility analysis of the optimized HEN was conducted using the heuristic method. The “downstream paths” approach was used to determine the variables for building a flexible HEN. The flexible HEN was improved to meet actual production requirements with a sufficient area margin and number of heat exchangers. This work proposes a feasible strategy to reduce the energy consumption, equipment costs, and operating costs of the syngas-to-methanol process by performing a flexibility analysis of an optimized HEN.

2.Methodology and simulation

2.1 Pinch technology

Pinch technology, proposed by Linnhoff and Hindmarsh (1983), is a methodology for minimizing the energy consumption of chemical processes. Pinch technology has been widely used in the design and construction of new HENs as well as in the energy-saving modification of old equipment. The pinch analysis mechanism, based on thermodynamics and computational topology, analyzes the distribution of energy flows and the temperature in a chemical process. The objectives of pinch analysis are to find the heat integration bottleneck, namely “the pinch,” and provide a solution to reduce the energy consumed by a given process (Binosi and Papavassiliou 2009; Tan et al. 2014). The Aspen Energy Analyzer software was utilized to compute and optimize the HEN according to pinch technology.

2.2 Flexible analysis

Generally, the design and optimization of HENs are based on stable operating conditions. However, uncertain factors affect the stability of HENs in actual production. Uncertain factors (e.g., feed status fluctuation) can cause the operation of the HEN to deviate from the optimal state or even cease altogether. Thus, in the design and optimization of the HEN, it is necessary to meet the flexibility requirements under fluctuating operating conditions, instead of simply considering stable operating conditions. The concept of the flexible index was first proposed by Swaney and Grossmann in 1987 to measure the maximum deviation of uncertain factors from stable operating conditions in the feasible region; a strict mathematical model was also established (Swaney and Grossmann 1985). The flexibility analysis includes mathematical programming and heuristic methods (Grossmann and Kravanja 1995; Linnhoff and Kotjabasakis 1986). The mathematical programming method is frequently used to study the feasibility of HENs in a specific disturbed domain under certain constraints, while the heuristic method is mainly applied in the detailed structural analysis of HENs to determine the feasibility in the entire disturbed domain based on the experience in chemical engineering (Grossmann and Kravanja 1995; Linnhoff and Kotjabasakis 1986). The heuristic method is much simpler and more convenient than the mathematical programming method, which is widely used in the field of chemical engineering, particularly in relatively short chemical processes. Moreover, a feasible and precise strategy for HENs can be obtained by the heuristic method, whereas the mathematical programming method is only used to test the feasibility of existing HENs under different disturbance conditions. In the heuristic method, the concept of downstream paths was developed to find a balance between the flexibility and economics of HENs. Downstream paths have been used as a convenient approach for determining the disturbed and controlled variables of HENs during flexibility analysis (Linnhoff and Kotjabasakis 1986). A path is defined as the uninterrupted connection between any two nodes in the grid diagram of a HEN. A downstream path is defined as the directed path that is always in the same direction as the streams through which it passes (Zhu et al. 1996). Thus, the disturbance occurring in the HEN can only affect the controlled variables along with the downstream paths between the disturbed and controlled variables.

2.3 Process simulation

The technical route adapted in this study is illustrated in Fig. S1 in Supplementary Material. The process simulation was performed using Aspen Plus V10. The original HEN, traditional optimal HEN, and flexible HEN were designed and optimized using the Aspen Energy Analyzer V10 based on the theories of pinch technology and heuristic methods. The entire syngas-to-methanol process was divided into two parts for the simulation: methanol synthesis and methanol purification, as shown in Fig. 1. The information of all blocks and streams are listed in Table S1 in Supplementary Material. The Lurgi low-pressure methanol synthesis process was adopted to simulate the methanol synthesis process. The synthetic reaction was carried out at 250 °C and 5 MPa in a Lurgi reactor, modeled by combining the RPLUG and RSTOIC reactor blocks in Aspen Plus. RPLUG and RSTOIC were used to model the two main reactions and six side reactions, respectively. The reaction equations and kinetic parameters are listed in Table S2 in Supplementary Material. The equipment parts in which methanol synthesis and methanol purification processes are performed were connected by stream PCF-IN. In the methanol purification process, small quantities of lightweight components (e.g., dimethyl ether and methyl formate) dissolved in crude methanol were removed in the predistillation column (Block PC). Double-effect distillation, which combined a high-pressure column (Block HC) and an atmospheric column (Block CC), was applied to remove impurities, such as water, ethanol, and n-butanol. The purity of the obtained methanol product reached 99.98 wt%. The compositions and volume fractions of syngas and crude methanol are listed in Table S3 and S4 in Supplementary Material, respectively.

Fig. 1
figure 1

Simulation flowsheet of syngas-to-methanol

3.Design and optimization of HEN

3.1 Determination of energy targets

The datasheet of the process simulation obtained by Aspen Plus was imported directly into the Aspen Energy Analyzer to design the HEN. The temperature and enthalpy of each stream participating in heat exchange are listed in Table 1. Through the computations of the Aspen Energy Analyzer, the total cost index was maintained at a low level when the minimum heat transfer temperature difference (ΔTmin) ranged between 4 and 11 °C (Fig. S2 in Supplementary Material). Here, the value of ΔTmin was set to 11 °C because the empirical value of ΔTmin is generally within the range of 10–20 °C by considering factors such as utilities and heat exchanger equipment cost, heat exchange media, heat transfer coefficient, and operation flexibility (Kang and Liu 2019). The composite curves and grand composite curves are shown in Fig. S3 and S4 in Supplementary Material, respectively. The pinch-point temperatures of the hot and cold streams were 135 °C and 124 °C, respectively, from the composite curves (see Fig. S3 in Supplementary Material) when ΔTmin was set to 11 °C. The minimum hot utility (QH, min) target was 1.460 × 108 kJ/h, while the minimum cold utility (QC, min) target was 2.396 × 108 kJ/h. According to Fig. S4 in Supplementary Material, the average temperature to process the pinch point was confirmed as 129.5 °C when the corresponding enthalpy was 0 kJ/h.

Table 1 Datasheet of streams involving heat exchange

Steam No.

Steam name

Tin (°C)

Tout (°C)

Enthalpy (kJ/h)

H1

REA-OUT_To_28

250.0

40.0

1.52 × 108

H2

HCW-OUT_To_CCF-IN

129.2

80.3

5.86 × 106

H3

HCD-OUT_To_49

122.3

121.8

1.46 × 108

H4

24_To_16

122.3

40.0

5.58 × 106

H5

27_To_W-WATWER

107.3

40.0

1.80 × 106

H6

To Condenser@CC_TO_CCD-OUTDuplicate

74.9

71.0

1.46 × 108

H7

PCD-OUT_To_S1-IN

73.3

64.0

1.55 × 107

H8

CCD-OUT_To_31

71.0

40.0

2.22 × 106

H9

S1D-OUT_To_S2-IN

64.0

40.0

1.71 × 106

C1

10_To_REA-IN

124.0

230.0

5.51 × 107

C2

To Reboiler@HC_TO_HCW-OUTDuplicate

126.9

129.2

1.46 × 108

C3

21_To_HCF-IN

81.7

127.0

8.76 × 106

C4

52_To_CCW-IN

107.3

110.9

1.46 × 108

C5

To Reboiler@PC_TO_PCW-OUT

81.1

81.7

1.85 × 107

C6

23_To_PCF-IN

20.4

75.0

8.65 × 106

3.2 Original HEN and traditional optimal HEN

The original HEN of the syngas-to-methanol process generated by the Aspen Energy Analyzer is shown in Fig. S5 in Supplementary Material. The original HEN includes two process-process heat exchangers (white matches), nine coolers (blue matches), and five heaters (red matches). The original HEN was optimized under the assumption of stable operating conditions to obtain the traditional optimal HEN, as shown in Fig. S6 in Supplementary Material. The optimized HEN contains seven process-process heat exchangers, seven coolers, and two heaters.

3.3 Flexible HEN

3.3.1 Flexibility analysis of HEN

The disturbed and controlled variables of the HEN in the syngas-to-methanol process were determined using the downstream paths approach. Disturbed variables fluctuate owing to changes in the external environment, equipment operation, production load, etc. In the actual production process, the temperatures of several streams were determined as the disturbed variables, i.e., the feedstock stream and outlet stream of the reactor; the mixing stream of crude methanol with water, and the bottom outlet stream of the predistillation column (Rashid et al. 2011). Here, θ is denoted as the disturbed variables; θ1 is the temperature of the feedstock stream C1; θ2 is the temperature of the outlet stream (H1) of the reactor; θ3 is the temperature of the mixing stream (C6) of crude methanol with water, and θ4 is the temperature of the bottom outlet stream (C3) of the predistillation column. Controlled variables are defined as those that are strictly manipulated to maintain stability and ensure safety in the actual production process. In the syngas-to-methanol system, the inlet stream temperature of the reactor must be strictly controlled to ensure the stability of the production load and feedstock conversion rate. In addition, control of the feed stream temperature of the pressured column and the atmospheric column is required owing to the sensitive inlet temperature requirements for double-effect distillation. Therefore, the temperatures of (1) the reactor inlet stream and (2) the feed streams of the pressured column and atmospheric column were set as the controlled variables. Here, z is the controlled variable; z1 is the temperature of the reactor inlet stream (C1); z2 is the temperature of the feed stream (C3) of the pressured column, and z3 is the temperature of the feed stream (H2) of the atmospheric column. Figure 2 shows the positions of the disturbed and controlled variables in the grid diagram of the traditional optimal HEN.

Fig. 2
figure 2

Disturbed variables and controlled variables in the grid diagram of traditional optimal HEN

A method for determining downstream paths was developed using the node adjacency matrix that represents the connection of units along the direction of the HEN streams (Zhu et al. 1996). Units, splitting, and mixing points as well as various variables in the HEN were all considered as nodes. In the adjacency matrix, both the row and column headings represent the abovementioned nodes. The (i, j)th entry of the matrix is filled with "*" if node j can directly reach node i in the direction of the stream; otherwise, it is left blank. The nodes of the coolers or heaters in the matrix, which were not adjacent to other nodes, were deleted for the sake of simplicity (Table 2). In this case, the determination of the downstream paths between the disturbed and controlled variables can be transformed into the determination of the downstream paths between the two nodes, i and j. The detailed approach for finding the downstream path is described as follows. (Step I) θi (i = 1, 2, …, I) and zj (j = 1, 2, …, J) represent the number of disturbed and controlled variables of the HEN, respectively. (Step II) S is defined as a set of nodes, and the serial number of these nodes should start from i = 1. (Step III) For the last node p in S, the node q that can be directly reached from node p along with the direction of the stream must be confirmed, i.e., if the (p, q)th entry is filled with "*" in the node adjacency matrix, then add q to set S. Step III will stop if S contains all the nodes of the HEN, or the (p, q)th entry is blank in the node adjacency matrix for each q. If node zj is included in set S, then there is a downstream path between node θi and node zj; otherwise, let i = i + 1 and go back to Step II until i > I. For the traditional optimal HEN, the results of identifying the downstream paths between all disturbed and controlled variables are given in Table 3; each downstream path is plotted in Fig. 3.

Table 2 Node adjacency matrix of traditional optimal HEN

j\i

1

2

3

4

5

6

7

8

9

HE1

HE2

CE1

CE2

θ1

θ2

θ3

θ4

z1

z2

z3

1

*

        

*

    

*

     

2

*

*

        

*

         

3

  

*

    

*

        

*

   

4

   

*

                

5

    

*

  

*

            

6

     

*

*

             

7

      

*

        

*

    

8

 

*

     

*

            

9

  

*

 

*

   

*

           

HE1

         

*

   

*

      

HE2

          

*

         

CE1

        

*

  

*

        

CE2

      

*

     

*

       

θ1

             

*

      

θ2

              

*

     

θ3

               

*

    

θ4

                

*

   

z1

*

                

*

  

z2

  

*

               

*

 

z3

     

*

             

*

Table 3 Downstream paths between disturbed and controlled variables

Number of downstream paths

Initial node

Process nodes

Final node

1

θ1

HE1, 1

z1

2

θ1

HE1, 1, 2, 8, 3

z2

3

θ2

1

z1

4

θ2

1, 2, 8, 3

z2

5

θ3

7, 6

z3

6

θ4

3

z2

Fig. 3
figure 3

Downstream paths between disturbed variables and controlled variables

3.3.2 Flexibility design of HEN

For an HEN with temperature fluctuations, heat recovery reaches a maximum when the inlet temperatures of all hot streams take the upper limit of fluctuations, while the inlet temperatures of all cold streams take the lower limit. With this assumption, it is possible to obtain the maximum design margin, and subsequently meet the requirements for various operating conditions in actual production. The stream-dividing principle was used to design a flexible HEN (Zhihong and Ben 1999). The details are as follows. (i) For hot stream i, its outlet temperature is denoted as Ti, out, while the inlet temperature is in the range of Ti, minTi, max. Hot stream i can be divided into the high-temperature stream and low-temperature stream. In the high-temperature stream, the inlet and outlet temperatures are Ti, max and Ti, min, respectively. In the low-temperature stream, the inlet and outlet temperatures are Ti, min and Ti, out, respectively. (ii) Similarly, for cold stream j, its outlet temperature is referred to as Tj, out, and the inlet temperature is in the range of Tj, minTj, max. Cold stream j can be divided into the high-temperature stream and low-temperature stream. In the high-temperature stream, the inlet and outlet temperatures are Ti, max and Tj, respectively. In the low-temperature stream, the inlet temperature is Ti, min, and the outlet temperature is Tj, max. The fluctuation ranges of the disturbed variables were determined according to the actual production, as listed in Table 4 (Zamora and Grossmann 1997). θN was set as the reference value. Δθ+ and Δθ are the expected deviations.

Table 4 Reference values of disturbed variables and their expected deviation

Disturbance variable

θN (°C)

Δθ+ (°C)

Δθ (°C)

θ1 = TC1

124.0

10

10

θ2 = TH1

250.0

10

10

θ3 = TC6

20.4

10

0

θ4 = TC3

82.7

5

5

The process streams with disturbed variables were divided according to the dividing principle; the data are listed in Table 4. The segmented streams were treated as independent streams. Then, the data of all streams shown in Table 5 were imported into the Aspen Energy Analyzer to design the flexible HEN. The final flexible HEN was obtained after merging the segmented streams and relaxing the energy load. The detailed procedure is as follows. (i) The segmented process streams were first merged to obtain the network structure shown in Fig. 4. (ii) Then, the load loops were relaxed along the direction of the downstream paths of the merged streams. (iii) Finally, the flexible HEN was obtained by adjusting the parameters of each heat exchanger, as shown in Fig. 5. The merging of segmented streams was conducted based on the design results under the condition of maximum heat recovery. Energy relaxation can reduce the number of heat exchangers in the entire process. Consequently, the final flexible HEN is optimized.

Table 5 Datasheet of process streams obtained by the dividing principle

New No.

Old No.

Tin (°C)

Tout (°C)

q (kJ/h)

h1

H1g

260

240

9.72 × 106

h2

H1d

240

40

1.48 × 108

h3

H2

129.2

80.3

5.86 × 106

h4

H3

122.3

121.8

1.46 × 108

h5

H4

122.3

40

5.58 × 106

h6

H5

107.3

40

1.80 × 106

h7

H6

74.9

71

1.46 × 108

h8

H7

73.3

64

1.55 × 107

h9

H8

71

40

2.22 × 106

h10

H9

64

40

1.71 × 106

c1

C1g

134

230

4.99 × 107

c2

C1d

114

134

1.02 × 107

c3

C2

126.8

129.8

1.46 × 108

c4

C3g

86.7

127

7.85 × 106

c5

C3d

76.7

86.7

1.76 × 106

c6

C4

107.3

110.9

1.46 × 108

c7

C5

81.1

81.7

1.85 × 107

c8

C6g

30.4

75

7.13 × 106

c9

C6d

20.4

30.4

1.52 × 106

Fig. 4
figure 4

Grid diagram of flexible HEN with merging the segmented streams

Fig. 5
figure 5

Final flexible HEN

4.Results and discussion

From Fig. S5 in Supplementary Material, the strategy of the original HEN employed the excessive number of 16 utilities, which can be attributed to two reasons. One is the situation in which the cooler crosses two pinch points and the heater crosses the process pinch point, which violates the rule that prohibits a cross-pinch. The other is the presence of three heaters below the pinch point and a cooler above the pinch point, which violates another rule of pinch technology: no heating utility below the pinch point and no cooling utility above the pinch point. It can be concluded that the original HEN contains massive energy that is expected to be recycled. Compared with the original HEN scenario, the traditional optimal HEN significantly improves the heat recovery capacity of the system; the number of heat exchangers is also significantly reduced. Although two process-process heat exchangers cross the pinch point, there is no cooling and heating utility above and below the pinch point, respectively (Fig. S6 in Supplementary Material). In the flexible HEN (Fig. 4), a total of six downstream paths were obtained. The results of the flexibility analysis can be summarized as follows. (i) For downstream paths 1 and 2 with the disturbance variable θ1 as the initial node, it is possible to stop θ1 from affecting the controlled variables z1 and z2 by adjusting the usage of the heating utility because these two paths pass through node HE1. (ii) For downstream paths 3 and 6, there is only one node of a heat exchanger between the disturbance and control variables, with the former having a direct effect on the latter. Some differences between paths 3 and 6 were observed. Both variables on downstream path 3 were located in different streams. In contrast, the variables on downstream path 6 were in the same stream. (iii) For downstream paths 4 and 5, there are many nodes on each but no utilities node; hence, the fluctuation of the disturbance variable will be transferred to the controlled variable along the downstream path direction. In summary, the controlled variables will be affected as long as the disturbance variables of downstream paths 3–6 fluctuate. Thus, it is difficult for this regular optimized HEN to meet the flexibility requirements in a practical scenario.

Table 6 shows the comparison of major parameters of the original HEN and flexible HEN. The total energy consumption of the flexible HEN decreased by 13.89% compared with that of the original HEN, indicating a high energy recycling efficiency. The capital cost, operating cost, and total cost in the flexible HEN scenario decreased by 39.60%, 11.53%, and 20.82%, respectively.

Table 6 Parameters of original HEN and flexible HEN

Parameter

Original HEN

Flexible HEN

Reduction (%)

Heating utilities (kJ/h)

1.819 × 108

1.508 × 108

17.12

Cooling utilities (kJ/h)

2.756 × 108

2.432 × 108

11.76

Number of shells

55

37

32.73

Total area (m2)

2.066 × 104

1.191 × 104

42.35

Capital cost (Cost)

5.099 × 106

3.080 × 106

39.60

Operating cost (Cost/s)

0.1053

0.0931

11.53

Total cost (Cost/s)

0.1573

0.1246

20.82

Table 7 shows a comparison of the major parameters of the traditional optimal HEN and flexible HEN. The cooling and heating utilities in the flexible HEN scenario increased slightly at increments of 2.72% and 1.13%, respectively, compared with the regular optimized HEN. The number and total area of the heat exchangers increased by 15.63% and 29.35%, respectively. This illustrates that the essence of the HEN flexibility design in our work is to increase the area margin and appropriately set up alternative heat exchangers. However, the operating cost in the flexible HEN strategy decreased by 1.69%, implying that increasing the HEN flexibility can reduce the frequency and difficulty of operations. The total cost slightly increased (4.08%). In the long term, the extended lifespan of heat exchangers will offset the partial increase in capital costs stemming from the additional heat exchangers and area margin.

Table 7 Parameters of traditional optimal HEN and flexible HEN

Parameter

Regular optimized HEN

Flexible HEN

Reduction (%)

Heating utilities (kJ/h)

1.468 × 108

1.508 × 108

− 2.72

Cooling utilities (kJ/h)

2.405 × 108

2.432 × 108

− 1.13

Number of shells

32

37

− 15.63

Total area (m2)

9209

1.191 × 104

− 29.35

Capital cost (Cost)

2.445 × 106

3.080 × 106

− 25.97

Operating cost (Cost/s)

0.0947

0.0931

1.69

Total cost (Cost/s)

0.1197

0.1246

− 4.08

5.Conclusions

In this study, the overall HEN of a syngas-to-methanol process was designed and optimized based on pinch technology under stable operating conditions. Flexibility analysis of the traditional optimal HEN was implemented using the downstream paths method. Six downstream paths between the disturbance variables and control variables were determined by the modified procedure with the node adjacency matrix. The results show that two-thirds of the downstream paths cannot meet flexibility requirements, indicating that the optimized design of the HEN under stable operating conditions is not flexible enough to resist disturbances in actual production. The flexible HEN was realized under the condition of maximum heat recovery by dividing and subsequent merging of streams. Although the total cost of the flexible HEN is 4.08% higher than that of the traditional optimal HEN, the flexible HEN has a sufficient area margin and extra alternative heat exchangers to avoid production disturbances and maintain its stability and safety.

References

[1] Bessa LCBA, Batista FRM, Meirelles AJA (2012) Double-effect integration of multicomponent alcoholic distillation columns. Energy 45:603–612. https://doi.org/10.1016/j.energy.2012.07.038
[2] Binosi D, Papavassiliou J (2009) Pinch technique: theory and applications. Phys Rep 479:1–152. https://doi.org/10.1016/j.physrep.2009.05.001
[3] Chien IL, Teng Y-P, Huang H-P, Tang YT (2005) Design and control of an ethyl acetate process: coupled r eactor/column configuration. J Process Control 15:435–449. https://doi.org/10.1016/j.jprocont.2004.07.003
[4] Cui C, Sun J, Li X (2017) A hybrid design combining double-effect thermal integration and heat pump to the methanol distillation process for improving energy efficiency. Chem Eng Process 119:81–92. https://doi.org/10.1016/j.cep.2017.06.003
[5] Douglas AP, Hoadley AFA (2006) A process integration approach to the design of the two- and three-column methanol distillation schemes. Appl Therm Eng 26:338–349. https://doi.org/10.1016/j.applthermaleng.2005.07.001
[6] Galadima A, Muraza O (2015) From synthesis gas production to methanol synthesis and potential upgrade to gasoline range hydrocarbons: a review. J Nat Gas Sci Eng 25:303–316. https://doi.org/10.1016/j.jngse.2015.05.012
[7] Gao D, Qiu X, Zhang Y, Liu P (2018) Life cycle analysis of coal based methanol-to-olefins processes in China. Comput Chem Eng 109:112–118. https://doi.org/10.1016/j.compchemeng.2017.11.001
[8] Grossmann IE, Kravanja Z (1995) Mixed-integer nonlinear programming techniques for process systems engineering. Comput Chem Eng 19:189–204. https://doi.org/10.1016/0098-1354(95)00072-A
[9] Kang L, Liu Y (2019) Synthesis of flexible heat exchanger networks: a review. Chin J Chem Eng 27:1485–1497. https://doi.org/10.1016/j.cjche.2018.09.015
[10] Li Z, Gao D, Chang L, Liu P, Pistikopoulos EN (2010) Coal-derived methanol for hydrogen vehicles in China: energy, environment, and economic analysis for distributed reforming. Chem Eng Res Des 88:73–80. https://doi.org/10.1016/j.cherd.2009.07.003
[11] Linnhoff B, Hindmarsh E (1983) The pinch design method for heat exchanger networks. Chem Eng Sci 38:745–763. https://doi.org/10.1016/0009-2509(83)80185-7
[12] Linnhoff B, Kotjabasakis E (1986) Downstream paths for operable process design. Chem Eng Prog 5:263–281
[13] Payet L, Thery Hétreux R, Hétreux G, Bourgeois F, Floquet P (2018) Flexibility assessment of heat exchanger networks: from a thorough data extraction to robustness evaluation. Chem Eng Res Des 131:571–583. https://doi.org/10.1016/j.cherd.2017.11.036
[14] Rashid SRA, Ibrahim UK, Alauddin SM (2011) Retrofit design of heat exchanger network (HEN) on synthesis and purification unit of methanol plant. In: 2011 IEEE Colloquium on Humanities, Science and Engineering Research (CHUSER 2011) Dec 5–6 2011, pp 33–36. https://doi.org/10.1109/CHUSER.2011.6163746
[15] Riaz A, Zahedi G, Klemeš JJ (2013) A review of cleaner production methods for the manufacture of methanol. J Clean Prod 57:19–37. https://doi.org/10.1016/j.jclepro.2013.06.017
[16] Sun J, Wang F, Ma T, Gao H, Wu P, Liu L (2012) Energy and exergy analysis of a five-column methanol distillation scheme. Energy 45:696–703. https://doi.org/10.1016/j.energy.2012.07.022
[17] Swaney RE, Grossmann IE (1985) An index for operational flexibility in chemical process design. Part I Formulation and theory. AIChE J 31:621–630. https://doi.org/10.1002/aic.690310412
[18] Tan YL, Ng DKS, El-Halwagi MM, Foo DCY, Samyudia Y (2014) Floating pinch method for utility targeting in heat exchanger network (HEN). Chem Eng Res Des 92:119–126. https://doi.org/10.1016/j.cherd.2013.06.029
[19] Zamora JM, Grossmann IE (1997) A comprehensive global optimization approach for the synthesis of heat exchanger networks with no stream splits. Comput Chem Eng 21:S65–S70. https://doi.org/10.1016/S0098-1354(97)87480-7
[20] Zhihong L, Ben H (1999) Synthesis of flexible heat exchanger networks with stream no-splitting(I) based on ranges of stream supply temperatures. J Chem Ind Eng 50:317–325. https://doi.org/10.3321/j.issn:0438-1157.1999.03.005
[21] Zhu J, Han Z, Rao M, Chuang KT (1996) Identification of heat load loops and downstream paths in heat exchanger networks. Can J Chem Eng 74:876–882. https://doi.org/10.1002/cjce.5450740609

About this article

Cite this article

Liu, J., Zhang, P., Xie, Q. et al. Flexibility analysis and design of heat exchanger network for syngas-to-methanol process.Int J Coal Sci Technol 8, 1468–1478 (2021).
  • Received

    29 September 2020

  • Revised

    06 January 2021

  • Accepted

    25 March 2021

  • Issue Date

    December 2021

  • DOI

    https://doi.org/10.1007/s40789-021-00426-4

  • Share this article

    Copy to clipboard

For Authors

Explore