基于Aspen Plus醋酸乙烯精餾塔的模擬優(yōu)化

陳 果 王 煤 程 橙 涂開慧

（四川大學(xué)化學(xué)工程學(xué)院，成都 610065）

研究與開發(fā)

基于Aspen Plus醋酸乙烯精餾塔的模擬優(yōu)化

陳 果 王 煤 程 橙 涂開慧

（四川大學(xué)化學(xué)工程學(xué)院，成都 610065）

利用Aspen Plus軟件，選擇NRTL活度系數(shù)方程和Hayden-O′connell逸度系數(shù)方程的熱力學(xué)模型，應(yīng)用系統(tǒng)中的RadFrac精餾模塊對(duì)醋酸乙烯精餾塔進(jìn)行模擬，模擬值與實(shí)際值基本吻合。討論了進(jìn)料位置、回流比、塔頂側(cè)線采出量等參數(shù)對(duì)精餾分離精度與能耗的影響，提出優(yōu)化方案為：進(jìn)料板為第62塊，回流比為32，側(cè)線采出質(zhì)量流量為37.5 t/h。此參數(shù)下，重新進(jìn)行計(jì)算，結(jié)果顯示，塔頂冷凝器和塔釜再沸器的熱流量分別降低了15.5%和16.9%，塔頂側(cè)線采出液中醋酸乙烯和塔釜采出液中醋酸的質(zhì)量分?jǐn)?shù)分別上升了0.4%和0.13%。

Aspen Plus；醋酸乙烯；醋酸；精餾塔；模擬；優(yōu)化

Aspen Plus化工模擬系統(tǒng)是美國(guó)麻省理工學(xué)院于20世紀(jì)70年代后期開發(fā)的大型化工模擬軟件，20世紀(jì)80年代初由美國(guó)Aspen Plus技術(shù)公司推向市場(chǎng)。Aspen Plus用嚴(yán)格的方法進(jìn)行單元和全過程計(jì)算，為企業(yè)提供了準(zhǔn)確的單元操作模型，該系統(tǒng)功能齊全、規(guī)模龐大，目前已在世界范圍內(nèi)廣泛使用參考[1]。

在醋酸乙烯（VAC）生產(chǎn)工藝中，VAC精餾塔在整個(gè)流程中占據(jù)著重要的位置，它直接影響著VAC的質(zhì)量和合成后期聚合物的好壞，以及整個(gè)生產(chǎn)過程的經(jīng)濟(jì)效益。對(duì)VAC精餾及精餾塔的分析研究前人已經(jīng)做過一些，但主要是采用不同的化工模擬軟件或者是僅從個(gè)別操作條件的改變對(duì)分離精度的影響，而對(duì)不同操作參數(shù)和工藝條件的改變對(duì)VAC精餾塔分離精度和能耗的影響尚未進(jìn)行系統(tǒng)的分析和研究，特別是塔頂側(cè)線采出量對(duì)分析精度的影響。本文中以增加VAC產(chǎn)量和節(jié)約能源為目的，通過Aspen Plus模擬系統(tǒng)對(duì)VAC精餾塔進(jìn)行了詳細(xì)的模擬計(jì)算，對(duì)VAC精餾塔的操作參數(shù)進(jìn)行了靈敏度分析，以便為工業(yè)生產(chǎn)中VAC精餾塔的設(shè)計(jì)提供基礎(chǔ)數(shù)據(jù)。

1 VAC精餾工藝

某廠由天然氣乙炔路線生產(chǎn)的粗VAC經(jīng)脫氣、粗餾等系統(tǒng)除去溶解乙炔氣及大部分乙醛、丙酮及水等輕組分后作為VAC精餾塔一進(jìn)料，塔釜用再沸器加熱，塔一塔頂氣相進(jìn)入塔二底部，塔二的釜液經(jīng)回流泵供塔一回流。在塔二中部側(cè)線液相采出精VAC，塔一塔底采出粗醋酸（HAC）去HAC精餾塔，從塔二塔頂出來(lái)的氣相經(jīng)冷凝部分回流入塔。進(jìn)料及產(chǎn)品組成見表1。

表1 進(jìn)料及產(chǎn)品組成Tab 1 Composition of feed and product

2 模型的建立及模擬結(jié)果

2.1 熱力學(xué)方法的選取

在用Aspen Plus進(jìn)行模擬計(jì)算時(shí)，熱力學(xué)模型的選取尤為重要，其正確與否直接影響計(jì)算的物理性能的準(zhǔn)確程度與計(jì)算結(jié)果的精確度。在分離技術(shù)中，相平衡是精餾操作的基礎(chǔ)，其汽-液平衡關(guān)系計(jì)算如下[2]：

