[go: up one dir, main page]

CN110006738B - A method for evaluating rock brittleness based on stress-strain curve and scratch testing - Google Patents

A method for evaluating rock brittleness based on stress-strain curve and scratch testing Download PDF

Info

Publication number
CN110006738B
CN110006738B CN201910153061.4A CN201910153061A CN110006738B CN 110006738 B CN110006738 B CN 110006738B CN 201910153061 A CN201910153061 A CN 201910153061A CN 110006738 B CN110006738 B CN 110006738B
Authority
CN
China
Prior art keywords
rock
stress
value
strain
peak
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
CN201910153061.4A
Other languages
Chinese (zh)
Other versions
CN110006738A (en
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.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201910153061.4A priority Critical patent/CN110006738B/en
Publication of CN110006738A publication Critical patent/CN110006738A/en
Application granted granted Critical
Publication of CN110006738B publication Critical patent/CN110006738B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0019Compressive
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058Kind of property studied
    • G01N2203/0069Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0075Strain-stress relations or elastic constants
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/0202Control of the test
    • G01N2203/0212Theories, calculations
    • G01N2203/0218Calculations based on experimental data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/025Geometry of the test
    • G01N2203/0256Triaxial, i.e. the forces being applied along three normal axes of the specimen
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/06Indicating or recording means; Sensing means
    • G01N2203/067Parameter measured for estimating the property
    • G01N2203/0676Force, weight, load, energy, speed or acceleration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/06Indicating or recording means; Sensing means
    • G01N2203/067Parameter measured for estimating the property
    • G01N2203/0682Spatial dimension, e.g. length, area, angle

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention provides a rock brittleness evaluation method based on a stress-strain curve and a scratch test, which comprises the following steps: evaluating the rock from the ground stress environment and the pore pressure of the rock to obtain the stratum condition of the rock; selecting rock mechanics experiments which accord with the environment of the rock sample to be tested, including conventional triaxial compression experiments or triaxial compression experiments considering pore pressure, and scratch tests; acquiring a pre-peak strain energy density value, a crack initiation stress value, a peak stress value, a residual stress value, and strain magnitude and crack line density corresponding to the peak stress value and the residual stress value; and determining the value range of the effective stress coefficient alpha according to the pore pressure of the tested rock sample, and substituting the pre-peak strain energy value, the initiation stress value, the peak stress value, the residual stress value, the strain values corresponding to the residual stress value and the crack line density into a brittleness index calculation formula to obtain the brittleness index of the rock to be tested. The brittleness evaluation method improves the accuracy and the applicability of rock brittleness evaluation.

Description

一种基于应力应变曲线和划痕测试的岩石脆性评价方法A method for evaluating rock brittleness based on stress-strain curve and scratch testing

技术领域technical field

本发明涉及深层岩石力学技术领域,具体地,涉及一种基于应力应变曲线和划痕测试的岩石脆性评价方法。The invention relates to the technical field of deep rock mechanics, in particular to a rock brittleness evaluation method based on stress-strain curve and scratch test.

背景技术Background technique

经典观点认为,岩石在破裂之前没有或很少发生永久变形为脆性破裂。对于深层岩石来说,其一般具有较高脆性,同时常伴随着高温高压的复杂环境。脆性是岩石很重要的一种性质,其评价对岩石工程有着重要的指导意义。例如在石油工程中,脆性是评价储层地质力学特征及水力压裂裂缝扩展评价的重要指标;而在复杂应力条件下的深部岩体工程中,岩体的脆性是影响岩爆等工程灾害的重要内部因素。目前,国内外学者关于岩石的脆性尚未形成统一的评价标准,如公开号为CN108827774A的中国发明专利申请公开的一种煤岩脆性评价方法,通过建立幂函数分布的岩石损伤本构模型,并考虑煤岩力学特性和割理、裂隙系统分布特征,该方法较为复杂且只适用于煤岩;公开号为CN106547034A的中国发明专利申请公开的一种基于测井资料快速计算岩石脆性的方法,由于测井数据的分析和解释会带来一定误差,其测试结果可靠性存疑。由于岩石的脆性与其力学特性密切相关,因此,对岩石所处环境下的力学特性进行分析可准确评价岩石脆性。The classic view is that rock undergoes little or no permanent deformation to brittle fracture before fracture. For deep rocks, they generally have high brittleness and are often accompanied by a complex environment of high temperature and high pressure. Brittleness is an important property of rock, and its evaluation has important guiding significance for rock engineering. For example, in petroleum engineering, brittleness is an important index for evaluating the geomechanical characteristics of reservoirs and the evaluation of hydraulic fracturing crack propagation; while in deep rock mass engineering under complex stress conditions, the brittleness of rock mass is an important factor affecting engineering disasters such as rock bursts. important internal factors. At present, scholars at home and abroad have not formed a unified evaluation standard for rock brittleness. For example, a coal-rock brittleness evaluation method disclosed in the Chinese invention patent application with publication number CN108827774A is established by establishing a power function distribution rock damage constitutive model, and considering Coal rock mechanical properties and cleat, fracture system distribution characteristics, this method is relatively complex and only suitable for coal rock; a method for fast calculation of rock brittleness based on logging data disclosed in Chinese Patent Application Publication No. CN106547034A. The analysis and interpretation of well data will bring certain errors, and the reliability of the test results is questionable. Because the brittleness of rock is closely related to its mechanical properties, the analysis of the mechanical properties of the rock in the environment can accurately evaluate the brittleness of the rock.

公开号为CN106908322A的中国发明专利申请提供了一种基于全应力—应变曲线的岩石脆性指数评价方法,该方法包括:从地应力、温度、水压等方面对待测岩样所处地质环境进行评价;选取符合待测岩样地质环境的岩石力学实验,包括单轴压缩实验,考虑温度的三轴压缩实验,考虑水压的三轴压缩实验;获取实验过程中各特征应力值以及其对应的应变大小;根据所提出的脆性指数计算方法,将特征应力值及应变带入进行计算,得到待测岩石的脆性指数。但该评价方法不适用于裂缝性岩石的脆性评价、准确性及可靠性低,且目前没有针对裂缝性岩石的脆性科学的评价方法。The Chinese invention patent application with publication number CN106908322A provides a rock brittleness index evaluation method based on the total stress-strain curve, the method includes: evaluating the geological environment where the rock sample to be measured is located from the aspects of in-situ stress, temperature, water pressure, etc. ;Select rock mechanics experiments that conform to the geological environment of the rock sample to be tested, including uniaxial compression experiments, triaxial compression experiments considering temperature, and triaxial compression experiments considering water pressure; obtain each characteristic stress value and its corresponding strain during the experiment According to the proposed calculation method of brittleness index, the characteristic stress value and strain are brought into the calculation to obtain the brittleness index of the rock to be tested. However, this evaluation method is not suitable for the brittleness evaluation of fractured rocks, and the accuracy and reliability are low, and there is currently no scientific evaluation method for the brittleness of fractured rocks.

发明内容SUMMARY OF THE INVENTION

为了解决现有技术中的目前没有针对裂缝性岩石的脆性科学的评价方法的技术问题,本发明提供一种基于应力应变曲线和划痕测试的岩石脆性评价方法,该方法基于岩石三轴压缩实验中全应力应变曲线的全过程和划痕测试进行分析,可以便捷准确地得到岩石的脆性评价结果。In order to solve the technical problem in the prior art that there is currently no scientific evaluation method for fractured rock brittleness, the present invention provides a rock brittleness evaluation method based on stress-strain curve and scratch test, the method is based on rock triaxial compression test The whole process and scratch test of the medium full stress-strain curve can be analyzed, and the brittleness evaluation results of the rock can be obtained conveniently and accurately.

为了实现上述目的,本发明提供一种基于应力应变曲线和划痕测试的岩石脆性评价方法,包括依次执行的如下步骤:In order to achieve the above object, the present invention provides a method for evaluating rock brittleness based on stress-strain curve and scratch test, including the following steps performed in sequence:

(1)从岩石所处的地应力环境和孔隙压力对岩石进行评价,得到岩石所处地层条件;(1) Evaluate the rock from the in-situ stress environment and pore pressure where the rock is located, and obtain the stratum conditions where the rock is located;

(2)选取符合待测岩样所处环境的岩石力学实验,包括常规三轴压缩实验或考虑孔隙压力的三轴压缩实验,和,划痕测试;(2) Select rock mechanics experiments that conform to the environment of the rock samples to be tested, including conventional triaxial compression experiments or triaxial compression experiments considering pore pressure, and scratch tests;

(3)获取实验过程中峰前应变能密度值,以及起裂应力值、峰值应力值和残余应力值及其所对应的应变大小;根据划痕测试结果得到裂缝线密度DF(3) obtaining the pre-peak strain energy density value in the experimental process, as well as the crack initiation stress value, the peak stress value and the residual stress value and their corresponding strain sizes; obtain the crack line density D F according to the scratch test result;

(4)由测试岩样的孔隙压力,确定有效应力系数α的取值范围,之后根据脆性指数计算方法,将峰前应变能值,起裂应力值、峰值应力值、残余应力值及其对应的应变值和裂缝线密度带入进行计算,得到待测岩石的脆性指数;(4) Determine the value range of the effective stress coefficient α from the pore pressure of the test rock sample, and then calculate the pre-peak strain energy value, the crack initiation stress value, the peak stress value, the residual stress value and their corresponding values according to the calculation method of the brittleness index. The strain value and fracture linear density are brought into the calculation to obtain the brittleness index of the rock to be tested;

所述的脆性指数B的计算方法为:The calculation method of the brittleness index B is:

B=B1+B2 B=B 1 +B 2

其中,in,

Figure BDA0001982024980000021
Figure BDA0001982024980000021

Figure BDA0001982024980000022
Figure BDA0001982024980000022

B2=log(DF+1)B 2 =log( DF +1)

其中:in:

σp为峰值应力,εp为峰值应变,σr为残余应力,εr为残余应变,α为有效应力系数,PP为孔隙压力,UD为峰前应变能密度,ε为从原点到峰值点这一段曲线对应的应变值,σ为与ε相应的曲线上的应力值,V为岩心试样体积,DF为裂缝线密度,由划痕测试得到。σ p is the peak stress, ε p is the peak strain, σ r is the residual stress, ε r is the residual strain, α is the effective stress coefficient, P P is the pore pressure, U D is the pre-peak strain energy density, ε is the The strain value corresponding to this section of the curve at the peak point, σ is the stress value on the curve corresponding to ε, V is the core sample volume, and D F is the crack line density, obtained from scratch testing.

优选的是,对于致密性岩石,岩石力学实验须选用常规三轴压缩实验和划痕测试,在三轴压缩实验条件设定中需要根据地层条件设定围压。Preferably, for tight rocks, conventional triaxial compression experiments and scratch tests must be selected for rock mechanics experiments, and confining pressures need to be set according to formation conditions in the setting of triaxial compression experimental conditions.

上述任一方案优选的是,对于裂缝性岩石,岩石力学实验须选用考虑孔隙压力的三轴压缩实验和划痕测试。In any of the above solutions, it is preferable that, for fractured rocks, the rock mechanics experiment must select triaxial compression experiment and scratch test considering pore pressure.

上述任一方案优选的是,常规三轴压缩实验及考虑孔隙压力的三轴压缩实验中,根据地层条件设定围压为地层深度/100Mpa,如地层深度为4000米,那么相应的围压设定为4000/100MPa=40MPa。In any of the above schemes, it is preferable that, in the conventional triaxial compression experiment and the triaxial compression experiment considering pore pressure, the confining pressure is set as formation depth/100Mpa according to the formation conditions. If the formation depth is 4000 meters, then the corresponding confining pressure is set. Set as 4000/100MPa=40MPa.

上述任一方案优选的是,孔隙压力为零,脆性指数计算方法中有效应力系数α取0。In any of the above solutions, preferably, the pore pressure is zero, and the effective stress coefficient α in the brittleness index calculation method is zero.

上述任一方案优选的是,孔隙压力不为零,脆性指数计算方法中有效应力系数α取0.3-0.5。In any of the above solutions, preferably, the pore pressure is not zero, and the effective stress coefficient α in the calculation method of the brittleness index is 0.3-0.5.

上述任一方案优选的是,对于低孔低渗储层,孔隙度在10-15%之间,渗透率5-50×10-3μm2,有效应力系数取0.3。In any of the above solutions, preferably, for a low-porosity and low-permeability reservoir, the porosity is between 10-15%, the permeability is 5-50×10 -3 μm 2 , and the effective stress coefficient is 0.3.

上述任一方案优选的是,对于孔隙度和渗透率较高的储层,比如疏松砂岩储层,孔隙度达到>25%,渗透率>500×10-3μm2,这种有效应力系数可以取0.4-0.5。In any of the above solutions, it is preferred that, for reservoirs with high porosity and permeability, such as unconsolidated sandstone reservoirs, the porosity reaches >25% and the permeability is >500×10 -3 μm 2 . This effective stress coefficient can be Take 0.4-0.5.

本发明的基于应力应变曲线和划痕测试的岩石脆性评价方法,是一种针对岩石所处地质环境并基于三轴压缩实验全应力—应变曲线的脆性指数综合计算方法,不仅适用于致密性岩石的脆性评价,同时还适用于裂缝性岩石的脆性评价,准确性高,可靠性强,适用性广。The rock brittleness evaluation method based on the stress-strain curve and the scratch test of the present invention is a comprehensive calculation method of the brittleness index based on the full stress-strain curve of the triaxial compression experiment according to the geological environment where the rock is located, and is not only suitable for compact rocks It is also suitable for the brittleness evaluation of fractured rocks, with high accuracy, strong reliability and wide applicability.

附图说明Description of drawings

图1为根据本发明的基于应力应变曲线和划痕测试的岩石脆性评价方法一优选实施例的流程示意图。FIG. 1 is a schematic flowchart of a preferred embodiment of a method for evaluating rock brittleness based on stress-strain curves and scratch testing according to the present invention.

图2为峰前应变能密度计算方法示意图。Figure 2 is a schematic diagram of the calculation method of the pre-peak strain energy density.

图3为根据本发明的基于应力应变曲线和划痕测试的岩石脆性评价方法一优选实施例的岩石力学实验的应力-应变曲线。3 is a stress-strain curve of a rock mechanics experiment of a preferred embodiment of the method for evaluating rock brittleness based on stress-strain curve and scratch test according to the present invention.

图4为图3所示实施例的划痕测试结果。FIG. 4 is a scratch test result of the embodiment shown in FIG. 3 .

图5为根据本发明的基于应力应变曲线和划痕测试的岩石脆性评价方法另一优选实施例的岩石力学实验的应力-应变曲线。5 is a stress-strain curve of a rock mechanics experiment of another preferred embodiment of the method for evaluating rock brittleness based on stress-strain curve and scratch test according to the present invention.

图6为图5所示实施例的划痕测试结果。FIG. 6 is a scratch test result of the embodiment shown in FIG. 5 .

图7为本根据本发明的基于应力应变曲线和划痕测试的岩石脆性评价方法又一优选实施例的的岩石力学实验的应力-应变曲线。7 is a stress-strain curve of a rock mechanics experiment according to another preferred embodiment of the method for evaluating rock brittleness based on stress-strain curve and scratch test according to the present invention.

图8为图7所示实施例的划痕测试结果。FIG. 8 is a scratch test result of the embodiment shown in FIG. 7 .

具体实施方式Detailed ways

下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。以下实施例中,用于三轴压缩实验的岩石试样为直径25mm,长度50mm的圆柱体,用于划痕测试的试验岩石试样长度均为200mm的岩心柱。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. The described embodiments are only some, but not all, embodiments of the present invention. In the following embodiments, the rock samples used for the triaxial compression test are cylinders with a diameter of 25 mm and a length of 50 mm, and the test rock samples used for the scratch test are all core columns with a length of 200 mm.

实施例1Example 1

对不同埋深的致密碳酸盐岩进行脆性指数评价,方法步骤如图1所示,具体如下:To evaluate the brittleness index of tight carbonate rocks with different burial depths, the method steps are shown in Figure 1, and the details are as follows:

(一)从岩石所处的地应力环境和孔隙压力对岩石进行评价:对致密碳酸盐岩所处地质环境进行评价:埋深分布广泛,从地表到地下4000米,岩石极其致密,内部几乎不包含流体;(1) Evaluation of the rock from the in-situ stress environment and pore pressure of the rock: Evaluation of the geological environment of the tight carbonate rock: The burial depth is widely distributed, from the surface to the underground 4000 meters, the rock is extremely dense, and the interior is almost does not contain fluids;

(二)选取符合待测岩样所处环境的岩石力学实验:由于岩样所处地质环境中,主要影响因素为地应力,因此选用常规三轴压缩实验,实验数量为4,设定围压分别为10Mpa、20MPa、30MPa和40MPa。对于极其致密的储层,有效应力系数取0。通过RTR-1500三轴实验系统获得4块试样完整的应力应变曲线,如图3所示;通过划痕测试仪获得致密碳酸盐岩的划痕测试结果,如图4所示;(2) Select a rock mechanics experiment that matches the environment of the rock sample to be tested: Since the geological environment where the rock sample is located, the main influencing factor is the in-situ stress, so the conventional triaxial compression experiment is selected, the number of experiments is 4, and the confining pressure is set They are 10MPa, 20MPa, 30MPa and 40MPa respectively. For extremely tight reservoirs, the effective stress coefficient is 0. The complete stress-strain curves of 4 samples were obtained by the RTR-1500 triaxial experimental system, as shown in Figure 3; the scratch test results of the dense carbonate rock were obtained by the scratch tester, as shown in Figure 4;

(三)获取实验过程中峰前应变能密度值,以及起裂应力值、峰值应力值和残余应力值及其所对应的应变大小;根据划痕测试结果得到裂缝线密度DF:峰前应变能密度与应力应变曲线的关系如图2所示,本实施例的应力应变结果如图3所示,结合图2和图3,可以得出四块岩样的应力-应变曲线获得岩样各项应力值、应变值和峰前应变能密度如表1所示。(3) Obtain the pre-peak strain energy density value during the experiment, as well as the crack initiation stress value, the peak stress value and the residual stress value and their corresponding strain values; according to the scratch test results, obtain the crack line density DF : pre-peak strain The relationship between energy density and stress-strain curve is shown in Figure 2, and the stress-strain results of this embodiment are shown in Figure 3. Combining Figures 2 and 3, the stress-strain curves of the four rock samples can be obtained. The term stress value, strain value and pre-peak strain energy density are shown in Table 1.

表1致密碳酸盐岩各项应力值、应变值和峰前应变能密度值Table 1 Various stress values, strain values and pre-peak strain energy density values of tight carbonate rock

Figure BDA0001982024980000041
Figure BDA0001982024980000041

Figure BDA0001982024980000051
Figure BDA0001982024980000051

同时,根据图4,致密碳酸盐岩的裂缝数量为4条,折算为裂缝线密度DF为20条/m。At the same time, according to Fig. 4, the number of fractures in the tight carbonate rock is 4, which is converted into a linear fracture density DF of 20 fractures/m.

(四)由测试岩样的孔隙压力,确定有效应力系数的取值范围,本实施例中,由于岩石是致密碳酸盐岩,几乎不包含流体,孔隙压力接近零,所以有效应力系数取值为0,之后根据脆性指数计算方法,将峰前应变能值,起裂应力值、峰值应力值、残余应力值及其对应的应变值和裂缝线密度带入进行计算,得到待测岩石的脆性指数;(4) Determine the value range of the effective stress coefficient based on the pore pressure of the test rock sample. In this embodiment, since the rock is a tight carbonate rock, it hardly contains fluid, and the pore pressure is close to zero, so the effective stress coefficient is taken as a value. is 0, then according to the calculation method of brittleness index, the pre-peak strain energy value, crack initiation stress value, peak stress value, residual stress value and its corresponding strain value and crack linear density are brought into the calculation to obtain the brittleness of the rock to be tested. index;

所述的脆性指数B的计算公式为:The calculation formula of the brittleness index B is:

B=B1+B2 B=B 1 +B 2

Figure BDA0001982024980000052
Figure BDA0001982024980000052

Figure BDA0001982024980000053
Figure BDA0001982024980000053

B2=log(DF+1)B 2 =log( DF +1)

经过计算,得到致密碳酸盐岩在10MPa、20MPa、30MPa和40MPa下的脆性指数为2.9578、2.7464、2.6643和2.3475。由此可知,对于致密碳酸盐岩,在不同围压条件下,其脆性程度相对大小为:10MPa>20MPa>30MPa>40MPa。After calculation, the brittleness indices of tight carbonate rock at 10MPa, 20MPa, 30MPa and 40MPa are 2.9578, 2.7464, 2.6643 and 2.3475. It can be seen that for tight carbonate rock, under different confining pressure conditions, the relative brittleness degree is: 10MPa>20MPa>30MPa>40MPa.

实施例2Example 2

对某地区某一低孔低渗储层(孔隙度在10-15%之间,渗透率5-50*10-3μm2)的泥页岩进行脆性指数评价,如图1所示,本发明实施例提供的脆性指数计算方法,步骤如下:The brittleness index was evaluated for the shale in a low-porosity and low-permeability reservoir (porosity between 10-15%, permeability 5-50* 10-3 μm 2 ) in a certain area, as shown in Fig. 1. The brittleness index calculation method provided by the embodiment of the invention, the steps are as follows:

(一)对泥页岩所处地质环境进行评价:埋深在4000米到4200米,岩石具有低孔隙度低渗透率特性,区域孔隙压力差异性大,压力系数从1.1到1.5不等。(1) Evaluate the geological environment of mud shale: the burial depth is 4000 meters to 4200 meters, the rock has the characteristics of low porosity and low permeability, and the regional pore pressure varies greatly, and the pressure coefficient ranges from 1.1 to 1.5.

(二)由于岩样所处地质环境中,主要影响因素为地应力和孔隙压力,因此选用三轴压缩实验,实验数量为4,设定围压为80MPa,孔隙压力分别为45MPa,50MPa,55MPa和60MPa。对于低孔低渗储层,有效应力系数取0.3。通过RTR-1500三轴实验系统获得4块试样完整的应力应变曲线,如图5所示;(2) The main influencing factors are in-situ stress and pore pressure in the geological environment where the rock sample is located. Therefore, the triaxial compression experiment was selected, the number of experiments was 4, the confining pressure was set to 80MPa, and the pore pressures were 45MPa, 50MPa, and 55MPa respectively. and 60MPa. For low-porosity and low-permeability reservoirs, the effective stress coefficient is taken as 0.3. The complete stress-strain curves of 4 samples were obtained through the RTR-1500 triaxial experimental system, as shown in Figure 5;

(三)根据图5,四块岩样的应力-应变曲线获得岩样各项应力值、应变值和峰前应变能密度如表2所示。(3) According to Figure 5, the stress-strain curves of the four rock samples obtained various stress values, strain values and pre-peak strain energy density of the rock samples, as shown in Table 2.

表2某一低孔低渗储层的泥页岩各项应力值、应变值和峰前应变能密度值Table 2 Various stress values, strain values and pre-peak strain energy density values of shale in a low-porosity and low-permeability reservoir

Figure BDA0001982024980000061
Figure BDA0001982024980000061

同时,根据图6,本实施例中的低孔低渗储层的泥页岩的裂缝数量为10条,折算为裂缝线密度DF为50条/m。Meanwhile, according to FIG. 6 , the number of fractures in the shale of the low-porosity and low-permeability reservoir in this embodiment is 10, which is converted into a linear fracture density D F of 50 fractures/m.

(四)由测试岩样的孔隙压力,确定有效应力系数的取值范围,本实施例中,由于孔隙压力不为零,对于低孔低渗储层,有效应力系数取0.3。然后,将峰前应变能值,起裂应力值、峰值应力值、残余应力值及其对应的应变值和裂缝线密度带入进行计算,得到待测岩石的脆性指数;(4) Determine the value range of the effective stress coefficient according to the pore pressure of the test rock sample. In this embodiment, since the pore pressure is not zero, for low-porosity and low-permeability reservoirs, the effective stress coefficient is 0.3. Then, the pre-peak strain energy value, the crack initiation stress value, the peak stress value, the residual stress value and its corresponding strain value and crack linear density are brought into the calculation to obtain the brittleness index of the rock to be tested;

脆性指B的计算方法为:The calculation method of brittleness index B is:

B=B1+B2 B=B 1 +B 2

Figure BDA0001982024980000062
Figure BDA0001982024980000062

Figure BDA0001982024980000063
Figure BDA0001982024980000063

B2=log(DF+1)B 2 =log( DF +1)

得到泥页岩在80MPa围压条件下,不同孔隙压力45MPa,50MPa,55MPa和60MPa下的脆性指数为3.6157,3.5357,3.3509和3.0108。由此可知,对于泥页岩,在80MPa围压条件下,不同孔隙压力条件下,其脆性程度相对大小为:孔压45MPa>孔压50MPa>孔压55MPa>孔压60MPa。The brittleness indices of shale under the confining pressure of 80MPa and different pore pressures of 45MPa, 50MPa, 55MPa and 60MPa are 3.6157, 3.5357, 3.3509 and 3.0108. It can be seen that for mud shale, under the condition of 80MPa confining pressure and different pore pressures, the relative brittleness degree is as follows: pore pressure 45MPa>pore pressure 50MPa>pore pressure 55MPa>pore pressure 60MPa.

实施例3Example 3

对致密碳酸盐岩进行脆性指数评价,方法步骤如图1所示,具体如下:The brittleness index evaluation of tight carbonate rock, the method steps are shown in Figure 1, and the details are as follows:

(一)对该碳酸盐岩所处地质环境进行评价:该段碳酸盐岩较为致密,所处地下环境中无流体存在,其埋深约4000m;(1) Evaluate the geological environment where the carbonate rock is located: the carbonate rock in this section is relatively dense, there is no fluid in the underground environment, and its burial depth is about 4000m;

(二)根据岩石所处环境进行三轴压缩实验的条件设置:由于岩样所处地质环境中,无流体存在,因此无需考虑流体孔隙压力的影响,其主要影响因素为地应力,因此选用常规三轴压缩实验。(2) Setting the conditions for the triaxial compression experiment according to the environment where the rock is located: Since there is no fluid in the geological environment where the rock sample is located, there is no need to consider the influence of fluid pore pressure. The main influencing factor is in-situ stress, so conventional Triaxial compression experiment.

该例中实验数量为1,根据岩石埋深设定围压为40MPa。该例中由于无流体存在,空隙压力为零,因此有效应力系数取0。通过RTR-1500三轴实验系统获得该试样完整的应力应变曲线,如图7所示。In this example, the number of experiments is 1, and the confining pressure is set to 40 MPa according to the burial depth of the rock. In this example, since there is no fluid, the void pressure is zero, so the effective stress coefficient is taken as zero. The complete stress-strain curve of the sample was obtained by the RTR-1500 triaxial experimental system, as shown in Figure 7.

对其进行划痕测试,得到实验结果如图8所示(三)获取实验过程中峰前应变能密度值,以及起裂应力值、峰值应力值和残余应力值及其所对应的应变大小;根据划痕测试结果得到裂缝线密度DF:峰前应变能密度与应力应变曲线的关系2所示,本实施例的应力应变结果如图7所示,结合图2和图7,可以得出岩样的应力-应变曲线获得岩样各项应力值、应变值和峰前应变能密度如表3所示。Scratch test is carried out on it, and the experimental results are shown in Figure 8. (3) Obtain the pre-peak strain energy density value during the experiment, as well as the crack initiation stress value, peak stress value and residual stress value and their corresponding strain values; According to the scratch test results, the crack line density DF is obtained: the relationship between the pre-peak strain energy density and the stress-strain curve is shown in Figure 2. The stress-strain results of this embodiment are shown in Figure 7. Combining Figures 2 and 7, it can be obtained From the stress-strain curve of the rock sample, the stress values, strain values and pre-peak strain energy density of the rock samples are obtained as shown in Table 3.

表3致密碳酸盐岩各项应力值、应变值和峰前应变能密度值Table 3 Various stress values, strain values and pre-peak strain energy density values of tight carbonate rock

Figure BDA0001982024980000071
Figure BDA0001982024980000071

同时,根据图8,该实验岩心段长度为20cm,根据连续强度实验结果可以观察到,该碳酸盐岩的裂缝数量为5条,因此,折算为线密度为25条/m。At the same time, according to Figure 8, the length of the experimental core section is 20cm. According to the continuous strength test results, it can be observed that the number of fractures in the carbonate rock is 5, so the linear density is converted to 25/m.

(四)由测试岩样的孔隙压力,确定有效应力系数的取值范围,本实施例中,由于岩石是致密碳酸盐岩,无流体存在,孔隙压力为零,所以有效应力系数取值为0,之后根据脆性指数计算方法,将峰前应变能值,起裂应力值、峰值应力值、残余应力值及其对应的应变值和裂缝线密度带入进行计算,得到待测岩石的脆性指数;(4) Determine the value range of the effective stress coefficient from the pore pressure of the test rock sample. In this embodiment, since the rock is a tight carbonate rock, there is no fluid, and the pore pressure is zero, the effective stress coefficient is taken as 0, and then according to the calculation method of brittleness index, the pre-peak strain energy value, crack initiation stress value, peak stress value, residual stress value and its corresponding strain value and crack linear density are brought into the calculation to obtain the brittleness index of the rock to be tested. ;

带入脆性指数B的计算方法,Bring in the calculation method of the brittleness index B,

B=Bl+B2 B=B l +B 2

Figure BDA0001982024980000081
Figure BDA0001982024980000081

Figure BDA0001982024980000082
Figure BDA0001982024980000082

B2=log(DF+1)B 2 =log( DF +1)

得到该碳酸盐岩在40MPa下的脆性指数为2.278。The brittleness index of the carbonate rock at 40MPa is 2.278.

关于本发明提供的脆性指数评价方法的准确性和可靠性,与现有技术中的评价方法进行了对比,结果表明,现有技术中的评价方法,皆会存在满足一定特征的一些应力应变曲线,导致其无法准确评价其脆性,会出现不同曲线却得到同样的结果等异常结论;而本发明中的脆性评价方法,可以避免出现异常结论,准确进行评价。因此,本发明的方法准确性更好、可靠性更强。The accuracy and reliability of the brittleness index evaluation method provided by the present invention are compared with the evaluation methods in the prior art. The results show that the evaluation methods in the prior art all have some stress-strain curves that meet certain characteristics. , which makes it impossible to accurately evaluate its brittleness, and there will be abnormal conclusions such as different curves but the same results are obtained; and the brittleness evaluation method in the present invention can avoid abnormal conclusions and accurately evaluate. Therefore, the method of the present invention is more accurate and more reliable.

以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The specific embodiments described above further describe the purpose, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above-mentioned specific embodiments are only specific embodiments of the present invention, and are not intended to limit the scope of the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims (8)

1.一种基于应力应变曲线和划痕测试的岩石脆性评价方法,包括依次执行的如下步骤:1. A method for evaluating rock brittleness based on stress-strain curve and scratch testing, comprising the following steps performed in sequence: (1)从岩石所处的地应力环境和孔隙压力对岩石进行评价,得到岩石所处地层条件;(1) Evaluate the rock from the in-situ stress environment and pore pressure where the rock is located, and obtain the stratum conditions where the rock is located; (2)选取符合待测岩样所处环境的岩石力学实验,包括三轴压缩实验和划痕测试,所述三轴压缩实验包括常规三轴压缩实验或考虑孔隙压力的三轴压缩实验;(2) Selecting rock mechanics experiments that conform to the environment where the rock samples to be tested are located, including triaxial compression experiments and scratch tests, where the triaxial compression experiments include conventional triaxial compression experiments or triaxial compression experiments considering pore pressure; (3)获取实验过程中峰前应变能密度值,以及起裂应力值、峰值应力值和残余应力值及其所对应的应变大小;根据划痕测试结果得到裂缝线密度DF(3) obtaining the pre-peak strain energy density value in the experimental process, as well as the crack initiation stress value, the peak stress value and the residual stress value and their corresponding strain sizes; obtain the crack line density D F according to the scratch test result; (4)由测试岩样的孔隙压力,确定有效应力系数α的取值范围,之后根据脆性指数计算方法,将峰前应变能值,起裂应力值、峰值应力值、残余应力值及其对应的应变值和裂缝线密度代入进行计算,得到待测岩石的脆性指数;(4) Determine the value range of the effective stress coefficient α from the pore pressure of the test rock sample, and then calculate the pre-peak strain energy value, the crack initiation stress value, the peak stress value, the residual stress value and their corresponding values according to the calculation method of the brittleness index. The strain value and fracture linear density are substituted into the calculation to obtain the brittleness index of the rock to be tested; 所述的脆性指数B的计算方法为:The calculation method of the brittleness index B is: B=B1+B2 B=B 1 +B 2 其中,in,
Figure FDA0002486772810000011
Figure FDA0002486772810000011
Figure FDA0002486772810000012
Figure FDA0002486772810000012
B2=log(DF+1)B 2 =log( DF +1) 其中:in: σp为峰值应力,εp为峰值应变,σr为残余应力,εr为残余应变,α为有效应力系数,PP为孔隙压力,UD为峰前应变能密度,ε为从原点到峰值点这一段曲线对应的应变值,σ为与ε相应的曲线上的应力值,V为岩心试样体积,DF为裂缝线密度,由划痕测试得到。σ p is the peak stress, ε p is the peak strain, σ r is the residual stress, ε r is the residual strain, α is the effective stress coefficient, P P is the pore pressure, U D is the pre-peak strain energy density, ε is the The strain value corresponding to this section of the curve at the peak point, σ is the stress value on the curve corresponding to ε, V is the core sample volume, and D F is the crack line density, obtained from scratch testing.
2.如权利要求1所述的评价方法,其特征在于,对于致密性岩石,岩石力学实验须选用常规三轴压缩实验和划痕测试,在三轴压缩实验条件设定中需要根据地层条件设定围压。2. The evaluation method as claimed in claim 1 is characterized in that, for compact rock, conventional triaxial compression experiment and scratch test must be selected for rock mechanics experiment, and in the setting of triaxial compression experiment conditions, it is necessary to set according to formation conditions. Fixed confining pressure. 3.如权利要求1所述的评价方法,其特征在于,对于裂缝性岩石,岩石力学实验须选用考虑孔隙压力的三轴压缩实验和划痕测试。3. The evaluation method according to claim 1, characterized in that, for fractured rock, triaxial compression test and scratch test considering pore pressure must be selected for rock mechanics test. 4.如权利要求2或3所述的评价方法,其特征在于,常规三轴压缩实验或考虑孔隙压力的三轴压缩实验中,根据地层条件设定围压为地层深度/100Mpa。The evaluation method according to claim 2 or 3, characterized in that, in a conventional triaxial compression experiment or a triaxial compression experiment considering pore pressure, the confining pressure is set as formation depth/100Mpa according to formation conditions. 5.如权利要求4所述的评价方法,其特征在于,孔隙压力为零,脆性指数计算方法中有效应力系数α取0。5 . The evaluation method according to claim 4 , wherein the pore pressure is zero, and the effective stress coefficient α in the brittleness index calculation method is zero. 6 . 6.如权利要求4所述的评价方法,其特征在于,孔隙压力不为零,脆性指数计算方法中有效应力系数α取0.3-0.5。6 . The evaluation method according to claim 4 , wherein the pore pressure is not zero, and the effective stress coefficient α in the brittleness index calculation method is 0.3-0.5. 7 . 7.如权利要求6所述的评价方法,其特征在于,对于孔隙度在10-15%之间,渗透率5-50×10-3μm2的储层,有效应力系数取0.3。7 . The evaluation method according to claim 6 , wherein, for a reservoir with a porosity of 10-15% and a permeability of 5-50×10 -3 μm 2 , the effective stress coefficient is 0.3. 8 . 8.如权利要求6所述的评价方法,其特征在于,对于孔隙度达到>25%,渗透率>500×10-3μm2的储层,有效应力系数取0.4-0.5。8 . The evaluation method according to claim 6 , wherein, for a reservoir with a porosity of >25% and a permeability of >500×10 −3 μm 2 , the effective stress coefficient is 0.4-0.5. 9 .
CN201910153061.4A 2019-02-28 2019-02-28 A method for evaluating rock brittleness based on stress-strain curve and scratch testing Active CN110006738B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910153061.4A CN110006738B (en) 2019-02-28 2019-02-28 A method for evaluating rock brittleness based on stress-strain curve and scratch testing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910153061.4A CN110006738B (en) 2019-02-28 2019-02-28 A method for evaluating rock brittleness based on stress-strain curve and scratch testing

Publications (2)

Publication Number Publication Date
CN110006738A CN110006738A (en) 2019-07-12
CN110006738B true CN110006738B (en) 2020-08-18

Family

ID=67166265

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910153061.4A Active CN110006738B (en) 2019-02-28 2019-02-28 A method for evaluating rock brittleness based on stress-strain curve and scratch testing

Country Status (1)

Country Link
CN (1) CN110006738B (en)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110864966B (en) * 2019-10-12 2022-08-05 中国石油天然气股份有限公司 Compressibility evaluation method suitable for fractured rock
CN110726608B (en) * 2019-11-08 2022-04-05 西南石油大学 A shale brittleness evaluation method based on energy evolution of stress-strain curve
CN112881184A (en) * 2019-11-29 2021-06-01 中国石油化工股份有限公司 Gravel brittleness index determination method and system based on whole-course stress strain
CN111238931B (en) * 2019-12-30 2023-08-22 长江大学 Shale brittleness index evaluation method based on energy evolution
CN111238939A (en) * 2020-03-17 2020-06-05 中联煤层气有限责任公司 Method and device for determining rock brittleness index
CN111504780B (en) * 2020-04-07 2021-02-02 武汉大学 Method and device for determining rock softening curve in bonding crack model
CN111504779B (en) * 2020-04-07 2021-04-16 武汉大学 Method and device for determining rock softening curve by using brittleness index
CN111983194A (en) * 2020-08-31 2020-11-24 长安大学 Oil and gas reservoir rock fracturing experimental analysis method
US12050297B2 (en) 2020-09-11 2024-07-30 Saudi Arabian Oil Company Method and system for determining energy-based brittleness
CN112461662B (en) * 2020-11-26 2022-05-06 西南石油大学 Novel method for quantitatively representing rock brittleness based on stress-strain curve
CN112505085B (en) * 2021-02-05 2021-04-09 西南石油大学 Method for measuring porosity effective stress coefficient based on nuclear magnetic resonance
CN113051727B (en) * 2021-03-13 2022-12-23 西华大学 Brittleness evaluation method and system based on rock pre-peak crack initiation and post-peak stress characteristics
CN114813303B (en) * 2022-02-28 2025-01-28 陕西延长石油(集团)有限责任公司 A calculation method for rock brittleness index of tight sandstone reservoirs
CN115356223B (en) * 2022-10-20 2022-12-20 中国矿业大学(北京) Device and method for measuring shale brittleness index continuous section based on high-temperature and high-pressure scratches

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104406849B (en) * 2014-11-21 2017-01-04 中国石油天然气股份有限公司 Method and device for predicting brittleness of reservoir rock
CN104730596B (en) * 2015-01-25 2017-11-07 中国石油大学(华东) A kind of discrete fractures modeling method based on multiple dimensioned factor constraints
CN106094052A (en) * 2016-06-01 2016-11-09 中国地质大学(武汉) A kind of development degree of micro cracks in oil recognition methods of dense dolomite reservoir
JP6724761B2 (en) * 2016-06-08 2020-07-15 日本製鉄株式会社 Hydrogen embrittlement evaluation apparatus, hydrogen embrittlement evaluation method, and test piece used therein
CN106908322B (en) * 2017-02-23 2019-04-23 成都理工大学 An evaluation method of rock brittleness index based on total stress-strain curve

Also Published As

Publication number Publication date
CN110006738A (en) 2019-07-12

Similar Documents

Publication Publication Date Title
CN110006738B (en) A method for evaluating rock brittleness based on stress-strain curve and scratch testing
Wang et al. Permeability evolution of coal subjected to triaxial compression based on in-situ nuclear magnetic resonance
Peng et al. Numerical investigation of brittleness effect on strength and microcracking behavior of crystalline rock
Bohnsack et al. Stress sensitivity of porosity and permeability under varying hydrostatic stress conditions for different carbonate rock types of the geothermal Malm reservoir in Southern Germany
Hou et al. Effect of water and nitrogen fracturing fluids on initiation and extension of fracture in hydraulic fracturing of porous rock
Kamali-Asl et al. Elastic, viscoelastic, and strength properties of Marcellus Shale specimens
Liu et al. Experimental study on the difference of shale mechanical properties
Al-Shayea Comparing reservoir and outcrop specimens for mixed mode I–II fracture toughness of a limestone rock formation at various conditions
Ye et al. Laboratory investigation of fluid flow and permeability evolution through shale fractures
Zheng et al. Assessment on anisotropy degree and fracture modes for weakly anisotropic sandstone using the acoustic emission technique
Ma et al. Nuclear magnetic resonance experiment on the influence of confining pressure on spontaneous imbibition of water in coal
Lu et al. Experimental study of dynamic permeability changes in coals of various ranks during hydraulic fracturing
Chen et al. Stress sensitivity of multiscale pore structure of shale gas reservoir under fracturing fluid imbibition
Ye et al. Failure behavior of the Poorman schist and its fractures from EGS collab stimulation site
Goyal et al. A comparative study of monotonic and cyclic injection hydraulic fracturing in saturated tight rocks under triaxial stress
Hou et al. Mechanical properties and brittleness of shale with different degrees of fracturing-fluid saturation
Shi et al. Effect of the intermediate principal stress on the evolution of mudstone permeability under true triaxial compression
Sun et al. Pore structure evolution of mudstone caprock under cyclic load-unload and its influence on breakthrough pressure
Zeng et al. Experimental Determination of Geomechanical and Petrophysical Properties of Jackfork Sandstone-A Tight Gas Formation
Zhao et al. Geomechanical properties of laminated shale and bedding shale after water absorption: A case study of the Chang 7 shale in Ordos basin, China
Ren et al. Insight into the pore structures and its impacts on movable fluid in tight sandstones
Wei et al. Experimental study on the anisotropic characteristics and engineering application of tight sandstone
Liu et al. A new method for wellbore stability evaluation based on fractured carbonate reservoir rock breaking degree
Chavez et al. Field and experimental brittleness (toughness) determination of Vaca Muerta shale
Mighani et al. Stress dependency of rock tensile strength

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