一种岩屑-岩心-井眼-储层多尺度的页岩储层三维可压裂性评价方法
【技术领域】
[0001] 本发明属于页岩油气开发技术领域，尤其涉及一种岩肩-岩心-井眼-储层多尺 度的页岩储层三维可压裂性评价方法。
【背景技术】
[0002] 随着我国国民经济持续快速增长，能源需求急剧增加，石油天然气资源供需矛盾 突出。我国页岩气资源十分丰富，根据2012年国土资源部油气中心最新研究成果表明，我 国页岩气可采资源量为25万亿方，发展潜力巨大。随着美国Barnett、Eagle Ford和我国 涪陵等典型页岩气田的成功开发，页岩气有望成为重要的接替能源。然而，我国页岩气的开 发普遍存在着单井成本高、整体产量低、压后效果参差不齐等瓶颈问题，再加上国际油价的 大幅下跌，使得我国页岩气资源的开发面临着投资成本难以回收的难题，严重制约页岩气 产业的发展。
[0003] 美国的页岩气储层地貌平坦、埋藏较浅、储层物性均匀，通常通过对比水平井井眼 各位置的可压裂性来优选射孔簇位置和压裂层段，压后效果较好。然而，四川地区地质构造 剧烈、储层地质条件复杂，页岩储层内部物性和非均质性差异较大，现有的脆性指数评价方 法仅属于井眼尺度，无法实现整个储层脆性的预测，若水平井井眼轨迹位于脆性较差的储 层或钻遇弹性模量较高和泊松比较小的脆性隔层，采用现有方法优化的射孔位置和压裂层 段将难以得到体积压裂所要求的裂缝网络，从而使页岩气井的产能达不到预期目标甚至部 分井压后并不见气，导致前期钻井和大规模水力压裂的巨额投资难以得到回报。为了优化 页岩气井长水平段轨迹的设计，使1000_2000m的长水平段均位于储层物性和可压裂性较 好的位置，必须建立整个储层的三维可压裂性模型。
【发明内容】
[0004] 本发明的目的在于克服现有技术的缺点，提供一种岩肩-岩心-井眼-储层多尺 度的页岩储层三维可压裂性评价方法，建立页岩储层的三维可压裂性模型。结合地质甜点、 地应力和天然裂缝发育情况，此模型可用于指导页岩气水平井的井眼轨迹设计，为射孔簇 位置和压裂层段的优化提供理论支持，以形成最大化的裂缝网络，提高最终压裂效果，尽快 回收页岩气井的投资成本。
[0005] 本发明的目的通过以下技术方案来实现，一种岩肩-岩心-井眼-储层多尺度的 页岩储层三维可压裂性评价方法，包括下列步骤：
[0006] S1、采用矿物含量脆性指数评价经验公式，计算页岩取心位置处的矿物含量脆性 指数；
[0007] S2、建立内摩擦角、I型和II型裂缝断裂韧性与岩石力学特征参数之间的关系式；
[0008] S3、建立综合考虑矿物含量、弹性参数、内摩擦角、临界应变能释放率和断裂韧性 的页岩可压裂性评价模型，计算钻井位置处储层的可压裂性纵向曲线；
[0009] S4、运用支持向量机算法，得到储层可压裂性与弹性参数之间的聚类分析关系，并 采用单井可压裂性纵向曲线验证可压裂性与弹性参数之间的聚类分析模式；
[0010] S5、应用聚类分析模式到储层三维弹性参数数据体，建立基于岩肩-岩心-井 眼-储层多尺度的页岩储层三维可压裂性模型。
[0011] 进一步地，所述步骤S3包括以下子步骤：
[0012] S31、采用弹性参数脆性指数评价经验公式，计算取心位置处的弹性参数脆性指数 EEn;
[0013] S32、建立综合考虑所述步骤Sl中的矿物含量脆性指数Bn、弹性参数脆性指数EEn和所述步骤S2中的内摩擦角0、临界应变能释放率Ge和断裂韧性的页岩可压裂性评价模 型：
[0014] F/, = u'x B1^1 -(I - ir)X G( " X Since "
[0015] FL = ir X 5,.,, -l· (I -- \ν) χ Κ" " χ K111 " χ 5in^ "
[0016] 式中，FIp FI2分别为考虑页岩应变能释放率和断裂韧性的可压裂性评价指数，无 量纲；
[0017] w为数值范围0~1的权重系数，无量纲；
[0018]
为综合矿物含量脆性指数Bn和弹性参数脆性指数EE n的脆性表达 式，Bn n为综合矿物含量和弹性参数的页岩脆性指数，无量纲；B n、EEn均无量纲；
[0019]
为临界应变能释放率权重表达式，Ge n为砂岩储层的应变能 释放率权重，无量纲；Ge、Ge__、Gejllin分别为砂岩储层的临界应变能释放率、砂岩储层的最大 临界应变能释放率和最小临界应变能释放率，单位为N/m ;
[0020]
为内摩擦角权重表达式，n为砂岩储层的内摩 擦角权重，无量纲；吵、&'鮮^min分别为砂岩储层计算位置处内摩擦角的正弦 值、砂岩储层的最大内摩擦角的正弦值和最小内摩擦角的正弦值，无量纲；
[0021]
为I型断裂韧性权重表达式，Kie n为砂岩储层的I型断裂 韧性权重，无量纲；Κιε、Κκ__、Kiejllin分别为砂岩储层计算位置处的I型断裂韧性、砂岩储层 的最大I型断裂韧性和最小I型断裂韧性，单位为MPa · m1/2;
[0022]
为II型断裂韧性权重表达式，KIK n为砂岩储层的II 型断裂韧性权重，无量纲；Kik、Kik__、Klljllin分别为砂岩储层计算位置处的II型断裂韧性、 砂岩储层的最大Π 型断裂韧性和最小II型断裂韧性，单位为MPa · m1/2;
[0023] S33、利用所建立的综合多因素可压裂性评价模型，结合页岩气井的测井数据和所 述步骤S2中的内摩擦角、I型和II型断裂韧性与岩石力学特征参数之间的关系式，计算页 岩气钻井位置储层的可压裂性纵向曲线特征。
[0024] 进一步地，所述步骤S4包括：
[0025] 采用支持向量机算法，考虑到地震反演储层弹性参数的可靠性，选择纵波阻抗，横 波阻抗，泊松比，拉梅参数与剪切模量之比作为聚类分析训练的弹性参数，训练页岩储层可 压裂性与页岩弹性参数之间的分类关系，找出储层弹性参数与可压裂性之间的聚类分析模 式，并采用所述步骤S3中单井可压裂性纵向曲线验证可压裂性与弹性参数之间的聚类分 析模式。
[0026] 进一步地，所述步骤S5包括：
[0027] 对叠前地震数据体进行预处理以确保数据体质量并得到三维储层弹性参数数据 体，利用所述步骤S3所建立的单井可压裂性纵向曲线，通过应用所述步骤S4中的聚类分析 模式到三维储层弹性参数数据体，建立基于岩肩-岩心-井眼-储层多尺度的页岩储层三 维可压裂性模型。
[0028] 本发明具有以下优点：
[0029] 1、建立了页岩气储层的三维可压裂性评价模型，此模型可真实准确地量化储层不 同空间位置的可压裂性，结合地质甜点，使页岩气水平井总是钻遇高产、高可压裂性层位， 在一定程度上弥补地质条件认识不清楚的不足，提高页岩储层体积压裂的效果，最大化页 岩气井的产能，减小页岩气井投资成本的回收时间；
[0030] 2、运用所建立的模型可评价单井钻井位置处的可压裂性，优选射孔簇位置和压裂 层段；
[0031] 3、运用所建立的模型可使井眼轨迹总是位于可压裂性高的储层内部，钻头在高可 压裂性储层钻进时钻速得到提高，缩短钻井周期，节约钻井成本。
【附图说明】
[0032] 图1是本发明一种岩肩-岩心-井眼-储层多尺度的页岩储层三维可压裂性评价 方法的流程图；
[0033] 图2是本发明中井眼钻井位置处储层的矿物含量脆性指数曲线图；
[0034] 图3是本发明中井眼钻井位置处储层的弹性参数脆性指数曲线图；
[0035] 图4是本发明中综合考虑矿物含量、弹性参数、内摩擦角、临界应变能释放率的钻 井位置处储层的可压裂性纵向曲线图；
[0036] 图5是本发明中基于储层弹性参数与可压裂性之间的聚类分析模式的纵向剖面 图；
[0037] 图6是本发明中所建立的页岩储层三维可压裂性垂向截面图；
[0038] 图7是本发明中所建立的页岩储层三维可压裂性横向截面图。
【具体实施方式】
[0039] 下面结合附图对本发明做进一步的描述，但本发明的保护范围不局限于以下所 述。
[0040] 一种岩肩-岩心-井眼-储层多尺度的页岩储层三维可压裂性评价方法，包括下 列步骤：
[0041] S1、采用X-射线衍射仪等测试设备，开展目标区块钻井岩肩和取样岩心碎片的 矿物组分测试，采用Wang和Gale的矿物含量脆性指数Bn评价经验公式（Wang, F. P.，and J. F. ff. Gale. Screening criteria for shale-gas systems:Gulf Coast Association of Geological Societies Transactions,v. 59,p. 779-793, 2009)：
[0042] Bn= (W石英、方解石+W白云岩）/W总质量 ⑴计算钻井位置或取心位置处页岩的矿 物含量脆性指数Bn，无量纲，式中，为石英和方解石的质量，Wteg为白云石质量，单 位为Kg 为矿物的总质量，单位为Kg;结合校正后的矿物含量曲线，计算井眼钻井位 置处储层的矿物含量脆性指数曲线见图2所示；
[0043] S2、采用高温高压三轴岩石力学测试系统，测试页岩岩心的单轴和三轴岩石力学 强度参数，包括弹性模量E、泊松比V、内摩擦角t等，采用巴西圆盘实验测定页岩试样的抗 拉强度、I型和II型断裂韧性MP κ IIC，建立内摩擦角f、I型和II型断裂韧性与岩石力 学特征参数之间的关系式：
[0044]
(2)
[0045] 式中，σ 致密砂岩的抗拉强度，单位为MPa ; σ n为裂缝面法向围压，单位为MPa ; P