CN108918678B - 利用非线性超声检测物品劣化参数的方法及装置 - Google Patents
利用非线性超声检测物品劣化参数的方法及装置 Download PDFInfo
- Publication number
- CN108918678B CN108918678B CN201810502786.5A CN201810502786A CN108918678B CN 108918678 B CN108918678 B CN 108918678B CN 201810502786 A CN201810502786 A CN 201810502786A CN 108918678 B CN108918678 B CN 108918678B
- Authority
- CN
- China
- Prior art keywords
- coefficient
- calibration sample
- wavelet
- ultrasonic
- equation
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Signal Processing (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Acoustics & Sound (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明实施例提供了一种利用非线性超声检测物品劣化参数的方法及装置,可以对待测物品进行非线性超声检测,获得所述待测物品输出的超声信号;对所述超声信号进行小波变换,获得多个频段对应的小波系数;根据所述小波系数确定所述待测物品的非线性超声系数;将所述非线性超声系数和小波系数作为自变量带入与所述待测物品的材料对应的第一方程中,获得所述待测物品的第一劣化参数的取值。由于本发明使用了非线性超声检测待测物品的劣化参数,因此无需对待测物品进行破坏。同时,本发明对超声信号进行小波变换,具有多分辨率和局部分析的优势,提高了检测的准确性和稳定性。
Description
技术领域
本发明涉及超声检测领域,特别涉及利用非线性超声检测物品劣化参数的方法及装置。
背景技术
承受高温、高压工况的工件大量采用金属材料,如发电厂、化工厂的高温高压管道、管件、阀门等。随着服役时间延长,金属材料逐渐老化,导致金属工件的安全性和稳定性降低。
通过检测金属工件的劣化参数可以确定金属工件的老化程度,进而可以及时对老化较严重的金属工件进行及时的更换。现有的劣化参数检测方法一般为取样检测法,即:通过破坏性方法(如割管、整体剖切等)对待测金属工件取样,然后通过长时高温蠕变试验(试验周期数千到数万小时)、短时高温拉伸试验等试验手段对取样的工件的劣化参数进行检测。
由于现有技术需要对待测金属工件取样,因此可能会损坏待测金属工件。同时,除管道部件外的其他承受高温、高压的部件(如联箱、阀体、缸体等)通常不允许破坏性取样。当然,除金属工件外,其他材料的工件也可能存在相同问题。
因此,如何在不对工件进行损坏的前提下实现劣化参数的检测仍是本领域一个亟待解决的技术难题。
发明内容
有鉴于此,本发明提供一种利用非线性超声检测物品劣化参数的方法及装置,以通过非线性超声检测物品的劣化参数。
为了实现上述发明目的,本发明提供以下技术方案:
一种利用非线性超声检测物品劣化参数的方法,包括:
对待测物品进行非线性超声检测,获得所述待测物品输出的超声信号;
对所述超声信号进行小波变换,获得多个频段对应的小波系数;
根据所述小波系数确定所述待测物品的非线性超声系数;
将所述非线性超声系数和所述小波系数作为自变量带入与所述待测物品的材料对应的第一方程中,获得所述待测物品的第一劣化参数的取值,其中,所述第一方程的因变量为所述第一劣化参数的取值且所述第一方程与所述第一劣化参数对应。
可选的,在所述对待测物品进行非线性超声检测前,所述方法还包括:
获得材料与所述待测物品相同的多个定标样品,其中,各定标样品的第一劣化参数的取值已知;
对每一个定标样品:对该定标样品进行非线性超声检测,获得该定标样品输出的超声信号,对该定标样品输出的超声信号进行小波变换,获得该定标样品的多个频段对应的小波系数,根据该定标样品的多个频段对应的小波系数确定该定标样品的非线性超声系数;
根据各定标样品的第一劣化参数的取值、小波系数以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数的取值,其中,所述待拟合方程为第一劣化参数与非线性超声系数、小波系数的关联方程;
将所述待拟合方程中的未知系数修改为所确定的取值,获得所述第一方程。
可选的,所述小波系数为Wi,其中,i为小波系数的编号,i为自然数,所述根据该定标样品的多个频段对应的小波系数确定该定标样品的非线性超声系数,包括:
根据公式
β=W2/W1 2
γ=W3/W1 3
计算得到该定标样品的非线性超声系数γ和β,其中,W1为对该定标样品进行非线性超声检测时使用的发射端超声换能器的中心频率F1所在的频段对应的小波系数,W2为频率F2所在的频段对应的小波系数,W3为频率F3所在的频段对应的小波系数,其中,F2=2F1,F3=3F1。
可选的,所述待拟合方程为:
可选的,所述根据各定标样品的第一劣化参数的取值、小波系数以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数的取值,包括:
根据各定标样品的第一劣化参数的取值以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数k1和k2的取值;
对每个定标样品,根据公式:
根据各定标样品的残差以及小波系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的fi的取值。
一种利用非线性超声检测物品劣化参数的装置,包括:信号获得单元、小波系数获得单元、超声系数获得单元和劣化参数获得单元,
所述信号获得单元,用于对待测物品进行非线性超声检测,获得所述待测物品输出的超声信号;
所述小波系数获得单元,用于对所述超声信号进行小波变换,获得多个频段对应的小波系数;
所述超声系数获得单元,用于根据所述小波系数确定所述待测物品的非线性超声系数;
所述劣化参数获得单元,用于将所述非线性超声系数和所述小波系数作为自变量带入与所述待测物品的材料对应的第一方程中,获得所述待测物品的第一劣化参数的取值,其中,所述第一方程的因变量为所述第一劣化参数的取值且所述第一方程与所述第一劣化参数对应。
可选的,所述装置还包括:样品获得单元、样品检测单元、系数拟合单元和方程获得单元,
所述样品获得单元,用于在所述信号获得单元对待测物品进行非线性超声检测前,获得材料与所述待测物品相同的多个定标样品,其中,各定标样品的第一劣化参数的取值已知;
所述样品检测单元,用于对每一个定标样品:对该定标样品进行非线性超声检测,获得该定标样品输出的超声信号,对该定标样品输出的超声信号进行小波变换,获得该定标样品的多个频段对应的小波系数,根据该定标样品的多个频段对应的小波系数确定该定标样品的非线性超声系数;
所述系数拟合单元,用于根据各定标样品的第一劣化参数的取值、小波系数以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数的取值,其中,所述待拟合方程为第一劣化参数与非线性超声系数、小波系数的关联方程;
所述方程获得单元,用于将所述待拟合方程中的未知系数修改为所确定的取值,获得所述第一方程。
可选的,所述小波系数为Wi,其中,i为小波系数的编号,i为自然数,所述样品检测单元具体设置为:
对每一个定标样品:对该定标样品进行非线性超声检测,获得该定标样品输出的超声信号,对该定标样品输出的超声信号进行小波变换,获得该定标样品的多个频段对应的小波系数,根据公式
β=W2/W1 2
γ=W3/W1 3
计算得到该定标样品的非线性超声系数γ和β,其中,W1为对该定标样品进行非线性超声检测时使用的发射端超声换能器的中心频率F1所在的频段对应的小波系数,W2为频率F2所在的频段对应的小波系数,W3为频率F3所在的频段对应的小波系数,其中,F2=2F1,F3=3F1。
可选的,所述待拟合方程为:
可选的,所述系数拟合单元包括:初次拟合子单元、残差确定子单元和二次拟合子单元,
所述初次拟合子单元,用于根据各定标样品的第一劣化参数的取值以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数k1和k2的取值;
所述残差确定子单元,用于对每个定标样品,根据公式:
所述二次拟合子单元,用于根据各定标样品的残差以及小波系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的fi的取值。
通过上述技术方案,本发明实施例提供的一种利用非线性超声检测物品劣化参数的方法及装置,可以对待测物品进行非线性超声检测,获得所述待测物品输出的超声信号;对所述超声信号进行小波变换,获得多个频段对应的小波系数;根据所述小波系数确定所述待测物品的非线性超声系数;将所述非线性超声系数和所述小波系数作为自变量带入与所述待测物品的材料对应的第一方程中,获得所述待测物品的第一劣化参数的取值。由于本发明使用了非线性超声检测待测物品的劣化参数,因此无需对待测物品进行破坏。同时,本发明对超声信号进行小波变换,具有多分辨率和局部分析的优势,提高了检测的准确性和稳定性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。
图1为本发明实施例提供的一种利用非线性超声检测物品劣化参数的方法的流程图;
图2为本发明实施例提供的另一种利用非线性超声检测物品劣化参数的方法的流程图;
图3为本发明实施例提供的一种利用非线性超声检测物品劣化参数的装置的结构示意图;
图4为本发明实施例提供的另一种利用非线性超声检测物品劣化参数的装置的结构示意图。
具体实施方式
本发明公开了一种利用非线性超声检测物品劣化参数的方法及装置,本领域技术人员可以借鉴本文内容,适当改进工艺参数实现。特别需要指出的是,所有类似的替换和改动对本领域技术人员来说是显而易见的,它们都被视为包括在本发明。本发明的方法及应用已经通过较佳实施例进行了描述,相关人员明显能在不脱离本发明内容、精神和范围内对本文所述的方法和应用进行改动或适当变更与组合,来实现和应用本发明技术。
下面结合实施例,进一步阐述本发明:
如图1所示,本发明实施例提供的一种利用非线性超声检测物品劣化参数的方法,可以包括:
S100、对待测物品进行非线性超声检测,获得所述待测物品输出的超声信号;
其中,待测物品可以为金属材料或其他材料的物品,本发明在此不做限定。
具体的,本发明可以通过非线性超声设备对待测物品进行超声检测,该非线性超声设备包括:发射端超声换能器和接收端超声换能器。其中,发射端超声换能器将非线性超声信号发射到待测物品上,该非线性超声信号在待测物品上进行传导,然后被接收端超声换能器接收。
可选的,发射端超声换能器的中心频率可以位于区间[2.5,10]MHz内。
可选的,接收端超声换能器的中心频率可以为发射端超声换能器的中心频率的两倍。
S200、对所述超声信号进行小波变换,获得多个频段对应的小波系数;
其中,小波变换(Wavelet Transform,WT)是一种变换分析方法,它继承和发展了短时傅立叶变换局部化的思想,同时又克服了窗口大小不随频率变化等缺点,能够提供一个随频率改变的“时间-频率”窗口,是进行信号时频分析和处理的理想工具。它的主要特点是通过变换能够充分突出问题某些方面的特征,能对时间(空间)频率的局部化分析,通过伸缩平移运算对信号(函数)逐步进行多尺度细化,最终达到高频处时间细分,低频处频率细分,能自动适应时频信号分析的要求,从而可聚焦到信号的任意细节,解决了傅立叶变换的困难问题。
传统的傅里叶变换方法无法观测信号的局部频域特性,所以在处理长度有限且非稳态的的信号方面存在一定的不足,容易出现伪信号,影响测量的精度。本发明采用小波分析的方法进行信号的时频转换,具有多分辨率和局部分析的优势,克服了傅里叶变换的不足,对于一定时间区间内的信号进行聚焦,从而得到多分辨率的频域信号,有效地提取非线性超声信号中的二次谐波和三次谐波信号,并且利用这些特征信号与样品的性能劣化参数之间建立关联方程,提高了测量的准确性和稳定性。
具体的,本发明采用的小波变换方法可以具体为mallat算法或者小波包分解算法。其中,本发明可以通过公式
计算小波系数,其中,W(a,b)为小波系数,t为积分变量,f(t)为特征光谱图中的谱线强度数据构成的原函数,ψa,b(t)为小波函数,ψ*为小波函数的共轭函数,a为小波函数的缩放参数,b为小波函数的平移参数。
在实际应用中,本发明可以获得发射端超声换能器的中心频率对应的小波系数,本发明还可以获得该中心频率的二倍频以及三倍频分别对应的小波系数。
S300、根据所述小波系数确定所述待测物品的非线性超声系数;
具体的,本发明可以根据公式
β=W2/W1 2
γ=W3/W1 3
计算得到该待测物品的非线性超声系数γ和β,其中,W1为对该待测物品进行非线性超声检测时使用的发射端超声换能器的中心频率F1所在的频段对应的小波系数,W2为频率F2所在的频段对应的小波系数,W3为频率F3所在的频段对应的小波系数,其中,F2=2F1,F3=3F1。S400、将所述非线性超声系数和所述小波系数作为自变量带入与所述待测物品的材料对应的第一方程中,获得所述待测物品的第一劣化参数的取值,其中,所述第一方程的因变量为所述第一劣化参数的取值且所述第一方程与所述第一劣化参数对应。
所述第一劣化参数可以为:持久强度、蠕变强度、析出物数量、晶粒度或硬度。
其中,第一方程可以为:
其中,C为第一劣化参数,k1、k2和fi为第一方程的系数。
本发明实施例提供的一种利用非线性超声检测物品劣化参数的方法,可以对待测物品进行非线性超声检测,获得所述待测物品输出的超声信号;对所述超声信号进行小波变换,获得多个频段对应的小波系数;根据所述小波系数确定所述待测物品的非线性超声系数;将所述非线性超声系数和所述小波系数作为自变量带入与所述待测物品的材料对应的第一方程中,获得所述待测物品的第一劣化参数的取值。由于本发明使用了非线性超声检测待测物品的劣化参数,因此无需对待测物品进行破坏。同时,本发明对超声信号进行小波变换,具有多分辨率和局部分析的优势,提高了检测的准确性和稳定性。
在实际应用中,本发明可以通过对定标样品检测来确定该第一方程的系数。具体的,在第一方程的系数未知时,该第一方程为待拟合方程。如图2所示,本发明实施例提供的另一种利用非线性超声检测物品劣化参数的方法,在步骤S100前,还可以包括:
S001、获得材料与所述待测物品相同的多个定标样品,其中,各定标样品的第一劣化参数的取值已知;
在实际应用中,本发明可以分别针对每种材料进行劣化参数对应的待拟合方程中未知系数的确定。
S002、对每一个定标样品:对该定标样品进行非线性超声检测,获得该定标样品输出的超声信号,对该定标样品输出的超声信号进行小波变换,获得该定标样品的多个频段对应的小波系数,根据该定标样品的多个频段对应的小波系数确定该定标样品的非线性超声系数;
其中,步骤S002进行非线性超声检测时所使用的非线性超声的频率可以和步骤S100非线性超声检测时所使用的非线性超声的频率相同。
进一步,步骤S002进行小波变换后获得小波系数对应的各频段可以和步骤S200进行小波变换后获得小波系数对应的各频段。
例如:步骤S002进行小波变换后获得小波系数对应的各频段分别为:(0,3]MHz、(3,6]MHz和(6,9]MHz,则步骤S200进行小波变换后获得小波系数对应的各频段也分别为:(0,3]MHz、(3,6]MHz和(6,9]MHz。
其中,所述小波系数可以为Wi,其中,i为小波系数的编号,i为自然数,根据该定标样品的多个频段对应的小波系数确定该定标样品的非线性超声系数的具体过程可以包括:
根据公式
β=W2/W1 2
γ=W3/W1 3
计算得到该定标样品的非线性超声系数γ和β,其中,W1为对该定标样品进行非线性超声检测时使用的发射端超声换能器的中心频率F1所在的频段对应的小波系数,W2为频率F2所在的频段对应的小波系数,W3为频率F3所在的频段对应的小波系数,其中,F2=2F1,F3=3F1。
其中,本发明实施例进行非线性超声检测时使用的发射端超声换能器的中心频率可以位于区间[2.5,10]MHz内。
S003、根据各定标样品的第一劣化参数的取值、小波系数以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数的取值,其中,所述待拟合方程为第一劣化参数与非线性超声系数的关联方程;
其中,待拟合方程可以为:
具体的,数据拟合的方式可以有多种,如最小二乘拟合。
在实际应用中,步骤S003可以具体包括:
根据各定标样品的第一劣化参数的取值以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数k1和k2的取值;
对每个定标样品,根据公式:
根据各定标样品的残差以及小波系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的fi的取值。
其中,本发明在确定k1和k2的取值时可以首先忽略残差,即待拟合方程变为:C=k1β+k2γ,然后根据各定标样品的第一劣化参数的取值以及非线性超声系数进行数据拟合确定k1和k2的取值。在确定k1和k2的取值后,本发明就可以确定各定标样品的残差,进而对残差进行二次拟合,从而确定fi的取值。
S004、将所述待拟合方程中的未知系数修改为所确定的取值,获得所述第一方程。
为方便理解,下面举例说明:
假设有n个材料与所述待测物品相同的定标样品,定标样品的编号分别为1至n,各定标样品的第一劣化参数的取值分别为C1至Cn。
对每一个定标样品:对该定标样品进行非线性超声检测,获得该定标样品输出的超声信号,对该定标样品输出的超声信号进行小波变换,获得该定标样品的多个频段对应的小波系数Wi,其中,i为小波系数的编号,该编号与定标样品的编号相同。具体的,每个定标样品的Wi可以为多个,分别与不同的频段对应。根据该定标样品的多个频段对应的小波系数Wi确定该定标样品的非线性超声系数βi和γi。
至此,可以得到材料与所述待测物品相同的n个定标样品的第一劣化参数的取值:C1至Cn,也可以得到上述n个定标样品中每个定标样品的小波系数W1至Wm,也可以得到上述n个定标样品的非线性超声系数:β1至βn,以及γ1至γn。其中,m为每个定标样品的小波系数的个数。
这样,通过上述C1至Cn、β1至βn、γ1至γn就可以进行数据拟合,得到公式
C=k1β+k2γ
中的未知系数k1和k2的取值。
但上述公式忽略了残差,使用A代表残差,则可以通过A=C-k1β-k2γ计算得到各定标样品的残差,例如:A1=C1-k1β1-k2γ1。因此,本发明可以获得各定标样品的残差:A1至An。
这样,本发明就可以确定待拟合方程
中各未知系数的取值,从而确定第一劣化参数与非线性超声系数、小波系数的关联方程。在获得该关联方程后,就可以通过图1所示方法获得待测物品的第一劣化参数的取值。
与上述方法相对应,本发明实施例还提供了一种利用非线性超声检测物品劣化参数的装置。
如图3所示,本发明实施例提供的一种利用非线性超声检测物品劣化参数的装置,可以包括:信号获得单元100、小波系数获得单元200、超声系数获得单元300和劣化参数获得单元400,
所述信号获得单元100,用于对待测物品进行非线性超声检测,获得所述待测物品输出的超声信号;
其中,待测物品可以为金属材料或其他材料的物品,本发明在此不做限定。
具体的,本发明可以通过非线性超声设备对待测物品进行超声检测,该非线性超声设备包括:发射端超声换能器和接收端超声换能器。其中,发射端超声换能器将非线性超声信号发射到待测物品上,该非线性超声信号在待测物品上进行传导,然后被接收端超声换能器接收。
可选的,发射端超声换能器的中心频率可以位于区间[2.5,10]MHz内。
可选的,接收端超声换能器的中心频率可以为发射端超声换能器的中心频率的两倍。
所述小波系数获得单元200,用于对所述超声信号进行小波变换,获得多个频段对应的小波系数;
其中,小波变换(Wavelet Transform,WT)是一种变换分析方法,它继承和发展了短时傅立叶变换局部化的思想,同时又克服了窗口大小不随频率变化等缺点,能够提供一个随频率改变的“时间-频率”窗口,是进行信号时频分析和处理的理想工具。它的主要特点是通过变换能够充分突出问题某些方面的特征,能对时间(空间)频率的局部化分析,通过伸缩平移运算对信号(函数)逐步进行多尺度细化,最终达到高频处时间细分,低频处频率细分,能自动适应时频信号分析的要求,从而可聚焦到信号的任意细节,解决了傅立叶变换的困难问题。
传统的傅里叶变换方法无法观测信号的局部频域特性,所以在处理长度有限且非稳态的的信号方面存在一定的不足,容易出现伪信号,影响测量的精度。本发明采用小波分析的方法进行信号的时频转换,具有多分辨率和局部分析的优势,克服了傅里叶变换的不足,对于一定时间区间内的信号进行聚焦,从而得到多分辨率的频域信号,有效地提取非线性超声信号中的二次谐波和三次谐波信号,并且利用这些特征信号与样品的性能劣化参数之间建立关联方程,提高了测量的准确性和稳定性。
具体的,本发明采用的小波变换方法可以具体为mallat算法或者小波包分解算法。其中,本发明可以通过公式
计算小波系数,其中,W(a,b)为小波系数,t为积分变量,f(t)为特征光谱图中的谱线强度数据构成的原函数,ψa,b(t)为小波函数,ψ*为小波函数的共轭函数,a为小波函数的缩放参数,b为小波函数的平移参数。
在实际应用中,本发明可以获得发射端超声换能器的中心频率对应的小波系数,本发明还可以获得该中心频率的二倍频以及三倍频分别对应的小波系数。
所述超声系数获得单元300,用于根据所述小波系数确定所述待测物品的非线性超声系数;
具体的,超声系数获得单元300可以设置为根据公式
β=W2/W1 2
γ=W3/W1 3
计算得到该待测物品的非线性超声系数γ和β,其中,W1为对该待测物品进行非线性超声检测时使用的发射端超声换能器的中心频率F1所在的频段对应的小波系数,W2为频率F2所在的频段对应的小波系数,W3为频率F3所在的频段对应的小波系数,其中,F2=2F1,F3=3F1。所述劣化参数获得单元400,用于将所述非线性超声系数和所述小波系数作为自变量带入与所述待测物品的材料对应的第一方程中,获得所述待测物品的第一劣化参数的取值,其中,所述第一方程的因变量为所述第一劣化参数的取值且所述第一方程与所述第一劣化参数对应。
所述第一劣化参数可以为:持久强度、蠕变强度、析出物数量、晶粒度或硬度。
其中,第一方程可以为:
其中,C为第一劣化参数,k1和k2为第一方程的系数。
本发明实施例提供的一种利用非线性超声检测物品劣化参数的装置,可以对待测物品进行非线性超声检测,获得所述待测物品输出的超声信号;对所述超声信号进行小波变换,获得多个频段对应的小波系数;根据所述小波系数确定所述待测物品的非线性超声系数;将所述非线性超声系数和所述小波系数作为自变量带入与所述待测物品的材料对应的第一方程中,获得所述待测物品的第一劣化参数的取值。由于本发明使用了非线性超声检测待测物品的劣化参数,因此无需对待测物品进行破坏。同时,本发明对超声信号进行小波变换,具有多分辨率和局部分析的优势,提高了检测的准确性和稳定性。
在实际应用中,本发明可以通过对定标样品检测来确定该第一方程的系数。具体的,在第一方程的系数未知时,该第一方程为待拟合方程。如图4所示,在图3所示装置基础上,本发明实施例提供的另一种利用非线性超声检测物品劣化参数的装置还可以包括:样品获得单元001、样品检测单元002、系数拟合单元003和方程获得单元004,
所述样品获得单元001,用于在所述信号获得单元对待测物品进行非线性超声检测前,获得材料与所述待测物品相同的多个定标样品,其中,各定标样品的第一劣化参数的取值已知;
在实际应用中,本发明可以分别针对每种材料进行劣化参数对应的待拟合方程中未知系数的确定。
所述样品检测单元002,用于对每一个定标样品:对该定标样品进行非线性超声检测,获得该定标样品输出的超声信号,对该定标样品输出的超声信号进行小波变换,获得该定标样品的多个频段对应的小波系数,根据该定标样品的多个频段对应的小波系数确定该定标样品的非线性超声系数;
其中,样品检测单元002进行非线性超声检测时所使用的非线性超声的频率可以和信号获得单元100进行非线性超声检测时所使用的非线性超声的频率相同。
进一步,样品检测单元002进行小波变换后获得小波系数对应的各频段可以和小波系数获得单元200进行小波变换后获得小波系数对应的各频段相同。
例如:步骤S002进行小波变换后获得小波系数对应的各频段分别为:(0,3]MHz、(3,6]MHz和(6,9]MHz,则步骤S200进行小波变换后获得小波系数对应的各频段也分别为:(0,3]MHz、(3,6]MHz和(6,9]MHz。
其中,所述小波系数为Wi,其中,i为小波系数的编号,i为自然数,所述样品检测单元002可以具体设置为:
对每一个定标样品:对该定标样品进行非线性超声检测,获得该定标样品输出的超声信号,对该定标样品输出的超声信号进行小波变换,获得该定标样品的多个频段对应的小波系数,根据公式
β=W2/W1 2
γ=W3/W1 3
计算得到该定标样品的非线性超声系数γ和β,其中,W1为对该定标样品进行非线性超声检测时使用的发射端超声换能器的中心频率F1所在的频段对应的小波系数,W2为频率F2所在的频段对应的小波系数,W3为频率F3所在的频段对应的小波系数,其中,F2=2F1,F3=3F1。其中,本发明实施例进行非线性超声检测时使用的发射端超声换能器的中心频率可以位于区间[2.5,10]MHz内。
所述系数拟合单元003,用于根据各定标样品的第一劣化参数的取值、小波系数以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数的取值,其中,所述待拟合方程为第一劣化参数与非线性超声系数、小波系数的关联方程;
其中,所述系数拟合单元003可以包括:初次拟合子单元、残差确定子单元和二次拟合子单元,
所述初次拟合子单元,用于根据各定标样品的第一劣化参数的取值以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数k1和k2的取值;
所述残差确定子单元,用于对每个定标样品,根据公式:
所述二次拟合子单元,用于根据各定标样品的残差以及小波系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的fi的取值。
其中,本发明在确定k1和k2的取值时可以首先忽略残差,即待拟合方程变为:C=k1β+k2γ,然后根据各定标样品的第一劣化参数的取值以及非线性超声系数进行数据拟合确定k1和k2的取值。在确定k1和k2的取值后,本发明就可以确定各定标样品的残差,进而对残差进行二次拟合,从而确定fi的取值。
所述方程获得单元004,用于将所述待拟合方程中的未知系数修改为所确定的取值,获得所述第一方程。
其中,待拟合方程可以为:
具体的,数据拟合的方式可以有多种,如最小二乘拟合等。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中的各个实施例均采用相关的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于***实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
Claims (4)
1.一种利用非线性超声检测物品劣化参数的方法,其特征在于,包括:
对待测物品进行非线性超声检测,获得所述待测物品输出的超声信号;
对所述超声信号进行小波变换,获得N个频段对应的小波系数;其中,N≥3且N为正整数,所述N个频段至少包括中心频率的频段、二倍频率的频段和三倍频率的频段;
根据所述小波系数确定所述待测物品的非线性超声系数γ和β;其中,β为根据二倍频率计算得到的系数,γ为根据三倍频率计算得到的系数;
将所述非线性超声系数和所述小波系数作为自变量带入与所述待测物品的材料对应的第一方程中,获得所述待测物品的第一劣化参数的取值,其中,所述第一方程的因变量为所述第一劣化参数的取值且所述第一方程与所述第一劣化参数对应;
其中,所述第一方程的获取过程包括:
获得材料与所述待测物品相同的多个定标样品,其中,各定标样品的第一劣化参数的取值已知;
对每一个定标样品:对该定标样品进行非线性超声检测,获得该定标样品输出的超声信号,对该定标样品输出的超声信号进行小波变换,获得该定标样品的多个频段对应的小波系数,根据该定标样品的多个频段对应的小波系数确定该定标样品的非线性超声系数;
根据各定标样品的第一劣化参数的取值以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数k1和k2的取值;
根据各定标样品的残差以及小波系数进行数据拟合,确定与所述定标样品的材料对应
的待拟合方程中的fi的取值;其中,所述待拟合方程为:,C为第
一劣化参数,k1、k2和fi为所述待拟合方程的未知系数,m为小波系数的个数,为
残差;其中,Wi为小波系数,i为小波系数的编号,i为自然数;fi为与Wi对应的残差拟合系数;
将所述待拟合方程中的未知系数修改为所确定的取值,获得所述第一方程。
2.根据权利要求1所述的方法,其特征在于,所述根据该定标样品的多个频段对应的小波系数确定该定标样品的非线性超声系数γ和β,包括:
根据公式
β= W2/ W1 2
γ= W3/ W1 3
计算得到该定标样品的非线性超声系数γ和β,其中,W1为对该定标样品进行非线性超声检测时使用的发射端超声换能器的中心频率F1所在的频段对应的小波系数,W2为频率F2所在的频段对应的小波系数,W3为频率F3所在的频段对应的小波系数,其中,F2=2F1,F3=3F1。
3.一种利用非线性超声检测物品劣化参数的装置,其特征在于,包括:样品获得单元、样品检测单元、系数拟合单元、方程获得单元、信号获得单元、小波系数获得单元、超声系数获得单元和劣化参数获得单元,
所述信号获得单元,用于对待测物品进行非线性超声检测,获得所述待测物品输出的超声信号;
所述小波系数获得单元,用于对所述超声信号进行小波变换,获得多个频段对应的小波系数;
所述超声系数获得单元,用于根据所述小波系数确定所述待测物品的非线性超声系数;
所述劣化参数获得单元,用于将所述非线性超声系数和所述小波系数作为自变量带入与所述待测物品的材料对应的第一方程中,获得所述待测物品的第一劣化参数的取值,其中,所述第一方程的因变量为所述第一劣化参数的取值且所述第一方程与所述第一劣化参数对应;
所述样品获得单元,用于在所述信号获得单元对待测物品进行非线性超声检测前,获得材料与所述待测物品相同的多个定标样品,其中,各定标样品的第一劣化参数的取值已知;
所述样品检测单元,用于对每一个定标样品:对该定标样品进行非线性超声检测,获得该定标样品输出的超声信号,对该定标样品输出的超声信号进行小波变换,获得该定标样品的多个频段对应的小波系数,根据该定标样品的多个频段对应的小波系数确定该定标样品的非线性超声系数;
所述系数拟合单元,用于根据各定标样品的第一劣化参数的取值、小波系数以及非线
性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数的
取值,其中,所述待拟合方程为:,C为第一劣化参数,k1、k2和fi为所
述待拟合方程的未知系数,m为小波系数的个数,为残差;其中,Wi为小波系数,i为
小波系数的编号,i为自然数;fi为与Wi对应的残差拟合系数;
所述方程获得单元,用于将所述待拟合方程中的未知系数修改为所确定的取值,获得所述第一方程;
其中,所述系数拟合单元包括:初次拟合子单元、残差确定子单元和二次拟合子单元,
所述初次拟合子单元,用于根据各定标样品的第一劣化参数的取值以及非线性超声系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的未知系数k1和k2的取值;
所述二次拟合子单元,用于根据各定标样品的残差以及小波系数进行数据拟合,确定与所述定标样品的材料对应的待拟合方程中的fi的取值。
4.根据权利要求3所述的装置,其特征在于,所述样品检测单元具体设置为:
对每一个定标样品:对该定标样品进行非线性超声检测,获得该定标样品输出的超声信号,对该定标样品输出的超声信号进行小波变换,获得该定标样品的多个频段对应的小波系数,根据公式
β= W2/ W1 2
γ= W3/ W1 3
计算得到该定标样品的非线性超声系数γ和β,其中,W1为对该定标样品进行非线性超声检测时使用的发射端超声换能器的中心频率F1所在的频段对应的小波系数,W2为频率F2所在的频段对应的小波系数,W3为频率F3所在的频段对应的小波系数,其中,F2=2F1,F3=3F1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810502786.5A CN108918678B (zh) | 2018-05-23 | 2018-05-23 | 利用非线性超声检测物品劣化参数的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810502786.5A CN108918678B (zh) | 2018-05-23 | 2018-05-23 | 利用非线性超声检测物品劣化参数的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108918678A CN108918678A (zh) | 2018-11-30 |
CN108918678B true CN108918678B (zh) | 2021-01-29 |
Family
ID=64402746
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810502786.5A Active CN108918678B (zh) | 2018-05-23 | 2018-05-23 | 利用非线性超声检测物品劣化参数的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108918678B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111141822A (zh) * | 2020-01-06 | 2020-05-12 | 大连理工大学 | 一种基于非线性扭波的hp型炉管高温组织劣化无损评价方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104458915A (zh) * | 2013-09-22 | 2015-03-25 | 南京大学 | 一种风塔焊缝非线性超声检测方法 |
CN104764805A (zh) * | 2015-03-19 | 2015-07-08 | 北京工业大学 | 一种用于有机材料热氧老化评价的非线性超声检测方法 |
-
2018
- 2018-05-23 CN CN201810502786.5A patent/CN108918678B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104458915A (zh) * | 2013-09-22 | 2015-03-25 | 南京大学 | 一种风塔焊缝非线性超声检测方法 |
CN104764805A (zh) * | 2015-03-19 | 2015-07-08 | 北京工业大学 | 一种用于有机材料热氧老化评价的非线性超声检测方法 |
Non-Patent Citations (1)
Title |
---|
Wavelet based harmonics decomposition of ultrasonic signal in assessment of plastic strain in aluminum;Amir Mostavi et al.;《Measurement》;20171231;第106卷;第66-78页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108918678A (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103837492B (zh) | 一种基于近红外光谱技术的猕猴桃膨大果无损检测方法 | |
Song et al. | Wavelet based noise suppression technique and its application to ultrasonic flaw detection | |
CN103134767B (zh) | 一种红外光谱校正鉴定白酒品质的方法 | |
WO2007013068A3 (en) | System and methods for non-destructive testing of tubular systems | |
CN108918678B (zh) | 利用非线性超声检测物品劣化参数的方法及装置 | |
CN109340586B (zh) | 一种供水管道泄露的检测方法及*** | |
CN108508091B (zh) | 检测物品劣化参数的方法及装置 | |
CN100405384C (zh) | 粗晶材料超声检测时频分析处理方法 | |
CN103926324B (zh) | 一种超声表面波检测主汽管道蠕变损伤的方法 | |
CN114660105A (zh) | 土壤重金属Cd含量反演方法、***、介质、计算机设备 | |
CN107528646B (zh) | 一种基于宽带频谱的干扰信号识别及提取方法 | |
CN106126879B (zh) | 一种基于稀疏表示技术的土壤近红外光谱分析预测方法 | |
CN106092961A (zh) | 一种快速检测农产品中化学成分的方法 | |
WO2007132509A1 (ja) | 共振現象を利用した超音波探査方法 | |
Chen et al. | Noise-robust adaptive feature mode decomposition method for accurate feature extraction in rotating machinery fault diagnosis | |
US5696691A (en) | Self-adjusting statistical noise analyzer with interference suppression | |
CN105891148A (zh) | 利用近红外光谱测定废水厌氧处理过程中vfa浓度的方法 | |
CN108107111A (zh) | 一种耐热钢部件非线性超声波检测方法 | |
CN105954206B (zh) | 紫玉米叶片花色苷含量的测量方法及*** | |
CN211477587U (zh) | 在线安全预警装置和故障诊断*** | |
CN111010193B (zh) | 基于压电导纳的数据压缩重构及结构损伤识别方法、设备 | |
CN110532635A (zh) | 一种基于时域的管道泄漏检测算法 | |
US20160087660A1 (en) | Method and system for extending dynamic range of receiver by compensating for non-linear distortion | |
TW202326155A (zh) | 一種用於量測待測物的複介電係數的裝置、以及用於複合介電質的時域多重反射訊號的量測裝置及其量測方法 | |
CN110596249A (zh) | 燃气调压器的故障检测方法及装置 |
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. 2096, Beijing New Energy Technology Research Institute Patentee after: Guoneng boiler and Pressure Vessel Inspection Co.,Ltd. Address before: No. 2096, Beijing New Energy Technology Research Institute Patentee before: GUODIAN BOILER AND PRESSURE VESSEL INSPECTION Co.,Ltd. |
|
CP01 | Change in the name or title of a patent holder |