1、2020/1/8,岩体结构数字识别与力学参数表征,于庆磊 东北大学,2,研究背景与现状,提 纲,1,3,1、研究背景与现状,岩体的结构特性岩体是经过地质作用过的、由岩块和地质 结构面（节理、裂隙、断层等）所组成、 具有一定结构特征，赋存于一定地质环境 （地应力、地下水、地温等）中的地质体。,不同规模的节理、裂隙、断层等结构面,4,1、研究背景与现状,在结构面的五级分类方案中，采矿、水利及交通等工程岩体开挖过程中揭露并需要现场测量的结构面主要包括三级： 级：小规模的次级断层，规模较大，数量较少，和开挖临空面相互交切组合，主要制约局部危岩体的崩落、滑动方式 级：原生的具有一定分布规律和统计意义的节
2、理面、卸荷裂隙、风化裂隙等 级：连接性极差的微小节理或隐微裂面 级和级数量大、随机展布，主要弱化岩体的力学性质,5,不同级别的结构面的空间分布和成组特征，直接影响岩体的工程特性，导致了岩体工程特性参数如变形、渗透性、强度等的各向异性、非均匀性及离散性，节理岩体特性参数的确定一直是计算岩体力学的难点,1、研究背景与现状,6,1、研究背景与现状,岩体的力学特性,7,岩体渗流力学研究方法大结构：主干裂隙系统数量少、规模大、起主导控制作用采用离散网络方法 小裂隙：裂隙岩块系统数目多、分布广、提供较大渗流空间采用连续介质方法 从宏观上来描述其对渗流场的贡献 采用张量分析计算方法 然后，通过两系统间的水量
3、交换建立联系,1、研究背景与现状,8,1、研究背景与现状,岩体力学、 水力学参数,数学模型数值计算 方法,求解结果,岩体力学问题求解,不可靠,无论是哪种模型，其根本的问题是如何确定参数（张有天）,9,1、研究背景与现状,数值计算在岩石工程应用中一直存在“声誉高、信誉低”的问题，制约其在岩石工程应用的技术“瓶颈”问题为：,数据不完备的复杂地质系统与理论严密的力学模 型之间相互脱节，岩体结构地学描述的多尺度定 性方法与力学分析的精细定量方法衔接不好,岩体参数和数值计算的可靠性取决于对 岩体结构面测量和表征的完备程度及数值模型中对控制结构面 的描述程度,关键,测量和表征岩体结构面,10,结构面测量现
4、状,精测线法和窗口法 传统方法：即通过皮尺和罗盘人工现场逐一接触测量结构面信息 缺点：采集的信息有限，有些高陡岩体不可能全面接触，使得测量数据的代表性受现场条件的限制.低效、费力、耗时，不安全、甚至难以接近实体和不能满足现代快速施工的要求,1、研究背景与现状,11,钻孔定向取芯技术和孔内照相技术 缺点：获取岩体结构面信息规模小、应用效果不佳， 该方法可作为工程地质调查的辅助手段,1、研究背景与现状,12,岩体结构参数三维非接触测量系统3GSM 数字摄影测量技术,从奥地利startup公司引进的一套3G软件和测量产品JointMetriX3D和ShapeMetriX3D是一个全新的、代表当今最高
5、水平的岩体几何参数三维不接触测量系统 该系统由一个可以进行高分辨率立体摄像的照相机、进行三维图像生成的模型重建软件和对三维图像进行交互式空间可视化分析的分析软件包组成。,jointmetrix3d成像系统,shapemetrix3d成像系统,1、研究背景与现状,13,通过标定的高分辨率的图像照相扫描系统，获取包括: 岩体不连续面的几何参数如间距、面的连通率和方向等信息。由此可推导出岩体的标识参数、所有的分析都是三维、不接触并以出处的坐标给出。 对该信息进行数字化建模，建立岩体三维力学模型，在联想深腾1800并行计算机上进行渗流、应力等物理过程计算，得到结果导入GOCAD系统中，在我们购进的大型
6、虚拟现实系统中演示，有力的支撑了数值模拟和虚拟现实过程研究，并使得该中心具有更好的工程应用前景。,14,3GSM测量原理,立体图像合成原理：两个图像上相应点P(u,v)组成三维空间物体点P(x,y,z)（露天边坡）,软件系统对不同角度的图像进行一系列的技术处理如基准标定、像素点匹配、图像变形偏差纠正等，实现物体表面真三维模型重构，在计算机可视化屏幕上从任何方位观察三维实体图像,15,3GSM的软件系统,16,3GSM的软件系统,使用电脑鼠标进行交互式操作来实现每个结构面个体的识别、定位、拟合、追踪以及几何形态信息参数（产状、迹长、间距、断距等）的获取，并进行纷繁复杂结构面的分级、分组、几何参数
7、统计。,17,数字摄影测量3GSM技术创建了一个实时的地质信息交流和反馈 环境，提高地质纪录任务的效率，降低不完整信息和信息丢失 的可能性，大大的帮助地质工作者区分鉴定地质特征，在已完 成工作面节理图像的基础上预测没完成工作面上弱面的位置和 方向。但是，该技术不和数值力学计算结合起来，其工作潜力 还远远没有发挥出来。,利用详实的结构面信息，开展力学分析与计算,1、研究背景与现状,18,研究思路,测量数据不完备的复杂岩体结构系统与理论严密的力学模型之间相互脱节是岩体力学应用的关键问题。,19,借助于先进的3GSM三维岩体不接触测量技术，获取一系列真实描述岩体宏观结构的数字图像，提取节理几何形态空
8、间分布信息；在此基础上，对于级别较低的节理面，基于考虑节理统计分布的几何损伤理论、Hoek-Brown系统和辅以现场试验数据修正的裂隙样本法，建立岩体力学、水力学参数表征方法；对级别较高的结构面按节理单元处理，实现岩体结构参数（几何形态和力学性质）数字信息与力学分析的精细定量计算有机衔接，建立能反映控制结构面真实空间分布、充分考虑各级别结构面作用的三维岩体破坏过程分析数值模型，为岩体工程结构稳定性分析和破坏机理研究提供更加真实、可靠的手段。,研究思路,20,2、结构面数字摄影测量与三维形态重构,结构面数字摄影测量,合成三维图,21,2、结构面数字摄影测量与三维形态重构,边坡岩体中节理分布情况,
9、第一组节理参数统计,第二组节理参数统计,两组节理赤平极射投影,22,2、结构面数字摄影测量与三维形态重构,第一组结构面几何参数(部分),第二组结构面几何参数(部分),23,第一组,第二组,迹长(负指数),倾向(正态),倾角(对数正态),迹长(对数正态),倾向(正态),倾角(正态),2、结构面数字摄影测量与三维形态重构,24,2、结构面数字摄影测量与三维形态重构,结构面几何参数概率模型,25,2、结构面数字摄影测量与三维形态重构,一组结构面的间距s和断距d的确定,26,2、结构面数字摄影测量与三维形态重构,三维形态重构,在3GSM三维岩体不接触测量数字图像系统基础上，对空间分布不同级别结构面的进
10、行准确识别，实现岩体结构几何形态空间分布的数字表征，利用可视化编程平台Visual C+和AUTOCAD08系统开发结构面空间展布可视化三维模型（AUTOCAD格式），得到不同级别、组别结构面几何形态参数。,27,28,2、结构面数字摄影测量与三维形态重构,生成的101010 m3三维裂隙内部分布,生成的101010 m3裂隙网络,对于统计型节理，假定节理面形状为圆盘状，即Baecher圆盘模型，利用Mont-Carlo方法重构三维裂隙网络,29,2、结构面数字摄影测量与三维形态重构,A-A剖面,B-B剖面,30,2、结构面数字摄影测量与三维形态重构,31,3、岩体力学、水力学参数表征,岩体力
11、学参数表征,现场原位试验法：获得参数最为准确，但时间长、费用高 等效折减法：对于工程岩体参数表征和力学分析计算，一般基于传统的测线法采集岩体结构面信息数据，进行岩体参数的等效折减：采用岩体完整性系数修正法、费先科方法、Gergi方法、分形维理论、蒙特卡洛方法等统计方法并结合RMR系统、Q系统和Hoek-Brown系统，同时参考室内岩块强度指标来确定岩体强度指标。,32,3、岩体力学、水力学参数表征,岩体力学参数表征,基于几何损伤力学理论确定岩体的初始损伤张量，由此确定岩体初始弹性模量 基于Hoek-Brown公式的m，s参数（反映岩体软硬程度和结构面分布密度）和GSI（地质体强度指数），结合详
12、实的岩体结构测量信息，定量估计岩体强度,对于结构面（原生统计结构面）和级结构面（开挖扰动次生结构面），其平均间距小于研究区域模型边界尺寸的1/20时，按照等效连续介质考虑：,33,岩体初始弹性模量,由几何损伤张量理论，含节理岩体的初始损伤由下式计算：,式中，l 节理面最小间距；V 研究的含节理岩体体积；N 节理的数量；ak 第K条节理的面积；nk 第K条节理的法矢量；,通过上式计算三维条件下裂隙岩体损伤张量，根据式,Em=（I-D）Er,式中：Em岩体初始弹性模量；Er岩块弹性模量。,34,岩体强度定量估算,推广后的Hoek-Brown准则为：,式中： 岩体破坏时的最大主应力岩体破坏时的最小主
13、应力；组成岩体完整岩块的单轴抗压强度；岩体的Hoek-Brown常量；组成岩体的完整岩块的Hoek-Brown常数。,35,对于GSI25的岩体：,对于GSI25的岩体：,GSI为岩体的地质强度指标。,岩体强度定量估算,36,岩体强度定量估算,表1 由岩石类型所决定的Hoek-Brown常量mI,表2岩体结构特征定量描述的Jv表示,37,38,岩体渗透系数确定,现场测试法 单孔压水试验、三段压水试验、修正渗透张量压水试验、三孔交叉压水试验 由于裂隙岩体渗透性的离散程度大，现场测试法的结果具有尺寸效应 单孔压水试验得出的渗透系数是无方向性的标量，其理论基础是各向同性的孔隙介质，故所得K值不能反映
14、裂隙岩体渗透性的非均匀性和各向异性 三段压水试验、修正渗透张量的压水试验等方法，试验成本较高且均具有自身的局限性,应用广泛的是现场测试法和裂隙样本法,39,岩体渗透系数确定,裂隙样本法 现场收集和统计研究岩区的裂隙网络的倾向、倾角、开度、分布、粗糙度、连通度等，然后依据单裂隙渗流理论和张量分析方法获得裂隙网络的渗透系数及主渗透方向等 裂隙样本法在理论上己比较完善，不足之处是由于裂隙的隙宽、隙壁粗糙度、填充物以及连通度难以准确测量而未能考虑这些因素的影响,两者相结合，采用裂隙样本法确定岩体的渗透张量，并辅以现场测试法的数据进行修正，既能较好反映原位地质环境裂隙岩体渗透性又能消除压水试验尺寸效应，
15、克服单孔压水试验不能反映渗透各向异性的不足 该技术途径不但适合工程岩体宏观渗流分析的特点，而且能够较为深刻地认识裂隙渗流的本质。,40,裂隙样本法确定渗透张量,渗流区域裂隙网络示意图,如图所示，二维渗流区域MN和MN为定水头边界，MM和NN为零流量边界，由下式得其裂隙网络渗流数学模型为,式中，未知量为内节点和零流量边界节点的水头H；方程式中q包含有H项，其中第j线元的流量由单一裂隙渗流公式求得：,MN，MN之间的等效渗透系数为：,41,裂隙样本法确定渗透张量,对每隔15度方向的单元网络施加上图所示的水头边界条件，并用流量等效原则计算出水流流入（或流出）研究区域内的流量，由式可得网络整体的渗透系
16、数。据对称性，一般选取6个方向的研究网络，可得不同方向的渗透系数，进而求出网络单元的渗透张量。,渗透系数张量图,42,范各庄应用实例,开滦范各庄煤矿是我国自行设计、施工的一座大型现代化矿井，矿区的水文地质条件极其复杂。含煤地层属石炭二叠系，主要可采煤层为下二叠统大苗庄组的5煤层、7煤层、9煤层和上石炭统的赵各庄组12煤层。其中12煤层位于本区煤系地层的下部，为复杂结构的厚煤层，煤层厚度1.058.32m，平均厚度为3.54m，在煤炭开采过程中底板砂岩破坏诱发突水事故频繁。根据开滦范各庄煤矿的复杂地质条件，12煤底板含水砂岩层内含有平均压力为2MPa的承压水，精确查明该岩体的力学、水力学参数对于
17、分析底板含水砂岩层破断突水机理具有重要的科学意义。,范各庄煤矿砂岩岩体条件,43,岩体结构测量,本次采用3GSM设备于08年10月7日在-530m水平12煤底板砂岩巷道中进行。,左视图,右视图,44,系统根据像素点对应合成得到掌子面三维表面图,在合成的三维图上，根据主要的节理裂隙的分布情况，对其进行分组，不同颜色代表不同的组，主要的裂隙分布有两组，即右图中的黄色组和红色组。,岩体结构测量,45,倾伏角概率模型统计表,节理面概率模型统计,46,评价范各庄煤矿12煤底板突水危险性，需要分析砂岩岩体下部奥陶系灰岩层中承压水沿砂岩岩体垂直方向的裂隙向上渗透，计算砂岩岩体垂直方向的渗透系数十分关键。所以
18、，为了简化计算，建立砂岩岩体垂直方向二维X-Z方向裂隙网络模型，开展X-Z方向渗透张量分析研究。,模拟生成裂隙网络的区域为16m16m 选取测定渗透张量的区域边长分别为 3m、5m、8m、9m、10m 固定区域的中心点，逆时针旋转矩形，每隔15计算一次裂隙网络的渗透系数，最终求得的渗透张量,范各庄砂岩裂隙网络模拟及渗透张量分析,47,范各庄砂岩裂隙网络模拟及渗透张量分析,0,0,45,45,完整裂隙网络,连通裂隙网络,48,渗透张量（10-6 m/s）,右图直观地反映出网络渗透特性受控于结构面几何分布规律。计算结果表明，随着研究域的增大，渗透系数主值逐渐减小，当研究区域为9m左右时，渗透张量的
19、主值已很接近，9m与10m计算的各方向的渗透系数的最大误差为5.79%。可见，10m10m可作为渗透张量的表征单元体，得到的最大渗透系数Kmax=1.4210-6 m/s，最小渗透系数Kmin=6.6710-7 m/s，方向角为45，Kmax/Kmin=2.13。,不同尺度时的裂隙渗透系数主值,范各庄砂岩裂隙网络模拟及渗透张量分析,49,当用Hoek-Brown准则估计节理化岩体强度指标与力学参数时，需用3个基本参数 (1) 组成岩体的完整岩块的单轴抗压强度c ； (2) 组成岩体的完整岩块的Hoek-Brown常数m1（查表） (3) 岩体的地质强度指标GSI（查表）,范各庄煤矿砂岩节理化岩
20、体强度计算,50,完整岩块的单轴抗压强度=175MPa，掌子面面积6.2m2；节理条数：22；Jv=4，GSI=60，砂岩m1=19 由 得 mb=4.553,摩擦角 40.195； 内聚力 c=10.148MPa 岩体的单轴抗拉强度 tm= -0.451MPa 节理化岩体抗压强度 cm=43.717MPa,范各庄煤矿砂岩节理化岩体强度计算,最终求得：,51,渗流场计算,整个模型由泥岩、12煤、隔水泥岩层、含水砂岩层、奥陶系灰岩层组成。,52,渗流场计算,53,渗流场计算,=0水压力和速度场矢量,=18水压力和速度场矢量,54,渗流场计算,=90水压力和速度场矢量,=144水压力和速度场矢量,
21、55,5、前景与展望,岩体参数和数值计算的可靠性取决于对岩体结构面测量和表征的完备程度及数值模型中对控制结构面的描述程度,3GSM为检测岩体结构提供了有效工具,根据岩结构与工程规模，控制结构面以节理单元表示，统计结构面的作用在参数表征过程中体现，建立充分反映各级工程岩体结构面作用的数值模型,借助大规模并行计算工具计算岩体稳定性,56,2008-06-10,汇报结束 谢谢大家！,57,什么是博士？,58,假设人类所有的知识，就是一个圆。圆的内部代表已知，圆的外部代表未知,59,读完小学，你有了一些最基本的知识,60,读完中学，你的知识又多了一点,61,读完本科，你不仅有了更多的知识，而 且还有了
22、一个专业方向,62,读完硕士，你在专业上又前进了一大步,63,进入博士生阶段，你大量阅读文献，接 触到本专业的最前沿,64,你选择边界上的一个点，也就是一个非常 专门的问题，作为自己的主攻方向,65,你在这个点上苦苦思索，也许需要好几年,66,终于有一天，你突破了这个点,67,你把人类的知识向前推进了一步，这 时你就成为博士了,68,现在你就是最前沿，其他人都在你身后,69,但是，不要陶醉在这个点上，不要把整 张图的样子忘了,70,美国犹他大学的助理教授Matt Might，用这组图解释，博士学位到底是什么意思。他说，每年都有新生的入学教育，但是有些观点语言说不清楚，不如画图。 老子说“大道至简“，可是真的要很简单地表达出来，却是非常难的一件事。,71,
