石脑油催化裂解结构化反应器反应特性的CFD模拟
许敏, 刘辉*, 王红霞
北京化工大学,化工资源有效利用国家重点实验室,能源环境催化北京市重点实验室, 北京100029
通讯联系人:刘 辉。

作者简介:许 敏,1992年生,女,山东省菏泽市人,在读硕士研究生。

摘要

以石脑油为原料,采用催化裂解六集总动力学模型,建立描述结构化反应器内催化裂解的反应器数学模型,并利用CFD软件对结构化反应器内的石脑油催化裂解性能进行数值模拟。通过改变孔道直径、反应器长度以及反应器内温度、气体入口速率考察反应器结构尺寸和反应条件对目标产物乙烯、丙烯的收率及石脑油转化率的影响。结果表明,反应器孔道直径的增加,目标产物收率减小,反应器长度20 mm时反应完全,升高反应温度和增大入口速率均有利于目标产物的生成。在入口温度680 ℃和入口速率0.4 m·s-1条件下,石脑油转化率92%,乙烯收率19.3%,丙烯收率23.1%。而在相同反应条件下的固定床反应器中乙烯收率10.3%,丙烯收率13.3%,石脑油转化率80.0%。

关键词: 石油化学工程; 石脑油; 集总动力学模型; 催化裂解; 结构化反应器; 数值模拟
中图分类号:TE624.4;O643.12    文献标志码:A    文章编号:1008-1143(2018)05-0127-06
CFD simulation of naphtha catalytic cracking reaction in monolith reactors
Xu Min, Liu Hui*, Wang Hongxia
State Key Laboratory of Chemical Resource Engineering,Beijing Key Laboratory of Energy Environmental Catalysis,Beijing University of Chemical Technology, Beijing 100029,China
Abstract

Mathematic model for naphtha catalytic cracking in monolith reactor is established using six lumped kinetics by CFD simulations. Effects of reactor size and reaction conditions on yield of ethylene, propylene and naphtha conversion are investigated by changing channel diameter,reactor length,reaction temperature and gas inlet velocity. The results show that increase in channel diameter leads to decrease of products yield.Reaction completes at channel length around of 20 mm.Increasing in reaction temperature and entrance velocity are beneficial to production of target products. The monolith reactor can achieve a propylene yield up to 23.1% and ethylene yield to 19.3% with a naphtha conversion being 92% at 680 ℃ and 0.4 m·s-1; this is in contrast to a typical fixed reactor under the same reaction conditions, in which case the ethylene yield is 10.3%, the propylene yield is 13.3%, and the conversion rate of naphtha was 80.0%.

Keyword: petrochemical engineering; naphtha; lump kinetics; catalytic cracking; monolith reactor; numerical simulation

丙烯和乙烯为石油化工基础原料, 主要来自蒸气裂解(SC)装置和催化裂化(FCC)装置[1]。近年来对乙烯和丙烯的需求不断增大[2], 如何通过改进反应器技术提高乙烯和丙烯收率显得尤为重要。传统催化裂解工艺中通常使用提升管反应器、固定床反应器和流化床反应器。其中, 以石脑油为原料, 在固定床反应器中催化裂解的烯烃收率较高, 但反应温度降低幅度不大, 难以从根本上克服蒸气裂解工艺能耗高投资大的局限[3, 4]; 研究机构开发的石脑油催化裂解制烯烃技术[5]均采用流化床反应工艺。

结构化反应器具有众多相互平行的、规则的直通孔道, 催化活性组分被制成极薄的涂层结构[(10~150) μ m]负载于孔道内壁; 由于其所具有的特殊蜂窝孔道结构, 故具有低压力降、床层分布均匀、比表面积大、无催化剂磨损、放大简单、操作灵活和易分离等优点, 能够克服常规颗粒型固定床反应器的问题[6]

本文建立石脑油催化裂解结构化反应器数学模型, 采用CFD软件对石脑油催化裂解反应进行不同反应条件下的模拟计算, 并基于石脑油转化率、乙烯和丙烯收率模拟结果, 与相同反应条件下的固定床反应器中催化裂解结果进行对比, 以期为此类新反应器的应用提供依据。

