CN106228559A - 基于分解‑重构建立裂隙岩体三维数值模型的方法 - Google Patents
基于分解‑重构建立裂隙岩体三维数值模型的方法 Download PDFInfo
- Publication number
- CN106228559A CN106228559A CN201610600268.8A CN201610600268A CN106228559A CN 106228559 A CN106228559 A CN 106228559A CN 201610600268 A CN201610600268 A CN 201610600268A CN 106228559 A CN106228559 A CN 106228559A
- Authority
- CN
- China
- Prior art keywords
- rock
- reference point
- numerical model
- crack
- block
- 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
Classifications
-
- 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
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/08—Indexing scheme for image data processing or generation, in general involving all processing steps from image acquisition to 3D model generation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30184—Infrastructure
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
- G06T2219/20—Indexing scheme for editing of 3D models
- G06T2219/2008—Assembling, disassembling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Remote Sensing (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Processing Or Creating Images (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种基于分解‑重构建立裂隙岩体三维数值模型的方法,包括如下步骤:1)在岩石表面贴上用于坐标识别的参考点;2)利用三维扫描仪对岩体表面进行扫描,获得裂隙岩体的三维表面形态的测量数据,利用逆向工程软件完成对扫描数据的拼合,获得岩体的三维数值模型,记录参考点的坐标;3)沿着岩体内部的天然贯通裂隙将块体分解,分解后获得含有内部裂隙面的各块体,保留块体上的参考点;4)利用三维扫描仪对分解的块体进行逐个扫描,同步骤2)获得三维表面数值模型,记录扫描的参考点坐标;5)利用分解之前参考点的坐标对获得的三维数值模型进行坐标转化拼接,即获得重构的裂隙岩体三维数值模型。本发明有效模拟岩体内部真实的裂隙几何形态。
Description
技术领域
本发明涉及一种基于对裂隙岩体的分解-重构建立裂隙岩体三维数值模型的方法,适用于内部存在贯通裂隙的天然岩体。
背景技术
天然岩体的形成需要经历一个复杂而又漫长的过程,由于不同的地质成因,岩体在经历多次构造运动的改造作用,发生了构造变形与破裂,破坏了岩体的连续性和完整性。以上作用会在其内部形成众多节理、裂隙等不连续构造,其裂隙岩体内部富含各种缺陷,包括微裂纹、孔隙以及节理裂隙等宏观非连续面,在地下内部复杂的水环境下,这些岩体裂隙严重影响着岩体的力学特性和渗透特性。
随着经济的发展,边坡加固、地铁、隧道、水利水电等岩土工程正在快速的建设中,而岩体内部的裂隙对工程岩体的稳定性有着十分重要的影响,国内外很多岩土工程的事故都是由于贯通裂隙造成的。针对此问题,国内外众多学者对岩体的力学与渗流问题进行了广泛的讨论和理论研究。一般基于裂隙长度、开度、频率分布等参数,描述岩体内裂隙的几何形态,这些参数都以其各自的方式影响着裂隙的力学和渗流特性。
目前基于裂隙三维数值模型对裂隙岩体的研究较多,其优点在于1)可以模拟岩体内部裂隙的发育情况,较为真实的反映出天然岩体的裂隙;2)数值模型在电脑操作,便于对模型的数值模拟分析;3)可以改变裂隙的各项几何参数,分析各项参数的影响。但目前裂隙岩体的三维数值模型多基于对内部裂隙的抽象模拟,并不能完全反映岩体内部裂隙的真实发育情况。
发明内容
为了克服现有的裂隙三维数值模型不能完全模拟岩体内部真实的裂隙几何形态的不足,本发明提供了一种有效模拟岩体内部真实的裂隙几何形态的基于分解-重构建立裂隙岩体三维数值模型的方法。
本发明所采用的技术方案是:
一种基于分解-重构建立裂隙岩体三维数值模型的方法,包括如下步骤:
1)在岩石表面贴上用于坐标识别的参考点:对岩体试样的表面进行清理,贴上参考点;
2)构建岩石块体表面的三维数值模型:利用三维扫描仪对岩体表面进行扫描,获得裂隙岩体的三维表面形态的测量数据,利用逆向工程软件完成对扫描数据的拼合,获得岩体的三维数值模型,记录参考点的坐标;
3)岩石块体的分解:沿着岩体内部的天然贯通裂隙将块体分解,分解后获得含有内部裂隙面的各块体,保留块体上的参考点;
4)构建各个岩石块体的三维数值模型:利用三维扫描仪对分解的块体进行逐个扫描,同步骤2)获得三维表面数值模型,记录扫描的参考点坐标;
5)岩体三维数值模型的重构:利用分解之前参考点的坐标对获得的三维数值模型进行坐标转化拼接,即获得重构的裂隙岩体三维数值模型。
所述的参考点为特制的可用于坐标识别的参考点,三维扫描仪扫描之后建立三维数值模型,参考点可被识别和标记出。
所述步骤2)中,所述的三维扫描仪在扫描过程中,收集的数据即为点云的三维坐标,因此,只要参考点可被三维扫描仪识别,即可记录参考点的坐标。
所述步骤1)中,所述的贴参考点需要注意的是,在岩体表面贴参考点,不得覆盖表面裂隙,且参考点为了便于块体的重构拼接,不能规则分布,在分解块体时,不得对参考点进行任何的移动、撕毁。
优选的,所述三维扫描仪为高精度结构光三维扫描仪,亦可是其他高精度三维扫描仪。
所述逆向工程软件及参考点坐标转化拼接,需要用到数据采集软件、数据分析软件和点云处理软件。
更进一步,获得的三维数值模型包含了岩体内部各条裂隙的表面形态和空间坐标数据。
由于采用了上述技术方案,本发明所取得有益效果为:
1)本发明中建立的裂隙岩体三维模型由真实的裂隙岩体所得,可完全还原岩体中裂隙的三维空间分布和表面形态。裂隙岩体三维模型可用于数值模拟,和实验结果作对比分析。可基于实测数据改变裂隙长度、开度、粗糙度等参数以及模型尺寸,实施多参数分析。
2)本发明中运用参考点的相对坐标转化拼接分解的块体,此方法在建立裂隙岩体三维模型中第一次试用。坐标转化的优点是可将分解后散乱的块体按原状拼接,且拼接后得到的模型可还原真实的内部裂隙面发育情况。
3)利用高精度三维扫描仪,裂隙面的模拟效果好,具有很好的操作性和实用性。
附图说明
图1为贴了参考点的裂隙岩体以及沿裂隙分解块体的式样。
图2为扫描前裂隙岩体的三维数值模型。
图3为重构后的裂隙岩体三维数值模型。
其中,1、参考点,2、分解的块体。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1~图3,一种基于分解-重构建立裂隙岩体三维数值模型的方法,包括如下步骤:
1)现场采集天然裂隙岩体试样,对岩体试样的表面进行清理,贴上参考点(如附图1)。参考点如附图1上所示,不能规则分布,不得覆盖裂隙。
2)利用三维扫描仪对岩体表面进行扫描,获得裂隙岩体各个表面的测量数据,利用逆向工程软件完成对扫描数据的拼合,获得裂隙岩体的三维表面数值模型,记录参考点的坐标。即附图2所示的三维模型,对岩体表面的扫描只能显示表面裂缝的情况,不能反映内部裂隙的发育。
3)沿着岩体内部的天然裂隙将块体分解,分解后获得含有内部裂隙面的各块体,保留块体上的参考点。
4)利用三维扫描仪对分解的块体进行逐个扫描,同步骤2)获得三维数值模型,记录扫描的参考点坐标,利用分解之前参考点的坐标对获得的三维数值模型进行坐标转化拼接,即获得重构的三维数值模型。附图3为重构后裂隙岩体的三维数值模型,可以精确反映内部裂隙的表面形态和空间分布情况。
所述的参考点为特制的可用于坐标识别的参考点,三维扫描仪扫描之后建立三维数值模型,参考点可被识别和标记出。
所述步骤2)中,所述的,三维扫描仪在扫描过程中,收集的数据即为点云的三维坐标,因此,只要参考点可被三维扫描仪识别,即可记录参考点的坐标。
所述步骤1)中,所述的贴参考点需要注意的是,在岩体表面贴参考点,不得覆盖表面裂隙,且参考点为了便于块体的重构拼接,不能规则分布,在分解块体时,不得对参考点进行任何的移动、撕毁。
优选的,所述三维扫描仪为高精度结构光三维扫描仪,亦可是其他高精度三维扫描仪。
所述逆向工程软件及参考点坐标转化拼接,需要用到数据采集软件、数据分析软件和点云处理软件。
更进一步,获得的三维数值模型包含了岩体内部各条裂隙的表面形态和空间坐标数据。
Claims (5)
1.一种基于分解-重构建立裂隙岩体三维数值模型的方法,其特征在于:所述方法包括如下步骤:
1)在岩石表面贴上用于坐标识别的参考点:对岩体试样的表面进行清理,贴上参考点;
2)构建岩石块体表面的三维数值模型:利用三维扫描仪对岩体表面进行扫描,获得裂隙岩体的三维表面形态的测量数据,利用逆向工程软件完成对扫描数据的拼合,获得岩体的三维数值模型,记录参考点的坐标;
3)岩石块体的分解:沿着岩体内部的天然贯通裂隙将块体分解,分解后获得含有内部裂隙面的各块体,保留块体上的参考点;
4)构建各个岩石块体的三维数值模型:利用三维扫描仪对分解的块体进行逐个扫描,同步骤2)获得三维表面数值模型,记录扫描的参考点坐标;
5)岩体三维数值模型的重构:利用分解之前参考点的坐标对获得的三维数值模型进行坐标转化拼接,即获得重构的裂隙岩体三维数值模型。
2.如权利要求1所述的基于分解-重构建立裂隙岩体三维数值模型的方法,其特征在于:所述步骤2)中,所述三维扫描仪在扫描过程中,收集的数据即为点云的三维坐标,因此,只要参考点可被三维扫描仪识别,即可记录参考点的坐标。
3.如权利要求1或2所述的基于分解-重构建立裂隙岩体三维数值模型的方法,其特征在于:所述步骤1)中,在岩体表面贴参考点,不得覆盖表面裂隙,且参考点为了便于块体的重构拼接,不能规则分布,在分解块体时,不得对参考点进行任何的移动、撕毁。
4.如权利要求1或2所述的基于分解-重构建立裂隙岩体三维数值模型的方法,其特征在于:所述步骤5)中,逆向工程软件及参考点坐标转化拼接,需要用到数据采集软件、数据分析软件和点云处理软件。
5.如权利要求4所述的基于分解-重构建立裂隙岩体三维数值模型的方法,其特征在于:所述步骤5)中,获得的三维数值模型包含岩体内部各条裂隙的表面形态和空间坐标数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610600268.8A CN106228559B (zh) | 2016-07-26 | 2016-07-26 | 基于分解-重构建立裂隙岩体三维数值模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610600268.8A CN106228559B (zh) | 2016-07-26 | 2016-07-26 | 基于分解-重构建立裂隙岩体三维数值模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106228559A true CN106228559A (zh) | 2016-12-14 |
CN106228559B CN106228559B (zh) | 2019-07-12 |
Family
ID=57533429
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610600268.8A Active CN106228559B (zh) | 2016-07-26 | 2016-07-26 | 基于分解-重构建立裂隙岩体三维数值模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106228559B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106803281A (zh) * | 2017-02-17 | 2017-06-06 | 武汉大学 | 一种开挖边坡岩块信息提取与三维重构方法 |
CN106844845A (zh) * | 2016-12-15 | 2017-06-13 | 长江水利委员会长江科学院 | 考虑弯扭贡献效应的三维时效破裂模型构建方法 |
CN107463740A (zh) * | 2017-07-27 | 2017-12-12 | 中南大学 | 考虑中间主应力效应的岩石类材料真三轴试验数值模拟方法 |
CN110415167A (zh) * | 2019-08-02 | 2019-11-05 | 山东科技大学 | 一种基于数字图像技术的粗糙面裂隙生成方法及试验系统 |
CN112085825A (zh) * | 2020-09-21 | 2020-12-15 | 华北科技学院 | 一种三维岩石裂隙的分形量化方法 |
CN112211646A (zh) * | 2020-09-16 | 2021-01-12 | 山东大学 | 一种适合隧道工程岩体结构精细化重构的方法及系统 |
CN112763288A (zh) * | 2020-12-22 | 2021-05-07 | 清华大学 | 基于3d打印的任意裂缝组合的岩体模型制作方法 |
CN114299239A (zh) * | 2022-02-10 | 2022-04-08 | 煤炭科学研究总院有限公司 | 一种对岩体天然裂隙结构的张量确定方法及设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1793780A (zh) * | 2005-12-05 | 2006-06-28 | 天津大学 | 基于拼接靶的大型三维形体形貌测量拼接方法和装置 |
CN102222359A (zh) * | 2011-05-24 | 2011-10-19 | 中国石油天然气股份有限公司 | 一种岩心三维孔隙结构重构方法 |
CN103114848A (zh) * | 2013-01-18 | 2013-05-22 | 西南石油大学 | 一种基于岩心测量的地层裂缝空间重构方法 |
CN104121864A (zh) * | 2014-08-11 | 2014-10-29 | 卢渊 | 一种岩石错动裂缝面的间隙评价方法 |
CN104616350A (zh) * | 2015-02-09 | 2015-05-13 | 西南石油大学 | 缝洞型碳酸盐油藏三维物理模型建立方法 |
CN104819874A (zh) * | 2015-03-19 | 2015-08-05 | 绍兴文理学院 | 一种岩石相似表面形态模型结构面制作方法 |
-
2016
- 2016-07-26 CN CN201610600268.8A patent/CN106228559B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1793780A (zh) * | 2005-12-05 | 2006-06-28 | 天津大学 | 基于拼接靶的大型三维形体形貌测量拼接方法和装置 |
CN102222359A (zh) * | 2011-05-24 | 2011-10-19 | 中国石油天然气股份有限公司 | 一种岩心三维孔隙结构重构方法 |
CN103114848A (zh) * | 2013-01-18 | 2013-05-22 | 西南石油大学 | 一种基于岩心测量的地层裂缝空间重构方法 |
CN104121864A (zh) * | 2014-08-11 | 2014-10-29 | 卢渊 | 一种岩石错动裂缝面的间隙评价方法 |
CN104616350A (zh) * | 2015-02-09 | 2015-05-13 | 西南石油大学 | 缝洞型碳酸盐油藏三维物理模型建立方法 |
CN104819874A (zh) * | 2015-03-19 | 2015-08-05 | 绍兴文理学院 | 一种岩石相似表面形态模型结构面制作方法 |
Non-Patent Citations (6)
Title |
---|
CAICHU XIA 等: "Numerical method for estimating void spaces of rock joints and the evolution of void spaces under different contact states", 《JOURNAL OF GEOPHYSICS AND ENGINEERING》 * |
LANRU JING: "Block system construction for three-dimensional discrete element models of fractured rocks", 《INTERNATIONAL JOURNAL OF ROCK MECHANICS AND MINING SCIENCES》 * |
LIANGCHAO ZOU 等: "Roughness decomposition and nonlinear fluid flow in a single rock fracture", 《INTERNATIONALJOURNALOFROCKMECHANICS&MININGSCIENCES》 * |
刘泉声 等: "多场耦合作用下岩体裂隙扩展演化关键问题研究", 《岩土力学》 * |
朱贵娜: "复杂岩体三维岩石块体系统构建方法及程序研制", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑(月刊) 》 * |
黄曼 等: "系列尺度岩石结构面相似表面模型制作的逆向控制技术研究", 《岩土力学》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106844845A (zh) * | 2016-12-15 | 2017-06-13 | 长江水利委员会长江科学院 | 考虑弯扭贡献效应的三维时效破裂模型构建方法 |
CN106803281A (zh) * | 2017-02-17 | 2017-06-06 | 武汉大学 | 一种开挖边坡岩块信息提取与三维重构方法 |
CN107463740A (zh) * | 2017-07-27 | 2017-12-12 | 中南大学 | 考虑中间主应力效应的岩石类材料真三轴试验数值模拟方法 |
CN107463740B (zh) * | 2017-07-27 | 2020-09-08 | 中南大学 | 考虑中间主应力效应的岩石类材料真三轴试验数值模拟方法 |
CN110415167A (zh) * | 2019-08-02 | 2019-11-05 | 山东科技大学 | 一种基于数字图像技术的粗糙面裂隙生成方法及试验系统 |
CN110415167B (zh) * | 2019-08-02 | 2020-03-06 | 山东科技大学 | 一种基于数字图像技术的粗糙面裂隙生成方法及试验系统 |
CN112211646A (zh) * | 2020-09-16 | 2021-01-12 | 山东大学 | 一种适合隧道工程岩体结构精细化重构的方法及系统 |
CN112211646B (zh) * | 2020-09-16 | 2021-10-29 | 山东大学 | 一种适合隧道工程岩体结构精细化重构的方法及系统 |
CN112085825A (zh) * | 2020-09-21 | 2020-12-15 | 华北科技学院 | 一种三维岩石裂隙的分形量化方法 |
CN112085825B (zh) * | 2020-09-21 | 2021-12-14 | 华北科技学院 | 一种三维岩石裂隙的分形量化方法 |
CN112763288A (zh) * | 2020-12-22 | 2021-05-07 | 清华大学 | 基于3d打印的任意裂缝组合的岩体模型制作方法 |
CN114299239A (zh) * | 2022-02-10 | 2022-04-08 | 煤炭科学研究总院有限公司 | 一种对岩体天然裂隙结构的张量确定方法及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN106228559B (zh) | 2019-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106228559A (zh) | 基于分解‑重构建立裂隙岩体三维数值模型的方法 | |
CN102222359B (zh) | 一种岩心三维孔隙结构重构方法 | |
Kufner et al. | Deep India meets deep Asia: Lithospheric indentation, delamination and break-off under Pamir and Hindu Kush (Central Asia) | |
CN110415167B (zh) | 一种基于数字图像技术的粗糙面裂隙生成方法及试验系统 | |
CN106289901B (zh) | 一种基于分解-重构制作裂隙岩体模型的方法 | |
CN106501090B (zh) | 用于水力压裂模拟实验的裂缝表征方法 | |
CN108507842A (zh) | 一种岩石三维自然结构面的制作方法 | |
CN104268940B (zh) | 基于ct扫描图像的mems结构重构与检测方法 | |
CN102855658A (zh) | 基于气象雷达基数据三维重构方法 | |
CN101364311A (zh) | 大规模城市仿真中快速自动建立模型的方法 | |
CN104819874A (zh) | 一种岩石相似表面形态模型结构面制作方法 | |
Xiao et al. | Implementation of manifold coverage for 3D rock fracture network modeling and its application in rock permeability prediction | |
CN109376395A (zh) | 一种建筑设计管理方法及系统 | |
CN107392873A (zh) | 物体虚拟修复方法、物体修复方法及物体虚拟展示系统 | |
Jia et al. | Exploration on the characteristics of 3D crack network expansion induced by hydraulic fracturing: a hybrid approach combining experiments and algorithms | |
CN105574930B (zh) | 一种汉字立体化数据采集与处理系统 | |
CN104121864B (zh) | 一种岩石错动裂缝面的间隙评价方法 | |
CN104318619A (zh) | 面向无损检测的自适应压缩感知的重建方法 | |
CN104881564A (zh) | 结构面粗糙度系数尺寸效应概率密度函数模型的构建方法 | |
CN109559377B (zh) | 一种利用经纬度多边形切割生成三维地图模型的方法 | |
CN114925567A (zh) | 岩石力学参数获取方法及其装置、计算机设备、存储介质 | |
CN105717021A (zh) | 一种定量确定裂隙岩体渗透各向异性的方法 | |
Tian | 3D modeling and digital preservation of ancient architectures based on autoCAD and 3Dmax | |
CN103218853A (zh) | 一种作物单根可变形建模方法 | |
Verdiani et al. | E. Stone, an archive for the Sardinia monumental witnesses |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |