CN110929599A - 基于小波分析的岩体结构面轮廓曲线形态分离方法 - Google Patents

基于小波分析的岩体结构面轮廓曲线形态分离方法 Download PDF

Info

Publication number
CN110929599A
CN110929599A CN201911085984.7A CN201911085984A CN110929599A CN 110929599 A CN110929599 A CN 110929599A CN 201911085984 A CN201911085984 A CN 201911085984A CN 110929599 A CN110929599 A CN 110929599A
Authority
CN
China
Prior art keywords
structural
wavelet
structural surface
rock mass
profile
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
Application number
CN201911085984.7A
Other languages
English (en)
Other versions
CN110929599B (zh
Inventor
雍睿
林杭
谢世杰
杜时贵
刘文连
王秀庆
郑荣跃
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ningbo University
Original Assignee
Ningbo University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ningbo University filed Critical Ningbo University
Priority to CN201911085984.7A priority Critical patent/CN110929599B/zh
Publication of CN110929599A publication Critical patent/CN110929599A/zh
Application granted granted Critical
Publication of CN110929599B publication Critical patent/CN110929599B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/148—Wavelet transforms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

一种基于小波分析的岩体结构面轮廓曲线确定方法，包括以下步骤：(1)采用轮廓曲线仪对结构面表面进行定向测量；(2)选取小波基函数；(3)最大分解尺度J，计算低频细节信号A1‑AJ和高频细节信号D1‑DJ；(4)将AJ作为宏观几何轮廓；(5)将D1‑DJ中每一级别的微观粗糙度高度分布频率与分布频率的信息绘制成图，依据矿物颗粒大小分布规律，将符合高斯分布，统计平均值接近0小波变换所对应的P个D之和为所求微观粗糙度；(6)将D1‑Dp叠加在一起为分离所得微观粗糙度；(7)将原始结构面轮廓曲线除去宏观几何轮廓与微观粗糙度，剩下的就是表面起伏轮廓。本发明适用于定量确定不同类型结构面的表面轮廓形态要素。

Description

基于小波分析的岩体结构面轮廓曲线形态分离方法
技术领域
本发明属于工程技术领域，涉及一种基于小波分析的岩体结构面轮廓曲线分离方法，特别是本发明在岩体结构面表面形态三级划分的基础上，根据不同级别表面形态数值化处理后的频率、幅值的差异，提供了一种基于小波分析的合理分离、重构结构面表面形态的技术方案，解决了结构面起伏表面形态多级定量描述问题，显著提高了结构面表面形态描述的准确性。
背景技术
结构面的存在是岩体介质区别于其他介质的根源，更是影响岩体力学特性的重要因素。结构面对岩体的变形与破坏起着控制作用，而结构面表面形态又是结构面力学性质的重要影响因素，现有技术通常一般按平整度(Unevenness)和起伏度(Waviness)两级形态来描述岩体结构面表面形态模型，名称界线不明确，且没有说明采用二级形态的确定方法，不便于实际应用。因此，杜时贵(2005)将结构面表面形态划分为3级：宏观几何轮廓、表面起伏形态、微观粗糙度。该结构面形态划分模式充分考虑了表面形态的力学机理，形态分级的界线概念明确。本发明主要通过小波分析来实现结构面表面形态三级要素(宏观几何轮廓、表面起伏形态、微观粗糙度)的分离与重构，确定各级表面形态的性质与特征，为岩体结构面表面形态的定量描述提供科学依据，具有重要的工程实用价值。
结构面表面起伏形态不仅包含各种频率信号，还包括几乎没有任何周期或频率特征的突起等局部特征。小波分析在时域和频域上同时具有良好的局部化能力，多分辨率的特性使其可以提取信号任意细节的信息。这一特性恰恰符合了表面形貌微凸体多尺度分布的特点，使其在表面微观形貌重构方面得到了广泛应用。一些代表性的研究如下：
徐磊(2010)截取人工劈裂方法制备的花岗岩结构面二维轮廓线进行几何组成分析，将其分为非稳态倾向、稳态一阶起伏度、稳态二阶粗糙度三个组成部分。并利用小波分析实现了稳态一阶起伏度、稳态二阶粗糙度的分离。
潘学哉(2014)在精确测试岩块试样失稳断裂结构面表面形貌并进行数值化的基础上，定义了能量和信噪比，抽取了岩块试样断裂面不同测段进行研究，结果表明小波变换能判别不同加载模式下的岩石断裂面形貌特征差异特征。
荆君涛(2014)在原有测量数据基础上进行小波分析，将表面粗糙度从原有测量数据中分离出，减少了数据丢失。重构后的微观形貌表面均方根偏差小，为解决表面完整性综合评价提供了一种新的技术途径。
李艳娜(2015)对Mallat算法进行了优化，依据二代提升小波变换理论构建了表面形貌分离的数学模型，通过对模型的仿真实现，成功实现了表面形貌的分离。然而其提取的精确度仍有待进一步深入研究。
以上研究，尚未明确小波分析结构面表面起伏特征的分解尺度的确定方法，更没有提及结构面宏观几何轮廓、表面起伏形态、微观粗糙度的分离与重构的具体方法。由于不同级别表面形态对结构面力学性质的影响存在显著差异，为了研究结构面宏观几何轮廓、表面起伏形态、微观粗糙度对结构面抗剪强度的贡献，迫切需要提出结构面轮廓曲线的定量化分离方法。
发明内容
为了克服传统方法无法量化分析岩体结构面各级表面形态的不足，保证岩体结构面表面形貌特征三级要素分离的精确性，本发明提出一种基于小波分析的岩体结构面轮廓曲线分离的方法，适用于定量确定不同类型结构面的表面轮廓形态要素。
本发明解决其技术问题所采用的技术方案是：
一种基于小波分析的岩体结构面轮廓曲线确定方法，包括以下步骤：
(1)采用轮廓曲线仪对结构面表面进行定向测量，基于图像灰度数据提取技术将结构面测量图像进行坐标数据转化处理；
(2)选取小波基函数，选取的小波基函数必须同时满足对称性、正交性、较短的支撑和较高的消失矩，结合具体结构面表面形貌实际进行确定；
(3)根据结构面表面形貌原始测量数据及工程需要，确定最大分解尺度J，计算每次分解通过小波低通分解滤波器分解出的低频细节信号A1-AJ，计算每次分解通过小波高通分解滤波器分解出的高频细节信号D1-DJ；
(4)将小波变换的最大分解尺度J所对应的AJ作为宏观几何轮廓，它反映了采用峰谷包络线来表征结构面表面起伏形态；
(5)将步骤(4)分离出的D1-DJ中每一级别的微观粗糙度高度分布频率与分布频率的信息绘制成图，依据构成结构面的矿物颗粒大小分布规律，将符合高斯分布，统计平均值接近0小波变换所对应的P个D之和即为所求微观粗糙度，该级别结构面表面起伏形态具体体现了矿物颗粒或细小矿物晶体在结构面表面的分布排列特征；
(6)将步骤(5)确定的D1-Dp叠加在一起即为分离所得微观粗糙度；
(7)将原始结构面轮廓曲线除去步骤(4)、(6)分离出的宏观几何轮廓与微观粗糙度，剩下的就是表面起伏轮廓，至此，结构面表面形态三级要素均已分离得出。
进一步，所述步骤(2)中，选用小波基函数sym5。
本发明的有益效果主要表现在：能够准确的分离出岩体结构面三级表面形态，实现了对三级表面形态的统计、整合、重构，排除了实际表面随机现象对表面特征的影响。本发明为结构面轮廓曲线分离提供了一个科学高效的方法，能够有效解决现有结构面形貌分离中存在的不足与问题。
附图说明
图1为岩体结构面原始几何表面。
图2为基于小波分析的岩体结构面各级二阶粗糙度。
图3为基于小波分析的岩体结构面宏观几何轮廓。
图4为基于小波分析的岩体结构面微观粗糙度。
图5为基于小波分析的岩体结构面表面起伏轮廓。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
下面结合附图对本发明作进一步描述。
参照图1～图5，一种基于小波分析的岩体结构面轮廓曲线分离方法，包括以下步骤：
(1)采用轮廓曲线仪对结构面表面进行定向测量，基于图像灰度数据提取技术将结构面测量图像进行坐标数据转化处理；
(2)选取小波基函数sym5，对结构面轮廓曲线进行二维分离；
(3)根据结构面表面形貌原始测量数据及工程需要，确定最大分解尺度J＝[log2N]＝12，计算每次分解通过小波低通分解滤波器分解出的低频细节信号A1-A12，计算每次分解通过小波高通分解滤波器分解出的高频细节信号D1-D12，见图2；
(4)将小波变换的最大分解尺度12所对应的A12作为宏观几何轮廓，它反映了采用峰谷包络线来表征结构面表面起伏形态，见图3。
(5)将步骤(4)分离出的D1-D12中每一级别的微观粗糙度高度分布频率与分布频率的信息绘制成图，依据构成结构面的矿物颗粒大小分布规律，将符合高斯分布，统计平均值接近0小波变换所对应的D1、D2、D3、D4之和作为微观粗糙度，见图4，该级别结构面表面起伏形态具体体现了矿物颗粒或细小矿物晶体在结构面表面的分布排列特征；
(6)将步骤(5)确定的D1-Dp叠加在一起即为分离所得微观粗糙度；
(7)将原始结构面轮廓曲线除去步骤(4)、(6)分离出的宏观几何轮廓与微观粗糙度，剩下的就是表面起伏轮廓，见图5。至此，结构面表面形态三级要素均已分离得出。
虽然结合附图对发明的具体实施方式进行了详细地描述，但不应理解为对本专利的保护范围的限定。在权利要求书所描述的范围内，本领域技术人员不经创造性劳动即可做出的各种修改和变形仍属本专利的保护范围。

Claims (2)

1.一种基于小波分析的岩体结构面轮廓曲线确定方法，其特征在于，所述方法包括以下步骤：
(1)采用轮廓曲线仪对结构面表面进行定向测量，基于图像灰度数据提取技术将结构面测量图像进行坐标数据转化处理；
(2)选取小波基函数，选取的小波基函数必须同时满足对称性、正交性、较短的支撑和较高的消失矩，结合具体结构面表面形貌实际进行确定；
(3)根据结构面表面形貌原始测量数据及工程需要，确定最大分解尺度J，计算每次分解通过小波低通分解滤波器分解出的低频细节信号A1-AJ，计算每次分解通过小波高通分解滤波器分解出的高频细节信号D1-DJ；
(4)将小波变换的最大分解尺度J所对应的AJ作为宏观几何轮廓，它反映了采用峰谷包络线来表征结构面表面起伏形态；
(5)将步骤(4)分离出的D1-DJ中每一级别的微观粗糙度高度分布频率与分布频率的信息绘制成图，依据构成结构面的矿物颗粒大小分布规律，将符合高斯分布，统计平均值接近0小波变换所对应的P个D之和即为所求微观粗糙度，该级别结构面表面起伏形态具体体现了矿物颗粒或细小矿物晶体在结构面表面的分布排列特征；
(6)将步骤(5)确定的D1-Dp叠加在一起即为分离所得微观粗糙度；
(7)将原始结构面轮廓曲线除去步骤(4)、(6)分离出的宏观几何轮廓与微观粗糙度，剩下的就是表面起伏轮廓，至此，结构面表面形态三级要素均已分离得出。
2.如权利要求1所述的基于小波分析的岩体结构面轮廓曲线确定方法，其特征在于，所述步骤(2)中，选用小波基函数sym5。
CN201911085984.7A 2019-11-08 2019-11-08 基于小波分析的岩体结构面轮廓曲线形态分离方法 Active CN110929599B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911085984.7A CN110929599B (zh) 2019-11-08 2019-11-08 基于小波分析的岩体结构面轮廓曲线形态分离方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911085984.7A CN110929599B (zh) 2019-11-08 2019-11-08 基于小波分析的岩体结构面轮廓曲线形态分离方法

Publications (2)

Publication Number Publication Date
CN110929599A true CN110929599A (zh) 2020-03-27
CN110929599B CN110929599B (zh) 2023-05-26

Family

ID=69852464

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911085984.7A Active CN110929599B (zh) 2019-11-08 2019-11-08 基于小波分析的岩体结构面轮廓曲线形态分离方法

Country Status (1)

Country Link
CN (1) CN110929599B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113049403A (zh) * 2021-03-02 2021-06-29 宁波大学 考虑形貌频谱特征的结构面法向卸荷剪切损伤试验方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103149348A (zh) * 2013-03-05 2013-06-12 天津大学 基于小波多尺度分析的先进陶瓷磨削表面损伤的评价方法
US20150347854A1 (en) * 2014-04-25 2015-12-03 Huntington Ingalls Incorporated System and Method for Using Augmented Reality Display in Surface Treatment Procedures
CN110227968A (zh) * 2019-05-16 2019-09-13 中国工程物理研究院激光聚变研究中心 基于小波变换的磁流变去除函数寄生条纹消除方法及装置
CN110288700A (zh) * 2019-06-26 2019-09-27 东北大学 一种岩质边坡结构面自动分组及位移预测方法
CN110362955A (zh) * 2019-07-25 2019-10-22 四川大学 岩质高边坡稳定性分析三维地质力学模型试验方法及应用

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103149348A (zh) * 2013-03-05 2013-06-12 天津大学 基于小波多尺度分析的先进陶瓷磨削表面损伤的评价方法
US20150347854A1 (en) * 2014-04-25 2015-12-03 Huntington Ingalls Incorporated System and Method for Using Augmented Reality Display in Surface Treatment Procedures
CN110227968A (zh) * 2019-05-16 2019-09-13 中国工程物理研究院激光聚变研究中心 基于小波变换的磁流变去除函数寄生条纹消除方法及装置
CN110288700A (zh) * 2019-06-26 2019-09-27 东北大学 一种岩质边坡结构面自动分组及位移预测方法
CN110362955A (zh) * 2019-07-25 2019-10-22 四川大学 岩质高边坡稳定性分析三维地质力学模型试验方法及应用

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HANG LIN.ET.: "An empirical statistical constitutive relationship for rock joint shearing considering scale effect" *
侯钦宽等: "结构面粗糙度统计测量最小样本数确定方法" *
孙盛玥等: "天然岩石节理双阶粗糙度分形特征研究" *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113049403A (zh) * 2021-03-02 2021-06-29 宁波大学 考虑形貌频谱特征的结构面法向卸荷剪切损伤试验方法

Also Published As

Publication number Publication date
CN110929599B (zh) 2023-05-26

Similar Documents

Publication Publication Date Title
CN102944174B (zh) 一种三维激光点云数据的去噪与精简方法及系统
Zahouani et al. Fractal models of surface topography and contact mechanics
Du et al. A shearlet-based separation method of 3D engineering surface using high definition metrology
CN110807781B (zh) 一种保留细节与边界特征的点云精简方法
CN112150430B (zh) 一种利用岩石细观结构数字图像的数值分析方法
Wang Bilateral recovering of sharp edges on feature-insensitive sampled meshes
CN101082484B (zh) 一种二维表面粗糙度评定中建立轮廓基准线的方法
CN100461217C (zh) 一种复杂性测度的图像纹理分割方法
CN103149348B (zh) 基于小波多尺度分析的先进陶瓷磨削表面损伤的评价方法
CN104268940B (zh) 基于ct扫描图像的mems结构重构与检测方法
CN110929599A (zh) 基于小波分析的岩体结构面轮廓曲线形态分离方法
CN110908001A (zh) 一种大地电磁测深信号的重构方法及系统
CN114707555A (zh) 崩岸前兆次声波信号降噪处理方法、特征识别方法及装置
Li et al. Robust estimation in Gaussian filtering for engineering surface characterization
CN113049403B (zh) 考虑形貌频谱特征的结构面法向卸荷剪切损伤试验方法
CN109671015A (zh) 基于dem分形特征的坡度尺度变换方法
Stachowiak et al. Classification of tribological surfaces
CN117173345A (zh) 基于分形理论和蒙特卡罗方法的粗糙表面形貌模拟方法
OUYANG et al. Edge detection in pavement crack image with beamlet transform
Ayala-Landeros et al. Correlation between roughness (Ra) and fractal dimension (D) using artificial vision systems for on-site inspection
CN107239629B (zh) 一种岩石结构面实验室合理尺寸确定的分形维数分析方法
CN113049326B (zh) 结构面形貌频谱成分定量可控的制样方法
Podsiadlo et al. Hybrid fractal-wavelet method for characterization of tribological surfaces—a preliminary study
CN106682237A (zh) 一种山体边界自动提取方法
Podulka Application of image processing methods for the characterization of selected features and wear analysis in surface topography measurements

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