一种确定岩石初始损伤程度的方法

技术领域

本发明涉及岩石损伤确定方法，特别涉及一种确定岩石初始损伤程度的方法。

背景技术

岩石在漫长复杂的成岩过程中及外界各种复合扰动如温度、湿度、高地应力作用下，其内部或多或少会产生孔隙、裂隙、节理弱面、断层等原生缺陷，致使岩石具有初始损伤。相较于理想完整无损岩块，这些初始缺陷或弱结构面在岩石基体内部交错分布，不仅改变了其宏观变形力学参数如弹性模量、抗压强度等，更对其损伤演变过程及本构关系产生显著影响。因此，确定岩石初始损伤程度对深刻理解岩石变形破坏过程具有重要意义。

目前，针对岩石损伤演化方面的研究，大多假设岩石初始状态无损伤，这显然与实际情况不符。因此，亟需研制出一种确定岩石初始损伤程度的方法。

发明内容

经研究，发明人发现：从岩石应力-应变曲线上看，正是这些原生初始缺陷的存在导致应力-应变曲线前期的非线性变形，且岩石初始损伤程度越高，非线性特征越强烈，弹性模量也相应降低。从能量演化角度看，岩石压缩变形破坏过程实质上是外界能量输入、可释放弹性能积聚及能量耗散的一系列复杂演化过程。从岩石破裂过程的细观角度看，其变形破坏过程又是内部原生微缺陷压缩闭合、新裂纹萌生及新旧微裂纹不断发展并最终汇通形成宏观裂缝的过程。而耗散能的演化过程恰好代表并反映了岩石内部微缺陷及损伤程度演变规律。故岩石初始压密阶段由于初始缺陷压缩闭合产生的摩擦作用会消耗部分能量，该阶段耗散的能量反应了岩石初始损伤程度。基于此，本发明提出了一种确定岩石初始损伤程度的方法。

为实现本发明目的，本发明采用如下技术方案：

本发明提供了一种确定岩石初始损伤程度的方法，包括如下步骤：

步骤100：将工程现场取回的岩样加工成试样，对试样进行单轴压缩试验，获得该试样的应力-应变曲线及弹性模量

步骤200：根据已建立的理论本构模型得到该试样的单轴压缩条件下应力-应变全过程曲线；

步骤300：根据步骤200中得到的理论应力-应变全过程曲线，结合能量计算公式可得到总能量、弹性能和耗散能演化曲线；

步骤400：根据步骤300中的耗散能演化曲线可确定岩石初始压密阶段的耗散能

，

其中，

即可确定岩石初始损伤程度。

根据本发明的一实施方式，其中，所述试样是将所述岩样按照国际岩石力学学会标准加工成直径为50mm、高度为100mm的标准圆柱体。

根据本发明的一实施方式，其中，所述弹性模量

根据本发明的一实施方式，其中，所述理论本构模型为：

，

式中，

为单轴方向上的主应力，MPa；/>

为裂纹闭合应变，根据应力-应变曲线弹性阶段反向延长线与横坐标轴的交点确定；

为裂纹闭合阶段的等效弹性模量，根据裂纹应变-应力曲线拟合结果求得，MPa；

E

D

为裂纹闭合应力/>

根据本发明的一实施方式，其中，所述能量计算公式为单轴压缩条件下的能量公式：

，

，

，

，

式中，

为岩石总共吸收的能量，MJ/m

为岩石变形过程中内部储存的弹性应变能，MJ/m

为加载过程中释放耗散的能量，MJ/m

为单轴方向上的主应力，MPa；/>

E

本发明中的一个实施例具有如下优点或有益效果：

本发明的确定岩石初始损伤程度的方法，包括：将工程现场的岩样加工成试样，对试样进行单轴压缩试验，获得该试样的应力-应变曲线及弹性模量

附图说明

通过参照附图详细描述其示例实施方式，本发明的上述和其它特征及优点将变得更加明显。

图1是根据一示例性实施方式示出砂岩单轴压缩过程中能量演化示意图；

图2是根据一示例性实施方式示出基于变模量衰减法获得损伤变量示意图；

图3是根据一示例性实施方式示出损伤演化曲线及拟合结果示意图；

图4是根据一示例性实施方式示出砂岩单轴压缩应力-应变曲线拟合结果示意图；

图5是根据一示例性实施方式示出砂岩理论能量演化全过程的示意图。

具体实施方式

现在将参考附图更全面地描述示例实施方式。然而，示例实施方式能够以多种形式实施，且不应被理解为限于在此阐述的实施方式；相反，提供这些实施方式使得本发明将全面和完整，并将示例实施方式的构思全面地传达给本领域的技术人员。图中相同的附图标记表示相同或类似的结构，因而将省略它们的详细描述。

