岩体边坡三维模型及块体滑落分析方法技术

一种岩体边坡三维模型及块体滑落分析方法，属于岩土工程、矿山开采、道路工程、水利工程等岩体工程的实体结构空间建模建立和分析领域。按如下步骤进行：步骤一、获取工程岩体现场结构原始数据；步骤二、工程岩体原始数据的处理与提取；步骤三、构建工程岩体模型；步骤四、识别块体及稳定性分析；步骤五、结果显示；本发明专利技术的优点：通过该方法实现了工程岩质边坡的空间模型建立，关键块体的识别，考虑工程实施过程中快速辨别新结构面，并能自动搜索关键块体，统计出可移动块体数量和几何信息，并能初步计算分析关键块体的可动性。

【技术实现步骤摘要】

本专利技术属于岩土工程、矿山开采、道路工程、水利工程等岩体工程的实体结构空间 建模建立和分析领域，特别涉及一种。
技术介绍
边坡稳定性问题一直是岩土工程的一个重要研究内容。随着国民经济建设的蓬勃 发展，铁路、公路、水利水电、矿山等工程的高陡边坡越来越多，崩塌的问题更加突出。为了 防止崩塌地质灾害的发生，有计划、有重点地对边坡进行整治，对边坡崩塌稳定性进行分析 就变得更加迫切。大量的研究表明，崩塌现象的发生绝大部分与岩体结构面有关。结构面是地质历 史发展过程中在岩体内形成的具有一定方向、一定规模、一定形态和特性的地质界面，由于 结构面的存在，不仅破坏了岩体的完整性，而且直接影响岩体的力学性质和应力分布形态， 因此研究结构面对岩体稳定性的影响至关重要。从某种意义上说，能否较合理的分析结构 面对岩体稳定性的影响，成为衡量该方法是否适用于岩体稳定性分析的重要标准。块体理 论以其独特的计算分析思路和功效在岩体稳定分析中占据重要地位。块体理论针对岩体在 结构面切割下的稳定性问题进行研究，其实质是应用几何方法(拓扑学和集合论)研究由 结构面和临空面组成的块体类型及可动性，根据块体的静力平衡条件，求出滑动力，评价它 的稳定性。传统的分析方法是基于有限元、离散元、边界元及其耦合等数值方法考虑，模拟过 程很难考虑岩体结构的真实特征和工程的扰动过程，目前，大多采取块体理论解决上述问 题，但是还存在一些不足，如不定位块体中计算块体的大小、考虑粘聚力作用的块体稳定性 计算等问题。
技术实现思路
为了解决现有技术的不足，本专利技术提供一种岩体边坡三维模型及块体滑落分析方 法，基于现场工程岩体结构特性和地质调查报告，建立反映岩质边坡真实结构特性空间三 维模型，实现了工程岩体结构和施工过程结合起来模拟真实岩体结构，本专利技术实现了块体 搜索与显示的方法。同时，通过对现场数据的对比以及现场观测数据的反馈，进一步分析了 工程岩体破裂与失稳的状态；在坚硬和半坚硬地层中，岩体被结构面切割成各种类型的空 间镶嵌块体，这些块体通常分为关键块体和一般块体，在自然状态下，这些块体空间处于静 力平衡状态。当进行边坡，隧道等人工开挖，或对岩体施加新的荷载后，使暴露在临空面上 的某些块体失去原始的静力平衡状态，因而造成某些块体首先沿着结构面滑移、失稳、进而 产生连锁反应，造成整个岩体工程破坏。我们称这种首先失稳的块体为“关键块体”。本发 明主要就是在建立岩质边坡空间模型之后，搜索出关键块体，然后分析它的稳定性，寻找出 影响整个边坡稳定性的可移动关键块体，为岩质边坡稳定性分析及其施工方案优化、加固 方案设计提供了一条新路径。一种，按如下步骤进行，如附图说明图1所示步骤一、获取工程岩体现场结构原始数据，包括工程岩体结构参数、工程岩体结 构面信息和工程岩体结构状态，其中工程岩体结构参数包括边坡的尺寸；工程岩体结构面 信息和工程岩体结构状态包括断层参数和节理参数。步骤二、工程岩体原始数据的处理与提取；通过岩体结构面的结构特点会产生多种原始数据，同时，在对这些数据进行处理过程 中，也会产生大量的中间数据，如何对这些数据进行科学的管理和有效的处理，以充分利用 和挖掘其中的有效信息，从而保证成果的准确可靠性，是一个关键问题，工程岩体裂隙发育 的特点决定了对工程岩体裂隙的研究必须建立在大量的统计分析基础之上，以实现对研究 区具体工程部位的工程岩体裂隙进行优势方位分析及配套参数提取；在现场采集的原始数据基础上，利用赤平投影和概率统计方法对结构面进行分组，整 理出1)岩体裂隙方位，以及每一条裂隙所在的方位、层位、工程部位、风化带信息，进行裂隙 优势方位分析；2)存储工程普遍测网法获得的每条裂隙的基本信息，包括每条裂隙对应的测点号、位 置和产状，进行绘制裂隙平面展布图的坐标信息；3)存储工程岩体测点的各条裂隙的产状及相对测网原点的局部几何坐标；4)在对普遍测网法实测数据有效管理的基础上，根据现场实测数据对各测点的岩体结 构参数通过蒙特卡罗模拟方法计算结构面随机变量，模拟结构面的密度、迹长、倾向、倾角 参数；结构面分组框架图如图2所示。在此基础上，为本专利技术提供的初始数据，边坡工程结构参数和结构面参数的数据 范围及精度要求根据前述的结构面几何数据对表征岩体结构指标进行有效的测试，包括结构面的优势 方位、间距、密度、连通率、各类迹线长、测网延时质量指标(RQD)值、裂隙系数、三维空间结 构模式、裂隙网络分形维数、损伤张量、渗透张量等十余个岩体结构参数，并具有对选择的 多个测点进行综合统计分析和结构面参数概率分布拟合及裂隙网络模拟计算，包括1)基本分析计算本文档来自技高网...

