彭小奇 ，宋國輝，宋彥坡 ，張建智，劉振國

(1. 中南大學 能源科學與工程學院，湖南 長沙，410083；2. 中南大學 信息科學與工程學院，湖南 長沙，410083；3. 湖南第一師范學院 信息科學與工程系，湖南 長沙，410205)

我國氧化鋁生產(chǎn)與工業(yè)發(fā)達國家相比，在能耗方面存在較大差距[1]。其中，蒸發(fā)工序擔負著平衡全廠水和堿的作用，能耗尤為顯著：其汽耗占總汽耗的30%～40%，成本占總成本的10%～12%[2]。因此，降低蒸發(fā)工序的能耗對降低氧化鋁生產(chǎn)成本具有重要意義。為提高蒸發(fā)工序的運行性能和用能效率，氧化鋁生產(chǎn)科研工作者進行了長期探索和研究，其主要成果可概括為2類：通過更新升級設備提高蒸發(fā)系統(tǒng)或單元的性能[3-4]；根據(jù)長期運行經(jīng)驗或工業(yè)試驗優(yōu)化操作參數(shù)或方式[5]。然而，關于蒸發(fā)工序能耗的熱力學分析研究報道較少或不夠深入。為降低蒸發(fā)工序的能耗，使用熱力學分析方法研究其能耗狀況非常必要。作為一種更科學的熱力學分析方法，分析已用于指導眾多部門的節(jié)能工作，如電力、化工和鋼鐵等[6-8]。但分析在有色冶金工業(yè)尤其是氧化鋁生產(chǎn)上的應用很少。郭沈[9]曾將分析應用于間接加熱的預脫硅過程，但未將物料的化學考慮在內。目前，尚無氧化鋁生產(chǎn)的主要中間物流，即工業(yè)鋁酸鈉溶液值計算式的報道。而溶液的值，尤其是化學的計算，需要使用溶液的密度、比熱容，各組分的活度因子和標準化學等性質，常常成為應用分析的一大障礙。目前，已有成熟的工業(yè)鋁酸鈉溶液密度和比熱容的計算方法[10-12]；Zhou等[13]和Li等[14]分別用實驗和理論方法研究了NaOH-NaAl(OH)4-H2O體系的活度因子，彭小奇等[15]發(fā)展了一種基于 Bromley模型的 NaOHNaAl(OH)4-Na2CO3-H2O體系活度因子的計算模型；鄭丹星等[16]修正了龜山-吉田環(huán)境模型并更新了元素的標準化學。本文作者將在此基礎上，推導并討論工業(yè)鋁酸鈉溶液的計算式；對中國鋁業(yè)中州分公司四效蒸發(fā)器-三級閃蒸器系統(tǒng)的進行全面深入分析，以便尋找系統(tǒng)的用能薄弱環(huán)節(jié)，確定能耗損失的部位、大小及原因，為進一步的節(jié)能改造、優(yōu)化指明方向。

1 流程簡介

中國鋁業(yè)中州分公司蒸發(fā)V組所采用的四效管式降膜蒸發(fā)器-三級閃蒸器蒸發(fā)系統(tǒng)如圖 1所示，其中除E4(末)效蒸發(fā)器無預熱器外，其他各效均有預熱器。新蒸汽進入 E1效，經(jīng)換熱后變?yōu)槔淠僖来谓?jīng)過1-1號、1-2號和1-3號等冷水罐回收部分余熱，其余冷凝水被排出體系，而相應3個冷凝水罐中產(chǎn)生的二次蒸汽分別用于加熱E1，E2和E3效蒸發(fā)器。蒸發(fā)原液首先分別進入E4和E3效蒸發(fā)器，然后依次經(jīng)過E2和E1效蒸發(fā)器以及S1，S2和S3級閃蒸器蒸發(fā)濃縮。E4效蒸發(fā)器產(chǎn)生的乏汽被水冷卻后排出體系，其余各效產(chǎn)生的二次蒸汽依次分別用于加熱前一效蒸發(fā)器和本效預熱器。三級閃蒸器產(chǎn)生的二次蒸汽分別供入E1，E2和E3效預熱器，與上一效蒸發(fā)器的出料進行混合預熱。

