+高级检索
网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

砂岩剪切破裂面形貌特征及粗糙度各向异性  PDF

  • 朱谭谭 1,2
  • 朱焕焕 1
  • 李昂 1
  • 宗西垒 1
  • 孙旭亮 1
  • 宿钟鸣 2
  • 何信 2
1. 长安大学 公路学院,陕西 西安 710064; 2. 黄土地区公路建设与养护技术交通行业重点实验室,山西 太原 030032

中图分类号: TU45

最近更新:2025-01-21

DOI: 10.16339/j.cnki.hdxbzkb.2025022

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

岩体破裂面的形态特征是影响其剪切滑移破坏的重要因素,法向应力是岩体剪切破裂面形态特征的重要影响因素.为研究法向应力对砂岩剪切破裂面形貌特征和粗糙度的影响,首先开展了砂岩双面剪切试验,获得不同法向应力作用下砂岩剪切破裂面,然后对直剪试验中形成的剪切破裂面进行三维光学扫描,最后对不同法向应力作用下砂岩剪切破裂面进行三维重构.分析了法向应力对剪切破裂面粗糙高度特征参数的影响,基于各向异性参数研究了砂岩剪切破裂面粗糙度各向异性特征,研究了砂岩剪切破裂面分形特征与粗糙度的关系.结果表明,砂岩剪切破裂面粗糙高度期望值、标准差、均值和最大起伏差都随法向应力的增大而增大;剪切破裂面粗糙度表现出显著的各向异性,随着与剪切方向夹角的增大,砂岩剪切破裂面粗糙度先减小后增大;法向应力为压应力时,剪切破裂面粗糙度在与剪切方向成75° ~ 90°时取得最小值;剪切破裂面粗糙度极差和各向异性参数随法向应力的增大逐渐增大,剪切破裂面各向异性特征增强;砂岩剪切破裂面分形维数与粗糙度呈幂函数关系.

岩体破裂面形态对破裂面的内摩擦角、黏聚力、剪胀性、峰值剪切强度和位移等具有重要影响.学者们对岩体破裂面形貌特征开展了大量研究工作.

Barton

1基于直剪试验结果,提出了考虑粗糙度的经验公式预测未填充岩石节理的抗剪强度.Barton2描述了粗糙度从0到20的10个标准剖面线,用以定量描述岩体节理面粗糙度,由于该方法简单有效,被国际岩石力学学会建议为岩石节理粗糙度评价标准方法.基于该研究成果,众多学者针对粗糙度和剪切强度之间的关系开展了大量研3-6.Reeves7提出了一种将岩石表面粗糙度参数与剪切强度联系起来的方法,该方法可以用来评估岩体破裂面粗糙度对剪切强度的影响.Özvan8选取了三种不同风化程度的花岗岩试样,评估了其表面粗糙度和风化程度对剪切强度的影响.Ban9提出了基于两个新的粗糙度参数的预估峰值剪切强度的方程,并通过试验数据验证了这种方法的可靠性.Tang10针对岩石节理形态和法向应力对其剪切强度的影响进行了研究,采用水泥砂浆材料制作了三组不同形态的岩石节理,并进行了恒定法向荷载剪切试验,研究结果表明剪切强度随节理粗糙度和法向应力的增大而增大,并提出了一种岩石节理峰值剪切强度预测准则.Belem11提出,使用三维平均倾角θs、表观各向异性度Ka、节理表面平均梯度Z2s和表面扭曲参数Ts等5个参数来描述破裂面的粗糙度,并基于这些参数预测剪切破裂面的力学行为.破裂面的物理力学性质具有显著的各向异性特征,Zhang12通过光学扫描和3D打印技术制备了与剪切岩石试样相同的混凝土试样,并沿不同方向进行剪切试验,研究发现破裂面在不同方向的剪切强度存在显著差异.Huang13用节理面正交方向的粗糙度权重比作为各向异性特征指标,提出了各向异性变化系数,推导得到了各向异性解析函数,并验证该函数能够反映各向异性分布趋势.陈世江14通过高清数码相机获得结构面粗糙度信息,应用地质统计学原理,提出了采用变异函数参数表示结构面粗糙度的方法,并应用该方法深入探讨了结构面形貌各向异性和尺寸效应特征.邹先坚15提出了一种钻孔图像结构面剖面线分析提取方法,精确还原了孔内原位岩石结构面形貌特征,并采用分形维数描述岩石结构面粗糙度特性和各向异性特征,用来判定岩体结构面的优势抗滑方向.Babanouri16通过剪切试验,研究了剪切峰值前后节理表面轮廓的破坏形式、面积和节理粗糙度的变化规律.周枝华17利用三维激光扫描技术精确测量了岩石节理表面形态,分析了节理面几何特性参数的分布规律.

