孫帆，程雙，于小芳，張新平，張春雷，周興貴

?

醋酸乙酯加氫合成乙醇反應(yīng)器的模型化

孫帆1,2，程雙1，于小芳1，張新平1，張春雷1，周興貴2

(1上海華誼集團(tuán)技術(shù)研究院，上海 200241；2華東理工大學(xué)化學(xué)工程聯(lián)合國(guó)家重點(diǎn)實(shí)驗(yàn)室，上海 200237）

根據(jù)實(shí)驗(yàn)室數(shù)據(jù)對(duì)醋酸乙酯加氫制乙醇宏觀動(dòng)力學(xué)方程進(jìn)行參數(shù)估計(jì)，得到宏觀動(dòng)力學(xué)模型。根據(jù)動(dòng)力學(xué)模型，同時(shí)考慮內(nèi)擴(kuò)散的影響以及物料、熱量和動(dòng)量衡算方程建立了醋酸乙酯加氫制乙醇固定床反應(yīng)器的一維擬均相模型，并分別得到不同氫酯比、不同控制溫度、不同壓力下的多種工況的模擬結(jié)果。模擬結(jié)果與實(shí)際生產(chǎn)數(shù)據(jù)吻合。該模型可用于指導(dǎo)工業(yè)反應(yīng)器的設(shè)計(jì)和操作優(yōu)化。

醋酸乙酯；固定床；反應(yīng)器；乙醇；動(dòng)力學(xué)；一維擬均相

引 言

乙醇是重要的化工原料和燃料。關(guān)于乙醇的制備技術(shù)，目前我國(guó)主要采用糧食發(fā)酵以及乙烯水合法[1]，也有采用天然氣為原料制備的報(bào)道[2-4]，國(guó)外近幾年有較多以天然氣為原料制備乙醇的專(zhuān)利報(bào)道[5-7]。近年來(lái)，醋酸酯加氫制備乙醇技術(shù)越來(lái)越受到人們的重視，一方面可以緩解目前國(guó)內(nèi)醋酸行業(yè)處于產(chǎn)能?chē)?yán)重過(guò)剩、市場(chǎng)持續(xù)低迷的困境，開(kāi)發(fā)其下游產(chǎn)業(yè)鏈[8]，另一方面由于石油的不可再生和石油產(chǎn)區(qū)的不穩(wěn)定性，燃料能源安全問(wèn)題在全球范圍內(nèi)引起了越來(lái)越多的關(guān)注，因此，大力發(fā)展燃料乙醇、乙醇柴油等以乙醇為原料的能源產(chǎn)業(yè)已經(jīng)成為國(guó)際燃料能源產(chǎn)業(yè)的一大亮點(diǎn)[9]。與其他技術(shù)相比，采用醋酸乙酯加氫合成乙醇，催化劑價(jià)格較低、單程轉(zhuǎn)化率高、擇性好且副產(chǎn)物少，可以大大降低后續(xù)分離的難度[10]。

加氫反應(yīng)器是醋酸乙酯加氫制乙醇生產(chǎn)中的重要裝置，目前國(guó)內(nèi)外文獻(xiàn)還沒(méi)有關(guān)于該反應(yīng)器模型的相關(guān)報(bào)道，僅有關(guān)于催化劑的文獻(xiàn)[11-13]。本研究首先根據(jù)實(shí)驗(yàn)室數(shù)據(jù)建立動(dòng)力學(xué)模型，然后根據(jù)動(dòng)力學(xué)模型，同時(shí)考慮內(nèi)擴(kuò)散的影響以及物料、熱量和動(dòng)量衡算方程建立了反應(yīng)器模型，模型結(jié)果與單管實(shí)際工況較吻合。

1 宏觀動(dòng)力學(xué)模型

原料氫氣（H2）和汽化的醋酸乙酯（EA）在固定床反應(yīng)器中與銅基催化劑接觸，進(jìn)行加氫反應(yīng)，生成乙醇（EtOH）。反應(yīng)過(guò)程可以描述為

