致密砂岩脆性指数测井评价方法——以鄂尔多斯盆地陇东地区长7段致密砂岩储集层为例

致密砂岩脆性指数测井评价方法——以鄂尔多斯盆地陇东地区长7段致密砂岩储集层为例李华阳;周灿灿;李长喜;李潮流;石玉江;王长胜【摘要】岩石的脆性指数是致密砂岩油气体积压裂设计中应考虑的重要因素之一.通过对岩石弹性参数法和矿物成分法两种脆性指数测井评价方法的分析,并应用岩石力学实验和全岩X衍射分析实验标定,建立了应用常规测井的多矿物模型确定致密砂岩脆性指数的方法.对鄂尔多斯盆地陇东地区长7段致密砂岩储集层实际资料处理结果表明,应用常规测井资料对岩石脆性指数进行评价,可以弥补由于阵列声波(或声波扫描)测井采集井数较少难以开展岩石脆性指数平面分布研究的不足,为致密砂岩油气“甜点”优选提供技术支持.【期刊名称】《新疆石油地质》【年(卷),期】2014(035)005【总页数】5页(P593-597)【关键词】致密砂岩储集层;脆性指数;多矿物模型;岩石力学【作者】李华阳;周灿灿;李长喜;李潮流;石玉江;王长胜【作者单位】中国石油勘探开发研究院,北京100083;中国石油勘探开发研究院,北京100083;中国石油勘探开发研究院,北京100083;中国石油勘探开发研究院,北京100083;中国石油长庆油田分公司勘探开发研究院,西安710065;中国石油长庆油田分公司勘探开发研究院,西安710065【正文语种】中文【中图分类】TE125.1;P631.824致密砂岩油气在中国分布广泛且潜力巨大，但由于储集层致密、岩性复杂，均需采用压裂改造以获得工业产能[1-2]，因此，可压性评价对于优选压裂井段具有重要意义。

文献［3］通过页岩脆性指数表征岩石的可压性，为致密储集层的可压性评价开辟了思路。

岩石的脆性指数是致密砂岩油气体积压裂设计中应考虑的重要因素之一，也是测井“七性关系”研究的重要内容之一[4]。

国内外针对非常规储集层岩石脆性指数开展了大量研究工作，在测井评价中主要采用岩石弹性参数计算法和岩石矿物组分计算法计算岩石脆性指数。

但由于对岩石脆性的物理意义及实验测量、岩石脆性指数与岩石弹性参数关系、岩石脆性指数与矿物组分关系的定量分析和研究不够，一定程度上影响了脆性指数在非常规储集层测井评价中的应用效果。

本文通过岩石力学实验和全岩X衍射分析实验标定，分析了岩石脆性指数与岩石弹性参数以及矿物成分之间的关系，建立了应用常规测井的多矿物模型确定致密砂岩脆性指数的方法，并通过实例应用，证明该方法可为优选致密砂岩储集层的有利压裂层段提供重要依据，为致密砂岩油气的工程“甜点”优选提供技术支持。

在非常规储集层测井评价中，通常应用岩石脆性表示储集层岩石的压裂难易程度，但当前尚无关于岩石脆性的明确定义及岩石物理测量方法，多用于定性描述中，并且不同地区不同脆性指数计算方法的计算结果差异很大。

一般认为，脆性指数与岩石矿物成分、岩石力学性质等密切相关，杨氏模量越大，泊松比越小，岩石脆性指数越高，越容易破裂形成复杂裂缝[5]。

脆性指数较高的储集层一般性质硬脆，受构造运动的影响，其天然裂缝发育，对压裂作业非常敏感，能迅速形成复杂的网状裂缝，对提高产能有很好的帮助，而脆性指数低的储集层一般形成简单的双翼型裂缝[6-8]。

目前，实验室研究岩石脆性的主要手段是应力—应变实验、静态弹性参数实验以及岩石矿物成分分析实验等[9-10]。

根据应力—应变实验分析，当应力由某一初始弹性态加载到峰值强度后，将发生突变而迅速跌落至残余强度面上。