目前对岩体破裂面形貌特征的研究主要针对压剪破裂面,然而,地下硐室围岩和高陡边坡等常发生拉剪破坏,岩石拉剪应力和压剪应力下破裂面破坏机制和形貌特征存在显著区别.由于拉剪试验技术上的困难,人们对拉剪应力下破裂面的粗糙度特征、分形特征、各向异性特征等的认识不足.基于此,本文首先采用自研拉剪装置开展直剪试验,制作了不同法向应力下的剪切破裂面,然后通过光学扫描对剪切破裂面进行了三维重构,研究了不同法向应力作用下剪切破裂面的形貌特征、粗糙高度特征参数、粗糙度特征及其各向异性、粗糙度和分形维数的关系等.

1 试验方案

试验采用的砂岩取自三峡库区,呈暗红色,主要矿物为石英、方解石、长石及少量磁铁矿和赤铁矿.通过切割、打磨等工序将砂岩加工成60 mm×60 mm ×120 mm的长方体试样.试验系统及辅助装置如图1所示.为研究砂岩剪切破裂面形貌特征,首先通过双面剪切试验得到不同法向应力条件下的剪切破裂面,然后采用Cronos三维光学扫描仪对剪切破裂面扫描重构.试验开始前,测得试样的平均抗拉强度为1.61 MPa.因此,在开展双面剪切试验时,法向拉应力设置-0.25 MPa、-0.5 MPa、-0.75 MPa、-1.0 MPa、-1.25 MPa和-1.5 MPa六个拉应力水平.为对比分析法向拉应力和法向压应力条件下剪切破裂面各向异性特征,设置9 MPa、6 MPa、3 MPa三个法向压应力水平.双面剪切试验工况及产生的剪切破裂面如表1所示.

fig

图1  试验系统及辅助装置

Fig.1  Test system and auxiliary devices

表1  双面剪切试验工况及对应的剪切破裂面
Tab. 1  Double-sided shear test conditions and shear fracture surfaces
试样法向应力/MPa破裂面编号试样法向应力/MPa破裂面编号
1# -1.50 T-1.5-1 6# -0.25 T-0.25-1
T-1.5-2 T-0.25-2
2# -1.25 T-1.25-1 7# 0 0-1
T-1.25-2 0-2
3# -1.00 T-1.0-1 8# 3.00 P-3-1
T-1.0-2 P-3-2
4# -0.75 T-0.75-1 9# 6.00 P-6-1
T-0.75-2 P-6-2
5# -0.50 T-0.5-1 10# 9.00 P-9-1
T-0.5-2 P-9-2

双面剪切试验过程中,首先施加法向应力至设定值,加载方式为应力控制模式,加载速率为0.05 kN/s,然后施加剪应力直至砂岩试样发生破坏.剪应力加载采用位移控制模式,加载速率为0.2 mm/min. Cronos三维光学扫描得到的剪切破裂面点云为无序点云,通过MATLAB编程对点云数据进行处理可以得到剪切破裂面有序点云数据.

图1中剪切破裂试样所示,在剪切试验过程中,试样除发生沿剪切面的剪切破坏,在试样表面也会发生局部剥落破坏.为了消除试样表面剥落对计算结果的影响,在计算剪切破裂面参数时,选择剪切破裂面中间区域作为研究范围.如图2所示,破裂面边缘各删除 2 mm宽度,即研究范围为56 mm×56 mm的正方形.

