周彩榮， 王青云

(鄭州大學(xué) 化工與能源學(xué)院 河南 鄭州 450001)

DOI: 10.13705/j.issn.1671-6841.2017089

己二酸二甲酯反應(yīng)精餾塔的模擬

周彩榮， 王青云

(鄭州大學(xué) 化工與能源學(xué)院 河南 鄭州 450001)

使用Aspen Plus軟件，對以己二酸和甲醇為原料，對甲苯磺酸為催化劑制備己二酸二甲酯的反應(yīng)精餾過程進行了模擬.以己二酸二甲酯的質(zhì)量分數(shù)為考察目標，通過Aspen Plus軟件中的設(shè)計規(guī)定、靈敏度分析、塔板和填料的設(shè)計與校核，確定優(yōu)化結(jié)果為：己二酸的轉(zhuǎn)化率為98.57%,塔釜內(nèi)己二酸二甲酯的質(zhì)量分數(shù)為77.4%，流率為66.26 kmol/h；理論塔板數(shù)為25，塔徑為0.535 m; 己二酸和甲醇分別由12級和14級塔板處進料，進料流率比和回流比分別為0.355和2.5；冷凝器和再沸器熱負荷分別為-1 153.79 kW和1 248.93 kW；塔內(nèi)最大液泛因子為0.072，塔段壓降為9.103 kPa，最大降液管液位/板間距為0.325.

己二酸二甲酯； 反應(yīng)精餾； 模擬

DOI: 10.13705/j.issn.1671-6841.2017089

0 引言

有機酸和有機醇進行的化學(xué)反應(yīng)屬于復(fù)雜的反應(yīng)體系，因酯化反應(yīng)與水解反應(yīng)同時存在，屬于可逆的化學(xué)平衡反應(yīng).要提高酯化反應(yīng)物的收率,主要的方法是蒸出酯或水，使反應(yīng)朝著酯化反應(yīng)的方向進行.目前常用的提高酯化率的方法[1-5]有：在體系中添加帶水劑、反應(yīng)-吸附耦合法、反應(yīng)-精餾耦合法.己二酸二甲酯(DMA)為無色清澈透明液體，沸點112 ℃，易溶于醇、醚類，不溶于水，主要用于溶劑、醫(yī)藥、造紙和纖維素的生產(chǎn).此外，它還是合成1,6-己二醇的中間體.己二酸二甲酯是重要的精細化工產(chǎn)品和化工原料，工業(yè)上比較成熟的合成方法是采用己二酸與甲醇在硫酸或硝酸催化作用下制備,但是由于硫酸的加入對合成設(shè)備具有較大的腐蝕，且不能滿足環(huán)保的要求[6].本文以對甲苯磺酸為催化劑，以產(chǎn)品中己二酸二甲酯的質(zhì)量分數(shù)為考察目標，用Aspen Plus軟件[7-8]中的設(shè)計規(guī)定、靈敏度分析、塔板和填料的設(shè)計與校核，確定該工藝的工藝參數(shù)，為反應(yīng)精餾耦合工藝制備己二酸二甲酯的工業(yè)化設(shè)計和實際操作提供借鑒.

A：甲醇進料口；B：己二酸+催化劑進料口； C：塔頂出料口； D：塔底產(chǎn)品出料口圖1 己二酸二甲酯的工藝流程Fig.1 Flowsheet of dimethyl adipate

1 反應(yīng)精餾合成己二酸二甲酯的基本工程數(shù)據(jù)

1.1 工藝流程

己二酸二甲酯的工藝流程如圖1所示.己二酸由反應(yīng)精餾塔上部加入，對甲苯磺酸由中上部加入，甲醇由反應(yīng)精餾塔下部加入.在反應(yīng)段內(nèi)，下降的己二酸物流與上升的甲醇物流逆流接觸，在催化劑的作用下發(fā)生反應(yīng).己二酸二甲酯和水由塔釜流出，未完全反應(yīng)的甲醇由塔頂蒸出.圖中精餾段主要是甲醇、己二酸二甲酯和己二酸的分離.在反應(yīng)段，己二酸和過量的甲醇在對甲苯磺酸的催化作用下開始反應(yīng)生成己二酸二甲酯和水，未反應(yīng)完的甲醇上升,己二酸二甲酯和水下降.提餾段主要是實現(xiàn)甲醇、己二酸二甲酯與水的分離.

