深层天然气扩散行为与碳同位素分馏模型研究进展.pdf

第5 4卷 第4期2 0 2 3年7月 太原理工大学学报J OUR NA L O F TA I YUAN UN I V E R S I T Y O F T E CHNO L OG Y V o l.5 4 N o.4 J u l.2 0 2 3 引文格式:董海裕,邹雨,王国建,等.深层天然气扩散行为与碳同位素分馏模型研究进展J.太原理工大学学报,2 0 2 3,5 4(4):5 9 9-6 0 8.D ON G H a i y u,Z OU Y u,WAN G G u o j i a n,e t a l.R e s e a r c h p r o g r e s s o n d i f f u s i o n b e h a v i o r a n d c a r b o n-i s o t o p e f r a c-t i o n a t i o n m o d e l s o f d e e p n a t u r a l g a sJ.J o u r n a l o f T a i y u a n U n i v e r s i t y o f T e c h n o l o g y,2 0 2 3,5 4(4):5 9 9-6 0 8.收稿日期:2 0 2 3-0 1-2 8;修回日期:2 0 2 3-0 4-1 9 基金项目:国家自然科学基金资助项目(4 1 8 0 2 1 9 1);山西省自然科学基金资助项目(2 0 1 9 0 1 D 1 1 1 0 5 0)第一作者:董海裕(1 9 9 7-),硕士研究生,(E-m a i l)d o n g h a i y u 0 7 5 9 l i n k.t y u t.e d u.c n 通信作者:刘东娜(1 9 8 2-),副教授,博士,主要从事煤地质学、显微岩石学及盆地沉积构造演化的研究,(E-m a i l)l i u d o n g n a t y u t.e d u.c n深层天然气扩散行为与碳同位素分馏模型研究进展董海裕1,邹 雨2,王国建2,刘东娜1,周 瑞1(1.太原理工大学 矿业工程学院,太原 0 3 0 0 2 4;2.中国石化石油勘探开发研究院 无锡石油地质研究所,江苏 无锡 2 1 4 1 2 6)摘 要:【目的】深层油气资源逐渐成为我国油气勘探与开发领域的重要目标。为适应发展需求,开展深层天然气扩散过程的可视化研究势在必行。【方法】以天然气分子为出发点,以深层环境为背景,结合多种影响因素,在总结现有扩散方式F i c k扩散、K n u d s e n扩散、表面扩散以及构型扩散的基础上,进行甲烷扩散行为模型及其碳同位素分馏模型的研究进展总结,评价了多种模型在表达深层天然气扩散行为与碳同位素分馏效应方面的适用性。同时认为,建立相应的数学表述模型是准确和全面理解扩散行为机制的关键。【结论】现有的耦合模型不仅能分析深层纳米孔隙中的天然气扩散行为,其高压低孔低渗条件下K n u d s e n扩散还会影响总扩散强度及碳同位素分馏程度。此外,在这些数学模型的基础上,可以结合地层埋藏演化史,进行气藏散失强度计算。研究分析旨在更精准、定量地刻画深层天然气扩散行为,服务于深层油气勘探与开发,对深部气藏保存与破坏评价具有指导意义。关键词:深层天然气;扩散行为;同位素分馏;数学模型中图分类号:T E 1 3 文献标识码:AD O I:1 0.1 6 3 5 5/j.c n k i.i s s n 1 0 0 7-9 4 3 2 t y u t.2 0 2 3.0 4.0 0 3 文章编号:1 0 0 7-9 4 3 2(2 0 2 3)0 4-0 5 9 9-1 0R e s e a r c h P r o g r e s s o n D i f f u s i o n B e h a v i o r a n d C a r b o n-i s o t o p e F r a c t i o n a t i o n M o d e l s o f D e e p N a t u r a l G a sD O N G H a i y u1,Z O U Y u2,WA N G G u o j i a n2,L I U D o n g n a1,Z H O U R u i1(1.C o l l e g e o f M i n i n g E n g i n e e r i n g,T a i y u a n U n i v e r s i t y o f T e c h n o l o g y,T a i y u a n 0 3 0 0 2 4,C h i n a;2.W u x i R e s e a r c h I n s t i t u t e o f P e t r o l e u m G e o l o g y,S i n o p e c P e t r o l e u m E x p l o r a t i o n a n d P r o d u c t i o n R e s e a r c h I n s t i t u t e,W u x i J i a n g s u 2 1 4 1 2 6,C h i n a)A b s t r a c t:【P u r p o s e s】I n C h i n a,t h e e x p l o r a t i o n a n d d e v e l o p m e n t o f d e e p o i l a n d g a s r e s o u r c e s h a v e g a i n e d s i g n i f i c a n t a t t e n t i o n i n r e c e n t y e a r s.I n o r d e r t o m e e t t h e d e m a n d s o f d e v e l o p m e n t,m a n y s c h o l a r s h a v e c o n d u c t e d r e s e a r c h o n t h e p r o c e s s o f d e e p n a t u r a l g a s d i f f u s i o n a n d h a v e m a d e n o t e w o r t h y a d v a n c e m e n t s.【M e t h o d s】A i m g t o e v a l u a t e t h e a p p l i c a b i l i t y o f d e e p n a t u r a l g a s d i f-f u s i o n b e h a v i o r a n d t h e c a r b o n i s o t o p e f r a c t i o n a t i o n e f f e c t,t h e s t u d y i s b a s e d o n a n a n a l y s i s o f F i c k d i f f u s i o n,K n u d s e n d i f f u s i o n,s u r f a c e d i f f u s i o n,a n d c o n f i g u r a t i o n d i f f u s i o n.E s t a b l i s h i n g a n a c c u r a t e m a t h e m a t i c a l m o d e l i s b e l i e v e d c r u c i a l i n c o m p r e h e n s i v e l y u n d e r s t a n d i n g t h e m e c h a n i s m o f d i f f u s i o n b e h a v i o r.【C o n c l u s i o n s】T h e c u r r e n t m o d e l s f o r c o u p l i n g c a n a n a l y z e n o t o n l y t h e b e-h a v i o r o f n a t u r a l g a s d i f f u s i o n i n d e e p n a n o m e t e r p o r e s,b u t a l s o h o w i t a f f e c t s t h e o v e r a l l i n t e n s i-t y o f d i f f u s i o n a n d t h e d e g r e e o f c a r b o n i s o t o p e f r a c t i o n a t i o n i n l o w-p e r m e a b i l i t y c o n d i t i o n s.F u r-t h e r m o r e,t h e s e m a t h e m a t i c a l m o d e l s a l l o w f o r t h e c a l c u l a t i o n o f g a s l o s s i n t e n s i t y w h e n c o m-b i n e d w i t h t h e h i s t o r y o f s t r a t a e v o l u t i o n.T h e o b j e c t i v e o f t h i s r e s e a r c h i s t o p r o v i d e a m o r e p r e-c i s e a n d q u a n t i t a t i v e r e p r e s e n t a t i o n o f d e e p g a s d i f f u s i o n b e h a v i o r.T h i s w i l l a i d i n t h e e x p l o r a t i o n a n d d e v e l o p m e n t o f d e e p o i l a n d g a s r e s e r v e s,a n d a l s o h a v e s i g n i f i c a n t i m p l i c a t i o n s f o r t h e e v a l u-a t i o n o f d e e p g a s r e s e r v o i r s i n t e r m s o f p r e s e r v a t i o n a n d d e s t r u c t i o n.K e y w o r d s:D e e p n a t u r a l g a s;D i f f u s i o n b e h a v i o r;I s o t o p e f r a c t i o n a t i o n;M a t h e m a t i c a l m o d e l 目前,随着浅至中层地层中大型油气藏的勘探难度逐渐增大,我国的油气勘探领域已经逐步转向深层和超深层资源的探寻1-4。随着地质理论和开发技术的发展,在塔里木盆地、四川盆地及周缘等地区深层-超深层油气勘探中取得的一系列重大突破5-8,近年来中国深层油气的新增储量呈现增长趋势9。对于深层地质背景中的天然气藏,由浓度驱使的甲烷扩散是重要的气体运移机制,其具有易发生性和和持续作用的特征,不仅影响着气藏的保存与散失评价1 0-1 2,还密切关系着深层油气勘探开发决策与部署1 3-1 5。深层油气藏通常表现为层系多、非均质性强、高温高压、以及微纳米孔隙发育的低孔低渗1,例如四川盆地及周缘奥陶系-志留系的深层页岩气藏,主要以埋深大于3 5 0 0 m,高温高压且变化范围广(6 0 1 2 0,3 51 0 0 MP a),以及基质渗透率低(0.1 11 0-30.0 91 0-3 m2)为特征1 6。在这种超临界状态下,某一空间位置不存在释放压力的过程,天然气扩散引起的运移是一个自发的物理过程1 7,即从高浓度区运移到低浓度区,于是多尺度的微纳米孔隙依旧能够提供天然气储集和运移通道1 8-1 9。因此,微纳米孔隙甲烷气体扩散行为的数学模型(扩散行为模型)研究是对天然气藏准确模拟的基础,有助于了解深层的油气成藏动力学。在过往的相关研究和建模工作中,学者们采用了以下方法:一种是分别针对多种扩散行为(如表面扩散、F i c k扩散、K n u d s e n扩散等)进行逐一分解模拟与研究2 0;另一种则是构建了能够描述总体扩散行为的耦合数学模型1 2-1 3,2 1。另外,甲烷碳同位素分馏特征也可以作为一种重要的地球化学示踪指标,不仅与气源成熟度、生烃史、生物化学作用等油气形成过程有关2 2,还与分子运移的物理过程密切相关2 3-2 4。尤其在近些年来,前人建立了并在持续完善的甲烷碳同位素分馏动力学模型,明确了天然气解析/生产过程中存在显著的碳同位素分馏2 5,并可用于有效区分不同赋存状态气体(吸附气、游离气)运移过程及阶段2 6-2 7。虽然该分馏模型尚未广泛用于描述深层天然气运移,其中分子扩散的同位素分馏研究也处于起步阶段2 8,但其同位素分馏模型可以潜在的与扩散行为模型 相 结 合,也 是 阐 明 耦 合 扩 散 过 程 的 关 键 途径2 9。因此,本文以天然气分子为出发点,以深层环境为背景,结合多种影响因素,进行甲烷扩散行为模型及其碳同位素分馏模型的研究进展总结,旨在更精准、定量地刻画深层天然气扩散行为,服务于深层油气勘探与开发。1 深层天然气扩散机制与测定我国东部及西部地区分别以3 5 0 0 m和4 5 0 0 m埋深界定深层环境3 0-3 1,在深层环境中,孔隙直径一般处于纳米级别,而天然气的扩散行为与努森数K n密切相关3 2。K n定义为气体分子平均自由程(m)与孔隙直径d(m)之比,见公式(1):K n=d.(1)基于分子行为特征以及K n值对扩散行为进行详细划分1 2,2 0如图1所示。1)F i c k扩散(K n值较小,分子平均自由程小于孔隙直径),或称为黏性流动,分子扩散主要表现为分子之间的碰撞,通常发生在大孔隙和高压条件下;2)K n u d s e n扩散(K n值较大,分子平均自由程大于孔隙直径),分子扩散主要表现为分子与孔隙壁面之间的碰撞;3)表面扩散,驱动力是化学势梯度,即吸附层的流体分子通过在吸附势场的作用下发生的跳跃和移动,对于强吸附剂或小孔,该扩散行为显著;4)构型扩散,吸附气体从有机质孔壁解吸后,有机质中溶解气体存在浓度梯度而发生的扩散。需要指出的是,根据K n值来区分F i c k和K n u d s e n扩散存在不同的界定标准,主要有两种划分方法:1)以K n1为标准划分出F i c k和K n u d s e n扩散1 3;2)以K n1 0为标准划分出F i c k和K n u d s e n扩散,而0.1K n1 0为介于两种扩散行为中间的过渡扩散3 2。不仅如此,也有学者提出,将F i c k和K n u d-006太 原 理 工 大 学 学 报 第5 4卷 s e n扩散行为视为一种整体的浓度扩散(用于区别于表面扩散),称为体相扩散1 2,3 3。Fick 扩散（Kn或 0.1）Knudsen 扩散（Kn或 10）表面扩散构型扩散图1 深层天然气扩散行为特征F i g.1 C h a r a c t e r i s t i c s o f d e e p g a s d i f f u s i o n b e h a v i o r扩散行为的作用强度通常使用扩散系数(m2/s)进行表征。在深层环境中,高温会略微增强所有扩散作用,但扩散系数不会产生数量级变化,而随着压力的增大,扩散系数趋于减小3 4-3 5,这可能是由于气体分子的热运动与温度有关3 6;然而随着压力的增大,即使是高温环境下,扩散系数依旧趋于减小1 1,1 3。不仅如此,低孔低渗条件在限制扩散作用的同时,还会显著限制总扩散量,与气藏保存与破坏密切相关,有待结合数学模型进行深入探讨。另外,在微纳米孔隙中普遍存在构型扩散与表面扩散,但是相比较扩散系数较高的体相扩散(通常为1 0-1 41 0-5 m2/s),构型扩散系数非常低(甚至低至1 0-2 2 m2/s)2 0,因此数学模型研究中通常没有考虑构型扩散作用(见下文)。目前可用于原位模拟天然气扩散及测定扩散系数的实验方法主要有三种:解析法、等压扩散法和压力衰减法。解析法2 6基于均质和均匀孔隙结构的假设,用于模拟油气、煤层气开发阶段的天然气扩散行为。等压扩散法1 1,1 8通过在密闭容器中施加压力与温度,模拟甲烷在岩石中的扩散过程,是最接近实际深层地质环境的方法。压力衰减法2 0通过对岩心注入甲烷并追踪其扩散行为来测定扩散系数,并能区分不同的扩散行为。这三种方法各有优缺点,适用于不同的场景。具体方法及适用范围如表1所示。表1 天然气扩散原位模拟实验方法与扩散系数测定T a b l e 1 E x p e r i m e n t a l m e t h o d s a n d d e t e r m i n a t i o n o f d i f f u s i o n c o e f f i c i e n t s f o r i n s i t u s i m u l a t i o n o f n a t u r a l g a s d i f f u s i o n实验方法扩散系数测定公式适用范围解析法单向模型MtM=1-62n=11n2e x p-D n22tr2p .双向模型MtM=33+1-62n=11n2e x p-Dan22tR2a +31-62n=11n2e x p-Din22tR2i .适用于开发阶段中吸附解吸过程研 究,不 适 于 深层扩散行为研究等压扩散法D=l nC0Ci A(V1+V2)t)LV1V2.适用于深层扩散行为研究压力衰减法MtM=821-e x p-D 24h2t .适用于不同扩散行为的解耦研究 注:1)解析法测定公式中的参数见文献3 7-3 8;2)等压扩散和压力衰减法中,D为扩散系数,m2/s;A为岩心横截面积,m2;L为岩心长度,m;t为时间;C0和 Ci分别为在初始时间和某一时间(岩心两端高压容器中的)浓度差;V1和V2分别为等压扩散实验中岩心两端高压容器的体积,m3;Mt和M分别为某一时间和总的甲烷扩散量,k g;h为构型扩散深度。2 体相扩散模型2.1 体相扩散的分子模型通过气体分子的性质来描述气体传输机制,称为分子模型,其主要是针对纳米多孔材料中的体相扩散行为,强调了K n u d s e n扩散的作用,可以准确模拟纳米孔中的各种物理机制3 9-4 0。对于三维孔隙中气体的体相扩散,分子位移的变化与时间成正比,这里以C o p p e n s和D a mm e r s模型为例4 0,可以由式(2)和式(3)进行表述:l i mt=6D t.(2)S()2-df.(3)r（3），t（3）r（2），t（2）r（n），t（n）r（1），t（1）（a）（b）tw（1）r（1），t（1）r（2），t（2）tw（0）t（n-1）tw（2）图2 分子模型中的体相扩散特征(修改自文献4 0)F i g.2 B u l k p h a s e d i f f u s i o n c h a r a c t e r i s t i c s i n t h e m o l e c u l a r m o d e l(m o d i f i e d f r o m l i t e r a t u r e4 0)106 第4期 董海裕,等:深层天然气扩散行为与碳同位素分馏模型研究进展式中:r(t)为从任意原点开始,作为时间t函数的分子位置;D为扩散系数,m2/s;df为孔隙表面的分形维数,介于2(有效平滑)和3(空间填充)之间。其中许多孔隙具有粗糙的内表面,在有限的孔隙空间内,扩散分 子 可 到 达 的 表 面 积 取 决 于 其 有 效 直 径(m),数学表述遵循幂定律,即公式(3).图2(a)示意性地展示了分子模型中的体相扩散,即嵌入多孔介质中的不均匀性。当粗糙度(统计上)自相似时,从分子尺度上,可以导出体相扩散系数的表达式3 9-4 0:D()=()df-2D0.(4)其中D0是光滑孔隙内壁中的扩散系数,其孔隙直径与原始粗糙孔隙相同。表面粗糙度导致体相扩散系数随着气体分子尺寸的增加而单调下降,因为较小的分子可以进入较大的表面积,并且更容易被粗糙面困住(如图2(a).考虑到沿孔壁粗糙面与光滑孔隙的对应关系,气体分子轨迹在光滑孔隙中需要相应延长。因此,可以将粗糙孔隙中气体扩散转化为光滑孔隙中的扩散(如图2(b),其中气体分子在粗糙孔隙某一空间位置及时间以r(i),t(i)表示,而在光滑孔隙中某一空间位置及时间以r(i),t(i)表示,而在孔壁捕获等待时间为t(i)w.因此,基于图2示意模型以及上述理论,可以基于计算机数学模拟进行微观拟合,能够对体相扩散进行精确的模型计算。但是,分子模型更多是分析了分子与孔壁的碰撞,没有考虑有机质孔壁对分子强吸附的影响,忽视了温压作用的影响4 0;并且,在实际天然气扩散研究中,该模型使用的分子建模技术需要大量的计算机资源和计算时间,因此其应用受到限制。2.2 体相扩散的宏观模型描述页岩纳米孔内气体运移机理的宏观模型有两种,第一种是在连续介质模型中修改获得的水动力模型1 5,4 1:该模型结合了含气储层中甲烷滞留和运移的相关机制,用于确定天然气的固有渗透和扩散作用,是根据统一的H a g e n-P o i s e u i l l e型方程,考虑了各种流动状态、气体特性和岩石物性,并基于非线性压力方程的数值解来实现的,如式(5)和式(6):u=-(p)+K(p)+f(p)-(p)D p x.(5)D=K(p)+(p)+(1-)Ka(p)-(1-)(qp(p).(6)式中:U为表观通量矢量,m/s;D为水动力模型扩散系数,m2/s;、K、f、分别为流体密度、固有渗透率、流动条件函数、流体黏度、以及孔隙的等温压缩系数,P a-1;p为气体压力,P a;q为每固体体积吸附的气体质量,k g/m3;x为水平流动方向上的距离,m;为孔隙度,%;为气体动态黏度,P as.该水动力模型通过对气体流体扩散运移进行微观表征来测定天然气扩散系数,并从模拟实验与模型中获取各种重要参数,与流体模拟实验具有更好的吻合度,并考虑了稀疏和滑移效应,可以解释压力与扩散系数密切关系,能够表达纳米孔隙中的气体传输过程1 5,4 1。另一种宏观模型为权重模型,是通过使用一定的权重系数来组合不同传递机制形成的1 2。相比较其他体相扩散的模型(例如分子模型、水动力模型等),权重模型的优势在于结合各种气体运移理论的优势,进行权重关系的计算,来获取更符合实际情况的计算模型。在前 人研究的基 础 上4 2-4 5,WU e t a l1 2和Z HONG e t a l1 3提出了一种合理可信的权重模型来计算F i c k扩散和K n u d s e n扩散的加权系数,并结合温度、压力、气体分子性质、孔隙结构、岩性等因素构建了 较为完善的 体 相 扩 散 系 数 计 算模型:Db u l k=bT2 2P2r+bT2 2P 2r38R T M+2r2r+bT2 2P bT3g.(7)式中,b为B o l t z m a n n常数,取值1.3 8 0 51 0-2 3 J/K;T为绝对温度,K;为气体分子碰撞直径,m;R为通用气体常数,取值8.3 1 4 J/m o l/K;M为气体摩尔质量,k g/m o l;g为气体黏度,P as.在该模型中,各种参数都可以通过实验测试获得,不存在模型拟合的复杂过程,能够进行直观的计算获得气体分子的体相扩散系数,便于推广至不同地质背景下的深层天然气扩散研究中。3 表面扩散模型3.1 表面扩散的跳跃模型跳跃模型是指固体孔壁吸附气体分子从吸附点206太 原 理 工 大 学 学 报 第5 4卷 跳跃到相邻吸附点,表现为活化的吸附气体分子,如图3所示。当被吸附的气体分子能够获得足够的能量从表面逃逸时,它就返回到气相;当跃迁到邻近的吸附位置,则发生表面扩散,适用于单层吸附的气体表面扩散3 3。需要注意的是,气体吸附/解吸是气体分子的相变,但对气体运移没有贡献。由于表面扩散的驱动力是化学势梯度,当假设不存在气体覆盖,吸附气体在表面的扩散系数可表示为3 3:D0s u r f a c e=8.2 91 0-7 T0.5 e x p-H0.8R T .(8)式中:D0s u r f a c e为不存在气体覆盖情况下的表面扩散,m2/s;H为等位吸附热,J/m o l.为了表征高压条件下的表面扩散,则需要考虑气体覆盖度的影响,利用动力学的方法可以推导出真实的表面扩散系数4 6:Ds u r f a c e=D0s u r f a c e(1-)+2(2-)+H(1-)(1-)22(1-+2)2.(9)式中:Ds u r f a c e为表面扩散系数,m2/s;为阻塞速度系数与正向速度系数 之比,无量 纲;H(1-)为H e a v i s i d e函数;为吸附气体覆盖率,可用L a n g-m u i r等温理论计算,%.该跳跃模型的计算结果能够与实验数据进行相互验证,证实了其可靠性,可以应用于纳米孔隙中吸附气体的表面扩散。空吸附位图3 表面扩散的跳跃模型F i g.3 J u m p m o d e l o f s u r f a c e d i f f u s i o n3.2 表面扩散的流体力学模型流体力学模型是假设吸附流体的表面扩散是由表面膜的黏性流动引起的,该理论仅适用于多层吸附的流体表面扩散。KA I NOUR G I AK I S e t a l4 7基于流体动力学,提出了适用于表面扩散的均质流体动力学模型,该模型是对非均质表面的扩展,可以表述为式(1 0)所示:PR1-VSVS C.(1 0)式中:PR为是气体相对渗透率,%;VS为吸附气体的体积,m3;VS C为系统渗滤阈值下的吸附气体体积,m3.该模型还考虑了各种孔隙连通性的三维规则晶格,如式(1 1)和式(1 2)4 7:PR(fb-fb c)t0.(1 1)PRL-tv.(1 2)式中:t0为网络维度的通用临界指数,对于三维晶格t0为1.8 7;fb为 进 行 毛 细 管 冷 凝 的 孔 隙 分数,%;fb c为渗流阈值,%;L为三维孔隙网络的大小,m.在使用不同组的孔隙材料和吸附质进行气体相对扩散实验,当在渗流阈值附近考虑时,无论多孔隙介质的拓扑结构和吸附质的性质如何,都表明扩散、吸附解吸行为普遍存在4 7。由于该模型只能够应用于多层吸附的流体表面扩散行为,因此在后续的推广应用中受到限制。4 耦合的扩散行为模型在深层页岩孔隙,因为受到孔隙形状、孔径大小、介质性质等因素的影响,不同类型的扩散行为(如体相扩散和表面扩散)往往是相互耦合。考虑不同类型扩散耦合的模型对于理解物质在孔隙介质中的传输机理和实现预测具有重要的意义。F i c k i a n理论认为表面扩散与体相扩散对气体的传输量是相互独立的1 2,因此深层天然气的扩散行为可以用耦合的扩散行为模型进行诠释,该模型综合考虑了体相扩散和表面扩散的影响。通过耦合的扩散行为模型,天然气在孔隙中总的扩散系数(Dt o t a l)可以简单的表述为式(1 3):Dt o t a l=Db u l k+Ds u r f a c e.(1 3)在深层环境中,孔隙通常为毛细、弯曲的,且在压力作用下会发生形变,会阻碍扩散行为。除了上述的分子模型,其他模型都需要进行相应的修正,可以通过式(1 4)和式(1 5)进行表述3 3,4 8:m b=.(1 4)r=r01+Kp-Lpp+pL.(1 5)式中:m b为纳米孔隙中气体扩散的修正系数,无量纲;为迂曲度,无量纲;r0和r分别为常压和施加压力下的孔隙半径,m;K为样品体积模量,P a;L为L a n g m u i r体积应变常数;pL为L a n g m u i r压力,P a;为B i o t系数,无量纲;p为孔隙与气体压力差,P a.不仅如此,气藏中存在吸附引起的溶胀响应,因此还需要对表面扩散进行近一步的修正,如式(1 6)所示3 3:306 第4期 董海裕,等:深层天然气扩散行为与碳同位素分馏模型研究进展m s=m b1-dMr -2-1 .(1 6)式中:m s为纳米孔隙中气体表面扩散的修正系数,无量纲;dM为气体分子直径,m.最终,将修正系数引入公式(1 3)中,可以获得修正后的总扩散系数(Dt o t a l):Dt o t a l=m bDb u l k+m sDs u r f a c e.(1 7)总扩散系数的计算方法为耦合的扩散行为模型,WU e t a l1 2基于各种材料及岩石(例如玻璃、活性炭、石墨烯、氧化铝、沸石等)纳米孔隙的扩散模拟实验,证实了耦合的扩散行为模型的可行性,适用于深层天然气扩散研究中。在岩石样品分析中,邹雨等1 1、Z HONG e t a l1 3也实现了该模型在页岩气分析中的应用。然而,在实际的深层气藏中,有机质的分布、孔隙大小及其非均质性、孔隙连通性在模型中缺乏考虑;另外,上述所有模型都是基于只有一种气体分子,且没有考虑实际深层孔隙中的水分或其他烃类,都需要进行相关研究与建立模型2 1,分子动力学模型研究尚待深入3 2。5 分子碳同位素分馏模型在深层天然气中,以浓度扩散为主,所以不存在压力驱动的渗流分馏,虽然在纳米孔隙中扩散过程中,也存在吸附解吸过程,但是在某一空间位置,不存在释放压力以及物性环境的显著变化,吸附解吸过程处于平衡状态,因此使用菲克第二定律,即公式(1 8)就能完整描述深层天然气扩散行为中的碳同位素分馏。随着天然气在纳米孔隙中扩散行为的推进,1 2C和1 3C甲烷分子在时间i的累积体积(Q)可以表示为公式(1 9),瞬时甲烷碳同位素1 3CCH4可以表示为公式(2 0)2 7:C#k t=-1xm xxmD#C#k x ,C*k t=-1xm xxmD*C*k x .(1 8)式中:#和*上标分别代表1 2C和1 3C甲烷分子;t为脱气时间,s;x为距离,m;m是几何因子(对于平面板材,m=0;对于圆柱体,m=1;对于球体,m=2);为吸附气体的覆盖率,%;D为扩散系数(包含体相和表面扩散),m2/s;Ck是干酪根结构孔隙内的气体浓度,m o l/m3.Q#(ti)=Skx=Hkx=0C#k0-C#k(x,ti)dxVm,Q*(ti)=Skx=Hkx=0C*k0-C*k(x,ti)dxVm,Q(ti)=Q#(ti)+Q*(ti).(1 9)1 3 CC H4(ti)=Q*(ti)-Q*(ti-1)Q#(ti)-Q#(ti-1)Rs t-1 1 0 0 0.(2 0)式中:Q为气体运移体积,L;Sk是干酪根中的有效开放面积,m2;Hk是干酪根的平均厚度,m.在该分子碳同位素分馏模型中,甲烷分子扩散系数(D)的获得是关键一步,可以基于该分馏模型进行建模拟合获得。首先,设置合适的初始和边界条件,包括基质或有机质孔隙中气体流动-扩散的初始和边界条件,例如初始压力p0(P a)、初始同位素值1 3C0()、饱和浓度Ck0(%)等;然后,选择并优化模型参数将其代入同位素分馏模型,以获得基质孔隙中压力或气体浓度的时空分布;最终当Q和Qt e s t(实验测试),以及1 3CC H4和1 3CC H4-t e s t(实验测试)之间的误差之和最小时,可以推导出目标参数(例如扩散系数),这其中需要使用计算机进行拟合,可以采用最小二乘法进行参数优化2 7。拟合获得的扩散系数需要实验测试和模拟,并需要大量计算机资源。因此,可以将上述扩散模型计算得到的扩散系数引入分馏模型中,从纳米孔隙扩散行为出发,建立新的耦合模型。6 现有模型分析与深层天然气保存评价应用 通 过 构 建 数 学 模 型,可 以 认 为F i c k扩 散、K n u d s e n扩散和表面扩散不是相互独立的,它们之间的关系由气体分子的平均自由程和孔径决定。在深层天然气藏中,任何一种单一扩散行为的改变都将影响整个扩散过程,同时深层条件的变化也可以通过该模型进行理论研究。基于分子模型4 0的研究强调了分子与孔壁碰撞的过程,突出了K n u d s e n扩散行为对气体运移的重要贡献。使用不同的分子模型进行计算机模拟,如连续的N a v i e r-S t o k e s方程和原子分子动力学方程建立的模型以及分子动力学建立的二元分子混合物模型4 9-5 0,都显示了纳米孔隙中分子扩散对气体运移的积极效果。宏观模型简化了以往繁琐的计算过程,尤其是对体相扩散计算的权重模型,能够更好地契合深层高温高压、低孔低渗环境。随后,WU e t a l1 2、Z HONG e t a l1 3基于权重模型揭示了有关扩散行为机制,提出了高压条件对应着低水平的扩散作用,这可能与高压条件下K n u d s e n扩散受到显著406太 原 理 工 大 学 学 报 第5 4卷 限制有关,如图4所示;而在相对较低的有效应力条件下,由于F i c k扩散和K n u d s e n扩散的贡献增加,甲烷的总扩散系数增大。对于表面扩散模型,许多学者已经证明,由跳跃模型导出的理论更适合于纳米孔隙中的表面扩散5 1-5 2,并且与体相扩散可以进行线性叠加(F i c k i a n理论)形成耦合的扩散行为模型1 2,最终纳米孔隙中甲烷分子的扩散行为如图4所示。rrrr增压增大压力引起的压缩形变Kundsen 扩散；Fick 扩散；13CH4分子；13CH4分子；受阻的 Kundsen 扩散；表面扩散图4 深层天然气甲烷分子的扩散行为模式图F i g.4 D i f f u s i o n b e h a v i o r p a t t e r n o f d e e p n a t u r a l g a s m e t h a n e m o l e c u l e s 基于分子碳同位素分馏模型,相关研究表明在深层环境中高压是影响同位素分馏最重要的因素(当 压 力 稳 定 不 变 时),同 位 素 分 馏 受 到 显 著 限制2 5。而只有在地层释放压力、开发过程中,非稳态下(气体压力浓度随时间变化)吸附/解吸会导致同位素显著的分馏2 7。S HE N e t a l2 8通过模型分析,研究结果表明深层环境中努森数K n与分馏正相关,即高压下K n u d s e n扩散受到限制,相应的K n值变小,产生的同位素分馏也较小,于是扩散行为中同位素分馏特征依旧可以使用图4进行表征。通过上述模型的建立与相关理论的研究,可以将获得扩散系数或同位素分馏应用至深层天然气保存评价中1 1,2 8。深层天然气的扩散系数(D)对应着其在深层储层中的迁移能力,这与盖层封闭能力密切相关5 3,5 4。尽管天然气扩散速度很慢,但在漫长的地质历史中,扩散过程导致的气藏损失也可能相当大5 5。当把天然气扩散运移的对象假设为相同的孔隙材料,累积耗散强度(QD/S)可用于评估气体损失1 1,5 6:QD/S=t0DSdCdZdtS=t0DdCdZdt.(2 1)DW=D7 2.1 8-1.0 7 1 6.(2 2)DW t=TWT0W 0DW.(2 3)式中:t为扩散时间,s;QD为在时间t扩散引起的甲烷散失的累积体积,m3;S为天然气藏储层面积,k m2;dC/dZ为浓度梯度,m-1;DW为饱和水岩样扩散系数,m2/s;T和T0分别为地层和实验温度,K;W和W 0分别为地