項(xiàng)順伯，黃 燕

（1. 廣東省石化裝備故障診斷重點(diǎn)實(shí)驗(yàn)室，廣東石油化工學(xué)院，廣東 茂名 525000；2. 廣東石油化工學(xué)院 化學(xué)工程學(xué)院， 廣東 茂名 525000）

工藝與裝備

無(wú)分流的管殼式換熱網(wǎng)絡(luò)優(yōu)化研究

項(xiàng)順伯1，黃 燕2

（1. 廣東省石化裝備故障診斷重點(diǎn)實(shí)驗(yàn)室，廣東石油化工學(xué)院，廣東 茂名 525000；2. 廣東石油化工學(xué)院 化學(xué)工程學(xué)院， 廣東 茂名 525000）

提出無(wú)分流管殼式換熱網(wǎng)絡(luò)的分級(jí)超結(jié)構(gòu)模型，得出優(yōu)化該換熱網(wǎng)絡(luò)的目標(biāo)函數(shù)和約束條件，采用位置加權(quán)的改進(jìn)粒子群算法進(jìn)行優(yōu)化舉例。結(jié)果表明，通過(guò)選擇合適的加權(quán)因子，位置加權(quán)的改進(jìn)粒子群算法能加快算法的收斂速度，更快地得到全局最優(yōu)解，從而得到最佳的管殼式換熱網(wǎng)絡(luò)優(yōu)化參數(shù)。

無(wú)分流; 管殼式換熱網(wǎng)絡(luò)；加權(quán)因子；改進(jìn)粒子群算法

換熱網(wǎng)絡(luò)是工業(yè)過(guò)程中能量回收和綜合利用的一種設(shè)備，換熱網(wǎng)絡(luò)很多，管殼式換熱網(wǎng)絡(luò)使用的范圍最為廣泛。由于管殼式換熱網(wǎng)絡(luò)的復(fù)雜結(jié)構(gòu)以及諸多變量間的非線性關(guān)系，其數(shù)學(xué)模型并不能準(zhǔn)確地建立，長(zhǎng)期以來(lái)，其優(yōu)化一直是一個(gè)研究的問(wèn)題。無(wú)分流的管殼式換熱網(wǎng)絡(luò)是指冷熱股流體在流體過(guò)程中除換熱外不存在分支流動(dòng)的流體，其優(yōu)化方法很多，傳統(tǒng)的方法是夾點(diǎn)分析法[1,2]和數(shù)學(xué)規(guī)劃法[3]。隨著人工智能技術(shù)的發(fā)展與應(yīng)用，其也應(yīng)用于換熱網(wǎng)絡(luò)的優(yōu)化中，如遺傳算法[4]、進(jìn)化算法[5]、粒子群算法[6]等，由于這類(lèi)算法自身的局限性，在優(yōu)化的搜索范圍內(nèi)，能比較容易地找到較優(yōu)解，但不能被證明是最優(yōu)解，因此，人工智能算法也在發(fā)展優(yōu)化中，如改進(jìn)的粒子群算法[8]。本文以管殼式換熱網(wǎng)絡(luò)綜合費(fèi)用最小為目標(biāo)，以改進(jìn)的粒子群算法進(jìn)行優(yōu)化舉例，在全局范圍內(nèi)，搜索出最優(yōu)的解，從而得到更優(yōu)的管殼式換熱網(wǎng)絡(luò)結(jié)構(gòu)。

1 無(wú)分流管殼式換熱網(wǎng)絡(luò)優(yōu)化模型

1.1 無(wú)分流管殼式換熱網(wǎng)絡(luò)分級(jí)超結(jié)構(gòu)模型

本文無(wú)分流的管殼式換熱網(wǎng)絡(luò)模型采用分級(jí)超結(jié)構(gòu)模型[7]，熱流體2股、冷流體3股，所有冷、熱流體之間的換熱組合成一個(gè)級(jí)，共分3級(jí)，其結(jié)構(gòu)圖如圖1所示。NH=2為熱流體數(shù)目，NC=3為冷流體數(shù)目，則圖1中的2股熱流體和3股冷流體之間的換熱組合共有NH·NC=6種，每種冷熱流體的組合需要一個(gè)換熱器，圖1所示的每一級(jí)換熱組合需要6個(gè)換熱器。圖中Hi(i= 1，2，…，NH) 表示熱流體編號(hào)，Cj(j= 1，2，…，NC) 表示冷流體編號(hào)，水平線及箭頭表示流體的流動(dòng)方向，冷、熱流體之間的換熱用2 個(gè)“○”以及它們之間的連線表示。換熱網(wǎng)絡(luò)的級(jí)數(shù)為NK，這里取NK=max(NH, NC)，這樣換熱網(wǎng)絡(luò)最大的換熱器數(shù)目為：熱、冷公用工程換熱器分別加在第KN級(jí)的后面和第1級(jí)的前面。