【技术保护点】
一种岩体边坡三维模型及块体滑落分析方法，按如下步骤进行：　　步骤一、获取工程岩体现场结构原始数据，包括工程岩体结构参数；工程岩体结构状态断层参数和节理参数；　　步骤二、工程岩体原始数据的处理与提取；　　在现场采集的原始数据基础上，利用赤平投影和概率统计方法，对结构面进行分组，整理出：　　１）岩体裂隙方位，以及每一条裂隙所在的方位、层位、工程部位、风化带信息，进行裂隙优势方位分析；　　２）存储工程普遍测网法获得的每条裂隙的基本信息，包括每条裂隙对应的测点号、位置和产状，进行绘制面对块体切割分解为对块体表面多边形的切割，进而分解为裂隙面对直线的切割，实现裂隙与多边形的切割，裂隙切割多边形；　　分有三种情况：　　Ａ如果多边形顶点全在裂隙面的上方，此时把这个多边形编号后放到上面的新块体中；　　Ｂ如果多边形顶点全在裂隙面的下方，此时把这个多边形编号后放到下面的新块体中；　　Ｃ如果多边形顶点分布在裂隙面的两侧，用裂隙平面把这个多边形切割成两个新多边形，编号后分别放入上下两个块体中；　　采用裂隙切割其他块体，直到切割完与裂隙相关的所有块体，然后进入下一个裂隙的切割过程；无论是虚拟结构面还是真实的结构面，在切割的方法中是不区别对待的，结构面把与其接触的块体一分为二，而与其不接触的块体保持原来的状态，依次加入结构面即虚结构面或真实结构面，每添加一个结构面，遍历一次现存的块体，判断块体是否与结构面相交，如果相交就进行切割，把新形成的块体放入块体数据链的尾端，删除母体，继续判断下一块体，不相交则直接进入下一块体的判断；　　２）　添加结构面：网格划分后，加入实测结构面和拟合出来的结构面，对现有的小块体进行再次切割；　　３）、消除网格，生成块体；　　首先判断相邻两个小单元面的交集是否与虚拟结构面有交集，如果有交集则合并两个小块，没有交集就不必合并，按照此方法，依次对存在的所有小块体进行处理，最后记录下合并后块体的数据；　　４）、确定块体的面积及体积　　首先把多边形分解成三角形，然后对各三角形的面积求和，得到多边形的面积；而多面体首先被分解成四面体，然后对所有的四面体体积求和，最后得到多面体的体积；　　如果多边形有ｎ个边，可以分解为ｎ－２个三角形，若三角形的三个顶点坐标分别为：Ｐ↓［１］（ｘ↓［１］，ｙ↓［１］，ｚ↓［１］），Ｐ↓［２］（ｘ↓［２］，ｙ↓［２］，ｚ↓［２］），Ｐ↓［３］（ｘ↓［３］，ｙ↓［３］，ｚ↓［３］），　　三角形面积...