fig

图2  试样表面剥落及剪切破裂面研究范围

Fig.2  Sample surface spalling and the research scope of shear fracture surface

2 试验结果与分析

2.1 剪切破裂面形貌特征

基于光学扫描得到的剪切破裂面点云数据,将破裂面高程最低点所在XOY平面(平行于试样端面)高程坐标设置为零,对不同法向应力作用下砂岩剪切破裂面进行三维重构,得到破裂面形貌特征如图3所示.由于篇幅原因,图3中相同法向应力下仅呈现了其中一个剪切破裂面.剪切破裂面最大高程范围为4.84~12.43 mm,最大值为12.43 mm,在法向应力为6 MPa时取得;最小值为4.84 mm,在法向应力为-0.25 MPa时取得.

fig

(a) σn=-1.5 MPa

(b) σn=-1.25 MPa

  

fig

(c) σn=-1.0 MPa

(d) σn=-0.75 MPa

  

fig

(e) σn=-0.5 MPa

(f) σn=-0.25 MPa

  

fig

(g) σn=0 MPa

(h) σn=3 MPa

  

fig

(i) σn=6 MPa

(j) σn=9 MPa

  

图3  剪切破裂面三维重构

Fig.3  3D reconstruction of shear fracture surface

粗糙高度是指破裂面上某高点距离基准面的高度,剪切破裂面起伏高度相对频率可以反映破裂面粗糙高度的分布规律.以z=0 mm对应的平面作为基准面,计算扫描点的高程分布频度.粗糙高度期望值μ、标准差Sq、均值Sa、最大起伏差Sp等粗糙高度特征参数能够定量表征剪切形成的破裂面的形态特征.其中,粗糙高度最大起伏差Sp是剪切破面上最高点到计算参考面的垂直距离.砂岩试样剪切破裂面粗糙高度特征参数与法向应力关系如图4所示.

fig

(a) 粗糙高度期望值

(b) 粗糙高度标准差

  

fig

(c) 粗糙高度均值

(d) 粗糙高度最大起伏差

  

图4  剪切破裂面粗糙高度特征参数

Fig.4  Variations of the roughness height indexes for the shear fracture surfaces

图4可知,随着法向应力的增大,粗糙高度期望值、最大起伏差、标准差和均值均表现为逐渐增大.这是因为除局部锯齿状起伏会影响破裂面粗糙高度外,破裂面整体起伏同样会对粗糙高度产生影响.随着法向应力的增大,破裂面局部锯齿状起伏减小,但破裂面整体起伏增大,导致粗糙高度随法向应力的增大逐渐增大.黄达

18通过砂岩剪切试验同样发现法向应力越大,剪切破裂面整体起伏越大;法向应力越小,剪切破裂面越平直,这与本试验结果相符.当法向应力为拉应力时,粗糙高度特征参数随法向应力的变化离散性较强,波动范围较大,规律性不明显.当法向应力为压应力时,粗糙高度特征参数与法向应力的相关性更为明显,粗糙度特征参数波动范围减小.

2.2 粗糙度各向异性

试样剪切破坏后,假定采用一个平行于法相应力方向的平面将试样切开,每个破裂面与该平面有一条交线,该交线即为破裂面的剖面线.Barton

2给出了10个剖面线来表示破裂面的粗糙度,该方法被国际岩石力学与岩石工程学会推荐作为岩体破裂面粗糙度标准曲线.

Tse

19对剖面线进行了数字化,提出了岩体破裂面粗糙度和破裂面等高线之间的关系方程:

JRC=32.2+32.47logZ2 (1)
Z2=iN-1Zi+1-Zi2N-1Δs2 (2)

