一种岩石结构面实验室合理尺寸确定的分形维数分析方法与流程

本发明涉及岩石结构面尺寸效应的定量化表征，特别涉及一种岩石结构面实验室合理尺寸确定的分形维数分析方法。

背景技术：

大量试验表明，岩石结构面的力学性质是存在尺寸效应的，而这一现象主要源于结构面粗糙度的尺寸效应。因此，在实验室如何利用结构面的粗糙信息规律确定合理的岩石试件尺寸具有重要的现实意义。目前表示结构面粗糙度的方法主要有统计参数表征法，分形维数描述法。分形几何是描述自然界不规则几何体的有效方法，因此，分形维数描述结构面粗糙度的成果较多。然而，在应用分形维数研究结构面的尺寸效应规律方面，目前仅局限于采用结构面某一或某几条剖面线的分形维数来描述，这样存在以偏概全的缺陷。因此，为了克服上述不足，本发明应用三维扫描仪获取岩石结构面形貌数据，采用改进投影覆盖法计算不同尺寸大小结构面的分形维数，分析其变化规律，给出其函数表达式，进而提出采用函数斜率倾角确定实验室内同类型结构面的合理尺寸的新方法。

技术实现要素：

本发明的目的是克服现有技术的不足，提供一种更为精确的岩石结构面实验室合理尺寸的确定方法。

本发明提供了一种岩石结构面实验室合理尺寸确定的分形维数分析方法，包括以下步骤：

根据需要按照不同尺寸划分方案，将结构面划分为大小不同的区域；

计算各区域结构面的分形维数值；

研究大小不同结构面分形维数的变化规律，给出各尺寸划分方案下结构面尺寸效应的函数关系式，用对数函数关系表示，所述的a、b分别为系数，d是分形维数，l0为基础结构面边长，l0＝32mm，l为所计算结构面的边长，l>l0；

比较各方案结构面尺寸效应函数关系式，以系数a的绝对值最大的函数关系式作为确定实验室内岩石结构面合理尺寸的计算依据；

当函数曲线某一点斜率倾角等于k°时，将该点值作为结构面的合理临界尺寸，亦即求解时，l的值，所述的k是评判指标，通过在实验室比较函数斜率较小时的两块岩石结构面的力学性能，验证后进行确定。

所述的计算各区域结构面的分形维数值的步骤具体为：

(1)获取该区域结构面的三维形貌数据，包括结构面各点的高度信息，所述高度信息是指该点与结构面内最低点的落差；

(2)尺度变量δ分别取值重复进行步骤(3)至(4)；

(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的绝对值最大的函数关系式作为确定实验室内岩石结构面合理尺寸的计算依据；

当函数曲线某一点斜率倾角等于k°时，将该点值作为结构面的合理临界尺寸，亦即求解时，l的值，所述的k是评判指标，通过在实验室比较函数斜率较小时的两块岩石结构面的力学性能，验证后进行确定。

具体实施方式二：本实施方式是对具体实施方式一所述的岩石结构面实验室合理尺寸确定的分形维数分析方法的进一步限定，如图2所示，所述的计算各区域结构面的分形维数的步骤为：

(1)获取该区域结构面的三维形貌数据，包括结构面各点的高度信息，所述高度信息是指该点与结构面内最低点的落差；

(2)尺度变量δ分别取值重复进行步骤(3)至(4)；

(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.岩石结构面实验室内合理尺寸确定

统计分析表1数据，其分形维数d与的关系为f(d)＝-0.117ln(l/l0)+2.4984；其发展趋势图见图5。图4五种方案结构面尺寸效应函数关系式见表2。

表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次抽样计算，若要获得更精确结果，可进行更多次抽样计算。

以上所述，仅为本发明的具体实施方式和实施例，本发明的保护范围并不局限于此，对于本领域普通技术人员来说，可以根据上述说明不加创造性地加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。