2 工業(yè)鋁酸鈉溶液的值計算

圖1 中國鋁業(yè)中州分公司蒸發(fā)工序流程圖Fig.1 Flow sheet of the evaporation process in Zhongzhou Branch of Chalco

2.1 物理

根據(jù)相關的計算方法[18-19]，考慮到偏摩爾性質和易用性，將工業(yè)鋁酸鈉溶液密度[10]和比熱容[11]的擬合式代入式(1)，得其比物理的計算式為：

式中：ρNK，ρA和ρNC分別為工業(yè)鋁酸鈉溶液中苛性堿、氧化鋁和碳酸堿的質量濃度，g/L；ρs為工業(yè)鋁酸鈉溶液的密度，kg/m3。

在蒸發(fā)工序運行條件下，工業(yè)鋁酸鈉溶液的密度為 1 200～1 400 kg/m3且|p-p0|＜0.5，故式(2)末項(即(p-p0)×103/ρs)的值相對較小。實測數(shù)據(jù)計算顯示，式(2)末項的值占式(2)值的-3%～2%，因此，可忽略式(2)末項，即忽略壓力對工業(yè)鋁酸鈉溶液物理的影響，且在熱工測試中可以省去對某些測試難度較大位置的壓力測試。這既可減少測試工作量，也有利于設備安全運行和壽命維護。

2.2 化學

式中：xi為溶液中第i種組分的摩爾分數(shù)；ei為溶液中第i種組分的標準化學，kJ/mol；fi為使用摩爾分數(shù)時溶液中第i種組分的活度因子；R為摩爾氣體常數(shù)，8.314 J/(mol·K)。

根據(jù)熱力學基本原理，將式(4)和(5)改寫為：

式中：mi為溶液中第i種組分的質量摩爾濃度，mol/kg；γi為使用質量摩爾濃度時溶液中第 i種組分的活度因子，可由 NaOH-NaAl(OH)4-Na2CO3-H2O體系的活度系數(shù)計算模型[15]獲得。

將RT0/1 000=2.479代入式(7)，并經(jīng)單位換算和簡化后，可得工業(yè)鋁酸鈉溶液比化學的計算式為：

式中 ：ρw為單位體積工業(yè)鋁酸鈉溶液所含水的質量，kg/m3。

以某次穩(wěn)定運行狀態(tài)下熱工測試所得數(shù)據(jù)為例，根據(jù)式(2)和(8)計算蒸發(fā)原液和各效(級)出料的比物理、比反應和比擴散，結果如表1所示。由表1可見：比反應遠大于其他兩者之和，因此，總對其變化非常敏感，在具有化學反應的其他氧化鋁生產(chǎn)工序的分析中，需注意成分分析的精度。而蒸發(fā)工序無化學反應，理論上溶質各組分的總質量不變，實測數(shù)據(jù)也近似守恒，故可認為蒸發(fā)料液的反應總量不變，可不計入反應項。在表1中，比擴散出現(xiàn)了正、負兩類值，這與環(huán)境模型的規(guī)定和活度因子有關，不影響分析結果的正確性。

表1 蒸發(fā)料液的工藝參數(shù)及其值Table1 Process parameters and exergy values of materials in evaporation process

表1 蒸發(fā)料液的工藝參數(shù)及其值Table1 Process parameters and exergy values of materials in evaporation process

注：E1～E4和S1～S3分別表示相應蒸發(fā)器和閃蒸器的出口料液。

?

3 蒸發(fā)工序的效率分析

式中：Eg為系統(tǒng)或單元的收益(系統(tǒng)或E4效蒸發(fā)器的收益是進出蒸發(fā)系統(tǒng)或E4效蒸發(fā)器的物料的增；其他單元的收益是相應單元的出料(包括增濃后的料液和二次蒸汽)與進料相比的増)；Ep為系統(tǒng)或單元的支付(系統(tǒng)或E1效蒸發(fā)器的支付是新蒸氣供入的；其他單元的支付是相應蒸發(fā)器或閃蒸器所供入的二次蒸汽的)。

