一种基于全应力‑应变曲线的岩石脆性指数评价方法

摘要

本发明公开了一种基于全应力‑应变曲线的岩石脆性指数评价方法，该方法包括：从地应力，温度，水压等方面对待测岩样所处地质环境进行评价；选取符合待测岩样地质环境的岩石力学实验，包括单轴压缩实验，考虑温度的三轴压缩实验，考虑水压的三轴压缩实验；获取实验过程中各特征应力值以及其对应的应变大小；根据所提出的脆性指数计算方法，将特征应力值及应变带入进行计算，得到待测岩石的脆性指数。本方法是一种针对岩石所处地质环境并基于岩石力学实验全过程的脆性指数综合计算方法，提高了岩石脆性评价的合理性与准确性。

说明书

技术领域

本发明涉及岩石力学技术领域，具体地，涉及一种基于全应力-应变曲线的岩石脆性指数评价方法。

背景技术

脆性是岩石在极小的形变条件下发生破坏，并释放大量能量，深部岩体一般具有较高脆性，同时常伴随着高地温高水压的恶劣环境。作为岩石极为重要的一种性质，脆性评价对岩体工程有着重要的指导意义，例如在高地应力下的深部岩体工程中，岩体的脆性是影响岩爆等工程灾害的重要内部因素；脆性指标也是油气藏工程领域中储层力学特征、井壁稳定性评价以及水力压裂效果评价的重要指标；在岩体开挖进程中，岩石的脆性也决定了TBM掘进速度和钻机的钻孔效率。目前，国内外学者关于岩石脆性尚未形成统一的评价标准，如专利号201410461329.8公开的一种基于煤岩工业组分的脆性指数确定方法，采用煤岩工业组分结合测井资料方法确定脆性指数，其方法繁琐且只能适用于煤气领域，专利号201310254628.X公开的一种应用声发射能量值测定岩石脆性指数的方法，由于岩石声发射能量值具有较大的随机性，其测试结果不稳定，具有较大误差，但由于岩石的脆性与其力学特性密切相关，因此，对岩体所处环境下的力学特性进行研究可实现对脆性指标的准确评价。

发明内容

本发明目的在于提供一种基于全应力-应变曲线的岩石脆性指数评价方法，该方法基于岩石单轴实验或三轴实验的应力-应变曲线全过程的分析，可提高岩石脆性评价的合理性与准确性。

为实现上述目的，本发明实施例提供一种基于全应力-应变曲线的岩石脆性指数评价方法，包括：

从地应力，温度，水压方面对待测岩样所处地质环境进行评价；

选取符合待测岩样地质环境的岩石力学实验，包括单轴压缩实验，考虑温度的三轴压缩实验，考虑水压的三轴压缩实验；

获取实验过程中各特征应力值以及其对应的应变大小，各特征应力值包括起裂应力、峰值应力、残余应力；

根据所提出的脆性指数计算方法，将特征应力值及应变带入进行计算，得到待测岩石的脆性指数。

所述的脆性指数计算方法为：

其中：σ_i为起裂应力；σ_p为峰值应力；σ_r为残余应力；ε_i为起裂应变；ε_p为峰值应变；ε_r为残余应变；t为三轴压缩实验所设置温度；t₂₀为外界环境温度，一般为20℃；D为三轴压缩实验所设置水压；k为系数，待测岩石为硬岩时取4，待测岩石为软岩时取2。

为了在实验中获取更好的峰后曲线，实验加载系统采用三段式控制，在实验开始时采用轴向位移控制，当轴向应力达到峰值强度70％-85％时，切换为环向位移控制，环向位移速率逐级递增，初始设定值为0.0025-0.0040mm/s，当应力达到峰值并开始下降至峰值强度的80％-90％时，将环向位移速率调整为初始值的1.5-2倍，当下降至峰值强度的50％-60％时，将环向位移速率调整为初始值的2.5-3倍，当下降至峰值强度的30％-40％时，调整为初始值的3.5-5倍，当下降至峰值强度10％-20％时，再次切换为轴向位移控制，来获得最终完整的峰后曲线。当试验岩样为硬岩时，在范围内取较大值；当试验岩样为软岩时，在范围内去取较小值。

结合岩石所处地质环境，当岩石无需考虑其他外部因素时，岩石力学实验选用单轴压缩实验，本发明所述的脆性指数计算方法采用B₁。

结合岩石所处地质环境，当岩石需要考虑温度因素时，一般伴随较高的地应力，岩石力学实验须选用考虑温度的三轴压缩实验，本发明所述的脆性指数计算方法采用B₂。

