CN101777083A - Heat exchange network optimization synthesis method based on Delta T-contribution value correction - Google Patents
Heat exchange network optimization synthesis method based on Delta T-contribution value correction Download PDFInfo
- Publication number
- CN101777083A CN101777083A CN200910155061A CN200910155061A CN101777083A CN 101777083 A CN101777083 A CN 101777083A CN 200910155061 A CN200910155061 A CN 200910155061A CN 200910155061 A CN200910155061 A CN 200910155061A CN 101777083 A CN101777083 A CN 101777083A
- Authority
- CN
- China
- Prior art keywords
- heat exchange
- unit
- heat
- network
- contribution margin
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Abstract
The invention relates to a heat exchange network optimization synthesis method based on Delta T-contribution value correction, which comprises the following steps: (1) heat exchange fluid heat transfer temperature difference contribution value Delta T is calculated and corrected; and (2) with thermodynamic analysis as the basis and the starting point, the synchronous integration of the energy target, area target, the cost target and the network structure design of a heat exchange network is realized through a mathematical model building method, and the specific process is as follows: (2.1) the corrected Delta T-contribution value corrects the temperature of all heat exchange units, a temperature enthalpy diagram of change pinches is drawn, and the pinch temperature and the energy target after correction are obtained; (2.2) a change balance compound curve is created on the temperature enthalpy diagram, an enthalpy interval is segmented, and the super-structure model of the heat exchange network is built; and (2.3) minimum total cost serves as the objective function, the super-structure is solved, and the optimal heat exchange network is obtained. The heat exchange network optimization synthesis method based on Delta T-contribution value correction can effectively prevent arbitrary temperature correction, and has good practicality and systematicness.
Description
Technical field
The present invention relates to a kind of heat exchanger network synthetic method, especially a kind of heat exchange network optimization synthesis method based on modified Δ T-contribution margin.
Background technology
Heat exchanger network is as an important subsystem of procedures system, be used for the available energy of recovery system, minimizing is to the demand of outside public work, be one of chief component of process industrials such as chemical industry, oil refining, for the energy utilization of production run with energy-saving and cost-reducingly have a crucial meaning.
The optimization of heat exchanger network comprehensively is exactly to find the solution satisfying each strand process stream by the initial temperature heating or be cooled under the prerequisite of intended target temperature, have the minimum equipment investment expense and a heat exchanger network structure of operation cost.
The synthetic method of heat exchanger network mainly is divided into folder point analysis method, mathematical programming approach and artificial intelligence approach at present.Most popular in engineering practice is the folder design method of points, and the folder point design method can be used as the important foundation of other optimization method of development based on thermodynamics.The heat exchanger network composition problem that has the heat exchanging fluid of different film heat-transfer coefficients for solution has the scholar to propose Δ T-contribution margin method, and Δ T-contribution margin can calculate with formula (3):
In formula (3), Δ T
jBe the temperature difference contribution margin of stream gang j, h
jBe the heat transfer coefficient of stream gang j, k and z are empirical parameters.For given parameter k, z, the change of every stream burst temperature equals its corresponding Δ T-contribution margin size.After determining the Δ T-contribution margin of each stream thigh, each stream burst revised temperature is:
Wherein, Δ T
H, iWith Δ T
C, jBe respectively the Δ T-contribution margin of hot-fluid thigh i and cold flow thigh j; T
H, iAnd T
C, jBe respectively initially the advancing of hot-fluid thigh i and cold flow thigh j (going out) mouthful temperature; T
* H, iAnd T
* C, jBe respectively that hot-fluid thigh i and cold flow thigh j are revised to advance (going out) mouthful temperature.
The solution procedure of Δ T-contribution margin method is similar to tradition folder point analysis method, and difference is that it is based on temperature T change or that change
* H, iAnd T
* C, j, rather than actual temperature T
H, iAnd T
C, jIn addition, during utilization Δ T-contribution margin method, a crucial step is to parameter k, the optimization of z, and the value of parameter k and z and unrestricted.This just causes existing heat exchange network optimization method based on Δ T-contribution margin to exist temperature difference correction randomness, heat exchanger network design to lack deficiencies such as systematicness, therefore develops energy-conservation significant to process industrial of more efficiently heat exchange network optimization and synthetic method.Still there is not relevant achievement at present based on the heat exchange network optimization synthesis method aspect of modified Δ T-contribution margin.
Summary of the invention
For the temperature difference correction randomness that overcomes existing heat exchanger network synthetic method, lack systematic deficiency, the invention provides and a kind ofly effectively avoid temperature difference correction randomness, have good practicality and systematic heat exchange network optimization synthesis method based on the correction of Δ T-contribution margin.
The technical solution adopted for the present invention to solve the technical problems is:
A kind of heat exchange network optimization synthesis method based on the correction of Δ T-contribution margin may further comprise the steps:
1), calculate heat exchanging fluid heat transfer temperature difference contribution margin Δ T, and contribution margin Δ T is revised, innovation representation is:
ΔT
H，i+ΔT
C，j≤T
H，i-T
C，j (1)
In the following formula (1): Δ T
H, iBe the heat transfer temperature difference contribution margin of hot fluid i, unit is ℃; Δ T
C, jBe the heat transfer temperature difference contribution margin of cold fluid j, unit is ℃; T
H, iBe the inlet temperature of hot fluid i, unit is ℃; T
C, jBe the inlet temperature of cold fluid j, unit is ℃;
2), based on thermodynamic analysis and starting point, adopt the method for mathematical modeling carry out the energy object of heat exchanger network, area target, cost objective and network structure design synchronously not comprehensively, detailed process is as follows:
(2.1) adopt revised Δ T-contribution margin that the temperature of each heat exchange stream thigh is revised, draw the warm enthalpy diagram that becomes the folder point, obtain revised folder point temperature and energy object;
(2.2) on warm enthalpy diagram, construct the weighing apparatus composite curve chart that flattens, cut apart the enthalpy interval, the superstructure model of structure heat exchanger network;
(2.3) be objective function with the total cost minimum, find the solution superstructure, obtain best heat exchanger network; The expression formula of objective function is:
TAC＝(Q
Hmin·C
H)+(Q
Cmin·C
C)+C
CAP (2)
In the following formula (2): TAC is the total annual charge in heat exchange networking, and unit is unit/year; Q
HminAnd Q
CminBe respectively the minimum thermal of trying to achieve in the described step (2.1), cold public work consumption, unit is J/; C
HAnd C
CBe respectively the heat of using in the heat exchanger network, the unit costs of cold public work, unit is unit/J; C
CAPBe the total annualized capital cost usefulness of network, unit is unit/year.
Further, in described step (2.3), when finding the solution the superstructure model, regulation MINIMUM HEAT-TRANSFER AREA constraint condition makes and rejects little load, small size heat interchanger automatically in the solution procedure of heat exchanger network superstructure, realizes the lax of heat exchanger network.
Further again, in described step (2.3), at the interval k of a certain given sub-enthalpy, if I bar hot-fluid thigh and J bar cold flow thigh are arranged, so in this sub-enthalpy interval, total (the individual coupling number of I * J), and the operative constraint equation number is (I+J-1), when I * J＞I+J-1, in determining this sub-enthalpy interval, need carry out inside son optimization during Load Distribution, the investment cost in the sub-for this reason enthalpy of its optimization aim interval minimizes.
Technical conceive of the present invention is: the systems approach that adopts thermodynamic analysis and mathematical computations to combine, at first,, traditional heat exchanging fluid heat transfer temperature difference contribution margin Δ T is revised based on thermodynamic analysis; Then, the method that adopts folder point analysis and mathematics to combine, the structural design of energy object, area target, cost objective and the network of heat exchanger network is comprehensive synchronously; At last, realized that by setting up minimum area constraint condition heat exchanger network relaxes and carrying out synchronously of optimizing.
Beneficial effect of the present invention is: 1, this method is by setting up temperature difference correction restrictive condition, avoided original method since the temperature difference modified value that temperature difference correction randomness is caused greater than the appearance of this unreasonable situation of actual heat transfer temperature difference; 2, the thought that combines with mathematics based on thermodynamics of this method adopts the multi-level optimization strategy, has better practicability and systematicness; 3, compare with classic method, the minimum area constraint has realized relaxing and optimization is comprehensive synchronously by setting up, and finds the solution the efficient height.
Description of drawings
Fig. 1 is the building-block of logic of this method.
Fig. 2 is the synoptic diagram of the heat exchanger network after optimizing.
Embodiment
Below in conjunction with accompanying drawing the present invention is further described.
With reference to Fig. 1, a kind of heat exchange network optimization synthesis method based on the correction of Δ T-contribution margin may further comprise the steps:
1), calculate heat exchanging fluid heat transfer temperature difference contribution margin Δ T, and contribution margin Δ T is revised, innovation representation is:
ΔT
H，i+ΔT
C，j≤T
H，i-T
C，j (1)
In the following formula (1): Δ T
H, iBe the heat transfer temperature difference contribution margin of hot fluid i, unit is ℃; Δ T
C, jBe the heat transfer temperature difference contribution margin of cold fluid j, unit is ℃; T
H, iBe the inlet temperature of hot fluid i, unit is ℃; T
C, jBe the inlet temperature of cold fluid j, unit is ℃;
2), based on thermodynamic analysis and starting point, adopt the method for mathematical modeling carry out the energy object of heat exchanger network, area target, cost objective and network structure design synchronously not comprehensively, detailed process is as follows:
(2.1) adopt revised Δ T-contribution margin that the temperature of each heat exchange stream thigh is revised, draw the warm enthalpy diagram that becomes the folder point, obtain revised folder point temperature and energy object;
(2.2) on warm enthalpy diagram, construct the weighing apparatus composite curve chart that flattens, cut apart the enthalpy interval, the superstructure model of structure heat exchanger network;
(2.3) be objective function with the total cost minimum, find the solution superstructure, obtain best heat exchanger network; The objective function expression formula is:
TAC＝(Q
Hmin·C
H)+(Q
Cmin·C
C)+C
CAP (2)
In the following formula (2): TAC is the total annual charge in heat exchange networking, and unit is unit/year; Q
HminAnd Q
CminBe respectively the minimum thermal of trying to achieve in the described step (2.1), cold public work consumption, unit is J/; C
HAnd C
CBe respectively the heat of using in the heat exchanger network, the unit costs of cold public work, unit is unit/J; C
CAPBe the total annualized capital cost usefulness of network, unit is unit/year.
In the present embodiment, at first,, find the solution heat transfer temperature difference contribution margin Δ T based on thermodynamic analysis, and heat exchanging fluid heat transferring temperature difference contribution margin Δ T additional corrections condition; Then, the method that adopts folder point analysis and mathematics to combine, the structural design of energy object, area target, cost objective and the network of heat exchanger network is comprehensive synchronously; At last, realized that by setting up minimum area constraint condition heat exchanger network relaxes and carrying out synchronously of optimizing.The specific implementation step is as follows:
(1) adopts Δ T-contribution margin that the temperature of each heat exchange stream thigh is revised, draw the warm enthalpy diagram that becomes the folder point, obtain revised folder point temperature and energy object.
(2) on warm enthalpy diagram, construct the weighing apparatus composite curve chart that flattens, cut apart the enthalpy interval, the superstructure model of structure heat exchanger network.
(3) be objective function with the total cost minimum, under the prerequisite that satisfies the constraint of heat balance constraint condition and temperature difference correction constraint and heat interchanger minimum area, find the solution superstructure.Especially, at the interval k of a certain given sub-enthalpy, if I bar hot-fluid thigh and J bar cold flow thigh are arranged, so in this sub-enthalpy interval, it is total that (the individual coupling number of I * J), and the operative constraint equation number is (I+J-1) is when I * J＞I+J-1, need carry out inside son optimization during Load Distribution in determining this sub-enthalpy interval, the investment cost in the sub-for this reason enthalpy of its optimization aim interval minimizes.
When (4) satisfying the condition of convergence, finish to calculate, export best heat exchanger network structure.
Heat interchanger expense formula is: 1300 * A
0.6$/yr; Hot public work expense: 80$/kWyr, cold public work expense: 20$/kWyr.
Table 1 instance data
In order to verify effect of the present invention, and compare with existing existing method both at home and abroad, logistics data such as table 1 that present embodiment adopts, exist in this problem 3 strands of hot-fluids (H1, H2, H3), 2 strands of cold flows (C1, C2), hot public H4, cold public C3.
The result of calculation of this method is as shown in table 2, heat exchanger network structure such as Fig. 2 of acquisition.Case study shows, under equal conditions, adopts the heat exchanger network structure that the inventive method obtained, and has littler heat interchanging area and annual total costs.And, have better practicability and systematicness because the thought that this method combines with mathematics based on thermodynamics adopts the multi-level optimization strategy.
Table 2 result of calculation
The foregoing description is used for the present invention that explains, rather than limits the invention, and in the protection domain of spirit of the present invention and claim, any modification and change to the present invention makes all fall into protection scope of the present invention.
Claims (3)
1. heat exchange network optimization synthesis method based on the correction of Δ T-contribution margin, it is characterized in that: described heat exchange network optimization synthesis method may further comprise the steps:
1), calculate heat exchanging fluid heat transfer temperature difference contribution margin Δ T, and contribution margin Δ T is revised, innovation representation is:
ΔT
H，i+ΔT
C，j≤T
H，i-T
C，j (1)
In the following formula (1): Δ T
H, iBe the heat transfer temperature difference contribution margin of hot fluid i, unit is ℃; Δ T
C, jBe the heat transfer temperature difference contribution margin of cold fluid j, unit is ℃; T
H, iBe the inlet temperature of hot fluid i, unit is ℃; T
C, jBe the inlet temperature of cold fluid j, unit is ℃;
2), based on thermodynamic analysis and starting point, adopt the method for mathematical modeling carry out the energy object of heat exchanger network, area target, cost objective and network structure design synchronously not comprehensively, detailed process is as follows:
(2.1) adopt revised Δ T-contribution margin that the temperature of each heat exchange stream thigh is revised, draw the warm enthalpy diagram that becomes the folder point, obtain revised folder point temperature and energy object;
(2.2) on warm enthalpy diagram, construct the weighing apparatus composite curve chart that flattens, cut apart the enthalpy interval, the superstructure model of structure heat exchanger network;
(2.3) be objective function with the total cost minimum, find the solution the superstructure optimization problem, obtain best heat exchanger network; The expression formula of objective function is:
TAC＝(Q
Hmin·C
H)+(Q
Cmin·C
C)+C
CAP (2)
In the following formula (2): TAC is the total annual charge in heat exchange networking, and unit is unit/year; Q
HminAnd Q
CminBe respectively the minimum thermal of trying to achieve in the described step (2.1), cold public work consumption, unit is J/; C
HAnd C
CBe respectively the heat of using in the heat exchanger network, the unit costs of cold public work, unit is unit/J; C
CAPBe the total annualized capital cost usefulness of network, unit is unit/year.
2. a kind of heat exchange network optimization synthesis method as claimed in claim 1 based on the correction of Δ T-contribution margin, it is characterized in that: in described step (2.3), when finding the solution the superstructure model, regulation MINIMUM HEAT-TRANSFER AREA constraint condition, make and reject little load, small size heat interchanger automatically in the solution procedure of heat exchanger network superstructure, realize the lax of heat exchanger network.
3. improved heat exchanger network Δ T-contribution margin as claimed in claim 1 or 2 is optimized synthetic method, it is characterized in that: in described step (2.3), at the interval k of a certain given sub-enthalpy, if I bar hot-fluid thigh and J bar cold flow thigh are arranged, so in this sub-enthalpy interval, total (the individual coupling number of I * J), and the operative constraint equation number is (I+J-1), when I * J＞I+J-1, need carry out inside son optimization during Load Distribution in determining this sub-enthalpy interval, the investment cost in the sub-for this reason enthalpy of its optimization aim interval minimizes.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200910155061A CN101777083A (en) | 2009-12-15 | 2009-12-15 | Heat exchange network optimization synthesis method based on Delta T-contribution value correction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200910155061A CN101777083A (en) | 2009-12-15 | 2009-12-15 | Heat exchange network optimization synthesis method based on Delta T-contribution value correction |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN101777083A true CN101777083A (en) | 2010-07-14 |
Family
ID=42513544
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN200910155061A Pending CN101777083A (en) | 2009-12-15 | 2009-12-15 | Heat exchange network optimization synthesis method based on Delta T-contribution value correction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101777083A (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102155860A (en) * | 2010-12-28 | 2011-08-17 | 浙江工业大学 | Method for constructing heat exchange network based on exergy consumption cost |
| CN102446299A (en) * | 2011-07-19 | 2012-05-09 | 北京三博中自科技有限公司 | Heat exchanger network analysis method for process industry |
| CN103914604A (en) * | 2012-12-31 | 2014-07-09 | 北京宜能高科科技有限公司 | Heat exchanger network transformation optimization method taking certain point temperature in maximized network as target |
| CN104793495A (en) * | 2015-04-15 | 2015-07-22 | 浙江工业大学 | Method for determining maximum heat recycling potential of heat exchange networks with non-isothermal phase-change fluid |
| CN104914719A (en) * | 2015-04-15 | 2015-09-16 | 浙江工业大学 | Method for determining minimum cooling common engineering requirement of heat exchange network containing non-isothermal phase-change fluid |
| CN107665280A (en) * | 2017-09-27 | 2018-02-06 | 浙江工业大学 | A kind of Retrofit of Heat Exchanger Networks optimization method based on performance simulation |
| CN107918280A (en) * | 2017-11-21 | 2018-04-17 | 华东理工大学 | A kind of NETWORK IN REFINERY Optimization Scheduling for pressing from both sides point method and being mixed with superstructure method |
| CN109708516A (en) * | 2018-11-14 | 2019-05-03 | 浙江工业大学 | A kind of old heat exchange equipment Efficient utilization method for improving industrial heat exchange network |
| CN112035962A (en) * | 2020-09-09 | 2020-12-04 | 中国航空工业集团公司沈阳飞机设计研究所 | Superstructure model-based optimization method for airborne heat exchange network |
-
2009
- 2009-12-15 CN CN200910155061A patent/CN101777083A/en active Pending
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102155860A (en) * | 2010-12-28 | 2011-08-17 | 浙江工业大学 | Method for constructing heat exchange network based on exergy consumption cost |
| CN102155860B (en) * | 2010-12-28 | 2012-11-14 | 浙江工业大学 | Method for constructing heat exchange network based on exergy consumption cost |
| CN102446299A (en) * | 2011-07-19 | 2012-05-09 | 北京三博中自科技有限公司 | Heat exchanger network analysis method for process industry |
| CN102446299B (en) * | 2011-07-19 | 2014-06-25 | 北京三博中自科技有限公司 | Heat exchanger network analysis method for process industry |
| CN103914604A (en) * | 2012-12-31 | 2014-07-09 | 北京宜能高科科技有限公司 | Heat exchanger network transformation optimization method taking certain point temperature in maximized network as target |
| CN103914604B (en) * | 2012-12-31 | 2017-02-08 | 北京宜能高科科技有限公司 | Heat exchanger network transformation optimization method taking certain point temperature in maximized network as target |
| CN104914719A (en) * | 2015-04-15 | 2015-09-16 | 浙江工业大学 | Method for determining minimum cooling common engineering requirement of heat exchange network containing non-isothermal phase-change fluid |
| CN104793495A (en) * | 2015-04-15 | 2015-07-22 | 浙江工业大学 | Method for determining maximum heat recycling potential of heat exchange networks with non-isothermal phase-change fluid |
| CN107665280A (en) * | 2017-09-27 | 2018-02-06 | 浙江工业大学 | A kind of Retrofit of Heat Exchanger Networks optimization method based on performance simulation |
| CN107918280A (en) * | 2017-11-21 | 2018-04-17 | 华东理工大学 | A kind of NETWORK IN REFINERY Optimization Scheduling for pressing from both sides point method and being mixed with superstructure method |
| CN107918280B (en) * | 2017-11-21 | 2020-09-04 | 华东理工大学 | Oil refinery hydrogen network optimization scheduling method with mixed pinch method and superstructure method |
| CN109708516A (en) * | 2018-11-14 | 2019-05-03 | 浙江工业大学 | A kind of old heat exchange equipment Efficient utilization method for improving industrial heat exchange network |
| CN112035962A (en) * | 2020-09-09 | 2020-12-04 | 中国航空工业集团公司沈阳飞机设计研究所 | Superstructure model-based optimization method for airborne heat exchange network |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101777083A (en) | Heat exchange network optimization synthesis method based on Delta T-contribution value correction | |
| CN106815661B (en) | Decomposition coordination scheduling method of combined heat and power system | |
| Varbanov et al. | Total Site targeting with process specific minimum temperature difference (ΔTmin) | |
| CN103558046B (en) | A kind of heat exchanger energy efficiency evaluation system | |
| CN102155860B (en) | Method for constructing heat exchange network based on exergy consumption cost | |
| CN104866702B (en) | A kind of method of utilization pivot decoupling computation hot-summer and cold-winter area building load | |
| Liu et al. | Simultaneous integrated design for heat exchanger network and cooling water system | |
| MY163918A (en) | A plant for regasification of lng | |
| CN203940504U (en) | The suction-type lithium bromide heat-exchange system of two path water heat supply simultaneously | |
| CN107665280A (en) | A kind of Retrofit of Heat Exchanger Networks optimization method based on performance simulation | |
| Pan et al. | Multi-plant indirect heat integration based on the Alopex-based evolutionary algorithm | |
| CN105576943B (en) | A kind of high voltage direct current converter valve heat recovery system and its implementation | |
| CN205037807U (en) | Utilize steam heat transfer system of steam ejector exhaust steam recovery | |
| Xia et al. | A novel method based on entransy theory for setting energy targets of heat exchanger network | |
| Rév et al. | Hidden and pseudo pinch phenomena and relaxation in the synthesis of heat-exchange networks | |
| CN212941540U (en) | A system for nitrotoluene melt crystallization | |
| CN105135396B (en) | The system and its recovery method of a kind of tetrafluoroethene cracking gas energy recovery | |
| CN207299138U (en) | A kind of big temperature difference heat exchange station of cascade utilization | |
| CN111402074B (en) | Comprehensive optimization method for mass energy of circulating water system | |
| CN204529906U (en) | A kind of modified version insulation ageing device | |
| Chena et al. | Large-scale MILP transshipment models for heat exchanger network synthesis | |
| CN204897408U (en) | Waste heat utilization system in rare nitric acid manufacture | |
| PIRES et al. | An energy integration tool for batch process | |
| Tian et al. | Multi-plant heat integration with shell-and-tube heat exchangers for multi-period operations | |
| Zhu et al. | Simulation and optimization of cooling water systems with closed wet cooling tower based on reduced order model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C12 | Rejection of a patent application after its publication | ||
| RJ01 | Rejection of invention patent application after publication |
Open date: 20100714 |