亚星游戏官网-www.yaxin868.com

专利名称：一种超声波钢管混凝土柱质量检测方法
技术领域：
本发明涉及一种超声波钢管混凝土柱质量检测方法及检测系统，尤其涉及一种经 过滤波处理和HHT变换的超声波钢管混凝土柱质量检测方法及检测系统。
背景技术：
传统的超声检测法主要有首波声时法、首波幅值法、波形识别法和首波频率法 等。工程中经常以首波声时法为主，同时综合考虑首波幅值、首波波形和首波频率来对钢管 混凝土的浇筑质量做出评价。实际检测结果表明，在对钢管混凝土进行质量检测时，超声检 测测得的超声波信号经常会由于含有噪声信号和波形畸变导致首波声时的难以辨读和首 波幅值、波形、频率的误判，最终导致钢管混凝土柱的损伤识别效果较差。

发明内容
本发明解决的技术问题是提供一种超声波钢管混凝土柱质量检测方法及超声波 钢管混凝土柱质量检测系统，克服现有技术的首波声时难以判读，波形受到噪声信号影响 以致特征不明显的缺点。本发明的技术方案是提供一种超声波钢管混凝土柱质量检测方法，包括检测钢 管混凝土柱质量的超声波检测仪，将超声波检测仪的超声波发生器和超声波接收器分别设 置在所述预检测钢管混凝土柱的两侧，超声波钢管混凝土柱质量检测方法包括如下步骤发射和接收超声波所述超声波发生器发射超声波，所述超声波接收器接收所述 超声波发生器发射的经钢管混凝土柱的超声波；滤波处理对接收的超声波信号进行滤波处理；HHT变换对经过滤波处理的超声波信号进行HHT (Hilbert HuangTranslate,简 称“HHT”)变换，即(经验模态分解与Hilbert谱分析)；质量检测根据HHT变换后的超声波信号进行钢管混凝土柱质量检测。本发明的进一步技术方案是在滤波处理步骤中，包括在滤波前先对超声波信号 末尾进行补零延拓。本发明的进一步技术方案是在滤波处理步骤中，补零长度设置为相位延迟时间 的长度。本发明的进一步技术方案是在滤波处理步骤中，还包括在滤波完成之后再在超 声波信号前端截掉相应的数据长度。本发明的进一步技术方案是在质量检测步骤中，包括钢管混凝土脱层位置的确定。本发明的进一步技术方案是在质量检测步骤中，包括钢管混凝土内含杂物位置 的确定。本发明的进一步技术方案是在质量检测步骤中，包括钢管混凝土中心孔洞位置 的确定。
本发明的技术方案是构建一种超声波钢管混凝土柱质量检测系统，包括检测钢 管混凝土柱质量的超声波检测仪，将超声波检测仪的超声波发生器和超声波接收器分别设 置在所述预检测钢管混凝土柱的两侧，还包括将所述超声波接收器接收的超声波进行滤波 处理的滤波单元，对所述滤波单元滤波后的超声波信号进行HHT变换的HHT变换单元，根据 HHT变换后的超声波信号检测钢管混凝土柱质量的质量检测单元。本发明的进一步技术方案是还包括对所述超声波接收器接收的超声波信号进行 预处理的预处理单元，所述预处理单元在滤波前先对超声波信号末尾进行补零延拓，在滤 波完成之后再在超声波信号前端截掉相应的数据长度。本发明的进一步技术方案是补零长度设置为相位延迟时间的长度。本发明的技术效果是本发明一种超声波钢管混凝土柱质量检测方法及检测系 统，通过对所述超声波接收器接收所述超声波发生器发射的经钢管混凝土柱的超声波进行 滤波处理后再进行HHT变换后，根据HHT变换后的超声波信号进行钢管混凝土柱质量评估。 本方法可靠性强，检测准确，方便实用。