式中:Z2为剖面线一阶导数均方根,Zi为破裂面采样点的高程坐标,Δs为剖面线上数据采样间隔距离(如图5所示).标准剖面线和由其得到的粗糙度可以很好地描述曲线的粗糙特征.然而,岩体破裂面是一个三维曲面,在破裂面上沿不同角度可以得到无数条剖面线,且各剖面线的粗糙度不同,因此,单一标准剖面线法不能描述岩体破裂面粗糙度的各向异性特征.为研究岩石剪切破裂面粗糙度各向异性特征,基于剪切破裂面三维扫描得到的点云数据,沿某一方向描绘出多条剪切破裂面剖面线,根据式(1)式(2)即可计算出每条剖面线的粗糙度,以某一方向所有剖面线粗糙度的平均值作为剪切破裂面在该方向上的粗糙度:

JRCα=1ni=1nJRCαi (3)
fig

图5  剪切破裂面粗糙度各向异性计算原理

Fig.5  Principle for anisotropy calculation of roughness

式中:JRCα为剪切破裂面在某方向上的粗糙度,JRCαi为在某方向上第i条剖面线的粗糙度.岩石剪切破裂面粗糙度各向异性计算原理如图5所示.其中α为剖面线在水平面上的投影与剪切方向之间的倾角,d为相邻两条剖面线之间的距离,∆s为剪切破裂面上数据点采样间隔,Zαmn为第m行和第n列测点的实测高度.当倾角不等于0°时,采用线性插值法得到目标位置的高度.

不同法向应力作用下,剪切破裂面粗糙度随α的变化规律如图6所示.在不同法向应力下,粗糙度整体上表现为先减小后增大.当法向应力大于0时,粗糙度基本呈轴对称分布,其对称轴的范围为75°~90°;当法向应力小于0时,粗糙度虽然随着倾角α的增大先减小后增大,但粗糙度并不呈轴对称分布.

fig

  

fig

  

fig

  

fig

  

fig

  

图6  剪切破裂面粗糙度随倾角α的变化规律

Fig.6  Changes in the JRC of shear fracture surfaces with the dip angle α

粗糙度极差为砂岩剪切破裂面不同方向上最大粗糙度和最小粗糙度的差值,可以用来定量描述岩体破裂面在不同方向上粗糙度差异.不同法向应力作用下,砂岩剪切破裂面粗糙度极差变化规律如图7(a)所示.随着法向应力的增大,砂岩剪切破裂面粗糙度极差整体上表现为逐渐增大,破裂面各向异性特征增强,特别是当法向应力为压应力时,各向异性随法向应力的增大而增大这一特征更为显著.剪切破裂面粗糙度极差与法向应力之间的关系可用线性方程表示:

JRCR=0.45σn+3.37 (4)
fig

图7  粗糙度各向异性指标变化规律

Fig.7  Variation of roughness anisotropy index

(a)粗糙度极差 (b)各向异性参数

各向异性参数DAC可以用来定量描述破裂面粗糙度的各向异性程

20,其取值范围为0~1,DAC越大,破裂面各向异性特征越显著;DAC越小,各向异性越不明显.各向异性参数可用下式计算:

DAC=1-1eCV (5)
CV=1n-1α=0nJRCα-JRC¯2JRC¯ (6)
JRC¯=1nα=0nJRCα (7)

式中:JRCα是倾角为α时破裂面的粗糙度,JRC¯为JRCα的平均值,CV为破裂面统计方向上粗糙度的变异系数.各向异性参数与法向应力的关系如图7(b)所示.由图可知,砂岩剪切破裂面各向异性参数随着法向应力的增大整体上呈逐渐增大趋势,这表示破裂面各向异性特征随法向应力的增大逐渐增强.各向异性参数与法向应力的关系可用式(8)表示:

DAC=0.056σn+1.910.32 (8)

2.3 分形特征

Xie

21提出了一种使用高阶凸起的平均基长和平均高度来定义分形维数的方法,分形维数可用式(9)计算:

D=log4log21+cosarctan2HL (9)
H=1Ni=1Nhi (10)
L=1Ni=1Nli (11)

其中:D为分形维数,hili分别为剪切破裂面剖面线第i个斜坡的高度和长度.

