四川盆地长宁地区海相页岩吸附气含量演化特征
——以N201井五峰组—龙马溪组一段为例
2021-04-06邹源红肖红纱陈义才彭文秋
聂 舟,衡 德,邹源红,肖红纱,陈义才,康 宇,彭文秋
1 四川长宁天然气开发有限责任公司勘探开发部;2 中国石油西南油气田公司气田开发管理部;3 成都理工大学能源学院
0 前 言
上奥陶统五峰组—下志留统龙马溪组一段(简称龙一段)是四川盆地长宁地区页岩气的主力产层段,其吸附气含量占页岩总含气量的20%~85%,为页岩气长期稳产提供了重要的物质基础。在喜马拉雅期构造抬升过程中,长宁地区龙马溪组上覆地层剥蚀厚度一般达到15 00~3 500 m,页岩吸附气含量随着地层温度与压力的变化,经历了复杂的演化过程。近年来,前人对川南地区五峰组—龙马溪组页岩吸附气含量及其影响因素等方面进行了比较系统的研究,但是关于页岩吸附气含量演化过程的研究相对薄弱[1-3]。为了进一步明确长宁地区五峰组—龙马溪组页岩气的成藏演化规律,本文根据四川长宁天然气开发有限责任公司15口评价井的现场解吸、测井解释以及13个页岩样品的等温吸附实验资料,在采用盆地数值模拟软件恢复典型评价井埋藏史和热演化史的基础上,结合页岩吸附气含量主要影响因素的相关性分析,建立了吸附气含量的多元回归数学模型,定量分析了五峰组—龙一段页岩自侏罗纪大量生气以来,页岩吸附气含量随埋藏深度、热演化程度变化的演化特征。
1 地质概况
长宁地区位于四川盆地南部(图1a),四川长宁天然气开发有限责任公司在该区的页岩气矿权登记区块面积约为7 000 km2。长宁地区地表大部分出露二叠系、三叠系及侏罗系,在长宁背斜核部,寒武系及奥陶系出露于地表(图1b);五峰组—龙一段埋藏深度一般在1 500~4 500 m 之间。长宁背斜长轴总体呈NWW—SEE 向,东北翼地层倾角为50°~65°,西南翼地层倾角为30°~45°。长宁背斜在早燕山期初具雏形,喜马拉雅期构造挤压使其地层变形更加剧烈,同时伴随持续隆升,奥陶系及其上覆地层遭受不同程度的风化剥蚀[4-6]。
目前,矿权区内完钻评价井有20余口。页岩气主力勘探开发层段为五峰组—龙一段下亚段(图2),页岩厚度一般为35~50 m,岩性为碳质页岩和黑色页岩,页岩总有机碳含量(TOC)一般在2%~4%之间。上覆龙一段上亚段厚度为120~150 m,岩性为灰色页岩及粉砂质泥岩,TOC一般为0.5%~2%。
图2 长宁地区五峰组—龙一段综合柱状图Fig.2 Comprehensive column of Wufeng Formation-Longmaxi Member 1 in Changning area
2 页岩现今吸附气含量
页岩吸附气是指吸附在页岩干酪根和无机矿物颗粒表面的天然气。吸附气含量通常以每吨页岩中所含吸附气在标准状态下(0℃,101.325 kPa)的体积来代表。目前,页岩吸附气含量难以通过密封高压取心直接测量,主要通过等温吸附、现场解吸、测井解释等方法获取。
2.1 等温吸附实验含气量
页岩等温吸附是页岩气解吸的逆过程。根据页岩吸附机理和等温吸附实验,可以定量分析页岩的吸附能力,明确吸附气含量和压力的关系。关于页岩气吸附机理的研究,许多学者基于页岩样品特征建立了不同的页岩气吸附模型,如Langmuir模型、BET 模型、SDR 模型等[7-12]。国外的勘探实践表明,北美页岩气吸附特征大多服从Langmuir等温吸附。Langmuir 吸附模型的基本假设条件是表面均匀光洁的页岩固体形成单分子层吸附,并且被吸附的气体分子之间没有相互作用力,吸附平衡时处于一种动态平衡。根据Langmuir 等温吸附曲线,可计算页岩在某一温度不同地层压力条件下的吸附气含量:
公式(1)中:Vg为地层条件下的吸附气含量,m3/t;P为地层孔隙气体压力,MPa;VL为Langmuir 体积(吸附量达到饱和时的吸附气含量),m3/t;PL为Langmuir压力(VL的一半所对应的压力),MPa。
长宁地区五峰组—龙一段页岩样品在地层温度下的等温吸附特征曲线(图3)表明:一般在0~20MPa 的压力增长区间内,页岩吸附气含量随之快速增大;当压力超过20MPa 继续增加时,页岩吸附气含量增长速率放缓直至达到饱和。
长宁地区五峰组—龙一段5口井13个页岩样品的Langmuir 压力(PL)在4.28~38.50 MPa 之间(表1),其中10 个样品的Langmuir 压力低于15MPa,表明甲烷气体容易被页岩吸附;Langmuir 体积(VL)有3个样品大于4 m3/t,7 个样品在2~4 m3/t 之间,3 个样品为1.2~2 m3/t,总体上反映页岩具有较强的吸附能力。根据公式(1)计算的结果,五峰组—龙一段页岩在现今地层条件下的吸附气含量(Vg)在0.58~3.08 m3/t之间,平均值为1.80 m3/t。
图3 长宁地区五峰组—龙一段页岩甲烷等温吸附曲线Fig.3 Methane isotherm adsorption curve of the shale of Wufeng Formation-Longmaxi Member 1 in Changning area
表1 长宁地区五峰组—龙一段页岩等温吸附实验数据与回归计算结果Table 1 Isothermal adsorption experimental data and regression calculation results of the shale of Wufeng Formation-Longmaxi Member 1in Changning area
2.2 现场解吸含气量
页岩现场解吸的总含气量包括解吸气、残留气和损失气等3部分。解吸气量采用排水集气的原理测定,即将岩心样品装入解吸装置后,利用水浴加热至地层温度,在大气压力下收集释放的吸附气。然后在密闭的球磨罐中将样品破碎计量滞留在岩心中的残余气量[13-14]。现场解吸时,由于钻井取心过程中地层温度、压力的变化,岩心样品在装入解吸罐之前必然有部分吸附气损失。损失气量的大小,主要受取心方式、测定方法、逸散时间等影响,可通过岩心颗粒的球形扩散原理估算,或者选取岩心装罐开始后10~15个解吸点数据进行曲线拟合推算损失气量[15]。
通过对长宁地区15 口评价井共344 个页岩样品现场解吸测定结果的统计,五峰组—龙一段页岩解吸总含气量在0.48~4.19 m3/t 之间,平均值为1.78 m3/t,主要分布区间为1.0~2.5 m3/t,总含气量超过2.5 m3/t的样品约占13%(图4)。
图4 长宁地区五峰组—龙一段页岩解吸气含量统计直方图Fig.4 Histogram of desorbed gas content of the shale of Wufeng Forma⁃tion-Longmaxi Member 1 in Changning area
2.3 测井解释含气量
测井解释页岩吸附气含量是以Langmuir 等温吸附为理论依据,利用等温吸附实验数据和大量测井数据,优选拟合回归地层温度、压力、总有机碳含量等关键参数,建立基于常规测井资料的计算模型而获得的吸附气含量[16-18]。测井解释法不仅能够估算吸附气含量,而且通过页岩孔隙度、含气饱和度及地层压力等参数可以同时计算页岩游离气含量。近年来,以等温吸附和体积模型为基础的页岩含气量测井解释模型不断发展和完善,在国内外页岩气开采中得到广泛应用[19-20]。
长宁地区15 口评价井五峰组—龙一段测井解释页岩吸附气含量在1.04~2.33 m3/t 之间,平均值为1.72 m3/t。页岩吸附气含量在不同评价井之间的差异较小(图5),但是吸附气含量占总含气量的比例变化较大,在19.6%~85.2%之间,平均为43%;注意到当页岩吸附气含量在总含气量中所占比例超过50%时,总含气量一般低于4 m3/t。
图5 长宁地区五峰组—龙一段测井解释页岩含气量统计直方图Fig.5 Histogram of shale gas content by logging interpretation of Wufeng Formation-Longmaxi Member 1 in Changning area
对比等温吸附法、现场解吸法、测井解释法的结果可见,五峰组—龙一段页岩吸附气含量主要分布在1.0~2.5 m3/t 之间,3 种方法获得的吸附气含量虽然存在一定差异,但是总体上比较相近,平均值在1.72~1.80 m3/t之间。
3 页岩吸附气含量影响因素及多元回归
页岩吸附气含量的影响因素较多,其中内在因素包括有机碳含量、矿物组成、有机质成熟度,外部因素包括地层压力与地层温度[1]。
3.1 页岩总有机碳含量
有机碳是页岩生气的物质基础。页岩中有机质存在大量的纳米孔隙,可产生巨大的比表面积,为吸附气提供了吸附载体。对长宁地区5口井13个等温吸附实验结果的Langmuir甲烷气体吸附体积VL(单位重量页岩饱和吸附气量)、地层条件下的吸附气含量Vg分别与TOC作相关性分析(图6),五峰组—龙一段吸附气含量与TOC均呈良好的正相关关系,表明有机碳含量的增高能够显著增强页岩吸附气体的能力。
图6 长宁地区五峰组—龙一段页岩吸附气含量与TOC相关性Fig.6 Correlation between adsorbed gas content and TOC of the shale of Wufeng Formation-Longmaxi Member 1 in Changning area
3.2 页岩无机矿物组成
长宁地区五峰组—龙一段页岩无机矿物主要由石英、黏土矿物组成。黏土矿物属于层状硅酸盐矿物,比表面较柱状石英的大,在相同条件下吸附气体的能力高于石英。通过长宁地区5口井13个等温吸附实验结果的Langmuir甲烷气体吸附体积VL与页岩黏土矿物、石英含量的相关性分析,可见吸附气含量与黏土矿物含量之间存在一定的正相关性(图7a),与石英含量之间存在微弱的负相关性(图7b)。
图7 长宁地区五峰组—龙一段页岩吸附气含量与矿物组成相关性Fig.7 Correlation between adsorbed gas content and mineral composition of the shale of Wufeng Formation-Longmaxi Member 1 in Changning area
长宁地区五峰组—龙一段页岩吸附气含量与石英含量之间的相关性,可能受到石英来源的影响。薄片观察发现,五峰组—龙一段页岩的石英存在两种主要来源:一是盆外的陆源输入,二是盆内生物成因和化学成因。陆源输入的石英含量增加往往降低页岩的有机质含量。在深水还原环境发育的富有机质页岩中,当盆内来源的石英较发育时,虽然在一定程度上会降低页岩岩石的吸附能力,但是页岩中富含的有机质可使吸附能力整体增加。
3.3 地层温度与压力
气体在固体表面的吸附作用是一个放热的过程,即温度升高会抑制气体的吸附。温度升高会使气体分子的热运动更加活跃,促使其脱离吸附介质的束缚,解吸而转化为自由游离态。当压力一定时,页岩岩样原地吸附能力与温度呈负相关关系。压力的增大能够有效减小气体分子吸附所需的结合能,进而让气体分子更容易被固体介质吸附[21]。
长宁地区五峰组—龙一段页岩在不同地质时期埋藏深度变化幅度大,尤其是晚白垩世以来,随着喜马拉雅期构造抬升和上覆地层的大幅度剥蚀,地层温度和压力显著降低。由于地层温度和地层压力对页岩吸附性能具有相互抑制的作用,因此页岩吸附气含量在不同地质时期的演化主要取决于地层温度和地层压力的综合影响。
3.4 页岩吸附气含量多元回归
为了定量反映不同地质时期吸附气含量随地层温度与压力同时发生变化的演变过程,根据长宁地区五峰组—龙一段页岩等温吸附实验结果,通过多元回归分析建立了吸附气含量与TOC、Ro、地层温度及地层压力的回归公式:
公式(2)中:V为吸附气含量,m3/t;T为地层温度,℃;P 为地层压力,MPa;Ro为镜质组反射率,%;TOC 为有机碳含量,%。
把通过公式(2)计算的吸附气含量与长宁地区五峰组—龙一段不同样品的等温吸附实测结果(Vg)进行对比分析(表1),发现存在两种情况:①当页岩TOC≤2%和埋藏深度≤1 300 m 时,通过公式(2)计算的结果与实测结果的误差较大,原因一是TOC 低于2%的样品点较少,二是当埋藏深度较浅时,页岩吸附气保存条件变差,页岩气散失作用强烈;②当页岩TOC 大于2%和埋藏深度超过1 300 m 时,通过公式(2)计算的结果与等温吸附实测结果的误差值一般低于0.25 m3/t,其中N201 井五峰组和龙一段页岩计算的吸附气含量分别为2.81 m3/t和2.84 m3/t,与等温吸附实验结果的误差值分别为-0.12 m3/t和0.24 m3/t。总体而言,五峰组—龙一段优质页岩在距离长宁背斜核部较远、埋藏深度较大、吸附气保存条件较好的地区,通过回归公式计算的吸附气含量具有较高的可靠性。
4 典型井页岩热演化与吸附气含量演化
N201井在五峰组—龙一段获得高产页岩气,页岩TOC 平均值为3.54%,吸附气含量为1.67 m3/t,总含气量为6.47 m3/t。该井位于长宁背斜构造的南翼(图1b),地表出露下三叠统嘉陵江组,龙一段顶界埋深为2 495.2 m,五峰组—龙一段页岩热演化史和吸附气含量演化具有一定代表性。
4.1 页岩热演化史
页岩热演化史主要取决于页岩有机质经历的沉积埋藏史和古地温。长宁地区志留系龙马溪组沉积以来,在加里东、东吴、印支及喜马拉雅构造运动期间发生过多次地层剥蚀,其中喜马拉雅期地层剥蚀强烈[22-24]。根据钻井地层厚度趋势法估算,长宁背斜核部在喜马拉雅期的地层剥蚀厚度达4 500~5 500 m,在背斜翼部地层剥蚀厚度一般为1 500~3 000 m(图1b)。
根据烃源岩Ro、不同期次包裹体均一温度对古热流进行恢复:四川盆地在加里东期之前的热状态较为稳定,热流值较低,海西期热流开始逐渐增大,印支期及其以后古热流逐渐降低[25-26]。采用盆地数值模拟软件恢复长宁地区N201 井的古埋深和古地温(图8):五峰组—龙一段在二叠纪时的埋藏深度为1 500~2 500 m,地温梯度为(4.2~5.4)℃/100 m,地层温度为95~130 ℃;三叠纪—白垩纪,埋藏深度和地层温度持续增加,但是地温梯度逐渐降低到(2.8~3.3)℃/100 m;白垩纪末期,埋藏深度和地层温度分别达到5 300 m 和193 ℃;古近纪以来,随着燕山期—喜马拉雅期构造抬升作用与上覆地层的剥蚀,埋藏深度和地层温度分别逐渐降低到现今的2 495 m和93℃(图8)。
4.2 页岩吸附气含量演化
长宁地区五峰组—龙一段页岩有机质类型为腐泥型,在成熟阶段有机质以生油为主,伴生的天然气主要溶解在液态烃中而呈油相,在高—过成熟阶段液态烃、干酪根裂解而形成气相产物。根据N201 井五峰组—龙一段页岩有机质热演化史恢复结果(图8),目的层在二叠纪—三叠纪处于成熟阶段,侏罗纪及其以后进入高—过成熟阶段。以白垩纪末期最大埋深为界,页岩吸附气含量演化可划分为富集和解吸两个阶段。
页岩吸附气富集阶段发生在侏罗纪—白垩纪。侏罗纪—白垩纪,五峰组—龙一段页岩埋藏深度持续增加,页岩有机质经过生油高峰开始大量生气,地层温度的不断升高虽然不利于烃类气体的吸附,但是页岩孔隙流体压力的增加和有机质成熟度的升高有利于气体被吸附。N201 井五峰组—龙一段页岩现今地层压力为49.87 MPa,按照理想气体估算,在白垩纪末期最大埋深时的地层压力为68.94 MPa。通过公式(2)计算,页岩吸附气含量随着埋藏深度增加逐渐升高,白垩纪末期达到4.42 m3/t。
页岩吸附气的解吸阶段发生在古近纪—第四纪。古近纪以来,随着上覆地层被大幅度剥蚀,五峰组—龙一段页岩的地层温度和地层压力逐渐降低到现今状态,页岩Ro基本保持不变。地层温度的降低虽然会导致吸附气含量的增加,但是在低地温梯度背景下地层压力的降低更有利于页岩吸附气的解吸。通过公式(2)计算,页岩吸附气含量随着埋藏深度降低而逐渐降低到1.67 m3/t。
图8 N201井埋藏史、热史及五峰组—龙一段页岩吸附气含量演化史Fig.8 Burial history,thermal history and evolution history of shale adsorbed gas content of Wufeng Formation-Longmaxi Member 1 of Well N201
5 结 论
(1)根据等温吸附法、现场解吸法、测井解释法的综合分析结果,长宁地区五峰组—龙一段页岩吸附气含量主要分布在1.0~2.5 m3/t 之间,但是吸附气含量占总含气量的比例变化较大。
(2)侏罗纪以来,在低地温梯度背景下,五峰组—龙一段页岩吸附气含量经历了随埋藏深度增加而升高、在抬升过程中随埋深变浅而降低的演化过程。喜马拉雅期构造抬升幅度较大,不利于页岩气的保存,五峰组—龙一段页岩吸附气含量降低。
(3)随着将来页岩等温吸附实验的温度、压力上限的提高,以及国内外深层页岩气生产数据的补充,对吸附气含量的多元回归可进一步改善。