CN1696687A - 粗晶材料超声检测时频分析处理方法 - Google Patents
粗晶材料超声检测时频分析处理方法 Download PDFInfo
- Publication number
- CN1696687A CN1696687A CNA2005100117622A CN200510011762A CN1696687A CN 1696687 A CN1696687 A CN 1696687A CN A2005100117622 A CNA2005100117622 A CN A2005100117622A CN 200510011762 A CN200510011762 A CN 200510011762A CN 1696687 A CN1696687 A CN 1696687A
- Authority
- CN
- China
- Prior art keywords
- time
- frequency
- signal
- frequency image
- exp
- 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
- 239000000463 material Substances 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000012360 testing method Methods 0.000 title description 5
- 238000010183 spectrum analysis Methods 0.000 title description 3
- 230000007547 defect Effects 0.000 claims abstract description 26
- 238000001514 detection method Methods 0.000 claims abstract description 12
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 4
- 230000021615 conjugation Effects 0.000 claims description 3
- 238000011002 quantification Methods 0.000 claims description 3
- SYHGEUNFJIGTRX-UHFFFAOYSA-N methylenedioxypyrovalerone Chemical compound C=1C=C2OCOC2=CC=1C(=O)C(CCC)N1CCCC1 SYHGEUNFJIGTRX-UHFFFAOYSA-N 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 4
- 230000000694 effects Effects 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 description 26
- 238000001228 spectrum Methods 0.000 description 21
- 230000002950 deficient Effects 0.000 description 8
- 229910000963 austenitic stainless steel Inorganic materials 0.000 description 6
- 230000007812 deficiency Effects 0.000 description 6
- 238000013459 approach Methods 0.000 description 3
- 206010061619 Deformity Diseases 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 208000037656 Respiratory Sounds Diseases 0.000 description 1
- 229910001069 Ti alloy Inorganic materials 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
一种粗晶材料超声检测的时频分析处理方法,属于超声无损检测技术领域。先获得优质时频图像,再时频图像缺陷信息提取,最后缺陷信息的A型显示。本发明共由三个技术步骤组成,因此简称为“三步法”。该方法具有更强的小缺陷发现能力和更好的信噪比增强效果,且克服了分离谱技术的参数敏感性问题,比传统的分离谱技术具有更强的缺陷发现能力,可以检测粗晶材料中更微小的缺陷,具有非常好的信噪比增强效果,且由于摒弃了分离谱技术的非线性统计处理,避免了分离谱技术的参数敏感性问题,在航空航天、核工业、石油化工等大量应用粗晶材料的重要部门具有广阔的应用前景。
Description
技术领域
本发明涉及一种粗晶材料超声检测的时频分析处理方法,属于超声无损检测技术领域。
背景技术
粗晶材料(如奥氏体不锈钢和钛合金)由于在核工业和航空航天工业中的大量应用而备受瞩目,其安全性检测十分重要。但由于粗大晶粒对超声波的强烈散射作用,粗晶材料的超声反射信号伴有严重的材料噪声,致使对粗晶材料内部缺陷(如裂纹)的检测十分困难。粗晶材料的超声检测是无损检测领域的一项重大技术难题。
目前,分离谱技术是粗晶材料超声检测的主导技术,该技术于1982年由美国V.L.Newhouse博士发明,一直沿用至今。分离谱技术存在两大本质技术缺陷。其一,分离谱技术对小缺陷的发现能力不足,目前的水平最高只能检测到直径为1.5mmd的缺陷,对于更小的缺陷则难以发现。这一技术缺陷的根本原因是由于分离谱对超声信号进行的是恒定时频分辨率的观测,无法凸现缺陷信号与噪声信号在低频端的频谱差异。分离谱技术的另一技术缺陷是,它对相关处理参数十分敏感,性能不稳定,人工经验依赖性强,无法实现对超声信号的自动化分析和处理。这一技术缺陷的根本原因是分离谱采用非线性统计处理对噪声加以抑制,处理参数的变化可能同时造成对缺陷信号的抑制。因此,分离谱技术缺乏实用化基础,目前仍无法在生产中获得推广应用。
近年来,有学者将连续小波分析技术引入超声信号处理中。连续小波分析实现对超声信号的可变时频分辨率观测,可以更好的凸现缺陷信号与噪声信号在低频段的频谱差异。因此,连续小波分析技术具有小缺陷检测的潜在能力。但当前连续小波分析技术进行噪声抑制的方法是基于门限阈值处理的能量消减法,由于门限阈值处理在抑制噪声信号的同时消减了大量的缺陷信息,因此连续小波降噪技术会造成严重的信号失真,而且不具备缺陷信息的增强和放大能力,不能满足粗晶材料超声检测的需要。分离谱技术仍在粗晶材料超声检测中处于主导地位。
发明内容
针对分离谱技术和连续小波降噪技术存在的问题,本发明提出了一种更适合粗晶材料超声检测的时频分析处理方法。
粗晶材料超声检测时频分析处理方法由三个步骤组成:
第一步:获得优质时频图像:以Q值优化(Q=1/4.8)的高斯小波对超声反射信号进行常Q时频分析,获得高质量的超声信号时频图像;
优化的高斯母小波为
其中,优化参数Bm=0.58(即Q=1/4.8),fc=1。
采用优化后的高斯连续小波对粗晶材料超声反射信号x(t)按下式进行时频分解,这样便可以获得由分解数据WT(τ,a)所构成的超声信号时频图像。
式中,*号代表复数取共轭;a为尺度因子;τ为时移变量。
第二步:时频图像缺陷信息提取,即从时频图像中提取材料内部缺陷的量化信息:以匹配追踪算法实现对时频图像中各行子带信号的自适应时频分解,从而获得缺陷信号在各子带内的幅度Ak、脉冲宽度σk、时间中心tk、频率中心fk等重要的量化信息;
第三步:缺陷信息的A型显示:统计子带缺陷信号在时间中心tk出现的概率,并将不同tk处的概率值以曲线形式加以表达,即获得缺陷信息的A型显示。
本发明的第一步获得优质时频图像,继承了连续小波分析对缺陷信号与噪声信号在低频段频谱差异的凸现能力。同时,本发明实现了对连续小波分析关键参数(Q值)的优化,使对超声信号的观测具有了更高的时频分辨率,优化了超声信号的时频图像。
本发明的核心是其第二步,时频图像缺陷信息提取,即从时频图像中提取材料内部缺陷的量化信息。这一步摒弃了分离谱技术的非线性统计方法和连续小波降噪技术的门限阈值处理方法,而是采用匹配追踪算法,从超声信号的时频图像中提取缺陷信息,这样可以更完整地保留缺陷信息并加以量化,而不像门限阈值处理那样将低于门限阈值的缺陷信息和噪声一并删除。
本发明的第三步缺陷信息的A型显示,利用第二步获得的缺陷量化信息,以A型显示方式对缺陷加以表达,符合超声检测领域工程显示方式的要求,便于操作人员的理解和判断。
由于本发明共由三个技术步骤组成,因此可以简称为“三步法”。该方法具有更强的小缺陷发现能力和更好的信噪比增强效果,且克服了分离谱技术的参数敏感性问题,比传统的分离谱技术具有更强的缺陷发现能力,可以检测粗晶材料中更微小的缺陷,具有非常好的信噪比增强效果,且由于摒弃了分离谱技术的非线性统计处理,避免了分离谱技术的参数敏感性问题,具有很强的实用性能。本发明是粗晶材料超声检测实用而方便的新型处理技术,在航空航天、核工业、石油化工等大量应用粗晶材料的重要部门具有广阔的应用前景。
附图说明
图1是一种粗晶奥氏体不锈钢获得优质时频图像的处理结果;其中1(a)为原始的低信噪比超声反射信号,1(b)为优化高斯小波分析的时频图像;
图2是对图1所示奥氏体不锈钢时频图像缺陷信息提取后重构而成的时频图像;其中2(a)为原始的低信噪比超声反射信号;2(b)为重构而成的时频图像;
图3是缺陷信息A型显示的对比图;其中3(a)为原始的低信噪比超声反射信号,3(b)为本发明(“三步法”)处理结果的A型显示,3(c)为采用传统的分离谱技术所获得的结果。
具体实施方式
附图所用检测对象为粗晶奥氏体不锈钢(牌号为OCr18Ni9Ti),利用传统的分离谱技术无法检测其内部含有的相距5mm的两个φ0.5mm横通孔缺陷。。
按本发明的“三步法”原理,具体实施过程如下:
第一步:获得优质时频图像。
采用连续小波变换对粗晶材料超声反射信号x(t)进行时频分解:
式中,*号代表复数取共轭;a为尺度因子;τ为时移变量,优化的高斯母小波为
其中,优化参数Bm=0.58(即Q=1/4.8),fc=1。这样便可以获得由分解数据WT(τ,a)所构成的超声信号时频图像,如图1所示。
图1是对粗晶奥氏体不锈钢获得优质时频图像的处理结果;其中1(a)为原始的低信噪比超声反射信号,1(b)为优化高斯小波分析的时频图像;也即是本发明第一步的处理结果。
第二步:时频图像缺陷信息提取。
采用匹配追踪算法对获得的时频图像进行缺陷信息提取。匹配追踪是语音信号处理领域已公开的时频分析技术,已经有相当成熟的快速算法。匹配追踪只需预先确定一个完备的函数集合(不需要其它技术参数),即可以快速的将待处理信号分解为相应函数集中与信号最匹配的基函数的线性组合。对于超声信号,本发明选择高斯函数集,其基函数形式如下
这样,便可以将时频图像中各子带超声信号fa(t)分解为少数几个高能量信号原子的线性组合
分解程度可根据信号的残余能量设定,通常当信号残余能量小于原始信号能量的1%时,便可认为残余为噪声,终止匹配追踪过程。由此,便可以获得缺陷信号在各子带内的幅度Ak、脉冲宽度σk、时间中心tk、频率中心fk等量化信息。
图2为利用本发明第二步所提取的信息重构而成的时频图像,即是对一种粗晶奥氏体不锈钢时频图像缺陷信息提取重构而成的时频图像;其中2(a)为原始的低信噪比超声反射信号;2(b)为重构而成的时频图像。采用匹配追踪算法对图1中时频图像所含缺陷信息的提取,可以获得缺陷信号在各子带内的幅度Ak、脉冲宽度σk、时间中心tk、频率中心fk等量化信息。图2是对这些信息的综合反映。由图可见,在高尺度范围(即低频范围)内,噪声信号被很好地消除掉了,同时缺陷信息得到了完整的保留。
第三步:缺陷信息的A型显示。统计子带缺陷信号在时间中心tk出现的概率,并将不同tk处的概率值以曲线形式加以表达,即获得缺陷信息的A型显示。
图3是缺陷信息A型显示的对比图。图3(a)为原始的低信噪比超声反射信号;
图3.(b)为本发明第三步的最终处理结果。由图可见,本项技术发明可以帮助检测人员非常清晰地对缺陷加以判别。
图3.(c)为采用传统的分离谱技术所获得的结果。由图可见,分离谱的处理结果不足以使检测人员对缺陷的存在和位置做出明确的判定,对缺陷信息的增强和定位能力都不及本发明。
Claims (1)
1、一种粗晶材料超声检测的时频分析处理方法,其特征在于,该方法由三个步骤组成:
第1步:获得优质时频图像:以Q值优化(Q=1/4.8)的高斯小波对超声反射信号进行常Q时频分析,获得高质量的超声信号时频图像;
优化的高斯母小波为
其中,优化参数Bm=0.58(即Q=1/4.8),fc=1;
采用优化后的高斯连续小波对粗晶材料超声反射信号x(t)按下式进行时频分解,这样便可以获得由分解数据WT(τ,a)所构成的超声信号时频图像;
式中,*号代表复数取共轭;a为尺度因子;τ为时移变量;
第2步:时频图像缺陷信息提取,即从时频图像中提取材料内部缺陷的量化信息:以匹配追踪算法实现对时频图像中各行子带信号的自适应时频分解,从而获得缺陷信号在各子带内的幅度Ak、脉冲宽度σk、时间中心tk、频率中心fk等重要的量化信息;
第2.1步:选择高斯函数集,其基函数形式如下
第2.2步:将时频图像中各子带超声信号fa(t)分解为少数几个高能量信号原子的线性组合
分解程度可根据信号的残余能量设定,通常当信号残余能量小于原始信号能量的1%时,便可认为残余为噪声,终止匹配追踪过程;由此,便可以获得缺陷信号在各子带内的幅度Ak、脉冲宽度σk、时间中心tk、频率中心fk等量化信息;
第3步:缺陷信息的A型显示:统计子带缺陷信号在时间中心tk出现的概率,并将不同tk处的概率值以曲线形式加以表达,即获得缺陷信息的A型显示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100117622A CN100405384C (zh) | 2005-05-23 | 2005-05-23 | 粗晶材料超声检测时频分析处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100117622A CN100405384C (zh) | 2005-05-23 | 2005-05-23 | 粗晶材料超声检测时频分析处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1696687A true CN1696687A (zh) | 2005-11-16 |
CN100405384C CN100405384C (zh) | 2008-07-23 |
Family
ID=35349506
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100117622A Expired - Fee Related CN100405384C (zh) | 2005-05-23 | 2005-05-23 | 粗晶材料超声检测时频分析处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100405384C (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101109732B (zh) * | 2007-08-08 | 2010-06-30 | 哈尔滨工业大学 | 基于模糊平面特征的超声无损检测回波信号分类方法 |
CN101726545B (zh) * | 2009-12-23 | 2011-11-23 | 哈尔滨工业大学 | 异种材料扩散焊界面超声信号时频幅度特征和时频相位特征提取方法 |
CN105181805A (zh) * | 2015-09-30 | 2015-12-23 | 中国计量学院 | 一种基于时反music的多滤波超声成像方法 |
CN105361855A (zh) * | 2016-01-11 | 2016-03-02 | 东南大学 | 一种有效获取脑磁图信号中事件相关磁场信息的方法 |
CN105954360A (zh) * | 2016-06-03 | 2016-09-21 | 河北省电力建设调整试验所 | 一种20Cr1Mo1VNbTiB高温螺栓粗晶超声波检验方法 |
CN106815552A (zh) * | 2016-12-09 | 2017-06-09 | 云南航天工程物探检测股份有限公司 | 基于时频分析的数字信号后处理方法 |
CN111855802A (zh) * | 2020-07-28 | 2020-10-30 | 武汉大学 | 一种消除激光超声行进波的缺陷可视化成像方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI431271B (zh) | 2010-07-15 | 2014-03-21 | Ind Tech Res Inst | 瑕疵檢測系統及方法 |
CN102346171B (zh) * | 2010-07-30 | 2013-09-04 | 财团法人工业技术研究院 | 瑕疵检测***及方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4224035A1 (de) * | 1992-07-21 | 1994-01-27 | Siemens Ag | Ultraschallprüfverfahren |
WO2004034231A2 (en) * | 2002-10-11 | 2004-04-22 | Flint Hills Scientific, L.L.C. | Intrinsic timescale decomposition, filtering, and automated analysis of signals of arbitrary origin or timescale |
-
2005
- 2005-05-23 CN CNB2005100117622A patent/CN100405384C/zh not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101109732B (zh) * | 2007-08-08 | 2010-06-30 | 哈尔滨工业大学 | 基于模糊平面特征的超声无损检测回波信号分类方法 |
CN101726545B (zh) * | 2009-12-23 | 2011-11-23 | 哈尔滨工业大学 | 异种材料扩散焊界面超声信号时频幅度特征和时频相位特征提取方法 |
CN105181805A (zh) * | 2015-09-30 | 2015-12-23 | 中国计量学院 | 一种基于时反music的多滤波超声成像方法 |
CN105361855A (zh) * | 2016-01-11 | 2016-03-02 | 东南大学 | 一种有效获取脑磁图信号中事件相关磁场信息的方法 |
CN105954360A (zh) * | 2016-06-03 | 2016-09-21 | 河北省电力建设调整试验所 | 一种20Cr1Mo1VNbTiB高温螺栓粗晶超声波检验方法 |
CN106815552A (zh) * | 2016-12-09 | 2017-06-09 | 云南航天工程物探检测股份有限公司 | 基于时频分析的数字信号后处理方法 |
CN106815552B (zh) * | 2016-12-09 | 2020-03-17 | 云南航天工程物探检测股份有限公司 | 基于时频分析的数字信号后处理方法 |
CN111855802A (zh) * | 2020-07-28 | 2020-10-30 | 武汉大学 | 一种消除激光超声行进波的缺陷可视化成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100405384C (zh) | 2008-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1696687A (zh) | 粗晶材料超声检测时频分析处理方法 | |
Song et al. | Wavelet based noise suppression technique and its application to ultrasonic flaw detection | |
CN105223482B (zh) | 局部放电特高频信号波形的小波分解二值去噪方法 | |
CN101799916A (zh) | 基于贝叶斯估计的生物芯片图像小波去噪方法 | |
CN114799610B (zh) | 一种基于傅里叶逆变换及自编码器的焊接质量实时检测方法及*** | |
CN117434153B (zh) | 基于超声波技术的道路无损检测方法及*** | |
CN115420806A (zh) | 一种基于神经网络及图像融合的无损超声检测方法 | |
CN111579634B (zh) | 一种铁磁性材料表面残余应力的抽样评估方法 | |
CN115655455A (zh) | 一种基于自适应噪声变换和随机共振的机械故障诊断方法 | |
Zhang et al. | Denoising of gamma-ray signals by interval-dependent thresholds of wavelet analysis | |
Yumnam et al. | Empirical mode decomposition based techniques for imaging of shallow delamination in concrete using impact echo | |
Wang et al. | The optimal wavelet threshold de-nosing method for acoustic emission signals during the medium strain rate damage process of concrete | |
Prajna et al. | Efficient harmonic regeneration noise reduction‐based Wiener filter for acoustic emission signal detection | |
Yao et al. | Research on wavelet denoising method based on soft threshold in wire rope damage detection | |
CN109932697B (zh) | 基于压缩感知的微弱目标检测方法 | |
CN112183260A (zh) | 一种基于总变差降噪与rqa的单向阀故障诊断方法 | |
Lan et al. | Improved wavelet packet noise reduction for microseismic data via fuzzy partition | |
CN115753107A (zh) | 基于非线性量子信息熵的微弱滚动轴承故障诊断方法 | |
CN108254333A (zh) | 基于红外光声光谱的含氮类农药残留快速检测方法 | |
CN112788057A (zh) | 一种基于FSWT时频分布的LDoS攻击检测方法 | |
Yu et al. | The research on the application of electronic nose in discriminate the rice varieties | |
CN117972536B (zh) | 一种脉冲分类方法及*** | |
Li et al. | A novel denoising approach based on multiresolution thresholding and its NDT application | |
El Tokhy et al. | Automatic detection algorithm of defects in casting radiography images based on cepstral coefficients | |
CN114384156B (zh) | 一种基于信息熵理论的声发射信号特征提取与裂纹损伤识别方法 |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080723 |