根据光学扫描得到的点云数据和式(9)可以计算出每条剖面线的分形维数.不同法向应力下砂岩试样剪切破裂面分形维数与粗糙度的关系如图8所示.在不同的法向应力下,剪切破裂面的分形维数呈非线性增加.当剖面线粗糙度趋近于0时,分形维数趋近于1,这是因为粗糙度JRC=0时,剪切破裂面是一个光滑平面,式(9)中斜坡平均高度H=0,D=1.

fig

  

fig

  

fig

  

fig

  

fig

  

图8  粗糙度与分形维数的关系

Fig. 8  Relationship between roughness and fractal dimension

图8所示,粗糙度和分形维数之间的关系可以用式(12)来表示:

JRC=a(D-1)b (12)

式中:ab是拟合参数.法向应力在-1.5~9 MPa范围内,参数a的取值范围为65.90~82.57,参数b的取值范围为0.43~0.49.

图9(a)所示,参数ab存在强线性关系:

a=281.47b-55.13 (13)

式(13)代入式(12),砂岩剪切破裂面粗糙度与分形维数的关系可以改写为:

JRC=(281.47b-55.13)(D-1)b (14)

参数ab随法向应力的增大整体表现为逐渐增大.如图9(b)所示,由于砂岩材料是非连续非均质材料,当法向应力为-1 MPa和-0.25 MPa时,两个参数明显偏离了数据的整体变化趋势,可认为这两组数据异常,在曲线拟合时忽略这两组数据.参数ab与法向应力的关系可用对数函数表示:

a=69.78+3.59lnσn+2.5 (15)
b=0.44+0.014ln(σn+2.5) (16)
fig

(a)  参数ab关系曲线

fig

(b)  参数ab与法向应力关系

图9  拟合参数ab变化规律

Fig.9  Changes in fitting parameters a and b

3 结 论

为研究砂岩剪切破裂面形貌特征和粗糙度各向异性特征,开展了不同法向应力作用下砂岩剪切试验,并对剪切破裂面进行了三维扫描和重构,研究了砂岩剪切破裂面的形貌特征、粗糙度特征及其各向异性、粗糙度和分形维数关系等,可以得到以下主要结论:

1)砂岩试样剪切破裂面粗糙高度期望值、最大起伏差、标准差和均值随法向应力的增大整体逐渐增大.

2)剪切破裂面的粗糙度具有显著的各向异性特征.随着与剪切方向夹角的增大,剪切破裂面粗糙度先减小后增大.法向应力为压应力时,粗糙度基本呈轴对称分布,对称轴的范围为75°~90°.

3)粗糙度与分形维数之间的关系可用JRC=

aD-1b描述,其中参数ab随法向应力的增大呈非线性增大,ab呈强线性关系.

参考文献

1

BARTON NReview of a new shear-strength criterion for rock joints[J].Engineering Geology197374): 287-332 [百度学术] 

2

BARTON NCHOUBEY VThe shear strength of rock joints in theory and practice[J]. Rock Mechanics1977101): 1-54 [百度学术] 

3

LI BMO Y YZOU L Cet al. Influence of surface roughness on fluid flow and solute transport through 3D crossed rock fractures[J].Journal of Hydrology2020582124284 [百度学术] 

4

MA S MSINGER GCHEN S Het al. Objective-driven solid-surface-roughness characterization for enhanced nuclear-magnetic-resonance petrophysics[J]. SPE Journal2021265):2860-2879 [百度学术] 

5

金磊磊魏玉峰黄鑫基于节理面三维形貌的岩石节理抗剪强度计算模型[J].岩土力学20204110): 3355-3364 [百度学术] 

JIN L LWEI Y FHUANG Xet alShear strength calculation model of rock joints based on three-dimensional morphology of joint surface[J].Rock and Soil Mechanics20204110):3355-3364(in Chinese) [百度学术] 

6

YANG Z YLO S CDI C CReassessing the joint roughness coefficient (JRC) estimation using Z2[J].Rock Mechanics and Rock Engineering2001343): 243-251 [百度学术] 

7

REEVES M JRock surface roughness and frictional strength[J].International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts1985226): 429-442 [百度学术] 

8

