甲烷水合物常压分解

t= 2
t/ min
0 2 8 17 20 23 34 39 42 56 62 67 77 84
V r/ L
0 7. 419 18. 821 32. 111 36. 115 38. 470 47. 279 49. 741 50. 283 54. 444 56. 769 56. 968 58. 614 58. 927
∃ m(t)
t= 0
%
∃ V total =
Vr
t= 0
∃%
m( t) =
t= 0
P V tot al RT
式中, V r 为湿式流量计测量的气体体积( L ) ; V total为气体体积记录仪中收集到的总体积( L ) 。体积由实验测
定, 进而将 nH/ n0 表示出来( 实验数据见表 1) , 得出的( nH/ n0) t 曲线见图 2。
实验装置中设置缓冲罐, 容积为 12 L, 最高工作压力为 15 M Pa。水合物形成过程中的耗气量由流量计 测定, 其重复精度均为量程的 2% 。流量显示仪表可同时显示瞬间流量和累积流量, 采用数据采集系统进 行流量、温度、压力的采集。此外, 分解气体体积由湿式气体记录仪记录。 1. 2 过程
实验前用蒸馏水把反应釜清洗 2 次, 并用实验气体进行吹扫, 然后抽真空。在反应釜中注入 300 g 左右 的蒸馏水。开启水浴, 调节恒温水浴的温度至实验设定温度后, 开启缓冲罐, 给反应釜充气, 使反应釜中压力 达到实验设定压力。开始进行气体流量的计量, 气体流量计显示为零时表明反应结束。迅速调节水浴至分 解反应需要温度, 待温度稳定后开启放气阀, 直至系统压力为常压, 关闭放气阀。开启通向湿式气体流量计 的气体阀, 记录不同时刻分解得到的气体的体积值, 直至湿式流量计的示数 5 m in 不变为止。
Key words: methane hydrate; kinetics; dissociation; active energy
天然气水合物是由一种或几种烃类气体在一定的温度和压力下和水作用生成的一种非化学计量的笼型 晶体化合物。形成天然气水合物的主要成分是甲烷, 甲烷体积分数超过 99. 9% 的天然气水合物通常称为甲 烷水合物[ 1 3] 。1 体积水合物可以储存 150~ 180 体积标准状态的天然气[ 4 6] , 利用这一性质进行天然气水合 物储存和运输更安全而且更有效率[ 7] 。水合物储运技术主要包括制备、处理、储存运输和分解等过程。由 于气体水合物的分解涉及到气、液、固 3 相, 是一种较为复杂的物理化学过程, 到目前为止, 气体水合物分解 动力学的研究还不够完善。K amat h 等[ 8] 使用加热法首先量化了甲烷和丙烷水合物的热分解速率, 认为气体 水合物的分解受界面传热控制, 和流体泡核沸腾现象具有一定的相似性。U llerich 等[ 9] 也对甲烷水合物的 热分解进行研究, 假定分解过程产生的水直接被甲烷气体携带离开固体表面, 水合物分解可以被认为是一个 移动的界面消融问题。K im 等[ 10] 使用降压法在半连续搅拌釜式反应器中研究了甲烷水合物的减压分解现 象。刘犟等[ 11] 采用降压法研究了封闭体系内, 甲烷水合物在活性炭中的分解动力学。孙长宇等[ 12] 测定了 CH4 水合物在不同温度、压力下的分解动力学数据, 建立了甲烷水合物分解动力学模型。而水合物运输过程 中涉及到常压分解, 但对该过程研究的报道还不多。
( 6)
第 28 卷 第 7 期
田 龙, 等: 甲烷水合物常压分解
25
为表达 nP , 给出假设 5) , 即所有实验中水合物分解具有相同的颗粒直径 D 0( !m) , D 0 不随分解过程而变。
在以上的假设下, 颗粒数可以被计算出, 即
nP =
nH/
H VP =
nH/
H(
D
3 0
/
6)
( 7)
2 分解动力学模型的建立
甲烷水合物分解过程包括: 颗粒表面的晶格破坏; 客体分子( 甲烷) 由表面脱附。以上分解过程由于发生
于固体表面, 而不是固体内部。随着分解过程的进行, 颗粒缩小, 甲烷气体从固体表面释放出来, 产生的气体
随后进入气相主体, 满足
dnH/ d t = k d A s( f e - f )
t= 4
t / min
0 2 8 17 20 23 34 39 42 56
V r/ L
0 8. 171 23. 082 37. 083 41. 857 43. 456 51. 112 52. 783 53. 070 56. 022
收稿日期: 2006 02 13. 基金项目: 国家自然科学基金( 50176051) . 作者简介: 田 龙( 1979 ) , 男 , 硕士生. E mail: tianlong@ ms. giec. ac. cn
24
武 汉理工大学学报
2006 年 7 月
因此, 针对甲烷水合物常压分解试验提出了常压分解动力学模型, 并结合实验数据, 给出了相应的分解 动力学模型及活化能数据, 通过与其他人测定的活化能数据的对比, 验证了模型的可靠性。
摘 要: 在容积为 1 L 的高压反应釜中 制备甲 烷水合 物, 采 用常压 恒温分 解法测定 了甲 烷水合 物分 解过 程的实 验数
据, 并提出了以微分方程表达的宏观分解动力学模型, 同时计算出甲烷水合物分解的活化能 为 76. 93 kJ/ mol。 并与其他
人测定的活化能数据进行了对比, 验证了动力学模型的可靠性。
式( 10) 表明甲烷水合物的分解速率与球形因子、水合物密度和颗粒直径以及推动力有关。
当 t= 0 时, nH= n0, 给出
ln nH/ n0 = - k∀t
( 11)
推动力可表示为
∀f = f e- f
( 12)
式中, 甲烷的 逸度 f 由 SRK 方 程得出; 三 相平 衡压 力 f e 由 甲烷 水合 物生 成经验 公式 lg P = 1. 415+
nH = n0- [ m( t) - m(0)]
nH/ n 0 = 1 - [ m ( t ) - m ( 0) ] / n0
%
令
n0=
∃ m( t),当
t= 0
t=
0 时,
m( t) =
m ( 0) =
0, m ( t ) =
PV r/ RT , 则
nH/ n 0 =
1-
PV r/ RT -
%
m ( 0)
第 28 卷 第 7 期 2006 年 7 月
武汉理工 大学学报
JOURNAL OF WUHAN UNIVERSITY OF TECHNOLOGY
Vol. 28 No. 7 Jul. 2006
甲烷水合物常压分解
田 龙1, 2, 樊栓狮1, 郝文峰1, 3
( 1. 中国科学院广州天然气水合物研究中心, 广州 510640; 2. 中国科学院研究生院, 北京 100039; 3. 大连理工大学化工学院, 大连 116012)
( 1. Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences, Guangzhou 510640, China; 2. Graduate School, Chinese A cademy of Sciences, Beijing 100039, China; 3. School of Chemical Eng ineering, Dalian U niversit y of
将式( 7) 代入式( 6) 得
A s = 6/ H D 0 nH
( 8)
将式( 8) 代入式( 1) 得
- d nH/ dt = k∀nH
( 9)
从式( 9) 可以看出甲烷水合物的常压分解速率与未分解的水合物物质的量成线性关系
k∀ = kd 6/ H D 0 ( f e - f )
( 10)
式中, k∀为表观分解速率常数, m in- 1。
( 1)
式中, nH 为反应器中剩余物质的量( mol) ; k d 为本征动力学反应常数( mol!m- 2!MPa- 1!min- 1) ; A s 为分解
表面积( m2) ; f e 为三相线逸度( Pa) ; f 为气相逸度( Pa) 。
为求解分解速率方程, 作出如下假设: 1) 水合物组成不变, 单位体积的水合物含有的甲烷数量已知, 即 H
no rmal pressure. Exper imental data of the dissociatio n process w er e automatically recorded throug h a data collector. A new kinet ic model for methane hydrate dissociat ion w as proposed, and the activation energy of the reaction, 76. 93 kJ/ mol, w as calculat ed. M oreover, reliability of the kinetics model was validated by comparison with others.
T echnolog y, Dalian 116012, China)
Abstract: M ethane hydrates, w hich were prepared in a high pressure reactor of 1 L , w er e disscociated under isothermal and
表 1 湿式流量计中气体 量随时间的变化
t= 0
t / min
0 2 8 17 20 23 34 39 42 56 62 67 77 84
Vr/ L
0 5. 563 12. 933 24. 734 29. 433 30. 849 39. 954 40. 857 41. 206 46. 663 48. 307 49. 274 54. 538 55. 698
1实验
1. 1 装置 实验装置主要包括反应釜、恒温水浴、温度与压力测
量仪表、流量计和数据采集系统等( 见图 1) 。装置的核心 是高压反应釜, 容积为 1 L , 最大工作压力 20 MPa, 工作 温度范围- 15~ 100 。采用无级调速永磁旋转搅拌装 置, 转速调节范围 0~ 1 500r/ min。反应釜的温度由恒温 水浴控制, 恒温水浴控制精度为 0. 01 。反应釜内的 温度由 2 个 Pt 100 铂电阻测量, 压力调节阀前后的压力由 2 个 0. 25 级精度的压力表测量, 压力调节阀前的压力表 测量范围为 0~ 25 MP a, 压力调节阀后的压力表测量 范 围为 0~ 10 MPa。
( mol/ m3) 已知。2) 反应器中颗粒具有相同的体积 V P( m3) , 同一时刻分解速率相同, 且是分层分解。因此
nH = H nP V P
( 2)
式中, nP 为水合物颗粒数。3) 颗粒随分解过程均匀变化, 且分解过程温度基本不变。4) 颗粒是非球形的, 定 义颗粒具有相同的直径。则
VP = D3/ 6
关键词: 甲烷水合物; 恒压; 分解; 动力学
中图分类号: T Q 026. 7
文献标志码: A
文章编号: 1671 4431( 2006) 07 0023 04
Methane Hydrate Dissociation at Normal Pressure
T IA N L ong1, 2, FA N Shuan shi 1, H A O Wen f eng 1, 3
( 球体积公式)
( 3)
表面积
A P = D2/ =
6 V P / 2/ 3
( 4)
式中, 为颗粒偏心因子。结合方程式( 2) 和式( 4) 给出单个颗粒表面积
AP =
1/ 3 6 nH 2/ 3 HnP
( 5)
对于整个表面积
A s = nP A P = ( nP ) 1/ 3 6nH / H 2/ 3/
0. 041 7 # ( t + 0. 01t 2) 和 SRK 方程法计算得出。
3 结果与讨论
3. 1 分解气体量随时间的变化
Βιβλιοθήκη Baidu
对在反应釜内形成的甲烷水合物进行分解实验研究, 对实验结果进行数据处理。下面对反应器中甲烷
气体和气体体积记录仪进行物料衡算。
反应器中: 初始态为 n 0( mol) , t 时刻为 nH( mol) ; 计量器( 气体体积记录仪) 中: 初始态为 m ( 0) ( mol) , t 时刻为 m ( t ) ( mol) 。