萬義群，崔國民，方大俊，彭富裕

（上海理工大學(xué)新能源科學(xué)與工程研究所，上海 200093）

換熱網(wǎng)絡(luò)的綜合及優(yōu)化是實(shí)現(xiàn)能量綜合利用的重要環(huán)節(jié)，其在節(jié)能和降低產(chǎn)品成本方面有著非常重要的意義。迄今國內(nèi)外已開發(fā)出許多針對換熱網(wǎng)絡(luò)的優(yōu)化綜合方法，這些方法基本可以分為夾點(diǎn)技術(shù)法[1-2]、數(shù)學(xué)規(guī)劃法[3]和啟發(fā)式方法[4-6]三大類。

換熱網(wǎng)絡(luò)的優(yōu)化過程中常用的目標(biāo)函數(shù)是換熱網(wǎng)絡(luò)的年綜合費(fèi)用，但從優(yōu)化的結(jié)果來看，換熱網(wǎng)絡(luò)的綜合費(fèi)用并不是完全與換熱效能相一致的，這是由于以綜合費(fèi)用為目標(biāo)的優(yōu)化方法并沒有從提高總體換熱效能的角度出發(fā)，獲得的換熱網(wǎng)絡(luò)只是在綜合費(fèi)用上達(dá)到一定程度的最優(yōu)。從本質(zhì)上講，提高換熱網(wǎng)絡(luò)的換熱效能是解決換熱網(wǎng)絡(luò)優(yōu)化問題的最根本途徑，因為換熱效能的提高能夠使得能量得到有效地綜合利用，如果要獲得更優(yōu)的換熱網(wǎng)絡(luò)結(jié)構(gòu)，必須從換熱效能的角度為換熱網(wǎng)絡(luò)提出新的評價指標(biāo)及新的優(yōu)化思路，從根本上改進(jìn)換熱網(wǎng)絡(luò)結(jié)構(gòu)。

過增元院士等[7-8]從換熱器沿程溫差與換熱效能的關(guān)系出發(fā)，提出了兩股流換熱器的溫差均勻性優(yōu)化原則，依據(jù)該原則對換熱器進(jìn)行改進(jìn)能夠提高換熱器的換熱效能；文獻(xiàn)[9]將溫差均勻性優(yōu)化原則應(yīng)用在多股流換熱器上，建立了多股流換熱器設(shè)計的溫差均勻性因子，能夠用于衡量多股流換熱器通道排列的優(yōu)劣，從而能指導(dǎo)多股流換熱器的優(yōu)化設(shè)計；文獻(xiàn)[10-12]從換熱網(wǎng)絡(luò)整體強(qiáng)化的角度出發(fā)提出了換熱網(wǎng)絡(luò)的溫差均勻性優(yōu)化原則，建立了換熱網(wǎng)絡(luò)的溫差均勻性因子，試圖為換熱網(wǎng)絡(luò)的換熱效能提供一種新的評價指標(biāo)，并能指導(dǎo)換熱網(wǎng)絡(luò)的優(yōu)化設(shè)計。

從場協(xié)同理論[13-14]的角度來看，將換熱器沿程的對流換熱系數(shù)以及換熱器的固壁導(dǎo)熱系數(shù)、固壁均勻程度都考慮進(jìn)來，那么換熱器的熱流邊界條件越均勻則其換熱效能就越高[15-16]；而換熱網(wǎng)絡(luò)是一個多流體，多物性的復(fù)雜換熱流場，因此換熱網(wǎng)絡(luò)的溫差均勻性因子并不能作為換熱效能絕對的評價指標(biāo)。如果將各流體的換熱系數(shù)都考慮進(jìn)來從而建立換熱網(wǎng)絡(luò)的熱流均勻性因子，則其將更適合于作為換熱網(wǎng)絡(luò)效能的評價指標(biāo)。

鑒于此，本文在換熱器網(wǎng)絡(luò)的溫差均勻性因子及場協(xié)同理論的基礎(chǔ)上，提出并建立換熱網(wǎng)絡(luò)的熱流均勻性因子，著重驗證及分析換熱網(wǎng)絡(luò)整體換熱量與熱流均勻性因子之間的關(guān)系，試圖從能量有效利用的角度為換熱網(wǎng)絡(luò)的換熱效能建立新的評價指標(biāo)，并為換熱器網(wǎng)絡(luò)的綜合及優(yōu)化提供新的方法和思路。

