收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

岩石微孔隙中气体吸附、链状分子运移的计算模拟及其油气地质意义

杨侃  
【摘要】:随着常规油气资源勘探难度的增加,煤层气、页岩气等非常规油气资源在能源结构中的地位越来越重要。认识和了解煤层和页岩孔隙中甲烷和二氧化碳的吸附特征及岩石形变规律,对于煤层气、页岩气的资源评价和开发有着十分重要的意义。本文在借鉴前人岩石孔隙表征结果的同时,采用氮气探针气体吸附法分析了4个煤岩样品和11个页岩样品的孔隙结构的特征,发现纳米尺度(50 nm)孔隙较为丰富。鉴于纳米级孔隙所独有的纳米尺寸效应及其对气体吸附的重要作用,本研究针对这类孔隙-气体的相互作用开展了一系列计算模拟工作,淬火固体密度泛函理论(Quenched solid density functional theory, QSDFT)和DPD方法是本文采用的主要研究方法。 首先,本文采用淬火固体密度泛函理论(Quenched solid density functional theory, QSDFT)研究了地质条件下甲烷和二氧化碳在煤层中的吸附行为。研究的重点是吸附过程中温度,压力以及孔径分布对于煤岩变形的影响。研究结果包括,0.5到50nm的煤层孔在298K和360K下100MPa的压力范围内的吸附量,密度分布曲线以及溶剂化压力曲线。研究发现甲烷吸附导致的溶剂化压力对孔径的变化曲线存在两种主要类型:类型Ⅰ表现为在整个压力范围内的单调膨胀。这种类型在最小的孔中1.3σff(0.5nm)中最为典型。这类孔的特点是不能容纳超过一层的甲烷分子。类型Ⅱ表现为在低压下的收缩,而后膨胀。类型Ⅱ在能容纳整数倍层数甲烷分子(2-6层)的孔中最为常见。360K下的二氧化碳吸附导致的溶剂化压力与外压的依存关系也存在两种类型的变形行为:第Ⅰ类形变表现为在全压力范围内的单调膨胀,但在极高压下(100MPa)又会表现出收缩或者继续膨胀。该类行为在最小的孔(0.5nm)中较为典型,该类孔不能容纳超过一层的二氧化碳分子,见0.5nm的密度分布曲线。第Ⅱ类变形行为在低压下表现出收缩,而后膨胀。高压下,该类变形行为也可以出现或膨胀或收缩的现象。第Ⅱa类行为出现在0.6和0.9nm的孔中。所有0.9nm的孔都表现出第Ⅱb的形变。通过与文献中甲烷和二氧化碳吸附导致形变的实验数据的对比,验证了QSDFT模型的有效性。通过对一个假定模型盆地中甲烷和二氧化碳吸附行为的研究,建立了甲烷和二氧化碳吸附量、溶剂化压力与煤层孔结构以及埋藏深度之间的联系,并对二氧化碳对甲烷的驱替行为进行了初步的探讨。研究发现煤层的变形对于埋藏深度有很强的依赖性,在不同的埋藏深度膨胀或者收缩的发生依赖于煤样品孔径分布的情况。通过对地质条件下煤层中二氧化碳驱替甲烷的模拟,探讨了增强的甲烷开采方案(ECBM)中二氧化碳储量以及二氧化碳注入过程中煤层孔隙结构变化的情况。研究发现,微孔在驱替过程中的应变变化比较大,而介孔的变化较小。不同微孔应变达到极值的深度也不尽相同。对于介孔来说,5nm的孔表现出最小的区别(0.6%)。而对于微孔来说,体相应变的差异在1nm孔径100m的情况下可以达到1.7%。这么大的应变差异将有可能导致储层渗透率的降低。 其次,本文采用同样的方法研究了甲烷和二氧化碳在页岩中的吸附储量及由吸附研究的变形行为,相应地获得了甲烷和二氧化碳吸附量、溶剂化压力与页岩孔结构以及埋藏深度之间的联系。深度浅的时候,微孔比介孔表现出更大的单位吸附能力。小于1nm的孔的甲烷吸附能力随深度增加而快速增加,到埋藏为1500m时吸附量达到8 mmol/cm3左右的极大值。而2nm和5nm的孔的甲烷吸附量在到达6000m之前一直在持续增加,其对吸附量的贡献在3000m深度之后超过微孔。本文发现甲烷的吸附导致了三种不同的应变-深度关系。第一种,在整个深度范围内都表现为先膨胀后收缩,例如0.5nm和0.7nm的孔。在达到膨胀极值之前(0.5nm的孔1000m下的0.11%,0.7nm的孔为2000m下的0.12%),体相应变随埋藏深度快速增长,之后收缩。第二种,1nm的孔,在极浅部随深度增加表现为收缩,至100m达到一个极小值后缓慢增长,最后在6000m时膨胀至0.11%。第三种,对于≥2 nm的孔,在整个深度范围内其体积应变的变化≤0.04%,与微孔相比其影响可以忽略不计。在将二氧化碳封存到页岩中时,当深度小于1000m的时候,微孔中的所吸附的二氧化碳占主导地位,但当深度大于1000m时,介孔中的存储将变得更加重要。如果页岩层有大量的≤1nm的孔,那么其在6000m的深度之内的形变都有较大变化,有可能对二氧化碳注入过程中页岩孔结构产生影响。 通过对甲烷和二氧化碳在煤层及页岩中的吸附研究,本文发现由流体分子导致的纳米孔隙的吸附变形行为可能存在一定的规律。为了进一步认识形变与流体分子的吸附行为的关系,本文采用QSDFT的方法定量得研究了表面粗糙程度对流体吸附和颗粒之间的相互作用的影响。表面的不均一性对于固体的润湿性和粘着性质有着很重要的影响。然而,大多数计算固体与流体相互作用的理论和模拟方法都假设存在一个标准的平滑的Gibbs流体固体分界面。这个假设导致了靠近固体表面的流体相出现人为的分层现象。具体的表现就是密度分布曲线的明显波动。这种流体密度的成层性导致了溶剂化压力以离孔壁的距离为函数发生波动,但是溶剂化压力的这种明显的波动性却在试验中极少见到。在QSDFT中,表面粗糙成分采用粗糙度参数来定量表示,该粗糙度参数代表了固体表面松散层的厚度。通过对LJ流体在平行板状孔中吸附(模拟了氮气在石墨中的吸附)的例子,本文发现在饱和状态,溶剂化压力随孔径有规律的波动。这种波动随表面粗糙度的增加(从0增加到1/2倍的分子直径)而快速减小。溶剂化压力的波动性与流体密度在匹配孔中流体密度的最大值以及不匹配孔中流体密度的最小值有很好的对应。后者以松散排列和流体的密度减小为特征。该分析表明分子尺度上的表面粗糙度对局限流体的成层排列以及成层性都有很大的抑制作用。另外,本文计算了在吸附过程中溶剂化压力随外压的变化情况,其变化随孔径和粗糙度的不同可以分为几种不同的机制。由于溶剂化压力决定了多孔介质的吸附变形性质,溶剂化压力随外压变化的非单调性很好的对应了石墨等微孔材料的在不同压力下的收缩膨胀性质。 第三,本文采用三维DPD的模拟方法研究了链状大分子在不同外力场的作用下穿越狭窄孔道时链长的行为,以期对链状分子的运移行为进行深入理解。研究发现外力与平均穿越时间存在着反比关系,链长对穿越时间有很大影响,其在良溶剂下的指数关系为1.56,在不良溶剂下的指数关系为1.76。这与τtrans~N1+v的理论计算结果很接近。模拟观测表明穿越过程中的平均回旋半径主要受其初始状态的平衡值的影响。在穿越过程结束的时候,其指数关系为1,这与其受相同外力在无限空间中的运动行为一致。这可能是受到的外力力场强度大以及带电电荷之间相互作用较强造成的(在带电链状大分子穿越静电力场的情况下)


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 白正华,王先彬;C-O-H体系高温高压流体相组成的热力学研究[J];矿物学报;1998年03期
2 杨荣兴,殷辉安;含流体相变质反应的热力学计算[J];桂林工学院学报;1992年02期
3 柏天宝,A.F.Koster,VanGroos,赵斌;Pb,Zn,Mn,Zr,Ti,Cr,Co,Ni在安山岩熔体与流体相之间分配的初步实验研究[J];科学通报;1995年23期
4 B.J.Wood;J.V.Walther;刘永祥;;变质作用期间的流体流动及流体—岩石比率的应用[J];世界地质;1991年02期
5 李福春,朱金初,漆亮,饶冰,潘根兴;富氟花岗岩体系岩浆流体内稀土元素演化规律的实验研究[J];高校地质学报;2002年01期
6 胡晓燕;尚林波;毕献武;胡瑞忠;樊文苓;陈佑纬;;锡在流体和花岗质硅酸盐熔体间分配行为的实验研究[J];矿物岩石地球化学通报;2007年04期
7 龙国清;邓宏文;朱继田;杨希冰;王立峰;;浅层气藏流体相控属性参数建模及应用[J];海相油气地质;2007年04期
8 彭省临;;铅、锌流-熔分配实验结果及其在矿床成因研究中的应用[J];地质学报;1988年04期
9 刘俊来;上部地壳岩石流动与显微构造演化——天然与实验岩石变形证据[J];地学前缘;2004年04期
10 Hans Keppler;Peter J.Wyllie;张家富;;流体在岩浆活动过程中对铀、钍搬运和分离的作用[J];世界核地质科学;1992年03期
11 Y.Tatsumi ,李增田;俯冲岩石圈释放的流体相的化学特征及岛弧岩浆成因:高压实验和天然岩石的证据[J];地质地球化学;1988年03期
12 柏天宝;;浅析初始样品对元素分配的影响[J];矿物岩石地球化学通报;1988年03期
13 阙梅登;花岗岩浆铝碱饱和程度与矿化的关系[J];地质科技情报;1991年01期
14 K.J.Dunn;林桦;;在流体饱和的多孔圆柱体中的低频声衰减[J];勘探地球物理进展;1987年06期
15 赵国春;CO_2流体相在紫苏花岗岩形成中的作用[J];地质地球化学;1994年06期
16 曾志刚;花岗岩-流体相互作用的实验研究进展[J];地质地球化学;1997年04期
17 Robert D.Stoll;周根祥;;沉积物中的声波[J];水声译丛;1979年01期
18 B.Charoy ,E.Gouzalez-Partida ,思南;墨西哥坦克斯戈-瓜纳华托州块状硫化物和金银矿脉生成时的流体相分析[J];地质地球化学;1985年12期
19 P.J.Pollard ,戴爱华,吴志勤;富氟与富硼含锡体系的差异演化[J];地球与环境;1988年10期
20 魏进家,姜培正,胡春波;气固两相流中气体湍流强度的研究[J];西安交通大学学报;1998年01期
中国重要会议论文全文数据库 前10条
1 李慧;尚林波;王水龙;;铜在熔体流体间分配实验的研究[A];中国矿物岩石地球化学学会第13届学术年会论文集[C];2011年
2 刘俊来;关会梅;;含流体相条件下白云岩的低温流动——破裂-溶解-结晶过程的意义[A];全国包裹体及地质流体学术研讨会论文摘要[C];2002年
3 胡晓燕;毕献武;尚林波;胡瑞忠;樊文苓;;锡在晶体-熔体-流体间分配行为及其对锡成矿的制约[A];中国矿物岩石地球化学学会第11届学术年会论文集[C];2007年
4 张钋;李幼铭;;油气储集层的广义反射透射系数及意义[A];1994年中国地球物理学会第十届学术年会论文集[C];1994年
5 高俊;Reiner Klemd;苏文;熊贤明;;洋壳俯冲释放的富高场强元素(Ti-Nb-Ta)流体——来自西天山榴辉岩相含金红石分凝体和脉体的证据[A];2006年全国岩石学与地球动力学研讨会论文摘要集[C];2006年
6 姜元勇;徐曾和;;局部非热平衡条件下填充床中的可压缩流动与气固反应[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
7 王尚旭;牟永光;;双相介质分界面上波的反射与透射[A];1993年中国地球物理学会第九届学术年会论文集[C];1993年
8 韩其玉;王旭辉;刘仲一;;孔隙介质震电场研究[A];1995年中国地球物理学会第十一届学术年会论文集[C];1995年
9 刘晓峰;解习农;姜涛;张成;;东营凹陷超压含盐流体的特征及其石油地质意义[A];第三届全国沉积学大会论文摘要汇编[C];2004年
10 丛曰祥;;陕西商洛地区宽坪岩群中剪切带与流体包裹体[A];中国地质科学院地质研究所文集(25)[C];1993年
中国博士学位论文全文数据库 前10条
1 李丽妍;方阱流体气液相平衡性质和临界非对称性的计算机模拟研究[D];华东理工大学;2014年
2 王世芳;嵌有类分形树状分叉网络的多孔介质中非牛顿流体渗流特性研究[D];华中科技大学;2012年
3 董伟;尕斯库勒E_3~1油藏相控建模与剩余油分布研究[D];成都理工大学;2008年
4 周斌;水库诱发地震时空演化特征及其动态响应机制研究[D];中国地震局地质研究所;2010年
5 曹秀玲;三河尖矿深井高温热害资源化利用技术[D];中国矿业大学(北京);2010年
6 赵伟;东营凹陷古近系沙四上亚段滩坝砂体固体—流体相互作用与有效储层预测[D];中国地质大学(北京);2011年
7 吴家龙;超临界流体辅助的聚合物系脱挥基础研究[D];浙江大学;2002年
8 贺征;气粒热流场中微粒特性的动力学分析[D];哈尔滨工程大学;2009年
9 杨侃;岩石微孔隙中气体吸附、链状分子运移的计算模拟及其油气地质意义[D];南京大学;2011年
10 蒋国豪;氟、氯对热液钨、铜成矿的制约[D];中国科学院研究生院(地球化学研究所);2004年
中国硕士学位论文全文数据库 前10条
1 李婷;咸阳城区地压型热储流体同位素水文地球化学空间分布及其指示意义[D];长安大学;2012年
2 刘鑫;超临界流体相行为研究及其在涂料工业中的应用[D];北京化工大学;2001年
3 杨敏;留416断块沙三上段低渗透油藏储层建模与非均质评价[D];成都理工大学;2009年
4 杨丹;固体火箭发动机气—固两相流的数值模拟[D];哈尔滨工程大学;2006年
5 张引弟;室内颗粒物行为的数值分析[D];武汉科技大学;2007年
6 刘晓俐;高金属含量固体火箭发动机的相关问题研究[D];哈尔滨工程大学;2006年
7 张峰光;低渗透煤层瓦斯流固耦合理论的研究[D];太原理工大学;2007年
8 于康成;固体火箭发动机潜入喷管背壁区流动研究[D];哈尔滨工程大学;2007年
9 任晔;单裂隙岩体渗流与传热耦合的解析解与参数敏感度分析[D];北京交通大学;2009年
10 杨光;岩样动电耦合系数测量及相关力学问题研究[D];哈尔滨工业大学;2010年
中国重要报纸全文数据库 前2条
1 ;盲人有望拥有“机器眼”[N];计算机世界;2002年
2 庄士鹏;贝根管道 创造新生活模式[N];中国房地产报;2003年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978