為便于工程設(shè)計(jì)應(yīng)用，宏觀動(dòng)力學(xué)方程模型一般用冪函數(shù)形式，反應(yīng)速率常數(shù)采用修正的阿倫尼烏斯方程[14]，反應(yīng)動(dòng)力學(xué)方程形式如下

采用非線性最小二乘法對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行參數(shù)估值[15]，實(shí)驗(yàn)結(jié)果和參數(shù)擬合可見(jiàn)文獻(xiàn)[16]，經(jīng)計(jì)算得到動(dòng)力學(xué)參數(shù)為

01.585 Pa-1·s-1,6.056×104J·mol-1,

2 反應(yīng)器模型

2.1 內(nèi)擴(kuò)散的影響

內(nèi)擴(kuò)散效率因子計(jì)算式為

氣相反應(yīng)物A在固相中的有效擴(kuò)散系數(shù)由式(6)～式(9)計(jì)算[17]

所用催化劑的物理參數(shù)為：粒徑0.003 m，長(zhǎng)度0.003 m，顆粒密度1.6×103kg·m-3，顆粒平均孔徑260×10-10m，顆粒孔隙率0.61，曲節(jié)因子1.35。

2.2 數(shù)學(xué)模型

建立醋酸乙酯加氫合成乙醇固定床反應(yīng)器模型時(shí)，采用一維擬均相平推流模型[18]，做如下假設(shè)：①流體在反應(yīng)器中溫度、濃度徑向均一，僅沿軸向變化；②流體與催化劑在任一與流體流動(dòng)方向垂直的橫截面處的溫度、反應(yīng)物濃度相同；③流體在床層中以平推流流動(dòng)。

在反應(yīng)器床層高度上取微元d進(jìn)行物料、熱量和動(dòng)量衡算，得到模型數(shù)學(xué)表達(dá)式如下。

物料衡算方程為

H2和EA在反應(yīng)過(guò)程中的摩爾分?jǐn)?shù)均用生成乙醇的摩爾分?jǐn)?shù)表示

熱量衡算方程為

動(dòng)量衡算采用厄根（Ergun）方程[19]

氣體混合物的混合規(guī)則采用Kay規(guī)則[20]，該規(guī)則提出的虛擬臨界參數(shù)為

混合氣體的密度、比熱容等參數(shù)也均采用此混合規(guī)則。

3 模型計(jì)算結(jié)果

在模擬計(jì)算時(shí)，采用Runge-Kutta積分算法對(duì)反應(yīng)器模型沿軸向積分求解反應(yīng)結(jié)果。在不同工況條件下運(yùn)行模型，得到模型的模擬結(jié)果，與實(shí)際工況下的乙酯轉(zhuǎn)化率、反應(yīng)器熱點(diǎn)溫度和熱點(diǎn)位置等主要指標(biāo)進(jìn)行對(duì)比。

3.1 不同氫酯比下的模擬結(jié)果

在壓力為5.5 MPa、控制溫度為204℃時(shí)，對(duì)不同氫酯比下的5組工況（氫酯比分別為7.5、10、13、20和30時(shí)）進(jìn)行模擬，乙酯轉(zhuǎn)化率模擬值與實(shí)際值的對(duì)比如圖1所示。

圖1 不同氫酯比時(shí)的乙酯轉(zhuǎn)化率

氫酯比分別為7.5、13、20和30時(shí)的床層溫度分布模擬值與實(shí)際值的對(duì)比如圖2所示。其中，模擬轉(zhuǎn)化率與實(shí)際轉(zhuǎn)化率的平均相對(duì)誤差為0.32%，模擬熱點(diǎn)溫度與實(shí)際熱點(diǎn)溫度的平均相對(duì)誤差為0.20%。從模擬結(jié)果可以看出，隨著氫酯比的增大，反應(yīng)器出口的乙酯轉(zhuǎn)化率逐漸增大，熱點(diǎn)溫度逐漸降低，熱點(diǎn)位置略微提前。