1 換熱網(wǎng)絡(luò)熱流均勻性因子的建立

1.1 換熱網(wǎng)絡(luò)的溫差均勻性因子

文獻(xiàn)[11]從換熱網(wǎng)絡(luò)的換熱效能出發(fā)，提出了換熱網(wǎng)絡(luò)溫差均勻性因子的概念，并通過對換熱網(wǎng)絡(luò)溫差場的分析，建立了基于換熱網(wǎng)絡(luò)離散溫差場的溫差均勻性因子，其建立的溫差均勻性因子表達(dá)式如式（1）。

式中，n為換熱器的個數(shù)；ΔTi為第i個換熱器的對數(shù)平均溫差。

通過具體的實(shí)例分析表明[11]：換熱網(wǎng)絡(luò)溫差均勻性因子φT越小，則整個換熱網(wǎng)絡(luò)的溫差場越均勻，換熱效能越高。但該表達(dá)式?jīng)]有將換熱網(wǎng)絡(luò)中各流體的換熱系數(shù)考慮進(jìn)來，在實(shí)際工程中所面對的換熱網(wǎng)絡(luò)優(yōu)化問題總是包含不同換熱系數(shù)的流體，對于多流體，多物性的換熱網(wǎng)絡(luò)復(fù)雜流場，溫差均勻性因子的適用性受到了很大的限制，因此本研究將建立更適合于評價換熱效能的換熱網(wǎng)絡(luò)的熱流均勻性因子。

1.2 換熱網(wǎng)絡(luò)的熱流均勻性因子

根據(jù)對流換熱過程的場協(xié)同理論，恒壁溫的邊界條件將抑制對流換熱，而恒熱流的換熱邊界條件，將最大限度地發(fā)揮對流換熱的作用，據(jù)此建立了換熱器熱流均勻性因子的概念[15]，將換熱器沿軸向劃分為N個微元換熱體，每個微元的熱流密度為qi(i=1，2，··，N)，則單體換熱器的熱流均勻性因子φq,s可表示為式（2）。

式中，N為單個換熱器的微元體個數(shù)；qi為第i個微元體的熱流密度。

本研究借鑒單體換熱器的熱流均勻性因子表達(dá)式，將復(fù)雜換熱網(wǎng)絡(luò)中每個換熱器作為一個換熱單元體劃分為n個單元，建立換熱網(wǎng)絡(luò)的熱流均勻性因子φq,net，表達(dá)式如式（3）、式（4）。

綜上所述，胎心監(jiān)護(hù)、臍動脈血流聯(lián)合血清HIF-1α水平檢測對胎兒窘迫的診斷有重要的臨床價值，其指標(biāo)表達(dá)異常可直接關(guān)乎新生兒的結(jié)局，臨床對于臨產(chǎn)的孕婦應(yīng)建議進(jìn)行聯(lián)合檢測，盡早發(fā)現(xiàn)是否出現(xiàn)胎兒窘迫，積極采取有效措施，將新生兒傷害風(fēng)險降至最低。

式中，n為換熱器的個數(shù)；Qi、Ai、qi分別為第i個換熱器的換熱量、換熱面積和熱流密度。

考慮到換熱網(wǎng)絡(luò)中單個換熱器的換熱量均不相同，其對熱流均勻性因子的影響及貢獻(xiàn)不同，如果所有的換熱器都按相同的權(quán)重計算，那么就很難準(zhǔn)確反映出這一網(wǎng)絡(luò)的換熱特點(diǎn)，因此通過加權(quán)的方式增加換熱能力強(qiáng)的換熱器在熱流均勻性因子中的作用。計算各換熱器的換熱量Qi(i=1，2，··，n)，取其中的最小值Qmin=min｛Q1，Q2，··，Qn｝，將各換熱器的換熱量與Qmin的比值mi作為各換熱器的權(quán)重值，則熱流均勻性因子按換熱量加權(quán)的計算公式為式（5）、式（6）。

式中，n為換熱器的個數(shù)；Qi、mi、qi分別為第i個換熱器的換熱量、權(quán)重值和熱流密度。

