页岩基质解吸-扩散-渗流耦合实验及数学模型

页岩基质解吸-扩散-渗流耦合实验及数学模型
郭为;胡志明;左罗;高树生;于荣泽;曾博
【摘要】采用川南地区龙马溪组页岩样品,设计了页岩基质解吸-扩散-渗流耦合物理模拟实验,揭示了页岩基质气体流动特征以及压力传播规律.推导了页岩气解吸-扩散-渗流耦合数学模型并且利用有限差分法对数学模型进行数值求解,与实验结果相比较表明该数学模型能够很好地描述气体在页岩基质中的流动规律.同时对页岩基质气体流动的影响因素进行了分析,认为页岩基质的渗透率、扩散系数、解吸附常数等因素均能影响页岩基质气体的流量和压力传播规律,在页岩气藏的开发过程中需要考虑这些参数的影响,该数学模型为页岩气井产能计算提供了更准确的计算方法.
【期刊名称】《力学学报》
【年(卷),期】2015(047)006
【总页数】7页(P916-922)
【关键词】页岩基质;解吸附;扩散;渗流;实验;数学模型
【作者】郭为;胡志明;左罗;高树生;于荣泽;曾博
【作者单位】中国石油勘探开发研究院廊坊分院,河北廊坊065007;中国石油勘探开发研究院廊坊分院,河北廊坊065007;中国石油勘探开发研究院廊坊分院,河北廊坊065007;中国石油勘探开发研究院廊坊分院,河北廊坊065007;中国石油勘探开发研究院廊坊分院,河北廊坊065007;国家能源局页岩气研发(实验)中心,河北廊坊065007;中国石油勘探开发研究院廊坊分院,河北廊坊065007;国家能源局页岩气研发(实验)中心,河北廊坊065007
【正文语种】中文
【中图分类】TE132.2
页岩气是指赋存于以富有机质页岩为主的储集岩系中的非常规天然气[1].据预测，
世界页岩气资源量为456.24×1012m3，主要分布于北美、中亚和东亚的中国、
中东和北非、拉丁美洲以及欧洲的前苏联等地区，其中，北美的页岩气储量为108.79×1012m3，占总资源量的23.84%.中国南方海相页岩地层是页岩气主要富集地区，此外，松辽、鄂尔多斯、吐哈、准噶尔等陆相沉积盆地页岩地层也具有页岩气富集的基础和条件[2].中国页岩气技术可采资源量估算约为26×1012m3.页岩气的渗流机理研究是指导页岩气产量预测的重要依据，因此页岩气的渗流机理是页岩气开发的基础.国外学者针对页岩的渗流机理进行了一定量的研究且取得了一些
认识，目前普遍认为页岩气流动过程为：首先是裂缝中的气体被采出，裂缝和基质的压力差使得基质表面的吸附气体解吸然后通过裂缝网络流向井筒，该过程气体的流动遵循达西定律[37]；基质内部气体在浓度差作用下扩散到表面后通过裂缝流入井筒，扩散过程用Fick定律描述[814].页岩气的基质含有大量的纳米级孔隙，孔
隙里面游离气占着一定的比例，气体从基质到裂缝中的流动由Fick定律描述欠妥，除此之外，页岩气流动的物理模拟实验目前相关文献相对较少，本研究针对这些问题，以物理模拟实验为基础研究气体在页岩基质中的流动规律.
1.1 实验样品和实验方法
页岩样品取自川南地区龙马溪组，取样深度1408m.实验设备采用中石油勘探开发研究院廊坊分院自主研发的页岩气解吸--扩散--渗流耦合实验装置，实验流程图见图1.实验步骤如下：(1)所取样品均为现场取芯页岩样品，页岩岩心取上来后用吸
水纸擦拭干净并测量长度、直径和重量，然后放入全直径岩心夹持器中.(2)给岩心
夹持器加轴压和围压，压力均为40MPa；(3)打开阀门1和2给岩心夹持器的两端
均接入装有甲烷气体的中间容器来饱和岩心，中间容器用ISCO泵稳定气体压力
29MPa，饱和岩心的目的是为了恢复地层条件下页岩气的存储状态；饱和岩心的
时间为240d，此时给甲烷气瓶提供压力源的ISCO泵液量3d内不再变化，可以
认为页岩岩心已经接近完全饱和状态.(4)岩心饱和完成后关闭阀门1和2并且撤掉高压气瓶，开始解吸--扩散--渗流耦合实验，设定出口压力为2MPa，打开阀门2，同时计量岩心夹持器两端的压力、气体流量和时间.(5)实验进行37d后为了研究不同工作制度下对气体流动规律的影响，撤掉回压阀，出口端的压力变为大气压.
1.2 实验结果
图2为入口压力和出口压力随着时间变化的关系曲线，从图中可以看出开始产气后，页岩基质中压力传播速度很慢，100d后压力波及到入口边界，压力传播到入口端后，入口端的压力下降非常缓慢.图3为日产气量与时间的变化关系曲线，实
验初期气体的产量大并且递减速度非常快，而在实验后期气体的产量小但递减速度非常缓慢，这与页岩气井的实际生产中产量递减规律一致(图4).
考虑均质的页岩基质，页岩基质中气体的流动包括渗流项以及扩散项，因此其总流动速度vT是达西流动速度vD和扩散流动速度vK之和
达西流动项可以用达西定理描述
k为基质渗透率，m2，µ为气体黏度，Pa·s，p为气体压力，Pa.扩散项可以用Kn 扩散定律描述
D为 Kn扩散系数，m2/s；C为气体的质量浓度，kg/m3；ρg为气体密度，
kg/m3.页岩气以吸附和游离两种状态存在于页岩储层中，以一维情况为例，在页
岩气储层中，假设页岩体积单元为Vr，m3，页岩岩石的密度为ρr，kg/m3，页
岩储层中气体符合真实气体状态方程，气体摩尔分子量Mg，mol/kg，则游离气
和吸附气在页岩体积单元中的质量浓度分别为
式中，Cfre，Cads分别为游离气和参与解吸吸附气的质量浓度，kg/m3；mfre，
mads分别为游离气和吸附气的质量，kg；Vr，Vg分别为页岩储层单元体积和气体体积，m3；ρr页岩储层单元岩石密度，kg/m3；ρa为标况条件下气体密度，kg/m3；Vdes页岩参与解吸吸附气量，m3/kg；SW为含水饱和度，无量纲量；φ为页岩孔隙度，无量纲量.
国内外学者普遍认为页岩的等温解吸曲线与等温吸附曲线是可逆的，并且使用Langmuir模型进行解吸曲线的计算；郭为等2012年利用大样量页岩等温吸附解吸实验得出了页岩等温吸附曲线和解吸曲线不重合，且解吸曲线应该用解吸式进行计算的结论，因此页岩解吸附量Vdes的计算按照解吸式模型进行计算[15].
式中，Vd,pd和c均为解吸常数，Vd为页岩解吸过程样最大吸附容量，m3/kg；pd为吸附速度、解吸速度与吸附热综合函数，Pa−1；c为匮乏压力下的残余吸附量，m3/kg.
气体的密度可以由真实气体的状态方程得到
将式(2)，式(3)代入式(1)可得
气体流动的连续性方程为
qd表示参与解吸的总吸附气量，Q表示源汇项，表示气体的产出或者流入
Vstd为标准状态下气体的摩尔体积，22.4×10−3m3/mol.
将式(4)～式(9),式(11)～式(13)代入式(10)，可得
考虑一维平面流动，方程(14)可以变为
方程(15)只能进行数值计算，利用有限差分法进行差分.考虑时间向前差分的显式差分，方程左边第1项进行离散
令
将方程左边第2项进行离散
方程右边第1项离散
方程右边第1项进行离散，令
将式(16)～式(23)代入式(15)就能得到时间向前差分的有限差分方程.从图 3可看到入口处压力变化非常缓慢，并且在100d后压力传播到岩心的另一端，因此在计算中假设为定压边界，压力为29MPa；产气端压力的边界条件和初始条件为：t=0d 时，出口压力为29MPa；t≤37d时，产气端压力为2MPa；t>37d时，出口压力为0.1MPa；岩心直径为10.5cm，长度L为14.2cm，将岩心按步长∆χ=0.02cm 进行等分，时间步长为∆t=2s，出口处第一个网格的流量就是气体的产量Q.
利用“Matlab”软件进行编程计算，可以求得日产气量和压力随着时间的变化关系，表1为编程计算时各个物理量取值，表格中的扩散系数D为拟合值，其他各个参数均为平行样品实验测量结果.
图5为日产气量随着时间变化的数值计算结果和实验结果对比图，从图中可以看到数值计算结果和实验结果比较吻合；初始阶段，日产气量下降快，到了后期日产气量下降变缓.图6为压力传播随着时间变化的关系图，从图中可以得知，压力的传播过程非常缓慢，100d时压力传播到入口处，这与实验结果一致，靠近出口端压力下降幅度大，远离出口端压力下降幅度小.
3.1 吸附气体的影响
为了对比研究各个影响因素对气体产量的影响，研究对象依然为实验所用的岩心，在计算不同参数对气体产量的影响时，产气端的压力恒定设为0.1MPa.吸附气体占页岩气很大一部分比例，吸附气体如何影响页岩气的产气量和压力传播的规律对于页岩气藏的有效开发非常重要.图7为有吸附气体和无吸附气体的日产量数值计算结果，从图中可以见到，当无吸附气体时，气体的产量下降非常快；初始阶段有无吸附气体对日产气量影响不大，但是到了后期，有吸附气体的日产量远远大于无吸附气体的日产量，吸附气体在生产后期能够维持日产气量.图8为有吸附气体和无吸附气体时压力的传播过程对比图，从图中可知，无吸附气体时气体的压力传播过程变快，吸附气体的存在可以延缓压力的传播，这是由于吸附气体的解吸能够补
充孔隙中气体的压力.
解吸常数pd和Vd也会影响气体的流动规律，图9为不同Vd条件下，日产气量
变化规律，Vd越大日产气量越大，且产量递减越缓慢.图10为不同pd条件下日
产气量变化规律，pd越大日产气量越大，pd对日产气量的影响程度不及Vd.
3.2 渗透率和扩散系数的影响
气体在页岩中的流动为渗流和扩散之和，渗透率和扩散系数势必会影响页岩的日产气量.图11为不同渗透率条件下的日产气量，渗透率越大日产量越高，图12为不同扩散系数条件下的日产气量，扩散系数越大，产气量越高.
3.3 孔隙度的影响
页岩中游离气的含量与孔隙度有关，游离气也会影响气体的流动规律，图13为不同孔隙度条件下的日产气量图，页岩孔隙度越高，游离气含量越多，产气量越大.
因此在页岩气藏开发中，孔隙度是非常重要的一项指标，直接关系到页岩的日产气量.
页岩基质解吸--扩散--渗流耦合实验结果显示气体的初期产量大，下降快，到了后期产量下降缓慢，这种现象与实际生产过程中页岩气井产量一致；压力在页岩中传播缓慢.页岩气的解吸--扩散--渗流耦合数学模型能够描述气体在页岩基质中的流动规律；吸附气体会影响气体的产量和压力传播速度，气体的解吸常数会影响气体的产量，Vd越大日产气量越大，pd越大日产气量越大；页岩基质的渗透率和扩散
系数也影响气体的流动规律，渗透率越大日产量越高，扩散系数越大，产气量越高；页岩基质孔隙度影响游离气含量，页岩孔隙度越高，游离气含量越多，产气量越大；在页岩气藏的开发过程中需要考虑气藏的解吸附能力、扩散系数、渗透率和孔隙度等物性参数，得到的数学模型为页岩气的产量预测提供了更准确的计算方法.
【相关文献】
1 张东晓,杨婷云.页岩气开发综述.石油学报,2013,34(4):792-801(Zhang Dongxiao,Yang Tingyun.An overview of shale-gas production.Acta Petrolei Sinica,2013,34(4):792-801(in Chinese))
2 李玉喜,张金川．我国非常规油气资源类型和潜力.国际石油经济，2011,3:61-67(Li Yuxi,Zhang Jinchuan.Unconventional gas &oil resource type and potential of China.International Petroleum Economy,2011,3:61-67(in Chinese))
3 Carlson NE,Mercer J.Devonian shale gas production:mechanisms and simple
models.Journal of Petroleum Technology,1991,43(4): 476-482
4 Warren JE,Root PJ.The behavior of naturally fractured reservoirs.Old SPE
Journal,1963,3(3):245-255
5 Kuuskpaa VA,Wiicks DE,Thurber JL.Geologic and reservoir mechanisms controlling gas recovery from the Antrim Shale.SPE Annual Technical Conference and Exhibition,1992
6 Carlso ES.Characterization of devonian shale gas reservoirs using coordinated single
well analytical models.SPE Eastern Regional Meeting,1994
7 Schepers K,Gonzalez R,Koperna G,et al.Reservoir modeling in support of shale gas tin American and Caribbean Petroleum Engineering Conference,2009
8 段永刚,魏明强,李建秋等.页岩气藏渗流机理及压裂井产能评价.重庆大学学报,2011,34(4):62-
66(Duan Yonggang,Wei Minqiang,Li Jianqiu,et al.Shale gas seepage mechanism and fractured wells’production evaluation.Journal of ChongQing University, 2011,34(4):62-
66(in Chinese))
9 李建秋,曹建红,段永刚等.页岩气井渗流机理及产能递减分析.天然气勘探与开发,2011,34(2):34-
37(Li Jianqiu,Cao Jianhong, Duan Yonggang,et al.Seepage mechanism and productivity decline of shale-gas reservoir.Natural Gas Eχporation&development, 2011,34(2):34-37(in Chinese))
10 李清泉,王新海,尹虎等.页岩气藏数值模拟及井底压力动态分析.东北石油大学学
报,2013,37(1):91-96(LiQinquan,WangXinhai,Yin Hu,et al.Simulation of shale-gas reservoirs and dynamic analysis of bottom hole pressure.Journal of Northeast Petroleum University,2013,37(1):91-96(in Chinese))
11 任建华,蔺景德,张亮等.页岩气藏吸附特征及其对产能的影响.新疆石油地质,2013,34(4):441-
444(Ren Jianhua,Lin Jingde, Zhang Liang,et al.Adsorption characteristics and impact on productivity of shale gas reservoir.Xinjiang Petroleum Geology,2013, 34(4):441-444(in Chinese))
12 李亚洲,李勇明,罗攀等.页岩气渗流机理与产能研究.断块油气田,2013,20(2):186-190(Li Yazhou，Li Yongming，Luo Pan，et al.Study on seepage mechanism and productivity of shale gas.Fault-Block Oil&Gas Field，2013，20(2)：186-190(in Chinese))
13 邓佳,朱维耀,刘锦霞等.考虑应力敏感性的页岩气产能预测模型.天然气地球科学,2013,24(3):456-460(Deng Jia,Zhu Weiyao, Liu Jinxia,et al.Productivity prediction model of shale gas
considering stress sensitivity.Natural Gas Geoscience,2013,24(3):456-460(in Chinese))
14 任飞，王新海，谢玉银等.考虑滑脱效应的页岩气井底压力特征.石油天然气学报,2013,35(3)：124-126(Ren Fei,Wang Xinhai,Xie Yuyin,et al.Characteristics of bottom pressure of shale-gas wellconsideringslippagee ff ect.JournalofOilandGasTechnology, 2013,35(3):124-126(in Chinese))
15 郭为，熊伟，高树生等.温度对页岩等温吸附/解吸特征影响研究.石油勘探与开发，2013,40(4)：481-485(GuoWei,XiongWei,Gao Shusheng,et al.Impact of temperature on the isothermal adsorption/desorption charac teristics of shale gas.Petroleum Eχploration and Development,2013,40(4)：481-485(in Chinese))。