乙醇为抑制剂时磺酸树脂催化异丁烯叠合反应的理论研究

石油炼制与化工PETROLEUM PROCESSING AND PETROCHEMICALS
2021年2月第52卷第2期
基础研究
乙醇为抑制剂时磺酸树脂催化异丁烯
叠合反应的理论研究
刘静】，肖之敏】，洪学思2,周晓龙1
（1.华东理工大学化工学院,上海200237； 2.中国石化上海石油化工研究院）
摘要：基于密度泛函理论，采用B3LYP/6-31G（d,p）计算方法对乙醇为抑制剂、磺酸树脂为催化剂的异丁烯叠合反应体系进行理论研究。

乙醇与异丁烯的醚化反应及异丁烯叠合反应的活化能分别为62.8kj/mol和167.8kj/mol,与异丁烯叠合反应相比，醚化反应的反应能垒较低，反应会优先发生。

乙醇在苯磺酸上的物理吸附能为一71.3kj/mol,异丁烯在苯磺酸上的物理吸附能为一44.8kj/mol,物理吸附态的异丁烯容易被磺酸基上的H质子化形成烷氧基中间体,烷氧态化学吸附的吸附能为一80.2kj/mol。

通过对比乙醇和异丁烯在苯磺酸上的吸附及两个反应的活化能，分析了乙醇在异丁烯叠合反应体系中的作用机理。

关键词：乙醇磺酸树脂异丁烯叠合活化能吸附能
我国具有丰富的G资源，催化裂化、蒸汽裂解、煤制烯烃等过程都会产生大量的Q烃，主要包括-丁烯、2-丁烯、异丁烯、异丁烷、正丁烷等［1］。

异丁烷和正丁烯可以通过直接烷基化反应来生产高辛烷值汽油添加组分一一三甲基戊烷；异丁烯可以与甲醇通过醚化反应生产甲基叔丁基醚（MTBE）⑷。

MTBE也曾作为高辛烷值添加组分加入到汽油池中，但是后来发现MTBE会污染地下水，欧洲部分国家已经禁止使用MTBE阪，导致大批MTBE生产装置闲置。

因此，如何合理利用生产MTBE的原料异丁烯是目前石化行业面临的一个重要课题。

探究C4烃类的转化规律，提高Q烃的高效利用具有广阔的应用前景。

异丁烯叠合反应的主要产物是二异丁烯（DIB）DIB经加氢后可得到异辛烷，其研究法辛烷值为100,且具有良好的抗爆性，是一种优良的高辛烷值清洁汽油调合组分［7］。

而且，将MTBE 的生产装置稍加改造即可用于异丁烯叠合反应,如此既可提高异丁烯的利用价值，又可降低DIB 的生产成本。

异丁烯叠合反应属于酸催化反应，遵守碳正离子机理。

在反应过程中，一分子异丁烯首先与催化剂的酸中心作用，形成叔碳正离子；叔碳正离子再与另一分子异丁烯叠合，形成C8碳正离子；该碳正离子既能脱去氢质子生成DIB,也能继续与异丁烯分子反应生成C12或C16碳正离子，进而生成三聚体或四聚体。

基于该反应机理,在异丁烯叠合反应过程中不可避免地会产生三聚体或四聚体。

三聚体加氢产物的沸点为170〜180c,处于汽油馏程（终馏点低于205c）的末端，其加入量需要严格控制，否则会影响汽油的蒸发性能；而四聚体加氢产物的沸点为230〜250c，已超出了汽油的馏程范围，不能添加到汽油中［10］。

因此,如果要得到高品质的汽油添加组分，必须抑制三聚体或四聚体的生成，提高二聚体的选择性。

有研究表明，在异丁烯叠合反应中添加醇类（甲醇、乙醇或叔丁醇）可以显著提高二聚体的选择性［1113］。

本课题组也曾研究过以磺酸树脂为催化剂时，乙醇对混合C4叠合反应体系中异丁烯和1-丁烯转化率以及产物分布的影响，结果表明，添加乙醇可以显著提高二聚体的选择性以及二聚体中三甲基戊烯的含量，有利于提高产品质量［14］,但未曾看到关于该反应体系的理论计算研究。

本研究基于密度泛函理论，采用B3LYP/ 6-31G（d,p）的计算方法对乙醇为抑制剂、磺酸树脂为催化剂的异丁烯叠合反应体系进行理论研究。

基于反应机理，搜索反应过程中所涉及到的关键过渡态，推测其反应路径，并计算乙醇与异丁烯的醚化反应及异丁烯叠合反应的活化能，通过活化能的大小比较不同反应的优先发生顺序。