本研究最終建立的換熱網(wǎng)絡(luò)的熱流均勻性因子φ'q,net表征的是一個換熱網(wǎng)絡(luò)內(nèi)部熱流均勻性的程度。從式（5）可以看出，熱流均勻性因子越接近1則換熱網(wǎng)絡(luò)內(nèi)部熱流密度越均勻，其換熱性能也越好；該因子能否和換熱網(wǎng)絡(luò)的總換熱量之間形成對應(yīng)關(guān)系，能否作為換熱網(wǎng)絡(luò)換熱效能的評價指標(biāo)，需要通過具體的實(shí)例進(jìn)行驗證及分析。

2 熱流均勻性因子的實(shí)例驗證及分析

2.1 算例一

算例一采用15股流換熱器網(wǎng)絡(luò)[17-20]來檢驗并分析換熱網(wǎng)絡(luò)中的熱流均勻性因子與換熱性能的關(guān)系，該換熱網(wǎng)絡(luò)由8股熱流體和7股冷流體組成，流股的初始參數(shù)見表1，文獻(xiàn)[19]中優(yōu)化的結(jié)構(gòu)如圖1所示，圖中用兩個“●”及其相連的豎直短線表示換熱器單元，單個“●”表示公用工程換熱單元。本研究將在保持換熱單元總面積相同的條件下與在改變換熱網(wǎng)絡(luò)結(jié)構(gòu)的條件下驗證分析換熱網(wǎng)絡(luò)的換熱性能和熱流均勻性因子之間的關(guān)系，其中換熱性能以換熱網(wǎng)絡(luò)總體換熱量的大小作為衡量指標(biāo)（不包括冷熱公用工程的熱負(fù)荷）。

表1 算例1的流體參數(shù)

2.1.1 保持換熱單元總面積相同條件下的熱流均勻性因子

采用圖1所示的換熱網(wǎng)絡(luò)結(jié)構(gòu)，隨機(jī)改變各換熱器單元的面積，計算30組熱流均勻性因子，規(guī)定保持換熱面積的總和等于3000 m2，以此來考察整個換熱網(wǎng)絡(luò)的換熱量與熱流均勻性因子之間的關(guān)系。將換熱器面積的總和保持不變則可以削弱換熱器面積大小對整個網(wǎng)絡(luò)換熱性能的影響，在這種情況下考察的換熱網(wǎng)絡(luò)熱流均勻性因子更具有參考價值。

熱流體數(shù)目為NH=8，冷流體數(shù)目為NC=7，換熱網(wǎng)絡(luò)的級數(shù)為NK=2，隨機(jī)獲取換熱器單元面積并進(jìn)行模擬計算，必須滿足如式（7）～式（9）的約束條件。

式中，TH,i,out、TC,j,out分別為熱流體和冷流體的出口溫度；TH,i、TC,j分別為熱流體和冷流體的目標(biāo)溫度。

換熱單元面積約束：

式中，ACU,i、AHU,j和Ai,j,k分別為冷公用工程單元、熱公用工程單元和換熱器單元的面積。

換熱單元換熱量約束：

式中，QCU,i、QHU,j和Qi,j,k分別為冷公用工程單元、熱公用工程單元和換熱器單元的換熱量。

首先考察總換熱量與是否考慮換熱器權(quán)重得到的熱流均勻性因子之間的關(guān)系，所得結(jié)果如圖2所示，圖中φq,net為采用式（3）沒有考慮加權(quán)求出的熱流均勻性因子，而φ'q,net為式（5）考慮換熱量權(quán)重求得的熱流均勻性因子，Q為換熱網(wǎng)絡(luò)的總換熱量。從圖2中可以看出，總體上熱流均勻性因子能與總換熱量之間形成較好的對應(yīng)關(guān)系：熱流均勻性因子越小，則對應(yīng)的整個換熱網(wǎng)絡(luò)的換熱量越大，即換熱網(wǎng)絡(luò)換熱性能越好。比較圖中φq,net與φ'q,net可以看出，φ'q,net與Q之間形成的對應(yīng)關(guān)系更好，由此可以說明考慮各換熱器換熱量的權(quán)重計算得到的熱流均勻性因子能更準(zhǔn)確地反映出換熱網(wǎng)絡(luò)的換熱性能。