圖2 蒸發(fā)系統(tǒng)及各單元的效率Fig.2 Exergy efficiencies of evaporation system and its units

4 四效蒸發(fā)器的損失分析

4.1 外部損失

在3次測試工況下，蒸發(fā)系統(tǒng)及各效冷凝水和末效乏汽外排形式的外部損系數(shù)rEX見表2。

表2 冷凝水和末效乏汽的損系數(shù)Table2 Exergy loss coefficients of condensed waters and exhausted vapor

表2 冷凝水和末效乏汽的損系數(shù)Table2 Exergy loss coefficients of condensed waters and exhausted vapor

注：rEXV(E4)為E4效蒸發(fā)器外排乏汽引起的外部損系數(shù)；rEXW(Ei)為 Ei效蒸發(fā)器外排冷凝水引起的外部損系數(shù)；rEX(SYS)為整個蒸發(fā)系統(tǒng)的外部損系數(shù)。

?

測試數(shù)據(jù)表明：1-3號和4號冷凝水罐排水的平均溫度約為95 ℃，末效冷凝水和乏汽的平均溫度約為60 ℃，而蒸發(fā)原液的平均溫度約為49 ℃?？梢?，1-3號和4號冷凝水罐的排水與蒸發(fā)原液間溫差較大；另外，雖然末效乏汽溫度較低，但其相變潛熱較大。因此，用冷凝水和乏汽預熱蒸發(fā)原液是一條值得探索的節(jié)能途徑。

4.2 內部損失

表3 四效蒸發(fā)器的傳熱損系數(shù)Table3 Exergy loss coefficients due to heat transfer in four-effect evaporations

表3 四效蒸發(fā)器的傳熱損系數(shù)Table3 Exergy loss coefficients due to heat transfer in four-effect evaporations

?

E1，E2和E3等三效蒸發(fā)器中預熱器內混合過程造成的混合損失均較小，系統(tǒng)的總混合損失相應較小。在3次測試工況下，蒸發(fā)系統(tǒng)的混合損系數(shù)rMT分別為0.023，0.016和0.030。

5 結論

(4) 為了實現(xiàn)蒸發(fā)工序的節(jié)能，建議加強冷凝水和二次蒸汽的余熱回收利用，優(yōu)化蒸發(fā)系統(tǒng)操作參數(shù)及傳熱溫差分布，并改進預熱器的使用方式。

[1] LIU Li-ru, LU A-ye, Lu Zhong-wu, et al. Analysis of the overall energy intensity of alumina refinery process using unit process energy intensity and product ratio method[J]. Energy, 2006,31(8/9): 1167-1176.

[2] 李續(xù)平, 周在毅. 氧化鋁蒸發(fā)節(jié)能降耗途徑的探討[J]. 輕金屬, 2006(8): 37-39.

LI Xu-ping, ZHOU Zai-yi. Discussion about energy saving for alumina evaporation process[J]. Light Metals, 2006(8): 37-39.

[3] Green J. Opportunities to reduce energy intensity in the aluminum industry[EB/OL]. [2002-03-26]. http://www.eere.energy.gov/evision/pdfs/presentation s/j_green.pdf.

[4] 孫克萍, 先晉聰, 宋強, 等. 我國氧化鋁業(yè)中蒸發(fā)裝置技術進步及效能分析[J]. 有色金屬, 2004(3): 43-45.

SUN Ke-ping, XIAN Jin-chong, SONG Qiang, et al. Technical progress and analysis of efficiency and energy on the evaporation installations to Chinas alumina industry[J]. Nonferrous Metallurgy, 2004(3): 43-45.

[5] 楊新民, 吳紅應, 梅永峰. 氧化鋁生產(chǎn)中降低蒸汽消耗的生產(chǎn)實踐[J]. 輕金屬, 2008(9): 28-30.

