CN108982650A - 一种钢轨顶面伤损漏磁数据的自适应滤波方法 - Google Patents
一种钢轨顶面伤损漏磁数据的自适应滤波方法 Download PDFInfo
- Publication number
- CN108982650A CN108982650A CN201810812431.6A CN201810812431A CN108982650A CN 108982650 A CN108982650 A CN 108982650A CN 201810812431 A CN201810812431 A CN 201810812431A CN 108982650 A CN108982650 A CN 108982650A
- Authority
- CN
- China
- Prior art keywords
- signal
- magnetic leakage
- leakage signal
- magnetic
- rail surface
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/72—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
- G01N27/82—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
Landscapes
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
本发明提出了一种钢轨顶面伤损漏磁数据的自适应滤波方法,属于高速铁路钢轨高速漏磁巡检技术领域。该方法包括如下步骤:步骤1,利用垂直于列车行进方向排成一行的路霍尔传感器采集钢轨顶面伤损漏磁信号,进行信号放大处理后存储,转化为数字电压信号;步骤2,利用小波滤波算法对步骤1得到的钢轨缺陷漏磁信号进行第一次滤波;步骤3,根据钢轨伤损漏磁信号构建虚拟通道,作为步骤4中的参考信号;步骤4,基于自适应滤波算法,对钢轨漏磁信号进行第二次滤波。本发明使自适应滤波算法成功运用于钢轨顶面伤损漏磁数据的滤波,钢轨顶面伤损检测时更易于特征提取,检出率更高,滤波效果明显,并提高了数据处理的效率精度。
Description
技术领域
本发明涉及一种钢轨顶面伤损漏磁数据的自适应滤波方法,属于高速铁路钢轨高速漏磁巡检技术领域。
背景技术
随着无损检测在铁路钢轨检测领域运用越来越广泛,漏磁检测技术因其传感器结构简单、检测灵敏度高、可实现非接触式检测,被广泛运用于诸如钢轨等铁磁性机件表面裂纹的无损检测。钢轨漏磁检测设备在采集和传输过程中由于机车振动、信号漂移、工频干扰、白噪声干扰等等情况的发生,使得其采集到的钢轨顶面缺陷漏磁信号都带有噪声、振动等干扰,而且干扰信号与缺陷信号频谱重叠,幅值相近,因此找到合适的滤波方法实现对漏磁信号的有效滤波具有重要意义。目前常用的滤波方法有频域滤波、卡尔曼滤波和自适应滤波,由于干扰信号与缺陷信号的特点,频域滤波和卡尔曼滤波都不适用于该漏磁信号的滤波。自适应滤波无需知道信号的先验统计特性,可以用于滤除漏磁信号中的干扰。然而自适应滤波需要不含有缺陷信息的通道的数据作为参考信号才能对漏磁信号进行有效滤波,而在实际钢轨巡检中,阵列传感器中每个传感器都有可能测到缺陷信息,若想利用自适应滤波方法,需要构建出合理的虚拟通道作为参考信号,才能对多通道阵列传感器采集到的漏磁信号进行最优滤波。
发明内容
针对现有技术的缺陷,本发明提出了一种钢轨顶面伤损漏磁数据的自适应滤波方法,使自适应滤波算法成功运用于钢轨顶面伤损漏磁数据的滤波,钢轨顶面伤损检测时更易于特征提取,检出率更高,滤波效果明显,并提高了数据处理的效率精度。
本发明为解决其技术问题采用如下技术方案:
一种钢轨顶面伤损漏磁数据的自适应滤波方法,包括如下步骤:
步骤1:漏磁信号采集
利用垂直于列车行进方向排成一行的n路霍尔传感器采集钢轨顶面伤损漏磁信号,进行信号放大处理后存储,并转化为数字电压信号;单路霍尔传感器采集到的漏磁信号为其中:为第i路通道第1个采样点,为第i路通道第2个采样点,为第i路通道第2个采样点,为第i路通道第j个采样点,为第i路通道第m个采样点,m为每个霍尔传感器的数据点数。由n路霍尔传感器采集到的漏磁信号构成漏磁信号矩阵X=(x1,x2,x3,...,xi,...,xn),其中,x1为第1路通道采集到的数据向量,x2为第2路通道采集到的数据向量,x3为第3路通道采集到的数据向量,xi为第i路通道采集到的数据向量,xn为第n路通道采集到的数据向量;
步骤2:小波去噪
基于小波滤波算法对钢轨漏磁信号进行第一次滤波,由漏磁信号矩阵X经过小波去噪后得到漏磁信号矩阵
步骤3:构建虚拟通道
根据步骤2所得漏磁信号找出的每行中绝对值最小的元素,记为其中j=1,2,3...m,为的第j行第1列的元素,为第j行第2列的元素,为第j行第3列的元素,为第j行最后一列的元素,为第j行第i列的元素;然后将找出的m个构成虚拟通道作为自适应滤波中的参考信号;其中:是第1行中绝对值最小的元素,是第2行中绝对值最小的元素,是第j行中绝对值最小的元素,是最后一行中绝对值最小的元素;
步骤4:基于自适应滤波算法,对钢轨漏磁信号进行第二次滤波
将步骤3中所得虚拟通道x0作为自适应滤波的参考信号,将步骤2中所得漏磁信号的每一个通道作为自适应滤波的输入信号,对13路漏磁信号中每个通道xi逐一进行单通道线滤波,最后实现整个漏磁信号矩阵的面滤波。
步骤1中,所述排成一行的n路霍尔传感器,排列长度需覆盖钢轨顶面宽度,n值由式子向上取整计算得出,其中D为钢轨顶面宽度,d为传感器宽度,s为伤损在垂直于列车行进方向上的分辨率。
步骤2的具体步骤如下:根据公式SNR为漏磁信号的信噪比,RMSE为漏磁信号的均方根误差,x(j)为原始纯净漏磁信号,为去噪后漏磁信号,m为某通道漏磁信号的采样点数,分别计算钢轨顶面漏磁检测数据的信噪比和均方根误差,由此确定小波去噪参数,通过计算选用的小波基是sym6小波,并进行10层小波分解。
步骤3中,所述虚拟通道x0由漏磁信号每行中绝对值最小的元素组成,其中虚拟通道x0第j行的元素是漏磁信号第j行中绝对值最小的元素,式中j=1,2,3,...,m。
本发明的有益效果如下:
1,本发明能够在缺陷信号频谱和噪声频谱重叠,和难以获取缺陷信号与干扰信号的先验统计特性的情况下,有效抑制噪声。根据实验计算得知,该自适应滤波方法能够使噪声强度降低81.44%,滤波效果明显。
2,该自适应滤波方法构建出了一个只含干扰信息而不含有缺陷信息的虚拟通道,作为参考信号,从而实现针对多通道阵列传感器采集到的漏磁信号的自适应滤波。因为利用漏磁信号的特点——检测出缺陷的信号幅值大于同一时刻未检测出缺陷的信号幅值;同时在自然缺陷中极少存在类似钢轨接缝的缺陷,即极少存在每路传感器都能检测到的缺陷,所以取每一时刻霍尔阵列传感器中幅值最小者构成虚拟通道,则该通道仅含干扰信息而不含有缺陷信息,即满足自适应滤波参考信号的选取原则,将该虚拟通道作为参考信息对漏磁信号进行滤波可实现最优滤波。
附图说明
图1(a)是原始带噪信号多通道三维图;图1(b)是自适应滤波后多通道效果图。
图2是滤波前后单通道对比图。
图3是钢轨顶面伤损漏磁数据自适应滤波方法流程图。
具体实施方式
以下结合附图及实例对本发明进行进一步详细说明。此处所描述的具体实例仅用于解释本发明,并不用于限定本发明。
参照图1,钢轨顶面伤损漏磁数据的自适应滤波方法应用于钢轨顶面缺陷漏磁数据的滤波,滤波方法包括以下步骤:
步骤1:漏磁信号采集。根据(其中D为钢轨顶面宽度,d为传感器宽度,s为裂纹在垂直于列车行进方向上的分辨率)计算确定n=13,因此采用13路霍尔传感器采集钢轨顶面伤损漏磁信号,并将信号放大处理后存储;单路霍尔传感器采集到的漏磁信号为其中:为第1时刻阵列传感器采集到的数据,为第2时刻阵列传感器采集到的数据,为第3时刻阵列传感器采集到的数据,为第j时刻阵列传感器采集到的数据,为最后时刻阵列传感器采集到的数据,共有4401个数据采样点数),13路霍尔传感器采集到的漏磁信号构成漏磁信号矩阵X=(x1,x2,x3,...,xi,...,x16),其中:x1为第1个霍尔传感器采集到的数据,x2为第2个霍尔传感器采集到的数据,x3为第3个霍尔传感器采集到的数据,xi为第i个霍尔传感器采集到的数据,x13为第13个霍尔传感器采集到的数据;
步骤2:小波去噪。基于小波滤波算法对钢轨漏磁信号进行第一次滤波,滤除步骤1中所得信号X的细节信息得到漏磁信号矩阵该步骤根据公式 其中:SNR为漏磁信号的信噪比,RMSE为漏磁信号的均方根误差,x(j)为原始纯净漏磁信号,为去噪后漏磁信号,m为某通道漏磁信号的采样点数,分别计算钢轨顶面漏磁检测数据的信噪比和均方根误差,由此确定小波去噪参数,即选用的小波基是sym6小波,进行10层小波分解,软阈值法对钢轨顶面漏磁检测数据进行滤波。
步骤3:构建虚拟通道。根据步骤2所得漏磁信号找出的每行中绝对值最小的元素,记为其中j=1,2,3...,4401,为的第j行第1列的元素,为第j行第2列的元素,为第j行第3列的元素,为为第j行最后一列的元素,为中第j行第i列的元素;最后将找出的4401个数据构成虚拟通道作为自适应滤波中的参考信号,其中:为第1行中绝对值最小的元素,为第2行中绝对值最小的元素,为第j行中绝对值最小的元素,为最后一行中绝对值最小的元素;
步骤4:基于自适应滤波算法,对钢轨漏磁信号进行第二次滤波。将步骤3中所得虚拟通道x0作为自适应滤波的参考信号,从第1路漏磁信号开始,逐一对13路漏磁信号xi进行单通道线滤波,从而实现整个漏磁信号矩阵的面滤波;
结果发现,之前几乎被干扰信号“淹没”的八条缺陷信号在进行自适应滤波之后完整浮现出来了,并且更加易于缺陷的识别与查找,同时定量分析了滤波前后纯噪声能量的变化率,发现噪声强度降低了81.44%,显著提高了漏磁信号的信噪比。
Claims (4)
1.一种钢轨顶面伤损漏磁数据的自适应滤波方法,其特征在于,包括如下步骤:
步骤1:漏磁信号采集
利用垂直于列车行进方向排成一行的n路霍尔传感器采集钢轨顶面伤损漏磁信号,进行信号放大处理后存储,并转化为数字电压信号;单路霍尔传感器采集到的漏磁信号为其中:为第i路通道第1个采样点,为第i路通道第2个采样点,为第i路通道第2个采样点,为第i路通道第j个采样点,为第i路通道最后一个采样点,m为每个霍尔传感器采集的数据点数;由n路霍尔传感器采集到的漏磁信号构成漏磁信号矩阵X=(x1,x2,x3,...,xi,...,xn),其中,x1为第1路通道采集到的数据向量,x2为第2路通道采集到的数据向量,x3为第3路通道采集到的数据向量,xi为第i路通道采集到的数据向量,xn为最后一路通道采集到的数据向量;
步骤2:小波去噪
基于小波滤波算法对钢轨顶面漏磁信号进行第一次滤波,由漏磁信号矩阵X经过小波去噪后得到漏磁信号矩阵
步骤3:构建虚拟通道
根据步骤2所得漏磁信号找出的每行中绝对值最小的元素,记为其中j=1,2,3...m,为的第j行第1列的元素,为第j行第2列的元素,为第j行第3列的元素,为第j行最后一列的元素,为第j行第i列的元素;然后将找出的m个构成虚拟通道作为自适应滤波中的参考信号;其中:r1 0是第1行中绝对值最小的元素,是第2行中绝对值最小的元素,是第j行中绝对值最小的元素,是最后一行中绝对值最小的元素;
步骤4:基于自适应滤波算法,对钢轨漏磁信号进行第二次滤波将步骤3中所得虚拟通道x0作为自适应滤波的参考信号,将步骤2中所得漏磁信号的每一个通道作为自适应滤波的输入信号,对13路漏磁信号中每个通道xi逐一进行单通道线滤波,最后实现整个漏磁信号矩阵的面滤波。
2.根据权利要求1所述的一种钢轨顶面伤损漏磁数据的自适应滤波方法,其特征在于,步骤1中所述排成一行的n路霍尔传感器,排列长度需覆盖钢轨顶面宽度,n值由式子向上取整计算得出,其中D为钢轨顶面宽度,d为传感器宽度,s为伤损在垂直于列车行进方向上的分辨率。
3.根据权利要求1所述的一种钢轨顶面伤损漏磁数据的自适应滤波方法,其特征在于,步骤2的具体步骤如下:根据公式SNR为漏磁信号的信噪比,RMSE为漏磁信号的均方根误差,x(j)为原始纯净漏磁信号,为去噪后漏磁信号,m为单通道漏磁信号的采样点数,分别计算钢轨顶面漏磁检测数据的信噪比和均方根误差,由此确定小波去噪参数,通过计算选用的小波基是sym6小波,并进行10层小波分解。
4.根据权利要求1所述的一种钢轨顶面伤损漏磁数据的自适应滤波方法,其特征在于,步骤3中,所述虚拟通道x0由漏磁信号每行中绝对值最小的元素组成,其中虚拟通道x0第j行的元素是漏磁信号第j行中绝对值最小的元素,式中j=1,2,3,...,m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810812431.6A CN108982650A (zh) | 2018-07-23 | 2018-07-23 | 一种钢轨顶面伤损漏磁数据的自适应滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810812431.6A CN108982650A (zh) | 2018-07-23 | 2018-07-23 | 一种钢轨顶面伤损漏磁数据的自适应滤波方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108982650A true CN108982650A (zh) | 2018-12-11 |
Family
ID=64549371
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810812431.6A Pending CN108982650A (zh) | 2018-07-23 | 2018-07-23 | 一种钢轨顶面伤损漏磁数据的自适应滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108982650A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112034036A (zh) * | 2020-10-16 | 2020-12-04 | 中国铁道科学研究院集团有限公司 | 钢轨漏磁信号滤波方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030020469A1 (en) * | 2001-07-17 | 2003-01-30 | Gopichand Katragadda | Transverse crack detection in rail head using low frequency eddy currents |
CN102832908A (zh) * | 2012-09-20 | 2012-12-19 | 西安科技大学 | 基于小波变换与变步长lms自适应滤波的信号降噪方法 |
CN103997315A (zh) * | 2014-05-16 | 2014-08-20 | 东北大学 | 一种管道漏磁信号自适应滤波装置及方法 |
CN106353397A (zh) * | 2016-09-06 | 2017-01-25 | 中国铁道科学研究院 | 一种车载的钢轨顶面伤损漏磁检测的设备以及系统 |
CN106645384A (zh) * | 2016-09-23 | 2017-05-10 | 东北大学 | 一种管道漏磁内检测器数据的自适应滤波方法 |
-
2018
- 2018-07-23 CN CN201810812431.6A patent/CN108982650A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030020469A1 (en) * | 2001-07-17 | 2003-01-30 | Gopichand Katragadda | Transverse crack detection in rail head using low frequency eddy currents |
CN102832908A (zh) * | 2012-09-20 | 2012-12-19 | 西安科技大学 | 基于小波变换与变步长lms自适应滤波的信号降噪方法 |
CN103997315A (zh) * | 2014-05-16 | 2014-08-20 | 东北大学 | 一种管道漏磁信号自适应滤波装置及方法 |
CN106353397A (zh) * | 2016-09-06 | 2017-01-25 | 中国铁道科学研究院 | 一种车载的钢轨顶面伤损漏磁检测的设备以及系统 |
CN106645384A (zh) * | 2016-09-23 | 2017-05-10 | 东北大学 | 一种管道漏磁内检测器数据的自适应滤波方法 |
Non-Patent Citations (2)
Title |
---|
W. CHEN ET AL.: "Rail crack recognition based on Adaptive Weighting Multi-classifier Fusion Decision", 《MEASUREMENT》 * |
丁顺意 等: "钢轨顶面伤损漏磁数据的自适应滤波方法", 《2018远东无损检测新技术论坛论文集》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112034036A (zh) * | 2020-10-16 | 2020-12-04 | 中国铁道科学研究院集团有限公司 | 钢轨漏磁信号滤波方法及装置 |
CN112034036B (zh) * | 2020-10-16 | 2023-11-17 | 中国铁道科学研究院集团有限公司 | 钢轨漏磁信号滤波方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104076404B (zh) | 运用多通道相干抑制地磁背景噪声的磁异常探测方法 | |
CN102636577A (zh) | 一种基于小波变换提取轮辋超声波探伤信号的方法 | |
CN106885846B (zh) | 树木缺陷检测装置及检测方法 | |
CN103344989B (zh) | 可控震源地震记录中脉冲噪声干扰的分析方法 | |
CN103997315B (zh) | 一种管道漏磁信号自适应滤波装置及方法 | |
CN106501372B (zh) | 基于小波包分析道岔裂纹的监测定位方法 | |
CN102269814A (zh) | 一种用于基桩质量智能检测的方法 | |
CN111103621A (zh) | 一种主动源共成像点叠加多道面波分析方法 | |
CN105629300B (zh) | 提高复杂构造偏移数据信噪比的方法 | |
Zhang et al. | Data analysis technique to leverage ground penetrating radar ballast inspection performance | |
CN110389377B (zh) | 基于波形互相关系数相乘的微震偏移成像定位方法 | |
Lee et al. | Resolution analysis of shallow marine seismic data acquired using an airgun and an 8-channel streamer cable | |
CN108797662A (zh) | 立柱下桩基长度无损检测方法及装置 | |
Chambers | Using DAS to investigate traffic patterns at Brady Hot Springs, Nevada, USA | |
CN103645248A (zh) | 一种基于超声相速度的高温合金晶粒度评价方法 | |
CN108982650A (zh) | 一种钢轨顶面伤损漏磁数据的自适应滤波方法 | |
CN108254788A (zh) | 一种地震波初至拾取方法及系统 | |
Zhang et al. | Railway traffic monitoring with trackside fiber-optic cable by distributed acoustic sensing Technology | |
CN112034036B (zh) | 钢轨漏磁信号滤波方法及装置 | |
CN108375789B (zh) | 联合采集地震数据的同步匹配方法 | |
CN101957458A (zh) | 一种处理超低频数据的方法 | |
CN106771928A (zh) | 一种局部放电脉冲初至时刻在线拾取方法 | |
CN108387640A (zh) | 一种基于二维多次覆盖超声反射波检测方法 | |
CN104267094A (zh) | 铁磁性构件的脉冲漏磁响应信号分离方法 | |
McCauley et al. | Migratory patterns and estimated population size of pygmy blue whales (Balaenoptera musculus brevicauda) traversing the Western Australian coast based on passive acoustics. |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181211 |
|
RJ01 | Rejection of invention patent application after publication |