基于灰色关联与层次分析的脆性指数预测方法——以准噶尔盆地吉木萨尔凹陷芦草沟组致密储层为例.pdf

第 47 卷第 4 期物 探 与 化 探Vol.47,No.4 2023 年 8 月GEOPHYSICAL&GEOCHEMICAL EXPLORATION Aug.,2023doi:10.11720/wtyht.2023.1242刘庆,张镇,杨帅,等.基于灰色关联与层次分析的脆性指数预测方法 以准噶尔盆地吉木萨尔凹陷芦草沟组致密储层为例J.物探与化探,2023,47(4):944-953.http:/doi.org/10.11720/wtyht.2023.1242Liu Q,Zhang Z,Yang S,et al.Method for brittleness index prediction based on grey correlation and analytic hierarchy process:A case study of the tight res-ervoirs in the Lucaogou Formation of the Jimusaer Sag,Junggar BasinJ.Geophysical and Geochemical Exploration,2023,47(4):944-953.http:/doi.org/10.11720/wtyht.2023.1242基于灰色关联与层次分析的脆性指数预测方法 以准噶尔盆地吉木萨尔凹陷芦草沟组致密储层为例刘庆1,张镇2,杨帅1,李枫凌1(1.中国石油 新疆油田分公司 油田技术服务分公司,新疆 克拉玛依 834000;2.中国石油 新疆油田分公司 克拉玛依红山油田有限责任公司,新疆 克拉玛依 834000)摘 要:吉木萨尔凹陷致密储层因储层物性差,需要水力压裂才能获得工业油流。开展对应的岩石力学特性研究及脆性评价可为水力压裂工艺提供参考依据。本文通过三轴力学试验获得吉木萨尔致密地层的力学特性。开展脆性指数的测井敏感性相关性分析,确定对岩石脆性潜在敏感的测井参数,依据灰色关联理论,确定敏感参数的初始序列并对筛选的参数归一化处理,继而将筛选出的敏感参数对脆性指数潜在敏感性进行量化关联分析,获得各参数对脆性指数的关联度大小,确定各关联度的次序。在此基础上,基于层次分析法构造敏感参数的两两比较判断矩阵,确定权重向量,构建了脆性指数与敏感参数间的函数关系模型,从而建立了新的脆性指数预测模型,将测井模型与基于力学特性的脆性模型、室内测试的脆性指数进行比较。研究表明:吉木萨尔致密地层具有较强的脆性特征,基于全过程应力应变曲线表征的综合脆性指数符合实际岩石的脆性特征:灰色关联法确定的对脆性指数敏感的测井参数的关联度大小的顺序依次为自然伽马(GR)、电阻率(Rt)、密度()、中子(CNL)、声波时差(T);层次分析法构建的新预测模型中上述敏感参数的权重系数分别是 0.33、0.22、0.18、0.16、0.11;将该法应用于吉木萨尔芦草沟组,预测的脆性指数与实验室测试的相对比,能够较真实地反映实际地层的脆性,吻合度高。试油结果表明,米采油指数与脆性指数呈正比,脆性指数高的层段,压裂后的产能高。因此,新方法为脆性指数预测提供了新的思路,为储层的压裂改造提供了参数指导。关键词:灰色关联;层次分析;脆性指数;致密储层;芦草沟组中图分类号:P631.4 文献标识码:A 文章编号:1000-8918(2023)04-0944-10收稿日期:2022-05-20;修回日期:2023-03-21基金项目:国家自然科学基金项目(51974332)第一作者:刘庆(1987-),男,新疆油田公司地质勘探专业工程师,本科学历,主要研究方向为测井、地质勘探。Email:dxlq 通讯作者:张镇(1990-),男,新疆油田公司油气田开发(地质)专业工程师,工商管理硕士,地质勘探学士。Email:hs_zhz 0 引言为了满足国家经济的高速发展,对石油能源的需求逐年增加。而国内常规油气增产难度越来越大,已无法满足工业需求。致密油气资源潜力大,分布广泛,近年来已越来越受到高度重视1-3。致密油气储层因非均质性强、孔隙度与渗透率极低等特征,加大了开采的难度,致使其获得高效的产能面临巨大挑战。为了突破该瓶颈,水平井分段体积压裂进行储层改造应用越来越广泛。而岩石脆性是评价地层可压性进而制定有效的体积压裂方案的重要因素。脆性指数在钻探工艺、煤层气开采、地下工程中也逐渐受到重视。在致密储层开采中,岩石的脆性指数反映了岩石破裂的难易程度,对于脆性并没有一个统一的定义。通常认为,脆性指数高的地层,内部的天然裂缝网络越发育,能够形成复杂的裂缝网络。目前获取脆性指数的方法总体分为两类,即室内测试与测井预测。室内测试法主要依据岩石力学试验,如单轴、4 期 刘庆等:基于灰色关联与层次分析的脆性指数预测方法 以准噶尔盆地吉木萨尔凹陷芦草沟组致密储层为例三轴试验,静态弹性参数测试及矿物测试分析4-6。尽管室内测试可以获得精确的脆性指数,但测试流程繁琐、依赖钻孔取心、测试费用昂贵,并且无法获得连续的脆性指数剖面。相比之下,测井预测方法的流程简单,花费的物力较少,获得的数据也能够真实接近岩石的脆性特征。测井预测法主要分为矿物组分法和弹性参数法。矿物组分法利用测井方法计算地层的石英、方解石等脆性矿物的含量,进而获得脆性指数。而弹性参数法则通过纵、横波时差转换和密度测井资料,获得杨氏模量和泊松比,然后得到脆性指数。矿物组分法受组分矿物计算的精度的影响,存在一定的局限性。弹性参数法则需要横波资料,然而油田现场并非所有井都有横波测井资料,因此该方法也难以全面推广。灰色关联法本质上是用有限数列解决无限空间问题,以获得系统内各因素的关联度,从而可知各因素对目标函数的重要程度,进而避免人为确定各因素指标权重的主观性7-9。近年来,灰色关联法在致密砂岩油藏压后产量预测、油藏开发方案优选方面有研究,而在脆性指数预测方面研究较少。层次分析法(AHP)是一种把复杂问题中的各因素通过划分为相互联系的有序层次且支持多准则决策的方法,在可压裂性评价方面应用较多10-11,目前在脆性指数预测方面也逐步有少量应用。而将灰色关联理论与层次分析法结合预测脆性指数目前还没有相关研究。本文以准噶尔盆地吉木萨尔凹陷芦草沟组地层为例,基于常规测井资料,首先将测井参数与脆性指数进行相关性分析,从定性的角度筛选出对脆性指数敏感的参数;紧接着通过灰色关联法确定敏感参数的初始序列,从定量的角度获得潜在敏感参数与脆性指数的关联度大小,获得对脆性指数敏感程度的排序;进而基于层次分析法构造两两判断矩阵,确定权重向量,构建脆性指数与敏感测井参数的函数关系模型;最终建立新的脆性指数预测模型,以实现对脆性指数的连续预测。将上述方法应用于实际井中,将预测的结果与实验室测试结果对比,并借助水力压裂过程中的微震监测结果辅助分析,证实该方法可以精确地预测脆性指数,为致密储层水力压裂提供参数依据。1 研究区地质概况及储层特征吉木萨尔凹陷位于准噶尔盆地东南方位,南、北方向分别以三台断裂、吉木萨尔断裂为界,西部方向以西地断裂为界,东部方向以奇台凸起为界,呈箕状凹陷(图 1),面积约为 1 278 km2。区内二叠系芦草沟组厚度大,平均厚度约为 300 m 左右,整体呈南厚北薄、西厚东薄的趋势,分布较为广泛,分布面积为900 km2左右。芦草沟组地层自下而上共分为芦草沟组一段(P2l1)和二段(P2l2)两套致密型砂泥岩正旋回储盖组合,表现为源储一体、近源成藏、纵向上整体含油的特征12-13。芦草沟组地层平均孔隙度为 14.3%,平均渗透率为 0.1 mD(图 2),为典型的致密储层。图 1 吉木萨尔凹陷构造位置Fig.1 Structural location map of Jimusaer Sag图 2 芦草沟组地层物性分布Fig.2 Physical property distribution of Lucaogou Formation549物 探 与 化 探47 卷 2 现有脆性评价方法统计分析目前脆性指数评价方法可分为 5 大类,共有 22种方法(表 1),总体来讲,主要是通过实验室分析与测井预测两种途径获得14-17。BI1 BI16主要通过岩石力学实验获取,依据岩石压裂后的特征及压裂过程中的测量数据进行评价,尽管这种方式会比较准确,但对于产层较多、较薄的致密储层,对每个小层的岩心样品进行实验测量需要昂贵费用与较长周期,同时实验获得的结果不连续,这种方法存在较大局限。BI17BI18既可通过实验测试也可利用测井资料预测获得,BI19是通过测井资料预测获得,BI20BI22可通过 XRD 测试或者矿物含量测井获得。虽然测井预测方法存在一定多解性,但因测井数据具有连续性强、分辨率高的特点,实际过程中可以将多解性减少到可接受的范围。因此,相比于室内试验,这类预测方法具有经济、实效与精度高的优势。表 1 现有脆性指数评价方法统计Table 1 Statistics of brittleness index evaluation method原理分类公式变量说明获取方法基干硬度或坚固性BI1=(Hm-H)/KH 为硬度,GPa;Hm为微观硬度,GPa;K 为体枳模量,GPaBI2=H/KICKIC为断裂韧性,MPam1/2BI3=HE/K2ICE 为静态杨氏模量,GPaBI4=qcq 为直径小于 0.6 mm 碎屑百分比,%;c为抗压强度,MPaBI5=c/dc 为裂纹长度,m;d 为韦氏测试特定载荷下贯入尺寸,mBI6=tdec/tinctdec为平均载荷减少时间,s;tinc为平均载荷增加时间,sBI7=Fmax/PFmax为试件所受最大载荷,kN;P 为相应的贯入深度,mm硬度测试硬度和韧性测试陶质材料测试兽氏冲击实絵贯入实验基于强度比值BI8=c/tBI9=(c-t)/(c+t)BI10=ct/2BI11=(ct)0.5/2c为抗压强度,MPa;t为抗拉强度,MPa单轴抗压測试和巴西劈裂实验基于全应力应变特征BI12=(pr)/pBI13=r/tBI14=Wr/WtBI15=ux100%BI16=(p-r)/pBI17=/4+/2BI18=sinp为剪切强度峰值,MPa;r为残余剪切强度,MPar为可恢复应变,无量纲;t为总应变,无量纲Wr为可恢更应变能,J;Wt为总应变能,Jux为不可恢复轴向应变,无量纲p为应变峰值,无量纲;r为残余应变,无量纲 为内摩擦角,rad应力应变测试应力应变测试或声波测井数据基于弹性力学参数BI19=(En+vn)/2En为归一化杨氏模量,无量纲;vn为归一化泊松比,无量纲密度与声波测井基于岩石矿物组分BI20=Wqtz/WotBI21=(Wqtz+Wdot)/WotBI22=(WQFM+WCar)WotWqtz为石英含量,%;Wot为矿物总量,%Wdot为白云岩含量,%WQFM为硅酸盐岩含量,%;WCar为脆性碳酸盐岩含量,%实验室 XRD 测试或矿物含量测井3 致密储层脆性指数实验分析采用美国生产的 MTS 全伺服三轴岩石力学测试系统,对致密储层标准岩样开展岩石力学实验,标准圆柱岩样的尺寸为 50 mm25 mm。加载过程中,围压预定值设置为 40 MPa,当围压达到了预设值后,开展轴向加载至岩石破坏,加载速率为 0.001 mm/s,当岩石破坏后,轴向加载速率改变为 0.001 5 mm/s,加载过程中对应力、应变数据进行记录。本文采用基于三轴力学实验中的全过程应力应变曲线特征来表征岩石的脆性特征,该方法既能够反映峰前力学特性,也能够反映峰后对应的特征。该方法中峰值应变值被用来表征脆性破坏的难易程度,而峰后曲线形态特征被用来表征脆性强弱,方法流程如图 3 所示。图 3 基于全过程应力应变的脆性指数计算方法Fig.3 Calculation method of brittleness index based on stress-strain in whole process 该方法将脆性指数分为两部分,即峰值应变指数 BIP、峰后形态指数 BIS,可用下式来表达:BIe=BIP+BIS,(1)式中,BIe为实验室方法获得的脆性指数。峰值应变指数 BIP求取方法如下:BP=P,(2)649 4 期 刘庆等:基于灰色关联与层次分析的脆性指数预测方法 以准噶尔盆地吉木萨尔凹陷芦草沟组致密储层为例BIP=BP-NM-N,(3)式中:BP为峰值应变参变量;P为峰值应变量;M、N分别对应单轴压缩实验中的峰值应变极值,M=max BP(),N=min BP(),根据测试的岩样数量决定,吉木萨尔凹陷芦草沟组取值为 0.02、0.001 8。峰后形态指数 BIS由公式CPBP=PP-r()Pr-P()=P-rEPr-P(),(4)BIS=CPBP+kCPBP+,(5)计算获得。式中:P、r对应峰值强度与残余强度;P、r对应峰值及残余应变量;EP为峰值强度点对应的割线模量;CPBP对应峰后应力应变曲线形态参数;、k、为标准化系数,通常取值范围为 01,其取值取决于岩石类型与测试环境。此次实验采用的岩样来自于新疆吉木萨尔凹陷芦草沟组,取值情况为:=0.031,k=0.102,=0.058。式(2)(5)中其他的参数取值如图 4、图 5、表 2 所示。由图 6 可以看出,吉木萨尔凹陷致密地层表现出很强的脆性特征,在围压一定、轴向压力达到峰值时,泥晶白云岩受压后的裂纹呈现张性单霹裂特征,泥页岩呈现为共轭状裂缝面,砂屑白云岩为单剪切面破坏,该剪切面贯穿整个样品,由此看出,这 3 类图 4 脆性指数各参数取值示意Fig.4 Schematic diagram of brittleness index parameters图 5 部分岩样的三轴应力应变曲线Fig.5 Triaxial stress-strain curves of some rock samples表 2 利用实验方法获得的岩样综合脆性指数Table 2 Comprehensive brittleness index of rock samples obtained by experimental method编号岩性P/MPar/MPaP/%r/%BIPBISBIeX1泥晶白云岩228.3121.91.31.70.6150.2590.875X2泥页岩178.6115.31.11.50.5050.1880.693X3砂屑白云岩201.988.70.60.80.2310.2820.512X4白云质粉砂岩162.389.30.50.90.1760.1330.308图 6 岩样被压缩后的形态Fig.6 Morphology of rock sample after compression岩性具有较强的脆性特征。相比于其他岩性,白云质砂岩脆性特征相对较弱,其剪切破坏并未贯穿样品。经过计算后的综合脆性指数可以看出,泥晶白云岩的脆性指数最高,为 0.875;泥页岩的综合脆性指数为 0.693,砂屑白云岩为 0.512,白云质粉砂岩为 0.308。因此,基于全过程应力应变的综合脆性指数可以表征实际的储层脆性特征,可以作为一种较精确的实验室测量值。4 脆性指数与测井参数的关联程度4.1 脆性指数敏感因素分析4.1.1 测井参数与脆性指数相关性分析目前对声波时差、密度测井与脆性指数之间的关系研究的较多,主要是由于这两个参数与杨氏模量、泊松比有关,而杨氏模量与泊松比常用于脆性指数计算。但脆性指数与其他测井参数,如自然伽马、电阻率、中子孔隙度之间的关系很少被研究。文中选取 J13 井、J25 井的自然伽马、声波时差、密度、电阻率、中子孔隙度与实验室所测定的脆性指数做交749物 探 与 化 探47 卷 会,探索这几个测井参数与脆性指数之间的关系。自然伽马主要反映泥质含量,图 7a 表明自然伽马与脆性指数呈负相关关系,自然伽马值越高,泥质含量越高,脆性越小。研究表明,脆性指数与孔隙度具有很好的相关关系。图 7b 表明中子值增大,脆性指数降低,中子值与脆性指数呈负相关关系。图 7c 表明,脆性指数随密度值增大而增大。图 7d 表明,脆性指数随声波时差的增大而减小。图 7e 表明,电阻率与脆性指数呈正相关关系,随着电阻率的升高,脆性指数增大。这主要由于地层越致密,孔隙连通性就越弱,导电性就越弱。而图 7c 表明地层密度越大,脆性指数越高,因此,可以得出电阻率与脆性指数呈正相关关系。综上,可以看出脆性指数对自然伽马、电阻率、声波时差、中子孔隙度、密度敏感。图 7 实验室测定的脆性指数与测井参数交会Fig.7 Cross plot of logging parameters and brittleness index measured in laboratory4.1.2 灰色关联理论的敏感参数关联程度确定利用灰色关联法确定敏感测井参数与脆性指数之间的关联度,从而判别各因素对脆性指数敏感程度的大小,该理论的基本步骤如下:1)确定初始序列。所研究问题的因变量和多个自变量因素构成了初始序列。因变量数据构成参考数列;每个自变量因素的多个样本数据构成一个比较数列。X0=X0(1),X0(2),X0(3),X0(n),(6)Xi=Xi(1),Xi(2),Xi(3),Xi(n),(i=1,2,m),(7)式中:X0为参考数列;Xi为比较数列;i 为比较数列849 4 期 刘庆等:基于灰色关联与层次分析的脆性指数预测方法 以准噶尔盆地吉木萨尔凹陷芦草沟组致密储层为例的编号;m 为自变量因素的个数;n 为样品的个数。2)样本数据归一化。因变量、自变量数据所组成的数列量纲存在差异,无法直接进行比较,采用极值变换法进行归一化处理。以归一化后的脆性指数为参考数列,以归一化后的各敏感测井参数为比较序列。X-ik()=Xik()-min Xik()max Xik()-min Xik(),(i=1,2,m;k=1,2,n)。(8)3)计算差序列、最大差、最小差。样本数据无量纲化后,计算归一化后的第 i 个比较数列中各样本参数值与参考数列中对应的参数值的绝对差值,公式如下:0ik()=X-ik()-X-0k(),(9)计算所获得的绝对差值的最大、最小值分别为max=maxi=1,2,mmaxk=1,2,n0ik(),(10)min=mini=1,2,mmink=1,2,n0ik()。(11)4)确定关联系数。依据式(9)对绝对差值阵进行变换:0ik()=min+max0ik()+max,(12)获得关联系数矩阵:011()021()0m1()012()022()0m2()01n()02n()0mn(),(13)式(12)中:为分辨系数,0 1,表征各关联系数之间差异性的大小,越小,关联系数间的差异越明显,通常取为 0.5。0ik()是比较数列第 k 个元素 Xi(k)与参考数列对应的元素 X0(k)的相对差值,称为比较数列 Xi对参考数列 X0在第 k 个元素上的关联系数。通常,0ik()为不超过 1 的正数,0ik()越小,0ik()越大,体现了第i个比较数列Xi与参考数列在第k 个样本数据上的关联程度。5)确定关联度。单个关联系数所代表的信息比较分散,不具有比较性,不能反映各敏感因素对脆性指数的关联程度。通常用 n 个关联系数的平均值表征数列之间的关联强度大小:0i=1nnk=10ik(),(14)式中,0i表示第 i 个比较数列与参考数列的关联程度的大小。关联度越大,比较数列与参考数列的变化越一致。6)确定权重及关联度排序。将关联度进行归一化,获得各比较数列关联度的权重 Wi:Wi=0imi=10i,(15)进而确定各敏感因素与脆性指数之间关联程度的相对大小,依据权重数值大小排序,就可以反映各个敏感因素之间的敏感程度。按照上述步骤,确定各敏感参数与脆性指数之间的关联度大小(表 3),可以看出脆性指数与自然伽马的关联度最大,由此可以看出脆性指数与黏土含量的关系最为密切。表 3 各因素与脆性指数的关联度及排序Table 3 Correlation degree and order of each factor and brittleness index因素自然伽马GR/API电阻率 Rt/(m)密度/(gcm-3)中子 CNL/0.01声波时差T/(ss-1)关联度0.8610.7520.7060.5650.543排序123454.2 脆性指数预测模型建立层次分析法(AHP)是一种把复杂问题中的各因素通过划分为相互联系的有序层次且支持多准则决策的方法。目前主要应用于确定综合评价的权重系数10-11。基于标准化的测井参数,依据波叠加原理构建目的层段脆性指数测井预测模型。依据 5 个归一化的参数,构建两两比较判断矩阵,计算权重向量及开展一致性检验,从而实现对影响脆性指数的因素综合评价。1)测井参数标准化。因各测井参数量纲不同,对评价会产生影响。因此,采用极差标准化方法,对敏感测井参数归一化。依据式(16),对 T、GR、Rt、CNL 等 5 项对储层脆性指数敏感的参数进行归一化处理,归一化后的参数量纲一致,且各参数间的相关性不变。由图2 可知,和 Rt为正向指标参数,脆性指数的大小与正向指标值的变化成正比,T、CNL、GR为负向指标参数,脆性指数的大小与负向指标值的变化成反比。正向指标归一化公式如下:Zi=Xi-min Xi()max Xi()-min Xi(),(16)式中:Zi为各参数归一化后的值,无量纲;Xi为测井参数原始值;Xmax、Xmin为目的层段各测井参数最大值、最小值。负向指标归一化公式:949物 探 与 化 探47 卷 Zi=max Xi()-Ximax Xi()-min Xi(),(17)式中:Zi为各参数归一化后的值,无量纲;Xi为测井参数原始值;Xmax、Xmin为目的层段各测井参数最大值、最小值。2)脆性指数测井预测模型建立。基于灰色关联分析的结果及敏感参数归一化,利用叠加原理构建了脆性指数测井预测模型,从而开展对致密储层的脆性评价。BI=ni=1si i,(18)式中:BI 为脆性指数,%;si为各参数的权值,无量纲;i为各参数,无量纲;n 为测井参数个数。因此式(18)可演化为BI=s1 GRN+s2 Rt,N+s3 N+s4 CNLN+s5 TN,(19)式中:BI 为脆性指数,%;N 为标准化后的参数;其他参数物理意义同上。3)权系数确定。脆性评价涉及多个测井参数,通过灰色关联法已经确定出各测井参数与脆性指数的关联度大小,在此基础上,结合层次分析法确定出各敏感测井参数在脆性指数预测模型中的权重。层次分析法(AHP)是一种层次权重决策分析方法,该方法将复杂问题进行分解,把目的问题分解为若干个影响因素,每个因素根据属性不同继续分解,从而形成层次结构,并建立判断矩阵确定各因素的权重。结合灰色关联确定的各敏感测井参数与脆性指数关联度大小,对参数进行两两对比并赋值,根据 9点标度法,即范围1,9中的任何整数值和相关倒数值,对各参数分配的值表示每个元素相对于其他元素的偏好强度。元素取值如表 4 所示。将5 项敏感参数两两比较后获得判断矩阵(表 5)。表 4 对比矩阵标度及含义Table 4 Comparison matrix scale and its meaning标度含义1表示两个因素相比,具有相同的重要性3表示两个因素相比,一个因素比另一个因素稍微重要5表示两个因素相比,一个因素比另一个因素明显重要7表示两个因素相比,一个因素比另一个因素强烈重要9表示两个因素相比,一个因素比另一个因素极端重要2,4,6,8上述两相邻判断重要性的中值倒数因素 i 与因素 j 相比,对应 aij,则因素 j 与因素 i 相比,对应 aji=1/aij表 5 判断矩阵Table 5 Judgment matrix元素GR/APIRt/(m)/(gcm-3)CNL/0.01T/(ss-1)GR1335/35/3Rt1/3135/35/31/31/3135/3CNL3/53/51/313T3/53/53/51/31 具体步骤为:1)建立判断矩阵AW=maxW,(20)式中:A 为权重比矩阵;W 为权重向量;max为最大特征值。2)利用层次分析理论中的和积法确定权重向量 W=0.33,0.22,0.18,0.16,0.11,即各参数的权重系数分别为 0.33、0.22、0.18、0.16、0.11。可得脆性指数预测模型为BI=0.33GRN+0.22Rt,N+0.18N+0.16CNLN+0.11TN。(21)3)利用判断矩阵的特征值对判断矩阵进行一致性检验:IC=max-nn-1,(22)ICR=ICIR,(23)式中:n 为选取的元素项数;IC为一致性指标;IR为随机一致性指标,其取值见表 6,本文取值为 1.12;ICR为一致性比例。当 ICR 0.1 时,认为判断矩阵一致性符合要求,否则,应对判断矩阵作适当修正。该判断矩阵计算的 IC=0.106 3,ICR=0.094 9resistivity(Rt)density()neutron(CNL)sonic interval transit time(T),which had a weight coefficient of 0.33,0.22,0.18,0.16,and 0.11,respectively in the new prediction model.The prediction method proposed in this study was applied to the Lucaogou Formation in the Jimusaer Sag.Compared with that determined through laboratory tests,the brittleness index predicted can reflect the actual brittleness of the formation,exhibiting a high consistency.As shown by the re-sults from well tests,the productivity index of oil was proportional to the brittleness index,and a higher brittleness index was associated a high production capacity after fracturing.Therefore,the new method provides a new approach to brittleness index prediction and guides the parameter selection for the fracturing of reservoirs.Key words:grey correlation;analytic hierarchy process;brittleness index;tight reservoir;Lucaogou Formation(本文编辑:叶佩)359