基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法技术

本发明专利技术公开了基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法，具体为：收集深部各向异性储层岩石的物理力学参数、工程地质环境参数和双边缺口巴西圆盘几何参数；根据几何参数，绘制模拟不同层理角度真实剪切断裂的岩体计算域；对计算域添加固体力学场、历史应变场以及相场等多个物理场，建立相互耦合的各向异性真实剪切相场模型，并采用交错方案求解耦合，保持位移加载速率不变，重复上述步骤，直至岩石完全剪切断裂。该方法充分考虑岩石的各向异性特征以及工程地质环境特征，定量地表征了各向异性岩石真实剪切断裂模式，为评价深部地下工程施工过程中岩石的力学参数提供了重要参考。供了重要参考。供了重要参考。

【技术实现步骤摘要】
基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法


[0001]本专利技术属于岩土工程
，具体涉及基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法。

技术介绍

[0002]随着经济的不断发展，地下工程建设逐渐向更大深度扩展，高围压下的工程地质环境导致储层岩石更容易发生剪切破坏。例如，在非常规油气钻探和深部隧道开挖过程中，井壁以及洞室附近的围岩容易因应力释放而发生剪切破坏和失稳。同时，由于大部分地层是由不同岩性岩石组成的层状结构，储层岩石普遍具有各向异性特征，造成剪切裂缝扩展行为复杂。因此，了解各向异性岩石的剪切破裂机制对于地下能源的有效开发和地下岩石工程的安全至关重要。
[0003]大量文献表明，各向异性岩石的剪切破裂机制在很大程度上尚未被探索，主要是由于难以进行基于剪切的断裂试验。现有标准室内剪切裂纹扩展试验中通常会观察到裂纹出现一定角度的扭结，这种扭结裂纹的形成是由拉应力驱动的。因此，大部分实验无法准确了解岩石的真实剪切破裂机制。尽管已有报道提出的穿通剪切(PTS)试验和剪切箱试验可以捕捉到岩石的真实剪切断裂形态，但由于其装置的复杂型和额外要求，仅能获取有限的实验数据。因此，受限于室内试验的苛刻条件，开展各向异性岩石真实剪切断裂扩展的数值模拟是有必要的。
[0004]双边缺口巴西圆盘(DNBD)的径向压缩作为一种测试真实剪切断裂的新兴方法，不仅操作简单，还可以确定任何方向上相对于各向异性方向的真实剪切断裂韧性。然而，目前基于双边缺口巴西圆盘的剪切断裂扩展模型尚未考虑岩石的各向异性和深部工程地质环境对岩石真实剪切裂纹扩展所带来的影响，存在一定局限性。

技术实现思路

