CN113188488A - 锚杆长度无损测量方法 - Google Patents
锚杆长度无损测量方法 Download PDFInfo
- Publication number
- CN113188488A CN113188488A CN202110473896.5A CN202110473896A CN113188488A CN 113188488 A CN113188488 A CN 113188488A CN 202110473896 A CN202110473896 A CN 202110473896A CN 113188488 A CN113188488 A CN 113188488A
- Authority
- CN
- China
- Prior art keywords
- data
- length
- reflected wave
- wave
- anchor rod
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B15/00—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明涉及一种锚杆长度无损测量方法,用以解决现有技术中存在的最锚杆无损检测数据处理分析时,无法准确获得长度信息、滤波器不稳定和技术员主观影响大的问题。本发明实施例首先对原始数据进行数据叠加,随后进行降频处理、滤波,如果滤波阶数n<20则恢复原采集频率,如果滤波阶数n大于20则调整下采集参数。然后截取反射波有效信息对整体波形数据进行滑动互相关数据处理得到相似系数波形。最后根据相似系数波型图确定反射波初始位置,计算锚杆长度。本发明可以有效解决无法准确获得长度信息、滤波器不稳定和人员对处理结果主观影响的问题。
Description
技术领域
本发明属于工程测量领域,尤其涉及一种锚杆长度无损测量方法。
背景技术
锚杆支护是通过围岩内部的锚杆改变围岩本身的力学状态,在巷道周围形成一个整体而又稳定的岩石带,利用锚杆与围岩共同作用,达到维护巷道稳定的目的。锚杆的力学作用主要有悬吊作用 、组合梁作用 、 组合拱作用、减跨作用 、加固作用。
锚杆不但支护效果好,且用料省、施工简单、有利于机械化操作、施工速度快。但是锚杆长度达不到设计长度,锚固密实度差会造成锚杆支护巷道失稳,形成安全隐患。
常见的锚杆检测技术有拉拔承载力试验和锚杆无损检测试验。锚杆无损检测试验具有非破坏性和高效的优点,但同时也存在如下技术问题:
1.当锚杆较短时,直达波与反射波相互干涉,无法获得准确的长度信息;
2.为了获取分辨率较高的波形图,检测人员往往采用较高的采样频率,而锚杆往往较短,所以当采样频率较高,样本点数较少时会导致巴特沃斯滤波器滤波阶数过高,当滤波阶数n大>20时则会造成滤波器不稳定;
3.锚杆无损检测获取的波形数据处理与分析主观影响大,需要技术员的强大知识储备与丰富经验。根据《锚杆锚固质量无损检测技术规程》JGJ/T 182和《水利水电工程锚杆无损检测规程》DL/T 5424均对波形的判定没有明确规定,执行力较差,检测人员很难根据规程规范的要求对锚杆长度进行判断。尤其当锚杆没有严格按照规范施工时,锚杆锚固密实度差,会产生大量绕射波和反射波形成干扰波,无法区分。
发明内容
本发明的目的是提供一种锚杆长度无损测量方法,以解决现有技术中锚杆长度无法准确获取的问题。该方法步骤如下:
锚杆长度无损测量方法,包括:
S100,基于应力波法获取多道锚杆检测的原始数据;
S200,对各道原始数据进行叠加;
S300,对叠加数据降频处理后进行巴特沃斯滤波;
S400,将滤波后数据的频率恢复为原始数据的频率;
S500,截取S400处理后数据中反射波的有效信息;
S600,输入反射波有效信息对S400处理后数据的叠加波形进行滑动互相关处理,并获取相似系数;
S700,利用相似系数波形图获取反射波初始位置;
S800,计算锚杆长度。
优选地,在步骤S100和S200之间,还包括:
对原始数据进行相位反转、初至校准和错误数据剔除。
优选地,步骤S300中还具有:
降频处理后求取滤波阶数n,如果n>20则需要调整下采集参数;如果n<20,则进行所述巴特沃斯滤波。
优选地,步骤S500的截取方法包括:
截取数据处理后可以分辨的完整反射波,并截取反射波波形前段的两到三个波峰或波谷。
优选地,步骤S600的利用滑动互相关算法对反射波形进行扫描的方法包括:
利用截取的反射波信息从整段波形的第一个采样点开始计算等个数采样点与反射波采样点的互相关值。
优选地,S700中反射波初始位置的获取方法包括:
基于相似系数波形图出现的多个波峰,第一个波峰处相似系数值为1,选取第一个波峰后且相似系数最接近1的点,其相对应的采样点即为反射波初始的位置,即可根据该位置计算反射波初始的时间。
优选地,S800中锚杆长度计算方法包括:
-利用室内模拟试验或现场试验确定杆体波速和杆系波速;
-确定 杆端外露长度;
-利用杆体速度、杆系速度、杆端外露长度、初至时间和反射波初试时间计算锚杆长度。
本发明的有益效果:
本发明对原始数据依次进行叠加、滤波、恢复频率、截取、滑动互相关处理的方式,获取反射波初始位置,以计算锚杆长度,可以有效解决无法准确获得长度信息、滤波器不稳定和人员对处理结果主观影响的问题。
附图说明
图1是本发明锚杆长度无损测量方法的流程示意图(在单次运行中右列位于左列下方);
图2是本发明锚杆长度无损测量方法的6道原始数据;
图3是本发明锚杆长度无损测量方法的未处理直接叠加后的波形图;
图4是本发明锚杆长度无损测量方法的经过相位反转、和初至校准后的叠加波形图及频谱特征;
图5是本发明锚杆长度无损测量方法的经过降频处理后的波形图及频谱特征;
图6是本发明锚杆长度无损测量方法的经滤波处理后的波形图及频谱特征;
图7是本发明锚杆长度无损测量方法的恢复原有采集频率后的波形图;
图8是本发明锚杆长度无损测量方法的相似系数波形图。
具体实施方式
本实施例中锚杆长度无损测量方法是由数据处理方法和数据分析方法两部分组成的。数据处理方法旨在于对应力波锚杆无损检测数据进行除噪,一方面剔除错误数据,另一方面为数据分析提供更多的高质量数据样本。数据分析方法旨在依据处理后数据进行量化分析,更加准确的计算出锚杆长度。
数据处理方法,包括:
-通过应力波锚杆无损检测获取原始数据;
-将多道原始数据进行线性叠加;
-将叠加后的数据进行傅里叶变化,获取频谱;
-将叠加后的数据进行降频处理;
-将降频后数据通过巴特沃斯滤波器进行滤波;
-恢复原有数据采集频率。
上述数据处理方法,由于数据采集时获取的多道数据存在相位相反、初至时间不同、信息错误等问题,所以在叠加前应排除原始数据中位相反、初至时间不同、信息错误等问题。为此,对通过应力波锚杆无损检测获取原始数据进行数据处理,包括:
-对相位相反的原始数据进行相位翻转;
-对初至时间不同的原始数据在时间轴进行移动;
-对信息错误的数据进行剔除。
上述数据处理方法,当锚杆较短时,由于采样频率/采样点数比值过大,导致巴特沃斯滤波器阶数n过大,导致滤波器不稳定,所以在滤波前应排除滤波器不稳定的问题。为此,对通过对叠加后的数据进行数据处理,包括:
-将叠加后数据进行巴特沃斯滤波器处理,计算叠加后数据滤波阶数n;
-当滤波阶数n大于20时,对叠加后数据进行降频处理(控制滤波阶数小于20即可,滤波阶数大于20滤波器容易崩溃);
-对降频后数据进行巴特沃斯滤波器处理,计算叠加后数据滤波阶数n,且滤波阶数n小于20。
上述数据处理方法,由于叠加数据进行降频处理,滤波后会导致数据分辨率降低,所以应排除数据分辨率较低的问题。为此,对降频处理的叠加数据进行数据处理,包括:
-对滤波后数据进行恢复原有采集频率。
数据分析方法,该数据分析程序是在matlab中实现,调整参数就是改变程序中的数值,包括:
-对原始数据进行数据处理;
-截取发射波信息;
-利用滑动互相关算法对整段波形进行扫描;
-获取相似系数波形图;
-确定反射波位置;
-计算锚杆长度。
上述数据分析方法,由于锚杆较短直达波和反射波相互干涉,或锚杆前段灌浆不密实,不密实区域绕射波或反射波和直达波相互干涉,造成直达波难以完整区分,所以应排除杆底反射波或不密实处绕射波或反射波进入截取的波段。为此,对反射波截取的数据分析包括:
-截取数据处理后可以分辨的完整直达波(通过excel数据截取例如第25~90个采样点的数据);
-截取数据处理后波形前段的两到三个波峰或波谷。
利用滑动互相关算法对整段波形进行扫描,包括:
-提取截取的直达波信息;
-利用截取的直达波信息从整段波形的第一个采样点开始计算等个数采样点与直达波采样点的互相关值,即为相似系数i0;
-逐一计算相似系数,获得相似系数i0、i1、i2、。。。in;
-绘制相似系数波形图。
上述数据分析方法,当不能截取完整直达波时,由于仅截取部分直达波信息,相似系数波形图会出现多个波峰,且第一个波峰处相似系数值为1,应选取第一个波峰后且相似系数最接近1的点,其相对应的采样点即为反射波初始的位置,即可根据相应位置计算反射波初始的时间。
上述数据分析方法,其中锚杆长度的计算方法是:
-利用室内模拟试验或现场试验确定杆体波速和杆系波速;
-确定杆端外露长度;
-利用杆体速度、杆系速度、杆端外露长度、初至时间和反射波初试时间计算锚杆长度。
依据上述数据处理方法和数据分析方法,可总结出如下锚杆无损检测数据处理分析方法,如图1所示,具体包括以下步骤:
-采用应力波锚杆无损检测获取多道原始数据。
-对多道原始数据进行叠加;
-对叠加数据进行降频处理并进行巴特沃斯滤波;
-滤波后恢复原有采集频率;
-对数据处理后的数据截取反射波有效信息;
-输入反射波有效信息对叠加波形进行滑动互相关处理,并获取相似系数;
-利用相似系数波形图获取反射波初始位置;
-计算锚杆长度。
这里需要说明的是,在多道原始数据叠加前应,应对相位相反、初至位置不同、信息错误的原始数据进行相位反转、初至校准和错误数据提出的处理。
在一种具体的实施方式中,如图2所示,第3道数据初至相对其他道有延迟,第6道数据相位相对其他道相反,不经处理直接叠加后如图3所示,对第3道数据初至校准后,第6道数据相位反转后叠加数据如图4所示,处理后数据比不经处理叠加后叠加数据波形明显更加清晰。
这里需要说明的是,在进行巴特沃斯数字滤波时,应注意滤波阶数n应小于20,当滤波阶数n大于20时,应重新调整下采样参数,时滤波阶数n小于20。
在一种具体的实施方式中,如图7所示,经巴特沃斯数字滤波后波形更加清晰,高频信号得到有效压制。
这里需要说明的是,相似系数波形图中会出现多个波峰,且第一个波峰相似系数值为1,应选择第一个波峰后相似系数最接近1的采样点为反射波起始位置。
在一种具体实施方式中,如图8所示,相似系数波形图第一个波峰值为1,且第四个波峰值为0.98最接近1,该点为第653个采样点,即第653个采样点为反射波初始位置。利用杆体波速、杆系波速、锚杆外露长度、采样间隔和采样点位置即可计算锚杆长度。
以上仅就本发明应用较佳的实例做出了说明,但不能理解为是对权利要求的限制,本发明的结构可以有其他变化,不局限于上述结构。总之,凡在本发明的独立权利要求的保护范围内所作的各种变化均在发明的保护范围内。
Claims (7)
1.锚杆长度无损测量方法,其特征在于,包括:
S100,基于应力波法获取多道锚杆检测的原始数据;
S200,对各道原始数据进行叠加;
S300,对叠加数据降频处理后进行巴特沃斯滤波;
S400,将滤波后数据的频率恢复为原始数据的频率;
S500,截取S400处理后数据中直达波的有效信息;
S600,输入反射波有效信息对S400处理后数据的叠加波形进行滑动互相关处理,并获取相似系数;
S700,利用相似系数波形图获取反射波初始位置;
S800,计算锚杆长度。
2.根据权利要求1所述的锚杆长度无损测量方法,其特征在于,在步骤S100和S200之间,还包括:
对原始数据进行相位反转、初至校准和错误数据剔除。
3.根据权利要求1所述的锚杆长度无损测量方法,其特征在于,在步骤S300中还具有:
降频处理后求取滤波阶数n,如果n>20则需要调整下采集参数;如果n<20,则进行所述巴特沃斯滤波。
4.根据权利要求1所述的锚杆长度无损测量方法,其特征在于,步骤S500的截取方法包括:
截取数据处理后可以分辨的完整反射波,并截取反射波波形前段的两到三个波峰或波谷。
5.根据权利要求1所述的锚杆长度无损测量方法,其特征在于,步骤S600的利用滑动互相关算法对反射波形进行扫描的方法包括:
利用截取的反射波信息从整段波形的第一个采样点开始计算等个数采样点与反射波采样点的互相关值。
6.根据权利要求1所述的锚杆长度无损测量方法,其特征在于,S700中反射波初始位置的获取方法包括:
基于相似系数波形图出现的多个波峰,第一个波峰处相似系数值为1,选取第一个波峰后且相似系数最接近1的点,其相对应的采样点即为反射波初始的位置,即可根据该位置计算反射波初始的时间。
7.根据权利要求1所述的锚杆长度无损测量方法,其特征在于,S800中锚杆长度计算方法包括:
-利用室内模拟试验或现场试验确定杆体波速和杆系波速;
-确定杆端外露长度;
-利用杆体速度、杆系速度、杆端外露长度、初至时间和反射波初试时间计算锚杆长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110473896.5A CN113188488B (zh) | 2021-04-29 | 2021-04-29 | 锚杆长度无损测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110473896.5A CN113188488B (zh) | 2021-04-29 | 2021-04-29 | 锚杆长度无损测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113188488A true CN113188488A (zh) | 2021-07-30 |
CN113188488B CN113188488B (zh) | 2023-07-14 |
Family
ID=76980844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110473896.5A Active CN113188488B (zh) | 2021-04-29 | 2021-04-29 | 锚杆长度无损测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113188488B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20230120870A1 (en) * | 2021-10-18 | 2023-04-20 | Institute Of Rock And Soil Mechanics, Chinese Academy Of Sciences | Dual Channel Nondestructive Testing Method for Rock Bolt and Related Devices |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183226A (zh) * | 2011-03-11 | 2011-09-14 | 石家庄铁道大学 | 基于多源信息融合的锚杆无损检测方法 |
CN102269814A (zh) * | 2011-04-21 | 2011-12-07 | 长江大学 | 一种用于基桩质量智能检测的方法 |
CN102296645A (zh) * | 2011-06-02 | 2011-12-28 | 朱德兵 | 一种桩基质量检测方法及其装置 |
CN102890117A (zh) * | 2012-02-15 | 2013-01-23 | 长江水利委员会长江科学院 | 一种长锚索锚固质量检测方法及装置 |
CN104264722A (zh) * | 2014-09-30 | 2015-01-07 | 云南航天工程物探检测股份有限公司 | 成桥基桩无损检测方法 |
-
2021
- 2021-04-29 CN CN202110473896.5A patent/CN113188488B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183226A (zh) * | 2011-03-11 | 2011-09-14 | 石家庄铁道大学 | 基于多源信息融合的锚杆无损检测方法 |
CN102269814A (zh) * | 2011-04-21 | 2011-12-07 | 长江大学 | 一种用于基桩质量智能检测的方法 |
CN102296645A (zh) * | 2011-06-02 | 2011-12-28 | 朱德兵 | 一种桩基质量检测方法及其装置 |
CN102890117A (zh) * | 2012-02-15 | 2013-01-23 | 长江水利委员会长江科学院 | 一种长锚索锚固质量检测方法及装置 |
CN104264722A (zh) * | 2014-09-30 | 2015-01-07 | 云南航天工程物探检测股份有限公司 | 成桥基桩无损检测方法 |
Non-Patent Citations (1)
Title |
---|
郭庆: "声频应力波法检测土遗址加固用玻璃纤维锚杆锚固质量初探", 《文物保护与考古科学》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20230120870A1 (en) * | 2021-10-18 | 2023-04-20 | Institute Of Rock And Soil Mechanics, Chinese Academy Of Sciences | Dual Channel Nondestructive Testing Method for Rock Bolt and Related Devices |
US11796513B2 (en) * | 2021-10-18 | 2023-10-24 | Institute Of Rock And Soil Mechanics, Chinese Academy Of Sciences | Dual channel nondestructive testing method for rock bolt and related devices |
Also Published As
Publication number | Publication date |
---|---|
CN113188488B (zh) | 2023-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Meng et al. | Detecting bridge dynamics with GPS and triaxial accelerometers | |
CN105067239B (zh) | 基于扫频激励振动的梁裂纹故障检测装置及方法 | |
CN109459131A (zh) | 一种旋转机械多通道振动信号的时频特征提取方法及装置 | |
CN102425409B (zh) | 一种井斜和工具面测量装置及垂直钻井设备 | |
CN112985639B (zh) | 基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法 | |
Li et al. | Time-varying parameter identification of bridges subject to moving vehicles using ridge extraction based on empirical wavelet transform | |
CN112597577B (zh) | 桥梁索力高精度动态测量方法、装置和计算机设备 | |
CN113188488A (zh) | 锚杆长度无损测量方法 | |
CN111337171A (zh) | 一种应用于临界折射纵波应力检测的声时差测量方法 | |
CN113295310A (zh) | 一种基于应变刚度代表值的桥梁损伤判定方法 | |
CN104897777A (zh) | 基于Burg算法的自回归谱外推技术提高TOFD检测纵向分辨率的方法 | |
CN111736222A (zh) | 单炮数据信噪比确定方法及装置 | |
CN112269192A (zh) | 一种快速自适应的动态北斗监测实时解算去噪方法 | |
CN113189641B (zh) | 一种两道多模式瑞利波地下探测系统及方法 | |
Giurgiutiu et al. | Comparison of short-time fourier transform and wavelet transform of transient and tone burst wave propagation signals for structural health monitoring | |
Yan et al. | Adaptive synchroextracting transform and its application in bearing fault diagnosis | |
CN110333506B (zh) | 一种提取索力测量雷达的拉索位置参数的方法 | |
CN105138823A (zh) | 一种基于自相关函数的生理信号质量检测方法 | |
CN110568081B (zh) | 一种古建筑损伤测试装置及测试方法 | |
CN110261911B (zh) | 地震数据时差、相位差自动识别及校正方法 | |
CN114200519B (zh) | 一种无井区窄入射角道集获得密度的反演方法 | |
JP5335182B2 (ja) | 巨視的探査理論に基づいたコンクリート電柱地中部の広帯域超音波探査方法 | |
CN111307030A (zh) | 顾及有色噪声影响的gb-rar高铁桥梁形变信息估计模型 | |
CN112667952A (zh) | 一种结构动态位移非积分重构方法 | |
EP3156914B1 (en) | A method of performing real time decomposition of a signal into components |
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 | ||
CP01 | Change in the name or title of a patent holder |
Address after: No.109 Jinshui Road, Jinshui District, Zhengzhou City, Henan Province, 450000 Patentee after: Jianghe Anlan Engineering Consulting Co.,Ltd. Address before: No.109 Jinshui Road, Jinshui District, Zhengzhou City, Henan Province, 450000 Patentee before: Jianghe Engineering Inspection Co.,Ltd. |
|
CP01 | Change in the name or title of a patent holder |