1 反应器模型
1.1 物理模型

选取常用的单孔道模型[7, 8], 分为气相(孔道内的气体流动)和催化剂涂层两个区域, 两者通过流-固表面相耦合。反应器长度L=30 mm, 内半径R1=0.475 mm, 涂层厚度0.025 mm。鉴于涂层厚度相对流动通道较薄, 因此忽略涂层厚度, 将涂层中的反应视为表面反应。通道为圆管, 可只对纵向x和径向r进行轴对称性二维定态模拟。单孔道反应器结构如图 1所示。

图1 单通道反应器结构Figure 1 Single channel configuration

1.2 控制方程[9]

对于反应通道内二维轴对称气相流动过程, 建立了包含连续性方程、动量方程、能量方程和组分质量方程的通用形式的控制方程。

aiuφix+1rraivφir=xΓiφix+1rrrΓiφir+Si(1)

式中, ϕ i为通用变量, u(i=2)为轴向速率, υ (i=3)为径向速率, T(i=4)为温度变量, CAi(i=5~15)为组分浓度, ai为参数, Si为广义源项, Γ i为广义扩散系数。通用方程各符号的形式如表1所示。

表1 通用方程各符号的形式 Table 1 The specific forms of the symbols in the general equation
1.3 边界条件[10]

(1) 反应通道入口

u=U0; v=0; CA=CA0; T=T0(2)

(2) 反应通道轴对称线

ur=vr=0; CAr=0; Tr=0(3)

(3) 反应通道出口

ur=vr=0; CAz=0; Tz=0(4)

(4) 壁面

u=0; v=0; CAr=0; Tw=T0(5)

1.4 集总反应动力学模型

催化裂解原料石脑油的性质如表2所示。

表2 催化裂解原料石脑油的性质[11] Table 2 Properties ofnaphtha

反应动力学采用六集总反应动力学模型[11], 结果如图2所示。

图2 石脑油催化裂解六集总动力学模型Figure 2 Six-lumped kinetic model for naphtha kinetic catalytic pyrolysis

图2中的反应均为一级反应。由于通道壁面催化剂涂层较薄, 且反应温度小于700 ℃可以忽略内部扩散的影响[12], 所以考虑壁面为表面反应。前述六集总动力学为体积反应速率(kmol· m-3· s-1), 因此要将其转换为表面反应速率(kmol· m-2· s-1)[13]。单通道中内表面积与涂层体积之比为:

paw=2πR1Lπ×L×R22-R12(6)

转换后的表面反应速率为:

-rSi=-rvipaw=Aiexp-EaiRTCA(7)

转换后的动力学模型参数如表3所示。

表3 六集总动力学模型的活化能和指前因子[11] Table 3 The pre-exponential factor and activation energy of the six lumps model

用有限体积差分法求解前述的耦合微分方程组, 使用CFD软件求解, 并对结果进行网格无关性检验。模拟所需基本条件如表4所示。

表4 模拟所需基本条件 Table 4 Basic simulation conditions
2 结果与讨论
2.1 反应物和产物分布

在不同入口速率和反应温度下进行预模拟计算, 确定入口速率0.4 m· s-1和反应温度680 ℃时反应为最优结果。反应物及产物沿反应器的轴向分布如图 3所示。

图3 反应物及产物沿反应器的轴向分布Figure 3 Distribution of Naphtha and the product along the reactor

由图3可见, 反应物及产物质量分数沿轴向呈指数变化, 乙烯和丙烯是主要产物。

随反应不断进行, 壁面丙烯浓度增加。

不同轴向位置丙烯径向分布如图4所示。由图4可见, 随距离增加, 产物丙烯的径向分布先增后趋向稳定不变。结合图3可知, 在L=20 mm时, 反应已经完全; 相应的停留时间约为0.1 s, 对停留时间分布模拟计算得出方差为0.035, 模型参数N=28.5, 表明通道内流型接近平推流。