表1 標準狀態(tài)下液態(tài)生成焓和標準熵

1.2 熱力學(xué)和動力學(xué)數(shù)據(jù)

反應(yīng)精餾塔的模擬與設(shè)計[9-11]所需的熱力學(xué)和動力學(xué)數(shù)據(jù)見表1[12]、表2 和表3[13].在50 ℃條件下合成己二酸二甲酯的正、逆反應(yīng)均服從二級動力學(xué)模型，正、逆反應(yīng)的表觀活化能分別為6.2×104J/mol、4.44×104J/mol，指前因子分別為1.382×108、1.003×105.

表2 不同溫度下反應(yīng)的焓變、Gibbs自由能變 和標準平衡常數(shù)值

表3 反應(yīng)速率常數(shù)隨溫度的變化

1.3 理論塔板數(shù)

采用估算法可快速求出反應(yīng)精餾塔的理論塔板數(shù)，具體步驟見文獻[14]. 設(shè)體系餾出物組分XD=0.99，釜液組分XW=0.01，求取雙組分溶液最少理論塔板數(shù)的芬斯克方程為

平均相對揮發(fā)度為

(2)

溫度方程式為

(3)

結(jié)合(1)～(3)得Nm=14.5，而真正的塔板數(shù)通常約為最少理論塔板數(shù)的2倍，所以本設(shè)計得出整個催化精餾塔的塔板數(shù)為29.

2 Aspen Plus軟件模擬與優(yōu)化

2.1 物性估算

在輸入數(shù)據(jù)步驟中，對己二酸二甲酯分別進行了物性查詢和物性估算，當(dāng)將所得結(jié)果帶入模擬中時，發(fā)現(xiàn)物性查詢存在參數(shù)缺失，物性估算滿足模擬要求，己二酸二甲酯的物性估算結(jié)果為：相對分子質(zhì)量174.197，相對密度1.067，常壓沸點230.73 ℃，凝固點10.29 ℃，偏心因子0.606，標準液體摩爾體積0.163 m3/kmol，標準生成焓-8.198×108J/mol，臨界值(溫度418.85 ℃，摩爾體積0.549 m3/kmol，壓縮因子0.24).

2.2 常壓己二酸二甲酯-甲醇-水三元體系相平衡

在設(shè)計反應(yīng)精餾塔時，要求反應(yīng)物甲醇大大過量，在反應(yīng)段中反應(yīng)物己二酸要全部轉(zhuǎn)化，因而在精餾段和提餾段涉及的是水、甲醇和己二酸二甲酯的三元體系.用NRTL物性方法對水-甲醇-己二酸二甲酯三元體系進行三元液液相平衡回歸分析，進行回歸分析時，需用到水-甲醇-己二酸二甲酯三元體系液液相平衡數(shù)據(jù)[15].根據(jù)文獻[15]數(shù)據(jù)，由Aspen Plus軟件模擬得到NRTL方程中的二元交互作用參數(shù)，結(jié)果見表4.

表4 NRTL方程下水、甲醇、己二酸二甲酯的二元交互作用參數(shù)

2.3 工藝模擬與優(yōu)化結(jié)果

2.3.1輸入?yún)?shù) 設(shè)計中有己二酸和甲醇兩股物流，兩者的物流參數(shù)如下：對于己二酸物流，進料溫度45 ℃，壓力2.026×105Pa，流量20 kmol/h；對于甲醇物流，進料溫度45 ℃，壓力2.026×105Pa，流量80 kmol/h.

因本設(shè)計為反應(yīng)精餾塔，在塔內(nèi)所發(fā)生的反應(yīng)為

C6H10O4+2CH3OH?C8H14O4+2H2O.

正反應(yīng)的動力學(xué)方程為

(4)

逆反應(yīng)的動力學(xué)方程為

(5)

