王 菊，鐘思青，張成芳，金照生

1.中國石油化工股份有限公司上海石油化工研究院，上海 201208；2.華東理工大學(xué)大型工業(yè)反應(yīng)器工程，教育部工程研究中心，上海 200237

50年代前絕大部分的乙烯是由乙醇脫水生產(chǎn)而來的。隨著石油化工的發(fā)展，大量廉價石油促使烴類熱裂解生產(chǎn)乙烯，以其經(jīng)濟上的優(yōu)勢取代了乙醇脫水法生產(chǎn)乙烯。但隨著生物技術(shù)的快速發(fā)展以及非糧食基生物質(zhì)乙醇制備技術(shù)的突破，使得乙醇來源日趨廣泛，乙醇法制乙烯工藝又得到了重視。

20世紀對乙醇脫水反應(yīng)的機理已進行了大量的研究[1-11]，目前關(guān)于乙醇制乙烯的研究主要集中在催化劑的研制[12]和工藝研究[13]兩大方面，而針對特定催化劑上反應(yīng)動力學(xué)研究又是工藝研究的基礎(chǔ)。至今得到廣泛應(yīng)用的是Hassan 提出了乙醇在固體酸堿催化劑上反應(yīng)的催化機理[14]，認為乙醇首先吸附在催化劑表面的酸堿中心上并形成吸附態(tài)化合物，然后吸附態(tài)中間產(chǎn)物脫水生成最終產(chǎn)物并恢復(fù)酸堿中心。由于所用催化劑的物理性質(zhì)和化學(xué)性質(zhì)的不同，對機理和動力學(xué)研究就很難達成一致，主要存在兩類爭議：（1）乙醇脫水機理是連串反應(yīng)機理、平行反應(yīng)機理和平行-連串反應(yīng)機理；（2）反應(yīng)機理模型是吸附于催化劑上的兩個鄰近乙醇分子之間進行的L-H 模型和一個氣相乙醇分子和一個催化劑上處于吸附態(tài)的乙醇分子之間進行的R-E 模型。本工作采用自研發(fā)的氧化鋁作為催化劑，進行乙醇脫水反應(yīng)性能的研究，依據(jù)平行-連串反應(yīng)機理建立包括副產(chǎn)物C4 的所選催化劑反應(yīng)體系的本征動力學(xué)模型，為催化劑的工業(yè)應(yīng)用和反應(yīng)裝置的工業(yè)放大提供數(shù)據(jù)參考。

1 實驗部分

1.1 實驗裝置

實驗裝置如圖1所示，采用?10 mm 不銹鋼-管制成的固定床積分反應(yīng)器，內(nèi)插?1 mm 的熱電偶套管。反應(yīng)器和預(yù)熱器由電熱絲加熱，溫度控制誤差不超過±1 ℃。預(yù)熱爐出口處至反應(yīng)器進口處和反應(yīng)器出口處至氣相色譜進口處流路用加熱帶保溫防止蒸汽冷凝。分析純的乙醇按照一定的水醇比（質(zhì)量比）混合后，通過SZB-2 雙柱塞微量計量泵進料?？瞻讓嶒灠l(fā)現(xiàn)，反應(yīng)器本身及里面裝填的惰性物料石英砂對反應(yīng)沒有催化作用。

1.2 內(nèi)外擴散的排除

為排除外擴散影響，采用不同催化劑的使用量，改變空速進行實驗。實驗表明，當反應(yīng)在常壓操作，反應(yīng)溫度為340 ℃的條件反應(yīng)下，空速大于0.1 mL/(h·g)，外擴散影響不存在。

內(nèi)擴散排除實驗：保持空速恒定為0.5 mL/(h·g)，常壓操作，反應(yīng)溫度為340 ℃的條件下，改變氧化鋁催化劑的顆粒度，測得乙醇轉(zhuǎn)化率隨催化劑粒度變化的曲線，實驗發(fā)現(xiàn)當催化劑粒度大于60 目后，乙醇轉(zhuǎn)化率受氧化鋁催化劑的粒度影響很小。實驗選用粒度為60～80 目的氧化鋁催化劑，稱取2 g 放入固定床反應(yīng)器管中。通入氮氣550 ℃活化2 h，最后降至所需反應(yīng)溫度。

圖1 實驗流程示意Fig.1 Experiment flow scheme

2 結(jié)果與討論

2.1 動力學(xué)模型建立

根據(jù)對所選催化劑上的乙醇脫水制乙烯反應(yīng)性能研究，乙醇脫水制乙烯可能存在下述的反應(yīng)：