【技术特征摘要】
一种岩体边坡三维模型及块体滑落分析方法，按如下步骤进行步骤一、获取工程岩体现场结构原始数据，包括工程岩体结构参数；工程岩体结构状态断层参数和节理参数；步骤二、工程岩体原始数据的处理与提取；在现场采集的原始数据基础上，利用赤平投影和概率统计方法，对结构面进行分组，整理出1)岩体裂隙方位，以及每一条裂隙所在的方位、层位、工程部位、风化带信息，进行裂隙优势方位分析；2)存储工程普遍测网法获得的每条裂隙的基本信息，包括每条裂隙对应的测点号、位置和产状，进行绘制裂隙平面展布图的坐标信息；3)存储工程岩体测点的各条裂隙的产状及相对测网原点的局部几何坐标；4)在对普遍测网法实测数据有效管理的基础上，根据现场实测数据对各测点的岩体结构参数通过蒙特卡罗模拟方法计算结构面随机变量，模拟结构面的密度、迹长、倾向、倾角参数；步骤三、构建工程岩体模型；1)、建立岩体模型将岩体几何体分成凸体和凹体两类，分别对凸体结构的岩体几何体和凹体结构的岩体几何体进行建模；凸体结构首先，把岩体几何体分解为多边形，依次添加多边形，多边形的建立是按同一顺序即逆时针或顺时针顺序，依次按顺序输入顶点坐标到文本文件中，按照岩体被划分后，任何两个相邻的岩体公共面必须是重合的原则将划分实际工程岩体划分为凸体；凹体结构按照相邻两个岩体几何体公共面重合的原则，把凹体分成多个凸体，并把重合面赋上虚拟面的标记值，将被分成的凸体按照凸体的建模规则建模；2)、建立数据结构先确定块体棱线的两个端点，放入线段存储结构中，再把形成的线段放入多边形存储结构中，最后把多边形放入块体的数据结构中；采用编程语言中自带的动态类对数据进行管理；3)、确定结构面信息的方法采用确定性结构面和非确定性结构面的方法确定性结构面方法通过步骤一得到采集的结构面数据，作为确定性结构面数据，倾向、倾角和迹长；非确定性结构面通过现场调查得到的结构面样本的信息，以蒙特卡罗方法为基本理论，假设块体在岩体中的分布规律，采用计算机自动生成随机数，并进一步根据样本信息模拟结构面的空间分布及力学参数，然后调用系统确定性分析模块数据库，自动搜索所有的随机可动块体；其特征在于步骤四、识别块体及稳定性分析把岩体模型划分为有限个网格，然后在依次加入结构面，对现有的小块体进行切割，最后再去掉网格，合并小的块体，识别出结构面切割得到的所有块体；具体方法如下1)、网格划分根据区域内结构面的密度，平均半径进行网格划分，首先确定区域的范围，然后设置三组相互垂直的虚拟结构面对岩体模型进行切割，最后记录下小单元的数据结构，网格的大小根据结构面的平均尺寸，密度而定，采用结构面的正交模型对岩体模型进行切割，虚拟结构面的间距即单元的边长，其半径大于模型尺寸；具体切割过程如下 切割过程为裂隙面对块体切割分解为对块体表面多边形的切割，进而分解为裂隙面对直线的切割，实现裂隙与多边形的切割，裂隙切割多边形；分有三种情况A如果多边形顶点全在裂隙面的上方，此时把这个多边形编号后放到上面的新块体中；B如果多边形顶点全在裂隙面的下方，此时把这个多边形编号后放到下面的新块体中；C如果多边形顶点分布在裂隙面的两侧，用裂隙平面把这个多边形切割成两个新多边形，编号后分别放入上下两个块体中；采用裂隙切割其他块体，直到切割完与裂隙相关的所有块体，然后进入下一个裂隙的切割过程；无论是虚拟结构面还是真实的结构面，在切割的方法中是不区别对待的，结构面把与其接触的块体一分为二，而与其不接触的块体保持原来的状态，依次加入结构面即虚结构面或真实结构面，每添加一个结构面，遍历一次现存的块体，判断块体是否与结构面相交，如果相交就进行切割，把新形成的块体放入块体数据链的尾端，删除母体，继续判断下一块体，不相交则直接进入下一块体的判断；2) 添加结构面网格划分后，加入实测结构面和拟合出来的结构面，对现有的小块体进...

【专利技术属性】
技术研发人员：王述红，张靖杰，郭牡丹，张娟霞，金日守，关永平，张航，杨昊，
申请(专利权)人：东北大学，
类型：发明
国别省市：89[中国|沈阳]