圖1 無(wú)分流管殼式換熱網(wǎng)絡(luò)分級(jí)超結(jié)構(gòu)圖Fig.1 Stage-wise superstructure figure of shell-and-tube heat exchanger network without stream splits

1.2 目標(biāo)函數(shù)

以管殼式換熱網(wǎng)絡(luò)年綜合總費(fèi)用ATC最小為優(yōu)化目標(biāo)函數(shù)優(yōu)化換熱網(wǎng)絡(luò)結(jié)構(gòu)，年綜合總費(fèi)用由設(shè)備投資總費(fèi)用和運(yùn)行總費(fèi)用兩部分組成，則最小年綜合總費(fèi)用如(1)式：

式中：前2項(xiàng)為運(yùn)行費(fèi)用，CCU為冷公用工程費(fèi)用系數(shù)，QCU,i為第i股熱流體與冷公用工程之間的換熱量；CHU為熱公用工程費(fèi)用系數(shù)，QHU,j第j股冷流體與熱公用工程之間的換熱量；(1)式的后三項(xiàng)為設(shè)備投資費(fèi)用，CC為單個(gè)冷卻器固定投資費(fèi)用，C0為冷卻器面積費(fèi)用系數(shù)，ACU,i為冷卻第i股熱流體的冷卻器換熱面積；CH為單個(gè)加熱器固定投資費(fèi)用，C1為加熱器面積費(fèi)用系數(shù)，AHU,j為加熱第j股冷流體的加熱器換熱面積；CE為單個(gè)換熱器固定投資費(fèi)用，C3為換熱器面積費(fèi)用系數(shù)，Aijk是第i股熱流體和第j股冷流體同時(shí)所經(jīng)過(guò)的第k個(gè)換熱器的表面積。

①約束條件[9]：

以Ti,in和Ti,out分別表示第i股熱流體的入口溫度和出口溫度，Tj,in和Tj,out分別表示第j股冷流體的入口溫度和出口溫度，Gi、Gj分別第i股熱流體和第j股冷流體的熱容流率，則

式中：ijkQ表示第i股熱流體與第j股冷流體在第k級(jí)交換熱量時(shí)換熱器的換熱量。

每級(jí)換熱器的熱平衡：

Ti,k,in和Ti,k,out分別表示第i股熱流體經(jīng)過(guò)第k級(jí)換熱器時(shí)的入口溫度和出口溫度；Tj,k,out和Tj,k,in分別表示第j股冷流體經(jīng)過(guò)第k級(jí)換熱器時(shí)的出口溫度和入口溫度。

②最小換熱溫差約束

式中dTijk、dTCU,i和dTHU,j分別表示分別表示每個(gè)換熱器、冷卻器和加熱器的換熱溫差，ΔTmin為最小允許傳熱溫差。

冷卻器換熱面積、加熱器換熱面積和換熱器的表面積分別由下式計(jì)算：

KCU,i、KHU,j和Kij分別為傳熱系數(shù)，TDCU,i、TDHU,j和TDijk分別是冷卻器、加熱器和換熱器的逆流對(duì)數(shù)平均溫差，由以下三式計(jì)算[11]：

③非負(fù)性約束

2 位置加權(quán)的改進(jìn)粒子群算法[8]

標(biāo)準(zhǔn)粒子群算法，簡(jiǎn)稱(chēng)PSO算法，算法源自鳥(niǎo)類(lèi)捕食行為的模擬，PSO算法具體表達(dá)：設(shè)在維搜索空間中，m個(gè)粒子構(gòu)成種群其中第i個(gè)粒子位置，速度vi={vi1, vi2,...,vin}T；粒子i經(jīng)過(guò)的最優(yōu)位置記為pi={pi1,pi2,...,pin}T，整個(gè)粒子群經(jīng)過(guò)的最優(yōu)位置記為