图1为本发明的流程图。
图2为本发明滤波前的超声波波形图。
图3为本发明滤波后的超声波波形图。
图4为本发明滤波前的超声波Hilbert能量谱图。
图5为本发明滤波后的超声波Hilbert能量谱图。
图6为本发明滤波前波形图对比。
图7为本发明无损与脱层的一阶固有模态波形对比图。
图8为本发明无缺陷位置的能量谱图。
图9为本发明脱层缺陷位置的能量谱图。
图10为本发明无损与含杂物的一阶固有模态波形对比图。
图11为本发明含杂物缺陷位置的能量谱图。
图12为本发明无损与内含孔洞的一阶固有模态波形对比图
图13为本发明内含孔洞缺陷位置的能量谱图。
图14本为本发明结构示意图。
具体实施例方式下面结合具体实施例，对本发明技术方案进一步说明。如图1所示，本发明的具体实施方式
是本发明提供一种超声波钢管混凝土柱质 量检测方法，包括检测钢管混凝土柱质量的超声波检测仪，将超声波检测仪的超声波发生 器和超声波接收器分别设置在所述预检测钢管混凝土柱的两侧，超声波钢管混凝土柱质量 检测方法包括如下步骤步骤100 发射和接收超声波，即，所述超声波发生器发射超声波，所述超声波接 收器接收所述超声波发生器发射的经钢管混凝土柱的超声波。步骤200 滤波处理，即，对接收的超声波信号进行滤波处理。具体处理过程如下
4超声波检测仪的原始信号为45KHz的窄带信号，超声波接收器接收到的超声波信号中混入 了低频和高频的噪声信号的干扰，本发明中具体实施例中，采用基于matlab软件的kaiser 窗带通数字滤波器来将超声信号中的噪声信号过滤掉，将带通数字滤波器的上下限频率设 为IOKHz和80KHz。在进行信号的滤波时，滤波得到的信号将会产生整体相延迟，同时导致 原始信号末尾的数据信息丢失，为了消除这一影响，在滤波前先对信号末尾进行补零延拓， 补零长度设置为相位延迟时间的长度，在滤波完成之后再在信号前端截掉相应的数据长 度，最终得到滤波前后的波形图及其相应的Hilbert能量谱图，图2为滤波前的超声波波形 图，图3为滤波后的超声波波形图，图4为滤波前的超声波Hilbert能量谱图，图5为滤波 后的超声波Hi Ibert能量谱图。由此可以看到，经过滤波后的波形明显变得光滑，不再有毛 刺，方便了首波声时的判读；同时由滤波前后的Hilbert能量谱图对比可以看出，带通数字 滤波器很好的过滤掉了低频噪声信号。步骤300 =HHT变换，即，对经过滤波处理的超声波信号进行HHT变换，即(经验模 态分解与Hilbert谱分析)。具体过程如下一、HHT变换方法。黄氏变换(Hilbert Huang Translate，简称HHT)新方法，即经验模态分解与 Hilbert谱分析方法。这一方法创造性地提出了固有模态函数(IntrinsicMode Function, 简称IMF)的新概念以及将任意信号分解为固有模态函数方法——经验模态分解法(the Empirical Mode Decomposition，简称EMD)，这样赋予了瞬时频率合理的定义、物理意义和 求法，初步建立了以瞬时频率为表征信号交变的基本量，以固有模态函数为基本时域信号 的新时频分析方法体系。作为一个崭新的时频分析方法，HHT完全独立于傅立叶变换。此 方法的关键部分是经验模态分解，经验模态分解能够进行非线性、非平稳信号线性化和平 稳化处理，保留信号本身特性。任何复杂的数据都可以在固有时间范围内被直接抽取能量， 从而产生一组固有模态函数集合，使它具有很好的Hilbert变换特性。随后通过Hilbert 变换，固有模态函数产生出瞬时频率和瞬时振幅，瞬时振幅在频率_时间平面上的分布就 是Hilbert谱，表现为能量-频率-时间三者之间的分布关系。这种分解同样可以看作是 固有模态函数的数据扩展。于是固有模态函数可以作为表述线性或非线性数据的基�。�它 是完整的并且几乎是正交的，最重要的是它具有很好的适应性。二、经验模态分解法。经验模态分解法的基本思路是用波动上、下包络的平均值去确定“瞬时平衡位 置”，进而提取出固有模态函数。上、下包络线是由三次样条函数对极大值点和极小值点进 行拟合得到的。经验模式分解过程的基本过程可概括如下(1)寻找信号x(t)所有局部极大值和局部极小值，为更好保留原序列的特性，局 部极大值定义为时间序列中的某个时刻的值，它只要满足既大于前一时刻的值也大于后一 时刻的值即可。局部极小值的提取同理，即该时刻的值满足既小于前一时刻的值也小于后 一时刻的值。使用三次样条函数进行拟合，获得上包络线xmaX(t)和下包络线Xmin(t)；(2)计算上、下包络线的均值 m(t) = [xfflax (t)+Xfflin (t)]/2 ；(3)用原信号x(t)减去均值m(t)，得到第一个组件h(t) =x(t)-m(t)；由于原始序列的差异，组件h(t)不一定就是一个IMF，如果h(t)不满足固有模态函数两个条件 整个数据范围内，极值点和过零点的数量相等或者相差一个；在任意点处，所有极大值点形 成的包络线和所有极小值点形成的包络线的平均值为零。就把h(t)当成原始信号，重复 (1)-(3)，直到满足条件为止，这时满足固有模态函数条件的h(t)作为一个IMF，令I1 (t)= h(t)，至此第一个IMF已经成功的提取了。由于剩余的r(t) =x (O-I1U)仍然包含具有更 长周期组件的信息，因此可以把它看成新的信号，重复上述过程，依次得到第二个12 (t)，第 三个I3(t)，…，当r(t)满足单调序列或常值序列条件时，终止筛选过程，可以认为完成了 提取固有模态函数的任务，最后的r(t)称为余项，它是原始信号的趋势项。由此可得x(t)
η
的表达式4/) = 1/;(0 + 〃(0即原始序列是由η个IMF与一个趋势项组成。如上所述，整个 ‘=1 ,
过程就像筛选过程，根据时间特性把固有模态函数从信号中提取出来。三、超声波信号的HHT变换。频率是个极其重要的物理量，定义为信号周期倒数，其物理含义显而易见。对于正 弦信号，它的频率为恒值。但是对于大部分信号，它的频率是随时间变化的函数，故提出瞬 时频率概念。瞬时频率即表征信号在局部时间点上瞬态频率特性，整个持续期上的瞬时频 率反映了信号频率的时变规律。
对于随机时间序列X⑴，对其进行Hilbert变换，可以得到Y(t)如下 汴)=丄沙(厂^1^)(ι)
其中，PV为柯西主值(Cauchy principal value)。该式表示Y (t)是X(t)与丄
的卷积。通过这个定义，X(t)和Y(t)组成了一个共轭复数对，于是可以得到一个解析信号 Z(t)如下Z (t) =X (t) +iY (t) = a (t) e1 θ (t)其中，
'a{t) = [x\t) + Y2(t)\
⑵
π1/2
θ{ ) = arctan
)、
Z^-S
3
κιν其中，a(t)为解析信号的幅值，θ (t)为解析信号的相位角。从理论上讲，虚部的定义方法有很多种。但是Hilbert变换为其提供了一个唯一 的虚部值，这就使得其结果成为一个解析函数。得到了相位，就可以得到瞬时频率，因为瞬 时频率就是相位导数。 对固有模态函数进行Hilbert变换后，可以用下面方式表示信号 由于rn信号中的能量很�。�此处将其忽略。由希尔伯特变换得出的振幅和频率都是时间的函数，如果用三维图形表达幅值、频率和时间之间的关系，或者把振幅用灰度的形式显示在频率-时间平面上，就可以得到 Hilbert谱H(w，t)，H(w, t)表示振幅，它是时间与频率的函数。如果将H (w，t)对时间积分，则得到Hilbert边际谱h (w)h(w) =S TH(w, t)dt(6)边际谱提供了对每个频率的总振幅的量测，表达了整个时间长度内累积的振 幅。另外，作为希尔伯特边际谱的附加结果，可以得到如式(7)定义的Hilbert瞬时能量 IE⑴IE (t) = / wH2(w, t)dw (7)瞬时能量提供了信号能量随时间的变换情况。将振幅的平方对时间积分，可以得 到 Hilbert 能量谱 ES (w)ES (w) = / tH2(w, t)dt (8)Hilbert能量谱提供了对于每个频率的能量的量测，表达了每个频率在整个时间 长度内所累积的能量。步骤400 质量检测，即，根据HHT变换后的超声波信号进行钢管混凝土柱质量检 测。具体过程如下以下实施例中，钢管混凝土柱中钢管直径为219mm，取超声波在混凝土中的传播速 度为4400m/s，以面波形式在钢管中传播的速度为2900m/s。则可算得超声波穿过核心混凝 土到达接收端的时间、=0. 219 + 4400 = 49. 8 μ s,以面波形式沿钢管传播到达接收端的 时间 t2 = 0. 219X π +2 + 2900 = 118. 6 μ s。一、钢管混凝土脱层位置的确定。图6为滤波前波形图对比，图7为无损与脱层的一阶固有模态波形对比图。对比 图6和图7可以看到结合滤波处理的固有模态分解结果很好的抑制了噪声的影响，分析 结果表明，在无缺陷位置，有较多能量的超声波直接穿过核心混凝土到达了接收端，可看 到在50 μ s左右开始有明显波形；在存在脱层缺陷的位置，只有很少的能量的超声波直接 穿过空气层入射到混凝土中，绝大多数能量都以面波形式绕钢管传播到达接收端，以至在 120 μ s左右处才开始有明显波形。图8为无缺陷位置的能量谱图，图9为脱层缺陷位置的能量谱图。对比图8和图9 可以看到脱层缺陷位置的能量谱图的能量幅值明显要比无缺陷位置的要高，且高频区的能 量幅值也比无缺陷位置的高频区能量幅值要大，这是因为大部分能量沿着钢管传播，衰减 得较少。二、钢管混凝土内含杂物位置的确定。对比图10中的两条波形曲线可以看到，无缺陷位置的超声波信号的在50μ s处开 始有明显的波形，而内含杂物位置的超声波信号则一直到80 μ s左右处才开始出现明显的 波形。对比图8和图11可以看到含杂物缺陷位置处的能量谱图的峰值明显要比无缺陷 位置处的能量谱图的峰值�。�这是因为超声波在杂物与混凝土的边界面上有一部分能量由 于反射而损失掉，没有直接传递到接收端。三、钢管混凝土中心孔洞位置的确定。
在中心孔洞直径约1. 5cm位置处采用对测法测得超声波信号，其一阶固有模态波 形如图12中虚线所示。对比图12中的两条波形曲线可以看到，内含孔洞位置的超声波信 号的一阶固有模态波形的首波到达时间比无损位置的稍微晚了一点。此外，可以看到两条 曲线之间在整个采样时间范围内都存在明显的相位差，内含孔洞位置的一阶固有模态波形 总比无损位置的要先到达波峰位置。图13为内含孔洞缺陷位置的能量谱图，对比图8和图13可见，无缺陷位置的能量 谱图与内含孔洞缺陷位置的能量谱图两者之间在形状和能量分布上差别不大，可见当孔洞 直径较小时，孔洞的存在对超声波信号的能量谱影响较�。�可能无法识别出区别。若孔洞较 大，则无缺陷位置的能量谱图与内含孔洞缺陷位置的能量谱图两者之间在形状和能量分布 上差别较大。本发明一种超声波钢管混凝土柱质量检测方法，通过对所述超声波接收器接收所 述超声波发生器发射的经钢管混凝土柱的超声波进行滤波处理后再进行HHT变换后，根据 HHT变换后的超声波信号进行钢管混凝土柱质量检测。本方法可靠性强，检测准确，方便实 用。如图14所示，本发明的
具体实施例方式构建一种超声波钢管混凝土柱质量检测系 统，包括检测钢管混凝土柱质量的超声波检测仪，将超声波检测仪的超声波发生器2和超 声波接收器3分别设置在所述预检测钢管混凝土柱1的两端，还包括将所述超声波接收器3 接收的超声波进行滤波处理的滤波单元5，对所述滤波单元5滤波后的超声波信号进行HHT 变换的HHT变换单元6，根据HHT变换后的超声波信号检测钢管混凝土柱1质量的质量检测 单元7。本发明的HHT变换单元6的具体处理过程见超声波钢管混凝土柱质量检测方法的 步骤300，S卩，HHT变换步骤中HHT变换的具体处理过程。本发明的质量检测单元7的具体处理过程见超声波钢管混凝土柱质量检测方法 的步骤400 质量检测，即，根据HHT变换后的超声波信号进行钢管混凝土柱质量检测。本发明的优选实施方式是本发明还包括对所述超声波接收器接收的超声波信号 进行预处理的预处理单元4，所述预处理单元4在滤波前先对超声波信号末尾进行补零延 拓，在滤波完成之后再在超声波信号前端截掉相应的数据长度。本发明的具体实施例中，补 零长度设置为相位延迟时间的长度。本发明一种超声波钢管混凝土柱质量检测系统，通过对所述超声波接收器接收所 述超声波发生器发射的经钢管混凝土柱的超声波进行滤波处理后再进行HHT变换后，根据 HHT变换后的超声波信号进行钢管混凝土柱质量检测。本方法可靠性强，检测准确，方便实 用。以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定 本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在 不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的 保护范围。
权利要求
一种超声波钢管混凝土柱质量检测方法，其特征在于，包括检测钢管混凝土柱质量的超声波检测仪，将超声波检测仪的超声波发生器和超声波接收器分别设置在所述预检测钢管混凝土柱的两侧，超声波钢管混凝土柱质量检测方法包括如下步骤发射和接收超声波所述超声波发生器发射超声波，所述超声波接收器接收所述超声波发生器发射的经钢管混凝土柱的超声波；滤波处理对接收的超声波信号进行滤波处理；HHT变换对经过滤波处理的超声波信号进行HHT变换，即(经验模态分解与Hilbert谱分析)；质量检测根据HHT变换后的超声波信号进行钢管混凝土柱质量检测。
2.根据权利要求1所述的一种超声波钢管混凝土柱质量检测方法，其特征在于，在滤 波处理步骤中，包括在滤波前先对超声波信号末尾进行补零延拓。
3.根据权利要求2所述的一种超声波钢管混凝土柱质量检测方法，其特征在于，在滤 波处理步骤中，补零长度设置为相位延迟时间的长度。
4.根据权利要求3所述的一种超声波钢管混凝土柱质量检测方法，其特征在于，在滤 波处理步骤中，还包括在滤波完成之后再在超声波信号前端截掉相应的数据长度。
5.根据权利要求1所述的一种超声波钢管混凝土柱质量检测方法，其特征在于，在质 量检测步骤中，包括钢管混凝土脱层位置的确定。
6.根据权利要求1所述的一种超声波钢管混凝土柱质量检测方法，其特征在于，在质 量检测步骤中，包括钢管混凝土内含杂物位置的确定。
7.根据权利要求1所述的一种超声波钢管混凝土柱质量检测方法，其特征在于，在质 量检测步骤中，包括钢管混凝土中心孔洞位置的确定。
8.一种应用超声波钢管混凝土柱质量检测方法的超声波钢管混凝土柱质量检测系统， 包括检测钢管混凝土柱质量的超声波检测仪，将超声波检测仪的超声波发生器和超声波接 收器分别设置在所述预检测钢管混凝土柱的两侧，其特征在于，还包括将所述超声波接收 器接收的超声波进行滤波处理的滤波单元，对所述滤波单元滤波后的超声波信号进行HHT 变换的HHT变换单元，根据HHT变换后的超声波信号检测钢管混凝土柱质量的质量检测单元。
9.根据权利要求8所述的超声波钢管混凝土柱质量检测系统，其特征在于，还包括对 所述超声波接收器接收的超声波信号进行预处理的预处理单元，所述预处理单元在滤波前 先对超声波信号末尾进行补零延拓，在滤波完成之后再在超声波信号前端截掉相应的数据 长度。
10.根据权利要求9所述的超声波钢管混凝土柱质量检测系统，其特征在于，补零长度 设置为相位延迟时间的长度。
全文摘要
本发明涉及一种超声波钢管混凝土柱质量检测方法，包括检测钢管混凝土柱质量的超声波检测仪，将超声波检测仪的超声波发生器和超声波接收器分别设置在所述预检测钢管混凝土柱的两侧，超声波钢管混凝土柱质量检测方法包括如下步骤发射和接收超声波所述超声波发生器发射超声波，所述超声波接收器接收所述超声波发生器发射的经钢管混凝土柱的超声波；滤波处理对接收的超声波信号进行滤波处理；HHT变换对经过滤波处理的超声波信号进行HHT变换，即(经验模态分解与Hilbert谱分析)；质量检测根据HHT变换后的超声波信号进行钢管混凝土柱质量检测。
文档编号G01N29/04GK101907606SQ20101022562
公开日2010年12月8日 申请日期2010年7月12日 优先权日2010年7月12日
发明者叶福相, 查晓雄 申请人:哈尔滨工业大学深圳研究生院