图4 不同轴向位置丙烯径向分布Figure 4 The radial distribution of C3H6 at different x axis

2.2 入口温度

在入口速率0.4 m· s-1条件下, 通过改变反应温度模拟催化裂解过程, 结果如图 5所示。由图5可见, 对吸热的催化裂解反应, 随着反应温度升高, 石脑油转化率增加, 乙烯和丙烯收率也明显增加, 680 ℃时石脑油转化率92.0%, 乙烯收率19.3%, 丙烯收率23.1%。从图3的产物分布中也可以看出。

图5 反应温度对石脑油转化率及乙烯、丙烯收率的影响Figure 5 Effect of reaction temperature on naphtha conversion, yield of C2H4and C3H6

2.3 入口速率

固定反应器孔径和长度, 改变反应器物料入口速率, 模拟不同反应温度条件下的反应过程。丙烯模拟结果如图6所示。由图6可见, 丙烯收率随着入口速率的增加先增后减, 入口速率0.4 m· s-1时具有最佳的反应结果。

图6 入口速率对丙烯收率的影响Figure 6 Effect of inlet velocity on the yield of C3H6 products

2.4 孔道直径

选取入口速率0.4 m· s-1, 改变反应器孔道直径, 模拟催化裂解过程, 表明孔道直径的增加降低了石脑油转化率, 降低了乙烯和丙烯收率。孔道直径变化对丙烯收率的影响如图7所示。由图7可见, d=1 mm时, 丙烯收率最高。

图7 反应器孔道直径对丙烯收率的影响Figure 7 Effect of reactor diameter on the yield of C3H6

2.5 与固定床反应器的对比

在入口温度680 ℃和入口速率0.4 m· s-1条件下, 将结构化反应器中石脑油催化裂解反应结果与相同反应条件下的固定床反应器中石脑油催化裂解反应结果[11]进行对比, 如表4所示。由表4可见, 结构化反应器中乙烯收率比固定床反应器高9.0%, 丙烯收率高9.8%。

表4 结构化反应器与固定床反应器反应结果对比 Table 4 Comparison of the results of the reaction between a monolith reactor and a fixed bed reactor
3 结 论

(1) 建立了结构化反应器单通道模型, 基于石脑油催化裂解六集总动力学模型数值模拟表明, 增大反应器直径, 降低了反应物转化率和目标产物收率。升高反应温度, 增大入口速率有利于提高反应转化率和产物收率。在本文模拟条件下, 通道长度20 mm时反应接近完全;

(2) 与传统固定床反应器相比, 结构化反应器能得到更高的石脑油转化率、乙烯和丙烯收率。在入口温度680 ℃和入口速率0.4 m· s-1条件下, 结构化反应器中乙烯收率19.3%, 丙烯收率23.1%, 石脑油转化率92%。相同反应条件下, 固定床反应器中乙烯收率10.3%, 丙烯收率13.3%, 石脑油转化率80.0%。

符号说明:

ai:方程(1)中的参数;

L:反应器长度, mm;

R1:反应器内半径, mm;

R2:反应器外半径, mm;

u:气体轴向线速度, m· s-1;

v:气体径向线速度, m· s-1;

k:热导率, W· (m· K)-1;

x:轴坐标, mm;

CA:摩尔浓度, mol· m-3;

rSi:第i个表面反应的反应速率, kmol· m-2· s-1;

rvi:第i个体积反应的反应速率, kmol· m-3· s-1;

Ai:第i个表面反应指前因子, m· s-1;

Eai:第i个表面反应活化能, kJ· mol-1;

DAB:组分扩散系数;

Si:广义源项;

T:反应温度, K;

P:系统压力, Pa;

M:组分摩尔质量, g· mol-1;

H:反应的焓变, J· mol-1;

Cp:比热容, J· (kg· K)-1;

ρ :气体密度, kg· m-3;

μ :气体黏度, Pa· s;

λ :热导率, W· (m· K)-1;

Γi:方程(1)中的扩散系数;

i:变量。

The authors have declared that no competing interests exist.