pg={pg1,pg2,...,pgn}T，粒子xi將按式（17）、式（18）改變速度和位置。

式中，ω為慣性權(quán)重，表示前一時(shí)刻的速度對(duì)當(dāng)前時(shí)刻速度的影響。ω較小時(shí)，粒子容易陷入局部最優(yōu)，當(dāng)ω較大時(shí)，算法收斂速度較慢。1c為粒子自身的加速度權(quán)重系數(shù)，代表粒子自身經(jīng)驗(yàn)的學(xué)習(xí)，2c為全局加速度權(quán)重系數(shù)，代表粒子向群體經(jīng)驗(yàn)的學(xué)習(xí)；1r、2r為介于[0,1]間的相互獨(dú)立的隨機(jī)數(shù)?？紤]粒子前一次位置對(duì)更新之后位置的影響，引入影響因子α，稱(chēng)之加權(quán)因子，對(duì)基本粒子群算法公式進(jìn)行更新，即在公式(18)后面加上

來(lái)進(jìn)行粒子位置的更新，把構(gòu)造出的改進(jìn)粒子群算法稱(chēng)為位置加權(quán)的改進(jìn)粒子群算法。式(18)和式(19)合并，得到

式（20）表明，對(duì)位置的加權(quán)因子α，相當(dāng)于對(duì)速度的約束因子(1-α)。這相當(dāng)于在搜索空間里面對(duì)搜索的方向未加限制，而對(duì)搜索的步長(zhǎng)加以約束。這樣能夠有效地避免在搜索過(guò)程中由于步長(zhǎng)過(guò)大而錯(cuò)過(guò)全局最優(yōu)解，減少搜索過(guò)程中的盲目性，能提高算法的搜索效率。

根據(jù)上述分析，位置加權(quán)的改進(jìn)粒子群算法的步驟如下：

（1）設(shè)置算法的參數(shù)和最大迭代次數(shù)；

（2）在滿(mǎn)足控制變量約束條件下隨機(jī)賦予粒子群中每個(gè)粒子初始速度和初始位置；

（3）對(duì)群體中的每個(gè)粒子，計(jì)算其適應(yīng)度值，找到每個(gè)粒子經(jīng)歷的最好位置和群體經(jīng)歷的最好位置；

（4）按式（20）對(duì)粒子位置進(jìn)行加權(quán)處理；

（5）利用式（17）、（18）和（20）更新粒子速度和位置；

（6）判斷運(yùn)算結(jié)束條件，若不滿(mǎn)足結(jié)束條件，則返回執(zhí)行步驟（3）；若滿(mǎn)足結(jié)束條件，則執(zhí)行步驟（7）；

（7）輸出最優(yōu)位置。

3 無(wú)分流管殼式換熱網(wǎng)絡(luò)優(yōu)化實(shí)例

本文的優(yōu)化實(shí)例源自文獻(xiàn)[10],該管殼式換熱網(wǎng)絡(luò)由6股熱流體、4股冷流體組成，換熱器器面積費(fèi)用系數(shù)為60A$/m2，熱公用工程費(fèi)用系數(shù)100$·kW-1·a-1，冷公用工程費(fèi)用系數(shù)為15$·kW-1·a-1，流體的初始參數(shù)如表1所示。

表1 優(yōu)化舉例流體數(shù)據(jù)表Table 1 Steam data table of optimization illustration

設(shè)minTΔ=20 K，采用位置加權(quán)的改進(jìn)粒子群算法進(jìn)行優(yōu)化，優(yōu)化中，最大迭代次數(shù)為2 000，粒子種群數(shù)為24，c1、c2均取為2，α取0.6，優(yōu)化結(jié)果如表2所示。

表2 舉例優(yōu)化結(jié)果Table 2 Result of optimization illustration

采用位置加權(quán)的改進(jìn)粒子群算法時(shí)，通過(guò)選擇合適的加權(quán)因子，能加快算法的收斂，更快地得到最優(yōu)解。本文的優(yōu)化結(jié)果與文獻(xiàn)[9]、文獻(xiàn)[10]進(jìn)行對(duì)比，通過(guò)比較發(fā)現(xiàn)，本文的優(yōu)化所得到的年綜合費(fèi)用最低，熱交換總面積也最小，這說(shuō)明本文采用的方法更能得到最優(yōu)的換熱網(wǎng)絡(luò)結(jié)構(gòu)。