式中，yi、xi分別為氣、液相中組分i的摩爾分?jǐn)?shù)，φiv為氣相混合物中組分i的逸度系數(shù)，fi*l為純組分i在混合物溫度下的液相逸度；γi為液相組分i的活度系數(shù)；p為體系的壓力。

2.1.1 活度系數(shù)模型

液相活度系數(shù)γi由活度系數(shù)模型計(jì)算，NRTL方程是以威爾遜提出的局部組成概念為基礎(chǔ)的活度系數(shù)方程，該模型能準(zhǔn)確模擬非理想溶液的VLE和LLE性質(zhì)[3]。根據(jù)本體系的特點(diǎn)，選用NRTL作為計(jì)算采用的活度系數(shù)模型。NRTL方程為[4-5]：

式中：Gji=exp(-ajiτji)為溶液相互作用自由能，i、j、k、m為不同組分。其中，τji為組分i-j之間作用能參數(shù)（τji=aji+(bji/T)+ejiln T+fjiT）；aji為溶液可調(diào)的特征函數(shù)，定在 0.20～0.47；aji=cji+dji(T-273.15 K)，為 NRTL方程二元參數(shù)。 模型需要的二元參數(shù) aji、bji、cji、dji、eji和fji等可通過VLE和（或）LLE的數(shù)據(jù)回歸來(lái)確定，在Aspen Plus數(shù)據(jù)庫(kù)中有很多自文獻(xiàn)及通過數(shù)據(jù)回歸得到的相關(guān)二元參數(shù)。

2.1.2 逸度系數(shù)模型

氣相逸度系數(shù)φiv由狀態(tài)方程計(jì)算。Aspen plus推薦選擇Hayden-O′Connell方程來(lái)計(jì)算中壓以下的氣相締合現(xiàn)象，該方程能預(yù)測(cè)含有極性組分或有締合成二聚體傾向的體系，適用于本體系[6]。

由第2 Virial系數(shù)關(guān)聯(lián)式計(jì)算逸度系數(shù)的方程為[7]：

Aspen Plus中包含活度系數(shù)NRTL方程和逸度系數(shù) Hayden-O′Connell方程的熱力學(xué)模型為NRTL-HOC，故本次模擬在 Aspen plus的Properties模塊中選擇NRTL-HOC模型。

2.2 精餾模塊的選取

RadFrac是一個(gè)嚴(yán)格模型，用于模擬所有類型的多級(jí)氣液平衡操作。這些操作包括一般精餾、吸收、再沸吸收、汽提、再沸汽提、萃取和共沸蒸餾等過程。對(duì)精餾過程，它適用于2相蒸餾體系、3相蒸餾體系、窄沸程和寬沸程體系，液相具有強(qiáng)非理想性的體系的計(jì)算[9]。因此，本文選用RadFrac精餾模塊。

2.3 模擬結(jié)果

在如表2所示的塔初始操作參數(shù)下，應(yīng)用Aspen plus，選擇Properties中的NRTL-HOC模型，利用RadFrac模塊對(duì)VAC精餾塔進(jìn)行精算，結(jié)果如表3。

表2 VAC精餾塔初始操作參數(shù)Tab 2 Initlal operation Parameters of vinyl acetate distillation

表3 模擬結(jié)果與實(shí)際操作結(jié)果Tab 3 The results of simulation and actual operation

從表3可以看出，VAC精餾塔側(cè)線采出VAC和塔底HAC的質(zhì)量分?jǐn)?shù)計(jì)算值與工廠生產(chǎn)值相比，相對(duì)誤差在5%以內(nèi)，故該模型的計(jì)算結(jié)果與實(shí)際的吻合性很好，可以用此模型進(jìn)行進(jìn)一步模擬分析；分離效果已經(jīng)較好，但如果進(jìn)一步提高VAC精餾塔的分離效果以及降低能耗，必須優(yōu)化VAC精餾塔的初始操作參數(shù)。

3 操作參數(shù)的優(yōu)化

ASPEN中的靈敏度分析工具可以方便地確定過程對(duì)關(guān)鍵操作變量和設(shè)計(jì)變量的響應(yīng)，即一個(gè)或多個(gè)流程變量變化對(duì)其他流程變量造成的影響，這是一種“如果…就”形式的研究工具[10]。