3.2 不同控制溫度下的模擬結(jié)果

在壓力為5.5 MPa、氫酯比為13～14時(shí)（模擬時(shí)氫酯比取為13.5），對(duì)不同控制溫度下的7組工況進(jìn)行模擬，乙酯轉(zhuǎn)化率模擬值與實(shí)際值的對(duì)比如圖3所示，床層溫度分布實(shí)際值和模擬結(jié)果如圖4所示。其中，模擬轉(zhuǎn)化率與實(shí)際轉(zhuǎn)化率的平均相對(duì)誤差為0.18%。從模擬結(jié)果可以看出，隨著進(jìn)口溫度的增大，反應(yīng)器出口的乙酯轉(zhuǎn)化率逐漸增大，熱點(diǎn)溫度逐漸增大，熱點(diǎn)位置略微提前。

圖3 不同進(jìn)口溫度時(shí)的乙酯轉(zhuǎn)化率

圖4 不同進(jìn)口溫度時(shí)的床層溫度分布

3.3 不同壓力下的模擬結(jié)果

在控制溫度為204℃、氫酯比為13時(shí)，對(duì)不同壓力下的3組工況進(jìn)行模擬，壓力分別為5.5、4.8、4 MPa時(shí)的乙酯轉(zhuǎn)化率模擬值與實(shí)際值的對(duì)比如圖5所示。

圖5 不同壓力下的乙酯轉(zhuǎn)化率

床層溫度分布模擬值與實(shí)際值的對(duì)比如圖6所示。其中，模擬轉(zhuǎn)化率與實(shí)際轉(zhuǎn)化率的平均相對(duì)誤差為0.34%，模擬熱點(diǎn)溫度與實(shí)際熱點(diǎn)溫度的平均相對(duì)誤差為0.47%。從模擬結(jié)果可以看出，隨著壓力的增大，反應(yīng)器出口的乙酯轉(zhuǎn)化率逐漸增大，熱點(diǎn)溫度基本不變，熱點(diǎn)位置略微提前。

圖6 不同壓力時(shí)床層溫度分布

4 結(jié) 論

（1）針對(duì)醋酸乙酯加氫制乙醇反應(yīng)體系建立了醋酸乙酯加氫合成乙醇固定床反應(yīng)器的一維擬均相模型。模擬計(jì)算結(jié)果表明，增大氫酯比、增加反應(yīng)器進(jìn)口溫度以及增加反應(yīng)器壓力均可以增大反應(yīng)器出口的乙酯轉(zhuǎn)化率，氫酯比和反應(yīng)器進(jìn)口溫度對(duì)反應(yīng)器熱點(diǎn)溫度影響較大，反應(yīng)器壓力則對(duì)熱點(diǎn)溫度影響較小。

（2）模型計(jì)算結(jié)果與單管裝置實(shí)際生產(chǎn)數(shù)據(jù)吻合度較好，偏差較小，表明所建立的反應(yīng)器模型及其他數(shù)據(jù)的計(jì)算和處理是合適的，所建模型可用于醋酸乙酯加氫合成乙醇固定床工業(yè)反應(yīng)器的模擬、設(shè)計(jì)和操作優(yōu)化。

符 號(hào) 說(shuō) 明

Am——催化劑顆粒外表面面積，m2 cp——混合氣體比熱容，J·mol-1·K-1 cpc——導(dǎo)熱油比熱容，J·mol-1·K-1 DA,Knu——Knudsen擴(kuò)散系數(shù)，m2·s-1 DA,mol——分子擴(kuò)散系數(shù)，m2·s-1 DA,pore——孔擴(kuò)散系數(shù)，m2·s-1 Deff,A——?dú)怏w反應(yīng)物A在固相的有效擴(kuò)散系數(shù)，m2·s-1 dpore——顆粒平均孔徑，m