参考文献
[1] 曹湘洪. 增产丙烯, 提高炼化企业盈利能力[J]. 化工进展, 2003, 22(9): 911-919.
Cao Xianghong. Enhance the economic gain ability of petroleum refinery industry by increasing propylene production[J]. Chemical Industry and Engineering Progress, 2003, 22(9): 911-919. [本文引用:1]
[2] 李振宇, 王红秋, 黄格省, . 我国乙烯生产工艺现状与发展趋势分析[J]. 化工进展, 2017, 36(3): 767-773.
Li Zhenyu, Wang Hongqiu, Huang Gexing, et al. Analysis of current situation and development trend of ethylene production technology in China[J]. Chemical Industry and Engineering Progress, 2017, 36(3): 767-773. [本文引用:1]
[3] 王建明. 催化裂解生产低碳烯烃技术和工业应用的进展[J]. 化工进展, 2011, 30(5): 911-917.
Wang Jianming. Progress in the production of low carbon olefin by catalytic cracking[J]. Chemical Industry and Engineering Progress, 2011, 30(5): 911-917. [本文引用:1]
[4] 魏晓丽, 毛安国, 张久顺, . 石脑油催化裂解反应特性及影响因素分析[J]. 石油炼制与化工, 2013, 44(7): 1-7.
Wei Xiaoli, Mao Anguo, Zhang Jiushun, et al. Study on reaction characteristics and influence factors of naphtha catalytic cracking[J]. Petroleum Processing and Petrochemicals, 2013, 44(7): 1-7. [本文引用:1]
[5] 刘剑, 孙淑坤, 张永军, . 石脑油催化裂解制低碳烯烃技术进展及其技术经济分析[J]. 化学进展, 2011, 29(11): 33-37.
Liu Jian, Sun Shukun, Zhang Yongjun, et al. Advance in catalytic cracking of naphtha to light olefins and techno-economic analysis[J]. Chemical Industry, 2011, 29(11): 33-37. [本文引用:1]
[6] 龙军, 邵潜, 贺振富, . 规整结构催化剂及反应器研究进展[J]. 化工进展, 2004, 23(9): 925-932.
Long Jun, Shao Qian, He Zhenfu, et al. Monolithic catalysts and monolithic reactors[J]. Chemical Industry and Engineering Progress, 2004, 23(9): 925-932. [本文引用:1]
[7] Wei J, Kuo J C W. A lumping analysis in monomolecular reaction systems[J]. Industrial & Engineering Chemistry Fundamentals, 1969, 8(1): 233-243. [本文引用:1]
[8] Guo W, Xiao W, Luo M. Comparison among monolithic and rand omly packed reactors for the methanol-to-propylene process[J]. Chemical Engineering Journal, 2012, 207-208: 734-745. [本文引用:1]
[9] Chen J, Yang H, Wang N, et al. Mathematical modeling of monolith catalysts and reactors for gas phase reactions[J]. Applied Catalysis A General, 2008, 345(1): 1-11. [本文引用:1]
[10] Liu H, Zhao J, Li C, et al. Conceptual design and CFD simulation of a novel metal-based monolith reactor with enhanced mass transfer[J]. Catalysis Today, 2005, 1505(3): 401-406. [本文引用:1]
[11] 吴小军. 石脑油催化裂解制乙烯丙烯反应动力学研究[D]. 上海: 华东理工大学, 2011.
Wu Xiaojun. kinetic model study on naphtha catalytic cracking for ethylene and propylene production[D]. Shanghai: East China University of Science and Technology, 2011. [本文引用:2]
[12] Hayes R E, Kolaczkowski S T. Mass and heat transfer effects in catalytic monolith reactors[J]. Chemical Engineering Science, 1994, 49(21): 3587-3599. [本文引用:1]
[13] Irani M, Alizadehdakhel A, Pour A N, et al. CFD modeling of hydrogen production using steam reforming of methane in monolith reactors: surface or volume-base reaction model?[J]. International Journal of Hydrogen Energy, 2011, 36(24): 15602-15610. [本文引用:1]