其次考察總換熱量與熱流均勻性因子和溫差均勻性因子之間的對應(yīng)關(guān)系，所得結(jié)果如圖3所示。圖中φ'q,net為熱流均勻性因子，φT,net為溫差均勻性因子，Q為換熱網(wǎng)絡(luò)的總換熱量。從圖3中可以看出，兩個因子總體上都能與換熱網(wǎng)絡(luò)的總換熱量之間形成較好的對應(yīng)關(guān)系，而熱流均勻性因子與總換熱量之間形成的對應(yīng)關(guān)系比溫差均勻性因子更佳，也就是說對于本文所取的各股流體換熱系數(shù)不同的換熱網(wǎng)絡(luò)來說，熱流均勻性因子更具適用性，能夠更準(zhǔn)確地反映出換熱網(wǎng)絡(luò)的換熱性能。

2.1.2 改變換熱網(wǎng)絡(luò)結(jié)構(gòu)條件下的熱流均勻性因子

圖2 換熱量與φq,net和φ'q,net的關(guān)系圖

圖3 換熱量與φ'q,net和φT,net的關(guān)系圖

上述內(nèi)容只是在同一個結(jié)構(gòu)中對熱流均勻性因子進(jìn)行討論，而換熱網(wǎng)絡(luò)的綜合及優(yōu)化總是針對不同的網(wǎng)絡(luò)結(jié)構(gòu)，不同結(jié)構(gòu)的換熱網(wǎng)絡(luò)總換熱量與本研究建立的熱流均勻性因子之間是否能形成較好的對應(yīng)直接關(guān)系到熱流均勻性因子能否作為換熱網(wǎng)絡(luò)換熱效能的評價指標(biāo)及能否指導(dǎo)換熱網(wǎng)絡(luò)的綜合及優(yōu)化等關(guān)鍵問題。本研究在圖1所示的結(jié)構(gòu)基礎(chǔ)上隨機(jī)改變各換熱器的冷、熱流體匹配，規(guī)定每股流體上至多匹配3個換熱器，同一位置至多匹配一個換熱器，以此取得30個不同的換熱網(wǎng)絡(luò)結(jié)構(gòu)，然后采用牛頓法對每個結(jié)構(gòu)進(jìn)行優(yōu)化，目前換熱網(wǎng)絡(luò)的優(yōu)化都以年綜合費(fèi)用為目標(biāo)，所以本研究以換熱網(wǎng)絡(luò)年綜合費(fèi)用F為目標(biāo)函數(shù)，換熱器單元面積A為變量，從而獲得30個換熱網(wǎng)絡(luò)的優(yōu)化結(jié)構(gòu)。約束條件如式（7）～式（9），目標(biāo)函數(shù)包括投資費(fèi)用和公用工程費(fèi)用，其表達(dá)式如式（10）。

式中，QCU,i、QHU,j分別為冷、熱公用工程換熱量；CCU、QHU分別為冷、熱公用工程費(fèi)用系數(shù)；ACU,i、AHU,j、Ai,j,k分別為冷公用工程單元、熱公用工程單元和換熱器單元的換熱面積；C1'、C2'、C3'分別為冷公用工程單元、熱公用工程單元和換熱器單元的面積費(fèi)用系數(shù)；C1、C2、C3分別為冷公用工程單元、熱公用工程單元和換熱器單元的固定投資費(fèi)用；B為面積費(fèi)用指數(shù)。

考察優(yōu)化后換熱網(wǎng)絡(luò)的總換熱量與其對應(yīng)的熱流均勻性因子之間的關(guān)系，所得結(jié)果如圖4所示。從圖4中可以看出，換熱網(wǎng)絡(luò)的總換熱量Q與熱流均勻性因子φ'q,net兩條曲線成相反趨勢，能夠形成換熱量越大，熱流均勻性因子越小的對應(yīng)關(guān)系，且在同一橫坐標(biāo)處換熱量有最大值時，熱流均勻性因子具有最小值。雖然存在一定的誤差但總體上能夠形成較好地對應(yīng)，由此可以說明本研究建立的換熱網(wǎng)絡(luò)熱流均勻性因子能作為優(yōu)化后換熱網(wǎng)絡(luò)換熱效能的評價指標(biāo)，從而能夠利用該熱流均勻性因子指導(dǎo)換熱網(wǎng)絡(luò)的綜合及優(yōu)化。