ds——顆粒當(dāng)量直徑，m E——活化能，J·mol-1 F——混合氣體流量，kg·m-2·s-1 Fc——導(dǎo)熱油流量，kg·m-2·s-1 H——反應(yīng)熱，J·mol-1 K——管程傳熱系數(shù)，J·m-2·K-1·s-1 k——反應(yīng)速率常數(shù)，Pa-1·s-1 k0——指前因子，Pa-1·s-1 l——床層位置，m M——組分相對(duì)分子質(zhì)量 p——?dú)怏w總壓，Pa pci——組分i的臨界壓力，Pa pcm——虛擬臨界壓力，Pa pEA——EA的分壓，Pa R——?dú)怏w常數(shù)，J·mol-1·K-1 Rem——Reynolds數(shù) S——反應(yīng)管內(nèi)截面積，m2

T——溫度，K Tc——導(dǎo)熱油溫度，K Tci——組分i的臨界溫度，K Tcm——虛擬臨界溫度，K t——反應(yīng)時(shí)間，s um——混合氣體質(zhì)量流速，kg·m-2·s-1 V——分子擴(kuò)散體積 yi——組分i的摩爾分?jǐn)?shù) α——H2的反應(yīng)級(jí)數(shù) β——EA的反應(yīng)級(jí)數(shù) εB——床層孔隙率 εp——催化劑顆?？紫堵?η——內(nèi)擴(kuò)散效率因子 ρB——催化劑床層密度，kg·m-3 ρg——混合氣體密度，kg·m-3 τp——曲節(jié)因子 ——Thiele模數(shù)

References

[1] Li Fang, Zhang Ke, Ma Hongfang, Zhang Haitao, Ying Weiyong, Fang Dingye. Thermodynamic analysis of different routes of ethanol synthesis [J]., 2013, 30 (3): 256-260

[2] Sun Lujing (孫陸晶). Alcohol fuel industry prospect [J].(化學(xué)工業(yè)), 2011, 29 (1): 8-12

[3] Wang Peng (王鵬)，Wang Xiangui (王憲貴)，Guo Zhanying (郭戰(zhàn)英). Research progress in producer gas to ethanol technology [J].(潔凈煤技術(shù)), 2010, 16 (1): 55-58

[4] Phillips J, Clausen E, Gaddy J. Biological production of ethanol from coal synthesis gas [J]., 1993, 40 (1): 559-571

[5] Johnston V, Chapman J, Chen L. Ethanol production from acetic acid utilizing a cobalt catalyst [P]: US, 7608744. 2009

[6] Wollrab R, Johnston V, Warner R. Process for producing ethanol by hydrogenation of acetic acid [P]: WO, 2011140468. 2011

[7] Johnston V, Kimmich B, Waal J. Ethanol production from acetic acid utilizing a molybdenum carbide catalyst [P]: US, 8350098. 2013

[8] Hu Yurong (胡玉容), Wang Ke (王科), Li Yang (李揚(yáng)), Fan Xin (范鑫), Yuan Xiaojin (袁小金), Xu Hongyun (許紅云). Simulation of synthesizing ethanol by hydrogenation of ethyl acetate [J].(化學(xué)工程師), 2011 (9): 25-27

[9] Yang Jie (楊潔), Ding Mingzhu (丁明珠), Li Bingzhi (李炳志), Yuan Yingjin (元英進(jìn)). Phospholipidome alterations ofin bioethanol fermentation process [J].(化工學(xué)報(bào)), 2012, 63 (6):1830-1835

[10] Wang Biao (王彪), Wang Xiting (王熙庭), Xu Guohui (徐國(guó)輝). Advances in technologies for hydrogenation of acetic acid and acetates to ethanol [J].(天然氣化工), 2013, 38 (3): 79-83

[11] Song Qingfeng (宋慶鋒)，Zhang Yong (張勇)，Zeng Qinghu (曾清湖). Techno-economic analysis of production process of syngas to ethanol [J].(工業(yè)催化), 2013, 21 (6): 17-21

