龙马溪—五峰组富有机质页岩三维岩石物理模板分析及“甜点”预测
2021-08-03檀文慧巴晶符力耘JosCarcione周欣
檀文慧， 巴晶 ， 符力耘， José M. Carcione， 周欣
1 河海大学地球科学与工程学院， 地球探测研究所， 南京 211100 2 中国石化石油物探技术研究院， 南京 211103 3 中国石油大学(华东)， 地球科学与技术学院， 山东青岛 266580 4 National Institute of Oceanography and Applied Geophysics (OGS), Trieste, Italy
0 引言
页岩气主要是指以吸附态和游离态赋存于富含有机质的低孔、低渗页岩储集系统中的一种天然气，已成为全球非常规油气勘探开发的重点领域(邹才能等，2014，2016).从我国油气资源分布来看，四川盆地及周缘的龙马溪—五峰组页岩成为国内页岩气勘探的首选层位之一，其储层富含有机质，埋藏深度适中、有机质演化程度高(邓继新等，2018).有机碳含量可以反映烃源岩的生烃能力，寻找优质页岩气储层不仅需要富含有机质的区域，同时也要寻找研究区内储层脆性高、裂缝发育、有利于压裂改造的层位(张晨晨等，2016).页岩储层经过压裂改造，可显著提高页岩气的产量，因此脆性评价与储层压裂品质息息相关(李庆辉等，2012；Smart et al., 2014；Lisle, 1997；Yang et al., 2015).一般认为富含有机质的脆性页岩，且孔隙度和微裂缝相对比较发育的储层，是页岩气“甜点区”勘探开发的首选目标，其关键评价指标及下限为压力系数大于1.3、TOC值大于3.0%、孔隙度大于3.0%及含气量大于3.0 m3/t(邹才能等，2014，2016).
岩石物理模型是连接测井数据、地震属性与储层参数的纽带，通过等效简化的方法将宏观声学响应和微观结构特征等实际情况联系起来，从而有效地预测弹性参数，为储层地震属性反演提供理论基础(Ba等，2017；Zhang等，2020，2021).Xu和White(1995)基于Kuster-Toksöz理论和微分等效介质模型(Differential effective medium，DEM)，将岩石孔隙分为砂岩孔隙和泥岩孔隙，并提出用纵横比来描述孔隙的几何形态，最终形成适用于砂泥岩储层的Xu-White模型.Guo等(2013)根据自洽模型(Self-consistent approximation，SCA)和Backus平均方法，提出了富含有机质的岩石物理模型，针对北美页岩分析了矿物成分、孔隙度等参数对弹性性质与脆性的影响.董宁等(2014)在SCA和DEM的基础上，引入Berryman三维孔隙形态及Brown-Korringa固体替换技术，建立适用于富有机质泥页岩的岩石物理模型.邓继新等(2015)根据岩石的矿物、孔隙度、TOC含量及微观结构特征对龙马溪组页岩进行地震岩石物理建模，研究结果可为该地区的测井及地震甜点预测提供理论依据.Yang等(2017)通过获取的四川盆地页岩样品的数据，总结出地震岩石物理特征与微观结构的联系，并加入了孔隙纵横比这一影响因素，合理运用SCA、DEM、Gassmann等多个模型，建立了适合页岩储层的岩石物理模型.印林杰等(2020)考虑微纳米孔隙和不同成熟度下干酪根的赋存状态，采用微纳米孔隙理论描述页岩微纳米孔隙特性，结合各向异性SCA-DEM模型、各向异性Eshelby-Cheng模型和Brown-Korringa固体替换方程等建立一种页岩储层岩石物理模型.
岩石物理模板(RPT)是从地质条件上进行约束的岩石物理模型，是预测岩性和流体的有利工具(Avseth et al., 2005)，可以直接反映储层不同参数(孔隙度、矿物等)和弹性参数之间的关系，能够定量地解释测井和地震数据(Perez et al., 2011；Guo et al., 2013；巴晶等, 2013；Pang et al., 2019).李宏兵等(2019)基于等效介质理论和Gassmann方程建立多重孔隙岩石的三维岩石物理模板，结合测井和地震数据校正模板，从而有效预测复杂储层孔隙扁度、孔隙度以及孔隙中的流体饱和度，提高了复杂储层孔隙参数的预测精度.Tan等(2020)考虑致密粉砂岩的矿物组分、孔隙度和微观结构特征，结合DEM和SCA模型建立致密粉砂岩岩石物理模型，构建多参数影响的弹性参数二维岩石物理模板，根据实际资料对致密储层的脆性矿物进行定量预测.庞孟强等(2020)考虑致密砂岩储层粒间孔隙和微裂隙的双重孔隙结构，采用Biot-Rayleigh双重孔隙方程，构建二维致密砂岩的岩石物理模板，结合实际资料对储层孔隙度和裂隙含量进行了定量预测，从而可以有效预测优质储层的分布区域.
页岩地震弹性性质不仅受矿物组分、有机碳含量(TOC)、孔隙度的影响，而且还与本身的微观结构有关.总体而言，由于其复杂性，目前页岩气“甜点”预测难度较大，难以精细刻画页岩气商业开采“甜点区”(赵文智等，2016；陈胜等，2017).本文根据丁山区龙马溪—五峰组的页岩样品的测井解释结果，分析其基本的物性和岩石物理特征，结合EIAS 理论模型构建对有机质和脆性敏感的三维岩石物理模板，将岩石裂隙、矿物组分、孔隙度和地震响应特征联系起来.利用页岩的测井数据对模板进行校正，基于校正后的三维岩石物理模板对工区地震数据进行应用，进而对页岩储层“甜点区”敏感参数进行定量预测.
1 工区基本特征
研究对象为丁山区志留系下统龙马溪组至奥陶系上统五峰组的页岩气储层，优质储层主要集中在龙马溪组底部-五峰组，厚度约为30 m.岩石类型以硅质页岩为主，其次为混合质页岩及少量粘土质页岩.丁山构造龙马溪组下段-五峰组地层为深水陆棚相沉积，笔石类生物化石丰富，机械压实过程中同时伴有生物成因的硅质胶结，造成岩石样品表现出高TOC 含量(邓继新等，2018).
1.1 岩芯样品特征
本次岩芯取自A井，页岩气储层孔隙度范围为0～7%，深度为3582.44～3729.41 m，采用X射线衍射实验定量分析岩芯的矿物组分，各个矿物组分含量随深度的变化关系见图1.随着深度的增加，岩石的石英矿物随之增大，含量范围为22.2%～69.1%，平均约为36%，黏土含量随着深度增大而减小，范围在13.2%～59.4%，平均成分约为39%，长石含量(斜长石和钾长石)平均约为8%，碳酸盐岩矿物含量(方解石和白云石)平均约为13%，黄铁矿含量平均约为2%，有机碳(TOC)含量范围为0.46%～5.86%，平均值为1.5%.图2是矿物组分分布的三元图，从图中可以看出矿物组分主要是由黏土和硅质矿物组成，成岩的过程中伴随着硅质胶结.页岩的孔隙空间主要由粒间孔隙、溶蚀裂缝以及有机质孔组成(邓继新等，2015).
图1 丁山区龙马溪组至五峰组页岩的矿物组分随着深度的变化规律Fig.1 Mineral components of the shale from Longmaxi to Wufeng Formations in Dingshan Area vary with the depth
图2 丁山区龙马溪至五峰组页岩矿物组分三元图Fig.2 Ternary of mineral components of the shale from Longmaxi to Wufeng Formations in Dingshan area
1.2 有机碳含量TOC
1.2.1 TOC与矿物组分之间的关系
图3是有机碳含量(TOC)与不同矿物组分的关系.图3a是TOC与石英含量之间的关系，随着TOC的增大，岩石石英矿物含量相应增大，具有较好的正相关性.前人研究认为石英属于脆性矿物(Jarvie et al., 2007)，因此可认为富有机质的层位脆性矿物含量高.图3b是TOC与碳酸盐岩矿物含量的变化关系，当TOC含量在0～1%之间增大时，碳酸盐岩矿物含量快速下降，随着TOC含量继续增大，碳酸盐岩矿物含量稳定在10%左右.图3c是TOC随黏土矿物含量的变化关系，随着TOC含量的增大，岩石的黏土矿物含量随之减少.同样一些学者研究发现硅质含量越高，页岩脆性越高，可压裂性好，越有利于裂隙的产生，裂隙发育对页岩气藏的开采效益和产量有直接的影响(Montgomery et al., 2005；Bowker, 2007；李登华等，2009；王淑芳等，2014).生物成因硅质重结晶可以形成石英(秦建中等，2010)，从而增加页岩的脆性.有机碳含量(TOC)对页岩气的富集是很重要的控制因素，决定了页岩生烃的能力，脆性矿物影响了岩石的压裂改造，岩石越脆，越容易压裂.富有机质、脆性矿物含量高的储层是页岩气勘探开发的主力层系.
图3 有机碳含量(TOC)与矿物组分的关系(a) TOC与石英的关系； (b) TOC 与碳酸盐岩的关系； (c) TOC与黏土的关系.Fig.3 The relation between TOC and mineral components(a) The relation between TOC and quartz content; (b) The relation between TOC and carbonate content; (c) The relation between TOC and clay content.
1.2.2 TOC与孔隙度、弹性参数的关系
关于TOC含量的预测，一般需要建立TOC含量与其敏感性参数之间的定量关系，然后利用此定量关系进行反演，从而获得TOC含量的平面分布(陈胜等，2017).如图4所示，TOC含量与孔隙度的关系中，TOC含量与孔隙度具有较好的正相关性，相关系数R2=0.83.该研究区二者之间的定量关系式为
TOC=0.85-0.198φ+0.144φ2，R2=0.83.
(1)
因此可以利用公式(1)，将得到的孔隙度转换成TOC含量.
图4 TOC和孔隙度的关系Fig.4 The relation between TOC and porosity
为了优选对有机碳含量敏感的储层参数，根据测井解释数据(井A)，如图5所示，图中色标是TOC含量，本文分析了目的层的有机碳含量与不同弹性参数和物性参数之间的关系.图5a是纵波阻抗与拉梅常数和密度的乘积的关系交会图，从图中可以看出λρ随纵波阻抗的增大而增大，λρ低(18～30 GPa·g·cm-3)和纵波阻抗低(9～10 g·cm-3·km·s-1)的情况下，对应TOC含量高；图5b显示的是孔隙度和密度的关系交会图，二者具有较好的负相关性，且孔隙度大、密度小则对应高的TOC含量；图5c是剪切模量与纵波阻抗的关系交会图，可以看出剪切模量随着纵波阻抗的增大而增大，高TOC含量对应着较高的剪切模量(13～18 GPa)；同样，图5d显示泊松比在0.18～0.22范围内，页岩气储层的TOC 含量高.
图5 (a) 纵波阻抗与拉梅常数和密度的乘积的交会图； (b) 孔隙度与密度的交会图； (c) 纵波阻抗与剪切模量的交会图； (d) 纵波阻抗与泊松比的交会图.色标是TOCFig.5 (a) Cross-plot between P-wave impedance and λ ρ; (b) Cross-plot between porosity and density; (c) Cross-plot between P-wave impedance and shear modulus; (d) Cross-plot between P-wave impedance and Poisson ratio. The color bar denotes TOC.
1.3 脆性指数
页岩的脆性评价方法没有统一的标准，目前我国主要是借鉴北美页岩气勘探开发中总结的经验公式，依据弹性参数定义的脆性指数来评价有利压裂层段(Rickman et al., 2008; Tan et al., 2018).杨氏模量是表征岩石在受力时维持裂缝扩展的能力，而泊松比则是反映应力的变化、岩石发生破坏时的强弱(李庆辉等，2012).Rickman等(2008)通过对北美页岩储层进行大量数据统计分析后提出了公式(2)—(4)，认为杨氏模量大和泊松比小的储层脆性高.
(2)
(3)
(4)
其中E是杨氏模量，υ是泊松比，Emax(Emin)是最大(最小)杨氏模量，υmax(υmin)是最大(最小)泊松比，最大值和最小值为测量数据中的极值.
同样为了分析对页岩脆性敏感的岩石物理参数，基于不同的弹性参数对页岩脆性进行敏感性分析.图6分别是剪切模量、泊松比与纵波阻抗的关系交会图，色标是脆性指数，结果显示高脆性对应高剪切模量，泊松比小对应岩石脆性高，剪切模量和泊松比能够很好地区分高脆性和低脆性储层.
图6 (a) 剪切模量与纵波阻抗的交会图，色标是TOC； (b) 泊松比和纵波阻抗的交会图，色标是脆性指数Fig.6 (a) Cross-plot between P-wave impedance and shear modulus, where the color bar denotes TOC; (b) Cross-plot between P-wave impedance and Poisson ratio, where the color bar denotes brittleness index
如图7所示，在TOC 含量与脆性指数的交会图中可知，当TOC>3%，BI>50%，储层属于优质页岩气储层，当2%
图7 TOC与脆性指数的交会图，色标是泊松比Fig.7 Cross-plot between TOC and brittleness index (The color bar is Poisson′s ratio)
2 页岩气“甜点区”敏感参数地震预测方法
页岩气储层“甜点区”是指具有勘探开发价值的页岩气富集区，主要表现为页岩储层厚度大，总有机碳(TOC)含量高，处于“生气窗”，含气量高，裂隙发育，脆性高易压裂等(袁书坤等，2014).本文首先在测井和地震资料的基础上，明确页岩储层的岩石物理特征，对“甜点区”敏感性参数进行优选，进行地震岩石物理建模.最后对优质储层的TOC、脆性以及孔隙度等核心参数进行定量预测，综合评价“甜点区”.
2.1 理论方法
对实际的含气储层进行岩石物理建模一般包括三个步骤(Xu and White，1996；Xu and Payne，2009)：(1)分析矿物组分的配比，计算岩石基质的弹性参数和密度；(2)计算岩石干骨架弹性模量与密度；(3)流体替换.
2.1.1 Voigt-Reuss-Hill理论
本文计算岩石基质的弹性模量采用Voigt-Reuss-Hill 平均法(Hill，1952)，即可得到理论的估测值，具体公式如下：
(5)
其中，i表示矿物的相或者孔隙空间的第i个介质，fi表示第i个介质的体积含量，Mi表示的是第i个介质的弹性模量.
2.1.2 等效嵌入体应力平均(EIAS) 理论模型
假设多孔岩石中含有球形孔隙和裂隙，饱和岩石的高频体积模量和剪切模量为(Endres and Knight, 1997)：
(6)
(7)
其中(Berryman, 1980a; Mavko et al., 2009)，
(8)
(9)
(10)
(11)
其中P1和Q1对应于球形孔，而P2和Q2近似对应硬币状裂隙.这些参数有更精确的表达式，适用于任意纵横比的扁状球体，包括球形孔隙和裂隙.
(12)
其中Tiijj和Tijij表示椭球状嵌入体内的应变与均匀远场应变场相关的张量(Berryman，1980b；Mavko等，2009).
当流体压力在整个岩石孔隙空间达到平衡状态时，低频饱和岩石的体积模量和剪切模量(Endres and Knight, 1997)为
(13)
(14)
此处
(15)
(16)
(17)
并且，
P0n=Pn(Kf=0), (n=1,2)，Q0n=Qn(Kf=0),
(n=1,2).
(18)
针对实际储层测井岩芯物性和岩石物理特征，具体的建模流程如下：
(1)首先针对实际储层测井岩芯的矿物组分，使用Voigt-Reuss-Hill边界理论估算页岩基质的体积模量和剪切模量，页岩不同矿物组分的弹性性质如表1所示.
表1 矿物组分性质 (Mavko et al., 2009)Table 1 Mineral component properties
(2)其次将岩石孔隙空间分为硬孔隙和裂隙，根据不同孔隙纵横比的大小表征岩石孔隙的大小，使用EIAS理论模型估算页岩干骨架的体积模量和剪切模量；
(3)如果岩石内部的孔隙填充流体，需要对岩石进行流体替换，可利用EIAS理论模型估算页岩饱和流体状态下的体积模量和剪切模量；
(4)最后根据构建的页岩岩石物理模型预测岩石的纵横波速度，结合岩石的体积密度，计算页岩气储层的弹性参数.
2.2 建模结果
2.2.1 不同影响因素下的岩石物理模板
为了识别高品质的页岩气储层，可基于不同影响因素下的岩石物理模板特征，分析页岩气储层的不同因素对“甜点区”敏感性弹性参数的影响.根据测井解释结果，可设页岩孔隙度的变化范围为0.5%～6.5%，裂隙纵横比变化范围为0.0006～0.0066，石英矿物体积分数从20%增加到70%，黏土矿物体积分数则相对从60%减少到10%，同时方解石矿物体积分数和长石矿物体积分数分别设为10%和7%，有机碳含量取平均值，体积分数约为3%，裂隙含量设为定值，占孔隙度的1.5%.
图8a是不同的孔隙度和矿物组分影响下剪切模量-密度-拉梅常数和密度的乘积的岩石物理模板，此时裂隙纵横比保持定值.结果显示随着孔隙度的增大，剪切模量、密度以及拉梅常数和密度的乘积都随之减小，随着石英含量的增大和黏土含量的相对减小，剪切模量和密度有增大的趋势，拉梅常数和密度的乘积变化不明显.图8b是不同的裂隙纵横比和矿物组分影响下，剪切模量-密度-拉梅常数和密度的乘积的岩石物理模板，此时孔隙度保持常数.随着裂隙纵横比的增大，剪切模量和拉梅常数和密度的乘积都随之增大，密度有略增大的趋势.图8c是不同的裂隙纵横比和孔隙度影响下，剪切模量-密度-拉梅常数和密度的乘积的岩石物理模板，此时矿物组分保持定值.随着孔隙度的增大，剪切模量、密度和拉梅常数和密度的乘积均减小，而随着裂隙纵横比的增大，三者均随之增大.
图8 不同影响因素下的弹性参数岩石物理模板(a) 在孔隙度和矿物组分的影响下岩石物理模板； (b) 在裂隙纵横比和矿物组分的影响下岩石物理模板； (c) 在孔隙度和裂隙纵横比的影响下岩石物理模板.Fig.8 The rock physics template by elastic parameters under different influencing factors(a) The effect of porosity and mineral component on the rock physics template; (b) The effect of crack aspect ratio and mineral component on the rock physics template; (c) The effect of porosity and crack aspect ratio on the rock physics template.
2.2.2 三维岩石物理模板
根据2.2.1页岩岩石物理模板设置的参数，结合不同影响因素下的弹性参数构建三维岩石物理模板，如图9所示.
图9a是在孔隙度、裂隙纵横比及矿物组分影响下的剪切模量-密度-泊松比的三维岩石物理模板.随着孔隙度的增大，剪切模量、密度及泊松比都具有减小的趋势.随着石英含量的增加，黏土含量的相对减小，剪切模量和密度随之增大，而泊松比减小.随着裂隙纵横比的增大，剪切模量、密度及泊松比都具有增大的趋势.图9b是在孔隙度、裂隙纵横比及矿物组分影响下的剪切模量-密度-拉梅常数和密度的乘积的三维岩石物理模板.随着孔隙度的增大，剪切模量、密度以及拉梅常数和密度的乘积都随之减小.随着石英含量的增大和黏土含量的相对减小，剪切模量和密度有增大的趋势，拉梅常数和密度的乘积变化不明显.随着裂隙纵横比的增大，剪切模量和拉梅常数和密度的乘积随之增大，密度变化不明显.
图9 不同影响因素下的三维岩石物理模板(a) 在孔隙度、裂隙纵横比及矿物组分影响下的剪切模量-密度-泊松比的岩石物理模板； (b) 在孔隙度、裂隙纵横比及矿物组分影响下的剪切模量-密度-拉梅常数和密度的乘积的岩石物理模板.Fig.9 3D rock physics template under different influencing factors(a) The effect of porosity, crack aspect ratio and mineral component on μ-ρ-υ rock physics template; (b)The effect of porosity, crack aspect ratio and mineral component on μ-ρ-λ ρ rock physical template.
3 测井数据校正三维岩石物理模板
本文针对页岩气储层建立了基于有机碳含量和脆性敏感的岩石物理模板，在利用图板对储层“甜点区”进行定量预测之前，需要利用测井数据对模板进行校正，从而获得符合实际储层的岩石物理模板，进而在图板的储层参数(文中指孔隙度、石英含量及裂隙纵横比)范围内，对储层相关参数进行定量的预测.
如图10给出校正后三维岩石物理模板与测井数据的对比.图10a和10b是剪切模量-密度-拉梅常数和密度的乘积岩石物理模板与测井数据的对比情况，色标是测井数据的孔隙度和有机质含量.对比测井数据和模板，结果显示二者的一致性较好，测井数据的孔隙度分布与模板的变化趋势一致，且高孔隙度的数据点对应高的有机质含量.同时高的有机碳含量具有较高的剪切模量和较低的拉梅常数和密度的乘积，以及低的密度，较高的石英含量.图10c是剪切模量-密度-泊松比岩石物理模板与测井数据的对比情况，色标是脆性指数.图中可知脆性较高的数据点泊松比较低，且对应高的有机质含量和孔隙度.因此基于该模板可对储层孔隙度、石英含量以及裂隙纵横比进行预测.
图10 三维岩石物理模板与测井数据(a) 色标为孔隙度； (b) 色标为有机碳含量； (c) 色标为脆性指数.Fig.10 3D rock physics template and log data(a) The color bar is porosity; (b) The color bar is TOC; (c) The color bar is brittleness index.
4 “甜点区”储层参数地震定量预测
本文对工区一条过A、B和C井的测线进行目的层孔隙度、裂隙纵横比及石英含量的定量预测.过A井目的层深度为3702～3731 m，厚度为28.9 m，过B井目的层深度为4332～4367 m，厚度为35.5 m，过C井目的层深度为2239～2272 m，厚度为32 m.首先提取反演得到的纵波阻抗、横波阻抗和密度，得到该测线的纵波阻抗、横波阻抗和密度二维剖面(图11)，根据纵波阻抗、横波阻抗和密度计算剪切模量、泊松比、拉梅常数和密度的乘积，得到该测线的剪切模量、泊松比、拉梅常数和密度的乘积二维剖面(图12).由图12(a—c)可知，目的层具有低拉梅常数和密度的乘积(18～30 GPa·g·cm-3)、低泊松比(0.18～0.22)和较高的剪切模量(13～18 GPa)的特点，这些参数在横向上连续性较好.
图11 目的层过井二维测线地震反演剖面(a) 横波阻抗反演结果； (b) 纵波阻抗反演结果； (c) 密度反演结果.Fig.11 Target layer profile of the seismic inversion results of crossing-well 2D test line(a) S-wave impedance inversion results; (b) P-wave impedance inversion results; (c) Density inversion results.
图12 目的层过井二维测线地震反演剖面(a) 拉梅常数和密度的乘积反演结果； (b) 剪切模量反演结果； (c) 泊松比反演结果.Fig.12 Targe layer profile of the seismic inversion results of crossing-well 2D test line(a) λ ρ inversion results; (b) μ inversion results; (c) υ inversion results.
基于前述提出的优质页岩储层参数三维岩石物理模板，结合储层地震反演获取的密度、剪切模量、泊松比、拉梅常数和密度的乘积，开展对比分析，基于模板的储层参数范围，进行孔隙度、裂隙纵横比和石英含量的预测.对数据点与模板边界差距较大的情况按非储层处理.如图13给出了三维岩石物理模板预测的储层孔隙度、裂隙纵横比和石英含量的二维剖面和由弹性参数得到的脆性指数二维剖面，结果显示目的层的孔隙度变化范围是3～7%之间，且过A和B井的目的层孔隙度要高于过C井目的层孔隙度.目的层的整体裂隙纵横比较低，连通性好，裂隙发育，而C井裂隙纵横比较差.整个目的层的石英矿物含量较高，且过A和B井的目的层石英含量要高于过C井目的层的石英含量，同时图13d可以看出过A和B井的目的层脆性指数要高于过C井目的层的脆性指数，目的层的整体脆性较高.综上所述，目的层横向连续性好，优质页岩储层具有高孔隙度、高脆性(脆性矿物含量高)、裂隙纵横比小的特点，且A和B井目的层要优于C井.
图13 三维岩石物理模板预测过A、B和C井的二维剖面(a) 孔隙度； (b) 裂隙纵横比； (c) 石英含量； (d) 脆性指数.Fig.13 The predicted 2D profile by 3D rock physics templates crossing the wells A, B and C(a) Porosity; (b) Crack aspect ratio; (c) Quartz content; (d) Brittleness index.
利用上述方法可以得到工区约180 km2的三维地震预测结果.如图14给出了工区龙马溪组底部-五峰组的中间切片三维页岩储层孔隙度、石英含量、脆性指数和裂隙纵横比的预测结果.结果显示B井所处的区域目的层孔隙度、脆性指数以及石英含量最高，A井附近储层的预测结果参数次之，C井储层预测结果最差.B井附近储层裂隙纵横比较小，A井附近储层裂隙纵横比次之，C井裂隙纵横比略大.图15是TOC含量的三维预测平面分布，通过前文得出的TOC含量和孔隙度之间的拟合关系换算得到.结果显示TOC 含量较为发育，B井属于优质气井，TOC含量高，A井次之.在产气报告中，A井属于高产井，初产气量为20.56(万方/天)，B井初产气量是10.5(万方/天)，C井初产量为3.3(万方/天).综合分析得出的结论，B井的储层优于A井附近的储层，而B井的初产气量低于A井的产气量，这可能是由于压裂工艺对页岩气产量的明显影响.整体而言，目的层预测结果与实际情况基本一致，高孔隙度、高脆性(脆性矿物含量高)、裂隙纵横比小可有效指示优质储层，工区的甜点有利区分布，主要集中于西北方向.
图14 “甜点区”关键参数平面分布(a) 孔隙度; (b) 石英含量; (c) 脆性指数; (d) 裂隙纵横比.Fig.14 The map distribution of key parameters of “sweet spot”(a) Porosity; (b) Quartz content; (c) Brittleness index; (d) Crack aspect ratio.
图15 TOC含量的三维平面分布Fig.15 3D map distribution of TOC content
5 结论
本文对研究工区的页岩气储层进行了岩石物理特征分析，在矿物组分和有机质碳含量的关系中，随着测井深度增加，石英矿物含量与有机质含量具有较好的正相关性，黏土矿物含量与有机质含量具有负相关性，即高有机碳含量的储层对应石英含量高、黏土含量低.在分析有机质和脆性的参数关系中，综合评价认为优质页岩气储层具有高的有机碳含量(>3%)、低的拉梅常数和密度的乘积(18～30 GPa·g·cm-3)、低泊松比(0.18～0.22)以及较高的剪切模量(13～18 GPa)，同时具有孔隙度高(3%～7%)、密度低(小于2560 kg·m-3)、储层脆性高(脆性指数>50%)的特征.
根据页岩气储层的物性和岩石物理特征，利用EIAS理论模型，分析矿物组分、孔隙度和裂隙纵横比对优质页岩气储层弹性参数的影响，进一步提出适合页岩气储层描述的三维岩石物理模板.利用测井数据对构建的三维岩石物理模板进行校正，从而可对优质页岩气储层的孔隙度、石英矿物和裂隙纵横比进行定量预测.
基于研究工区的储层地震数据，采用适用于优质页岩储层的三维岩石物理模板，实现了二维测线和三维地震资料的孔隙度、石英含量和裂隙纵横比的定量预测.结果显示，过三口井的目的层孔隙度变化与测井结果一致性好，且裂隙纵横比小、石英含量高、目的层横向连续性较好，预测结果与实际产气情况基本一致.根据实际地层岩石物理特征构建的三维岩石物理模板，可较好地表征页岩气“甜点”，从而结合地震参数反演估算矿物含量、孔隙度及裂隙参数，进而有效地指示优质页岩储层分布.