此外，还对比了乙醇和异丁烯的吸附过程，通过吸附
收稿日期:2020-09-07；修改稿收到日期：2020-10-30.
作者简介：刘静，博士研究生，主要从事碳四资源的综合利用研究工作。

通讯联系人：周晓龙，E-mail：xiaolong@ 。

34石油炼制与化工2021年第52卷
能及反应活化能的对比，分析乙醇在异丁烯叠合反应体系中的作用机理。

1计算模拟方法
采用量子化学B3LYP/6-31G(d,p)方法对反应物、过渡态和产物的几何构型进行优化，在结构的优化过程中，通过查看是否有且只有一个虚频来判断过渡态的正确性，进一步通过内禀反应坐标(IRC)计算推测反应路径，并计算醚化反应和叠合反应的活化能。

此外，根据式(1)计算出乙醇在磺酸基上的吸附能。

△E ads=E EtOH/Ph-SO3H—E P h-SO3H—E ei OH(1)式中:AE ads为乙醇在磺酸基上的吸附能，kj/mol E eioh为优化后乙醇的能量,kj/mol；E ph_so3h为优化后苯磺酸的能量，kj/mol；E EtOH/Ph-so3h为优化后吸附态的能量，kj/mol。

本研究所有的建模与计算均是在Gaussian09软件上完成。

2结果与讨论
21磺酸树脂模型的建立
在异丁烯叠合反应的实验研究中，本课题组曾用到Amberlyst-15和DH-2磺酸树脂［1-5。

它们均属于大孔磺酸树脂，是由磺化的苯乙烯和二乙烯基苯交联而成的一种聚合物，其中磺酸基是树脂的活性中心［16］。

因此，在计算吸附能和反应活化能之前，需要先建立合理的磺酸树脂模型。

如果构建相应的高聚物模型，计算结果会更加准确，但是由于该模型的原子数较多，会使得计算量非常大，耗费时间较长。

Ma等［17］用密度泛函理论研究Amberlyst-15树脂催化苯酚和烯烃烷基化反应的机理时，用苯磺酸代表磺酸树脂来建模进行计算，所得到的计算结果与实验结果基本匹配，这说明在进行磺酸树脂的相关理论计算时，可以选择苯磺酸作为模型。

这样不仅可以真实地反映不同物种与活性中心之间的相互作用，还可以缩短计算时间。

因此，在本研究中也选用苯磺酸作为磺酸树脂的模型，优化后的苯磺酸模型如图1所示。

图1优化后的苯磺酸模型
2.2醚化反应活化能的计算
乙醇与异丁烯的醚化反应以乙醇分子和异丁烯分子共同吸附在磺酸基上为起点，吸附后的反应物之间进行加成反应得到吸附态的产物，最后产物从催化剂活性位点上脱附得到产物，同时催化剂的活性中心被还原。

通过计算可知，乙醇与异丁烯生成乙基叔丁基醚(ETBE)的反应是协同反应。

图2为醚化反应过程的吸附态反应物、过渡态及产物的结构示意；图3为醚化反应过程的能量变化,其中ini,ts-ini,ts,pro分别表示起始反应物、中间体、过渡态及产物。

表1为反应物、吸附态、过渡态及产物的几何构型参数(键长)。

