CN108344385A - 一种节理三维形貌表征方法 - Google Patents
一种节理三维形貌表征方法 Download PDFInfo
- Publication number
- CN108344385A CN108344385A CN201810113931.0A CN201810113931A CN108344385A CN 108344385 A CN108344385 A CN 108344385A CN 201810113931 A CN201810113931 A CN 201810113931A CN 108344385 A CN108344385 A CN 108344385A
- Authority
- CN
- China
- Prior art keywords
- joint
- dimensional
- impact direction
- face
- point cloud
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000012512 characterization method Methods 0.000 claims abstract description 13
- 239000011435 rock Substances 0.000 claims description 27
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 claims description 9
- 230000001154 acute effect Effects 0.000 claims description 2
- 238000005325 percolation Methods 0.000 claims description 2
- 230000035515 penetration Effects 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
Landscapes
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Environmental & Geological Engineering (AREA)
- Dispersion Chemistry (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Length Measuring Devices By Optical Means (AREA)
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 (2)
1.一种节理三维形貌表征方法,其特征在于:以表观表面积率作为三维形貌表征的唯一核心形貌表征参数。
2.根据权利要求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 true CN108344385A (zh) | 2018-07-31 |
CN108344385B 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) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109948205A (zh) * | 2019-03-06 | 2019-06-28 | 东南大学 | 一种基于三维形态描述的节理面粗糙度计算方法 |
CN111649672A (zh) * | 2020-06-12 | 2020-09-11 | 东北大学 | 基于钻孔摄像技术的节理岩体体积节理数确定方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10148512A (ja) * | 1996-11-18 | 1998-06-02 | Nippon Telegr & Teleph Corp <Ntt> | 形状計測装置 |
CN103759677A (zh) * | 2014-01-27 | 2014-04-30 | 东北大学 | 基于三角面积比法度量岩体结构面三维粗糙度的方法 |
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
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10148512A (ja) * | 1996-11-18 | 1998-06-02 | Nippon Telegr & Teleph Corp <Ntt> | 形状計測装置 |
CN103759677A (zh) * | 2014-01-27 | 2014-04-30 | 东北大学 | 基于三角面积比法度量岩体结构面三维粗糙度的方法 |
CN103793938A (zh) * | 2014-03-05 | 2014-05-14 | 辽宁工程技术大学 | 一种山地表面构造方法 |
CN104008570A (zh) * | 2014-05-09 | 2014-08-27 | 中国地质大学(北京) | 一种矿山的双三维建模方法 |
Non-Patent Citations (1)
Title |
---|
王小江: "石结构面力学及水力特性实验研究", 《中国博士学位论文全文数据库》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109948205A (zh) * | 2019-03-06 | 2019-06-28 | 东南大学 | 一种基于三维形态描述的节理面粗糙度计算方法 |
CN109948205B (zh) * | 2019-03-06 | 2023-05-26 | 东南大学 | 一种基于三维形态描述的节理面粗糙度计算方法 |
CN111649672A (zh) * | 2020-06-12 | 2020-09-11 | 东北大学 | 基于钻孔摄像技术的节理岩体体积节理数确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108344385B (zh) | 2020-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
NL2030786B1 (en) | Method for identifying potential landslide hazard of reservoir bank based on rock mass degradation feature | |
CN106644206B (zh) | 一种用于浅埋隧道的围岩压力计算方法 | |
Teza et al. | Geometric characterization of a cylinder-shaped structure from laser scanner data: Development of an analysis tool and its use on a leaning bell tower | |
CN105526908A (zh) | 一种三维激光扫描和gps相结合的边坡监测方法 | |
CN107893437A (zh) | 基于远程无线传输技术的大型沉井基础施工实时监控系统 | |
CN108344385A (zh) | 一种节理三维形貌表征方法 | |
CN105928957A (zh) | 一种裂缝性碳酸盐岩三维数字岩心的构建方法 | |
Zhang et al. | Determination of the optimal threshold and length measurements for RQD calculations | |
CN111595403A (zh) | 一种基于点云测量技术的工程土方计量方法 | |
CN104573200B (zh) | 基于bim的内河码头桩基碰撞调整方法 | |
Gousie et al. | Constructing a DEM from grid-based data by computing intermediate contours | |
CN113536414B (zh) | 基于三维建模的岩质边坡稳定性分析方法、系统及介质 | |
KR101271402B1 (ko) | 침식기반 프랙탈기법을 이용한 하천지형 보간 방법 및 컴퓨터로 판독 가능한 기록매체 | |
CN109859301A (zh) | 一种岩石结构面粗糙度系数精细化表征方法 | |
CN109493425A (zh) | 一种矿山三维采空区实体模型的建立方法 | |
CN109712241A (zh) | 一种包含采空区的三维矿山实体模型的建立方法 | |
CN117313343A (zh) | 一种基于三维激光扫描仪爆破区域数字化模型获取方法 | |
CN106846481B (zh) | 一种地质剖面图的生成方法 | |
CN106920275A (zh) | 一种复杂属性边界三维矢量迭代方法及应用系统 | |
CN207894876U (zh) | 隧道围岩爆破损伤深度声波探测中的声测孔布局 | |
CN111596356B (zh) | 一种昔格达组地层岩质边坡的地震惯性力计算方法 | |
Roberson et al. | Application of borehole optical televiewing to investigating the 3-D structure of glaciers: implications for the formation of longitudinal debris ridges, midre Lovénbreen, Svalbard | |
Zeng et al. | Simple building reconstruction from LIDAR point cloud | |
CN112147619A (zh) | 基于声呐点云数据的桩体间距离迭代确定方法 | |
CN109858076A (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 |
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 |
|
CP01 | Change in the name or title of a patent holder |