在Aspen Plus軟件中，對于分離單元模擬，設(shè)計中選用RadFrac模塊，輸入?yún)?shù)為：理論塔板數(shù)29；全凝器為釜式再沸器；回流比2.5；塔頂產(chǎn)品與進料流率比0.337 5；己二酸進料位置是15級塔板上方進料(Above-stage)；甲醇進料位置是15級塔板上進料(On-stage)；塔內(nèi)每塊板上的液相持液量為0.3 L；再沸器液相持液量為1.0 L.

2.3.2運行模擬結(jié)果 Aspen Plus軟件的運行模擬結(jié)果如表5和表6所示.從表中數(shù)據(jù)可知，己二酸的轉(zhuǎn)化率為98.57%，從塔釜得到產(chǎn)品己二酸二甲酯，其質(zhì)量分數(shù)為77.4%.

表5 組分模擬結(jié)果

冷凝器與再沸器模擬結(jié)果見表7，其中塔釜產(chǎn)品流率為66.26 kmol/h，冷凝器熱負荷為-1 153.79 kW，再沸器熱負荷為1 248.93 kW.反應(yīng)精餾塔內(nèi)溫度分布曲線和液相質(zhì)量組成分布曲線如圖2和圖3所示，可以看出，在反應(yīng)精餾塔內(nèi)，精餾段塔內(nèi)溫度基本在64.5 ℃左右，在反應(yīng)段溫度變化幅度較大，溫度從上到下逐漸增加，在塔釜得到產(chǎn)品己二酸二甲酯.

3 靈敏度分析

設(shè)計中分別從以下幾個方面進行了靈敏度分析，以尋求對工藝條件的優(yōu)化[16].在此分析過程中，除了所考察因素的數(shù)值發(fā)生變化，其他所需因素的數(shù)值均固定不變，皆為2.3.1中所給數(shù)值.

表6 反應(yīng)精餾塔模擬結(jié)果

表7 冷凝器與再沸器模擬結(jié)果

圖2 反應(yīng)精餾塔內(nèi)溫度分布曲線Fig.2 Temperature profile in the reactive distillation column

圖3 反應(yīng)精餾塔內(nèi)液相質(zhì)量組成分布曲線Fig.3 Liquid composition profiles in the reactive distillation column

圖4 產(chǎn)品質(zhì)量分數(shù)隨進料流率比的變化曲線Fig.4 Product mass fraction with the change of flow ratio

3.1 產(chǎn)品質(zhì)量分數(shù)隨進料流率比的變化

圖4是產(chǎn)品質(zhì)量分數(shù)隨進料流率比的變化曲線.可以看出，進料流率比從0.335起開始逐漸增加，直至增加到0.354 8和0.354 9后開始逐漸下滑，所以最佳進料流率比選擇0.355.塔頂采出量的多少直接影響產(chǎn)品純度和裝置的產(chǎn)能，在產(chǎn)品純度滿足要求時應(yīng)盡量提高采出量以提高裝置的產(chǎn)能.

3.2 最小回流比隨理論塔板數(shù)的變化

理論塔板數(shù)增加，回流比會減少.當(dāng)回流比隨理論塔板數(shù)的增加而基本不變時，可認為回流比的值為最小回流比.在理論塔板數(shù)增大的過程中，可以假設(shè)進料位置與理論塔板數(shù)的比值是個定值.已知理論塔板數(shù)為29，最佳進料位置為15，所以此定值可取15/29=0.52.因進料位置隨理論塔板數(shù)而變，故進行靈敏度分析的同時需添加計算器(Calculator).計算器模塊也稱Fortran模塊，在此模塊中，可自行編寫Aspen Plus軟件可執(zhí)行的Fortran程序，將語句插入到流程計算中，運行模擬結(jié)果見表8.可以看出，隨著理論塔板數(shù)的增加，回流比不發(fā)生變化，所以本設(shè)計選擇回流比2.5為宜.

表8 回流比模擬結(jié)果

3.3原料進料位置對產(chǎn)品質(zhì)量分數(shù)的影響

