CN109577970B - 一种页岩储层i型断裂韧性的测井评价方法 - Google Patents

一种页岩储层i型断裂韧性的测井评价方法 Download PDF

Info

Publication number
CN109577970B
CN109577970B CN201811521748.0A CN201811521748A CN109577970B CN 109577970 B CN109577970 B CN 109577970B CN 201811521748 A CN201811521748 A CN 201811521748A CN 109577970 B CN109577970 B CN 109577970B
Authority
CN
China
Prior art keywords
sample
fracture toughness
ccnbd
time difference
shale
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811521748.0A
Other languages
English (en)
Other versions
CN109577970A (zh
Inventor
熊健
刘向君
梁利喜
石长有
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN201811521748.0A priority Critical patent/CN109577970B/zh
Publication of CN109577970A publication Critical patent/CN109577970A/zh
Application granted granted Critical
Publication of CN109577970B publication Critical patent/CN109577970B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种页岩储层Ⅰ型断裂韧性的测井评价方法，包括以下步骤：制备CCNBD试样；测量CCNBD试样的纵波速度和横波速度以获得试样的纵波时差和横波时差并计算试样的动态杨氏模量；将上述纵波时差、横波时差、动态杨氏模量输入页岩断裂韧性值的多元回归模型中得到页岩储层断裂韧性测井单井剖面。采用该方法，所需岩样数量少，数据测试量少，岩储层Ⅰ型断裂韧性的测井评价整体的耗时短且准确度高。

Description

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

Claims (5)

1.一种页岩储层I型断裂韧性的测井评价方法，其特征在于，包括以下步骤：
A、制备CCNBD试样，测量该CCNBD试样的几何参数，并计算该试样的体积密度ρb；
B、测量CCNBD试样的纵波速度和横波速度以获得试样的纵波时差Δtc和横波时差Δts并计算该试样的动态杨氏模量Ed；
C、将上述纵波时差Δtc、横波时差Δts、动态杨氏模量Ed输入页岩断裂韧性值的多元回归模型KIC中得到页岩储层断裂韧性测井单井剖面；
D、根据多元回归模型KIC的输出对I型断裂韧性进行评价；
所述动态杨氏模量Ed为：
Figure FDA0003386658630000011
其中，Ed为试样的动态杨氏模量，单位为GPa；Δtc、Δts分别为试样的纵波时差、横波时差，单位为us/m；ρb为试样的体积密度，g/cm3；ac为单位转换系数。
2.根据权利要求1所述的一种页岩储层I型断裂韧性的测井评价方法，其特征在于，所述多元回归模型KIC为：
KIC＝aΔtc+bΔts+cEd+d，其中，-0.002＜a＜0，-0.004＜b＜-0.002，-0.01＜c＜0.02，2＜d＜3。
3.根据权利要求1所述的一种页岩储层I型断裂韧性的测井评价方法，其特征在于，所述多元回归模型KIC为：
KIC＝-0.00112Δtc-0.00333Δts-0.00887Ed+2.8101。
4.根据权利要求1所述的一种页岩储层I型断裂韧性的测井评价方法，其特征在于，步骤A具体为：
制取岩样；
采样人字型切槽巴西圆盘对岩样进行加工，将岩样制备成直径为75mm、直径与厚度比例为5∶2的试样；
计算试样的体积密度。
5.根据权利要求1所述的一种页岩储层I型断裂韧性的测井评价方法，其特征在于，所述多元回归模型KIC的确认方法为：
按照步骤A的方法制备多组CCNBD试样；
挑选几何尺寸满足式(1)、式(2)的CCNBD试样并计算CCNBD试样的体积密度ρb，其中，式(1)为
Figure FDA0003386658630000021
式(2)为
Figure FDA0003386658630000022
其中，R为CCNBD试样半径，单位为mm；D为CCNBD试样直径，单位为mm；α0为无量纲初始裂纹长度，单位为mm；α1为无量纲最大切槽长度，单位为mm；αB为无量纲试样厚度，单位为mm；αs为无量纲切割刀具半径，单位为mm；a0为CCNBD试样初始裂纹长度，，单位为mm；a1为CCNBD试样最大切槽长度，单位为mm；B为CCNBD试样厚度，单位为mm；Ds为切割刀具直径，单位为mm；
获取试样的纵波时差、横波时差、动态杨氏模量；
对CCNBD试样进行加载试验，获取试样的载荷-位移曲线，并从该曲线中获取最大破坏载荷Pmax；
计算标准断裂韧性KIC1，
Figure FDA0003386658630000031
其中，标准断裂韧性KIC1的单位为MPa·m0.5；最大破坏载荷Pmax的单位为kN；D为试样直径，单位为cm；B为试样厚度，单位为cm；
Figure FDA0003386658630000032
为试样的无量纲临界应力强度因子，仅由岩样的几何参数α0、α1和αB决定；
根据标准断裂韧性KIC1的响应规律特征，建立页岩断裂韧性值与页岩岩石纵波时差、横波时差、动态杨氏模量多参数之间的多元回归模型KIC。
CN201811521748.0A 2018-12-11 2018-12-11 一种页岩储层i型断裂韧性的测井评价方法 Active CN109577970B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811521748.0A CN109577970B (zh) 2018-12-11 2018-12-11 一种页岩储层i型断裂韧性的测井评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811521748.0A CN109577970B (zh) 2018-12-11 2018-12-11 一种页岩储层i型断裂韧性的测井评价方法

Publications (2)

Publication Number Publication Date
CN109577970A CN109577970A (zh) 2019-04-05
CN109577970B true CN109577970B (zh) 2022-03-01

Family

ID=65928351

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811521748.0A Active CN109577970B (zh) 2018-12-11 2018-12-11 一种页岩储层i型断裂韧性的测井评价方法

Country Status (1)

Country Link
CN (1) CN109577970B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111898785A (zh) * 2019-05-05 2020-11-06 中国石油天然气集团有限公司 基于页岩的断裂韧性空间分布特征预测方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103256046A (zh) * 2013-04-28 2013-08-21 北京大学 非常规油气藏水平井全缝长压裂参数模拟的方法及装置
CN105277671A (zh) * 2014-06-18 2016-01-27 中国石油化工股份有限公司 一种确定页岩地层脆性指数的方法
CN108019205A (zh) * 2017-09-14 2018-05-11 中国石油天然气股份有限公司 一种确定储层可压裂性指数的方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103256046A (zh) * 2013-04-28 2013-08-21 北京大学 非常规油气藏水平井全缝长压裂参数模拟的方法及装置
CN105277671A (zh) * 2014-06-18 2016-01-27 中国石油化工股份有限公司 一种确定页岩地层脆性指数的方法
CN108019205A (zh) * 2017-09-14 2018-05-11 中国石油天然气股份有限公司 一种确定储层可压裂性指数的方法及装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
SUGGESTED METHOD FOR DETERMINING MODE I FRACTURETOUGHNESS USING CRACKED CHEVRON NOTCHED BRAZILIAN DISC(CCNBD) SPECIMENS;R. J. Fowell;《Int. J. Rock Mech. Min. Sci. & Geomech》;19950615;全文 *
利用测井资料预测深部地层岩石断裂韧性;金衍等;《岩石力学与工程学报》;20010730;第2-3节，图1-2 *
页岩储层I 型和II 型断裂韧性评价方法研究;陈建国等;《岩石力学与工程学报》;20150615;第2-4节，图1-4 *

Also Published As

Publication number Publication date
CN109577970A (zh) 2019-04-05

Similar Documents

Publication Publication Date Title
Landon et al. Nondestructive sample quality assessment of a soft clay using shear wave velocity
Diederichs et al. Measurement of spalling parameters from laboratory testing
Khaksar et al. Rock strength from core and logs, where we stand and ways to go
Rajesh Kumar et al. Prediction of uniaxial compressive strength, tensile strength and porosity of sedimentary rocks using sound level produced during rotary drilling
Mohammed Vipulanandan models to predict the mechanical properties, fracture toughness, pulse velocity and ultimate shear strength of shale rocks
Kolawole et al. Evaluation of geomechanical properties via scratch tests: Where are we and where do we go from here?
He et al. A comprehensive approach for fracability evaluation in naturally fractured sandstone reservoirs based on analytical hierarchy process method
Shi et al. A brittleness index evaluation method for weak-brittle rock by acoustic emission technique
CN110219644B (zh) 确定储层可压裂性指数值空间分布的方法
Tan et al. Brittleness characteristics of tight oil siltstones
CN112504838B (zh) 一种搭载于tbm的岩石力学综合试验及信息评价系统
Wen et al. Correlation of geo-mechanics parameters with uniaxial compressive strength and P-wave velocity on dolomitic limestone using a statistical method
Chen et al. Empirical relations between ultrasonic P-wave velocity, porosity and uniaxial compressive strength
CN112145165B (zh) 一种微裂缝-孔隙型储层动静态渗透率转换方法
Vinck et al. Advanced in situ and laboratory characterisation of the ALPACA chalk research site
Tang et al. An objective crack initiation stress identification method for brittle rock under compression using a reference line
Abugharara et al. Experimental investigation of the effect of shale anisotropy orientation on the main drilling parameters influencing oriented drilling performance in shale
CN109577970B (zh) 一种页岩储层i型断裂韧性的测井评价方法
CN110320571A (zh) 一种致密砂岩储层岩石脆性测井评价方法
Lou et al. Cracking mode analysis of crack initiation in rocks under uniaxial compression
CN116717231A (zh) 一种页岩气水平井岩石可钻性预测方法
Khaksar et al. Enhanced Rock Strength Modelling, Combining Triaxial Compressive Tests, Non-Destructive Index Testing and Well Logs
Bai Experimental research on rock drillability in the center of junggar basin
Aşcı et al. Correlation of physical and mechanical properties with ultrasonic pulse velocities of sandstones in Çenedağ, Kocaeli-Turkey
Wang et al. Use of acoustic emission to determine the effects of bedding and stress paths on micro-cracking evolution of anisotropic shale under cyclic loading tests

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant