一种岩石结构面多重分形参数的计算方法技术

本发明专利技术涉及一种对岩石结构面粗糙度进行定量化的多重分形参数计算方法，其步骤包括：获取岩石试件结构面的三维形貌数据；利用投影覆盖法将结构面划分成个小网格，生成随机数，根据随机数的奇偶性选择小网格中三角形的不同划分方案，计算每个小网格的面积和结构面的总面积、每个小网格的结构面的粗糙性概率、测度；利用最小二乘法拟合得到不同测距度对应的多重分形谱参数，计算广义分形维数、多重分形参数；重复抽样，使抽样结果满足3σ法则，给出结构面多重分形参数概率密度分布函数，并将此抽样结果的均值作为所研究结构面多重分形参数的精准计算值。

【技术实现步骤摘要】
一种岩石结构面多重分形参数的计算方法
本专利技术涉及一种岩石结构面粗糙度的定量化表征，特别涉及一种岩石结构面多重分形参数的计算方法。
技术介绍
岩石结构面的粗糙度对其力学性质有重要的影响。诸如结构面的剪切强度、节理裂隙的渗流特性均与结构面的粗糙度有着密切的关系。因此，岩石结构面粗糙度定量化表征是该领域研究的热点。多重分形参数表征是岩石结构面定量化描述的主要方法之一，而投影覆盖法又是计算岩石结构面多重分形参数的常用方法。2009年，采用向量分析法，对该方法中小网格面积的计算进行了改进，然而在计算时，由于每一个小网格中三角形划分方案不同，导致小网格表面积计算结果也不相同。如图3所示，小网格表面abcd的面积可由abc和acd两个三角形组成；也可由abd和bcd两个三角形组成。这样，小网格abcd采用不同的三角形划分方案，其多重分形参数计算结果是不同的，所以目前投影覆盖法依赖于划分方案，缺少一种精确地多重分形参数的计算方法。
技术实现思路
本专利技术的目的是克服上述不足，提供一种更为精确的岩石结构面多重分形参数的计算方法。本专利技术的一种岩石结构面多重分形参数的计算方法，该方法包括以下步骤：(1)获取岩石试件结构面的三维形貌数据，包括结构面各点的高度信息，所述高度信息是指该点与结构面内最低点的落差；(2)测度距q分别取值[-45,45]的整数，重复进行步骤(3)至(8)；(3)尺度变量δ分别取值重复进行步骤(4)至(6)；(4)将结构面划分成δ-1×δ-1个小网格，采用函数rand()生成随机数，判断所生成随机数的奇偶性，根据随机数的奇偶性选择小网格中三角形的不同划分方案，按照海伦公式计算每个小网格的面积Ai(δ)，其中i＝1,2,…,δ-1×δ-1，进而计算结构面的总面积(5)计算每个小网格的结构面粗糙性概率pi(δ)，(6)计算每个小网格的测度ui(q,δ)，进一步计算log(1/δ)以及(7)按照最小二乘法分别对不同尺度变量δ下的和的数据点对进行拟合，分别得到拟合函数的斜率，就是测度距为q时多重分形谱参数f(a)和a(q)的值；(8)根据测度距q对应的f(a)和a(q)值，计算测度距q对应的广义分形维数Dq，其中时，(9)计算结构面多重分形谱特征参数ΔD(q)和Δa(q)，ΔD(q)＝maxDq-minDq，Δa(q)＝maxa(q)-mina(q)；(10)重复上述(2)至(9)步，直至达到抽样次数要求，将抽样数据置入数据处理软件EXCEL中，利用函数NORMDIST求解累积概率值，通过累积概率值判定抽样结果是否满足3σ法则即抽样结果落在[μ-3σ，μ+3σ]范围内，若满足，则可给出结构面多重分形参数ΔD(q)和Δa(q)概率密度分布函数，并将此抽样结果的均值作为所研究结构面多重分形参数的精准计算值。本专利技术的有益效果是通过引入随机数，对投影覆盖法中小网格三角形划分方案进行随机选择，这样就使结构面的网格划分结果与实际结构形貌更加吻合，而且对随机数进行多次抽样计算，依据抽样结果满足3σ法则的原则，将抽样结果的均值作为所研究结构面多重分形参数的精准计算值，这样实现了对岩石结构面多重分形的精确计算，进而使岩石结构面粗糙度的定量化描述更加准确。附图说明图1是本专利技术中具体实施方式所述的岩石结构面多重分形参数计算方法的流程图。图2是本专利技术中对结构面小网格划分方案的示意图。图3是本专利技术中对结构面小网格中三角形划分方案示意图。图4是本专利技术具体实施例中的石膏结构面形貌图。图5是本专利技术具体实施例中石膏结构面多重分形参数首次抽样计算结果图。图6是本专利技术具体实施例中石膏结构面多重分形参数ΔD(q)频数直方图。图7是本专利技术具体实施例中石膏结构面多重分形参数Δa(q)频数直方图。图8是本专利技术具体实施例中石膏结构面多重分形参数ΔD(q)累积概率分布图。图9是本专利技术具体实施例中石膏结构面多重分形参数Δa(q)累积概率分布图。图10是本专利技术具体实施例中石膏结构面多重分形参数ΔD(q)概率密度函数分布图。图11是本专利技术具体实施例中石膏结构面多重分形参数Δa(q)概率密度函数分布图。具体实施方式具体实施方式：结合图1来说明本实施方式所述的一种岩石结构面多重分形参数的计算方法，包括以下步骤：(1)获取岩石试件结构面的三维形貌数据，包括结构面各点的高度信息，所述高度信息是指该点与结构面内最低点的落差；(2)测度距q分别取值[-45,45]的整数，重复进行步骤(3)至(8)；(3)尺度变量δ分别取值重复进行步骤(4)至(6)；(4)将结构面划分成δ-1×δ-1个小网格，采用函数rand()生成随机数，判断所生成随机数的奇偶性，根据随机数的奇偶性选择小网格中三角形的不同划分方案，按照海伦公式计算每个小网格的面积Ai(δ)，其中i＝1,2,…,δ-1×δ-1，进而计算结构面的总面积(5)计算每个小网格的结构面粗糙性概率pi(δ)，(6)计算每个小网格的测度ui(q,δ)，进一步计算log(1/δ)以及(7)按照最小二乘法分别对不同尺度变量δ下的和的数据点对进行拟合，分别得到拟合函数的斜率，就是测度距为q时多重分形谱参数f(a)和a(q)的值；(8)根据测度距q对应的f(a)和a(q)值，计算测度距q对应的广义分形维数Dq，其中时，(9)计算结构面多重分形谱特征参数ΔD(q)和Δa(q)，ΔD(q)＝maxDq-minDq，Δa(q)＝maxa(q)-mina(q)；(10)重复上述(2)至(9)步，直至达到抽样次数要求，将抽样数据置入数据处理软件EXCEL中，利用函数NORMDIST求解累积概率值，通过累积概率值判定抽样结果是否满足3σ法则即抽样结果落在[μ-3σ，μ+3σ]范围内，若满足，则可给出结构面多重分形参数ΔD(q)和Δa(q)概率密度分布函数，并将此抽样结果的均值作为所研究结构面多重分形参数的精准计算值。下面通过具体实施例，对本专利技术做进一步详细说明，应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用于限定本专利技术。如图4所示，本实施例采用石膏人工制作类岩石结构面，其规格长×宽为50mm×50mm，以下结合附图详细说明本专利技术的具体流程：1.结构面三维形貌信息数据获取采用三维扫描仪EinScan-S获取结构面形貌信息数据。采集的数据保存于Excel中，以便下一步计算分析。2.结构面多重分形参数计算按照图1计算步骤，对图4结构面的多重分形参数进行计算，本次共抽样计算了100次。每次计算，参数q的取值为：最小值-45，最大值45。图5是石膏结构面多重分形参数首次抽样计算结果图。结构面多重分形参数100次抽样计算结果汇总表，见表1。多重分形参数值落在某一区间的频数统计结果，见表2。表1石膏结构面多重分形参数100次抽样计算结果汇总表表2多重分形参数值落在某一区间的频数统计结果3.结构面多重分形参数抽样计算分布规律分析由表1数据，经分析计算可知：所计算结构面多重分形参数ΔD(q)的均值μ为0.847641，方差σ为0.000094；所计算结构面多重分形参数Δa(q)的均值μ为0.898478，方差σ为0.000098。由表2数据，可绘制多重分形参数区间-频数直方图；图6、图7分别为石膏结构面多重分形参数ΔD(q)频数分布直方图和石膏结构面多重分形本文档来自技高网...

