技术领域
[0001] 本
发明石油
天然气勘探技术领域,具体涉及一种页岩储层I型断裂韧性的测井评价方法。
背景技术
[0002] 随着我国
能源结构的调整和环境问题的日益突出,作为清洁能源的天然气越来越多的受到重视,特别是非常规油气资源,其中页岩气已成为全球非常规油气资源勘探与开发热点。2013年,美国能源信息署(EIA)发布了包含美国在内的42个国家的页岩气资源评价成果报告,指出全球页岩气的技术可采资源量为220.73×1012m3,其中中国页岩气的技术可12 3
采资源量为31.57×10 m ,显示出我国页岩气资源开发潜
力巨大的特点。压裂改造技术是页岩气资源的高效开发的关键技术之一,而页岩储
层压裂改造效果的好坏与储层的
岩石力学特征、储层性质等有较大关系。页岩储层进行压裂改造时,压裂过程中裂缝起裂和延伸是压裂的核心问题,页岩的脆性是评价压裂效果重要指标。针对页岩脆性评价方面,大量的研究学者已从室内评价、测井预测及其工程应用等方面取得了大量的认识成果,这些研究成果在页岩储层可压性评价方面起到重要作用。此外,压裂过程中裂缝起裂和延伸也受页岩断裂韧性的影响,当裂缝尖端
应力强度因子超过页岩断裂韧性值时,裂缝开始扩展。这说明了页岩岩石断裂韧性也可作为页岩储层可压性的评价参数,然而基于测井信息的页岩岩石断裂韧性预测方法相对较少。
[0003] 目前,关于岩石断裂韧性的测试方法,比如国际岩石力学协会推荐的岩石断裂韧性的测试方法,即1995年提出的人字型切槽巴西圆盘CCNBD,其中CCNBD方法的岩石试样体积小,且能承受较大的
载荷,便于加载,容易实现I型、II型及复合型加载。因此,人字型切槽巴西圆盘法是一种研究岩石断裂韧性的有效方法。但是,存在以下问题,首先,岩样的规格大致需求为直径为75mm、直径与厚度比例为5∶2,在实际
采样中,很少有岩样满足上述直径需求,其需要较多岩样,制备岩样耗时长;其次,在实验测试过程中,需要参数多,耗时长。
发明内容
[0004] 本发明为了解决上述技术问题提供一种页岩储层I型断裂韧性的测井评价方法。
[0005] 本发明通过下述技术方案实现:
[0006] 一种页岩储层I型断裂韧性的测井评价方法,包括以下步骤:
[0007] A、制备CCNBD试样,测量该CCNBD试样的几何参数,并计算该试样的体积
密度ρb;
[0008] B、测量CCNBD试样的纵
波速度和横波速度以获得试样的纵波时差Δtc和横波时差Δts并计算试样的动态
杨氏模量Ed;
[0009] C、将上述纵波时差Δtc、横波时差Δts、动态杨氏模量Ed输入页岩断裂韧性值的多元回归模型KIC中得到页岩储层断裂韧性测井单井剖面;
[0010] D、根据多元回归模型KIC的输出对I型断裂韧性进行评价。
[0011] 作为优选,所述多元回归模型KIC为:
[0012] KIC=aΔtc+bΔts+cEd+d,其中,-0.002<a<0,-0.004<b<-0.002,-0.01<c<0.02,2<d<3。
[0013] 作为优选,所述多元回归模型KIC为:
[0014] KIC=-0.00112Δtc-0.00333Δts-0.00887Ed+2.8101。
[0015] 本方案的页岩断裂韧性值仅与纵波时差Δtc、横波时差Δts、动态杨氏模量Ed相关,
数据采集时,仅需采集一
块岩样的上述数据,即可得出页岩断裂韧性值KIC,其所需岩样数量少,采集岩样耗时短;采集数据少,岩储层I型断裂韧性的测井评价整体的耗时短,评价准确度高,满足实际测试要求。
[0016] 本发明与
现有技术相比,具有如下的优点和有益效果:
[0017] 1、本发明的评价方法整体仅需采集一块岩样,所需岩样数量少,采集岩样耗时短;仅需测试岩样的
体积密度ρb、纵波时差Δtc、横波时差Δts的3个参数,数据测试量少,岩储层I型断裂韧性的测井评价整体的耗时短且准确度高。
附图说明
[0018] 此处所说明的附图用来提供对本发明
实施例的进一步理解,构成本
申请的一部分,并不构成对本发明实施例的限定。
[0019] 图1为CCNBD试样加工图。
[0020] 图2为图1中CCNBD试样的A-A剖视图。
[0021] 图3为17组试样的纵波时差、横波时差图。
[0022] 图4为17组试样的动态杨氏模量图。
[0023] 图5为页岩断裂韧性与弹性参数间的关系图。
[0024] 图6为页岩断裂韧性与纵波时差的关系图。
[0025] 图7为页岩断裂韧性与横波时差的关系图。
[0026] 图8为某地区的龙
马溪组页岩断裂韧性A井的评价图。
[0028] 图10为预测的页岩断裂韧性值与实测值的对比图。
具体实施方式
[0029] 为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
[0030] 实施例1
[0031] 如图9所示的一种页岩储层I型断裂韧性的测井评价方法,包括以下步骤:
[0032] A、制备CCNBD试样,测量该CCNBD试样的几何参数,并计算该试样的体积密度ρb;
[0033] B、测量CCNBD试样的纵波速度和横波速度以获得试样的纵波时差Δtc和横波时差Δts并计算试样的动态杨氏模量Ed;
[0034] C、将上述纵波时差Δtc、横波时差Δts、动态杨氏模量Ed输入页岩断裂韧性值的多元回归模型KIC中得到页岩储层断裂韧性测井单井剖面,
[0035] D、根据多元回归模型KIC的输出对I型断裂韧性进行评价。
[0036] 实施例2
[0037] 基于上述实施例的原理,本实施例公开一具体方法。
[0038] 确定多元回归模型KIC,其具体方法如下:
[0039] 从矿场取回页岩井下或露头岩样,对岩样进行岩性描述,在此
基础上,制备多组CCNBD岩样,本方案以17组岩样为例说明。
[0040] 岩石断裂韧性测试的实验样品按照国际岩石力学协会ISRM推荐的人字型切槽巴西圆盘CCNBD对岩样进行加工,将岩样制备为直径为75mm,直径与厚度比例为5∶2的试样,如图1、图2所示。
[0041] 挑选几何尺寸满足式(1)、式(2)的CCNBD试样并计算CCNBD试样的体积密度ρb,[0042] 其中,式(1)为
[0043] 式(2)为
[0044] 其中,R为CCNBD试样半径,单位为mm;D为CCNBD试样直径,单位为mm;α0为无量纲初始裂纹长度,单位为mm;α1为无量纲最大切槽长度,单位为mm;αB为无量纲试样厚度,单位为mm;.αs为无量纲切割刀具半径,单位为mm;α0为CCNBD试样初始裂纹长度,,单位为mm;α1为CCNBD试样最大切槽长度,单位为mm;B为CCNBD试样厚度,单位为mm;Ds为切割刀具直径,单位为mm。
[0045] 采用
超声波透射法对CCNBD试样进行纵波速度和横波速度进行测量,获取试样的纵波时差、横波时差,测量结果如图3所示。
[0046] 计算试样的弹性参数,包括动态杨氏模量Ed,
[0047]
[0048] 其中,为试样的动态杨氏模量,单位为GPa;、分别为试样的纵波时差、横波时差,单位为us/m;为试样的体积密度,g/cm3;为单位转换系数。17组岩样的计算结果如图4所示。
[0049] 采用MTS伺服刚性控制试验机对CCNBD试样进行加载试验,获取试样的载荷-位移曲线,并从该曲线中获取最大破坏载荷Pmax;
[0050] 结合试样的几何参数,计算标准断裂韧性,,
[0051] 其中,标准断裂韧性的单位为MPa.m0.5;最大破坏载荷Pmax的单位为kN;D为试样直径,单位为cm;B为试样厚度,单位为cm;为试样的无量纲临界应力强度因子,仅由岩样的几何参数、和决定;
[0052] 基于声波测试和断裂韧性测试结果,获取了页岩断裂韧性对声波属性参数的敏感性,如图5、6、7所示,其相关系数计算公式为
[0053]
[0054] 式中,r为相关系数,xi为第i个页岩样品断裂韧性响应参数实测值,为页岩断裂韧性响应参数所有实测值的平均值,yi为第i个页岩样品断裂韧性实测值,为页岩断裂韧性所有实测值的平均值。
[0055] 根据标准断裂韧性KIC1的响应规律特征,建立页岩断裂韧性值与页岩岩石纵波时差、横波时差、动态杨氏模量等多参数之间的多元回归模型KIC。
[0056] KIC=-0.00112Δtc-0.00333Δts-0.00887Ed+2.8101,其判定系数R2为0.8183。判定系数R2计算公式为
[0057]
[0058] 式中,yi为第i个页岩样品断裂韧性实测值,为页岩样品断裂韧性所有实测值的平均值, 为拟合公式的计算值。
[0059] 多次实验表明,多元回归模型KIC满足下式条件时均可实现对I型断裂韧性的可靠预测。KIC=aΔtc+bΔts+cEd+d,其中,-0.002<a<0,-0.004<b<-0.002,-0.01<c<0.02,2<d<3。
[0060] 上述多元回归模型KIC确定后,后期应用时:
[0061] 从矿场取回页岩井下或露头岩样,对岩样进行岩性描述,在此基础上,制备1组CCNBD岩样。
[0062] 岩石断裂韧性测试的实验样品按照国际岩石力学协会ISRM推荐的人字型切槽巴西圆盘CCNBD对岩样进行加工,将岩样制备为直径为75mm,直径与厚度比例为5∶2的试样。
[0063] 挑选几何尺寸满足式(1)、式(2)的CCNBD试样并计算CCNBD试样的体积密度ρb。
[0064] 采用
超声波透射法对CCNBD试样进行纵波速度和横波速度进行测量,获取试样的纵波时差、横波时差。
[0065] 计算试样的弹性参数,包括动态杨氏模量Ed,
[0066] 将上述纵波时差Δtc、横波时差Δts、动态杨氏模量Ed输入页岩断裂韧性值的多元回归模型KIC=-0.00112Δtc一0.00333Δts-0.00887Ed+2.8101中得到页岩储层断裂韧性的测井单井剖面。
[0067] 为了证明本方案方法的可行性,对其准确性进行了验证。具体的以13组实验数据予以证明,采用实测断裂韧性值与本方案的输出进行对比。引入相对误差E对其进行评价。
[0068]
[0069] 式中,Vexp,i为实验中第i组样品的实测断裂韧性值,单位为MPa.m1/2;Vcal,i为采用本方案输出的第i组样品的断裂韧性值,MPa.m1/2。
[0070] 结果表明本方案的多元回归模型KIC预测的页岩断裂韧性值和实测值较吻合,如图10所示,预测值的平均相对误差为8.82%,符合实际测试的使用要求。
[0071] 根据多元回归模型计算得到页岩储层断裂韧性测井单井剖面,如图8所述。页岩层段的断裂韧性值剖面随深度的变化趋势与脆性指数的关系存在相反的趋势,即脆性指数增加的井段,而断裂韧性呈减小的趋势。由于岩石脆性指数越高,岩石断裂韧性越低,其储层越易进行压裂改造。这说明了断裂韧性和脆性指数一样也可用于页岩储层可压性的评价,同时也说明了构建了页岩断裂韧性预测评价方法的可靠性。
[0072] 综上所述,本发明通过提供一种将断裂韧性测试、超声波测试等岩石物理测量结果和测井数据结合起来预测页岩储层I型断裂韧性的方法,从而实现了提供一种基于岩石物理测量结果和测井数据合理有效地预测页岩储层I型断裂韧性值的目的,通过本发明所提供的预测方法能够较为真实地描述出页岩储层断裂韧性特性。
[0073] 以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何
修改、等同替换、改进等,均应包含在本发明的保护范围之内。