2.2 算例二

算例二采用文獻(xiàn)[21-23]中的10股流換熱器網(wǎng)絡(luò)來檢驗及分析換熱網(wǎng)絡(luò)中的熱流均勻性因子與換熱性能之間的關(guān)系，該換熱網(wǎng)絡(luò)由6股熱流體和4股冷流體組成，流股的初始參數(shù)見表2，文獻(xiàn)[21]中優(yōu)化的結(jié)構(gòu)如圖5所示。

在圖5所示的結(jié)構(gòu)基礎(chǔ)上隨機(jī)改變各換熱器的冷、熱流體匹配，規(guī)定每股流體上至多匹配4個換熱器，同一位置至多匹配一個換熱器，以此獲得30個不同的換熱網(wǎng)絡(luò)結(jié)構(gòu)，以換熱網(wǎng)絡(luò)年綜合費(fèi)用F為目標(biāo)函數(shù)，換熱器單元面積A為變量，采用牛頓法對每個結(jié)構(gòu)進(jìn)行優(yōu)化，獲得30個換熱網(wǎng)絡(luò)的優(yōu)化結(jié)構(gòu)。

圖4 換熱量與熱流均勻性因子φ'q,net的關(guān)系圖

表2 算例2的流體參數(shù)

圖5 文獻(xiàn)[21]優(yōu)化的換熱網(wǎng)絡(luò)結(jié)構(gòu)

圖6 換熱量與熱流均勻性因子φ'q,net的關(guān)系圖

考察優(yōu)化后換熱網(wǎng)絡(luò)的總換熱量與其對應(yīng)的熱流均勻性因子之間的關(guān)系，所得結(jié)果如圖6所示。從圖6中可以看出，換熱網(wǎng)絡(luò)的總換熱量Q與熱流均勻性因子φ'q,net之間總體上能夠形成換熱量越大，熱流均勻性因子越小的對應(yīng)關(guān)系。通過對該算例中熱流均勻性因子的驗證，說明本研究所建立的換熱網(wǎng)絡(luò)的熱流均勻性因子能夠作為優(yōu)化后換熱網(wǎng)絡(luò)換熱效能的一種新的評價指標(biāo)，且該熱流均勻性因子能為換熱網(wǎng)絡(luò)的綜合及優(yōu)化提供一種新的思路。

3 結(jié)論

本研究從能量有效利用的角度出發(fā)，在換熱網(wǎng)絡(luò)的溫差均勻性因子及場協(xié)同理論的基礎(chǔ)上，提出并建立了換熱網(wǎng)絡(luò)的熱流均勻性因子，通過兩個具體的實(shí)例對換熱網(wǎng)絡(luò)熱流均勻性因子進(jìn)行驗證及分析，得到以下結(jié)論。

（1）根據(jù)各換熱器的換熱量進(jìn)行加權(quán)計算得到的熱流均勻性因子與總換熱量之間能形成更好的對應(yīng)關(guān)系，能夠更準(zhǔn)確地反映換熱網(wǎng)絡(luò)的換熱效能。

（2）換熱網(wǎng)絡(luò)的熱流均勻性因子比溫差均勻性因子更適于評價換熱網(wǎng)絡(luò)的換熱效能，使用的范圍更加廣泛。

（3）不同優(yōu)化結(jié)構(gòu)的換熱網(wǎng)絡(luò)熱流均勻性因子與總換熱量之間能夠形成較好的對應(yīng)關(guān)系。

總之，本研究所建立的換熱網(wǎng)絡(luò)的熱流均勻性因子總體上能夠反映換熱網(wǎng)絡(luò)的換熱效能，該因子能夠作為換熱網(wǎng)絡(luò)換熱效能的一種評價指標(biāo)，并能利用該因子指導(dǎo)換熱網(wǎng)絡(luò)的綜合及優(yōu)化。

[1] Linnhoff B，Hindmarsh E.The pinch design method for heat exchanger networks[J].Chemical Engineering Science，1983，38（5）：745-763.

[2] 趙野，孫琳，羅雄麟.多程換熱網(wǎng)絡(luò)綜合與夾點(diǎn)技術(shù)研究進(jìn)展[J].化工進(jìn)展，2012，31（8）：1685-1689.