原料己二酸、甲醇進料位置對產(chǎn)品質(zhì)量分數(shù)的影響結(jié)果如圖5、圖6所示.在一定的流率比(0.355)和回流比(2.5)時，產(chǎn)品己二酸二甲酯的質(zhì)量分數(shù)隨甲醇進料位置的下移而降低，當(dāng)進料位置大于14時，產(chǎn)品質(zhì)量分數(shù)又隨甲醇進料位置的下移而升高，但考慮到甲醇進料位置下移會使得再沸器負荷增加，因而選取甲醇進料位置在14級塔板處.由于己二酸的進料位置受甲醇進料位置的約束，工藝上要求己二酸進料位置要高于甲醇.因此，己二酸進料位置選取12級塔板處，即可以達到設(shè)計要求.

圖5 產(chǎn)品質(zhì)量分數(shù)隨甲醇進料位置的變化曲線Fig.5 Product mass fraction with the change of the feed-stage of CH3OH

圖6 產(chǎn)品質(zhì)量分數(shù)隨己二酸進料位置的變化曲線Fig.6 Product mass fraction with the change of the feed-stage of C6H10O4

3.4 塔板數(shù)對產(chǎn)品質(zhì)量分數(shù)和溫度的影響

塔板數(shù)對產(chǎn)品質(zhì)量分數(shù)和溫度的影響結(jié)果如圖7、圖8所示.由圖可知，當(dāng)塔板數(shù)大于25時，產(chǎn)品質(zhì)量分數(shù)和塔釜溫度基本上不發(fā)生變化，因而反應(yīng)精餾塔的塔板數(shù)不能少于25.

圖7 產(chǎn)品質(zhì)量分數(shù)隨塔板數(shù)的變化曲線Fig.7 Product mass fraction with the change of the number of stages

圖8 塔釜溫度隨塔板數(shù)的變化曲線Fig.8 Temperature of bottoms with the change of the number of stages

3.5 塔板和填料的設(shè)計與校核

為了確定反應(yīng)精餾塔中的最大液泛因子、塔段壓降、最大降液管液位/板間距，設(shè)計中需進行塔板和填料的設(shè)計與校核.其中設(shè)計結(jié)果要滿足最大液泛因子小于0.8；最大降液管液位/板間距為0.2～0.5.根據(jù)設(shè)計條件進行運行模擬的結(jié)果可知，本設(shè)計的最大液泛因子為0.072；塔段壓降為9.103 kPa；最大降液管液位/板間距為0.325. 故經(jīng)校核設(shè)計合格，按此條件設(shè)計的反應(yīng)精餾塔的塔徑為0.535 m.

4 結(jié)論

對己二酸二甲酯反應(yīng)精餾塔進行了設(shè)計，用Aspen Plus軟件對工藝過程進行模擬計算和驗證，并用NRTL方程預(yù)測水-甲醇-己二酸二甲酯三元體系的液液相平衡數(shù)據(jù).用靈敏度分析法對反應(yīng)精餾塔進行模擬優(yōu)化計算，得到理論塔板數(shù)為25，塔徑為0.535 m；甲醇進料位置14級塔板，己二酸進料位置12級塔板；回流比為2.5；進料流率比為0.355；塔釜產(chǎn)品流率為66.26 kmol/h，冷凝器熱負荷為-1 153.79 kW，再沸器熱負荷為1 248.93 kW；最大液泛因子為0.072，塔段壓降為9.103 kPa，最大降液管液位/板間距為0.325.相應(yīng)的己二酸的轉(zhuǎn)化率為98.57%，產(chǎn)品己二酸二甲酯的質(zhì)量分數(shù)可達到77.4%(塔底).塔底液相經(jīng)液液相分離器，除去水和未反應(yīng)的己二酸，得到純產(chǎn)品己二酸二甲酯，塔頂餾出液甲醇可循環(huán)使用.

[1] 肖傳豪，周彩榮. 己二酸二甲酯反應(yīng)-吸附耦合技術(shù)的動力學(xué)研究[J].化工時刊, 2011，25(5):8-12.

[2] DOHERTY M F, MALONE M F. Conceptual design of distillation systems[M]. New York: McGraw-Hill,2001.

[3] 金柳偉. 反應(yīng)精餾技術(shù)制備醋酸甲酯工藝的研究[D].浙江:浙江大學(xué)，2006.