YANG Xin-min, WU Hong-ying, MEI Yong-feng. The practice of reducing steam consumption during producing alumina[J].Light Metals, 2008(9): 28-30.

[6] Kamate S C, Gangavati P B. Exergy analysis of cogeneration power plants in sugar industries[J]. Applied Thermal Engineering, 2009, 29(5/6): 1187-1194.

[7] Antonio B A, Romildo P B, Luis S V, Exergetic analysis of distillation processes: A case study[J]. Energy, 2007, 32(7):1185-1193.

[8] Costa M M, Schaeffer R, Worrell E. Exergy accounting of energy and materials flows in steel production systems[J].Energy, 2001, 26(4): 363-384.

[9] 郭沈. 間接加熱脫硅系統(tǒng)的有效能分析[J]. 輕金屬, 2005(8):11-14.

GUO Shen. Exergy analysis of indirect heating desilication process[J]. Light Metals, 2005(8): 11-14,

[10] 楊重愚. 氧化鋁生產(chǎn)工藝學[M]. 北京: 冶金工業(yè)出版社,1982: 37-41.

YANG Zhong-yu. The technology of alumina production[M].Beijing: Metallurgical Industry Press, 1982: 37-41.

[11] 畢詩文, 于海燕. 氧化鋁生產(chǎn)工藝[M]. 北京: 冶金工業(yè)出版社, 2006: 23-27.

BI Shi-wen, YU Hai-yan. The technology of alumina production[M]. Beijing: Chemical Industry Press, 2006: 23-27.

[12] Reynolds J G , Carter R. Density model for sodium hydroxidesodium aluminate solutions[J]. Hydrometallurgy, 2007, 89(3/4):233-241.

[13] ZHOU Jun, CHEN Qi-yuan, LI Jie, et al. Isopiestic measurement of the osmotic and activity coefficients for the NaOHNaAl(OH)4-H2O system at 313.2 K[J]. Geochim Cosmochim Ac,2003: 67(18): 3459-3472.

[14] LI Xiao-bin, LU Wei-jun, LIU Gui-hua, et al. Activity coefficient calculation model for NaAl(OH)4-NaOH-H2O system[J].Transactions of Nonferrous Metal Society of China, 2005, 15(4):908-912.

[15] 彭小奇, 宋國輝, 宋彥坡, 等. NaOH-NaAl(OH)4-Na2CO3-H2O體系活度系數(shù)的計算模型[J]. 中國有色金屬學報, 2009, 19(7):1332-1337.

PENG Xiao-qi, SONG Guo-hui, SONG Yan-po, et al. Activity coefficient calculation model for NaOH-NaAl(OH)4-Na2CO3-H2O system[J]. Chinese Journal of Nonferrous Metals, 2009,19(7): 1332-1337.

[16] 鄭丹星, 武向紅, 鄭大山.函數(shù)的一致性基礎[J]. 化工學報,2002, 53(7): 673-679.

ZHENG Dan-xing, WU Xiang-hong, ZHENG Da-shan.Fundamental of thermodynamic consistence of exergy functions[J]. Journal of Chemical Industry and Engineering,2002, 53(7): 673-679.

[17] Tsatsaronis G. Definitions and nomenclature in exergy analysis and exergoeconomics[J]. Energy, 2007, 32(4): 249-253.

[18] Szargut J, Morris D R, Steward F R. Exergy analysis of thermal,chemical and metallurgical process[M]. New York: Hemisphere Publishing, 1988: 51-53.

[19] 宋國輝, 彭小奇, 宋彥坡, 等. 含相變過程的機械計算方法[J]. 熱科學與技術, 2009, 8(1): 64-68.

SONG Guo-hui, PENG Xiao-qi, SONG Yan-po, et al. Method for mechanical exergy calculation including phase change process[J]. Journal of Thermal Science and Technology, 2009,8(1): 64-68.

[20] Barin I, Gregor P G. Thermochemical data of pure substances[M]. New York: Wiley-VCH, 2004: 1117-1136.