[0005]本专利技术的目的是提供基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法，综合考虑了深部各向异性储层岩石的工程环境特征，从而准确预测岩石真实剪切断裂形态。
[0006]本专利技术所采用的技术方案是，基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法，具体按照以下步骤实施：
[0007]步骤1、收集深部各向异性储层岩石的物理力学参数、地下工程地质环境参数和双边缺口巴西圆盘几何参数；
[0008]步骤2、根据物理力学参数中获得的不同主方向下的剪切断裂韧性，计算出不同层理角度下岩石纯剪切断裂的临界能量释放率G
c
，公式如式(1)及式(2)所示：
[0009]K
IIc
＝K
IIc,1 sin2α1+K
IIc,2 cos2α1ꢀꢀꢀꢀꢀ
(1)；
[0010][0011]式中，K
IIc,1
为垂直于层理方向剪切断裂韧性；K
IIc,2
为沿层理方向剪切断裂韧性；E1为垂直于层理方向弹性模量；ζ＝E1/E2，E2为沿层理方向弹性模量；α1为层理角度；v
12
和v
23
分别为沿层理平面内和垂直层理平面内的泊松比；
[0012]步骤3、利用COMSOL软件绘制双边缺口巴西圆盘几何计算域，并利用COMSOL内置基矢坐标系表征岩石不同层理角度下的横观各向同性特征；
[0013]步骤4、对计算域分别添加固体力学场、历史应变场以及相场，建立相互之间的耦合关系，从而获得各向异性真实剪切相场模型；
[0014]步骤5、对各向异性真实剪切相场模型采用交错格式在分离耦合求解器上独立求解；
[0015]步骤6、保持位移加载速率不变，重复步骤5，直至岩石完全断裂后停止计算，即可获得不同层理角度时岩石的真实剪切断裂形态。
[0016]本专利技术的特点还在于，
[0017]步骤1中，物理力学参数包括：岩石沿层理方向和垂直层理方向的弹性模量、剪切模量、泊松比、剪切断裂韧性以及岩石的比奥系数；工程地质环境参数包括：岩石所受围压、孔隙压力以及原位应力；双边缺口巴西圆盘几何参数包括：圆盘的半径、预制缺口尺寸、加载角度。
[0018]步骤3中，绘制双边缺口巴西圆盘几何计算域的具体步骤为：首先以坐标原点为中心生成半径R为37.5mm的圆形以及94mm
×
1.8mm的矩形，随后以坐标形式建立两条插值曲线用以确定预制缺口长度，将所有几何对象通过布尔并集合并后删除两侧预制缺口区域，最后获得双边缺口巴西圆盘模型计算域；
[0019]建立内置基矢坐标系，如式(3)所示：
[0020][0021]式中，x和y分别表示全局坐标系x轴方向和y轴方向；x
′
和y
′
分别表示垂直于层理方向和沿层理方向；α表示沿层理方向与y轴方向的夹角，分别选取五个不同层理角度α＝0
°
、22.5
°
、45
°
、67.5
°
以及90
°
用以真实剪切模拟分析。
[0022]步骤4中，各向异性真实剪切相场模型包括固体力学方程、相场控制方程和历史应变方程；
[0023]其中，固体力学方程如式(4)所示：
[0024][0025]式中，σ＝[(1
‑
κ)(1
‑
φ)2+κ]￡:ε，0＜κ≤1用以避免数值模型的奇异性，φ为裂纹
萌生以及扩展的相场，￡为未退化的四阶刚度张量，ε为应变张量；f
b
和f
s
分别为岩石的体力和面力；n为方向矢量；
▽
为梯度算子；
[0026]相场控制方程如式(5)所示：
[0027][0028]式中，l0为长度尺寸参数；G
c
为临界能量释放率；B为各向异性梯度；弹性应变能的拉伸部分
[0029]其中，ε
ii
为线应变，ε
ij
(i≠j)为切应变；i＝1,2,3、j＝1,2,3；
[0030]▽
φ为相场梯度；C
l
为未退化的四阶刚度张量￡的Voight表示法，其表示为：
[0031][0032]式中，为刚度分量，其不同表达式如下：
[0033][0034][0035]式中，E
i
(i＝1,2,3)表示沿三个主要材料方向的弹性模量(e1,e2,e3)；G
ij
(i＝1,2,3、j＝1,2,3)是由e
i
和e
j
定义的正交平面内的剪切模量；v
ij
＝v
ji
E
i
/E
j
为正交平面中的泊松比；
[0036]历史应变方程如式(6)所示：
[0037][0038]式中，t为时间；x为二维域任意点的位置矢量；σ0为初始应力张量。
[0039]步骤5中，具体为：
[0040]将步骤1的深部各向异性储层岩石的物理力学参数、地下工程地质环境参数代入
各向异性真实剪切相场模型中，依据步骤3给定模型不同的层理角度，将圆本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法，其特征在于，具体为：收集深部各向异性储层岩石的物理力学参数、工程地质环境参数和双边缺口巴西圆盘几何参数；根据几何参数，绘制模拟不同层理角度真实剪切断裂的岩体计算域；对计算域添加多个物理场，建立相互耦合的各向异性真实剪切相场模型；采用交错格式在分离耦合求解器上独立求解；保持位移加载速率不变，重复上述步骤，直至岩石完全剪切断裂，即可获得不同层理角度时岩石的真实剪切断裂形态。2.根据权利要求1所述的基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法，其特征在于，具体按照以下步骤实施：步骤1、收集深部各向异性储层岩石的物理力学参数、地下工程地质环境参数和双边缺口巴西圆盘几何参数；步骤2、根据物理力学参数中获得的不同主方向下的剪切断裂韧性，计算出不同层理角度下岩石纯剪切断裂的临界能量释放率G
c
，公式如式(1)及式(2)所示：K
IIc
＝K
IIc,1 sin2α1+K
IIc,2 cos2α1ꢀꢀꢀꢀꢀ
(1)；式中，K
IIc,1
为垂直于层理方向剪切断裂韧性；K
IIc,2
为沿层理方向剪切断裂韧性；E1为垂直于层理方向弹性模量；ζ＝E1/E2，E2为沿层理方向弹性模量；α1为层理角度；v
12
和v
23
分别为沿层理平面内和垂直层理平面内的泊松比；步骤3、利用COMSOL软件绘制双边缺口巴西圆盘几何计算域，并利用COMSOL内置基矢坐标系表征岩石不同层理角度下的横观各向同性特征；步骤4、对计算域分别添加固体力学场、历史应变场以及相场，建立相互之间的耦合关系，从而获得各向异性真实剪切相场模型；步骤5、对各向异性真实剪切相场模型采用交错格式在分离耦合求解器上独立求解；步骤6、保持位移加载速率不变，重复步骤5，直至岩石完全断裂后停止计算，即可获得不同层理角度时岩石的真实剪切断裂形态。3.根据权利要求2所述的基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法，其特征在于，所述步骤1中，物理力学参数包括：岩石沿层理方向和垂直层理方向的弹性模量、剪切模量、泊松比、剪切断裂韧性以及岩石的比奥系数；工程地质环境参数包括：岩石所受围压、孔隙压力以及原位应力；双边缺口巴西圆盘几何参数包括：圆盘的半径、预制缺口尺寸、加载角度。4.根据权利要求2所述的基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法，其特征在于，所述步骤3中，绘制双边缺口巴西圆盘几何计算域的具体步骤为：首先以坐标原点为中心生成半径R为37.5mm的圆形以及94mm
×
1.8mm的矩形，随后以坐标形式建立两条插值曲线用以确定预制缺口长度，将所有几何对象通过布尔并集合并后删除两侧预制缺口区域，最后获得双边缺口巴西圆盘模型计算域；建立内置基矢坐标系，如式(3)所示：...

【专利技术属性】
技术研发人员：薛熠，梁鑫，杨博鹍，时旭阳，张云，滕腾，苏善杰，刘建强，周阳，
申请(专利权)人：西安理工大学，
类型：发明
国别省市：