(a)吸附态的反应物(b)过渡态
图2乙醇与异丁烯醚化反应过程涉及的结构示意(c)
产物
第2期刘静，等.乙醇为抑制剂时磺酸树脂催化异丁烯叠合反应的理论研究35
芋
PhSO 3H Js
V-°'f /-4.8T.
'、吸附/
\
'
Z
1
(匸0日・
3§
-67.6
,/脱去催化剂
\ts-prg^ ^36.4
图3醚化反应过程的能量变化
结合图2、图3及表1对醚化反应过程进行
分析。

首先乙醇和异丁烯共同吸附在苯磺酸上, 体系的能量降低67. 6 kj/mol,形成稳定的吸附 态，然后克服62. 8 kj/mol 的能垒形成过渡态。

在
反应过程中，首先，磺酸基上的H37进攻异丁烯双 键上富电子的C25,磺酸基上的O14—H37键长由
0.099 961 nm 不断拉长至0. 146 426 nm ,导致 O14—H37键断裂，释放出H37；异丁烯的C25=C28
双键键长由0. 134 526 nm 伸长至0. 141 016 nm ,
形成叔丁基碳正离子。

碳正离子很容易进攻共吸
附的乙醇分子上电负性较大的O24,使得O24与碳 正离子C28之间的距离不断缩短，进而形成共价键
O24—C28,键长为0. 250 058 nm ；与此同时，乙醇
分子上的 O24—H16 键长由 0096996nm 不断拉长
至0. 098 637 nm,导致键的断裂，释放出H16,该质
子与磺酸基上的O15结合，形成O15—H16新键, 树脂的活性中心被还原，同时形成了吸附在活性
中心上的醚化产物(ETBE )由过渡态到吸附态 产物，释放了 131. 6 kj/mol 的能量，这也说明乙醇 与异丁烯的醚化反应是一个放热过程。

最后，吸
附态的产物需要再从体系中获得80.9 kj/mol 的 能量后从活性中心上脱附。

表1醚化反应过程中反应物、吸附态、过渡态及产物的几何构型参数(键长)
nm
项目ini
ts-ini
ts ts-pro
pro
O24—H16
0.096 9960 0972230098637C25—C280.133 952
0 1345260141016O14—H37
0 099961
0146426
C25—H3701224840 1092390 109344O24—C2802500580 1465950 143947
O15—H16
0175438
0 102136
上完成。

图4和图5分别为烷氧基中间体的形2.3异丁烯叠合反应活化能的计算在进行异丁烯叠合反应的相关理论计算
时，遵循了 Eley-Rideal(E-R )反应机理。

该机理 认为异丁烯叠合反应的发生首先是异丁烯化学吸附在磺酸基上形成烷氧基中间体，然后再与
异丁烯分子反应形成DIB,反应在一个活性中心
成过程及异丁烯叠合反应过程中所涉及到的吸
附态反应物、过渡态和产物的结构示意，图6为 异丁烯叠合反应过程中的能量变化，表2为异 丁烯叠合反应过程中反应物、吸附态反应物、过
渡态及产物的几何构型参数。

(a)吸附态的异丁烯
⑹过渡态
图4烷氧基中间体形成过程涉及的结构示意
(c)
烷氧基
36石油炼制与化工2021年第52卷
@）反应物（b）过渡态（c）产物
图5异丁烯叠合反应过程涉及的结构示意
异丁烯分子首先被物理吸附在苯磺酸上，体系的能量降低44.8kj/mol,吸附态的异丁烯容易被磺酸基上的羟基H质子化，进而形成烷氧基中间体。

该过程需要克服74.9kj/mol的能垒。

根据优化后的分子模型可知，磺酸基上的28H到异丁烯两个双键碳（C16和C19）的距离分别为0.208711nm和0.231927nm,因而H28更容易进攻取代基较少的C16。

在过渡态的形成过程中，磺酸基上的H28进攻异丁烯双键上C16形成叔丁基碳正离子，该碳正离子的C19很容易进攻磺酸基上电负性较大的015。

C19和015的距离由0.237622nm缩短至0.148845nm,形成015—C19烷氧键；同时C16—C19的键长逐渐拉长至0.152642nm,碳碳双键变成单键，最终形成烷氧基中间体。

表2异丁烯叠合反应过程中反应物、吸附态反应物、过渡态及产物的几何构型参数（键长）nm
项目n tsl-ini ts1tsl-pro ts2-ini ts2ts2-pro pro
C16—C190.1335700.1343250.1410210.152642
014—H280.141741
H28—C160.2087110.1241010.109078
015—C190.2376220.1488450.1494410.377868
C29—C320.1338030.1388420.1514440.151569 C32—C370.1505450.1481870.1344770.133717 C19—C290.2058870.1558090.155839 H39—O130.2340100.098948
在进行异丁烯叠合反应时，游离的异丁烯分子先物理吸附在所形成的烷氧基中间体（s-pro）附近，体系的能量降低37.9kj/mol,得到较稳定的吸附态。

烷氧基中间体与吸附的异丁烯分子发生反应，需要克服167.8kj/mol的能垒得到ts2过渡态，然后释放166.4kj/mol的能量得到吸附态的产物（ts2-pro）,这也说明异丁烯叠合反应是放热反应。

根据过渡态的虚频振动方向，推测在s2过
第2期刘静，等.乙醇为抑制剂时磺酸树脂催化异丁烯叠合反应的理论研究37
渡态的形成过程中，烷氧中间体上的O15—C19键长由0.149441nm不断拉长至0.377868nm,烷氧键断裂形成碳正离子，碳正离子C19会进攻所吸附的另一分子异丁烯双键上的C29,使得异丁烯分子的C29—C32键长拉长至0.138842nm；同时C32—C37单键的键长缩短变成双键，H39与磺酸基的013成键，磺酸基活性中心被还原，并得到吸附在磺酸基上的叠合产物（DIB）最后，从体系中吸收40.0kj/mol的能量后，产物从活性中心上脱附。

对比醚化反应和异丁烯叠合反应，醚化反应的活化能为62.8kj/mol,异丁烯叠合反应的活化能为167.8kj/mol,前者的活化能比后者低105.0kj/mol,说明醚化反应比叠合反应更容易发生。

此外，叠合产物DIB从活性中心上脱附需要40.0kj/mol 的能量，醚化产物ETBE的脱附需要80.9kj/mol 的能量。

相对于DIB,ETBE需要获得更多的能量才能从活性中心上脱附，这可能是因为ETBE 中的O与磺酸基上的羟基H存在氢键相互作用，因此生成的ETBE可能并不能及时从活性中心上脱附使活性中心得到还原，而是仍占据在活性中心上，使有效的活性中心数量减少。

这可能也是引入乙醇后二聚选择性有所提高的一个原因。

2.4乙醇的作用机理分析
乙醇与异丁烯的醚化反应是可逆反应，当该反应达到平衡时，体系中仍会有未转化的乙醇存在。

体系中游离的乙醇与树脂的磺酸基之间由于存在氢键相互作用会吸附在磺酸基上，而这种吸附只是一种物理吸附。

根据吸附能的计算式，算出乙醇在磺酸基上的吸附能为一71.3kJ/mol。

由2.3节的计算结果可知，异丁烯物理吸附在苯磺酸上的吸附能为一44.8kj/mol。

物理吸附态的异丁烯由于自身的电负性强容易被磺酸基上的H质子化形成烷氧基中间体，烷氧态化学吸附的吸附能为一80.2kJ/mol。

根据吸附能的定义可知，吸附能为负值时，其绝对值越大，说明两者的吸附作用越强，所形成的吸附态越稳定，该物质越容易吸附。

虽然形成的烷氧基中间体的能量更低，但由物理吸附态转变为更稳定的化学吸附态的过程需要一定的反应时间，还要克服一定的反应能垒，因此，当乙醇存在时，优先吸附在磺酸基上的是乙醇，剩余的活性位被化学吸附态的异丁烯分子占据，乙醇的吸附能力比异丁烯强。

结合乙醇和异丁烯在磺酸基上的吸附能力以及醚化反应和叠合反应的活化能，对乙醇的作用机理进行详细的分析。

在磺酸树脂催化异丁烯叠合反应体系中，当添加乙醇时，由于乙醇的吸附能力较强，故其优先吸附在活性中心上。

当异丁烯分子也吸附在磺酸基附近时，由于醚化反应的活化能较低，异丁烯分子会优先与乙醇分子发生反应。

醚化反应是可逆反应，当该反应达到平衡时,异丁烯分子开始发生聚合反应。

由于此时树脂中的一部分活性中心及反应物中的一部分异丁烯已经参与了醚化反应，因此对于异丁烯聚合反应，有效的活性中心数量减少，反应物的实际浓度降低,使聚合反应速率降低，多聚反应受到抑制，从而提高了二聚体选择性。

当醚化反应达到平衡后，体系中仍有游离的乙醇分子，由于乙醇分子与磺酸基之间存在氢键相互作用，使磺酸基为叠合反应提供氢质子的能力减弱，也就是说树脂的酸强度降低，活性降低。

此外，通过计算还发现，醚化产物ETBE从活性中心上脱附所需要的能量高于叠合产物DIB的脱附，说明生成的ETBE可能并不能及时从活性中心上脱附使活性中心得到还原，而是仍占据着活性中心，这将更进一步减少有效活性中心的数量。

总的来说，在异丁烯叠合反应体系中，乙醇是通过减少有效活性中心的数量、减弱树脂的酸强度以及降低反应物异丁烯的实际浓度来提高二聚体的选择性的。

3结论
基于密度泛函理论，利用苯磺酸来模拟磺酸树脂，推测醚化反应和异丁烯叠合反应的反应路径。

醚化反应以乙醇分子和异丁烯分子共同吸附在苯磺酸上为起点，叠合反应以异丁烯形成烷氧基中间体为起点，醚化反应和叠合反应的活化能分别为62.8kj/mol和167.8kj/mol0在磺酸树脂催化异丁烯叠合反应体系中，当乙醇存在时，磺酸树脂的活性中心优先吸附乙醇，剩余的活性中心被化学吸附态的异丁烯所占据，且活化能较低的醚化反应优先发生。

因此，对于异丁烯叠合反应，反应物浓度的降低、有效活性中心数量的减少以及树脂酸强度的减弱均会使聚合反应速率降低，多聚反应发生的可能性降低，从而提高二聚体的选择性。

38石油炼制与化工2021年第52卷
参考文献
[1]何英华，詹海容,朱丽娜，等•混合G馏分综合利用技术进展
[].化工技术与开发，2013():4044
[2]Feller A,Guzman A,Zuazo I,et al.On the mechanism of
catalyzed isobutane/butene alkylation by zeolites]〕],journal ofCatalyss，2004，224(1)：80-93
[3]Zhang Honghua,Xu jiakuo,Tang Hongguo,et al.Isobutane/
2-butene alkylation reaction catalyzed by Cu modified RE-X zeolite[J].Industrial&Engineering Chemistry Research, 2019，58(22)：9690-9700
[]王进善，栾祥山，杨宗仁.甲基叔丁基醚(MTBE)生产技术的开发[]•齐鲁石油化工,991,9()：175-186
[5]Tsai MJ,Kolodziej R,Ching D,et al.Consider new technologies
toreplace MTBE[j]．Hydrocarbon Processing，2002，81
(2)：81-87
[]毕建国.烷基化油生产技术的进展[]•化工进展,2007,26
(7)：934-939
[7]温郎友，吴巍，刘晓欣.间接烷基化技术进展[]•当代石油
石化，2004,12():36-40
[]何奕工，舒兴田，龙军•正碳离子和相关的反应机理[]•石油学报(石油加工)2007,23():17
[9]AlcantaraR，AlcantaraE，CanoiraL，etal Trimerizationof
isobutene over Amberlyst-15catalyst[j].Reactive& FunctionalPolymers，2000，45(1)：19-27[10]Girolamo M D，Tagliabue L MTBEandalkylateco-production：
Fundamentals and operating experience[j]Catalysis
Today，1999，52(2)：307-319
[11]Honkela M L,Krause AOI.Influence of polar components
in the dimerization of isobutene[J].Catalysis Letters,2003,
87(3)：113-119
[12]Girolamo M D，Marchionna M Acidic and basic ion
exchangeresinsforindustrialapplications[j]journalof
MolecularCatalysisA：Chemical，2001，177：33-40
[13]OuniT，Honkela M，Kolah A，etal Isobutenedimerization
ina miniplant-scalereactor[j]ChemicalEngineeringand
Processing，2006，45：329-339
[14]Liujing，Ding Ning，Ge Yuena，etal Dimerization of
isobutene in C4mixtures in the presence of ethanol over acid
ion-exchangeresin DH-2[j]．Catalysis Le t ers，2019，149：
1277-1285
[15]Liu Jing,Ge Yuena,Song Yueqin,et al.Dimerization of iso­
butene on sodium exchanged Amberlyst-15resins[j]．CatalysisCommunications，2019，119：57-61
[16]Sharma M M．Somenovelaspectsofcationicion-exchange
resins as catalysts[J]•Reactive&Functional Polymers,
1995，26：3-23
[17]MaQ，Chakraborty D，Faglioni F，etal Alkylation of
phenol：A mechanistic view[j]journal of Physical
ChemistryA，2006，110(6)：2246-2252
THEORETICAL STUDY ON ISOBUTENE POLYMERIZATION OVER
SULFONIC ACID RESIN IN PRESENCE OF ETHANOL
Liu Jing1,Xiao Zhimin1,Hong Xuesi2,Zhou Xiaolong1
(1.East China University of Science and Technology,School of Chemical Engineering,Shanghai200237；
2.SINOPEC Shanghai Re s earch In s thute of P e trochemical Technology')
Abstract：Thepolymerizationofisobuteneoversulfonicacidresininthepresenceofethanolwas studied using DFT theory and B3LYP/6-31G(d,p)algorithm method.The activation energies of the etherificationreaction of ethanol and isobutylene and the polymerization reaction ofisobutyleneare 62.8kj/mol and167.8kj/mol,pared with the polymerization reaction of isobutylene, thereaction energy barrier ofthe etherification reaction is lower and the reaction wi l take place preferentially.The physical adsorption energy of ethanol on benzenesulfonic acid is—71.3kj/mol,and the physical adsorption energy of isobutylene on benzenesulfonic acid is—44.8kj/mol.The physically adsorbed isobutylene is easily protonated by the H of the sulfonic acid group to form alkoxy intermediate with chemical adsorption energy of—80.2kj/mol.Finally,the mechanism of ethanol in the isobutene polymerzatonsystem wasanalyzed.
Key Words：ethanol；sulfoncacdresn；sobutenepolymerzaton；actvatonenergy；adsorptonenergy。