CN103886613A - 一种岩石结构面形貌各向异性评价方法 - Google Patents
一种岩石结构面形貌各向异性评价方法 Download PDFInfo
- Publication number
- CN103886613A CN103886613A CN201410142106.5A CN201410142106A CN103886613A CN 103886613 A CN103886613 A CN 103886613A CN 201410142106 A CN201410142106 A CN 201410142106A CN 103886613 A CN103886613 A CN 103886613A
- Authority
- CN
- China
- Prior art keywords
- represent
- structural plane
- value
- structural
- rock
- 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
Images
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明一种岩石结构面形貌各向异性评价方法。其步骤包括:⑴获取岩石结构面三维形貌数据;⑵按照实验变异函数公式计算结构面分析方向上的点对,并绘制随变化的发展趋势图;⑶确定变异函数模型并进行函数模型拟合,进而确定模型参数变程a与基台C的值;⑷应用量纲分析方法,用H(由a和C表示)表征结构面各向异性特征;⑸按照步骤⑵、⑶和⑷计算结构面上各个方向上的值,以此来分析岩石结构面的各向异性特征。本发明提供的分析方法可以有效解决岩石结构面形貌各向异性定量化表征的问题,方便准确,该方法为进一步全面评价岩石结构面剪切强度提供依据。
Description
技术领域
本发明涉及一种岩石结构面形貌各向异性评价方法,属于岩石力学领域范畴。
背景技术
大量研究表明,结构面的力学行为与水力特性存在明显的各向异性,而这一特性很大程度上取决于结构面形貌的各向异性。因此,深入研究岩石结构面各向异性的特征有很重要的理论意义与实际工程应用价值。国内外的学者在岩石结构面各向异性的研究方面做了大量的工作,主要应用地理信息系统的方法、分形理论描述的方法、累计功率谱密度指数描述、结构面不同方向上剖线的粗糙度分析方法等进行研究,并取得了丰硕的成果。尽管国内外学者在这方面做了很大的努力,但因岩石结构面的复杂性,其各向异性特征定量化表征仍是一个难点。
发明内容
本发明的目的在于提供一种岩石结构面形貌各向异性评价方法,该方法能够定量评价岩石结构面各向异性特征,进而为全面评价岩石结构面剪切强度提供依据。
技术解决方案:
本发明首先获取岩石结构面三维形貌数据,然后应用地质统计学原理,计算结构面各个方向上的变异函数,通过分析变异函数特征,提出用变异函数参数(基台和变程)表示岩石结构面的各向异性特征的方法。
用该方法描述岩石结构面的各向异性特征的具体步骤如下:
(1)获取岩石结构面数字图像,用图像灰度值表示岩石结构面的三维形貌数据;
(2)以岩石结构面的三维数据为基础,按照实验变异函数计算公式 ,计算结构面分析方向上的点对,并绘制随变化的发展趋势图,其中,为变异函数值,代表结构面点处的起伏高度,表示结构面内点的编号,表示分析方向上两点的间距,表示分析方向上结构面内间隔为的点对数;
(4)为了便于比较分析结构面的各向异性特征,应用量纲分析、归一化处理的方法,用单一指标H表示结构面形貌各向异性特征,,其中表示基台,表示变程,表示结构面起伏高程的平均值,表示结构面分析方向长度的一半,H值能够表示结构面的各向异性特征;
(5)结构面各个方向的H值通过重复步骤(2)—(4)获得。
本发明步骤3)中变异函数模型包括:球状模型、或对数模型、或指数模型、或高斯模型。
本发明还包括用图像灰度值表示岩石结构面的三维形貌特征,具体步骤为:
①将获得的结构面图像转化为256灰度图;
②在灰度图中0代表最暗,255代表最亮,在岩石结构面中灰度最小值代表最低点,灰度最大值表示最高点,按照公式对结构面中的灰度值进行比例缩放,并使其单位转化为像素,从而使缩放后的灰度值表示结构面的实际尺寸,其中,表示结构面图像长度,单位像素,表示结构面实际长度,单位mm,表示图像中最大灰度值,表示图像中最小灰度值,表示结构面内最高点,单位mm,表示结构面内最低点,单位mm,为缩放比例系数;
本发明还包括:分析方向的确定,具体步骤为:
①图像中的像素点是以矩阵的形式存储的,首先在矩阵中建立平面直角坐标系;
②在直角坐标系中选择两个像素点;
③根据两点在直角坐标系中的位置计算方向角度,计算结果为确定的分析方向。
本发明优点:
因岩石结构面的复杂性,其各向异性特征定量化表征仍是一个难点。本发明在获取结构面三维形貌数据的基础上,利用地质统计学中的变异函数定量分析岩石结构面的各向异性特征,研究结果表明,该方法是分析结构面各向异性特征的一种非常有效地方法。
附图说明
图1为本发明岩石结构面形貌各向异性分析流程图;
图2为本发明石膏试件结构面形貌图;
图3为本发明图像灰度值表示岩石结构面三维形貌特征流程图;
图4为本发明结构面图像数据点分布方式及分析方向示意图;
图5 为本发明结构面27°方向上点对统计示意图;
图6为本发明结构面上爬坡数据示意图;
图8为本发明结构面各向异性特征图;
图9为本发明结构面方向指示图。
具体实施方式
本实施例采用石膏粉材料人工制作类岩石结构面,以下结合附图详细说明本发明的具体实施方式:
见图1是岩石结构面形貌各向异性分析流程图。
结构面图像获取与用灰度值表示结构面三维形貌
采用佳能60D数码相机获取石膏结构面形貌图(见图2),将该结构面形貌图像进行图像灰度转换,比例缩放,而后将处理后的数据进行存储。图像灰度值表示岩石结构面三维形貌特征流程图见图3。该过程中比例缩放确定过程具体步骤如下:
结构面分析方向上变异函数模型参数计算
按照图1步骤(4)首先确定结构面分析方向,如图4所示,首先在结构面图像像素点矩阵中建立平面直角坐标系;其次在直角坐标系中选择两个像素点,如像素点1和像素点2;最后根据1、2两点在直角坐标系中的位置计算方向角度,即,,;这样图4中所有斜直线方向的点将参与计算分析。
按照实验变异函数计算公式,计算该方向上的点对。见图5结构面27°方向上点对统计示意图,当为1时,在该方向上逐行追踪扫描图像中的每一点,并记录符合=1的点的值,将参与两点间高程差的计算,这样的点对在示例中共有20对,同样方法,逐步增大,相应地统计出点的对数,这样能够计算不同的对应的的值,进而绘制随变化的发展趋势图(见图7所示)。该过程主要为了说明点对的统计,并未考虑在该方向上各点的分布特点,为了区分同一方位正反方向的各向异性特征,计算时采用沿分析方向上所有爬坡数据参与变异函数的计算,爬坡数据示意图见图6,此过程可以通过比较相邻点数值的大小,分析其是否参与计算。
按照图1步骤(5)确定变异函数模型,根据图7中随变化的发展趋势,前半段数据点较陡,后半段数据点相对较缓,该特征符合二级球状模型变异函数曲线特点,故采取二级球状模型变异函数进行拟合,拟合结果见图7拟合曲线,相应变异函数模型参数见图7指示部分。二级球状模型变异函数拟合具体步骤如下:
下述公式为二级球状模型变异函数模型计算公式。
表1 二级球状模型变异函数计算参数统计表
3.结构面各方向上变异特征参数H计算及各向异性分析
表2 32个方向上H值统计表
方向 ° | 0 | 9 | 18 | 27 | 45 | 63 | 72 | 81 | 90 | 99 | 108 |
变异性指标 | 0.003791 | 0.004679 | 0.003756 | 0.007213 | 0.018005 | 0.013808 | 0.013966 | 0.030331 | 0.036064 | 0.039219 | 0.039757 |
方向 ° | 117 | 135 | 153 | 162 | 171 | 180 | 189 | 198 | 207 | 225 | 243 |
变异性指标 | 0.038763 | 0.017666 | 0.035505 | 0.040319 | 0.022346 | 0.026571 | 0.012770 | 0.011237 | 0.011614 | 0.020210 | 0.024482 |
方向 ° | 252 | 261 | 270 | 279 | 288 | 297 | 315 | 333 | 342 | 351 | |
变异性指标 | 0.019851 | 0.023860 | 0.028232 | 0.025217 | 0.032419 | 0.028017 | 0.022091 | 0.008788 | 0.006852 | 0.013866 |
按照图1步骤(8),分析该结构面32个方向上H值,根据表2数据绘制的结构面各向异性特征图(见图8),图9是相关方向指示。
由实施例结果可见,提出的参数H能够很好的定量化表示岩石结构面的各向异性特征,分析结果可为全面了解岩石结构面力学性质,从而进一步为准确全面评价工程岩石强度提供依据。
Claims (5)
1.一种岩石结构面形貌各向异性评价方法,其特征在于,方法步骤包括:
(1)获取岩石结构面数字图像,用图像灰度值表示岩石结构面的三维形貌数据;
(2)以岩石结构面的三维数据为基础,按照实验变异函数计算公式 ,计算结构面分析方向上的点对,并绘制随变化的发展趋势图,其中,为变异函数值,代表结构面点处的起伏高度,表示结构面内点的编号,表示分析方向上两点的间距,表示分析方向上结构面内间隔为的点对数;
(4)为了便于比较分析结构面的各向异性特征,应用量纲分析、归一化处理的方法,用单一指标H表示结构面形貌各向异性特征,,其中表示基台,表示变程,表示结构面起伏高程的平均值,表示结构面分析方向长度的一半,H值能够表示结构面的各向异性特征;
(5)结构面各个方向的H值通过重复步骤(2)—(4)获得。
4.根据权利要求1所述岩石结构面形貌各向异性评价方法,其特征在于,还包括:分析方向的确定,具体步骤为:
①图像中的像素点是以矩阵的形式存储的,首先在矩阵中建立平面直角坐标系;
②在直角坐标系中选择两个像素点;
③根据两点在直角坐标系中的位置计算方向角度,计算结果为确定的分析方向。
5.根据权利要求1所述一种岩石结构面形貌各向异性评价方法,其特征在于,步骤3)中变异函数模型包括:球状模型、或对数模型、或指数模型、或高斯模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410142106.5A CN103886613B (zh) | 2014-04-10 | 2014-04-10 | 一种岩石结构面形貌各向异性评价分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410142106.5A CN103886613B (zh) | 2014-04-10 | 2014-04-10 | 一种岩石结构面形貌各向异性评价分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103886613A true CN103886613A (zh) | 2014-06-25 |
CN103886613B CN103886613B (zh) | 2016-08-31 |
Family
ID=50955486
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410142106.5A Expired - Fee Related CN103886613B (zh) | 2014-04-10 | 2014-04-10 | 一种岩石结构面形貌各向异性评价分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103886613B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105203593A (zh) * | 2015-10-14 | 2015-12-30 | 东南大学 | 基于混合料细观结构特性的沥青和集料导热系数反演方法 |
CN105466790A (zh) * | 2015-11-10 | 2016-04-06 | 内蒙古科技大学 | 一种各向异性特征的岩石结构面剪切强度评估方法 |
CN106546189A (zh) * | 2016-11-14 | 2017-03-29 | 绍兴文理学院 | 一种系列尺寸岩体结构面起伏幅度各向异性的评价方法 |
CN106971425A (zh) * | 2017-03-14 | 2017-07-21 | 天津大学 | 基于多面体网络模型的岩石块体群精细建模与识别方法 |
CN107036561A (zh) * | 2016-11-14 | 2017-08-11 | 绍兴文理学院 | 基于中智数函数的结构面粗糙度各向异性的近似表达方法 |
CN108446431A (zh) * | 2018-02-06 | 2018-08-24 | 中国地质大学(武汉) | 岩石结构面剪切强度各向异性评价方法 |
CN110362857A (zh) * | 2019-06-04 | 2019-10-22 | 绍兴文理学院 | 岩石结构面各向异性分布特征的理论模型解析方法 |
CN110362859A (zh) * | 2019-06-06 | 2019-10-22 | 绍兴文理学院 | 考虑不同方向起伏贡献率的三维结构面表面形态构建方法 |
CN110543854A (zh) * | 2019-09-05 | 2019-12-06 | 广东水电二局股份有限公司 | 一种基于图像处理的岩土结构检测方法 |
CN110941012A (zh) * | 2019-11-29 | 2020-03-31 | 北京化工大学 | 基于混合智能算法的弹性矢量波场分离方法及装置 |
CN112903565A (zh) * | 2021-02-01 | 2021-06-04 | 核工业北京地质研究院 | 考虑岩石裂隙内部几何特征的渗透率测定方法 |
CN114862783A (zh) * | 2022-04-27 | 2022-08-05 | 天津大学 | 一种表征岩石结构面吻合程度的方法及设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102254095A (zh) * | 2011-07-01 | 2011-11-23 | 合肥工业大学 | 基于多维分形克里格方法的成矿异常提取方法 |
-
2014
- 2014-04-10 CN CN201410142106.5A patent/CN103886613B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102254095A (zh) * | 2011-07-01 | 2011-11-23 | 合肥工业大学 | 基于多维分形克里格方法的成矿异常提取方法 |
Non-Patent Citations (4)
Title |
---|
G.GRASSELLI 等: "Constitutive law for the shear strength of rock joints based on three-dimensional surface parameters", 《INTERNATIONAL JOURNAL OF ROCK MECHANICS & MINING SCIENCES》, vol. 40, no. 1, 31 January 2003 (2003-01-31), pages 25 - 40 * |
张征 等: "岩土参数的变异性及其评价方法", 《土木工程学报》, vol. 28, no. 6, 31 December 1995 (1995-12-31), pages 43 - 51 * |
游志城 等: "基于三维激光扫描技术的结构面抗剪强度", 《岩石力学与工程学报》, vol. 33, no. 1, 31 January 2014 (2014-01-31), pages 3003 - 3008 * |
胡小荣 等: "岩土体的非均质性及力学参数的条件模拟赋值", 《岩石力学与工程学报》, vol. 21, no. 1, 31 January 2002 (2002-01-31), pages 13 - 17 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105203593A (zh) * | 2015-10-14 | 2015-12-30 | 东南大学 | 基于混合料细观结构特性的沥青和集料导热系数反演方法 |
CN105203593B (zh) * | 2015-10-14 | 2018-02-02 | 东南大学 | 一种基于混合料细观结构特性的沥青和集料导热系数反演方法 |
CN105466790B (zh) * | 2015-11-10 | 2017-12-29 | 内蒙古科技大学 | 一种各向异性特征的岩石结构面剪切强度评估方法 |
CN105466790A (zh) * | 2015-11-10 | 2016-04-06 | 内蒙古科技大学 | 一种各向异性特征的岩石结构面剪切强度评估方法 |
CN107036561B (zh) * | 2016-11-14 | 2019-05-31 | 绍兴文理学院 | 基于中智数函数的结构面粗糙度各向异性的近似表达方法 |
CN106546189A (zh) * | 2016-11-14 | 2017-03-29 | 绍兴文理学院 | 一种系列尺寸岩体结构面起伏幅度各向异性的评价方法 |
CN107036561A (zh) * | 2016-11-14 | 2017-08-11 | 绍兴文理学院 | 基于中智数函数的结构面粗糙度各向异性的近似表达方法 |
CN106971425A (zh) * | 2017-03-14 | 2017-07-21 | 天津大学 | 基于多面体网络模型的岩石块体群精细建模与识别方法 |
CN108446431A (zh) * | 2018-02-06 | 2018-08-24 | 中国地质大学(武汉) | 岩石结构面剪切强度各向异性评价方法 |
CN110362857B (zh) * | 2019-06-04 | 2023-05-16 | 绍兴文理学院 | 岩石结构面各向异性分布特征的理论模型解析方法 |
CN110362857A (zh) * | 2019-06-04 | 2019-10-22 | 绍兴文理学院 | 岩石结构面各向异性分布特征的理论模型解析方法 |
CN110362859A (zh) * | 2019-06-06 | 2019-10-22 | 绍兴文理学院 | 考虑不同方向起伏贡献率的三维结构面表面形态构建方法 |
CN110362859B (zh) * | 2019-06-06 | 2023-05-16 | 绍兴文理学院 | 考虑不同方向起伏贡献率的三维结构面表面形态构建方法 |
CN110543854A (zh) * | 2019-09-05 | 2019-12-06 | 广东水电二局股份有限公司 | 一种基于图像处理的岩土结构检测方法 |
CN110941012A (zh) * | 2019-11-29 | 2020-03-31 | 北京化工大学 | 基于混合智能算法的弹性矢量波场分离方法及装置 |
CN112903565A (zh) * | 2021-02-01 | 2021-06-04 | 核工业北京地质研究院 | 考虑岩石裂隙内部几何特征的渗透率测定方法 |
CN114862783A (zh) * | 2022-04-27 | 2022-08-05 | 天津大学 | 一种表征岩石结构面吻合程度的方法及设备 |
CN114862783B (zh) * | 2022-04-27 | 2022-12-23 | 天津大学 | 一种表征岩石结构面吻合程度的方法及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN103886613B (zh) | 2016-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103886613A (zh) | 一种岩石结构面形貌各向异性评价方法 | |
CN103884312B (zh) | 一种岩石结构面形貌尺寸效应分析方法 | |
CN103926156B (zh) | 一种三维岩石结构面剪切力学行为多重分形评估方法 | |
Mou et al. | Numerical simulation of the effects of building dimensional variation on wind pressure distribution | |
CN101101612B (zh) | 一种模拟田面微地形空间分布状况的方法 | |
CN105466790B (zh) | 一种各向异性特征的岩石结构面剪切强度评估方法 | |
CN104834806A (zh) | 结构面粗糙度系数尺寸效应取样代表性评价方法 | |
Vaughan et al. | Using fractal analysis to compare the characteristic complexity of nature and architecture: re-examining the evidence | |
CN101976444A (zh) | 一种基于像素类型的结构类似性图像质量客观评价方法 | |
CN105678757A (zh) | 一种物体位移测量方法 | |
CN109598015B (zh) | 碎裂结构岩体碎裂程度分级评价方法 | |
Napolitano et al. | Minimizing the adverse effects of bias and low repeatability precision in photogrammetry software through statistical analysis | |
CN101183460B (zh) | 彩色图像背景杂波量化方法 | |
CN103366045A (zh) | 基于格子Boltzmann的流体可视化仿真方法 | |
CN110716998B (zh) | 一种精细尺度人口数据空间化方法 | |
CN103268423A (zh) | 一种地理现象多点模拟空间尺度选择的方法 | |
CN110610539A (zh) | 地层曲面构建方法和装置、设备及储存介质 | |
CN110907327A (zh) | 一种考虑等效效率优化的复杂裂隙网络等效方法 | |
Xue et al. | Bias estimation and correction for triangle-based surface area calculations | |
CN110222368A (zh) | 一种利用二维切片计算岩心三维孔隙度和渗透率的方法 | |
CN111191393B (zh) | 基于离散余弦变换的水动力型滑坡非确定性评价方法及装置 | |
Alghalandis | Stochastic Modelling of Fractures | |
Salah | Geostatistical analysis of groundwater levels in the south Al Jabal Al Akhdar area using GIS | |
Gorard | The possible advantages of the mean absolute deviation'effect'size | |
Lorenz | Fractal geometry of architecture: Implementation of the box-counting method in a CAD-software |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160831 Termination date: 20190410 |
|
CF01 | Termination of patent right due to non-payment of annual fee |