VAC精餾塔的分離效果可以用塔頂側(cè)線采出VAC含量，塔底HAC含量來(lái)判斷，含量越高，分離效果越好、能耗越低、經(jīng)濟(jì)效益越好。為了優(yōu)化分離效果降低能耗，本文以側(cè)線采出VAC含量，塔底HAC含量（質(zhì)量分?jǐn)?shù)）為指標(biāo)，在表2的初始參數(shù)下，利用Aspen Plus的Sensitivity模塊，探討原料進(jìn)料位置、回流比以及側(cè)線采出量等參數(shù)對(duì)分離精度和能耗的影響，來(lái)尋求最適宜的操作條件，以達(dá)到最大的經(jīng)濟(jì)效益。

3.1 進(jìn)料位置的影響

在相同的條件下，進(jìn)料位置的不同將對(duì)塔板的分離效果產(chǎn)生影響。不同的精餾效果導(dǎo)致塔頂和塔釜的組成產(chǎn)生改變，對(duì)再沸器和冷凝器負(fù)荷也會(huì)產(chǎn)生影響，因此存在最佳進(jìn)料位置。原料進(jìn)料位置對(duì)側(cè)線采出的VAC、塔底HAC含量和冷凝器和再沸器總熱流量的影響見圖1。

從圖1可以看出，當(dāng)進(jìn)料位置在55～65精餾效果最好，且冷凝器和再沸器總熱流量較小，而進(jìn)料在第62塊時(shí)進(jìn)料達(dá)到最佳效果，此時(shí)塔頂側(cè)線采出VAC的質(zhì)量分?jǐn)?shù)為0.999 6，塔底的HAC的質(zhì)量分?jǐn)?shù)為0.995 6，因此取最佳進(jìn)料位置為第62塊。

3.2 回流比的影響

在精餾塔中，回流比將直接影響每一層塔板的分離能力，因此對(duì)于精餾塔，無(wú)論是設(shè)計(jì)還是操作，回流比都是一個(gè)至關(guān)重要的參數(shù)?；亓鞅葘?duì)側(cè)線采出的VAC、塔底HAC含量和冷凝器、再沸器熱流量的影響見圖2。

從圖2可以看出，當(dāng)回流比達(dá)到32以上時(shí)，增大回流比對(duì)精餾效果的影響已經(jīng)不是很明顯，但由于冷凝器和再沸器的熱流量隨回流比的增加而增加也即意味能源費(fèi)的增加，因此過大的回流比是沒有必要的，綜合工藝和經(jīng)濟(jì)要求選擇回流比為32比較合適。

3.3 塔頂側(cè)線采出對(duì)精餾效果的影響

塔頂側(cè)線采出流量對(duì)側(cè)線采出的VAC、塔底HAC含量的影響見圖3。

從圖3可以看出，在塔頂側(cè)線采出質(zhì)量流量為37.5 t/h后，產(chǎn)品中VAC含量急劇下降，塔底HAC精餾效果趨于平穩(wěn)；而進(jìn)料液中VAC的質(zhì)量流量為39.0 t/h，采出量的增加，會(huì)增大冷凝器和再沸器的負(fù)荷，因此可以考慮在側(cè)線產(chǎn)品VAC含量較高時(shí)選擇合適的采出量。通過模擬計(jì)算可知，選擇塔頂側(cè)線采出量為37.5 t/h是最好的。

綜上所述，物料的進(jìn)料位置、回流比和塔頂（或塔底）的出料量是影響精餾效果的主要因素，選擇合適的操作條件才能滿足技術(shù)和經(jīng)濟(jì)2方面的要求，提出優(yōu)化方案為：原料進(jìn)料位置為第62塊板，回流摩爾比為32，側(cè)線采出質(zhì)量流量為37.5 t/h。

4 優(yōu)化參數(shù)方案下的模擬計(jì)算

在上述優(yōu)化參數(shù)下，對(duì)精餾塔模型進(jìn)行模擬計(jì)算，優(yōu)后技術(shù)指標(biāo)見表4。

表4 優(yōu)化后的技術(shù)指標(biāo)Tab 4 Comparison of technical index before and after optimization

將表4與表3對(duì)比，可以看出，塔頂冷凝器的熱流量從9.364 MW降低到7.908 MW，即降低了15.5%。塔釜再沸器的熱流量從9.618 MW降低到7.990 MW，即降低了16.9%。塔頂側(cè)線采出液中VAC的質(zhì)量分?jǐn)?shù)由0.999 2上升到0.999 6，即上升了0.4%。塔釜采出液中HAC的質(zhì)量分?jǐn)?shù)由0.994 1到0.995 3，上升了0.13%。說明對(duì)原工藝的優(yōu)化達(dá)到了預(yù)期的目的，優(yōu)化方案是成功的。