结合岩石所处地质环境，当岩石需要考虑水压因素时，一般伴随较高的地应力，岩石力学实验选用考虑水压的三轴压缩实验，本发明所述的脆性指数计算方法采用B₃。

本发明提供了一种基于全应力-应变曲线的岩石脆性指数计算方法，该方法模拟了岩样所处地质环境，同时也考虑到了岩石力学实验整个过程对岩样脆性指数的影响；相比于现有技术，本发明方法针对性强，操作简单，提高了岩石脆性评价的合理性与准确性。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单介绍。

图1为本发明的流程示意图；

图2为本发明实施例提供的四种岩样岩石力学实验的应力-应变曲线；

图3为本发明实施例提供的起裂应力值计算方法示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清除、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。

实施例一

对地表出露的大理岩、砂岩、粗颗粒花岗岩、细颗粒花岗岩进行脆性指数评价，如图1所示，本发明实施例提供的脆性指数计算方法，包括以下步骤：

对四种岩石所处地质环境进行评价：地表出露，地应力较小，且无温度、水压等影响；

由于岩样所处地质环境简单，选用单轴压缩实验，实验加载系统采用三段式控制，在实验开始时采用轴向位移控制，当轴向压力达到峰值强度70％时，切换为环向位移控制，环向位移速率逐级递增，初始设定值为0.0025mm/s，当曲线达到峰值并开始下降至峰值强度的90％时，将环向位移速率调整为0.005mm/s，当下降至峰值强度的60％时，将环向位移速率调整为0.01mm/s，当下降至峰值强度的40％时，调整为0.02mm/s，当下降至峰值强度10％时，再次切换为轴向位移控制，来获得最终完整的峰后曲线；

如图2所示，获得四种岩样在单轴压缩条件下的应力-应变曲线；

如图3所示，根据已有的技术方法，采用裂纹体积应变拐点法确定起裂应力值，四种岩样的应力-应变曲线获得岩样各特征应力值及应变值如表1。

表1实施例一的具体特征应力值及其应变值

带入单轴压缩的脆性指数计算方法：

最后，得到大理岩的脆性指数为B_大＝1.059，砂岩的脆性指数为B_砂＝11.799，粗颗粒花岗岩的脆性指数为B_粗＝9.54，细颗粒花岗岩的脆性指数为B_细＝11.858，因此可知，四种岩石的脆性程度为细颗粒花岗岩＞砂岩＞粗颗粒花岗岩＞大理岩。

实施例二

对20℃，40℃，60℃，90℃下大理岩进行脆性指数评价，如图1所示，本发明实施例提供的脆性指数计算方法，包括以下步骤：

对四种温度下岩石所处地质环境进行评价：高地温岩石一般深埋地下，地应力较大，无水压影响；

选用三轴压缩试验，由于地应力较大，围压选用20MPa；

获取实验过程中各特征应力值以及其对应的应变大小如表2所示，各特征应力值包括起裂应力，峰值应力，残余应力；

表2实施例二的具体特征应力值及其应变值

带入考虑温度三轴压缩的脆性指数计算方法：

最后得到20℃的脆性指数为B_20℃＝4.156，40℃的脆性指数为B_40℃＝9.716，60℃的脆性指数为B_60℃＝17.074，90℃的脆性指数为B_90℃＝35.747，因此可知，四种温度下大理岩的脆性程度为90℃＞60℃＞40℃＞20℃。

实施例三

对水压为0MPa，4MPa，8MPa，12MPa下大理岩进行脆性指数评价，如图1所示，本发明实施例提供的脆性指数计算方法，包括以下步骤：

对四种水压下岩石所处地质环境进行评价：高水压岩石一般深埋地下，地下水丰富，地应力较大，无温度影响；

选用三轴压缩试验，由于地应力较大，围压选用20MPa；

获取实验过程中各特征应力值以及其对应的应变大小如表3所示，各特征应力值包括起裂应力，峰值应力，残余应力；

表3实施例三的具体特征应力值及其应变值

带入考虑水压三轴压缩的脆性指数计算方法，其中k取4：

最后得到水压0MPa的脆性指数为B_0MPa＝4.271，4MPa的脆性指数为B_4MPa＝4.502，8MPa的脆性指数为B_8MPa＝5.040，12MPa的脆性指数为B_12MPa＝5.688，因此可知，四种水压下大理岩的脆性程度为12MPa＞8MPa＞4MPa＞0MPa。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。