[3] Yee T F，Grossmann I E.Simultaneous optimization for heat integration-II.Heat exchanger network synthesis[J].Computers and Chemical Engineering，1990，14（10）：1165-1183.

[4] 周聰海，黃智賢，邱挺，等.隨機(jī)搜索算法在換熱網(wǎng)絡(luò)優(yōu)化中的應(yīng)用進(jìn)展[J].化工進(jìn)展，2012，31（3）：495-501.

[5] Ravagnani M A S S，Silva A P，Arroyo P A，et al.Heat exchanger network synthesis and optimization using genetic algorithm[J].Applied Thermal Engineering，2005，25（7）：1003-1017.

[6] Lin B，Miller D C.Solving heat exchanger network synthesis problems with Tabu search[J].ComputersandChemical Engineering，2004，28（8）：1451-1464.

[7] 過增元，李志信，周森泉，等.換熱器中的溫差場均勻性原則[J].中國科學(xué)，1996，26（1）：25-31.

[8] 過增元.熱流體學(xué)[M].北京：清華大學(xué)出版社，1992.

[9] 呂巖巖，崔國民，郭佳，等.多股流換熱器設(shè)計的量綱1溫差均勻性優(yōu)化因子[J].化工學(xué)報，2007，58（10）：2476-2480.

[10] 崔國民，呂巖巖，張勤.換熱器網(wǎng)絡(luò)溫差場分析及溫差均勻性因子的建立[C]//中國工程熱物理學(xué)會傳熱傳質(zhì)學(xué)術(shù)會議論文集，南京，2006：1608-1611.

[11] 呂巖巖.場協(xié)同下的換熱器網(wǎng)絡(luò)整體強(qiáng)化技術(shù)研究[D].上海：上海理工大學(xué)，2009.

[12] 崔國民，呂巖巖，胡向柏.換熱網(wǎng)絡(luò)整體強(qiáng)化技術(shù)及其溫差均勻性因子[C]//中國工程熱物理學(xué)會傳熱傳質(zhì)學(xué)術(shù)會議論文集，鄭州，2008：454.

[13] H.哈肯.協(xié)同學(xué)：理論與應(yīng)用[M].楊炳弈譯.北京：中國科學(xué)技術(shù)出版社，1990.

[14] 過增元，黃素逸，等.場協(xié)同原理與強(qiáng)化傳熱新技術(shù)[M].北京：中國電力出版社，2004.

[15] 崔國民，盧洪波，李美玲.場協(xié)同條件下的換熱器恒熱流換熱研究[C]//中國工程熱物理學(xué)會傳熱傳質(zhì)學(xué)會議論文集，上海，2002：1029-1032.

[16] 過增元.對流換熱的物理機(jī)制及其控制：速度場與熱流場的協(xié)同[J]. 科學(xué)通報，2000，45（19）：2118-2122.

[17] Bjork Kaj-Mikael，Nordman Roger.Solving large-scale retrofit heat exchanger network synthesis problems with mathematical optimization methods[J].Chemical Engineering and Processing，2005，44（8）：869-876.

[18] 胡向柏，崔國民，涂惟民，等.換熱網(wǎng)絡(luò)填充函數(shù)法的全局優(yōu)化[J].化學(xué)工程，2011，39（1）：28-31.

[19] 張勤，崔國民，關(guān)欣.基于蒙特卡羅遺傳算法的換熱網(wǎng)絡(luò)優(yōu)化問題[J].石油機(jī)械，2007，35（5）：19-22.

[20] 胡向柏，崔國民，許海珠，等.換熱器優(yōu)化順序?qū)Q熱網(wǎng)絡(luò)全局優(yōu)化的影響[J].化工進(jìn)展，2012，31（5）：987-1003.

[21] Khorasany R M，F(xiàn)esanghary M.A novel approach for synthesis of cost-optimal heat exchanger networks[J].Comput.Chem.Eng.，2009，33（8）：1363-1370.

[22] Yerramsetty K M，Murty C V S.Synthesis of cost-optimal heat exchanger networks using differential evolution[J].Comput.Chem.Eng.，2008，32（8）：1861-1876.

[23] Ravagnani M A S S，Silva A P，Arroyo P A，et al.Heat exchanger network synthesis and optimization using genetic algorithm[J].Appl.Therm.Eng.，2005，25（7）：1003-1017.