CN107239629B - 一种岩石结构面实验室合理尺寸确定的分形维数分析方法 - Google Patents
一种岩石结构面实验室合理尺寸确定的分形维数分析方法 Download PDFInfo
- Publication number
- CN107239629B CN107239629B CN201710453575.2A CN201710453575A CN107239629B CN 107239629 B CN107239629 B CN 107239629B CN 201710453575 A CN201710453575 A CN 201710453575A CN 107239629 B CN107239629 B CN 107239629B
- Authority
- CN
- China
- Prior art keywords
- structural
- fractal dimension
- structural surface
- laboratory
- size
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
技术领域
本发明涉及岩石结构面尺寸效应的定量化表征,特别涉及一种岩石结构面实验室合理尺寸确定的分形维数分析方法。
背景技术
大量试验表明,岩石结构面的力学性质是存在尺寸效应的,而这一现象主要源于结构面粗糙度的尺寸效应。因此,在实验室如何利用结构面的粗糙信息规律确定合理的岩石试件尺寸具有重要的现实意义。目前表示结构面粗糙度的方法主要有统计参数表征法,分形维数描述法。分形几何是描述自然界不规则几何体的有效方法,因此,分形维数描述结构面粗糙度的成果较多。然而,在应用分形维数研究结构面的尺寸效应规律方面,目前仅局限于采用结构面某一或某几条剖面线的分形维数来描述,这样存在以偏概全的缺陷。因此,为了克服上述不足,本发明应用三维扫描仪获取岩石结构面形貌数据,采用改进投影覆盖法计算不同尺寸大小结构面的分形维数,分析其变化规律,给出其函数表达式,进而提出采用函数斜率倾角确定实验室内同类型结构面的合理尺寸的新方法。
发明内容
本发明的目的是克服现有技术的不足,提供一种更为精确的岩石结构面实验室合理尺寸的确定方法。
本发明提供了一种岩石结构面实验室合理尺寸确定的分形维数分析方法,包括以下步骤:
根据需要按照不同尺寸划分方案,将结构面划分为大小不同的区域;
计算各区域结构面的分形维数值;
研究大小不同结构面分形维数的变化规律,给出各尺寸划分方案下结构面尺寸效应的函数关系式,用对数函数关系表示,所述的a、b分别为系数,D是分形维数,L0为基础结构面边长,L0=32mm,L为所计算结构面的边长,L>L0;
比较各方案结构面尺寸效应函数关系式,以系数a的绝对值最大的函数关系式作为确定实验室内岩石结构面合理尺寸的计算依据;
所述的计算各区域结构面的分形维数值的步骤具体为:
(1)获取该区域结构面的三维形貌数据,包括结构面各点的高度信息,所述高度信息是指该点与结构面内最低点的落差;
(3)将结构面划分成δ-1×δ-1个小网格,采用函数rand()生成随机数,判断所生成随机数的奇偶性,根据随机数的奇偶性选择小网格中三角形的不同划分方案,按照海伦公式计算每个小网格的面积Ai(δ),其中i=1,2,…,δ-1×δ-1,进而计算结构面的总面积
(4)计算点对{ln(AT(δ)/AT0),ln(δ)},其中AT0表示结构面的直观面积,即投影面积,用结构面边长乘以边长来计算;
(5)按照最小二乘法分别对不同尺度变量δ下的{ln(AT(δ)/AT0),ln(δ)}的数据点对进行拟合,其斜率记为β,则结构面的分形维数为2-β;
(6)重复上述(2)至(5)步,直至达到抽样次数要求,将抽样数据置入数据处理软件EXCEL中,利用函数NORMDIST求解累积概率值,通过累积概率值判定抽样结果是否满足3σ法则,若满足,则可给出结构面分形维数D概率密度分布函数,并将此抽样结果的均值作为所研究结构面分形维数的精准计算值。
本发明的有益效果是提出了应用分形维数随岩石结构面大小的变化规律确定其实验室内合理尺寸的新方法,克服了通过计算岩石结构面二维剖面线分形维数来研究结构面尺寸效应片面性的缺点,使得岩石结构面的尺寸计算更为精确,为实验室内岩石结构面合理尺寸的确定提供了理论依据。
附图说明
图1为本发明岩石结构面实验室内合理尺寸确定分形维数分析流程图;
图2为本发明岩石各区域结构面分形维数计算流程图;
图3为本发明小网格中三角形划分方案图;
图4为本发明岩石结构面形貌及尺寸划分方案图;
图5为本发明岩石结构面分形维数随尺寸大小变化的发展趋势图(方案3)。
图6为本发明分形维数D所确定下的合理尺寸(214.49mm)
具体实施方式
具体实施方式一:结合图1来说明本实施方式所述的一种岩石结构面实验室合理尺寸确定的分形维数分析方法,包括以下步骤:
根据需要按照不同尺寸划分方案,将结构面划分为大小不同的区域;
计算各区域结构面的分形维数值;
研究大小不同结构面分形维数的变化规律,给出各尺寸划分方案下结构面尺寸效应的函数关系式,用对数函数关系表示,所述的a、b分别为系数,D是分形维数,L0为基础结构面边长,L0=32mm,L为所计算结构面的边长,L>L0;
比较各方案结构面尺寸效应函数关系式,以系数a的绝对值最大的函数关系式作为确定实验室内岩石结构面合理尺寸的计算依据;
具体实施方式二:本实施方式是对具体实施方式一所述的岩石结构面实验室合理尺寸确定的分形维数分析方法的进一步限定,如图2所示,所述的计算各区域结构面的分形维数的步骤为:
(1)获取该区域结构面的三维形貌数据,包括结构面各点的高度信息,所述高度信息是指该点与结构面内最低点的落差;
(3)将结构面划分成δ-1×δ-1个小网格,采用函数rand()生成随机数,判断所生成随机数的奇偶性,根据随机数的奇偶性选择小网格中三角形的不同划分方案,按照海伦公式计算每个小网格的面积Ai(δ),其中i=1,2,…,δ-1×δ-1,进而计算结构面的总面积
(4)计算点对{ln(AT(δ)/AT0),ln(δ)},其中AT0表示结构面的直观面积,即投影面积,用结构面边长乘以边长来计算;
(5)按照最小二乘法分别对不同尺度变量δ下的{ln(AT(δ)/AT0),ln(δ)}的数据点对进行拟合,其斜率记为β,则结构面的分形维数为2-β;
(6)重复上述(2)至(5)步,直至达到抽样次数要求,将抽样数据置入数据处理软件EXCEL中,利用函数NORMDIST求解累积概率值,通过累积概率值判定抽样结果是否满足3σ法则,若满足,则可给出结构面分形维数D概率密度分布函数,并将此抽样结果的均值作为所研究结构面分形维数的精准计算值。
下面通过具体实施例,对本发明做进一步详细说明,应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本实施例计算了一个天然红砂岩结构面的分形维数,该结构面规格长×宽为1024mm×1024mm,以下结合实施步骤详细说明本发明的具体实施方式:
1.结构面三维形貌信息数据获取
采用三维扫描仪EinScan-S获取结构面形貌信息数据。采集的数据保存于Excel中,以便下一步计算分析。
2.结构面分形维数D的计算
按照图1计算步骤,每个小网格三角形划分采用图3所示两种方案,在计算过程中随机选择其中一种,对图4结构面的分形维数D进行计算。本次结构面尺寸大小划分方案如图4所示,共进行5种不同划分方案下的分形维数计算。每种方案划分为5个尺寸不同的结构面,每个结构面共抽样计算了30次。表1是方案3结构面维数30次抽样计算结果汇总表。因数据较多,其余4种方案计算结果不再列表。
表1 方案3结构面维数30次抽样计算结果汇总表
3.岩石结构面实验室内合理尺寸确定
表2 各方案各尺寸分形维数的统计参数表
比较各关系式系数a的绝对值,方案3的函数关系式系数a的绝对值最大,为0.117。由此,用函数关系式f(D)=-0.117ln(L/L0)+2.4984作为确定实验室内该类型岩石结构面合理尺寸的计算依据。本实施例中通过和验证后得出判定评判指标K为1,求解曲线某一点斜率倾角等于1°时的L值,即求解得L的值为214.49mm。若用分形维数D评估,室内该类型结构面合理尺寸为214.49mm。
本实施例仅进行了图4各方案各尺寸结构面分形维数的30次抽样计算,若要获得更精确结果,可进行更多次抽样计算。
以上所述,仅为本发明的具体实施方式和实施例,本发明的保护范围并不局限于此,对于本领域普通技术人员来说,可以根据上述说明不加创造性地加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (2)
1.一种岩石结构面实验室内合理尺寸确定的分形维数分析方法,其特征在于,包括以下步骤:
根据需要按照不同尺寸划分方案,将结构面划分为大小不同的区域;
计算各区域结构面的分形维数值;
研究大小不同结构面分形维数的变化规律,给出各尺寸划分方案下结构面尺寸效应的函数关系式,用对数函数关系表示,所述的a、b分别为系数,D是分形维数,L0为基础结构面边长,L0=32mm,L为所计算结构面的边长,L>L0;
比较各方案结构面尺寸效应函数关系式,以系数a的绝对值最大的函数关系式作为确定实验室内岩石结构面合理尺寸的计算依据;
2.根据权利要求1所述的一种岩石结构面实验室内合理尺寸确定的分形维数分析方法,其特征在于,所述的计算各区域结构面的分形维数值的步骤为:
(1)获取岩石结构面的三维形貌数据,包括结构面各点的高度信息,所述高度信息是指该点与结构面内最低点的落差;
(3)将结构面划分成δ-1×δ-1个小网格,采用函数rand()生成随机数,判断所生成随机数的奇偶性,根据随机数的奇偶性选择小网格中三角形的不同划分方案,按照海伦公式计算每个小网格的面积Ai(δ),其中i=1,2,…,δ-1×δ-1,进而计算结构面的总面积
(4)计算点对{ln(AT(δ)/AT0),ln(δ)},其中AT0表示结构面的直观面积,即投影面积,用结构面边长乘以边长来计算;
(5)按照最小二乘法分别对不同尺度变量δ下的{ln(AT(δ)/AT0),ln(δ)}的数据点对进行拟合,其斜率记为β,则结构面的分形维数为2-β;
(6)重复上述(2)至(5)步,直至达到抽样次数要求,将抽样数据置入数据处理软件EXCEL中,利用函数NORMDIST求解累积概率值,通过累积概率值判定抽样结果是否满足3σ法则,若满足,则给出结构面的分形维数的概率密度分布函数,并将此抽样结果的均值作为所研究结构面的分形维数的精准计算值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710453575.2A CN107239629B (zh) | 2017-06-15 | 2017-06-15 | 一种岩石结构面实验室合理尺寸确定的分形维数分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710453575.2A CN107239629B (zh) | 2017-06-15 | 2017-06-15 | 一种岩石结构面实验室合理尺寸确定的分形维数分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107239629A CN107239629A (zh) | 2017-10-10 |
CN107239629B true CN107239629B (zh) | 2020-09-22 |
Family
ID=59986320
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710453575.2A Expired - Fee Related CN107239629B (zh) | 2017-06-15 | 2017-06-15 | 一种岩石结构面实验室合理尺寸确定的分形维数分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107239629B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110415283B (zh) * | 2019-07-03 | 2022-04-01 | 绍兴文理学院 | 解析岩体结构面各向异性尺寸效应特征的分形评估方法 |
CN110570401B (zh) * | 2019-08-19 | 2020-04-21 | 成都理工大学 | 一种基于图像识别的斜坡碎裂岩体信息获取方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120029888A1 (en) * | 2010-07-30 | 2012-02-02 | Jeng-Haur Horng | Topography shaping apparatus for forming surfaces of low friction coefficient |
CN103926156A (zh) * | 2014-04-10 | 2014-07-16 | 内蒙古科技大学 | 一种三维岩石结构面剪切力学行为多重分形评估方法 |
CN105701274A (zh) * | 2016-01-04 | 2016-06-22 | 中国矿业大学 | 一种岩土参数三维局部平均随机场样本的生成方法 |
-
2017
- 2017-06-15 CN CN201710453575.2A patent/CN107239629B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120029888A1 (en) * | 2010-07-30 | 2012-02-02 | Jeng-Haur Horng | Topography shaping apparatus for forming surfaces of low friction coefficient |
CN103926156A (zh) * | 2014-04-10 | 2014-07-16 | 内蒙古科技大学 | 一种三维岩石结构面剪切力学行为多重分形评估方法 |
CN105701274A (zh) * | 2016-01-04 | 2016-06-22 | 中国矿业大学 | 一种岩土参数三维局部平均随机场样本的生成方法 |
Non-Patent Citations (3)
Title |
---|
Direct Estimation of the Fractal Dimensions of a Fracture Surface of Rock;H.W.Zhou 等;《Surface Review and Letters》;20030216;第10卷(第5期);第751-762页 * |
基于多重分形特征的岩体结构面剪切强度研究;陈世江 等;《岩土力学》;20150310;第36卷(第3期);第703-710、718页 * |
岩体结构面粗糙度系数定量表征研究进展;陈世江 等;《力学学报》;20170331;第49卷(第2期);第239-256页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107239629A (zh) | 2017-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110409369B (zh) | 边坡开挖数字化施工与质量控制方法 | |
Shi et al. | Adaptive simplification of point cloud using k-means clustering | |
Zhou et al. | A comparison of fractal dimension estimators based on multiple surface generation algorithms | |
CN103701466A (zh) | 基于特征保留的散乱点云压缩算法 | |
CN103149348B (zh) | 基于小波多尺度分析的先进陶瓷磨削表面损伤的评价方法 | |
CN109598015B (zh) | 碎裂结构岩体碎裂程度分级评价方法 | |
CN107239629B (zh) | 一种岩石结构面实验室合理尺寸确定的分形维数分析方法 | |
CN102749053A (zh) | 基于三维可视化和蒙特卡罗方法的体积测量方法 | |
Chu et al. | Quantification of soil random roughness and surface depression storage: Methods, applicability, and limitations | |
CN107194114B (zh) | 一种岩石结构面多重分形参数的计算方法 | |
WO2019136870A1 (zh) | 一种利用平行光评价岩石结构面起伏程度的计算方法 | |
CN107292017B (zh) | 岩石结构面实验室合理尺寸确定的多重分形参数分析方法 | |
CN107292100B (zh) | 一种岩石结构面分形维数的计算方法 | |
CN109543236B (zh) | 岩体结构面粗糙度统计样本数确定方法 | |
CN109583003B (zh) | 基于面裂多边形的碎裂结构岩体碎裂程度量化方法 | |
CN112923900A (zh) | 一种基于rs技术的建筑物不均匀沉降的判断方法及装置 | |
CN112346139A (zh) | 一种重力数据多层等效源延拓与数据转换方法 | |
CN109614630B (zh) | 基于迹线节点的碎裂结构岩体碎裂程度量化方法 | |
CN113378909B (zh) | 一种岩石节理面粗糙度系数的分级表征方法、设备及介质 | |
CN104484900A (zh) | 基于Delaunay三角网的气象要素插值评估方法 | |
CN116882255B (zh) | 一种基于傅里叶级数随机生成多孔介质模型的方法及系统 | |
CN113239600B (zh) | 一种复杂岩体二维随机网络模型构建方法 | |
CN117132738B (zh) | 一种球面离散格网多尺度等距模式量化方法和系统 | |
CN117173345A (zh) | 基于分形理论和蒙特卡罗方法的粗糙表面形貌模拟方法 | |
CN112085825B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200922 Termination date: 20210615 |