[发明专利]一种煤焦油加氢氢耗的动力学计算方法有效
申请号: | 201310703104.4 | 申请日: | 2013-12-12 |
公开(公告)号: | CN103745088A | 公开(公告)日: | 2014-04-23 |
发明(设计)人: | 李冬;李稳宏;孙智慧;田盼盼;刘鑫;李斌;李学坤;范峥 | 申请(专利权)人: | 西北大学 |
主分类号: | G06F19/00 | 分类号: | G06F19/00;G06F17/50 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 710069 *** | 国省代码: | 陕西;61 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明涉及一种煤焦油加氢氢耗的动力学计算方法,属于煤焦油技术领域。所用动力学模型计算方法具体步骤为:(1)动力学模型基本假设;(2)进行轴向扩散分析;(3)动力学模型的建立;(4)求解动力学模型;(5)模型验证:为了验证模型的可靠性和预测能力,进行重复实验予以验证。本发明通过实验测定值与模型预测值的对比分析发现,煤焦油加氢化学氢耗计算模型的相对误差为2.10%,具有较好的预测性,在一定程度上能较好的反应煤焦油加氢过程氢气消耗的实际情况,对煤焦油加氢工艺设计具有一定的指导意义。 | ||
搜索关键词: | 一种 煤焦油 加氢 动力学 计算方法 | ||
【主权项】:
1.一种煤焦油加氢氢耗的动力学计算方法,其特征在于:所用动力学模型计算方法具体步骤为:(1)动力学模型基本假设:本模型的建立基于以下假设:①催化剂床层润湿完全,固定相与床层之间不发生相互位移;②气、液两相发生拟均相反应;③反应装置在恒温、恒压条件下运行;④试验在稳态条件下操作,催化剂活性不随时间改变;⑤反应速率常数与温度的关系符合阿伦尼乌斯公式;(2)进行轴向扩散分析:为了确保反应器在平推流状态下操作,必须考察能否忽略返混;用于估计可忽略返混影响的最小床层长度Lb计算规则:L b d p > 20 n p e z · ln 1 1 - x - - - ( 1 ) ]]> 式(1)中Lb为催化剂床层高度;dp为颗粒直径;n为反应级数;x为馏分转化率;Pez为Peclet值,根据雷诺数函数估算;(3)动力学模型的建立:(3a)加氢精制(加氢脱硫、加氢脱氮、加氢脱氧、烯烃饱和)氢耗模型的建立:依据等温平推流模型,在管式反应器中,加氢精制反应氢耗微分质量平衡方程为:(1-ε)rAdVR=-QLdcH (2)式(2)中ε为催化剂床层空隙率,rA为氢耗表观反应速率,VR为催化剂床层体积,QL为原料油体积流量,其值为催化剂床层总体积与空速的乘积,cH为当量氢浓度,也就是在反应条件下可能与氢气发生反应的化学键的浓度,单位μg-H2/g-oil;各类化学氢耗的表观反应速率表示为:rA=kv(cH)n (3)式(3)中n为反应级数,kv为反应速率常数;kv值与反应压力项PH2α和温度T有关,根据阿伦尼乌斯公式kv表示为:k v = k 0 exp ( - Eα / RT ) P H 2 α - - - ( 4 ) ]]> 加氢精制氢耗的速率表达式可写为:r A = k 0 exp ( - Eα / RT ) ( c H ) n P H 2 α - - - ( 5 ) ]]> 式(5)中k0为Arrhenius方程的指前因子,Eα为反应的表观活化能,J·mol-1,T为反应温度;若反应为非一级反应,即n≠1,将式(5)中的rA带入式(3)中并整理得:( c H ) 1 - n - ( c H 0 ) 1 - n = ( n - 1 ) ( 1 - ϵ ) ( V R / Q L ) k 0 exp ( - Eα / RT ) P H 2 α - - - ( 6 ) ]]> 否则为一级反应,n=1,有:ln c H 0 c H = ( 1 - ϵ ) ( V R / Q L ) k 0 exp ( - Eα / RT ) P H 2 α - - - ( 7 ) ]]> 求解式(6)和式(7)即可得到氢耗计算公式H chem = c H 0 - c H = c H 0 - ( ( n - 1 ) ( 1 - ϵ ) ( V R Q L ) k 0 exp ( - Eα RT ) P H 2 α - ( c H 0 ) 1 - n ) ( 1 / 1 - n ) ( n ≠ 1 ) c H 0 ( 1 - exp ( - ( 1 - ϵ ) ( V R Q L ) k 0 exp ( - Eα RT ) P H 2 α ) ( n = 1 ) - - - ( 8 ) ]]> (3b)加氢裂化和芳烃饱和氢耗模型的建立:采用二次模型对芳烃饱和氢耗和加氢裂化氢耗进行回归拟合,模型如下:H = α 0 + α 1 T HDT + α 2 T HDC + α 3 LHSV + α 4 P H 2 + α 11 T HDT 2 + α 22 T HDC 2 + α 33 LHSV 2 + α 44 P H 2 2 + α 12 T HDT T HDC + α 13 T HDT LHSV + α 14 T HDT P H 2 + α 23 T HDC LHSV + α 24 T HDC P H 2 + α 34 P H 2 LHSV - - - ( 9 ) ]]> (4)求解动力学模型:(4a)加氢精制氢耗模型的求解(4a.1)反应级数的确定:反应级数是由化学反应机制决定的,其值的准确性在一定程度上反映了动力学方程表达的客观性水平,对于化学氢耗动力学参数,若直接通过方程(8)求取比较困难,根据动力学原理,速率方程还可表达为:r A = - dc H dt = kc H 0 - - - ( 10 ) ]]> 对上式进行线性化处理有:ln r A = ln ( - dc H dt ) = ln k + n ln c H - - - ( 11 ) ]]> 在催化剂床层温度、反应压力和氢油体积比一定的情况下,反应速率仅与反应时间有关,又t=1/LHSV,因此对式(11)进行线性回归拟合即可得到反应级数n;(4a.2)其他动力学参数的求取:为了便于回归拟合以确定各化学氢耗的动力学参数,对式(6)(7)进行线性化处理有,若反应级数不是1级则:ln [ ( c H ) 1 - n - ( c H 0 ) 1 - n ] - ln [ ( n - 1 ) ( 1 - ϵ ) ( V R Q L ) ] = ln k 0 + α ln ( P H 2 ) - E α RT - - - ( 12 ) ]]> 若反应级数n=1,则线性化处理结果为:ln ( ln c H 0 c H ) - ln [ ( 1 - ϵ ) ( V R / Q L ) ] ln k 0 - Eα RT + α ln P H 2 - - - ( 13 ) ]]> (4b)加氢裂化和芳烃饱和氢耗模型的求解:将实验数据通过SPSS(Statistical Product and Service Solutions)软件对式(9)进行非线性拟合,拟合采用莱温伯格-麦夸特(Levenberg-Marquardt)法;(5)模型验证:为了验证模型的可靠性和预测能力,进行重复实验予以验证。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西北大学,未经西北大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201310703104.4/,转载请声明来源钻瓜专利网。
- 同类专利
- 专利分类
G06 计算;推算;计数
G06F 电数字数据处理
G06F19-00 专门适用于特定应用的数字计算或数据处理的设备或方法
G06F19-10 .生物信息学,即计算分子生物学中的遗传或蛋白质相关的数据处理方法或系统
G06F19-12 ..用于系统生物学的建模或仿真,例如:概率模型或动态模型,遗传基因管理网络,蛋白质交互作用网络或新陈代谢作用网络
G06F19-14 ..用于发展或进化的,例如:进化的保存区域决定或进化树结构
G06F19-16 ..用于分子结构的,例如:结构排序,结构或功能关系,蛋白质折叠,结构域拓扑,用结构数据的药靶,涉及二维或三维结构的
G06F19-18 ..用于功能性基因组学或蛋白质组学的,例如:基因型–表型关联,不均衡连接,种群遗传学,结合位置鉴定,变异发生,基因型或染色体组的注释,蛋白质相互作用或蛋白质核酸的相互作用
G06F 电数字数据处理
G06F19-00 专门适用于特定应用的数字计算或数据处理的设备或方法
G06F19-10 .生物信息学,即计算分子生物学中的遗传或蛋白质相关的数据处理方法或系统
G06F19-12 ..用于系统生物学的建模或仿真,例如:概率模型或动态模型,遗传基因管理网络,蛋白质交互作用网络或新陈代谢作用网络
G06F19-14 ..用于发展或进化的,例如:进化的保存区域决定或进化树结构
G06F19-16 ..用于分子结构的,例如:结构排序,结构或功能关系,蛋白质折叠,结构域拓扑,用结构数据的药靶,涉及二维或三维结构的
G06F19-18 ..用于功能性基因组学或蛋白质组学的,例如:基因型–表型关联,不均衡连接,种群遗传学,结合位置鉴定,变异发生,基因型或染色体组的注释,蛋白质相互作用或蛋白质核酸的相互作用