基于BAP的岩体结构面粗糙度与峰值抗剪强度研究

发布时间：2018-01-17 16:18

本文关键词：基于BAP的岩体结构面粗糙度与峰值抗剪强度研究　出处：《中国地质大学》2014年博士论文　论文类型：学位论文

【摘要】：岩体是自然界中具有一定结构特征的地质体,其内部发育有各种各样具有一定方向、规模和形态的结构面,如层面、节理、断层、裂隙、不整合接触面等等。结构面的存在使得岩体不同于其他材料,物理力学性质存在不连续性、各向异性、非均一性等特点,同时也降低了岩体的完整性和整体强度。岩体往往不是因为本身材料破坏发生失稳,而是沿着岩体内部的软弱结构面发生破坏。甚至可以说,岩体的稳定性在一定程度上是由这些结构面(特别是控制性结构面)所决定的。因此,开展结构面抗剪强度力学参数的研究是评价岩体稳定性的关键所在。而作为影响岩体结构面抗剪强度的重要因素之一,粗糙度与抗剪强度之间具有明显的相关性,结构面越粗糙,抗剪强度越大。基于粗糙度的结构面峰值抗剪强度估算模型研究一直是岩体力学与工程地质领域的重要研究课题,估算模型的成败在于能否寻找出合理的指标来表征岩体结构面的真实粗糙程度。然而,目前关于结构面粗糙度的研究仍然存在以下几个问题： (1)采集装备落后,试验周期长,不能满足岩体稳定性的快速评价要求。对于发展迅速的岩体破坏过程,如高速远程滑坡,还有一些重大地质灾害抢险工程,时间就是生命,需要研究学者能够在短时间内给出评价结果,并预测受灾范围,提出合理的工程治理方案。而恰恰由于测量仪器本身问题,却体现不出估算模型的优势。 (2)评价方法对采样方向、采样尺寸以及采样间距具有很大的依赖性,缺乏系统地开展粗糙度影响因素的研究及其相应合理的解决对策。 (3)粗糙度评价指标仅仅依据结构面的几何信息,忽视了结构面粗糙度的力学特征,结构面粗糙度与结构面的受力状态、受力方向是密切相关的。本论文研究目的就是借助先进的采样手段,系统分析采样参数(采样方向、采样尺寸以及采样间距)对结构面粗糙度评价结果的影响,进而提出一种快速精确的粗糙度评价方法,并力求在此基础上建立结构面峰值抗剪强度估算模型。以重庆武隆鸡尾山高速远程滑坡研究区域的灰岩结构面为例,本论文主要开展以下几方面的研究工作： (1)天然岩体结构面粗糙度数据的快速采集。LS05-1岩体结构面为散落在堆积区内的灰岩层面,非原位产状,采用ILRIS-36D三维激光扫描仪,获取野外大尺寸结构面粗糙度点云数据,这部分数据主要用于岩体结构面粗糙度影响因素、基于光亮面积百分比(BAP)的粗糙度表征方法和结构面峰值抗剪强度估算模型研究工作中,原始点云数据扫描尺寸(ROI)为4.7×2.5m,采样间距为0.8mm,共采集扫描点约1800多万个。JWO2-1BL、 JW02-2BL、JW03-1BL、JW03-2BL和JW03-3BR岩体结构面为采自于滑动面上的灰岩节理,运用VZ-400三维激光扫描仪,获取室内小尺寸结构面粗糙度点云数据,这部分数据主要服务于结构面峰值抗剪强度估算模型的研究工作,为数值模拟试验提供参数并验证计算结果,原始点云数据平均扫描尺寸(ROI)为0.181×0.362m,采样间距为1.0mm,每个结构面共采集65500左右个扫描点。对采集的原始点云数据进行预处理,主要包括点云数据解析、去噪和精简、几何模型建立三大部分,综合协调计算机运算能力、精度与处理时间等因素,最后LS05-1结构面精简成尺寸为4.70×1.14m,间距为0.01m,约54000个扫描点的点云数据；而JW02-1BL、JW02-2BL、JWO3-1BL、JW03-2BL和JW03-3BR结构面精简成平均尺寸为0.06×0.06m,间距保持不变仍为1.0mm,每块约3600个扫描点的点云数据,利用编写的Matlab程序,建立岩体结构面三维几何模型。 (2)岩体结构面粗糙度影响因素研究。以LS05-1岩体结构面预处理后的点云数据为主,Weierstrass-Mandelbrot分形随机函数产生为辅,分别从二维与三维角度,分析采样方向、采样尺寸以及采样间距对结构面粗糙度评价结果影响。分析过程中做到单一变量原则,其余两因素保持固定。以15。为间隔,选取从0。到345。共24个采样方向,采用分形理论计算各个方向下的结构面二维剖面分形维数D,对于三维结构面,采用改进的赤平投影极点图法求出相应方向下的法向量数量NV。选取9种不同尺寸(长度或面积)的岩体结构面,基于分形理论求出相应尺度下的分形维数D。以0.01m为间隔,选取从0.01m到1.00m共100种采样间距,采用分形理论计算相应间距下的分形维数D。并且针对每一种影响因素,提出切实可行的解决方法。 (3)基于光亮面积百分比(BAP)粗糙度表征方法的提出。从无充填低法向应力下的岩体结构面剪切破坏机制入手,以结构面中面向剪切方向的微小平面为研究对象,提出表征结构面粗糙度的新指标-BAP。基于Matlab软件平台,编写了光源模拟程序,并基于图像分割技术,给出BAP的定义和计算方法。对光源模拟与图像分割过程中涉及的两重要参数-光源入射角β和灰度阈值TL取值问题进行了探讨。通过描述结构面粗糙度的各向异性、尺寸效应以及间距效应,检验BAP法的适用性。 (4)基于BAP的结构面峰值抗剪强度估算模型的建立。基于Matlab软件平台编写XYZ23DEC程序,简便快捷地建立起结构面三维复杂离散元模型。通过与室内直剪试验结果对比,来调整和验证数值试验所需的岩块、节理力学参数。在此参数基础上,以抗剪强度的各向异性、尺寸效应以及对法向应力敏感性等特征为参照点,验证数值模型的可靠性。设置合理的采样尺寸与采样间距常量,考虑结构面的粗糙度BAP、各向异性(剪切方向α)、岩块强度(JCS)、结构面基本摩擦角(φb)以及法向应力(σ)等参数,开展60例结构面直剪数值试验,统计分析模拟结果,建立粗糙度评价指标BAP与结构面峰值抗剪强度之间的关系,提出后者估算模型。 (5)结构面峰值抗剪强度估算模型的应用。将该模型应用于重庆武隆鸡尾山高速远程滑坡岩体稳定性评价当中,结合研究区域具体的粗糙度BAP、法向应力σ、基本摩擦角φh以及壁面强度JCS参数,估算出滑动面抗剪力学参数(内摩擦角φ和黏聚力C)的变化区间。选择最小值、最大值、平均值三种组合,基于大型三维离散元软件3DEC,再现滑体初期破坏演化过程,分析其变形特征。在上述研究内容与方法的基础上,得出如下几点结论： (1)运用三维激光扫描技术可以实现岩体结构面粗糙度数据的快速采集目的。扫描时间与扫描间距、扫描范围直接相关。扫描间距越小,扫描范围越大,所花费的扫描时间会随之增加。除此之外,扫描时间还受到非技术因素的影响,如操作不当、电量供给、操作环境等等。扫描得到的原始点云数据需要对其进行预处理工作,点云数据的预处理是一个重复的过程,需要不断根据反馈进行修正与改善,预处理工作花费的时间很大程度上取决于工作人员的专业素质,操作越熟练所用时间越短。 (2)采样方向对岩体结构面粗糙度评价影响分析显示：不同采样方向上的粗糙度之间存在较大差异,解决采样方向引入误差的最直接最有效方法,就是在数据采集过程中,确保采样方向与研究方向保持一致,研究方向包括滑动方向、剪切方向、渗流方向等等,这些研究方向需要在地质调查阶段进行判别,如滑动方向可根据滑动面上的擦痕和阶步等现象进行确定。采样尺寸对岩体结构面粗糙度评价影响分析显示：不同采样尺寸的结构面其粗糙度也各不相同,粗糙度会随着采样尺寸的增加而减小,当尺寸增加到有效尺寸(3.5m和3.5×0.84m)时,结构面粗糙度基本趋于平稳不再下降,通过增加采样尺寸的方法可以有效消除粗糙度的尺寸效应,但是受目前测量仪器的量程限制,也可以通过增加测量数量来弥补测量尺寸的不足,同样可以达到有效消除尺寸效应的目的,该思路在当前阶段具有较强的实用价值。采样间距对岩体结构面粗糙度评价影响分析显示：不同采样间距的结构面具有不同的粗糙度,整体上随着采样间距的减小,更多细节信息被获取,相应的粗糙度也越来越大,当采样间距小到有效间距时,粗糙度上升趋势变缓,趋于稳定状态,同时当采样间距过大时,大量信息被忽略,使得岩体结构面不再具备分形特征。有效采样问距与采样尺寸呈正比关系,建议采样间距与采样尺寸的比值不应高于0.05,在此条件下,可以有效降低采样间距对粗糙度评价结果的影响。 (3)在无充填低法向应力作用下的岩体结构面剪切破坏过程中,面向剪切方向的微小平面主要发生剪切破坏,抗剪作用较大,结构面粗糙度表征方法的研究对象应集中在该接触部分。灰度阈值TL的选取至关重要,应结合结构面剪切破坏力学机制,遵循具体情况具体分析原则进行选取。模拟光源入射角β对BAP计算结果存在显著影响,BAP会随着入射角β的增大而增大,建议最佳入射角β取值范围为35。与70。之间。具体案例应用中,入射角β与灰度阈值TL的选取应该相互协调配合,给出合理组合。提出的BAP表征新方法可以很好地反映出岩体结构面粗糙度的各向异性、尺寸效应以及间距效应特性,该法是综合考虑结构面几何与力学因素下提出的粗糙度表征新方法,具有独特优势。 (4)结合室内岩体结构面直剪试验展开的数值模拟结果显示：数值试验可以很好地体现出结构面抗剪强度的各向异性、尺寸效应以及与法向应力密切相关等特征。用于图像分割中的灰度阈值TL,其数值选取受法向应力σ的影响。以重庆武隆鸡尾山高速远程滑坡堆积区的LS05-1结构面为例,在光源入射角β=35°的前提下,法向应力所对应的合理灰度阈值TL如下所示, 法向应力σ=0.5MPa时,灰度阈值TL=160; 法向应力σ=1.0MPa时,灰度阈值TL=155; 法向应力σ=1.5MPa时,灰度阈值TL=150; 法向应力σ=2.0MPa时,灰度阈值TL=145; 法向应力σ=2.5MPa时,灰度阈值TL=140。在大量数值模拟数据的基础上,经过多元统计分析,拟合出结构面峰值抗剪强度估算新模型,其形式如下, (5)通过27个案例估算出武隆鸡尾山高速远程滑坡滑动面9组黏聚力C和内摩擦角φ建议值,对比已出版文献数据,与各学者采用参数具有相似的统计规律,黏聚力C和内摩擦角φ的最小值组合为:φ=3.977°, C=0.443MPa,最大值组合为:φ=38.080°, C=1.917MPa,平均值组合为：(p=17.162°, C=1.094MPa,工况一(最小值组合)中提供的滑动面黏聚力C和内摩擦角φ比较符合实际情况。以工况一为基础,模拟滑坡破坏初期滑体的位移变化规律,数值模拟结果显示：滑体不同部位处的岩体具有不同的变形规律。滑体前缘和后缘块体具有较大的位移,前缘滑体沿T2裂缝走向滑动,后缘块体沿滑动面和T1裂缝交线方向滑动,滑体前缘块体由于临空面的存在较先发生破坏；相同位置处的滑体上部块体位移量比下部要大,上部块体先于下部块体发生破坏；滑体外侧(东侧)块体位移量要大于内侧位移量,尤其滑体后缘部分差异更为明显,不同步移动造成滑体破坏过程中发生旋转运动。通过与现场地质调查结果对比,证明了数值模拟结果的合理性,同时也证明了BAP结构面峰值抗剪强度估算模型得到了较为成功的应用。本论文主要创新成果有： (1)提出表征岩体结构面粗糙度的新指标BAP。它是一种三维评价方法,可以考虑粗糙度的各向异性、尺寸效应与间距效应等特征,并且结合了岩体结构面的剪切破坏机制；该法基于三维激光扫描试验获取的点云数据,所有计算过程实现了程序化,评价精度和速度得到大幅度提高。 (2)基于三维粗糙度评价新指标BAP,建立岩体结构面峰值抗剪强度估算新模型。该模型是在合理的采样尺寸与采样间距下,考虑了与结构面相关法向应力σ、壁面强度JCS、基本摩擦角φb、粗糙度BAP以及剪切方向等因素而提出的；结构相对来说比较简单,具有很强的力学含义,模型中涉及到的统计参数是在大量数值试验结果基础上拟合而得,具有一定的可靠性；由该模型的实际工程应用显示,适用性较强,能够估算出合理的结构面抗剪力学参数。
[Abstract]:......
【学位授予单位】：中国地质大学
【学位级别】：博士
【学位授予年份】：2014
【分类号】：TU45

【参考文献】