CN103884312B - 一种岩石结构面形貌尺寸效应分析方法 - Google Patents
一种岩石结构面形貌尺寸效应分析方法 Download PDFInfo
- Publication number
- CN103884312B CN103884312B CN201410141917.3A CN201410141917A CN103884312B CN 103884312 B CN103884312 B CN 103884312B CN 201410141917 A CN201410141917 A CN 201410141917A CN 103884312 B CN103884312 B CN 103884312B
- Authority
- CN
- China
- Prior art keywords
- structural plane
- structural
- face
- value
- 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.)
- Expired - Fee Related
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提供一种岩石结构面形貌尺寸效应测量分析方法。其步骤包括:⑴首先获取岩石结构面三维形貌数据;⑵按照实验变异函数公式计算结构面分析方向上的点对,并绘制随变化的发展趋势图;⑶确定变异函数模型,应用最小二乘法进行模型拟合,进而确定模型参数变程a与基台C;⑷应用量纲分析法,用<i>JRCv</i>(由a和C表示)表征结构面分析方向上结构面粗糙性特征;⑸按照步骤⑵、⑶和⑷计算结构面不同尺寸范围的各个方向上的<i>JRCv</i>值,以此分析结构面的尺寸效应特征。本发明提供的方法可以定量化表征结构面尺寸效应特征,分析结果为实验室或原位结构面力学试验选择合理尺寸试件提供依据。
Description
技术领域
本发明涉及一种岩石结构面形貌尺寸效应分析方法,属于岩石力学领域范畴。
背景技术
大量的研究表明,岩石结构面存在明显的尺寸效应特征。早在1977年Barton和Choubey研究岩石结构面剪切强度时,就对结构面尺寸效应进行过相关论述。之后,Barton和Bandi给出了结构面粗糙度系数JRC的尺寸效应修正公式。近年来,国内外诸多学者采用分形理论中的分形维数D和截距A研究了节理面的尺寸效应,得出了参数D和A随着结构面尺寸范围的增大而减小并趋于恒定的结论。尽管国内外学者做了大量研究,但因岩石结构面的复杂性,其尺寸效应特征定量化表征仍是一个难点。
发明内容
本发明的目的在于提供一种能够定量化描述岩石结构面尺寸效应特征的分析方法,从而进一步为实验室或原位结构面力学试验提供一种选择合理尺寸试件的依据。
技术解决方案:
本发明首先获取岩石结构面三维形貌数据,然后应用地质统计学原理,计算结构面各个方向上的变异函数,通过分析变异函数特征,提出用变异函数参数(基台和变程)定量化描述岩石结构面粗糙度(JRCv)的方法,最后应用该法分析岩石结构面的尺寸效应特征。
用该方法分析岩石结构面尺寸效应特征的具体步骤如下:
(1)获取岩石结构面数字图像,用图像灰度值表示岩石结构面的三维形貌特征;
(2)以岩石结构面的三维数据为基础,按照实验变异函数计算公式,计算结构面分析方向上的点对,并绘制随变化的发展趋势图,其中,为变异函数值,代表结构面点处的起伏高度,表示结构面内点的编号,表示分析方向上两点的间距,表示分析方向上结构面内间隔为的点对数;
(3)根据变异函数值随变化的发展趋势,确定变异函数模型,应用最小二乘法进行变异函数模型拟合,确定模型参数变程a与基台C的值;
(4)用JRCv表示结构面的粗糙度,,变程a表示分析方向结构面起伏的频率,而基台C表示结构面起伏度的大小;
(5)划分结构面尺寸范围,不同尺寸范围的结构面各个方向的粗糙度JRCv重复步骤(2)—(4)获得。
本发明步骤3)中变异函数模型包括:球状模型、或对数模型、或指数模型、或高斯模型。
本发明还包括用图像灰度值表示岩石结构面三维形貌特征,具体步骤为:
①将获得的结构面图像转化为256灰度图;
②在灰度图中0代表最暗,255代表最亮,在岩石结构面中灰度最小值代表最低点,灰度最大值表示最高点,按照公式对结构面中的灰度值进行比例缩放,并使其单位转化为像素,从而使缩放后的灰度值表示结构面的实际尺寸,其中,表示结构面图像长度,单位像素,表示结构面实际长度,单位mm,表示图像中最大灰度值,表示图像中最小灰度值,表示结构面内最高点,单位mm,表示结构面内最低点,单位mm,为缩放比例系数。
本发明还包括结构面的粗糙度JRCv定量化表达。变程很好的反映了变量的影响范围,相对于所研究结构面区域,在分析方向变程较大,说明沿该方向结构面起伏的频率相对较小,而基台反映了结构面起伏度大小的变化,值越大,结构面起伏度越大,因此,结构面的粗糙度与变程成反比,与基台成正比,变程的量纲为长度的二次方,基台的量纲为长度的一次方,因此与变程,基台的关系可表示成。
本发明优点:
由于岩石结构面形貌的复杂性,其尺寸效应的定量化描述仍是该领域的难点。本发明通过分析地质统计学中变异函数参数(基台和变程)的性质特征,采用量纲分析方法,提出用JRCv(由基台和变程表示)表征结构面的粗糙度。通过计算岩石结构面不同尺寸范围的各个方向上的JRCv值,来分析岩石结构面的尺寸效应特征。分析结果为实验室或原位结构面力学试验提供了一种选择合理尺寸试件的方法,从而进一步为准确全面评价岩石结构面剪切强度提供依据。
附图说明
图1为本发明岩石结构面形貌尺寸效应分析流程图;
图2为本发明石膏结构面形貌图;
图3为本发明图像灰度值表示岩石结构面三维形貌特征流程图;
图4为本发明结构面图像数据点分布方式及分析方向示意图;
图5为本发明结构面27°方向上点对统计示意图;
图6为本发明结构面上爬坡数据示意图;
图7为本发明27°方向上二级球状模型变异函数拟合及粗糙度参数JRCv计算结果图;
图8为本发明结构面尺寸范围划分图;
图9为本发明结构面尺寸效应特征图;
图10为本发明结构面方向指示图。
具体实施方式
实施例1
本实施例采用石膏粉材料人工制作类岩石结构面,以下结合附图详细说明本发明的具体实施方式。
见图1是用变异函数参数性质分析岩石结构面形貌尺寸效应流程图。
)结构面图像获取与用灰度值表示结构面三维形貌
采用佳能60D数码相机获取的石膏结构面形貌图(见图2),将该结构面形貌图像进行图像灰度转换,比例缩放,而后将处理后的数据进行存储。图像灰度值表示岩石结构面三维形貌特征流程图见图3。该过程中比例缩放确定过程具体步骤如下:
(1)经实际量测,该结构面最高点与最低点落差()为9mm;
(2)从结构面图像灰度值统计结果看,该结构面图像最大灰度值208,最小灰度值36,最大灰度差()为172;
(3)该结构面实际长度为50mm,结构面图像实际长度为1000像素;
(4)按照公式计算,比例系数为45/43。
)结构面分析方向上变异函数模型参数及粗糙度JRCv计算
按照图1步骤(4)首先确定结构面分析方向,如图4所示,首先在结构面图像像素点矩阵中建立平面直角坐标系;其次在直角坐标系中选择两个像素点,如像素点1和像素点2;最后根据1、2两点在直角坐标系中的位置计算方向角度,即,,;这样图4中所有斜直线方向的点将参与计算分析。
按照实验变异函数计算公式,计算该方向上的点对。见图5结构面27°方向上点对统计示意图,当为1时,在该方向上逐行追踪扫描图像中的每一点,并记录符合=1的点的值,将参与两点间高程差的计算,这样的点对在示例中共有20对,同样方法,逐步增大,相应地统计出点的对数,这样能够计算不同的对应的的值,进而绘制随变化的发展趋势图(见图7所示)。该过程主要为了说明点对的统计,并未考虑在该方向上各点的分布特点,为了区分同一方位正反方向的各向异性特征,计算时采用沿分析方向上所有爬坡数据参与变异函数的计算,爬坡数据示意图见图6,此过程可以通过比较相邻点数值的大小,分析其是否参与计算。
按照图1步骤(5)确定变异函数模型,根据图7中随变化的发展趋势,前半段数据点较陡,后半段数据点相对较缓,该特征符合二级球状模型变异函数曲线特点,故采取二级球状模型变异函数进行拟合,拟合结果见图7拟合曲线,相应变异函数模型参数见图7指示部分。二级球状模型变异函数拟合具体步骤如下:
下述公式为二级球状模型变异函数模型计算公式。
⑴先用前区数据点在拟合第一段球状模型函数,设,,,,,,得到,按照最小二乘法的原理和方法,得到,,的解。
⑵用后区的数据点在内拟合第二段球状模型,令,,,
,,,这样得到,同理按照最小二乘法的原理和方法,得到,,的解。
⑶进而能求得参数,,,,的值。
二级球状模型变异函数参数变程C=,基台a=,进而JRVc由获得。
按照上述方法,本实施例共计算了10个尺寸范围(见图8)(从100像素×100像素到1000像素×1000像素)32个方向的JRVc值,计算结果见表1
表1不同尺度范围的结构面上不同方向的JRCv值统计表
3)结构面形貌尺寸效应分析
见图9是根据表1数据绘制的结构面尺寸效应特征图,关于图中所示方向与图10对应,从图9中可以看出,结构面的粗糙度具有明显的尺寸效应,结构面尺寸较小时,JRCv值较大;随着尺寸的增大,JRCv值逐渐地减小;达到一定尺寸时(此尺度成为临界尺度),JRCv值趋于稳定。
由实施例结果可见,提出JRCv能够很好地表示岩石结构面的尺寸效应,分析结果为实验室或原位结构面力学试验提供了一种选择合理尺寸试件的手段,在准确全面评价结构面剪切强度实践中具有重要的理论意义与实用价值。
Claims (4)
1.一种岩石结构面形貌尺寸效应评价分析方法,其特征在于,方法步骤包括:
(1)获取岩石结构面数字图像,用图像灰度值表示岩石结构面的三维形貌特征;
(2)以岩石结构面的三维数据为基础,按照实验变异函数计算公式,计算结构面分析方向上的点对,并绘制随h变化的发展趋势图,其中,为变异函数值,Z(x)代表结构面x点处的起伏高度,i表示结构面内点的编号,h表示分析方向上两点的间距,N(h)表示分析方向上结构面内间隔为h的点对数;
(3)根据变异函数值随h变化的发展趋势,确定变异函数模型,应用最小二乘法进行变异函数模型拟合,确定模型参数变程a与基台C的值;
(4)用JRCv表示结构面的粗糙度,,变程a表示分析方向结构面起伏的频率,而基台C表示结构面起伏度的大小;
(5)划分结构面尺寸范围,不同尺寸范围的结构面各个方向的粗糙度JRCv重复步骤(2)—(4)获得。
2.根据权利要求1所述一种岩石结构面形貌尺寸效应评价分析方法,其特征在于,还包括用图像灰度值表示岩石结构面的三维形貌特征,具体步骤为:
①将获得的结构面图像转化为256灰度图;
②在灰度图中0代表最暗,255代表最亮,在岩石结构面中灰度最小值代表最低点,灰度最大值表示最高点,按照公式对结构面中的灰度值进行比例缩放,并使其单位转化为像素,从而使缩放后的灰度值表示结构面的实际尺寸,其中,Limage表示结构面图像长度,单位像素,L表示结构面实际长度,单位mm,hhmax表示图像中最大灰度值,hhmin表示图像中最小灰度值,hmax表示结构面内最高点,单位mm,hmin表示结构面内最低点,单位mm,k为缩放比例系数。
3.根据权利要求1所述一种岩石结构面形貌尺寸效应评价分析方法,其特征在于,还包括:结构面的粗糙度JRCv定量化表达,变程a反映了变量的影响范围,相对于所研究结构面区域,在分析方向变程a越大,说明沿该方向结构面起伏的频率越小,而基台C反映了结构面起伏度大小的变化,C值越大,结构面起伏度越大,因此,结构面的粗糙度JRCv与变程a成反比,与基台C成正比,变程a的量纲为长度的二次方,基台C的量纲为长度的一次方,因此JRCv与变程a,基台C的关系表示成。
4.根据权利要求1所述一种岩石结构面形貌尺寸效应评价分析方法,其特征在于,步骤3)中变异函数模型包括:球状模型、或对数模型、或指数模型、或高斯模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410141917.3A CN103884312B (zh) | 2014-04-10 | 2014-04-10 | 一种岩石结构面形貌尺寸效应分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410141917.3A CN103884312B (zh) | 2014-04-10 | 2014-04-10 | 一种岩石结构面形貌尺寸效应分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103884312A CN103884312A (zh) | 2014-06-25 |
CN103884312B true CN103884312B (zh) | 2016-05-18 |
Family
ID=50953337
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410141917.3A Expired - Fee Related CN103884312B (zh) | 2014-04-10 | 2014-04-10 | 一种岩石结构面形貌尺寸效应分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103884312B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105466790B (zh) * | 2015-11-10 | 2017-12-29 | 内蒙古科技大学 | 一种各向异性特征的岩石结构面剪切强度评估方法 |
CN105809692B (zh) * | 2016-03-10 | 2017-05-03 | 中国石油大学(华东) | 一种页岩结构的定量表征方法 |
CN107144243A (zh) * | 2017-04-19 | 2017-09-08 | 中国电力科学研究院 | 一种测量岩体结构面剪切破坏面积的方法和系统 |
CN107292017B (zh) * | 2017-06-15 | 2020-10-20 | 内蒙古科技大学 | 岩石结构面实验室合理尺寸确定的多重分形参数分析方法 |
CN109099880B (zh) * | 2018-07-12 | 2020-08-11 | 宁波大学 | 岩体结构面粗糙度系数全域搜索测量方法 |
CN109543280B (zh) * | 2018-11-19 | 2021-11-30 | 东南大学 | 一种基于形态修正的节理面粗糙度计算方法 |
CN112903431B (zh) * | 2021-01-18 | 2022-10-21 | 绍兴文理学院 | 节理岩体强度尺寸效应统计方法 |
CN112903434B (zh) * | 2021-01-19 | 2022-10-21 | 绍兴文理学院 | 节理岩体强度尺寸效应统计模型构建方法 |
CN112903433B (zh) * | 2021-01-19 | 2022-10-21 | 绍兴文理学院 | 节理岩体强度尺寸效应规律量化方法 |
CN113049403B (zh) * | 2021-03-02 | 2022-11-29 | 宁波大学 | 考虑形貌频谱特征的结构面法向卸荷剪切损伤试验方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1645049A (zh) * | 2004-12-15 | 2005-07-27 | 金华职业技术学院 | 典型岩体结构面粗糙度系数尺寸效应分维数确定方法 |
CN1645051A (zh) * | 2004-12-15 | 2005-07-27 | 金华职业技术学院 | 岩体结构面粗糙度系数尺寸效应有效长度确定方法 |
CN101038158A (zh) * | 2007-04-05 | 2007-09-19 | 同济大学 | 大尺寸岩石表面三维形貌的测量方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3704220B2 (ja) * | 1996-02-14 | 2005-10-12 | 増之 氏平 | 落石源の位置評定システム |
JP3892927B2 (ja) * | 1997-01-28 | 2007-03-14 | 太陽工業株式会社 | 擬岩パネルによる擬岩地形の設計方法及び設計支援システム |
-
2014
- 2014-04-10 CN CN201410141917.3A patent/CN103884312B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1645049A (zh) * | 2004-12-15 | 2005-07-27 | 金华职业技术学院 | 典型岩体结构面粗糙度系数尺寸效应分维数确定方法 |
CN1645051A (zh) * | 2004-12-15 | 2005-07-27 | 金华职业技术学院 | 岩体结构面粗糙度系数尺寸效应有效长度确定方法 |
CN101038158A (zh) * | 2007-04-05 | 2007-09-19 | 同济大学 | 大尺寸岩石表面三维形貌的测量方法 |
Non-Patent Citations (2)
Title |
---|
基于修正线粗糙度法的岩石节理粗糙度估值;陈世江等;《金属矿山》;20120630(第6期);第22-25页 * |
岩体结构面三维表面形貌的尺寸效应研究;徐磊等;《武汉理工大学学报》;20080430;第30卷(第4期);第103-105页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103884312A (zh) | 2014-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103884312B (zh) | 一种岩石结构面形貌尺寸效应分析方法 | |
CN103926156B (zh) | 一种三维岩石结构面剪切力学行为多重分形评估方法 | |
CN103886613A (zh) | 一种岩石结构面形貌各向异性评价方法 | |
Galantucci et al. | Advanced damage detection techniques in historical buildings using digital photogrammetry and 3D surface anlysis | |
Polidori et al. | Description of terrain as a fractal surface, and application to digital elevation model quality assessment | |
CN106248003B (zh) | 一种三维激光点云提取植被冠层聚集度指数的方法 | |
Tang et al. | Study on estimation method of rock mass discontinuity shear strength based on three-dimensional laser scanning and image technique | |
CN102841385B (zh) | 一种基于多重分形克里金法的局部地磁图构建方法 | |
CN105466790A (zh) | 一种各向异性特征的岩石结构面剪切强度评估方法 | |
CN108759699A (zh) | 一种大视场砌体结构材料三维全场变形的测量方法及系统 | |
CN105761310B (zh) | 一种天空可视域数字地图的模拟分析及图像显示方法 | |
Liu et al. | The detection effect of pavement 3D texture morphology using improved binocular reconstruction algorithm with laser line constraint | |
CN101178811A (zh) | 一种三步法结构光直线光条的图像特征提取方法 | |
Lai et al. | Surface roughness of rock faces through the curvature of triangulated meshes | |
CN107452059B (zh) | 一种三维插值建模方法、设备及其存储设备 | |
CN103268423A (zh) | 一种地理现象多点模拟空间尺度选择的方法 | |
CN102750729A (zh) | 一种使用近景普通数码相机的快速单木建模方法 | |
CN109493380A (zh) | 一种岩石节理面抗剪试验中不规则剪切面面积的计算方法 | |
Salah | Geostatistical analysis of groundwater levels in the south Al Jabal Al Akhdar area using GIS | |
Lorenz | Fractal geometry of architecture: Implementation of the box-counting method in a CAD-software | |
Alghalandis | Stochastic Modelling of Fractures | |
CN103606150B (zh) | 一种局部方差定量探测规则地物空间格局大小的方法 | |
CN107784171B (zh) | 一种空间模型构建方法 | |
KR100924929B1 (ko) | 유한분할을 통한 건물의 일조면적 분석방법 | |
ONIGA et al. | TERRESTRIAL LASER SCANNER SURVEYING VERSUS TOTAL STATION SURVEYING FOR 3D BUILDING MODEL GENERATION. |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160518 Termination date: 20170410 |