Structural Surface Roughness Based on UAV High Density Point Cloud Fractal Characteristics and Anisotropy
-
摘要: 为研究岩体结构面各向异性对粗糙度评价的影响，以藏东南某铁路察达工点高陡斜坡为研究对象，运用无人机综合摄影测量技术，提取研究区结构面高密度点云并剪裁结构面轮廓线，采用修正直边法与盒维数法求算粗糙度系数JRC与分形维数D，拟合JRC与D的新公式并利用数字化Barton标准线验证.选取压剪性和拉张性结构面各15个，运用新公式计算各采样方向的JRC.结果表明：压剪性结构面粗糙度各向异性规律显著，整体上JRC由剪切滑动方向至垂直剪切滑动方向递增，呈椭圆状分布；拉张性结构面粗糙度存在各向异性但无明显规律，JRC随采样角度变化波动较大，呈刺状分布.证明不同力学成因的结构面JRC各向异性存在差异，在评价粗糙度时应遵循不同采样规则.Abstract: To study the influence of anisotropy of rock mass structural plane on roughness evaluation, taking the rock mass on the high and steep slope of a railway Chada site in southeast Tibet as the research object, the integrated photogram technology of unmanned aerial vehicle was used to extract the high-density point cloud of the structural plane in the study area and cut the outline of the structural plane. The roughness coefficient JRC and fractal dimension D were calculated by the modified straight edge method and the box dimension method respectively. The new JRC and D formulas are fitted and verified by digital Barton standard lines. 15 structural planes formed under shear and tension were selected and the new formula was used to calculate the JRC in large quantities. The results show that the shear surface roughness is anisotropic, and the JRC is smaller along the shear direction, larger along the vertical shear direction, and the fluctuation with sampling angle is small. The tensile surface roughness also has anisotropy but no obvious rule, and JRC fluctuates greatly with sampling angle. It is proved that the JRC anisotropy of structural planes with different mechanical origin is also different, and different sampling rules should be followed when evaluating roughness.
-
0. 引言
结构面的存在会破坏岩体完整性并导致其强度降低（谷德振，1979），隧道、边坡等岩体工程的破坏大都是沿结构面发生的（张志飞等，2024），岩体结构面的研究是实际工程中分析岩体稳定性的基础工作（柴波等，2022）.其中，结构面粗糙度直接影响结构面的抗剪强度（王思敬，1990）.由于结构面的形貌特征存在显著的不均匀性与各向异性，导致其粗糙度也存在各向异性.因此，能够快速、大批量获取结构面信息并开展粗糙度各向异性研究能够节约工程时间、确保工程安全性，具有重要的工程应用价值.
结构面信息的获取方法主要有两种：接触式采集与非接触式采集.接触式采集主要运用针状轮廓仪（Barton and Choubey，1977）、简易纵剖面积仪（杜时贵，1992）等仪器.这些仪器成本低、便于携带和使用，但误差较大且测量范围较小；非接触式采集包括近景摄影测量法（王凤艳等，2012）、三维激光扫描法（葛云峰等，2017；罗泽军等，2022）等，这些方法具有较高的测量精度，但存在测量盲区和易受地形影响等缺点.无人机摄影测量技术日渐成熟，其具有操作简便、效率高、成本低、可快速建立三维模型（Zhao et al.，2023）等优势，在结构面形貌特征获取中可弥补传统方法的不足（Zhao et al.，2023），显著提升结构面信息采集的准确性和效率（杨泽等，2023；潘晓娟等，2024）.同时，将无人机摄影测量技术运用到岩体结构面信息获取上（崔溦等，2023），可实现大批量、大范围的采集（Song et al.，2024）.
Barton（1973）首次提出结构面粗糙度系数JRC（joint rounghness coefficient），用于表示结构面粗糙起伏程度对其抗剪强度的影响，并通过大量试验给出了10条结构面标准剖面线.此后，与这10条标准线进行对比的JRC确定方法被广泛运用在工程实例中，但这一方法存在较强的人为主观性.目前针对结构面粗糙度的评价方法主要分为二维和三维两种，其中二维方法主要有直边图解法（Barton and Bandis，1982）、修正直边法（杜时贵等，1996）、起伏度参数法（EL-Soudani，1978；Tse and Cruden，1979；Belem et al.，2000）和分形维数法（Mandelnorot，1982）等.分形维数法是用于描述一些高度不规则集合的方法，常用的二维曲线分形维数的计算方法有码尺法（杜时贵，1997）、节理分形模型h-L法（谢和平，1995）、变异函数分析法（Kulatilake et al.，1997）等，这些方法计算准确，但原理相对复杂；此外，标准差纵横比法、均方根法（Kulatilake et al.，1999）和基于数字处理的分形法（陈世江等，2012）等计算过程虽然相对复杂，但都为结构面粗糙度的分形计算提供了不同思路.相比于上述方法，盒维数法（Walsh and Watterson，1993）具有原理简单、结果准确、便于编程使用的优点，方便推广应用.
结构面粗糙度分析不仅局限于评价方法的研究，同一结构面在不同采样方向下粗糙度也会存在差异，即各向异性，亦是研究重点.杜时贵和唐辉明（1993）利用统计学原理，得到了节理和断层JRC的各向异性规律.Xie and Wang（1999）利用分形维数法描述了结构面的各向异性特征.游志诚等（2014）研究了不同粗糙度结构面抗剪强度参数的各向异性.葛云峰等（2016）改进了Grasselli法，通过计算结构面粗糙度的θmax*/（C+1）参数，研究了结构面的各向异性.洪陈杰等（2020）结合各向异性变异系数AVC3D对结构面粗糙度的各向异性进行了研究.王述红等（2023）通过不同采样方法，对结构面抗剪强度和粗糙度的各向异性进行了分析.以上研究均得出结构面粗糙度存在明显的各向异性，且采样方向的选择显著影响粗糙度评价结果，但对野外不同力学成因的天然结构面粗糙度各向异性尚未开展深入研究，对天然结构面粗糙度评价时采样方向如何选取也未给出科学建议.
为此，本文选取藏东南某铁路察达工点高陡岩质斜坡为研究区，运用无人机综合摄影测量技术，求算结构面轮廓线的分形维数D，并拟合出结构面粗糙度系数JRC与D的新公式，用于计算由压剪作用和拉张作用形成的结构面粗糙度系数，分析其粗糙度各向异性及分布规律，并给出野外实测结构面粗糙度时采样方向选取的建议.
1. 研究区概况
研究区藏东南某铁路察达工点，位于西藏自治区昌都市洛隆县腊久乡（图 1a），西邻念青唐古拉山脉东段，位于横断山脉中段，山脉走向为北西-南东向，地势险峻，地形起伏大.研究区发育怒江支流冻错曲，流向为自南向北，地势切割较深，形成冰川“U”型谷，整体属于高山峡谷地貌（图 1b）.该地区年平均温度7.6 ℃，年降水量约680 mm.
图 1 研究区地理位置、地形地貌与斜坡全貌a.研究区地理位置图；b.研究区地形地貌图；c.研究区高陡斜坡全貌图Fig. 1. Geographical location, landform and slope overall appearance of the study area研究区域位于冈底斯念青唐古拉陆块中的昂龙岗日-班戈-腾冲岩浆弧带.历经复杂的构造变质与变形作用，发育有丰富的褶皱、逆冲、推覆、伸展等断层结构.研究区内发育义俄-巴曲断层，为逆断层，产状为10°∠60°，在第四系地层未发现断裂活动迹象（郭长宝等，2024）.
研究区察达高陡岩质斜坡的坡顶和坡脚高程分别为4 065 m、3 683 m，相对高差382 m，斜坡上窄下宽，平均宽度达435 m.区内坡顶部位为古近系始新统宗白群紫红色砂砾岩，斜坡中下部分布石炭系上统二叠系下统来姑组深灰色砂板岩，斜坡底部局部有花岗闪长岩侵入，岩体结构主要为板状结构，岩体完整性为较破碎.斜坡两侧发育冲沟，坡脚附近崩积物主要有砂、泥和棱角状的块石等，块石直径大小不一，最大可达数米.斜坡中部发育两条深大断裂的次级小型逆断层，在坡面上延伸长度均超过200 m.受断裂活动影响，斜坡内各类结构面十分发育（图 1c）.
2. 点云数据采集与处理
2.1 结构面高密度点云获取
本研究为提取研究区高陡斜坡高密度点云，运用以无人机仿地飞行模式下的倾斜摄影测量为基础，无人机多角度贴近摄影测量为补充的综合摄影测量技术.倾斜摄影测量技术能够有效解决因地形高差巨大产生的影像分辨率波动大、影像模糊等问题，但由于高陡斜坡坡体高大复杂，在复杂块体及反倾面等位置仍会存在摄影盲区（杨泽等，2023）.为此，本研究在实际航摄中使用无人机多角度贴近摄影测量技术作为补充，具体步骤如下：（1）场地勘测.根据实际地形条件，选择无人机起降点并预估拍摄范围，为后续无人机航线设计提供重要参考；（2）低分辨率三维地形采集.运用无人机仿地飞行模式下的倾斜摄影技术.输入无人机拍摄范围，设置飞行模式为仿地模式，输入数字高程数据（DEM），确定方位角，设置无人机飞行速度，控制航向重叠率80%及以上，旁向重叠率60%及以上，进行无人机倾斜摄影，并建立三维模型；（3）多角度航线设计.根据研究区地形特征拟合飞行平面，并结合现场调查所得结构面发育情况在飞行平面上规划无人机多角度拍摄航线，基于步骤（2）建立的低分辨率三维模型进行检查并设计拍摄角度，进行多角度拍摄；（4）全景高分辨率图像采集.针对高陡斜坡重点区域复杂块体，基于大疆经纬M300 RTK旋翼无人机的高精度定位能力和较强的姿态控制能力，搭载禅思P1相机（详细参数见表 1），通过在现场控制无人机靠近坡面（距坡面约30~50 m）并调整相机拍摄角度，从不同位置、不同方向对高陡斜坡重点区域或关键块体进行多角度、全景式的综合摄影测量（图 2）.
表 1 无人机和相机相关参数Table Supplementary Table UAV and camera related parameters大疆M300-RTK无人机 禅思P1相机 无人机总重(kg) 6.3 重量(g) 800 最大飞行时长(min) 55 相机尺寸(mm) 198×166×129 最大可承受风速(m/s) 15 传感器面积(mm2) 35.9×24 RTK垂直精度 1.5 cm+1 ppm 镜头焦距(mm) 35 RTK水平精度 1.0 cm+1 ppm 工作温度(℃) -20~50 
图 2 多角度贴近摄影测量消除摄影盲区Fig. 2. Multi-angle proximity photogrammetry eliminates photographic blind spots通过上述步骤成功消除了传统无人机摄影中拉花、阴影等问题，实现复杂高陡岩质斜坡岩体结构面高精度影像数据获取，为获取研究区高密度点云，精确解译岩体结构面信息提供有力支持.在无人机综合摄影测量过程中，采用WGS84作为参考坐标系，作业面积0.587 km2，影像重叠度为航向85%、旁向70%.共飞行4个架次，耗时100 min，飞行总里程30.548 km，共获取影像2 435张.将无人机综合摄影测量获取的全部影像导入三维建模软件中并输入相关飞行参数，通过图像特征匹配、空中三角测量、三维模型生成等步骤，完成高陡岩质斜坡三维模的型构建（图 3）.在WGS84坐标系下生成三维网络，设置点云密度为1 mm，生成研究区高密度点云.
图 3 无人机摄影航线设计图与高精度三维模型a.无人机摄影测量航线设计图；b.研究区高精度三维模型Fig. 3. UAV photography route design diagram and high-precision 3D model2.2 结构面优势分组
在察达工点高陡岩质斜坡上选取200 m×100 m的矩形窗口作为研究区，斜坡整体产状为113°∠73°，斜坡岩体节理裂隙发育.基于高精度三维模型，通过EPS软件共解译结构面1 495条.解译对象包括线状、面状、阶梯式、混合式出露结构面（图 4），优先解译大尺寸、控制性结构面，忽略迹长不足10 cm的小尺寸结构面.解译时基于3点共面原理，在三维模型的目标结构面上提取3个及以上特征点坐标，进一步拟合平面，计算其法向量、获取结构面产状信息.对解译后的结构面采用QPSO-FCM全局优化算法（宋盛渊等，2015，2024）对其进行优势分组，划分优势组数为4组，其中第1组陡倾坡外、第2组陡倾且与坡面大角度相交、第3组陡倾且与坡面近直交、第4组缓倾且与坡面大角度相交，具体结果见图 5，表 2.
图 4 不同出露形式结构面解译示意图a.线状出露结构面；b.面状出露结构面；c.阶梯式出露结果；d.混合式出露结构面Fig. 4. Interpretation diagram of different exposure structures表 2 研究区结构面优势分组统计Table Supplementary Table Grouping statistics of structural plane dominance in the study area组号 平均倾向(°) 平均倾角(°) 结构面数 1 128 76 354 2 242 69 408 3 197 82 278 4 216 7 455 2.3 结构面轮廓线提取
针对传统的人工接触式测量法难以适应高陡斜坡地形环境的问题，本文基于高精度三维模型生成的高密度点云，精细解译结构面形貌信息，求算其粗糙度.在自然界中，岩体结构面的出露形式复杂多样，主要分为线状和面状结构面.
线状结构面出露形式为二维曲线，只能根据其轮廓线估算结构面粗糙度系数JRC.因此，本研究采用EPS软件中的三维测图模式直接在其表面测量线状结构面轮廓线的起伏特征点（图 6），并将提取的结构面轮廓特征点按照1 mm的间距进行插值，以保证结构面信息提取的精度.在实际测量中，所得到的结构面轮廓点云为三维坐标且并不在同一平面上.为方便JRC的计算，需将点云坐标由三维转化至二维.根据Zhao et al.（2023）提出的坐标转换方法，采用最小二乘法拟合轮廓线基准平面，将点云投影至该平面上并建立新的坐标系.在新坐系中所有特征点都位于同一个平面上，其z轴坐标可视为0.以实现点云数据从三维至二维的转化.
对于面状结构面，框选目标区域，识别并裁剪出面状结构面高密度点云，由于获取的点云文件是不规则的，需将其导入MATLAB中，使用Griddata函数对点云插值处理，求取点云规则化后（x，y）坐标对应的高程z，点云间距设置为Δx=Δy，实现点云的规则化处理（葛云峰等，2017）.将处理后的结构面点云按照合适的尺寸选择其中出露较好、完整的典型区域（图 7），并沿一定采样方向剪裁，提取结构面轮廓线.
在选取目标结构面时尽量选取新鲜的、风化程度较低、出露完整的结构面.在剪裁高密度点云时，选取该结构面上较为平整的、微风化的区域提取高密度点云，从而获得更加准确的结构面形貌信息，最大程度上减少风化作用对研究结果的影响.
根据JRC尺寸效应的研究，当结构面尺寸达到1 m时，结构面粗糙度的表征参数逐渐趋于稳定，当结构面尺寸达到1.0~1.2 m时，表征参数趋于稳定不再增加（葛云峰等，2016）.因此本文在研究区中分别运用上述方法提取出露完整的20条线状与30条面状结构面轮廓线，其尺寸均为1 m.
3. 粗糙度分形特征
3.1 粗糙度盒维数计算
盒维数法是Walsh提出用于描述不规则二维曲线分形维数的方法（Walsh and Watterson，1993），该方法在计算时需要构造一个边长为R的正方形盒子，使结构面轮廓线完整的包裹在其中.将此正方形盒子分割为边长为r的小正方形盒子，r依次取R/2、R/4、R/8、R/16、R/32等特定值，分别计算出不同r值的盒子中与结构面轮廓线相交的最小正方形盒子数（图 8），记为N（r）.带入公式求出结构面的分形维数D，如公式（1）：
D=−limr→0lnN(r)lnr. (1) 根据公式以及图像分析可得，分形维数D即lnN（r）和ln（r）的图像的斜率.一般情况下，二维曲线的起伏越大，其分形维数D就越大，即岩体结构面越粗糙，分形维数D就越大.
3.2 JRC计算
Barton于1982年提出粗糙度系数的直边法概念公式（Barton and Bandis，1982），并根据200多组结构面统计分析推导出确定粗糙度系数JRC值的经验公式.直边法是根据不同取样长度结构面表面轮廓线的起伏度确定粗糙度系数JRC，其原理见图 9.在Barton直边法的基础上，杜时贵等（1996）发现直边法中求得的结构面粗糙度系数JRC有可能大于20，且起伏度为毫米级的岩体结构面也同样适用直边法.为此，杜时贵改进了Barton的直边法，推导出修正直边法（杜时贵等，1996）的数学表达式如公式（2）：
JRC=49e29L0/450Lnarctan(8Ry/Ln), (2) 式中：L0表示结构面采样长度（cm）；Ln表示结构面轮廓线长度（cm）；Ry表示结构面的最大起伏度（cm）.
在使用杜时贵的修正直边法实测由点云构成的结构面轮廓线粗糙度时，由于结构面轮廓线并非水平放置，不能直接通过软件测量采样长度Ln和结构面的最大起伏度Ry，需进行一系列坐标变换，其具体步骤如下：（1）将剪裁后的三维结构面轮廓线点云经上文的坐标转换后输出为二维坐标形式；（2）将二维结构面轮廓线首尾相连并将连线设为结构面轮廓基准线，经过旋转、平移将基准线设定为x=0，并将其余点跟随基准线共同进行坐标变换；（3）在变换后的点云坐标中测量采样长度Ln和结构面的最大起伏度Ry.
由于Ln和Ry的测量操作复杂，而利用盒维数法仅需通过简单的编程，导入点云文件，即可计算出分形维数D，因此需要建立一个JRC与盒维数法求解的分形维数D的表达式.
3.3 新公式拟合与验证
本文分别利用盒维数法和修正直边法求解上文中提取的50条结构面轮廓线的分形维数D和粗糙度系数JRC，结果如表 3所示.
表 3 分形维数D与粗糙度系数JRC计算结果Table Supplementary Table Calculation results of fractal dimension D and roughness coefficient JRC编号 分形维数D JRC 编号 分形维数D JRC 1 1.010 6 5.339 2 26 1.017 9 11.836 7 2 1.011 5 5.715 3 27 1.018 2 12.024 6 3 1.012 2 5.966 7 28 1.0184 12.327 3 4 1.013 1 6.322 9 29 1.018 4 12.119 3 5 1.013 3 7.129 5 30 1.018 6 13.376 4 6 1.013 9 7.658 0 31 1.018 8 12.958 8 7 1.014 1 7.925 4 32 1.018 8 13.112 7 8 1.014 2 8.394 2 33 1.019 1 13.813 9 9 1.014 2 7.862 1 34 1.019 1 14.199 2 10 1.014 3 8.567 6 35 1.019 3 14.219 6 11 1.014 5 8.482 7 36 1.018 9 14.403 8 12 1.014 6 8.601 4 37 1.019 8 15.002 2 13 1.015 3 9.183 5 38 1.019 9 14.892 8 14 1.015 5 9.956 5 39 1.020 5 15.383 5 15 1.015 5 10.039 8 40 1.020 7 15.402 3 16 1.015 7 9.822 4 41 1.020 9 15.118 9 17 1.016 2 10.054 7 42 1.021 2 15.944 7 18 1.015 5 10.139 6 43 1.021 5 16.325 7 19 1.015 9 10.336 2 44 1.022 8 17.286 4 20 1.015 8 10.932 0 45 1.023 3 18.863 0 21 1.016 5 11.122 4 46 1.023 6 19.005 3 22 1.016 7 11.197 1 47 1.023 6 19.553 7 23 1.017 0 12.335 4 48 1.023 7 18.726 5 24 1.017 3 11.446 5 49 1.024 1 19.865 7 25 1.017 6 11.638 1 50 1.024 3 19.066 2 在计算分形维数D与粗糙度系数JRC值后，为建立二者之间的定量关系，拟采用如公式（3）形式的表达式：
JRC=a+b×D. (3) 拟合出新公式如公式（4），拟合优度R2=0.985 39，可见拟合结果较为理想（图 10）.
JRC=−1118.7011+1111.2214D. (4) 为验证新公式的准确性，将Barton标准剖面线原始图片进行修复，去除剖面线附近的杂点并使用增加像素点的方法补齐剖面线存在的断口.将处理后的图片导入Matlab中进行精细数字化处理，得到由离散点构成的标准剖面线.求解数字化Barton标准剖面线的分形维数D并带入新公式中计算JRC（图 11）.可见，利用新公式计算Barton标准剖面线的粗糙度系数JRC，所得结果均在合理范围内，证明新公式较为准确，具有可行性.
4. 结构面力学成因与各向异性
4.1 结构面力学成因
自然界中结构面的地质力学成因主要有压剪和拉张两种（伍法权和伍劼，2022）.本文的研究对象为研究区内压剪作用形成的压剪性结构面和拉张作用形成的拉张性结构面.
压剪性结构面是在地质应力作用下岩石发生压剪变形所形成的.当研究区内高陡岩体在巨大的挤压或挤扭作用下，岩体内部会产生剪切应力，使岩体的软弱面和断层面发生相对滑动，进而形成压剪性结构面（图 12a）；拉张性结构面是由于研究区岩体快速隆升，在地质应力作用下发生拉张变形所形成的.当岩体受到拉张作用时，岩体内部的裂隙、断层面会被拉开，形成拉张性结构面.同时，裂隙和断层面两侧的岩体也会受到张力而发生拉伸、脆性断裂等变形，形成拉张性结构面（图 12b）.
图 12 压剪性与拉张性结构面成因示意图a.压剪性结构面；b.拉张性结构面Fig. 12. Genetic diagram of compressive shear and tensile structural plane4.2 JRC各向异性
为研究粗糙度系数JRC的各向异性，根据研究区结构面优势分组结果，在缓倾且与坡面大角度相交的第4组中选取15个面状出露的压剪性结构面J1~J15；在陡倾坡外的第1组、陡倾且与坡面大角度相交的第2组、陡倾且与坡面近直交的第3组中各选取5个面状出露的拉张性结构面L1~L5、L6~L10、L11~15，提取高密度点云数据并导入点云处理软件中剪裁.为保证每条结构面轮廓线的取样长度均为1 m，再将1 m×1 m的正方形结构面剪裁成直径为1 m的圆形结构面.将剪裁后的圆形结构面在笛卡尔坐标系中从0°以10°的角度增量选取36个采样方向（图 13），剪裁结构面轮廓点云，并通过坐标转换得到结构面轮廓线（图 14）.将全部结构面轮廓线求算分形维数D并代入新公式中求出粗糙度系数JRC，结果如图 15所示.通过分析研究区压剪性与拉张性结构面JRC各向异性计算结果，计算每个结构面JRC各向平均值JRC、JRC方差σJRC、JRC各向极差RJRC并绘制成直方图（图 16a~16c），其中JRC极差RJRC计算公式如公式（5）：
RJRC=JRCmax−JRCmin, (5) 
图 15 压剪性与拉张性结构面粗糙度各向异性图a.压剪性结构面J1~J15；b.拉张性结构面L1~L15Fig. 15. Anisotropy of compressive shear and tensile structural planes
图 16 结构面JRC各向异性统计直方图a.JRC平均值；b.JRC方差；c.JRC极差；d.JRC平均变化量Fig. 16. Statistical histograms of structural plane JRC anisotropy式中：JRCmax与JRCmin分别表示该结构面各个方向上的JRC最大值与最小值.
计算每个结构面相邻两个采样方向上JRC变化量绝对值并求算平均变化量ΔJRC，绘制成直方图（图 16d），平均变化量ΔJRC计算公式如公式（6）：
ΔJRC=n∑i=1JRCi+1−JRCin−1, (6) 式中：JRCi+1与JRCi分别表示该结构面相邻两个采样方向上的JRC值；n为该结构面总采样条数.
由图 15可见，结构面粗糙度存在明显的各向异性.压剪性结构面存在JRC较大的优势方向与JRC较小的劣势方向，且优势方向与劣势方向夹角约为90°.根据研究区结构面优势分组结果，压剪性结构面所在的第4组倾向为216°与JRC劣势方向20°~200°较为接近，说明沿剪切滑动方向JRC较小，垂直于剪切滑动方向JRC较大，呈现椭圆状分布.拉张性结构面JRC在不同方向上也具有各向异性，但无明显的变化规律，不存在显著的优势方向与劣势方向，相邻采样方向上变化较大，呈现刺状分布.
由图 16可见，研究区压剪性结构面JRC平均值小于拉张性结构面，可见结构面在发生剪切滑动时，摩擦作用使结构面变得相对平直.研究区拉张性结构面不同采样方向上的JRC极差相较于压剪性更大，且拉张性结构面的方差均值为1.42，压剪性结构面方差均值为1.04，可见拉张性结构面各采样方向上JRC离散程度较大，压剪性结构面各采样方向上JRC离散程度较小.结合图 15和图 16分析，随着采样角度的变化，压剪性结构面JRC变化较小，且整体上由剪切滑动方向至垂直剪切滑动方向递增，拉张性结构面JRC变化较大且整体上呈现为无规律波动.
综上所述，压剪性与拉张性结构面由于其地质力学成因不同，JRC各向异性也存在差异，压剪性结构面的JRC各向异性具有一定规律，而拉张性结构面的JRC各向异性无明显规律，可证明结构面的粗糙度各向异性与其地质力学成因有关.为保障工程安全，在实际工程中JRC取值建议选取该结构面的最小值.因此，在野外实测结构面粗糙度时，应先根据地质力学成因对结构面进行分类.若为压剪性结构面，应依据其产状判别剪切滑动方向，沿滑动方向采样并计算JRC；若为拉张性结构面，应沿不同方向大量采样计算JRC.
5. 结论
本文以藏东南某铁路察达工点高陡岩质斜坡为研究对象，利用无人机综合摄影测量技术获取结构面高密度点云，进一步裁剪出50条结构面轮廓线，分别采用盒维数法和修正直边法求算其分形维数D和结构面粗糙度系数JRC，拟合新公式并验证，用以对压剪性和拉张性结构面进行JRC各向异性研究，得出以下结论：
（1）将修复后的数字化Barton标准剖面线代入拟合的新公式检验其准确度，所得结果均在合理范围内，表明新公式计算简便且结果可靠.
（2）结构面粗糙度的各向异性与其地质力学成因有关.压剪性结构面的粗糙度各向异性具有显著规律，且JRC沿剪切滑动方向较小为劣势方向，垂直于剪切滑动方向较大为优势方向，优势方向与劣势方向夹角约为90°.JRC由剪切滑动方向至垂直剪切滑动方向递增，呈椭圆状分布；拉张性结构面的粗糙度也具有各向异性，但无明显变化规律，不存在明显的优势方向与劣势方向，JRC随采样角度变化呈现无规律波动，整体上呈刺状分布.
（3）在野外实测结构面粗糙度时，应先对结构面进行地质力学成因分类，对于压剪性结构面可根据其产状判断剪切滑动方向，沿剪切滑动方向采样并计算JRC.对于拉张性结构面可沿不同方向大量采样并计算JRC.为保障安全，实际工程中建议选取该结构面JRC的最小值.
值得指出，本文所拟合的JRC和D的新公式与粗糙度各向异性的结论还需更多的工程实例与数据进行验证.在未来的研究中，将此计算方法应用于不同研究区域，深入分析不同力学成因结构面粗糙度的各向异性，以进一步验证公式、巩固结论，使其能够更好的推广应用，服务于工程实例.
-
图 1 研究区地理位置、地形地貌与斜坡全貌
a.研究区地理位置图；b.研究区地形地貌图；c.研究区高陡斜坡全貌图
Fig. 1. Geographical location, landform and slope overall appearance of the study area
图 2 多角度贴近摄影测量消除摄影盲区
Fig. 2. Multi-angle proximity photogrammetry eliminates photographic blind spots
图 3 无人机摄影航线设计图与高精度三维模型
a.无人机摄影测量航线设计图；b.研究区高精度三维模型
Fig. 3. UAV photography route design diagram and high-precision 3D model
图 4 不同出露形式结构面解译示意图
a.线状出露结构面；b.面状出露结构面；c.阶梯式出露结果；d.混合式出露结构面
Fig. 4. Interpretation diagram of different exposure structures
图 12 压剪性与拉张性结构面成因示意图
a.压剪性结构面；b.拉张性结构面
Fig. 12. Genetic diagram of compressive shear and tensile structural plane
图 15 压剪性与拉张性结构面粗糙度各向异性图
a.压剪性结构面J1~J15；b.拉张性结构面L1~L15
Fig. 15. Anisotropy of compressive shear and tensile structural planes
图 16 结构面JRC各向异性统计直方图
a.JRC平均值；b.JRC方差；c.JRC极差；d.JRC平均变化量
Fig. 16. Statistical histograms of structural plane JRC anisotropy
表 1 无人机和相机相关参数
Table 1. UAV and camera related parameters
大疆M300-RTK无人机 禅思P1相机 无人机总重(kg) 6.3 重量(g) 800 最大飞行时长(min) 55 相机尺寸(mm) 198×166×129 最大可承受风速(m/s) 15 传感器面积(mm2) 35.9×24 RTK垂直精度 1.5 cm+1 ppm 镜头焦距(mm) 35 RTK水平精度 1.0 cm+1 ppm 工作温度(℃) -20~50 
下载: 导出CSV
表 2 研究区结构面优势分组统计
Table 2. Grouping statistics of structural plane dominance in the study area
组号 平均倾向(°) 平均倾角(°) 结构面数 1 128 76 354 2 242 69 408 3 197 82 278 4 216 7 455
下载: 导出CSV
表 3 分形维数D与粗糙度系数JRC计算结果
Table 3. Calculation results of fractal dimension D and roughness coefficient JRC
编号 分形维数D JRC 编号 分形维数D JRC 1 1.010 6 5.339 2 26 1.017 9 11.836 7 2 1.011 5 5.715 3 27 1.018 2 12.024 6 3 1.012 2 5.966 7 28 1.0184 12.327 3 4 1.013 1 6.322 9 29 1.018 4 12.119 3 5 1.013 3 7.129 5 30 1.018 6 13.376 4 6 1.013 9 7.658 0 31 1.018 8 12.958 8 7 1.014 1 7.925 4 32 1.018 8 13.112 7 8 1.014 2 8.394 2 33 1.019 1 13.813 9 9 1.014 2 7.862 1 34 1.019 1 14.199 2 10 1.014 3 8.567 6 35 1.019 3 14.219 6 11 1.014 5 8.482 7 36 1.018 9 14.403 8 12 1.014 6 8.601 4 37 1.019 8 15.002 2 13 1.015 3 9.183 5 38 1.019 9 14.892 8 14 1.015 5 9.956 5 39 1.020 5 15.383 5 15 1.015 5 10.039 8 40 1.020 7 15.402 3 16 1.015 7 9.822 4 41 1.020 9 15.118 9 17 1.016 2 10.054 7 42 1.021 2 15.944 7 18 1.015 5 10.139 6 43 1.021 5 16.325 7 19 1.015 9 10.336 2 44 1.022 8 17.286 4 20 1.015 8 10.932 0 45 1.023 3 18.863 0 21 1.016 5 11.122 4 46 1.023 6 19.005 3 22 1.016 7 11.197 1 47 1.023 6 19.553 7 23 1.017 0 12.335 4 48 1.023 7 18.726 5 24 1.017 3 11.446 5 49 1.024 1 19.865 7 25 1.017 6 11.638 1 50 1.024 3 19.066 2
下载: 导出CSV
-
Barton, N., 1973. Review of a New Shear-Strength Criterion for Rock Joints. Engineering Geology, 7(4): 287-332. https://doi.org/10.1016/0013-7952(73)90013-6 Barton, N., Bandis, S., 1982. Effects of Block Size on the Shear Behavior of Jointed Rock. The 23rd US Symposium on Rock Mechanics, Berkeley, 739-760. Barton, N., Choubey, V., 1977. The Shear Strength of Rock Joints in Theory and Practice. Rock Mechanics and Rock Engineering, 10(1): 1-54. https://doi.org/10.1007/BF01261801 Belem, T., Homand-Etienne, F., Souley, M., 2000. Quantitative Parameters for Rock Joint Surface Roughness. Rock Mechanics and Rock Engineering, 33(4): 217-242. https://doi.org/10.1007/s006030070001 Chai, B., Shi, X. S., Du, J., et al., 2022. How to Realize Elaborated Analysis of Regional Rock Mass Structure? A Review and Idea. Earth Science, 47(12): 4629-4646(in Chinese with English abstract). Chen, S. J., Zhu, W. C., Zhang, M. S., et al., 2012. Fractal Description of Rock Joints Based on Digital Image Processing Technique. Chinese Journal of Geotechnical Engineering, 34(11): 2087-2092(in Chinese with English abstract). Cui, W., Gao, D. Y., Wang, X. H., et al., 2023. Identification of Rocky Ledge on Steep and High Slopes Based on Aerial Photogrammetry. Earth Science, 48(9): 3378-3388 (in Chinese with English abstract). Du, S. G., 1992. Simple Profile Instrument and Its Application on Studying Joint Roughness Coefficient of Rock. Geological Science and Technology Information, 11(3): 91-95(in Chinese with English abstract). Du, S. G., 1997. The Practicability of Fractal Methods on Estimating Rock Joint Roughness Coefficient. Earth Science, 22(6): 665-668(in Chinese with English abstract). Du, S. G., Chen, Y., Fan, L. B., 1996. Maethematical Expression of JRC Modified Straight Edge. Journal of Engineering Geology, 4(2): 36-43(in Chinese with English abstract). Du, S. G., Tang, H. M., 1993. A Study on the Anisotropy of Joint Roughness Coefficient in Rock Mass. Journal of Engineering Geology, 1(2): 32-42(in Chinese with English abstract). El-Soudani, S. M., 1978. Profilometric Analysis of Fractures. Metallography, 11(3): 247-336. https://doi.org/10.1016/0026-0800(78)90045-9 Ge, Y. F., Chen, Y., Wang, L. Q., et al., 2017. Estimation of 3D Joint Matching Coefficient of Rock Discontinuities Using Point Cloud Data Alignment Technology. Rock and Soil Mechanics, 38(11): 3385-3393(in Chinese with English abstract). Ge, Y. F., Tang, H. M., Wang, L. Q., et al., 2016. Anisotropy, Scale and Interval Effects of Natural Rock Discontinuity Surface Roughness. Chinese Journal of Geotechnical Engineering, 38(1): 170-179(in Chinese with English abstract). Gu, D. Z., 1979. Foundation of Rock Mass Engineering Geomechanics. Science Press, Beijing (in Chinese). Guo, C. B., Wu, R. A., Zhong, N., et al., 2024. Large Landslides along Active Tectonic Zones of Eastern Tibetan Plateau: Background and Mechanism of Landslide Formation. Earth Science, 49(12): 4635-4658(in Chinese with English abstract). Hong, C. J., Huang, M., Xia, C. C., et al., 2020. Study of Size Effect on the Anisotropic Variation Coefficient of Rock Joints. Rock and Soil Mechanics, 41(6): 2098-2109(in Chinese with English abstract). Kulatilake, P. H. S. W., Um, J., 1999. Requirements for Accurate Quantification of Self-Affine Roughness Using the Roughness-Length Method. International Journal of Rock Mechanics and Mining Sciences, 36(1): 5-18. https://doi.org/10.1016/S0148-9062(98)00170-3 Kulatilake, P. H. S. W., Um, J., Pan, G., 1997. Requirements for Accurate Estimation of Fractal Parameters for Self-Affine Roughness Profiles Using the Line Scaling Method. Rock Mechanics and Rock Engineering, 30(4): 181-206. https://doi.org/10.1007/BF01045716 Luo, Z. J., Zhang, Q. Z., Shi, Z. M., et al., 2022. Shear Creep Characteristics of Red Sandstone Discontinuities Considering Different Morphologies. Earth Science, 47(12): 4484-4497(in Chinese with English abstract). Mandelnorot, B. B., 1982. The Fractal Geometry of Nature. Freeman, New York. Pan, X. J., Zhang, W., Sun, Q., et al., 2024. Complex Structural Plane Distribution of High⁃Steep Rock Slope and Division of Statistical Homogeneous Zones. Earth Science, 49(9): 3334-3346(in Chinese with English abstract). Song, S. Y., Huang, D., Sui, J. X., et al., 2024. Analysis of Size Effect for Shear Characteristics of Rock Mass Based on 3D Fracture Network. Journal of Harbin Institute of Technology, 56(3): 9-18(in Chinese with English abstract). Song, S. Y., Wang, Q., Chen, J. P., et al., 2015. A Method for Multivariate Parameter Dominant Partitioning of Discontinuities of Rock Masses. Rock and Soil Mechanics, 36(7): 2041-2048(in Chinese with English abstract). Song, S. Y., Zhao, M. Y., Zhang, W., et al., 2024. Research on Fine Collection and Interpretation Methods of Discontinuities on High-Steep Rock Slopes Based on UAV Multi-Angle Nap-of-the-Object Photogrammetry. Bulletin of Engineering Geology and the Environment, 83(4): 142. https://doi.org/10.1007/s10064-024-03646-5 Tse, R., Cruden, D. M., 1979. Estimating Joint Roughness Coefficients. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 16(5): 303-307. https://doi.org/10.1016/0148-9062(79)90241-9 Walsh, J. J., Watterson, J., 1993. Fractal Analysis of Fracture Patterns Using the Standard Box-Counting Technique: Valid and Invalid Methodologies. Journal of Structural Geology, 15(12): 1509-1512. https://doi.org/10.1016/0191-8141(93)90010-8 Wang, F. Y., Chen, J. P., Yang, G. D., et al., 2012. Solution Models of Geometrical Information of Rock Mass Discontinuities Based on Digital Close Range Photogrammetry. Journal of Jilin University (Earth Science Edition), 42(6): 1839-1846(in Chinese with English abstract). Wang, S. H., Hou, Q. K., Yong, R., et al., 2023. The Effectiveness of Sample Selection Methods in Study of Shear Strength Anisotropy of Rock Joints. Engineering Mechanics, 40(1): 168-179(in Chinese with English abstract). Wang, S. J., 1990. Engineering Geomechanics Analysis of Dam Foundation Rock Mass. Science Press, Beijing (in Chinese). Wu, F. Q., Wu, J., 2022. Theory and Application of Statistical Rock Mechanics. Science Press, Beijing (in Chinese). Xie, H. P., 1995. Fractal Description of Rock Joints. Chinese Journal of Geotechnical Engineering, 17(1): 18-23(in Chinese with English abstract). Xie, H. P., Wang, J. N., 1999. Direct Fractal Measurement of Fracture Surfaces. International Journal of Solids and Structures, 36(20): 3073-3084. https://doi.org/10.1016/S0020-7683(98)00141-3 Yang, Z., Li, B. T., Song, S. Y., et al., 2023. Quantitative Analysis of Structural Surface Roughness Based on UAV Point Cloud and Improved RS Characterization Method. Journal of Harbin Institute of Technology, 55(11): 72-81(in Chinese with English abstract). You, Z. C., Wang, L. Q., Yang, Y. X., et al., 2014. Anisotropic Research on Shear Strength Parameters of Discontinuity Based on Three-Dimensional Laser Scanning Technology. Chinese Journal of Rock Mechanics and Engineering, 33(S1): 3003-3008(in Chinese with English abstract). Zhang, Z. F., Huang, M., Tang, Z. C., 2024. Peak Shear Strength Criterion for Discontinuities with Different Rock Types Based on Revisiting Frictional Angle. Earth Science, 49(8): 2826-2838(in Chinese with English abstract). Zhao, M. Y., Chen, J. P., Song, S. Y., et al., 2023. Proposition of UAV Multi-Angle Nap-of-the-Object Image Acquisition Framework Based on a Quality Evaluation System for a 3D Real Scene Model of a High-Steep Rock Slope. International Journal of Applied Earth Observation and Geoinformation, 125: 103558. https://doi.org/10.1016/j.jag.2023.103558 柴波, 史绪山, 杜娟, 等, 2022. 如何实现区域岩体结构精细化分析?综述与设想. 地球科学, 47(12): 4629-4646. 陈世江, 朱万成, 张敏思, 等, 2012. 基于数字图像处理技术的岩石节理分形描述. 岩土工程学报, 34(11): 2087-2092. 崔溦, 高德宇, 王轩毫, 等, 2023. 基于航空影像的高陡边坡突岩识别方法. 地球科学, 48(9): 3378-3388. 杜时贵, 1992. 简易纵剖面仪及其在岩体结构面粗糙度系数研究中的应用. 地质科技情报, 11(3): 91-95. 杜时贵, 1997. JRC分形估测方法的实用性. 地球科学, 22(6): 665-668. 杜时贵, 陈禹, 樊良本, 1996. JRC修正直边法的数学表达. 工程地质学报, 4(2): 36-43. 杜时贵, 唐辉明, 1993. 岩体断裂粗糙度系数的各向异性研究. 工程地质学报, 1(2): 32-42. 葛云峰, 陈勇, 王亮清, 等, 2017. 基于点云数据对齐技术的岩体结构面三维吻合度求取. 岩土力学, 38(11): 3385-3393. 葛云峰, 唐辉明, 王亮清, 等, 2016. 天然岩体结构面粗糙度各向异性、尺寸效应、间距效应研究. 岩土工程学报, 38(1): 170-179. 谷德振, 1979. 岩体工程地质力学基础. 北京: 科学出版社. 郭长宝, 吴瑞安, 钟宁, 等, 2024. 青藏高原东部活动构造带大型滑坡成灾背景与灾变机制. 地球科学, 49(12): 4635-4658. 洪陈杰, 黄曼, 夏才初, 等, 2020. 岩体结构面各向异性变异系数的尺寸效应研究. 岩土力学, 41(6): 2098-2109. 罗泽军, 张清照, 石振明, 等, 2022. 考虑表面形貌特征的岩体结构面蠕变特性. 地球科学, 47(12): 4484-4497. 潘晓娟, 张文, 孙琦, 等, 2024. 高陡岩质斜坡复杂结构面分布规律及统计均质区划分. 地球科学, 49(9): 3334-3346. 宋盛渊, 黄迪, 隋佳轩, 等, 2024. 基于三维裂隙网络的岩体剪切特性尺寸效应分析. 哈尔滨工业大学学报, 56(3): 9-18. 宋盛渊, 王清, 陈剑平, 等, 2015. 岩体结构面的多参数优势分组方法研究. 岩土力学, 36(7): 2041-2048. 王凤艳, 陈剑平, 杨国东, 等, 2012. 基于数字近景摄影测量的岩体结构面几何信息解算模型. 吉林大学学报(地球科学版), 42(6): 1839-1846. 王述红, 侯钦宽, 雍睿, 等, 2023. 结构面抗剪强度各向异性研究中试样获取方法有效性分析. 工程力学, 40(1): 168-179. 王思敬, 1990. 坝基岩体工程地质力学分析. 北京: 科学出版社. 伍法权, 伍劼, 2022. 统计岩体力学理论与应用. 北京: 科学出版社. 谢和平, 1995. 岩石节理的分形描述. 岩土工程学报, 17(1): 18-23. 杨泽, 李保天, 宋盛渊, 等, 2023. 基于无人机点云与改进RS表征法的结构面粗糙度定量分析. 哈尔滨工业大学学报, 55(11): 72-81. 游志诚, 王亮清, 杨艳霞, 等, 2014. 基于三维激光扫描技术的结构面抗剪强度? 参数各向异性研究. 岩石力学与工程学报, 33(增刊1): 3003-3008. 张志飞, 黄曼, 唐志成, 2024. 基于摩擦角修正的异性岩石结构面抗剪强度准则. 地球科学, 49(8): 2826-2838. 期刊类型引用(1)
1. 徐锟山. 地质找矿技术在资源勘查中的应用. 世界有色金属. 2025(05): 190-192 . 
百度学术其他类型引用(0)
-
下载:
百度学术
点击查看大图
计量
- 文章访问数: 100
- HTML全文浏览量: 90
- PDF下载量: 18
- 被引次数: 1