用语“一个”、“一”、“该”、“所述”用以表示存在一个或多个要素/组成部分/等；用语“包括”和“具有”用以表示开放式的包括在内的意思并且是指除了列出的要素/组成部分/等之外还可存在另外的要素/组成部分/等。

如图1-图5所示，图1是根据一示例性实施方式示出砂岩单轴压缩过程中能量演化示意图。图2是根据一示例性实施方式示出基于变模量衰减法获得损伤变量示意图。图3是根据一示例性实施方式示出损伤演化曲线及拟合结果示意图。图4是根据一示例性实施方式示出砂岩单轴压缩应力-应变曲线拟合结果示意图。图5是根据一示例性实施方式示出砂岩理论能量演化全过程的示意图。

发明人发现，耗散能的演化过程恰好代表并反映了岩石内部微缺陷及损伤程度演变规律。岩石初始压密阶段由于初始缺陷压缩闭合产生的摩擦作用会消耗部分能量，故该阶段耗散的能量反应了岩石初始损伤程度。基于此，本发明提出了一种确定岩石初始损伤程度的方法。

岩石在外界荷载作用过程中发生能量输入、可释放弹性应变能的储存以及能量的耗散释放。根据能量守恒定律，复杂应力状态下单位体积岩石的能量转化满足如下关系：

（1）/>

（2）

（3）

（4）

式中：

为岩石总共吸收的能量，单位MJ/m

为岩石变形过程中内部储存的弹性应变能，MJ/m

为加载过程中释放耗散的能量，单位MJ/m

、/>

（

为相应方向上的弹性应变；

为泊松比；

为卸载弹性模量，此处取弹性段峰值强度50%~60%处的弹性模量代替。

对于单轴压缩条件下，

（5）

（6）

（7）

式中：

按照上述计算方法，以室内砂岩单轴压缩试验结果为例，其能量演化过程如图1所示。

从图1中可以看出，岩石在受载前期（图中阶段Ⅰ）耗散能曲线先缓慢增加，后逐渐趋于稳定，这是由于岩石内部原生微缺陷、微裂纹受载压缩闭合需消耗少部分能量，当原生裂纹完全闭合时，耗散能曲线趋于稳定。岩石在裂纹完全闭合之前（图中阶段Ⅰ）消耗了少部分能量，记为

（8）

然而，根据式（8）计算得到的岩石初始损伤程度值受岩石最终破坏时刻总输入能量

岩石单轴压缩条件下理论分段本构模型如下：

（9）

式中：

为裂纹闭合应变，可根据应力-应变曲线弹性阶段反向延长线与横坐标轴的交点确定；

为裂纹闭合阶段的等效弹性模量，可根据裂纹应变-应力曲线拟合结果求得，单位MPa；

E

D

为裂纹闭合应力/>

基于变形模量衰减法获得岩石损伤变量的公式为：

（10）

式中：

为应力-应变曲线上任意一点与O＇连线的斜率；

E

根据式（10）及函数

根据式（9）建立的本构模型得到了砂岩单轴压缩应力-应变全过程理论曲线如图4所示。

根据图4中获得的砂岩理论应力-应变曲线可计算得到能量演化全过程如图5所示。根据图5可以得到砂岩最终破坏时刻总耗散能

综上所述，本发明提出了一种利用能量耗散和岩石理论应力-应变全过程曲线确定岩石初始损伤程度的方法，包括以下步骤：

步骤100：将工程现场取回的岩样按照国际岩石力学学会标准加工成直径50mm、高度100mm的标准圆柱体试样，对岩样进行单轴压缩试验，获得岩样应力-应变曲线及弹性模量E，E为岩石应力-应变曲线中弹性变形阶段曲线的斜率；

步骤200：根据已建立的理论本构模型得到岩石单轴压缩条件下应力-应变全过程曲线；

步骤300：根据步骤200中得到的理论应力-应变全过程曲线，结合上文能量计算公式可得到总能量、弹性能和耗散能演化曲线；

步骤400：根据步骤300中的耗散能演化曲线可确定岩石初始压密阶段的耗散能

故通过本发明能够深刻理解岩石变形破坏过程，以利于获得合理的岩石力学参数。

在本发明实施例中，术语“多个”则指两个或两个以上，除非另有明确的限定。术语“安装”、“连接”、“固定”等术语均应做广义理解，例如，“连接”可以是固定连接，也可以是可拆卸连接，或一体地连接。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明实施例中的具体含义。

本发明实施例的描述中，需要理解的是，术语“上”、“下”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明实施例和简化描述，而不是指示或暗示所指的装置或单元必须具有特定的方向、以特定的方位构造和操作，因此，不能理解为对本发明实施例的限制。

在本说明书的描述中，术语“一个实施例”、“一个优选实施例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或特点包含于本发明实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或实例。而且，描述的具体特征、结构、材料或特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

以上仅为本发明实施例的优选实施例而已，并不用于限制本发明实施例，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明实施例的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明实施例的保护范围之内。