4 結(jié) 論

本文提出了無(wú)分流管殼式換熱網(wǎng)絡(luò)的級(jí)超結(jié)構(gòu)模型，得出優(yōu)化該換熱網(wǎng)絡(luò)的目標(biāo)函數(shù)以及線性的約束條件，然后采用位置加權(quán)的改進(jìn)粒子群算法進(jìn)行優(yōu)化舉例，優(yōu)化結(jié)果表明本文的方法能更快更準(zhǔn)地得到無(wú)分流管殼式換熱網(wǎng)絡(luò)的相關(guān)參數(shù)，如其他方法相比，進(jìn)一步說(shuō)明本文的方法是有效的。

［1］Linnhoff B,Flower J R .Synthesis of Heat Exchanger Networks[J]. AIChE.,1978,24(4):633-642.

［2］Linnhoff B. The pinch design method for heat exchanger networ ks[J]. Chem Eng Sci, 1983,38(5).

［3］TGundersen, SDuvold, AHashem-iAhmady.An extended verticalMI LP modle for heat exchanger network synthesis[J]. Computers C hem Eng, 1996, 20:97-102.

［4］ALLEN B，SAVARD-GOGUEN M，GOSSELIN L． Optimizing heat exchanger networks with genetic algorithms for designing e ach heat exchanger including condensers[J]．Applied Thermal En gineering，2009，29( 16) : 3437-3444．

［5］YERRAMSETTY K M，MURTY C V S．Synthesis of cost-opti mal heat exchanger networks using differential evolution[J]．Com puters and Chemical Engineering，2008，32( 8) : 1861-1876．

［6］RAVAGNANI M, SILVA A P, BISCAIA J E C, et al. Optimal design of shell-and-tube heat ex changers using particle swarm optimization[J]. Industrial &Engineering Chemistry Research, 20 09, 48(6): 2927-2935.

［7］T.F.Yee,1.E.Grossmann. Simultaneous optimization model for he at integration-II: Heat exchanger network synthesis[J].Compute rs & Chemical Engineering,1990,14(10):1165-1184.

［8］朱童,李小凡,魯明文.位置加權(quán)的改進(jìn)粒子群算法[J].計(jì)算機(jī)工程與應(yīng)用, 2011,47(5):4-6.

［9］趙亮,霍兆義,李棟斌,等. 無(wú)分流換熱網(wǎng)絡(luò)同步優(yōu)化綜合[J].熱科學(xué)與技術(shù), 2013, 12(1):69-70.

［10］K HORASANYR M, FESANGHARY M. A novel approach for synthesis of cost-optimal heat exchanger networks[J]. Comput Chem Eng, 2009,33(8):1363-1370.

［11］T.F.Yeeetal.Simultaneous optimization models for heat integratio n-I: Area and energy targeting and modeling of multi-stream exchangers[J].Computers & Chemical Engineering, 1990, 14(10): 1151-1164.

Optimization Study of Shell-and-tube Heat Exchanger Network Without Stream Splits

XIANG Shun-bo1，HUANG Yan2
（1. Guangdong Provincial Key Laboratory of Petrochemical Equipment Fault Diagnosis，Guangdong University of Petrochemical Technology，Guangdong Maoming 525000；2. College of Chemical Engineering, Guangdong University of Petrochemical Technology, Guangdong Maoming 525000，China）

The stage-wise superstructure model of shell-and-tube heat exchanger network without stream spits was proposed, objective function and constraint conditions of the heat exchanger network were obtained, an optimization example was given using improved particle swarm optimization algorithm with position weighted. The result shows that, through choosing appropriate weighting factor, improved particle swarm optimization algorithm with position weighted can accelerate the algorithm’s convergence velocity, and can obtain global optimum solution, at last the best optimization parameters of hell-and-tube heat exchanger network can be obtained.

No stream spits; Shell-and-tube heat exchanger network; Weighting factor; Advanced PSO algorithm

TQ 018

A

1671-0460（2015）01-0145-03

廣東高校石化過(guò)程裝備故障診斷與信息化控制工程中心開(kāi)放基金項(xiàng)目( 512028).

2014-08-13

項(xiàng)順伯（1979-），男，安徽樅陽(yáng)人，講師，碩士，研究方向：從事人工智能方面的技術(shù)工作。E-mail：qingcheng33@163.com。