結(jié)合實驗數(shù)據(jù)分析和反應(yīng)體系的熱力學(xué)分析[15]，選擇反應(yīng)（1）、（2）、（3）和（4）為乙醇脫水反應(yīng)體系的四個主要反應(yīng)，所選催化劑上的乙醇脫水反應(yīng)呈現(xiàn)的規(guī)律更應(yīng)符合平行-連串機理。

實驗發(fā)現(xiàn)，系統(tǒng)中的水分壓對反應(yīng)是有影響的，為簡化系統(tǒng)和優(yōu)化計算參數(shù)速率，這里假設(shè)這些反應(yīng)都是不可逆反應(yīng)。采用冪指數(shù)速率方程形式，建立了如下半理論半經(jīng)驗的速率方程式：

各主要反應(yīng)物及產(chǎn)物的凈生成速率如下：

同時反應(yīng)體系中各物質(zhì)的速率也可表示為單位質(zhì)量催化劑上某一物質(zhì)的摩爾流量的變化：

式中ri為反應(yīng)速率，mol/(h·g)；Pi為分壓，MPa；Fi為摩爾流率，mol/h；為進口乙醇原料的總摩爾流量，mol/h；W為催化劑裝填量，g；XA為乙醇轉(zhuǎn)化率；Yi為收率。

聯(lián)立上述方程，可優(yōu)化求解出各反應(yīng)速率常數(shù)k1～4和反應(yīng)級數(shù)n1～8。

2.2 模型參數(shù)計算

反應(yīng)級數(shù)的初值按照一定的經(jīng)驗選取，四個反應(yīng)的速率常數(shù)的初值參考文獻[16]。利用MATLAB軟件，采用四五階Runge-Kutta 法對上述微分速率方程組求解，可得到不同反應(yīng)條件下的乙醇轉(zhuǎn)化率、乙烯收率、乙醚收率和C4 收率。同時以G為目標函數(shù)，采用Marquardt 算法對參數(shù)進行優(yōu)化，可求得目標函數(shù)達到最小值時的相應(yīng)的動力學(xué)參數(shù)。

通過對參數(shù)初值的不斷調(diào)整，全局優(yōu)化求出合適的反應(yīng)級數(shù)。然后把得到的反應(yīng)級數(shù)作為已知參數(shù)，對每個溫度下的模型優(yōu)化求出不同溫度下的反應(yīng)速率常數(shù)，結(jié)果見表1。

表1 不同溫度下的反應(yīng)速率常數(shù)Table 1 Rate constants at different temperatures

按照Arrhenius 方程：k=k0exp(-Ea/RT)，對所得的速率常數(shù)以lnk與T-1作圖，如圖2所示。通過對圖2中的數(shù)據(jù)點進行線性最小二乘法擬合，可得到各個反應(yīng)的表觀指前因子k0和活化能Ea。

對所得到的模型中的指前因子、活化能和反應(yīng)級數(shù)的小數(shù)點的個數(shù)進行修正，最終得到乙醇在所選氧化鋁催化劑上脫水反應(yīng)的本征動力學(xué)模型，如下所示：

圖2 lnk 與T-1 直線擬合Fig.2 Line polyfit scheme of lnk vs T-1

2.3 動力學(xué)模型的統(tǒng)計檢驗

將上述優(yōu)化所得的動力學(xué)參數(shù)代入相應(yīng)的速率方程，可計算不同反應(yīng)條件下的關(guān)鍵組分乙醇的轉(zhuǎn)化率和乙烯收率，計算值與實驗值比較見圖3。

圖3 乙醇轉(zhuǎn)化率及乙烯收率的殘差分析Fig.3 Residual analysis of conversion of ethanol and yield of ethylene

從圖3可以看出，除少數(shù)幾個點的誤差比較大，大多數(shù)都均勻分布在對角線附近，平均誤差在10%以內(nèi)，說明此模型及參數(shù)所得到計算值與實驗值吻合。

根據(jù)統(tǒng)計理論，對模型的適定性進行了檢驗，檢驗結(jié)果見表2。

表2 動力學(xué)模型的統(tǒng)計檢驗Table 2 Statistical test of kinetic model

統(tǒng)計結(jié)果表明：ρ2>0.9，F(xiàn)>>10FT。所以模型在該水平上是合適的。

3 結(jié) 論

對選用的氧化鋁催化劑上乙醇催化脫水制乙烯的反應(yīng)本征動力學(xué)特性進行了較為全面的實驗研究，建立了描述反應(yīng)生成物的半經(jīng)驗冪指數(shù)型的本征動力學(xué)模型。通過對乙醇催化脫水制乙烯的反應(yīng)動力學(xué)的研究，得到了以下的結(jié)論。