[4] SMITH L A. Catalyst system for separating isobutene from C4streams： US 4215011A[P].1980-01-25.

[5] 馬敬環(huán),劉家祺,李俊臺，等.反應(yīng)精餾技術(shù)的進展[J].化學(xué)反應(yīng)工程與工藝, 2003,19(1):1-8.

[6] 王海峰，周彩榮. 合成己二酸二甲酯的反應(yīng)-精餾耦合技術(shù)[J].天然氣化工(C1化學(xué)與化工)，2015，40(2)：21-24.

[7] 孫蘭義.化工流程模擬實訓(xùn):Aspen Plus教程[M].北京：化學(xué)工業(yè)出版社,2012.

[8] LUYBEN W L. Distillation design and control using Aspen simulation[M]. Hoboken: John Wiley amp; Sons, 2013.

[9] 安維中,林子昕，江月，等.考慮內(nèi)部熱集成的乙二醇反應(yīng)精餾系統(tǒng)設(shè)計與優(yōu)化[J]. 化工學(xué)報, 2013, 64(12):4634-4640.

[10] 項漢銀.反應(yīng)精餾合成乙酸乙酯的實驗研究與模擬[J].計算機與應(yīng)用化學(xué)，2009,26(6):789-793.

[11] 楊萬典，仲慶，王佳兵. 丁二酸二甲酯精制工藝的設(shè)計與模擬分析[J]. 化學(xué)工程，2015，43(12)：73-78.

[12] 楊超，魏曉霞，霍穩(wěn)周. 己二酸二甲酯合成體系的熱力學(xué)估算及分析[J].當(dāng)代化工,2011,40(12):1305-1308.

[13] 張曉娟,唐麗華,賈長英，等. 己二酸二甲酯合成反應(yīng)動力學(xué)研究[J].沈陽工業(yè)大學(xué)學(xué)報,2008,30(3):350-355.

[14] 秦玉楠. 沸點估算法求取理論塔板數(shù)[J]. 化工設(shè)計, 1990(2):24-25.

[15] HUNG S B， LIN H M， YU C C，et al. Liquid-liquid equilibria of aqueous mixtures containing selected dibasic esters and/or methanol[J]. Fluid phase equilibria, 2006, 248 (2): 174-180.

[16] RESTREPO J B，OLIVAR G，CARDONA C A. Bifurcation analysis of dynamic process models using Aspen Dynamics and Aspen Custom Modeler[J].Computers and chemical engineering, 2014, 62:10-20.

(責(zé)任編輯：孔 薇)

SimulationonDimethylAdipateReactionDistillationColumn

ZHOU Cairong, WANG Qingyun

(SchoolofChemicalEngineeringandEnergy,ZhengzhouUniversity,Zhengzhou450001,China)

Adipic acid and methanol as raw materials andp-toluenesulfonic acid as a catalyst, dimethyl adipate reaction distillation column design was optimized by using Aspen Plus software. With the mass fraction of dimethyl adipate as an investigation target, the sensitivity analysis and the design of the tray and packing were examined. The simulation results were as follows. The conversion rate of adipic acid was 98.57%, the mass fraction of dimethyl adipate was 77.4% at 66.26 kmol/h in the bottoms; theoretical number of stages was 25 with column diameter 0.535 m; the feed positions of adipic acid and methanol were the stages of 12th and 14th; the flow ratio was 0.355 at the reflux ratio 2.5; heat duty were -1 153.792 kW and 1 248.93 kW in the condenser and reboiler, respectively; the maximum flooding factor was 0.072, the section pressure drop was 9.103 kPa, and the maximum backup/tray spacing was 0.325.

dimethyl adipate; reactive distillation; simulation

2017-03-02

河南省教育廳科技攻關(guān)重點項目(13A530712).

周彩榮(1958—)，女，江蘇沭陽人，教授，主要從事精細有機化工、分離技術(shù)和工程基礎(chǔ)數(shù)據(jù)研究，E-mail：zhoucairong@zzu.edu.cn.

TQ053.5

A

1671-6841(2017)04-0087-06