[12] Yang Yong (楊勇), Wang Yangdong (王仰東), Liu Su (劉蘇). Mo-Co-K sulfide-based catalysts promoted by rare earth salts for selective synthesis of ethanol and mixed alcohols from syngas [J].(催化學(xué)報(bào)), 2007, 28 (12): 1028 -1030

[13] Pan Hui (潘慧), Bai Fenghua (白鳳華), Su Haiquan (蘇海泉). Research progress of Rh-based catalyst for conversion of syngas to ethanol [J].(化工進(jìn)展), 2010, 29 (2): 157-161

[14] Xia Xinyu (夏新宇), Zhang Wenzheng (張文正). A kinetic study on closed-system petroleum generation simulations under constant temperature and its significance [J].(石油勘探與開(kāi)發(fā)), 2000, 27 (5): 22-26

[15] Roger F, Sven L. Solving mathematical programs with complementarity constraints as nonlinear programs [J]., 2004, 19 (1): 15-40

[16] Yu Xiaofang (于小芳), Cheng Shuang (程雙), Sun Fan (孫帆), Zhang Xinping (張新平), Zhang Chunlei (張春雷), Tang Yi (唐頤). The macro-kinetics of ethyl acetate hydrogenation to ethanol [J].:(復(fù)旦學(xué)報(bào):自然科學(xué)版), 2014, 53 (1):77-83

[17] Zhang Xinping (張新平), Sui Zhijun (隋志軍), Zhou Xinggui (周興貴). Modeling of coke burning off regeneration for coked propane dehydrogenation catalyst [J].(化工學(xué)報(bào)), 2009, 60 (1): 163-167

[18] Ling Zeji (凌澤濟(jì)). Mathematical simulation of synthetic reactor of ethylene oxide [J].(江蘇化工), 2005, 33 (5): 54-57

[19] Xiong Wei (熊瑋), Bi Xuegong (畢學(xué)工), Zhou Guofan (周?chē)?guó)凡). Simulation study on gas pressure drop characteristics in the lower zone of blast furnace [J].(高校化學(xué)工程學(xué)報(bào)), 2007, 21 (2): 233-238

[20] An Chao (安超), Guo Xiao (郭肖), Zhang Yong (張勇), Li Gang (李剛), Ma Jun (馬軍). Calculation methods optimization of deviation factor for overpressure gas reservoir [J].:(重慶科技學(xué)院學(xué)報(bào): 自然科學(xué)版), 2010, 12 (5): 92-94

Modeling of ethyl acetate hydrogenation reactor

SUN Fan1,2, CHENG Shuang1, YU Xiaofang1, ZHANG Xinping1, ZHANG Chunlei1, ZHOU Xinggui2

(Technology Research Institute of Shanghai Huayi GroupShanghaiChinaState Key Laboratory of Chemical EngineeringEast China University of Science and TechnologyShanghaiChina

A kinetic model of synthesizing ethanol by hydrogenation of ethyl acetate was developed based on laboratory data, and then corrected by taking into consideration the pore diffusion effect. With the correction, the results could conform with actual operation. Based on the principles of material balance, energy balance, momentum balance and kinetics of the reactions, an one-dimensional pseudo-homogeneous model of the ethanol synthesis reactor by hydrogenation of ethyl acetate using fixed bed was established. The simulation results of a variety of conditions were obtained by the model with different H2/EA molar ratios, different temperatures and different pressures. The model could reflect the performance of the hydrogenation reactor and provide guidance for the design and operation optimization of industrial reactor.

ethyl acetate; fixed bed; reactor; ethanol; kinetics; one-dimensional pseudo-homogeneous

2014-08-11.

10.11949/j.issn.0438-1157.20141217

TQ 021.8

A

0438—1157（2015）02—0561—06

2014-08-11收到初稿，2014-09-30收到修改稿。

聯(lián)系人：張春雷，周興貴。第一作者：孫帆（1985—），男，博士。

Prof.ZHANG Chunlei, clzhang@shhuayitec.com; Prof. ZHOU Xinggui, xgzhou@ecust.edu.cn