a）乙醇催化脫水制乙烯的反應(yīng)體系是復(fù)雜的復(fù)合反應(yīng)體系，總體上屬于平行反應(yīng)體系帶有連串生成乙醚的反應(yīng)體系，即乙醇分子內(nèi)直接脫水生成乙烯，同時發(fā)生乙醇分子間脫水生成乙醚，再由乙醚發(fā)生歧化反應(yīng)生成乙烯。

b）生成乙烯的反應(yīng)的活化能Ea1（154 kJ/mol）和Ea3（162.6 kJ/mol）高于生成乙醚（142.9 kJ/mol）和C4（109.9 kJ/mol）的活化能，說明要想提高目的產(chǎn)物乙烯收率，該催化劑下宜采用適當提高反應(yīng)溫度。

c）乙醇生成乙烯的反應(yīng)級數(shù)小于生成混合C4 和乙醚的反應(yīng)級數(shù)，說明生成乙烯反應(yīng)的速率受乙醇濃度的影響相對于副產(chǎn)物乙醚和C4 要小，乙醇濃度低更有利于乙烯選擇性的提高，乙醇進料宜采用多股陸續(xù)加入加料方式。

[1]Rho S B.The conversion of ethanol to ethylene using a gas-phase catalytic reaction, including the effects of cyclic operation [J].International Chemical Engineering, 1984, 24 (3):567-577.

[2]Balandin A A, Balandin G V, Popov E I, et al.Use of radiocarbon for the investigation of the mechanism of the dehydration of ethanol over alumina [J].Russian Chemical Bulletin, 1958, 7 (2):218-220.

[3]Arai H, Take J I, Saito Y, et al.Ethanol dehydration on alumina catalysts:I.the thermal adsorption of surface compounds [J].Journal of Catalysis, 1967, 9 (2):146-153.

[4]Kn?zinger H, Bühl H, Ress E.The dehydration of alcohols over alumina virtue dependence of reaction direction on the substrate structure [J].Journal of Catalysis, 1968, 12 (2):121-128.

[5]Padmanabhan V R, Eastburn F J.Mechanism of ether from alcohols over alumina catalyst [J].Journal of Catalysis, 1972, 24 (1):88-91.

[6]Morávek V, Kraus M.Transient behavior of the system ethanol-diethyl ether-water-alumina [J].Journal of Catalysis, 1984, 87(2):452-460.

[7]Saito Y, Niiyama H.Reaction mechanism of ethanol dehydration on/in heterotopy compounds:analysis of transient behavior based on pseudo-liquid catalysis model [J].Journal of Catalysis, 1987, 106(2):329-336.

[8]Shinohara Y, Satozono H, Nakajima T, et al.Study of the interaction of ethanol with the broensted and Lewis acid sites on metal oxide surfaces using the DV-Xa method [J].Journal of Chemical Software, 1997, 4(2):41.

[9]Shinohara Yuji, Nakajima T, Suzuki S, et al.A computational chemical investigation of the dehydration and dehydrogenation of ethanol on oxide catalysts [J].Journal of Chemical Software, 1997, 4 (3):89.

[10]Shi B C, Dabbagh H A, Davis B H.Catalytic dehydration of alcohols, kinetic isotope effect for the dehydration oft-butanol [J].Topics in Catalysis, 2002, 18(3-4):59-265.

[11]Golay S, Doepper R, Renken A.Reactor performance enhancement under periodic operation for the ethanol dehydration overγ-alumina,a reaction with a stop-effect [J].Chemical Engineering Science, 1999, 54(20):4469-4474.

[12]索紅波, 蘇國東, 胡 燚, 等.工業(yè)氧化鋁催化劑上乙醇脫水動力學(xué) [J].化學(xué)反應(yīng)工程與工藝，2009, 25(5):458-462.Suo Hongbo, Su Guodong, Hu Yi, et al.Kinetics of ethanol dehydration over industrial catalystγ-Al2O3[J].Chemical Reaction Engineering and Technology, 2009, 25(5):458-462.

[13]余英哲.乙醇脫水制乙烯γ-Al2O3催化劑的分子模擬和實驗研究 [D].天津, 天津大學(xué), 2012.

[14]El-Salaam K M A, Hassan E A.Active surface centers in a heterogeneous CdO catalyst for ethanol decomposition [J].Surface Technology, 1982, 16(2):121-128.

[15]王 菊, 鐘思青, 謝在庫, 等.乙醇脫水制烯烴過程分析 [J].天然氣化工（C1 化學(xué)與化工), 2008, 33(5):27-32.Wang Ju, Zhong Siqing, Xie Zaiku, et al.Analysis for the process of ethanol dehydration to olefins [J].Natural Gas Chemical Industry,2008, 33(5):27-32.

[16]Kochar N K, Merims R, Padia A S.Ethylene from ethanol [J].Chemical Engineering Progress, 1981, 77(6):66-70.