ÖZVAN ADINÇER İACAR Aet alThe effects of discontinuity surface roughness on the shear strength of weathered granite joints[J].Bulletin of Engineering Geology and the Environment2014733): 801-813 [百度学术] 

9

BAN L RQI C ZLU C SA direction-dependent shear strength criterion for rock joints with two new roughness parameters[J].Arabian Journal of Geosciences20181116): 466 [百度学术] 

10

TANG Z CLIU Q SHUANG J HNew criterion for rock joints based on three-dimensional roughness parameters[J].Journal of Central South University20142112): 4653-4659 [百度学术] 

11

BELEM THOMAND-ETIENNE FSOULEY MQuantitative parameters for rock joint surface roughness[J].Rock Mechanics and Rock Engineering2000334):217-242 [百度学术] 

12

ZHANG Q ZWU C ZJANG B Aet al. Shear behavior of greenschist along foliation plane considering anisotropy[J]. International Journal of Earth Sciences20191085): 1755-1761 [百度学术] 

13

HUANG MHONG C JMA C Ret alCharacterization of rock joint surface anisotropy considering the contribution ratios of undulations in different directions[J].Scientific Reports2020101): 17117 [百度学术] 

14

陈世江朱万成刘树新岩体结构面粗糙度各向异性特征及尺寸效应分析[J].岩石力学与工程学报2015341):57-66 [百度学术] 

CHEN S JZHU W CLIU S Xet al. Anisotropy and size effects of surface roughness of rock joints[J]. Chinese Journal of Rock Mechanics and Engineering2015341): 57-66(in Chinese) [百度学术] 

15

邹先坚王益腾王川婴. 钻孔图像中岩石结构面三维形貌特征及优势抗滑方向研究[J]. 岩土力学202041增刊1): 290-298. [百度学术] 

ZOU X JWANG Y TWANG C Y. Three-dimensional morphological feature of rock structural surfaces and dominant anti-slip direction using the borehole images[J]. Rock and Soil Mechanics202041Sup.1): 290-298. (in Chinese) [百度学术] 

16

BABANOURI NKARIMI NASAB SModeling spatial structure of rock fracture surfaces before and after shear test:a method for estimating morphology of damaged zones[J].Rock Mechanics and Rock Engineering2015483): 1051-1065 [百度学术] 

17

周枝华杜守继岩石节理表面几何特性的三维统计分析[J].岩土力学2005268):1227-1232 [百度学术] 

ZHOU Z HDU S J. 3D statistic analysis of geometrical properties of a rock joint[J]. Rock and Soil Mechanics2005268): 1227-1232(in Chinese) [百度学术] 

18

黄达张永发朱谭谭. 砂岩拉-剪力学特性试验研究[J].岩土工程学报2019412): 272-276 [百度学术] 

HUANG DZHANG Y FZHU T Tet al. Experimental study on tension-shear mechanical behavior of sandstone[J]. Chinese Journal of Geotechnical Engineering2019412): 272-276(in Chinese) [百度学术] 

19

TSE RCRUDEN D M. Estimating joint roughness coefficients[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts1979165): 303-307 [百度学术] 

20

宋磊博江权李元辉不同采样间隔下结构面形貌特征和各向异性特征的统计参数稳定性研究[J].岩土力学2017384): 1121-1132 [百度学术] 

SONG L BJIANG QLI Y Het alStability of statistical parameters of discontinuities morphology and anisotropy based on different sampling intervals[J].Rock and Soil Mechanics2017384): 1121-1132(in Chinese) [百度学术] 

21

XIE H PPARISEAU W G. Fractai estimation of joint roughness coefficients[J]. Science in China (Series B)19943712): 1516-1524. [百度学术] 

作者稿件一经被我刊录用,如无特别声明,即视作同意授予我刊论文整体的全部复制传播的权利,包括但不限于复制权、发行权、信息网络传播权、广播权、表演权、翻译权、汇编权、改编权等著作使用权转让给我刊,我刊有权根据工作需要,允许合作的数据库、新媒体平台及其他数字平台进行数字传播和国际传播等。特此声明。
关闭