CN112132409A - 一种基于BQ和RQDt的巷道开挖不利方位求解方法 - Google Patents
一种基于BQ和RQDt的巷道开挖不利方位求解方法 Download PDFInfo
- Publication number
- CN112132409A CN112132409A CN202010892992.9A CN202010892992A CN112132409A CN 112132409 A CN112132409 A CN 112132409A CN 202010892992 A CN202010892992 A CN 202010892992A CN 112132409 A CN112132409 A CN 112132409A
- Authority
- CN
- China
- Prior art keywords
- rqd
- solving
- value
- network model
- follows
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
- G06F7/58—Random or pseudo-random number generators
- G06F7/588—Random number generators, i.e. based on natural stochastic processes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06395—Quality analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Geometry (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Economics (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Optimization (AREA)
- Marketing (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Computation (AREA)
- Game Theory and Decision Science (AREA)
- General Business, Economics & Management (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Architecture (AREA)
- Structural Engineering (AREA)
- Civil Engineering (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于BQ和RQDt的巷道开挖不利方位求解方法,属于巷道开挖方位领域,其步骤包括:(1)基于BQ指标的岩体质量计算;(2)岩体三维裂隙网络模型生成和剖切;(3)RQDt各向异性图绘制;(4)基于BQ反演的最佳阈值t求解方法;(5)空间RQDt求解方法;(6)巷道开挖不利方位求解方法。本发明将BQ理论、裂隙网络模型、广义RQD理论和圆锥曲线方程相结合,实现了RQDt的最佳阈值t求解、RQDt椭圆曲线和各向异性椭球体绘制,实现了最佳阈值t下的巷道开挖不利方位求解。本发明方法明确,所得最佳阈值t、椭圆曲线和各向异性椭球体能够准确反映出岩体的RQDt各向异性特征,得到的巷道开挖不利方位可用于指导巷道开挖施工应用。
Description
技术领域
本发明涉及一种基于BQ和RQDt的巷道开挖不利方位求解方法,特别的是本发明基于BQ指标、裂隙网络模型、广义RQD理论和圆锥曲线方程,反演计算出RQDt最佳阈值t的范围和最佳阈值t值,得到了最佳阈值t下的RQDt椭圆曲线和各向异性椭球体,得到了最佳阈值t下的巷道开挖不利方位,提供了一种基于BQ和RQDt的巷道开挖不利方位求解方法,属于巷道开挖方位领域。
背景技术
在自然界中,各向异性无处不在。岩体表现出非常明显的各向异性,其性质随观察角度和测量方向的变化而变化。岩体各向异性变化的起源可归结于地质起源,如沉积岩中独特的层理、玄武岩中独具的天然孔洞、片麻岩中的片麻状构造和板岩中的劈理,以及贯穿这些岩体的断层等。除了固有的地质起源造成的岩体各向异性外,一些其它因素也会影响岩体的各向异性,如岩体中节理组的产状和分布形态,岩体中的节理会造成岩体的应力各向异性,影响岩体的质量。
岩体的外观和质量强烈的依赖于岩体的各向异性程度,岩体的稳定性是由岩体质量所决定的。一些地质现象也已经证明了这一点。例如,在巷道或高速公路隧道开挖过程中,会沿着主要节理组倾角出现非垂直性的变形破坏,或者隧道中的某一块岩体出现楔形块体掉落,进而造成整个岩体的失稳或者巷道的坍塌,这种失稳或坍塌对巷道开挖是非常不利的,不仅影响工程施工进度,更重要的是会造成严重的人员伤害或者财产安全。因此,在巷道开采的前期,如何准确的获取岩体质量的各向异性程度,以及如何依据岩体质量的各向异性程度,选择巷道的开采方向和开挖走向,具有重要的研究价值和现实指导意义。而目前国内外的研究中,从岩体质量各向异性角度考虑巷道开挖方向的研究,尚不多见。
RQD是表征岩体质量的重要指标,RQD具有各向异性。如采用钻孔的方式获得RQD,不同的岩体部位得到的结果完全不同,RQD值结果将取决于方向,当钻孔方向与主要节理组平行时,得到较高的RQD值,当钻孔方向垂直于主要节理组方向时,会得到较低的RQD值。RQD的各向异性,从钻孔RQD到广义RQD的计算公式的拓展也可以看出。Deer在1964年提出了钻孔RQD的概念,由于钻孔RQD在应用时存在着如下两个缺点:对于不同工程规模岩体,100mm的阈值选取是否合理;钻孔钻探方向具有局限性,获得的RQD只能反映局部岩体情况,不能反映出RQD各向异性。因此,部分学者引入阈值t,提出了广义RQD的概念,即对任一间距阈值t,把沿某一测线方向大于t的间距之和与测线总长之比的百分数定义为广义RQD,用RQDt表示。广义RQD的引入,让RQD各向异性的求解成为了可能。
由RQDt概念可以知道,阈值t是广义RQD的一项特别重要的参数,是广义RQD能否真实反映出岩体质量的关键,但是阈值t具有任意性,不具有唯一性,因此寻找到能够表征出岩体质量的最佳阈值t,将具有非常重要的科学意义和研究意义。
国内外学者针对阈值t的研究,主要体现在以下两个方面:不同阈值t下的RQDt计算和最佳阈值t的研究。在不同阈值t下的RQDt计算方面,已有的研究主要是探讨阈值t的变化对RQDt值的影响,主要是通过模拟岩体裂隙网络模型,布置虚拟钻孔,以求解不同RQDt值,研究的载体是裂隙网络模型。但在这个方面的研究中,主要是为了研究RQDt的各向异性特征,对最佳阈值t的研究鲜有涉及。
在最佳阈值t的研究方面,有些学者也曾对最佳阈值t展开过探讨。如有的学者基于三维结构面网络模拟技术,运用分形理论分别计算其结构面分布的分形维数,通过不断改变广义RQD的阈值,得到不同阈值下的RQDt值,将各剖面的分形维数与广义RQD值对比分析,为准确选取广义RQD的最佳阈值提供依据。有的学者基于修改后的块度指数MBi,建立了35种假想的三维裂隙网络模型,测量了具有不同阈值的广义RQD值,确定广义RQD的最佳阈值。这两个方面的最佳阈值t的研究,适用条件都有一定的局限性,都是在特定理论和背景下的一种最佳阈值t求解方法,当背景或者模型变化了之后,其最佳阈值t将不再具有最佳性。而且,由于分形维数或者块度指数本身并不具备表征岩体质量的功能,通过分形维数或者块度指数获得的最佳阈值t,能否反映真实的岩体质量,还有值得商榷的地方。
因此,找到并提出一种最能反映出岩体质量的RQDt最佳阈值t的求解方法,具有非常重要的科学和研究意义,也是RQDt阈值t研究中亟需解决的问题。
RQD的各向异性直接影响岩体质量,RQD各向异性对岩体质量的影响机理,目前还尚未探索清楚。由于在阈值t方面,还尚未有学者给出最佳阈值t的求解计算方法,因此众多针对RQDt空间效应的研究,主要还是集中在三维RQDt的计算上,如开展不同阈值t下的三维RQDt计算,或者在此基础上研究不同方向不同位置上的RQDt值,或者建立三维空间云图。很少有学者针对最佳阈值t,开展最能反映岩体质量的RQDt空间效应研究。因此,得到最佳阈值t下的空间RQDt的求解方法具有重要的理论意义和研究意义。
鉴于此,本发明提出了一种基于BQ和RQDt的巷道开挖不利方位求解方法。
发明内容
为了实现巷道开挖的不利方位求解,本发明提供了一种基于BQ和RQDt的巷道开挖不利方位求解方法。本发明基于BQ指标、裂隙网络模型、广义RQD理论和圆锥曲线方程,反演计算出RQDt最佳阈值t的范围和最佳阈值t值,得到了最佳阈值t下的RQDt椭圆曲线和各向异性椭球体,得到了最佳阈值t下的巷道开挖不利方位,提供了一种基于BQ和RQDt的巷道开挖不利方位求解方法
为了解决上述技术问题,本发明提供如下的技术方案:
一种基于BQ和RQDt的巷道开挖不利方位求解方法,所述方法包括以下步骤:
(1)基于BQ指标的岩体质量计算,过程如下:
1.1:根据结构面参数计算岩体完整性系数,公式如下:
式中:Jv为岩体体积节理数,单位条/m3;
1.2:Jv计算公式如下:
式中:L1,L2,...,Ln为垂直于结构面测线长度;N1,N2,...,Nn为同组结构面数目;
1.3:根据岩石单轴抗压强度值和岩体完整性系数值,计算BQ值:
BQ=90+3Rc+250Kv
式中:Rc是岩石单轴抗压强度;Kv为岩体完整性系数;
1.4:在应用BQ计算公式过程中,遵循以下条件:
当Rc>90Kv+30时,以Rc=90Kv+30和Kv代入计算BQ值;
当Kv>0.04Rc+0.4时,以Kv=0.04Rc+0.4和Rc代入计算BQ值;
1.5:根据地下水、主要软弱结构面产状和天然应力影响对BQ进行修正,修正公式如下:
[BQ]=BQ-100(K1+K2+K3)
式中:K1为地下水影响修正系数;K2为软件结构面产状影响修正系数;K3为天然应力影响修正系数;
1.6:得到修正后的BQ岩体分级结果;
(2)岩体三维裂隙网络模型生成和剖切;
(3)RQDt各向异性图绘制;
(4)基于BQ反演的最佳阈值t求解方法;
(5)空间RQDt求解方法;
(6)巷道开挖不利方位求解方法。
进一步,所述步骤(6)中,巷道开挖不利方位求解方法的过程如下:
6.1:平面上巷道开挖不利方位求解方法,过程如下:
6.1.1:根据RQDt椭圆曲线,求解出曲线上最大值RQDtmax和最小值RQDtmin的值;
6.1.2:求解出RQDt的各向异性比值,各向异性比计算公式如下:
式中:RQDtmax为椭圆曲线上的RQDt最大值,RQDtmin为椭圆曲线上的RQDt最小值;
6.1.3:各向异性比越大,岩体的不利方位越明显;
6.1.4:判断巷道开挖方向与RQDtmin方向的关系,当巷道轴线与RQDtmin方向一致时,对围岩稳定非常不利;当巷道轴线与RQDtmin方向垂直时,对围岩稳定有利;
6.1.5:得到巷道开挖的不利方位和最优方位;
6.2:空间上巷道开挖不利方位求解方法,过程如下:
6.2.1:根据RQDt各向异性椭球体,求解出椭球体上最大值RQDtmax和最小值RQDtmin的值;
6.2.2:求解出RQDt的空间各向异性比值,空间各向异性比计算公式如下:
式中:RQDtmax为椭球体上的RQDt最大值,RQDtmin为椭球体上的RQDt最小值;
6.2.3:各向异性比越大,岩体的不利方位越明显;
6.2.4:判断巷道开挖方向与RQDtmin方向的关系,当巷道轴线与RQDtmin方向一致时,对围岩稳定非常不利;当巷道轴线与RQDtmin方向垂直时,对围岩稳定有利;
6.2.5:得到巷道开挖的空间不利方位和最优方位。
进一步,所述步骤(5)中,空间RQDt求解方法的过程如下:
5.1:平面RQDt求解方法
根据最佳阈值t下的36个方位的RQDt值,求解出任意角度的平面RQDt值,过程如下:
5.1.1:求解出最佳阈值t下的36个RQDt值;
5.1.2:以RQDt值沿x方向的分量RQDx作为横坐标,沿y轴方向的分量RQDy作为纵坐标,设定拟合圆锥曲线方程为:
式中:p1,p2,p3为拟合系数;
5.1.3:基于最小二乘法,应用Matlab软件编写拟合程序;
5.1.4:拟合得到RQDt的椭圆曲线;
5.1.5:通过该拟合椭圆曲线,可以得到平面上任意角度的RQDt值;
5.2:空间RQDt求解方法
根据拟合得到的RQDt椭圆曲线,求解出空间内任意方位任意角度的RQDt值,过程如下:
5.2.1:根据拟合得到的RQDt椭圆曲线,求解出椭圆的长轴所在位置;
5.2.2:以椭圆曲线的长轴作为旋转轴,利用Matlab软件编写程序,将椭圆曲线沿长轴旋转360°,得到RQDt的各向异性椭球体;
5.2.3:得到空间内任意方位任意角度的RQDt值。
进一步,所述步骤(4)中,基于BQ反演的最佳阈值t求解方法的过程如下:
4.1:基于BQ指标反演RQDt范围
结合BQ分级计算出的岩体质量级别,查找《岩石质量指标》表,反演确定出该岩体级别下,RQDt范围值;
4.2:最佳阈值t求解方法
4.2.1:在三维裂隙网络模型上,过中心点O,以任意角度剖切三个剖面,得到三个二维裂隙网络模型,导出二维裂隙网络模型及数据;
4.2.2:针对每一个二维裂隙网络模型,设置不同的阈值t,求解出不同阈值t下的RQDt值;
4.2.3:导出不同阈值t下的RQDt值,计算出RQDt均值;
4.2.4:以阈值t为横坐标,以RQDt的均值为纵坐标,绘制RQDt随阈值t变化的散点图;
4.2.5:根据散点图,设置拟合方程,拟合RQDt随阈值t变化的曲线图;
4.2.6:将反演出的RQDt范围值,带入到拟合出的RQDt随阈值t变化的曲线图中,结合函数方程和曲线图,求解出在该RQDt范围内阈值t的范围,共得到三组阈值t的范围;
4.2.7:针对三组阈值t的范围,取其范围的交集,作为最佳阈值t的范围;
4.2.8:以最佳阈值t范围的中点值作为最佳阈值t值,得到最佳阈值t;
4.2.9:输出最佳阈值t的范围和最佳阈值t值。
进一步,所述步骤(3)中,RQDt各向异性图绘制的过程如下:
3.1:RQDt求解计算
基于RQDt理论,确定出二维裂隙网络模型的中心点O,以每隔10°角度做射线,穿过岩体裂隙网络模型,求解出岩体36个方位的RQDt值,过程如下:
3.1.1:RQDt理论公式如下:
式中:xi表示沿某一测线方向的第i个大于给定阈值t的整段岩石或间距长度,RQDt代表对应阈值t的岩石质量指标,即阈值t下的RQDt值;
3.1.2确定二维裂隙网络模型的剖面中心点O、长度a和宽度b,以模型的左下角为坐标原点,水平向右为x轴,垂直向上为y轴,则中心点O坐标为:
边界方程为:
3.1.3:以O为起点,每隔α=10°角绘制一条测线,与裂隙网络模型相交,共绘制36条测线,测线长度L等于O点到裂隙网络模型边界的距离,用L0~L35表示,则测线方程为:
式中:s表示测线,α表示角度;
3.1.4:判断测线与边界的交点,设测线方程与边界方程的交点为(xa,ya),将测线方程与边界方程依次连立,判断测线是否与边界相交,原理如下:
3.1.5:求出测线与裂隙网络模型的边界方程交点(XZ,α,YZ,α);
3.1.6:确定测线所在区间,原理如下:
3.1.7:根据裂隙网络模型中每条节理的起点坐标(Xb,Yb)和终点坐标(Xc,Yc),建立相应的解析方程,定义节理方程如下:
3.1.8:求解第一条测线与各节理方程的交点,循环判断每个交点(Xj,Yj)范围,如果交点符合a<Xj<c且b<Yj<d,则记录该交点,遍历所有节理方程求出所有m个交点;
3.1.9:将记录的m个交点及测线起点坐标和终点坐标,按横坐标或纵坐标从小到大排序;
3.1.10:计算相邻交点的距离,公式如下:
x0=a/2
y0=b/2
xm+1=XZ,α
ym+1=YZ,α
3.1.11:输入一个阈值t;
3.1.12:循环比较di和t的大小,设初始lt=0,规则如下:
lt=lt+di,若di>t
3.1.13:求解每条测线对应的RQDt值,用mα表示,以及对应的测线起点与终点的距离lα,公式如下:
3.1.14:循环求出每条测线所对应的mα,获得36个方向上的RQDt值;
3.2:RQDt各向异性图绘制
根据求解出的36个RQDt值,绘制出RQDt的各向异性图,过程如下:
3.2.1:将36个RQDt值,按角度顺序依次排序;
3.2.3:依次连接36组射线的端点,若某一条射线上的RQDt为0,则取圆心;
3.2.4:绘制出RQDt的各向异性图;
3.2.5:输出RQDt各向异性图;
3.3:不同阈值t下RQDt各向异性图绘制
用于绘制不同阈值t下的岩体RQDt各向异性图,过程如下:
3.3.1:输入不同的阈值t,求解出对应的RQDt值;
3.3.2:将不同阈值t下的RQDt值,以角度为变量,以RQDt值为函数,绘制到相同坐标系下;
3.3.3:得到不同阈值t下的RQDt各向异性图;
3.3.4:输出不同阈值t下的RQDt各向异性图;
3.4:输出RQDt的值。
再进一步,所述步骤(2)中,岩体三维裂隙网络模型生成和剖切的过程如下:
2.1:随机数求解
基于蒙特卡洛方法和节理分布形态,求解随机数,过程如下;
2.1.1:伪随机数的产生,产生随机数的数学方法应满足以下条件:产生的随机数列应均匀分布在(0,1)区间;序列之间应无相关性;随机序列有足够长的重复周期,在计算机上产生的速度快,占有的内存空间小,具有完全可重复性;
2.1.2:给定分布下随机变量数值的产生,Monte Carlo方法是根据确立的结构面几何概率模型,再现服从这种模型的结构面网络模型,该方法的核心是抽样数的随机性,高品质随机数可以得到良好的模拟结果,即在(0,1)区间上生成均匀分布随机变量ri,利用这些均匀随机变量产生服从其他分布的随机数;
2.1.3:节理几何参数的密度函数有正态分布、对数正态分布、负指数分布、均匀分布四种;
2.1.4:根据求得的随机数,确定用于生成节理的基本几何参数;
2.2:岩体三维裂隙网络模型生成
根据确立的结构面几何概率模型,利用求得的随机数,确定出用于生成节理的基本几何参数,再现出服从这种模型的结构面网络模型,过程如下;
2.2.1:根据结构面数据自动统计结果和求得的随机数,将每组结构面的数据保存到一个文本文件中,用st.dat表示;
2.2.2:st.dat数据内容格式依次为:每条结构面圆盘中心点坐标(x,y,z),圆盘半径D,倾角DA,倾向DD,走向SD,厚度thin,法向方向NX,NY,NZ和节理分组;
2.2.3:为区分不同组别的结构面,对相同组的结构面圆盘赋予相同的颜色,用数列GID表示;
2.2.4:利用Matlab软件,编写程序,读取结构面数据文件st.dat,运行后自动生成岩体三维裂隙网络模型;
2.2.5:得到岩体三维裂隙网络模型;
2.3:二维裂隙网络模型切割
在三维裂隙网络模型的基础上,实现任意角度任意方位的二维裂隙网络模型切割,过程如下;
2.3.1:在三维裂隙网络模型上,结合Matlab软件编程工具,以三维裂隙网络模型中心点为中心,实现任意角度的剖面切割功能;
2.3.2:得到穿过中心点的任意角度的二维裂隙网络模型;
2.3.3:在三维裂隙网络模型上,结合Matlab软件编程工具,在三维裂隙网络模型任意位置上,实现任意角度和方位的剖面切割功能;
2.3.4:得到任意角度和方位的二维裂隙网络模型;
2.3.5:将切割剖面上的数据,保存到st1.dat文件中,此时剖面处于三维坐标系下,文件中数据格式自左向右依次为:节理中心点坐标(x,y,z),节理长度D,倾角DA,倾向DD,走向SD,厚度thin,法向方向NX,NY,NZ;
2.3.6:将三维坐标系转化成二维坐标系,并将二维剖面数据保存到st2.dat文件中,文件中数据格式自左向右依次为:节理中心点坐标(x,y),节理长度D,倾角DA,倾向DD,厚度thin,法向方向NX,NY,NZ;
2.4:裂隙网络模型数据输出
实现三维裂隙网络模型的数据输出和任意二维裂隙网络模型的数据输出,过程如下;
2.4.1:输出st.dat文件数据信息;
2.4.2:输出st1.dat文件数据信息;
2.4.3:输出st2.dat文件数据信息。
本发明具有以下有益效果:
1、本发明基于BQ指标、裂隙网络模型、广义RQD理论和圆锥曲线方程,提供了一种基于BQ和RQDt的巷道开挖不利方位求解方法;
2、本发明实现了岩体BQ指标分级;
3、本发明实现了三维裂隙网络模型生成和二维剖面模型生成;
4、本发明实现了RQDt求解计算和RQDt各向异性图的绘制;
5、本发明实现了RQDt的最佳阈值t求解;
6、本发明实现了RQDt椭圆曲线和各向异性椭球体的求解绘制;
7、本发明实现了巷道开挖的不利方位求解和最优方位求解;
8、本发明方法明确,得到的最佳阈值t能够准确反映出岩体质量,得到的RQDt椭圆曲线和各向异性椭球体能够准确反映出岩体的RQDt各向异性特征,得到的巷道开挖不利方位和最优方位可用于指导巷道开挖施工应用。
附图说明
图1是二维裂隙网络图。
图2是不同阈值t下的RQDt各向异性图。
图3是RQDt随阈值t的拟合曲线。
图4是RQDt椭圆曲线。
具体实施方式
下面参照附图对本发明做进一步说明。
参照图1~图4,一种基于BQ和RQDt的巷道开挖不利方位求解方法,包括以下步骤:
1)基于BQ指标的岩体质量计算,过程如下:
1.1:根据结构面参数计算岩体完整性系数,公式如下:
式中:Jv为岩体体积节理数,单位条/m3;
1.2:Jv计算公式如下:
式中:L1,L2,...,Ln为垂直于结构面测线长度;N1,N2,...,Nn为同组结构面数目;
1.3:根据岩石单轴抗压强度值和岩体完整性系数值,计算BQ值:
BQ=90+3Rc+250Kv (3)
式中:Rc是岩石单轴抗压强度;Kv为岩体完整性系数;
1.4:在应用BQ计算公式过程中,遵循以下条件:
当Rc>90Kv+30时,以Rc=90Kv+30和Kv代入计算BQ值;
当Kv>0.04Rc+0.4时,以Kv=0.04Rc+0.4和Rc代入计算BQ值;
1.5:根据地下水、主要软弱结构面产状和天然应力影响对BQ进行修正,修正公式如下:
[BQ]=BQ-100(K1+K2+K3) (4)
式中:K1为地下水影响修正系数;K2为软件结构面产状影响修正系数;K3为天然应力影响修正系数;
1.6:得到修正后的BQ岩体分级结果,岩体为III级岩体;
2)岩体三维裂隙网络模型生成和剖切,过程如下;
2.1:随机数求解
基于蒙特卡洛方法和节理分布形态,求解随机数,过程如下;
2.1.1:伪随机数的产生,产生随机数的数学方法应满足以下条件:产生的随机数列应均匀分布在(0,1)区间;序列之间应无相关性;随机序列有足够长的重复周期,在计算机上产生的速度快,占有的内存空间小,具有完全可重复性;
2.1.2:给定分布下随机变量数值的产生,Monte Carlo方法是根据确立的结构面几何概率模型,再现服从这种模型的结构面网络模型,该方法的核心是抽样数的随机性,高品质随机数可以得到良好的模拟结果,即在(0,1)区间上生成均匀分布随机变量ri,利用这些均匀随机变量产生服从其他分布的随机数;
2.1.3:节理几何参数的密度函数有正态分布、对数正态分布、负指数分布、均匀分布四种;
2.1.4:根据求得的随机数,确定用于生成节理的基本几何参数;
2.2:岩体三维裂隙网络模型生成
根据确立的结构面几何概率模型,利用求得的随机数,确定出用于生成节理的基本几何参数,再现出服从这种模型的结构面网络模型,过程如下;
2.2.1:根据结构面数据自动统计结果和求得的随机数,将每组结构面的数据保存到一个文本文件中,用st.dat表示;
2.2.2:st.dat数据内容格式依次为:每条结构面圆盘中心点坐标(x,y,z),圆盘半径D,倾角DA,倾向DD,走向SD,厚度thin,法向方向NX,NY,NZ和节理分组;
2.2.3:为区分不同组别的结构面,对相同组的结构面圆盘赋予相同的颜色,用数列GID表示;
2.2.4:利用Matlab软件,编写程序,读取结构面数据文件st.dat,运行后自动生成岩体三维裂隙网络模型;
2.2.5:得到岩体三维裂隙网络模型;
2.3:二维裂隙网络模型切割
在三维裂隙网络模型的基础上,实现任意角度任意方位的二维裂隙网络模型切割,过程如下;
2.3.1:在三维裂隙网络模型上,结合Matlab软件编程工具,以三维裂隙网络模型中心点为中心,实现任意角度的剖面切割功能;
2.3.2:得到穿过中心点的任意角度的二维裂隙网络模型;
2.3.3:在三维裂隙网络模型上,结合Matlab软件编程工具,在三维裂隙网络模型任意位置上,实现任意角度和方位的剖面切割功能;
2.3.4:得到任意角度和方位的二维裂隙网络模型,如图1所示;
2.3.5:将切割剖面上的数据,保存到st1.dat文件中,此时剖面处于三维坐标系下,文件中数据格式自左向右依次为:节理中心点坐标(x,y,z),节理长度D,倾角DA,倾向DD,走向SD,厚度thin,法向方向NX,NY,NZ;
2.3.6:将三维坐标系转化成二维坐标系,并将二维剖面数据保存到st2.dat文件中,文件中数据格式自左向右依次为:节理中心点坐标(x,y),节理长度D,倾角DA,倾向DD,厚度thin,法向方向NX,NY,NZ;
2.4:裂隙网络模型数据输出
实现三维裂隙网络模型的数据输出和任意二维裂隙网络模型的数据输出,过程如下;
2.4.1:输出st.dat文件数据信息;
2.4.2:输出st1.dat文件数据信息;
2.4.3:输出st2.dat文件数据信息;
3)RQDt各向异性图绘制,过程如下:
3.1:RQDt求解计算
基于RQDt理论,确定出二维裂隙网络模型的中心点O,以每隔10°角度做射线,穿过岩体裂隙网络模型,求解出岩体36个方位的RQDt值,过程如下:
3.1.1:RQDt理论公式如下:
式中:xi表示沿某一测线方向的第i个大于给定阈值t的整段岩石或间距长度,RQDt代表对应阈值t的岩石质量指标,即阈值t下的RQDt值;
3.1.2确定二维裂隙网络模型的剖面中心点O、长度a和宽度b,以模型的左下角为坐标原点,水平向右为x轴,垂直向上为y轴,则中心点O坐标为:
边界方程为:
3.1.3:以O为起点,每隔α=10°角绘制一条测线,与裂隙网络模型相交,共绘制36条测线,测线长度L等于O点到裂隙网络模型边界的距离,用L0~L35表示,则测线方程为:
式中:s表示测线,α表示角度;
3.1.4:判断测线与边界的交点,设测线方程与边界方程的交点为(xa,ya),将测线方程与边界方程依次连立,判断测线是否与边界相交,原理如下:
3.1.5:求出测线与裂隙网络模型的边界方程交点(XZ,α,YZ,α);
3.1.6:确定测线所在区间,原理如下:
3.1.7:根据裂隙网络模型中每条节理的起点坐标(Xb,Yb)和终点坐标(Xc,Yc),建立相应的解析方程,定义节理方程如下:
3.1.8:求解第一条测线与各节理方程的交点,循环判断每个交点(Xj,Yj)范围,如果交点符合a<Xj<c且b<Yj<d,则记录该交点,遍历所有节理方程求出所有m个交点;
3.1.9:将记录的m个交点及测线起点坐标和终点坐标,按横坐标或纵坐标从小到大排序;
3.1.10:计算相邻交点的距离,公式如下:
x0=a/2 (14)
y0=b/2 (15)
xm+1=XZ,α (16)
ym+1=YZ,α (17)
3.1.11:输入一个阈值t;
3.1.12:循环比较di和t的大小,设初始lt=0,规则如下:
lt=lt+di,若di>t (18)
3.1.13:求解每条测线对应的RQDt值,用mα表示,以及对应的测线起点与终点的距离lα,公式如下:
3.1.14:循环求出每条测线所对应的mα,获得36个方向上的RQDt值;
3.2:RQDt各向异性图绘制
根据求解出的36个RQDt值,绘制出RQDt的各向异性图,过程如下:
3.2.1:将36个RQDt值,按角度顺序依次排序;
3.2.3:依次连接36组射线的端点,若某一条射线上的RQDt为0,则取圆心;
3.2.4:绘制出RQDt的各向异性图;
3.2.5:输出RQDt各向异性图;
3.3:不同阈值t下RQDt各向异性图绘制
用于绘制不同阈值t下的岩体RQDt各向异性图,过程如下:
3.3.1:输入不同的阈值t,求解出对应的RQDt值;
3.3.2:将不同阈值t下的RQDt值,以角度为变量,以RQDt值为函数,绘制到相同坐标系下;
3.3.3:得到不同阈值t下的RQDt各向异性图;
3.3.4:输出不同阈值t下的RQDt各向异性图,如图2所示;
3.4:输出RQDt的值;
4)基于BQ反演的最佳阈值t求解方法,过程如下:
4.1:基于BQ指标反演RQDt范围
结合BQ分级计算出的岩体质量级别,岩体为III级岩体,查找《岩石质量指标》表,见表1,反演确定出该岩体级别下,RQDt范围值,RQDt范围在50%~75%
表1;
4.2:最佳阈值t求解方法
4.2.1:在三维裂隙网络模型上,过中心点O,以任意角度剖切三个剖面,得到三个二维裂隙网络模型,导出二维裂隙网络模型及数据;
4.2.2:针对每一个二维裂隙网络模型,设置不同的阈值t,求解出不同阈值t下的RQDt值;
4.2.3:导出不同阈值t下的RQDt值,计算出RQDt均值;
4.2.4:以阈值t为横坐标,以RQDt的均值为纵坐标,绘制RQDt随阈值t变化的散点图;
4.2.5:根据散点图,设置拟合方程y=aexp(bx),拟合RQDt随阈值t变化的曲线图,如图3所示;
4.2.6:将反演出的RQDt范围值,带入到拟合出的RQDt随阈值t变化的曲线图中,结合函数方程和曲线图,求解出在该RQDt范围内,阈值t的范围,共得到三组阈值t的范围,见表2
表2;
4.2.7:针对三组阈值t的范围,取其范围的交集,以0.076m~0.124m作为最佳阈值t的范围;
4.2.8:以最佳阈值t范围的中点值作为最佳阈值t值,得到最佳阈值t,最佳阈值t=0.1m;
4.2.9:输出最佳阈值t的范围和最佳阈值t值;
5)空间RQDt求解方法,过程如下
5.1:平面RQDt求解方法
根据最佳阈值t下的36个方位的RQDt值,求解出任意角度的平面RQDt值,过程如下:
5.1.1:求解出最佳阈值t下的36个RQDt值;
5.1.2:以RQDt值沿x方向的分量RQDx作为横坐标,沿y轴方向的分量RQDy作为纵坐标,设定拟合圆锥曲线方程为:
式中:p1,p2,p3为拟合系数;
5.1.3:基于最小二乘法,应用Matlab软件编写拟合程序;
5.1.4:拟合得到RQDt的椭圆曲线,如图4所示;
5.1.5:通过该拟合椭圆曲线,可以得到平面上任意角度的RQDt值;
5.2:空间RQDt求解方法
根据拟合得到的RQDt椭圆曲线,求解出空间内任意方位任意角度的RQDt值,过程如下:
5.2.1:根据拟合得到的RQDt椭圆曲线,求解出椭圆的长轴所在位置;
5.2.2:以椭圆曲线的长轴作为旋转轴,利用Matlab软件编写程序,将椭圆曲线沿长轴旋转360°,得到RQDt的各向异性椭球体;
5.2.3:得到空间内任意方位任意角度的RQDt值;
6)巷道开挖不利方位求解方法,过程如下:
6.1:平面上巷道开挖不利方位求解方法,过程如下:
6.1.1:根据RQDt椭圆曲线,求解出曲线上最大值RQDtmax和最小值RQDtmin的值;
6.1.2:求解出RQDt的各向异性比值,各向异性比计算公式如下:
式中:RQDtmax为椭圆曲线上的RQDt最大值,RQDtmin为椭圆曲线上的RQDt最小值;
6.1.3:各向异性比越大,岩体的不利方位越明显,计算得到的各向异性比为1.93;
6.1.4:判断巷道开挖方向与RQDtmin方向的关系,当巷道轴线与RQDtmin方向一致时,对围岩稳定非常不利;当巷道轴线与RQDtmin方向垂直时,对围岩稳定有利;计算得到RQDtmin的方向为190°,即当巷道开挖方向为190°时,为巷道开挖的不利方位,围岩稳定性处于最不利的状态,当巷道开挖方向为100°时,为巷道开挖的最优方位,围岩稳定性处于最优状态;
6.1.5:得到巷道开挖的不利方位和最优方位,巷道开挖的不利方位为190°,最优方位为100°;
6.2:空间上巷道开挖不利方位求解方法,过程如下:
6.2.1:根据RQDt各向异性椭球体,求解出椭球体上最大值RQDtmax和最小值RQDtmin的值;
6.2.2:求解出RQDt的空间各向异性比值,空间各向异性比计算公式如下:
式中:RQDtmax为椭球体上的RQDt最大值,RQDtmin为椭球体上的RQDt最小值;
6.2.3:各向异性比越大,岩体的不利方位越明显;
6.2.4:判断巷道开挖方向与RQDtmin方向的关系,当巷道轴线与RQDtmin方向一致时,对围岩稳定非常不利;当巷道轴线与RQDtmin方向垂直时,对围岩稳定有利;
6.2.5:得到巷道开挖的空间不利方位和最优方位。
Claims (6)
1.一种基于BQ和RQDt的巷道开挖不利方位求解方法,其特征在于,所述方法包括以下步骤:
(1)基于BQ指标的岩体质量计算,过程如下:
1.1:根据结构面参数计算岩体完整性系数,公式如下:
式中:Jv为岩体体积节理数,单位条/m3;
1.2:Jv计算公式如下:
式中:L1,L2,...,Ln为垂直于结构面测线长度;N1,N2,...,Nn为同组结构面数目;
1.3:根据岩石单轴抗压强度值和岩体完整性系数值,计算BQ值:
BQ=90+3Rc+250Kv
式中:Rc是岩石单轴抗压强度;Kv为岩体完整性系数;
1.4:在应用BQ计算公式过程中,遵循以下条件:
当Rc>90Kv+30时,以Rc=90Kv+30和Kv代入计算BQ值;
当Kv>0.04Rc+0.4时,以Kv=0.04Rc+0.4和Rc代入计算BQ值;
1.5:根据地下水、主要软弱结构面产状和天然应力影响对BQ进行修正,修正公式如下:
[BQ]=BQ-100(K1+K2+K3)
式中:K1为地下水影响修正系数;K2为软件结构面产状影响修正系数;K3为天然应力影响修正系数;
1.6:得到修正后的BQ岩体分级结果;
(2)岩体三维裂隙网络模型生成和剖切;
(3)RQDt各向异性图绘制;
(4)基于BQ反演的最佳阈值t求解方法;
(5)空间RQDt求解方法;
(6)巷道开挖不利方位求解方法。
2.如权利要求1所述的一种基于BQ和RQDt的巷道开挖不利方位求解方法,其特征在于,所述步骤(6)中,巷道开挖不利方位求解方法的过程如下:
6.1:平面上巷道开挖不利方位求解方法,过程如下:
6.1.1:根据RQDt椭圆曲线,求解出曲线上最大值RQDtmax和最小值RQDtmin的值;
6.1.2:求解出RQDt的各向异性比值,各向异性比计算公式如下:
式中:RQDtmax为椭圆曲线上的RQDt最大值,RQDtmin为椭圆曲线上的RQDt最小值;
6.1.3:各向异性比越大,岩体的不利方位越明显;
6.1.4:判断巷道开挖方向与RQDtmin方向的关系,当巷道轴线与RQDtmin方向一致时,对围岩稳定非常不利;当巷道轴线与RQDtmin方向垂直时,对围岩稳定有利;
6.1.5:得到巷道开挖的不利方位和最优方位;
6.2:空间上巷道开挖不利方位求解方法,过程如下:
6.2.1:根据RQDt各向异性椭球体,求解出椭球体上最大值RQDtmax和最小值RQDtmin的值;
6.2.2:求解出RQDt的空间各向异性比值,空间各向异性比计算公式如下:
式中:RQDtmax为椭球体上的RQDt最大值,RQDtmin为椭球体上的RQDt最小值;
6.2.3:各向异性比越大,岩体的不利方位越明显;
6.2.4:判断巷道开挖方向与RQDtmin方向的关系,当巷道轴线与RQDtmin方向一致时,对围岩稳定非常不利;当巷道轴线与RQDtmin方向垂直时,对围岩稳定有利;
6.2.5:得到巷道开挖的空间不利方位和最优方位。
3.如权利要求1所述的一种基于BQ和RQDt的巷道开挖不利方位求解方法,其特征在于,所述步骤(5)中,空间RQDt求解方法的过程如下:
5.1:平面RQDt求解方法
根据最佳阈值t下的36个方位的RQDt值,求解出任意角度的平面RQDt值,过程如下:
5.1.1:求解出最佳阈值t下的36个RQDt值;
5.1.2:以RQDt值沿x方向的分量RQDx作为横坐标,沿y轴方向的分量RQDy作为纵坐标,设定拟合圆锥曲线方程为:
式中:p1,p2,p3为拟合系数;
5.1.3:基于最小二乘法,应用Matlab软件编写拟合程序;
5.1.4:拟合得到RQDt的椭圆曲线;
5.1.5:通过该拟合椭圆曲线,可以得到平面上任意角度的RQDt值;
5.2:空间RQDt求解方法
根据拟合得到的RQDt椭圆曲线,求解出空间内任意方位任意角度的RQDt值,过程如下:
5.2.1:根据拟合得到的RQDt椭圆曲线,求解出椭圆的长轴所在位置;
5.2.2:以椭圆曲线的长轴作为旋转轴,利用Matlab软件编写程序,将椭圆曲线沿长轴旋转360°,得到RQDt的各向异性椭球体;
5.2.3:得到空间内任意方位任意角度的RQDt值。
4.如权利要求1所述的一种基于BQ和RQDt的巷道开挖不利方位求解方法,其特征在于,所述步骤(4)中,基于BQ反演的最佳阈值t求解方法的过程如下:
4.1:基于BQ指标反演RQDt范围
结合BQ分级计算出的岩体质量级别,查找《岩石质量指标》表,反演确定出该岩体级别下,RQDt范围值;
4.2:最佳阈值t求解方法
4.2.1:在三维裂隙网络模型上,过中心点O,以任意角度剖切三个剖面,得到三个二维裂隙网络模型,导出二维裂隙网络模型及数据;
4.2.2:针对每一个二维裂隙网络模型,设置不同的阈值t,求解出不同阈值t下的RQDt值;
4.2.3:导出不同阈值t下的RQDt值,计算出RQDt均值;
4.2.4:以阈值t为横坐标,以RQDt的均值为纵坐标,绘制RQDt随阈值t变化的散点图;
4.2.5:根据散点图,设置拟合方程,拟合RQDt随阈值t变化的曲线图;
4.2.6:将反演出的RQDt范围值,带入到拟合出的RQDt随阈值t变化的曲线图中,结合函数方程和曲线图,求解出在该RQDt范围内阈值t的范围,共得到三组阈值t的范围;
4.2.7:针对三组阈值t的范围,取其范围的交集,作为最佳阈值t的范围;
4.2.8:以最佳阈值t范围的中点值作为最佳阈值t值,得到最佳阈值t;
4.2.9:输出最佳阈值t的范围和最佳阈值t值。
5.如权利要求1所述的一种基于BQ和RQDt的巷道开挖不利方位求解方法,其特征在于,所述步骤(3)中,RQDt各向异性图绘制的过程如下:
3.1:RQDt求解计算
基于RQDt理论,确定出二维裂隙网络模型的中心点O,以每隔10°角度做射线,穿过岩体裂隙网络模型,求解出岩体36个方位的RQDt值,过程如下:
3.1.1:RQDt理论公式如下:
式中:xi表示沿某一测线方向的第i个大于给定阈值t的整段岩石或间距长度,RQDt代表对应阈值t的岩石质量指标,即阈值t下的RQDt值;
3.1.2确定二维裂隙网络模型的剖面中心点O、长度a和宽度b,以模型的左下角为坐标原点,水平向右为x轴,垂直向上为y轴,则中心点O坐标为:
边界方程为:
3.1.3:以O为起点,每隔α=10°角绘制一条测线,与裂隙网络模型相交,共绘制36条测线,测线长度L等于O点到裂隙网络模型边界的距离,用L0~L35表示,则测线方程为:
式中:s表示测线,α表示角度;
3.1.4:判断测线与边界的交点,设测线方程与边界方程的交点为(xa,ya),将测线方程与边界方程依次连立,判断测线是否与边界相交,原理如下:
3.1.5:求出测线与裂隙网络模型的边界方程交点(XZ,α,YZ,α);
3.1.6:确定测线所在区间,原理如下:
3.1.7:根据裂隙网络模型中每条节理的起点坐标(Xb,Yb)和终点坐标(Xc,Yc),建立相应的解析方程,定义节理方程如下:
3.1.8:求解第一条测线与各节理方程的交点,循环判断每个交点(Xj,Yj)范围,如果交点符合a<Xj<c且b<Yj<d,则记录该交点,遍历所有节理方程求出所有m个交点;
3.1.9:将记录的m个交点及测线起点坐标和终点坐标,按横坐标或纵坐标从小到大排序;
3.1.10:计算相邻交点的距离,公式如下:
x0=a/2
y0=b/2
xm+1=XZ,α
ym+1=YZ,α
3.1.11:输入一个阈值t;
3.1.12:循环比较di和t的大小,设初始lt=0,规则如下:
lt=lt+di,若di>t
3.1.13:求解每条测线对应的RQDt值,用mα表示,以及对应的测线起点与终点的距离lα,公式如下:
3.1.14:循环求出每条测线所对应的mα,获得36个方向上的RQDt值;
3.2:RQDt各向异性图绘制
根据求解出的36个RQDt值,绘制出RQDt的各向异性图,过程如下:
3.2.1:将36个RQDt值,按角度顺序依次排序;
3.2.3:依次连接36组射线的端点,若某一条射线上的RQDt为0,则取圆心;
3.2.4:绘制出RQDt的各向异性图;
3.2.5:输出RQDt各向异性图;
3.3:不同阈值t下RQDt各向异性图绘制
用于绘制不同阈值t下的岩体RQDt各向异性图,过程如下:
3.3.1:输入不同的阈值t,求解出对应的RQDt值;
3.3.2:将不同阈值t下的RQDt值,以角度为变量,以RQDt值为函数,绘制到相同坐标系下;
3.3.3:得到不同阈值t下的RQDt各向异性图;
3.3.4:输出不同阈值t下的RQDt各向异性图;
3.4:输出RQDt的值。
6.如权利要求1所述的一种基于BQ和RQDt的巷道开挖不利方位求解方法,其特征在于,所述步骤(2)中,岩体三维裂隙网络模型生成和剖切的过程如下:
2.1:随机数求解
基于蒙特卡洛方法和节理分布形态,求解随机数,过程如下;
2.1.1:伪随机数的产生,产生随机数的数学方法应满足以下条件:产生的随机数列应均匀分布在(0,1)区间;序列之间应无相关性;随机序列有足够长的重复周期,在计算机上产生的速度快,占有的内存空间小,具有完全可重复性;
2.1.2:给定分布下随机变量数值的产生,Monte Carlo方法是根据确立的结构面几何概率模型,再现服从这种模型的结构面网络模型,该方法的核心是抽样数的随机性,高品质随机数可以得到良好的模拟结果,即在(0,1)区间上生成均匀分布随机变量ri,利用这些均匀随机变量产生服从其他分布的随机数;
2.1.3:节理几何参数的密度函数有正态分布、对数正态分布、负指数分布、均匀分布四种;
2.1.4:根据求得的随机数,确定用于生成节理的基本几何参数;
2.2:岩体三维裂隙网络模型生成
根据确立的结构面几何概率模型,利用求得的随机数,确定出用于生成节理的基本几何参数,再现出服从这种模型的结构面网络模型,过程如下;
2.2.1:根据结构面数据自动统计结果和求得的随机数,将每组结构面的数据保存到一个文本文件中,用st.dat表示;
2.2.2:st.dat数据内容格式依次为:每条结构面圆盘中心点坐标(x,y,z),圆盘半径D,倾角DA,倾向DD,走向SD,厚度thin,法向方向NX,NY,NZ和节理分组;
2.2.3:为区分不同组别的结构面,对相同组的结构面圆盘赋予相同的颜色,用数列GID表示;
2.2.4:利用Matlab软件,编写程序,读取结构面数据文件st.dat,运行后自动生成岩体三维裂隙网络模型;
2.2.5:得到岩体三维裂隙网络模型;
2.3:二维裂隙网络模型切割
在三维裂隙网络模型的基础上,实现任意角度任意方位的二维裂隙网络模型切割,过程如下;
2.3.1:在三维裂隙网络模型上,结合Matlab软件编程工具,以三维裂隙网络模型中心点为中心,实现任意角度的剖面切割功能;
2.3.2:得到穿过中心点的任意角度的二维裂隙网络模型;
2.3.3:在三维裂隙网络模型上,结合Matlab软件编程工具,在三维裂隙网络模型任意位置上,实现任意角度和方位的剖面切割功能;
2.3.4:得到任意角度和方位的二维裂隙网络模型;
2.3.5:将切割剖面上的数据,保存到st1.dat文件中,此时剖面处于三维坐标系下,文件中数据格式自左向右依次为:节理中心点坐标(x,y,z),节理长度D,倾角DA,倾向DD,走向SD,厚度thin,法向方向NX,NY,NZ;
2.3.6:将三维坐标系转化成二维坐标系,并将二维剖面数据保存到st2.dat文件中,文件中数据格式自左向右依次为:节理中心点坐标(x,y),节理长度D,倾角DA,倾向DD,厚度thin,法向方向NX,NY,NZ;
2.4:裂隙网络模型数据输出
实现三维裂隙网络模型的数据输出和任意二维裂隙网络模型的数据输出,过程如下;
2.4.1:输出st.dat文件数据信息;
2.4.2:输出st1.dat文件数据信息;
2.4.3:输出st2.dat文件数据信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010892992.9A CN112132409A (zh) | 2020-08-25 | 2020-08-25 | 一种基于BQ和RQDt的巷道开挖不利方位求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010892992.9A CN112132409A (zh) | 2020-08-25 | 2020-08-25 | 一种基于BQ和RQDt的巷道开挖不利方位求解方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112132409A true CN112132409A (zh) | 2020-12-25 |
Family
ID=73848718
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010892992.9A Withdrawn CN112132409A (zh) | 2020-08-25 | 2020-08-25 | 一种基于BQ和RQDt的巷道开挖不利方位求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112132409A (zh) |
-
2020
- 2020-08-25 CN CN202010892992.9A patent/CN112132409A/zh not_active Withdrawn
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112150001A (zh) | 基于摄影测量、BQ、改进Mathews稳定图的围岩稳定性评价方法 | |
CN112132407A (zh) | 一种基于BQ反演最佳阈值t的空间RQDt求解方法 | |
CN112200423A (zh) | 一种基于bq、数值模拟的围岩稳定性动态评价方法 | |
CN112200419A (zh) | 基于激光扫描、BQ、改进Mathews稳定图的围岩稳定性评价方法 | |
CN112200426A (zh) | 基于激光扫描、bq、数值模拟的围岩稳定性动态评价方法 | |
CN112200417A (zh) | 基于摄影测量、BQ、RQDt、地应力的改进Mathews稳定图评价方法 | |
CN112132409A (zh) | 一种基于BQ和RQDt的巷道开挖不利方位求解方法 | |
CN112150002A (zh) | 基于激光扫描、BQ、RQDt、地应力的改进Mathews稳定图评价方法 | |
CN112149996A (zh) | 基于激光扫描、BQ、RQDt各向异性的改进Mathews稳定图评价方法 | |
CN112464514A (zh) | 一种基于摄影测量、RQD和RQDt的巷道开挖不利方位求解方法 | |
CN112132411A (zh) | 基于激光扫描、BQ、RQDt各向异性的Q各向异性求解方法 | |
CN112115403A (zh) | 一种基于BQ、最佳阈值t的RQDt各向异性求解方法 | |
CN112200431A (zh) | 基于激光扫描、bq、数值模拟的空区稳定性动态评价方法 | |
CN112132408A (zh) | 一种基于激光扫描和BQ反演最佳阈值t的空间RQDt求解方法 | |
CN112464509A (zh) | 一种基于RQD和RQDt的巷道开挖不利方位求解方法 | |
CN112200418A (zh) | 基于摄影测量、bq、数值模拟的围岩稳定性动态评价方法 | |
CN112113838A (zh) | 一种基于BQ、RQDt各向异性的Q各向异性求解方法 | |
CN112464511A (zh) | 一种基于RQD、最佳阈值t的RQDt各向异性求解方法 | |
CN112464516A (zh) | 一种基于激光扫描和RQD反演最佳阈值t的空间RQDt求解方法 | |
CN112115402A (zh) | 一种基于BQ反演的RQDt最佳阈值t求解方法 | |
CN112200422A (zh) | 一种基于BQ、RQDt各向异性的Mathews稳定图法的改进方法 | |
CN112200425A (zh) | 一种基于BQ、RQDt各向异性的改进Mathews稳定图评价方法 | |
CN112200430A (zh) | 一种基于BQ、RQDt、地应力的Mathews稳定图法的改进方法 | |
CN112132410A (zh) | 一种基于摄影测量和BQ反演最佳阈值t的空间RQDt求解方法 | |
CN112464503A (zh) | 一种基于RQD反演最佳阈值t的空间RQDt求解方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20201225 |