CN106644736A - 一种测量储层岩石的脆性指数的方法及装置 - Google Patents
一种测量储层岩石的脆性指数的方法及装置 Download PDFInfo
- Publication number
- CN106644736A CN106644736A CN201611192098.0A CN201611192098A CN106644736A CN 106644736 A CN106644736 A CN 106644736A CN 201611192098 A CN201611192098 A CN 201611192098A CN 106644736 A CN106644736 A CN 106644736A
- Authority
- CN
- China
- Prior art keywords
- brittleness index
- measured
- peak
- stress
- value
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
- G01N3/10—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces generated by pneumatic or hydraulic pressure
- G01N3/12—Pressure testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0019—Compressive
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0076—Hardness, compressibility or resistance to crushing
- G01N2203/0087—Resistance to crushing
Abstract
本申请公开了一种测量储层岩石的脆性指数的方法及装置。该方法包括以下步骤:对待测岩心进行破裂实验并获取待测岩心的应力值与应变值;基于获取的应力值和应变值生成应力应变关系曲线;基于应力应变关系曲线通过第一算法计算峰前脆性指数;基于应力应变关系曲线通过第二算法计算峰后脆性指数;基于峰前脆性指数以及峰后脆性指数通过第三算法计算储层岩石的脆性指数。该方法是一种综合的脆性指数判别方法,克服了之前参数多、标准不统一的难题,为致密储层岩石力学性质评价奠定坚实的基础。
Description
技术领域
本发明涉及石油勘探开发中的岩石力学评价领域,尤其涉及一种测量储层岩石的脆性指数的方法及装置。
背景技术
致密气、页岩气、煤层气、致密油等致密油气储存在致密储层中。随着开采技术的进步,这些致密油气具有较高的开发潜力,但由于致密储层孔隙度和渗透率都非常低,一般情况下无自然产能或者自然产能较低,需要进行规模压裂才能产出工业油气流。储层的脆性是寻找相对优质储层、遴选射孔改造层段以及设计压裂规模的重要基础参数。
现有的储层脆性的测量方法有很多,例如,1、基于地球物理资料的动态弹性参数法(参见2008年《SPE》,Rickman R和Mullen M J等著作的《A practical use of shalepetrophysics for stimulation design optimization all shale plays are notclones of the Barnett Shale》);2、基于储层的岩心分析或测井评价的岩性参数指示法,然而,上述的两类方法均存在多解性。
岩心的脆性测量方法还包括有基于破裂实验的应力应变关系法,该方法可以克服上述的多解性的问题。前人对基于破裂实验的应力应变关系法进行了大量的研究,主要有以下几种观点:1、H.Honda和Y.Sanada提出以硬度和坚固性差异表征脆性(参见杂志《Fuel》中文章《Hardness of coal》);2、V.Hucka和B.Das建议采用试样抗压强度和抗拉强度的差异表示脆性(参见杂志《International Journal of Rock Mechanics and MiningSciences and Geomechanics Abstracts》中文章《Brittleness determination of rocksby different methods》);3、A.W.Bishop认为应从标准试样的应变破坏试验入手,分析应力释放的速度进而表征脆性(参见论文集《Proceedings of the geotechnicalconference》中文章《Progressive failure with special reference to the mechanismcausing it》),但这些方法大多针对具体问题提出没有较高的通用性,且仅仅考虑了破裂后的应力应变的关系,没有考虑到破裂前的应力应变关系。
李庆辉(参见杂志《岩石力学与工程学报》中文章《页岩脆性的室内评价方法及改进》)提出基于全应力-应变的脆性评价方法。该方法考虑了破裂前和破裂后的两种应力应变关系曲线,但该脆性评价方法需要确定的参数较多,且没有统一的标准。因此如何使用较少的参数准确地测量储层脆性是当前亟待解决的问题。
发明内容
针对上述问题,本发明的目的在于提供一种测量储层岩石的脆性指数的方法及装置,该方法可以在获取较少的参数情况下测量储层的脆性指数。
为了实现上述目的,本发明提供的一种测量储层岩石的脆性指数的方法。该方法包括:
对待测岩心进行破裂实验并获取所述待测岩心的应力值与应变值;
基于获取的所述应力值和所述应变值生成应力应变关系曲线;
基于所述应力应变关系曲线通过第一算法计算峰前脆性指数;
基于所述应力应变关系曲线通过第二算法计算峰后脆性指数;
基于所述峰前脆性指数以及所述峰后脆性指数通过第三算法计算所述待测岩心的脆性指数。
优选的,所述对待测岩心进行破裂实验并获取待测岩心的应力值与应变值的步骤包括:
对所述待测岩心施加预定的围限压力;
对所述待测岩心施加预设的孔隙压力;
对所述待测岩心施加预设的所述围限压力以及预设的所述孔隙压力后,对所述待测岩心施加轴向压力并逐渐增大,同时获取所述待测岩心的所述应力值与所述应变值,直至待测岩心破裂。
优选的,所述围限压力的值根据待测岩心所在的深度和地层密度确定;所述孔隙压力的值根据所述待测岩心所在的深度和流体密度确定。
优选的,所述基于应力应变关系曲线通过第一算法计算峰前脆性指数的步骤包括:
获取所述应力应变关系曲线上达到最大屈服强度的应变值εy以及峰值强度的应变值εp;
基于获取的最大屈服强度的应变值εy以及峰值强度的应变值εp通过第一算法计算峰前脆性指数Bb。
优选的,所述第一算法为:
其中Bb为峰前脆性指数;εy为最大屈服强度的应变值,%;εp为峰值强度的应变值,%。
优选的,所述基于应力应变关系曲线通过第二算法计算峰后脆性指数的步骤包括:
获取所述应力应变关系曲线上达到峰值强度的应力值τp以及达到残余强度时的应力值τr;
基于获取的达到峰值强度的应力值τp以及达到残余强度时的应力值τr通过第二算法计算峰后脆性指数Ba。
优选的,所述第二算法为:
其中,Ba为峰后脆性指数;τp为达到峰值强度的应力值,单位为兆帕;τr为达到残余强度时的应力值,单位为兆帕。
优选的,所述基于峰前脆性指数以及峰后脆性指数通过第三算法计算储层岩石的脆性指数的步骤包括:
基于计算得到的峰前脆性指数Ba以及峰后脆性指数Bb通过所述第三算法计算待测岩心的脆性指数B;
其中,所述第三算法为:
B=(Bb+Ba)/2
B为储层岩石的脆性指数;Ba为峰后脆性指数;Bb为峰前脆性指数。
优选的,所述待测岩心的直径为2.5~2.6厘米之间;所述待测岩心的高度为4.5~5.5厘米之间。
本发明还提供了一种测量储层岩石的脆性指数的装置,包括:
破裂实验单元,用于对待测岩心进行破裂实验并获取待测岩心的应力值与应变值;
与所述破裂实验单元电连接的处理单元,用于获取所述待测岩心的应力值和应变值并基于获取到的应力值和应变值生成应力应变关系曲线;
与所述处理单元电连接的第一计算单元,用于基于应力应变关系曲线通过第一算法计算峰前脆性指数;
与所述处理单元电连接的第二计算单元,用于基于应力应变关系曲线通过第二算法计算峰后脆性指数;
分别与所述第一计算单元以及第二计算单元电连接的第三计算单元,用于获取峰前脆性指数以及峰后脆性指数并通过第三算法计算储层岩石的脆性指数。
优选的,所述待测岩心通过取芯钻机或钻井设备进行钻取。
优选的,所述破裂实验单元包括:
具有所述待测岩心的放置位的放置组件;
设置在所述放置组件上的围限压力施加组件,用于向所述待测岩心施加围限压力;
设置在所述放置组件上的孔隙压力施加组件,用于向所述待测岩心施加孔隙压力;
设置在所述放置组件上的轴向压力施加组件,用于向所述待测岩心施加轴向压力;
设置在所述放置组件上的数据获取组件,用于获取所述待测岩心的应力值和应变值。
本发明针对当前勘探热点页岩气和致密砂岩等致密储层,通过对待测的致密储层的岩心进行破裂实验分析获取待测岩心的应力应变关系曲线,再根据获取到的应力应变关系曲线分析曲线上峰前和峰后两个阶段并分别计算峰前脆性指数和峰后脆性指数,综合上述两个阶段的脆性指数后,进而得到储层岩石的脆性指数。该方法可以在获取较少的参数情况下测量储层的脆性指数,克服了之前参数多、标准不统一的难题,为致密储层岩石力学性质评价奠定坚实的基础。
附图说明
图1为本发明的实施例提供的测量储层的脆性的方法的流程图;
图2为本发明实施例中典型应力应变关系曲线示意图;
图3为本发明实施例中应力应变关系曲线关键点示意图;
图4a为本发明实施例中样品1的应力应变关系曲线图;
图4b为本发明实施例中样品2的应力应变关系曲线图;
图5为本发明实施例中测量储层脆性的装置中各单元的连接示意图;
图6为本发明实施例中测量储层脆性的装置的破裂实验单元中各部件的示意图。
附图说明:1、破裂实验单元;2、处理单元;3、第一计算单元;4、第二计算单元;5、第三计算单元。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
本发明公开了一种测量储层岩石的脆性指数的方法,参照图1所示,该方法包括:
S101:对待测岩心进行破裂实验并获取待测岩心的应力值与应变值。
在该步骤中,首先可以使用破裂实验的设备对待测岩心施加预定的围限压力。其中,围限压力的值可以根据待测岩心未钻取之前所在的深度和地层密度确定,即待测岩心在未被钻取之前所处在的储层对其的压力值。施加围限压力的目的在于模拟岩心在地下的环境。保证待测岩心在实验室中进行测量时,能保证该待测岩心处于较为真实的环境中,从而确保该方法获取到的数据的准确度。
对待测岩心施加预设的孔隙压力。其中,孔隙压力的值根据待测岩心所在的深度和流体密度确定,即孔隙压力为待测岩心处于储层中时,待测岩心内的孔隙中存在的油或者气对其的压力。施加孔隙压力值的目的在于模拟岩心在地下的环境,进而保证待测岩心在实验室中进行测量时,能保证该待测岩心处于较为真实的环境中,从而保证获取到的数据的准确度。其中,围限压力和孔隙压力的施加顺序没有具体的限制,即,可以先施加围限压力然后再施加孔隙压力,同样,也可以先施加孔隙压力再施加围限压力。
对待测岩心施加预设的围限压力以及预设的孔隙压力后,对待测岩心施加轴向压力并逐渐增大,直至待测岩心破裂。具体的,施加轴向压力的装置可以与待测岩心的端部相接触,进而对其施加一轴向的压力。
对待测岩心施加轴向压力的同时,操作人员可以获取待测岩心的应力值与应变值。具体的,操作人员可以在待测岩心以及向待测岩心施加轴向压力的装置上设置应力感应器以及应变感应器,进而获取在对待测岩心施加轴向压力时的应力值以及应变值。该技术为现有技术中获取应力与应变值的惯用手段,在此不再赘述。
在该步骤中,待测岩心可以通过取芯设备在储层中钻取获得,或者也可以通过钻井设备在钻井的过程中获得。为了保证待测岩心的尺寸能在实验室的设备上进行实验。待测岩心的直径可以为2.5~2.6厘米(cm)之间;待测岩心的高度为4.5~5.5厘米之间。
S102:基于获取的应力值和应变值生成应力应变关系曲线。
在该步骤中,如图2所示,该图为本发明实施例中根据应力值以及应变值生成的典型应力应变关系曲线示意图。该曲线可被划分为弹性区域、塑性区域和破裂区3个部分。
又如图3所示,该图为本发明实施例的应力应变关系曲线关键点示意图。通过该曲线可读出3个关键点坐标,分别是屈服强度的应力和应变(εy,τy),峰值强度的应力和应变(εp,τp)以及残余强度的应力和应变(εr,τr)。
例如,如图4a以及图4b所示,该图为本发明的实施例提供的两块样品的具体的应力应变关系曲线图。
图4a为样品1的应力应变关系曲线。参照图3标示的各个关键点坐标,可以根据该应力应变关系曲线读取屈服强度的应力和应变(εy,τy),峰值强度的应力和应变(εp,τp)以及残余强度的应力和应变(εr,τr)的值,进而可以根据上述的值计算样品1的脆性指数。其中,根据上述的步骤计算出的样品1的脆性指数值见下表1。
同样的,图4b为样品2的应力应变关系曲线。参照图3标示的各个关键点坐标,可以根据该应力应变关系曲线读取屈服强度的应力和应变(εy,τy),峰值强度的应力和应变(εp,τp)以及残余强度的应力和应变(εr,τr)的值,进而可以根据上述的值计算样品2的脆性指数。其中,根据上述的步骤计算出的样品2的脆性指数值见下表1。
S103:基于应力应变关系曲线通过第一算法计算峰前脆性指数。
在该步骤中,参照图3所示,首先获取应力应变关系曲线上达到最大屈服强度的应变值εy以及峰值强度的应变值εp,即,上述的屈服强度的应力和应变(εy,τy),峰值强度的应力和应变(εp,τp)中的横向坐标值。其中,屈服强度的应力是在岩石压缩过程中,当应力达到一定值时,应力有微小的增加而应变却急剧增长的现象,该应力值即为屈服压力。它可表征岩石弹性变形段的大小。
基于获取的最大屈服强度的应变值εy以及峰值强度的应变值εp通过第一算法计算峰前脆性指数Bb。其中,第一算法为:
其中Bb为峰前脆性指数;εy为最大屈服强度的应变值,%;εp为峰值强度的应变值,%。
S104:基于应力应变关系曲线通过第二算法计算峰后脆性指数。
在该步骤中,首先获取应力应变关系曲线上达到峰值强度的应力值τp以及达到残余强度时的应力值τr,即,上述曲线中的峰值强度的应力和应变(εp,τp)以及残余强度的应力和应变(εr,τr)的纵向坐标值。
基于获取的达到峰值强度的应力值τp以及达到残余强度时的应力值τr通过第二算法计算出峰后脆性指数Ba。其中,第二算法为:
其中,Ba为峰后脆性指数;τp为达到峰值强度的应力值,单位为兆帕(Mpa);τr为达到残余强度时的应力值,单位为兆帕(Mpa)。
S105:基于峰前脆性指数以及峰后脆性指数通过第三算法计算储层岩石的脆性指数。
在该步骤中,首先基于计算得到的峰前脆性指数Ba以及峰后脆性指数Bb通过第三算法计算储层岩石的脆性指数B;
其中,第三算法为:
B=(Bb+Ba)/2
其中,B为储层岩石的脆性指数;Ba为峰后脆性指数;Bb为峰前脆性指数。
参照表1所示,该表为19个样品通过该方法计算出的储层岩石的脆性指数B。
表1
编号 | τy | τp | τr | εy | εp | εr | B |
样品1 | 211.0 | 212.4 | 167.0 | 0.410 | 0.413 | 0.474 | 60.32 |
样品2 | 102.0 | 113.9 | 83.3 | 0.316 | 0.551 | 0.628 | 42.11 |
样品3 | 227.6 | 229.8 | 142.9 | 0.394 | 0.398 | 0.485 | 68.41 |
样品4 | 102.8 | 114.0 | 78.1 | 0.398 | 0.568 | 0.628 | 50.78 |
样品5 | 121.4 | 132.4 | 93.6 | 0.387 | 0.615 | 0.682 | 46.12 |
样品6 | 111.7 | 118.0 | 64.5 | 0.489 | 0.618 | 0.676 | 62.23 |
样品7 | 114.5 | 122.3 | 82.9 | 0.400 | 0.528 | 0.667 | 53.99 |
样品8 | 118.1 | 123.8 | 94.5 | 0.416 | 0.567 | 0.660 | 48.52 |
样品9 | 108.0 | 118.2 | 87.1 | 0.347 | 0.529 | 0.617 | 45.95 |
样品10 | 109.5 | 121.7 | 94.1 | 0.360 | 0.551 | 0.606 | 44.01 |
样品11 | 108.3 | 113.5 | 93.3 | 0.454 | 0.532 | 0.580 | 51.57 |
样品12 | 106.9 | 114.6 | 79.6 | 0.359 | 0.533 | 0.647 | 48.95 |
样品13 | 99.7 | 109.4 | 73.2 | 0.321 | 0.527 | 0.621 | 47.00 |
样品14 | 105.2 | 106.3 | 90.9 | 0.476 | 0.490 | 0.533 | 55.82 |
样品15 | 104.5 | 114.2 | 77.9 | 0.379 | 0.537 | 0.628 | 51.18 |
样品16 | 111.7 | 120.8 | 89.0 | 0.428 | 0.550 | 0.597 | 52.07 |
样品17 | 116.1 | 125.5 | 92.5 | 0.441 | 0.566 | 0.615 | 52.11 |
样品18 | 117.5 | 121.5 | 78.2 | 0.461 | 0.547 | 0.616 | 59.96 |
样品19 | 132.0 | 134.6 | 98.1 | 0.505 | 0.581 | 0.640 | 57.02 |
其中,表中第一列为编号,从样品1到样品19;第二列为τy,即达到屈服强度时的应力;第三列为τp,即达到峰值强度时的应力;第四列为τr,即达到残余强度时的应力;第五列为εy,即达到屈服强度时的应变;第六列为εp,即达到峰值强度时的应变;第七列为εr,即达到残余强度时的应变;第八列为B,即储层岩石的脆性指数。
通过表1不难看出:通过本发明提供的测量储层岩石的脆性指数的方法可以准确的测量出待测岩心的脆性指数。
本发明针对当前勘探热点页岩气和致密砂岩等致密储层,通过对待测岩心进行破裂实验分析获取待测岩心的应力应变关系曲线,再根据获取到的应力应变关系曲线分析曲线上峰前和峰后两个阶段并分别计算峰前脆性指数和峰后脆性指数,再综合两个阶段的脆性指数进而得到储层岩石的脆性指数。该方法可以在较少的参数情况下测量储层的脆性指数,克服了之前参数多、标准不统一的难题,为致密储层岩石力学性质评价奠定坚实的基础。
本发明还公开了一种测量储层岩石的脆性指数的装置,参照图5所示,包括:破裂实验单元1,用于对待测岩心进行破裂实验并获取待测岩心的应力值与应变值;与破裂实验单元1电连接的处理单元2,处理单元2用于基于获取的应力值和应变值生成应力应变关系曲线;与处理单元2电连接的第一计算单元3,用于基于应力应变关系曲线通过第一算法计算峰前脆性指数;与处理单元2电连接的第二计算单元4,用于基于应力应变关系曲线通过第二算法计算峰后脆性指数;分别与第一计算单元3以及第二计算单元4电连接的第三计算单元5,用于基于峰前脆性指数以及峰后脆性指数通过第三算法计算储层岩石的脆性指数。
其中,处理单元2、第一计算单元3、第二计算单元4以及第三计算单元5可以为集成芯片、第一计算单元3、第二计算单元4以及第三计算单元5的算法可以写入集成芯片内。
在一具体的实施方式中,参照图6所示,破裂实验单元1可以包括:具有待测岩心的放置位的放置组件;其中,放置组件可以为一能设置在实验室工作台上的金属支撑板。该金属支撑板上可以设置有能与待测岩心的直径相适配的放置位,放置位上可以设置有放置装置,例如,凹槽,夹持件等能将待测岩心固定设置在支撑板上的装置,进而使得向待测岩心施加各个方向的压力时,能保证待测岩心的稳定。
设置在放置组件上的围限压力施加组件,用于向待测岩心施加围限压力;其中,围限压力施加组件可以为液压组件。围限压力施加组件可以围设在待测岩心的周围。围限压力施加组件的围限压力值可以根据待测岩心未钻取之前所在的深度和地层密度确定,即,待测岩心在未被钻取之前所处在的储层对其的压力值。其中,施加围限压力值的目的在于模拟岩心在地下的环境。保证待测岩心在实验室中进行测量时,能保证该待测岩心处于较为真实的环境中,从而保证获取到的数据的准确度。
设置在放置组件上的孔隙压力施加组件,用于向待测岩心施加孔隙压力。其中,孔隙压力的值可以根据待测岩心所在的深度和流体密度确定,即,孔隙压力为待测岩心处于储层中时,待测岩心内的孔隙中存在的油或者气对其的压力。施加孔隙压力值的目的在于模拟岩心在地下的环境,进而保证待测岩心在实验室中进行测量时,能保证该待测岩心处于较为真实的环境中,从而保证获取到的数据的准确度。
孔隙压力组件可以为液压组件。孔隙压力组件的作用位置为待测岩心的空隙部,孔隙压力组件的压力流体可以进入待测岩心的孔隙内,使得该组件能还原待测岩心在地下时,存在的油或者气对其的压力。
设置在放置组件上的轴向压力施加组件,用于向待测岩心施加轴向压力。其中,轴向压力施加组件可以包括一具有能与待测岩心的端部相接触的压力头,当将待测岩心放置在放置组件上时,轴向压力施加组件的压力头可以与待测岩心的上端部相接触,轴向压力施加组件进而可以向待测岩心施加轴向的压力。
设置在放置组件上的数据获取组件,用于获取待测岩心的应力值和应变值。其中,数据获取组件可以包括应力感应器以及应变感应器。应力感应器以及应变感应器可以获取在对待测岩心施加轴向压力时的应力值以及应变值。
应该理解,以上描述是为了进行图示说明而不是为了进行限制。通过阅读上述描述,在所提供的示例之外的许多实施方式和许多应用对本领域技术人员来说都将是显而易见的。因此,本教导的范围不应该参照上述描述来确定,而是应该参照所附权利要求以及这些权利要求所拥有的等价物的全部范围来确定。出于全面之目的,所有文章和参考包括专利申请和公告的公开都通过参考结合在本文中。在前述权利要求中省略这里公开的主题的任何方面并不是为了放弃该主体内容,也不应该认为发明人没有将该主题考虑为所公开的发明主题的一部分。
Claims (12)
1.一种测量储层岩石的脆性指数的方法,其特征在于,包括:
对待测岩心进行破裂实验并获取所述待测岩心的应力值与应变值;
基于获取的所述应力值和所述应变值生成应力应变关系曲线;
基于所述应力应变关系曲线通过第一算法计算峰前脆性指数;
基于所述应力应变关系曲线通过第二算法计算峰后脆性指数;
基于所述峰前脆性指数以及所述峰后脆性指数通过第三算法计算所述待测岩心的脆性指数。
2.根据权利要求1所述的测量储层岩石的脆性指数的方法,其特征在于,所述对待测岩心进行破裂实验并获取待测岩心的应力值与应变值的步骤包括:
对所述待测岩心施加预定的围限压力;
对所述待测岩心施加预设的孔隙压力;
对所述待测岩心施加预设的所述围限压力以及预设的所述孔隙压力后,对所述待测岩心施加轴向压力并逐渐增大,同时获取所述待测岩心的所述应力值与所述应变值,直至待测岩心破裂。
3.根据权利要求2所述的测量储层岩石的脆性指数的方法,其特征在于,所述围限压力的值根据待测岩心所在的深度和地层密度确定;所述孔隙压力的值根据所述待测岩心所在的深度和流体密度确定。
4.根据权利要求1所述的测量储层岩石的脆性指数的方法,其特征在于,所述基于应力应变关系曲线通过第一算法计算峰前脆性指数的步骤包括:
获取所述应力应变关系曲线上达到最大屈服强度的应变值εy以及峰值强度的应变值εp;
基于获取的最大屈服强度的应变值εy以及峰值强度的应变值εp通过第一算法计算峰前脆性指数Bb。
5.根据权利要求4所述的测量储层岩石的脆性指数的方法,其特征在于,所述第一算法为:
其中Bb为峰前脆性指数;εy为最大屈服强度的应变值,%;εp为峰值强度的应变值,%。
6.根据权利要求1所述的测量储层岩石的脆性指数的方法,其特征在于,所述基于应力应变关系曲线通过第二算法计算峰后脆性指数的步骤包括:
获取所述应力应变关系曲线上达到峰值强度的应力值τp以及达到残余强度时的应力值τr;
基于获取的达到峰值强度的应力值τp以及达到残余强度时的应力值τr通过第二算法计算峰后脆性指数Ba。
7.根据权利要求6所述的测量储层岩石的脆性指数的方法,其特征在于,所述第二算法为:
其中,Ba为峰后脆性指数;τp为达到峰值强度的应力值,单位为兆帕;τr为达到残余强度时的应力值,单位为兆帕。
8.根据权利要求1所述的测量储层岩石的脆性指数的方法,其特征在于,所述基于峰前脆性指数以及峰后脆性指数通过第三算法计算储层岩石的脆性指数的步骤包括:
基于计算得到的峰前脆性指数Ba以及峰后脆性指数Bb通过所述第三算法计算待测岩心的脆性指数B;
其中,所述第三算法为:
B=(Bb+Ba)/2
B为储层岩石的脆性指数;Ba为峰后脆性指数;Bb为峰前脆性指数。
9.根据权利要求1所述的测量储层岩石的脆性指数的方法,其特征在于,所述待测岩心的直径为2.5~2.6厘米之间;所述待测岩心的高度为4.5~5.5厘米之间。
10.一种测量储层岩石的脆性指数的装置,其特征在于,包括:
破裂实验单元,用于对待测岩心进行破裂实验并获取待测岩心的应力值与应变值;
与所述破裂实验单元电连接的处理单元,用于获取所述待测岩心的应力值和应变值并基于获取到的应力值和应变值生成应力应变关系曲线;
与所述处理单元电连接的第一计算单元,用于基于应力应变关系曲线通过第一算法计算峰前脆性指数;
与所述处理单元电连接的第二计算单元,用于基于应力应变关系曲线通过第二算法计算峰后脆性指数;
分别与所述第一计算单元以及第二计算单元电连接的第三计算单元,用于获取峰前脆性指数以及峰后脆性指数并通过第三算法计算储层岩石的脆性指数。
11.根据权利要求10所述的测量储层岩石的脆性指数的装置,其特征在于,所述待测岩心通过取芯钻机或钻井设备进行钻取。
12.根据权利要求10所述的测量储层岩石的脆性指数的装置,其特征在于,所述破裂实验单元包括:
具有所述待测岩心的放置位的放置组件;
设置在所述放置组件上的围限压力施加组件,用于向所述待测岩心施加围限压力;
设置在所述放置组件上的孔隙压力施加组件,用于向所述待测岩心施加孔隙压力;
设置在所述放置组件上的轴向压力施加组件,用于向所述待测岩心施加轴向压力;
设置在所述放置组件上的数据获取组件,用于获取所述待测岩心的应力值和应变值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611192098.0A CN106644736B (zh) | 2016-12-21 | 2016-12-21 | 一种测量储层岩石的脆性指数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611192098.0A CN106644736B (zh) | 2016-12-21 | 2016-12-21 | 一种测量储层岩石的脆性指数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106644736A true CN106644736A (zh) | 2017-05-10 |
CN106644736B CN106644736B (zh) | 2019-11-08 |
Family
ID=58834774
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611192098.0A Active CN106644736B (zh) | 2016-12-21 | 2016-12-21 | 一种测量储层岩石的脆性指数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106644736B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108519281A (zh) * | 2018-02-27 | 2018-09-11 | 中国石油天然气股份有限公司 | 一种确定岩石的脆性指数的方法、装置及系统 |
CN108827774A (zh) * | 2018-06-23 | 2018-11-16 | 东北石油大学 | 煤岩脆性评价方法 |
CN109612840A (zh) * | 2018-12-27 | 2019-04-12 | 东北大学 | 用于获得脆性岩石峰后曲线和残留强度的实验装置及方法 |
CN110501758A (zh) * | 2019-07-30 | 2019-11-26 | 东北大学 | 一种砂砾岩储层纵向连续脆性指数预方法 |
CN110579400A (zh) * | 2019-09-25 | 2019-12-17 | 西南石油大学 | 一种脆性岩石微尺度强度及残余强度的测量计算方法 |
CN110926941A (zh) * | 2019-11-15 | 2020-03-27 | 长江大学 | 一种页岩脆性指数评价方法、装置及系统 |
CN111504779A (zh) * | 2020-04-07 | 2020-08-07 | 武汉大学 | 利用脆性指数确定岩石软化曲线的方法及装置 |
CN112461662A (zh) * | 2020-11-26 | 2021-03-09 | 西南石油大学 | 一种基于应力-应变曲线定量表征岩石脆性的新方法 |
CN113051727A (zh) * | 2021-03-13 | 2021-06-29 | 西华大学 | 基于岩石峰前起裂及峰后应力特征的脆性评价方法及系统 |
-
2016
- 2016-12-21 CN CN201611192098.0A patent/CN106644736B/zh active Active
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108519281A (zh) * | 2018-02-27 | 2018-09-11 | 中国石油天然气股份有限公司 | 一种确定岩石的脆性指数的方法、装置及系统 |
CN108519281B (zh) * | 2018-02-27 | 2020-08-11 | 中国石油天然气股份有限公司 | 一种确定岩石的脆性指数的方法、装置及系统 |
CN108827774A (zh) * | 2018-06-23 | 2018-11-16 | 东北石油大学 | 煤岩脆性评价方法 |
CN109612840A (zh) * | 2018-12-27 | 2019-04-12 | 东北大学 | 用于获得脆性岩石峰后曲线和残留强度的实验装置及方法 |
CN109612840B (zh) * | 2018-12-27 | 2021-03-09 | 东北大学 | 用于获得脆性岩石峰后曲线和残留强度的实验装置及方法 |
CN110501758A (zh) * | 2019-07-30 | 2019-11-26 | 东北大学 | 一种砂砾岩储层纵向连续脆性指数预方法 |
CN110579400A (zh) * | 2019-09-25 | 2019-12-17 | 西南石油大学 | 一种脆性岩石微尺度强度及残余强度的测量计算方法 |
CN110926941A (zh) * | 2019-11-15 | 2020-03-27 | 长江大学 | 一种页岩脆性指数评价方法、装置及系统 |
CN111504779A (zh) * | 2020-04-07 | 2020-08-07 | 武汉大学 | 利用脆性指数确定岩石软化曲线的方法及装置 |
CN112461662A (zh) * | 2020-11-26 | 2021-03-09 | 西南石油大学 | 一种基于应力-应变曲线定量表征岩石脆性的新方法 |
CN113051727A (zh) * | 2021-03-13 | 2021-06-29 | 西华大学 | 基于岩石峰前起裂及峰后应力特征的脆性评价方法及系统 |
CN113051727B (zh) * | 2021-03-13 | 2022-12-23 | 西华大学 | 基于岩石峰前起裂及峰后应力特征的脆性评价方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106644736B (zh) | 2019-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106644736A (zh) | 一种测量储层岩石的脆性指数的方法及装置 | |
Chen et al. | A semianalytical model for pressure-transient analysis of fractured wells in unconventional plays with arbitrarily distributed discrete fractures | |
CN106468172B (zh) | 一种超低渗砂岩油藏低阻储层测井解释方法 | |
Behmanesh et al. | Impact of distance-of-investigation calculations on rate-transient analysis of unconventional gas and light-oil reservoirs: new formulations for linear flow | |
Shimizu et al. | A study of the effect of brittleness on hydraulic fracture complexity using a flow-coupled discrete element method | |
Zhang et al. | Pressure–dependent fracture permeability of marine shales in the Northeast Yunnan area, Southern China | |
Zhang et al. | The stress sensitivity and porosity sensitivity of coal permeability at different depths: a case study in the Pingdingshan mining area | |
CN106153856A (zh) | 一种含裂缝泥页岩稳定性评价装置及方法 | |
CN105527652B (zh) | 一种岩石脆性的测井方法和装置 | |
Liu et al. | An improved capillary pressure model using fractal geometry for coal rock | |
CN105697002A (zh) | 一种用于识别煤系地层岩性的方法 | |
CN104865610A (zh) | 一种深部储层岩心的数值表征方法 | |
Zhang et al. | A study on effect of seepage direction on permeability stress test | |
CN107237626B (zh) | 一种井钻遇溶洞体积的确定方法 | |
CN107842394A (zh) | 大跨度地下硐室开采顶板稳定性的动弹模探测方法 | |
Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
CN109555515A (zh) | 地层坍塌压力确定方法和装置 | |
Esatyana et al. | Nanoindentation of shale cuttings and its application to core measurements | |
CN107728205A (zh) | 一种地层压力预测方法 | |
CN205982211U (zh) | 用于测试钻井液与岩石间压力传递的实验装置 | |
WO2022056194A1 (en) | Method and system for determining energy-based brittleness | |
CN106772674B (zh) | 一种裂缝性砂岩储层含油饱和度的计算方法 | |
Subrahmanyam | Evaluation of hydraulic fracturing and overcoring methods to determine and compare the in situ stress parameters in porous rock mass | |
CN113609702B (zh) | 一种页岩地层井壁天然裂缝扩展压力计算方法 |
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 |