岩石应力跌落的这一特征，被称为应力脆性跌落系数或脆性指数。

通过选取三轴实验应力—应变图中脆性较强的岩样（图1），本文应用岩石变形参数法来定义脆性指数BI，岩石变形参数表征脆性指数公式为式中a，b——分别为不可恢复应变和残余应变。

根据鄂尔多斯盆地某区块6口井长7段致密油层岩石力学实验结果，杨氏模量和泊松比与脆性指数具有很好的相关性（图2，图3），可较好表征岩石的脆性指数。

在页岩气脆性评价研究中，国外学者分析了不同矿物岩石的力学性质差异，认为矿物成分及其结构等对岩石脆性有直接影响，并建立了应用岩石矿物成分计算岩石脆性的测井评价模型。

由于不同学者对脆性矿物的认识差异（一般认为石英、方解石、长石、白云石等为脆性矿物，而黏土为塑性矿物），所建立的脆性计算模型也各有不同。

为分析致密砂岩储集层中各矿物组分与脆性指数关系，根据鄂尔多斯盆地某区块长7段的岩石力学实验和全岩X衍射分析实验结果发现（图4，图5），石英含量与岩石变形参数法确定的脆性指数之间有较好的正相关关系，而黏土含量与岩石变形参数法确定的脆性指数之间存在较好的负相关关系。

因此，可应用岩石矿物的石英含量来表征砂岩储集层岩石的脆性指数。

对于不同地区由于岩石成岩作用的差异性，需要通过岩石物理实验标定来建立相应的岩石矿物成分及含量与脆性指数之间的相关关系。

由鄂尔多斯盆地某区块的岩石力学实验结果（图2、图3）可知，根据实验测量的泊松比和杨氏模量与脆性指数之间的关系，可建立岩石弹性参数的脆性指数计算式式中BI——脆性指数，%.ΔEy——归一化处理后的杨氏模量，GPa；Δμ——泊松比相对值。

为了应用（2）式计算岩石脆性指数，杨氏模量和泊松比等岩石弹性参数的获取是关键。

为获取高精度的杨氏模量和泊松比等岩石弹性参数，测井采集系列中应配备高精度的密度测井和阵列声波测井（或声波扫描测井），同时，以岩石力学实验测量值刻度测井计算值。

岩石弹性参数的常用测定方法分为静态法和动态法。

静态弹性参数是通过岩石力学实验直接获取岩石属性参数，包括泊松比、杨氏模量、抗压强度、抗张强度、抗剪强度和内摩擦角等；动态弹性参数是通过测定岩样中的纵波和横波传播速度并结合体积密度转换后得到的，其计算公式如下：式中Δtp——纵波声波时差，μs/m；Δts——横波声波时差，μs/m；ρ——体积密度，g/cm3.需要注意的是，实验室测量值是静态弹性参数，测井计算值是动态弹性参数，二者之间存在较大的差异，而且这种差异随地层力学性质和实验条件的不同而不同，通常采用回归拟合关系来实现弹性参数的动静态转换[11]。

根据研究区岩石力学实验结果，泊松比和杨氏模量的转换公式分别为式中μs,μd——分别为静态和动态泊松比；Ey,s,Ey,d——分别为静态和动态杨氏模量，GPa.应用矿物成分计算岩石脆性指数的计算公式为式中Vqa，Vca，Vdo，Vcl——分别为石英、方解石、白云石和黏土的含量；K——地区修正系数。

在岩石脆性指数测井评价中，若有高精度的密度测井和阵列声波测井，可应用（2）式直接计算岩石脆性指数，应用实验室刻度进行岩石动、静态弹性参数转换是计算的关键。

但实际测井中，大多数井只测量了常规测井，无横波测井资料，难以应用岩石弹性参数法计算脆性指数。

根据前述岩石力学实验和全岩分析实验结果以及（9）式可知，矿物成分计算法在研究区储集层脆性指数计算中也具有较好的效果。

因此，针对常规测井资料，如何获得较准确的岩石矿物剖面是评价脆性指数的关键。

传统的砂泥岩剖面对于储集层划分和有效性评价具有很好的辅助作用，但利用自然伽马曲线计算得到的砂质和泥质主要反映粒度大小，难以反映矿物成分含量（石英、黏土等）。

因此，在致密油脆性指数评价中，为应用矿物成分法计算，需使用多矿物模型[12]对该地区的矿物组分进行精细处理。

多矿物模型是将复杂岩性储集层看成是由若干种骨架矿物、黏土矿物和孔隙流体等几部分组成，而地层的测井值就是这几部分的综合响应，其响应方程为其中χi——各矿物的密度、声波、中子、自然伽马等的测井响应值；Vi——各矿物的体积含量。

根据研究区矿物成分分析资料确定主要矿物组成，建立岩石物理解释模型，利用密度、声波、中子、伽马以及电阻率五条常规测井曲线建立测井响应方程组，采用最优化处理方法进行求解[13]。

图6中计算结果对比道为A井应用多矿物模型计算的岩石矿物组分质量分数与全岩X衍射分析结果对比（计算结果对比道中点数据为全岩X衍射分析结果，曲线为矿物的计算结果）具有很好的一致性，反映多矿物模型和组分骨架参数选择合理，计算结果可靠，可精细解释岩石的矿物剖面，为应用矿物成分法计算岩石脆性指数奠定了基础。

为分析两种脆性指数测井评价效果，选择研究区资料齐全的B井进行处理（图7）。

图7中第七道为岩石力学弹性参数计算的脆性指数剖面，第二道、第三道分别为应用最优化计算的岩石矿物组分剖面和应用单矿物模型计算的砂泥岩剖面，对比可看出二者有较大差异；第六道黄色曲线为应用砂质含量计算的脆性指数同岩石弹性参数法计算结果的对比，可看出砂质含量难以反映储集层岩石的脆性指数；而第五道中绿色曲线是根据多矿物模型计算结果应用矿物成分法计算岩石脆性指数，与应用岩石弹性参数法计算脆性指数一致性好，可有效指示纵向剖面上储集层岩石脆性指数的差异，为有利压裂层段选择提供依据。

基于矿物组分的脆性指数计算法实现了在无阵列声波测井资料情况下的岩石脆性指数评价，通过岩心刻度测井，具有较好的计算精度。

并且矿物成分计算法纵向上连续性好，井间对比方便，可在平面上进行规模化处理，便于脆性指数平面成图。

（1）脆性指数可以表征岩石压裂的难易程度，脆性指数越大，对压裂反应越灵敏，形成的裂缝越复杂，岩石的脆性指数是致密砂岩油气体积压裂设计中的重要因素之一，是测井“七性关系”研究的重要内容之一。

（2）通过岩石力学实验和全岩X衍射分析实验标定，分析了岩石脆性指数与岩石弹性参数以及矿物成分之间的关系，建立了鄂尔多斯盆地陇东地区长7段岩石弹性参数计算法和矿物成分计算法两种脆性指数测井评价方法。

（3）应用常规测井资料，通过岩心实验标定建立了多矿物最优化处理模型和基于矿物成分的脆性指数计算方法，与岩石弹性参数计算法得到的脆性指数一致性好，对鄂尔多斯盆地陇东地区长7段致密砂岩储集层实际资料处理结果表明该方法具有较好的应用效果，为优选致密砂岩储集层的有利压裂层段提供了重要依据。

（4）基于矿物成分计算法的脆性指数测井评价弥补了由于阵列声波（或声波扫描）测井采集井数较少难以开展岩石脆性指数平面分布研究的不足，可进行批量处理，便于脆性指数平面成图，为致密砂岩油气的工程“甜点”优选提供了技术支持。

【相关文献】［1］李斌,朱永铭,童孝华,等.靖安油田上三叠统长6储集层成岩作用［J］.新疆石油地质,2004,24（4）:292-295.Li Bin,Zhu Yongming,Tong Xiaohua,et al.Diagenesis of Chang⁃6 reservoir of Upper Triassic in Jing'an oilfield of Ordos basin［J］.Xin⁃jiang Petroleum Geology,2004,24（4）:292-295. ［2］邹才能，朱如凯，吴松涛，等.常规与非常规油气聚集类型、特征、机理及展望——以中国致密油和致密气为例［J］.石油学报，2012，33（2）：173-187.Zou Caineng,Zhu Rukai,Wu Songtao,et al.Types,characteristics, genesis and prospects ofconventional and unconventional hydrocar⁃bon accumulations:taking tight oil and tight gas in China as an in⁃stance［J］.Acta Petrolei Sinica,2012,33（2）173-184.［3］Chong K K,Grieser W V,Passman A.A completions guide book to shale⁃play development:a review of successful approaches towards shale⁃play stimulation in the last two decades［R］.SPE 133874, 2010.［4］赵政璋,杜金虎.致密油气［M］.北京:石油工业出版社,2013.Zhao Zhengzhang，Du Jinhu.Tight oil and gas［M］.Beijing:Petro⁃leum IndustryPress,2013.［5］Rickman R，Mullen M，Peter E，et al.A practical use of shale petro⁃physics for stimulation design optimization:all shale plays are not clones of the Barnett shale［R］.SPE 115258，2008.［6］Kahraman S,Altindag R.A brittleness index to estimate fracture toughness［J］.International Journal of Rock Mechanic&Mining⁃Science,2004（41）:343-348.［7］黄金亮，邹才能，李建忠，等.川南志留系龙马溪组页岩气形成条件与有利区分析［J］.煤炭学报，2012,37（5）:782-787.Huang Jinliang,Zou Caineng,Li Jianzhong,et al.Shale gas accumu⁃lation conditions and favorable zones of Silurian Longmaxi formation in south Sichuan basin，China［J］.Journal of China Coal Society, 2012,37（5）:782-787.［8］袁俊亮,邓金根,张定宇,等.页岩气储集层可压裂性评价技术［J］.石油学报，2013，34（03）：523-527Yuan Junliang,Deng Jingen,Zhang Dingyu,et al.Fracability evalu⁃ation of shale⁃gas reservoirs［J］.Acta Petrolei Sinica,2013,4（3）：23-527.［9］蒋裕强，董大忠，漆麟，等.页岩气储集层的基本特征及其评价［J］.天然气工业，2010,30（10）:7-12.Jiang Yuqiang,Dong Dazhong,Qi Lin,et al.Basic features and eval⁃uation of shale gas reservoirs［J］.Natural Gas Industry,2010,30（10）:7-12.［10］陈吉，肖贤明.南方古生界3套富有机质页岩矿物组分与脆性分析［J］.煤炭学报，2013，38（05）：822-826.Chen Ji,Xiao Xianming.Mineral composition and brittleness of three sets of Paleozoic organic⁃rich shales in China south area［J］. Journal of China Coal Society,2013,38（05）:822-826.［11］贺顺义,师永民,谢楠,等.根据常规测井资料求取岩石力学参数的方法［J］.新疆石油地质,2008,29（5）:662-664.He Shunyi,Shi Yongmin,Xie Nan,et al.The method for acquire⁃ment of conventional logging response⁃based lithomechanical pa⁃rameters［J］.Xinjiang Petroleum Geology,2008,29（5）:662-664.［12］杨宇，康毅力，康志宏.多矿物测井解释模型以及其在砂泥岩地层测井解释中的应用［J］.新疆地质，2006，24（04）：441-446.Yang Yu,Kang Yili,Kang Zhihong.Analysis of multi⁃mineral mod⁃el in log interpretation of detrital strata［J］.Xinjiang Geology, 2006,24（4）:441-446.［13］肖亮,毛志强,孙中春,等.最优化方法在复杂岩性储集层测井评价中的应用［J］.断块油气田,2011,18（3）:342-345.Xiao Liang,Mao Zhiqiang,Sun Zhongchun,et al.Application of op⁃timization method in log evaluation of complex lithologic reservoir［J］.Fault⁃Block Oil&Gas Field,2011,18（3）:342-345.。