5 結(jié)論

1)應(yīng)用Aspen Plus化工模擬系統(tǒng)中RadFrac精餾模型，選擇NRTL活度系數(shù)方程和Hayden-O′connell逸度系數(shù)方程的熱力學(xué)模型對(duì)VAC精餾塔進(jìn)行模擬計(jì)算，其計(jì)算結(jié)果與實(shí)際運(yùn)行結(jié)果基本吻合，設(shè)計(jì)的模型可以用來(lái)模擬該工藝過程。

2)通過Aspen Plus的靈敏度分析功能，尋找到了精餾塔操作條件中的最靈敏參數(shù)：進(jìn)料板為第62塊，回流摩爾比為32，側(cè)線采出質(zhì)量流量為37.5 t/h。在優(yōu)化參數(shù)條件下，塔頂冷凝器和塔釜再沸器的熱流量分別降低了15.5%和16.9%，塔頂側(cè)線采出液中VAC和塔釜采出液中HAC的質(zhì)量分?jǐn)?shù)分別上升了0.4%和0.13%，對(duì)實(shí)際操作有指導(dǎo)意義。

[1]謝揚(yáng),沈慶揚(yáng).ASPEN PLUS化工模擬系統(tǒng)在精餾過程中的應(yīng)用[J].化工生產(chǎn)與技術(shù),1999,16(3):17-22,28.

[2]陳鐘秀,顧飛燕,胡望明.化工熱力學(xué)[M].2版.北京:化學(xué)工業(yè)出版社,2001.

[3]鄭丹星.化工熱力學(xué)教程[M].北京:中國(guó)石化出版社,2000.

[4]Renon H,Prausnitz J M.Local compositions in thermodynamic excess functions for liquid mixtures.AIChE[J].1968,14(l):135-144.

[5]Schmidt K A G,MahamY,Mather A E.Use of the NRTL equationforsimultaneouscorrelationofvapour-liquidequilibria and exeess enthalpy[J].Journal of thermal analysis and calorimetry,2007,89(l):61-72.

[6]Hsieh C T,Ji W Y,Lin H M.MultiPhase equilibria for mixtures Containing water,acetic acid,propionic acid,methyl acetate and Methyl propionate[J].Fluid Phase Equilibria,2008,271(l/2):69-72.

[7]Lee Ming-Jer,Lai Cheng-Hai,Warg Tzu-Bing.Vaporliquid Equilibrium of mixtures containing adipic acid,glutcaric acid,Dimethyl adipate,dimethyl glutarate,methanol,and water[J].Journal Of Chemieal&Engineering Data,2007,52(4):1291-1296.

[8]Hayden J G,O′Connell J P.A generalized method for predicting Second virial coefficient[J].Journal of Ind Eng Chem Proees Des Dev,1975,14(3):209-216.

[9]ASPEN Tech Inc.ASPEN PLUS,Release 9[M].California:ASPEN Tech Inc,1995.

[10]屈一新.化工過程數(shù)值模擬及軟件[M].北京:化學(xué)工業(yè)出版社,2006.

Process simulation for rectification tower of Vinyl Acetate based on Aspen Plus

Chen Guo,Wang Mei,Cheng Cheng,Tu Kaihui

(School of Chemical Engineering,Sichuan University,Chengdu 610065)

Using Aspen Plus software,Thermodynamic models were chosen as NRTL equation for activity coefficient and Hayden-O′Connell equation for fugacity coefficient.The vinyl acetate rectification towers were simulated by using the RadFrac rectification modules of Aspen Plus chemical industry simulation system.The simulation values were basically consistent with the actual values.By Considering the effection of feeding position,reflux ratio and the side draw of overhead to the degree of separation and energy consumption.Optimization solution was proposed as:the feed stage is 62,the reflux ratio is 32,the side draw is 37.5 t/h.Using the optimized parameters,the simulation was executed again.The result shows that the condenser's heat duty declined by 15.5%,the heat duty of reboiler declined by 16.9%,the side-drew of vinyl acetate increased 0.4%,the acetic acid of the tower bottoms increased 0.13%.

Aspen Plus;vinyl acetate;rectification tower;simulation;optimization

TQ028.3.+1，TQ018

ADOI10.3969/j.issn.1006-6829.2011.05.005

2011-07-08；

2011－07－29