CN108344385B - 一种节理三维形貌表征方法 - Google Patents
一种节理三维形貌表征方法 Download PDFInfo
- Publication number
- CN108344385B CN108344385B CN201810113931.0A CN201810113931A CN108344385B CN 108344385 B CN108344385 B CN 108344385B CN 201810113931 A CN201810113931 A CN 201810113931A CN 108344385 B CN108344385 B CN 108344385B
- Authority
- CN
- China
- Prior art keywords
- joint
- dimensional
- plane
- stress direction
- triangle
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/30—Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/28—Measuring arrangements characterised by the use of optical techniques for measuring areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
Abstract
本发明公开了节理三维形貌表征方法,该方法是以表观表面积率作为三维形貌表征的唯一核心形貌表征参数,充分考虑了受力、渗流方向的影响,针对同一节理面不同的方向得到不同的粗糙度大小,简单明了,便于工程应用。
Description
技术领域
本发明涉及岩石节理形貌表征方法,尤其是涉及一种能够用于岩石节理抗剪强度及渗透特性评估的三维节理面粗糙度表征方法。
技术背景
近年来,随着水利水电、核废料储藏、石油矿山开采等大型工程的建设,工程岩体所处的地质环境越来越复杂。岩石节理抗剪强度及其渗透特性作为对岩体工程稳定性有着重要影响的因素越来越受到研究人员的重视,岩体中地应力与地下水的相互作用往往直接决定着工程结构的稳定与安全。
由于节理形貌对于节理抗剪强度、渗透特性的重要性,该方面的研究随着形貌表征方法的进步而不断发展。早年的研究成果大多是基于节理剖面的二维方法,主要分为三类:经验参数方法、统计参数方法及分形参数方法。其中,统计参数物理意义明确,运用统计参数表征节理粗糙度有利于节理本构的提出;经验参数更利于工程应用,节理粗糙度表征参数JRC是目前应用最广泛的节理粗糙度表征方法,并被国际岩石力学学会推荐使用。随着时代的发展,三维激光扫描仪的出现为研究者开展三维节理表面粗糙度表征提供了有力的测量工具,国内外学者通过对二维统计参数进行拓展展开了三维统计参数的研究。
然而,1,基于节理剖面的二维方法从二维参数出发描述节理粗糙度,低估了节理粗糙程度,与节理真实粗糙度有较大偏差,因此,会低估节理抗剪强度及渗透特性。2,目前的三维形貌表征方法大多是基于二维统计方法进行的扩展,大多需要多个统计参数对节理面粗糙程度进行表征,不利于工程应用。3目前的节理形貌表征方法无法体现节理面的各向异性,因此,无法应用到考虑剪切方向及渗流方向的工程应用中。
发明内容
本发明所要解决的技术问题是提供一种岩石节理三维形貌表征方法,提供一种基于三维激光扫描技术的三维形貌表征参数,即表观表面积率。
本发明还有一目的是解决现有形貌表征方法中表征参数多且片面,难以进行工程应用的技术问题。提供一种基于一个核心形貌表征参数即表观表面积率的三维形貌表征方法,表征方法简单明了,表征参数易于获得,有利于工程应用。
本发明再有一目的是解决现有形貌表征方法无法反应节理各向异性的技术问题;提供一种能够反应节理面各向异性,能够考虑工程中节理受力方向、节理中渗流方向的节理三维形貌表征方法。
本发明解决其技术问题所采用的技术方案是:
采用三维激光扫描仪对岩石节理进行三维点云重构,三维激光扫描仪自定义一个“大地坐标系”,即包含X、Y方向以及一个高程基准面H(该基准面初始高程h0小于或等于节理点最小高程值hmin)。从而生成的点云数据,具有在该“大地坐标系”下的三维坐标(X,Y,H)。
节理三维点云中,选择高程最大值点A1,其坐标值为(x,y,hmax)。将节理面受力(渗流)方向线定义为K1,受力(渗流)方向面定义为M1。
由高程最大值点A1向受力方向面M1作一条法线S1,从而得到无数个与受力方向面垂直的基于S1的法面MSi;每个法面MSi与方向面的交线为MJi,最终得到一条与节理受力方向线K1空间垂直的MJ1,从而得到通过MJ1且与K1垂直的法面MS1。
将法面MS1与受力方向面M1作为二维平面的X,Y坐标系的Xm,Ym,所有点云在该坐标系下的坐标为(Xmi,Ymi),定义受力方向为正值,受力方向反方向为负值。
将Xmi中大于零的点云选取出,重构TIN三角形,得到总面积T1;同理,将Xmi中小于零的点云选取出,重构TIN三角形,得到总面积T2。TIN全称不规则三角网。在TIN中,满足最佳三角形的条件为:尽可能的保证三角形的三个角都是锐角,三角形的三条边近似相等,最小角最大化。TIN是基于矢量的数字地理数据的一种形式,通过将一系列折点(点)组成三角形来构建,此时TIN非常有用,因为它们允许计算平面面积、表面积和体积。
定义表观表面积率SA=T1/T2。表观表面积率越大,表明节理粗糙度越高,节理的抗剪能力越强,节理的渗流能力越差。
由上述技术方案可知,本发明采用一个核心参数表征节理三维形貌,充分考虑了受力、渗流方向的影响,针对同一节理面不同的方向得到不同的粗糙度大小,简单明了,便于工程应用。
附图说明
图1是岩石节理扫描图。
图2是岩石节理面扫描点集图。
图3是表观表面积率获得示意图。
图中标记为:1-传感头,2-节理试样,3-位移台,4-控制器,5-高程最大值点,6-节理面受力方向,7-受力方向面,8-法线,9-垂直钱,10-法面。
具体实施方式
下面结合实施例对本发明作进一步说明,但不作为对本发明的任何限制。
实施例1:
对岩石节理进行点云采集。图1是一台三维非接触式高速、高精度激光扫描系统,该仪器由传感头1、控制器4、数据采集器等组成。岩体试样2放在可沿x、y方向移动的位移台3上,岩石结构面表面的起伏形状通过激光变位计转换成光信号,光信号再由CCD转换成电信号,结构面表面变化的起伏形貌就经由变化的电信号处理,最后由数据采集软件LK-Navigator以点集的形式存储。位移台带动岩块试样移动,移动的方向及路线配合扫描仪,即可获取结构面表面形貌上点的坐标。
利用三维激光扫描仪对岩石节理获取三维点云后,三维激光扫描仪自定义一个“大地坐标系”,即包含X、Y方向以及一个高程基准面H(该基准面初始高程h0小于或等于节理点最小高程值hmin)。从而生成的点云数据具有在该“大地坐标系”下的三维坐标(X,Y,H)。图2为一岩石节理扫描得到的表面形貌点集。
节理三维点云中,选择高程最大值点5定义为A1,将节理面受力方向线6定义为K1,受力方向面7定义为M1。
由高程最大值点5向受力方向面7作一条法线8,从而得到无数个与受力方向面7垂直的基于法线8的法面10;每个法面10与受力方向面7的交线为MJi,最终得到一条与节理受力方向线6空间垂直的垂直线9,从而得到通过垂直线9且与节理面受力方向线6垂直的法面10。
将法面10与受力方向面7作为二维平面的X,Y坐标系的Xm,Ym,所有点云在该坐标系下的坐标为(Xmi,Ymi),定义受力方向为正值,受力方向反方向为负值。
将Xmi中大于零的点云选取出,重构TIN三角形,得到总面积T1;同理,将Xmi中小于零的点云选取出,重构TIN三角形,得到总面积T2。
得到表观表面积率SA=T1/T2。
在工程现场,采取获得不同天然节理试样后进行扫描,得到不同节理的表观表面积数值,即可对岩石节理面的粗糙程度进行评估,进一步可以评估该岩石节理的抗剪强度、渗透特性。
以上只是本发明的具体应用范例,本发明还有其他的实施方式,凡采用等同替换或等效变换形成的技术方案,均落在本发明所要求的保护范围之内。
Claims (1)
1.一种节理三维形貌表征方法,其特征在于:以表观表面积率作为三维形貌表征的唯一核心形貌表征参数;该方法包括如下步骤:
步骤1,采用三维激光扫描仪对岩石节理进行三维点云重构;
步骤2,节理三维点云中,选择高程最大值点A1,将节理面受力方向线定义为K1,受力方向面定义为M1;
步骤3,由高程最大值点A1向受力方向面M1作一条法线S1,从而得到无数个与受力方向面垂直的基于S1的法面MSi;每个法面MSi与方向面的交线为MJi,最终得到一条与节理受力方向线K1空间垂直的MJ1,从而得到通过MJ1且与K1垂直的法面MS1;
步骤4,将法面MS1与受力方向面M1作为二维平面的X,Y坐标系的Xm,Ym,所有点云在该坐标系下的坐标为(Xmi,Ymi),定义受力方向为正值,受力方向反方向为负值;
步骤5,将Xmi中大于零的点云选取出,重构TIN三角形,得到总面积T1;同理,将Xmi中小于零的点云选取出,重构TIN三角形,得到总面积T2;
步骤6,在TIN中,满足最佳三角形的条件为:尽可能的保证三角形的三个角都是锐角,三角形的三条边近似相等,最小角最大化;
步骤7,定义表观表面积率SA=T1/T2;表观表面积率越大,表明节理粗糙度越高,节理的抗剪能力越强,节理的渗流能力越差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810113931.0A CN108344385B (zh) | 2018-02-05 | 2018-02-05 | 一种节理三维形貌表征方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810113931.0A CN108344385B (zh) | 2018-02-05 | 2018-02-05 | 一种节理三维形貌表征方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108344385A CN108344385A (zh) | 2018-07-31 |
CN108344385B true CN108344385B (zh) | 2020-10-23 |
Family
ID=62959647
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810113931.0A Active CN108344385B (zh) | 2018-02-05 | 2018-02-05 | 一种节理三维形貌表征方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108344385B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109948205B (zh) * | 2019-03-06 | 2023-05-26 | 东南大学 | 一种基于三维形态描述的节理面粗糙度计算方法 |
CN111649672B (zh) * | 2020-06-12 | 2021-06-25 | 东北大学 | 基于钻孔摄像技术的节理岩体体积节理数确定方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10148512A (ja) * | 1996-11-18 | 1998-06-02 | Nippon Telegr & Teleph Corp <Ntt> | 形状計測装置 |
CN103759677B (zh) * | 2014-01-27 | 2016-04-20 | 东北大学 | 基于三角面积比法度量岩体结构面三维粗糙度的方法 |
CN103793938A (zh) * | 2014-03-05 | 2014-05-14 | 辽宁工程技术大学 | 一种山地表面构造方法 |
CN104008570A (zh) * | 2014-05-09 | 2014-08-27 | 中国地质大学(北京) | 一种矿山的双三维建模方法 |
-
2018
- 2018-02-05 CN CN201810113931.0A patent/CN108344385B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108344385A (zh) | 2018-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111079217B (zh) | 一种基于bim的岩土工程综合勘察信息解译方法及系统 | |
CN104075944B (zh) | 基于普通土压力盒和菱形十二面体的三维土压力测试装置及装配计算方法 | |
CN103806478A (zh) | 一种地下连续墙变形的测量方法 | |
CN108344385B (zh) | 一种节理三维形貌表征方法 | |
CN110702881B (zh) | 岩土材料参数变异性结果的预测方法及其应用 | |
CN104613904A (zh) | 一种岩体结构面粗糙度系数的综合确定方法 | |
CN111595403B (zh) | 一种基于点云测量技术的工程土方计量方法 | |
CN105445789A (zh) | 基于多次反射折射波约束的三维菲涅尔体旅行时层析成像方法 | |
CN108896004A (zh) | 一种裂缝面粗糙度各向异性表征方法 | |
CN209960156U (zh) | 一种三维激光扫描用辅助支架 | |
CN113946980B (zh) | 一种大型露天矿开挖工程台阶土石方量计算方法及系统 | |
CN113536414B (zh) | 基于三维建模的岩质边坡稳定性分析方法、系统及介质 | |
CN102778689B (zh) | 一种宽弯线地震资料地下反射线建立方法 | |
CN110197013A (zh) | 基于Morphing的河床基岩面建模方法 | |
CN112950702A (zh) | 一种基于三维点云的煤堆体积计算方法 | |
CN117313343A (zh) | 一种基于三维激光扫描仪爆破区域数字化模型获取方法 | |
Sánchez Aparicio et al. | A comparative study between WMMS and TLS for the stability analysis of the San Pedro Church Barrel vault by means of the finite element method | |
Pan et al. | Assessment method of slope excavation quality based on point cloud data | |
CN107941194A (zh) | 一种获取和计算工程岩体代表性rqd值的方法 | |
CN106846481B (zh) | 一种地质剖面图的生成方法 | |
Hu et al. | Slope excavation quality assessment and excavated volume calculation in hydraulic projects based on laser scanning technology | |
Li et al. | Application of ANSYS 3D FEM in studies of surface deformation caused by pipe jacking | |
CN111596356B (zh) | 一种昔格达组地层岩质边坡的地震惯性力计算方法 | |
Jo et al. | Displacement analysis of five-story stone pagoda in Geumgolsan Mountain, Jindo, using terrestrial laser scanning | |
CN1959342A (zh) | 一种机械式微地形探测仪 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: 550002 No.27 Baoshan South Road, Nanming District, Guiyang City, Guizhou Province Patentee after: Guizhou water conservancy and Hydropower Survey, design and Research Institute Co.,Ltd. Address before: 550002 No.27 Baoshan South Road, Nanming District, Guiyang City, Guizhou Province Patentee before: GUIZHOU SURVEY & DESIGN Research Institute FOR WATER RESOURCES AND HYDROPOWER |