【技术保护点】
一种岩石结构面多重分形参数的计算方法，其特征在于，该方法包括以下步骤：(1)获取岩石试件结构面的三维形貌数据，包括结构面各点的高度信息，所述高度信息是指该点与结构面内最低点的落差；(2)测度距q分别取值[‑45,45]的整数，重复进行步骤(3)至(8)；(3)尺度变量δ分别取值

【技术特征摘要】
1.一种岩石结构面多重分形参数的计算方法，其特征在于，该方法包括以下步骤：(1)获取岩石试件结构面的三维形貌数据，包括结构面各点的高度信息，所述高度信息是指该点与结构面内最低点的落差；(2)测度距q分别取值[-45,45]的整数，重复进行步骤(3)至(8)；(3)尺度变量δ分别取值重复进行步骤(4)至(6)；(4)将结构面划分成δ-1×δ-1个小网格，采用函数rand()生成随机数，判断所生成随机数的奇偶性，根据随机数的奇偶性选择小网格中三角形的不同划分方案，按照海伦公式计算每个小网格的面积Ai(δ)，其中i＝1,2,…,δ-1×δ-1，进而计算结构面的总面积(5)计算每个小网格的结构面粗糙性概率pi(δ)，(6)计算每个小网格的测度ui(q,δ)，进一步计算log(1/δ)以及(7)按照最小二乘法分别...

【专利技术属性】
技术研发人员：陈世江，王超，郑文翔，
申请(专利权)人：内蒙古科技大学，
类型：发明
国别省市：内蒙古,15