[0058] 下面结合附图和实施例对本发明作进一步的说明。
[0059] 如图7所示，基于小波变换和希尔伯特变换的核电站松动部件定位方法，包括以下步骤：
[0060] 步骤(1)网格划分
[0061] 根据定位精度要求，在被撞击物上划分网格，定位精度要求越高，网格划分越小;网格长宽为 ，然后对划分好的网格编号。
[0062] 步骤(2)标定信号采集
[0063] 通过数据采集卡、电荷放大器和加速度传感器来实现标定信号的采集，将三个加速度传感器安装在被撞击物上，其中加速度传感器采集力锤敲击信号后，经电荷放大器进入数据采集卡；每个通道的标定信号 包括力锤敲击信号和本底噪声；
[0064] 所述的力锤敲击信号为力锤敲击不同网格内不同点所得的信号；
[0065] 步骤(3)标定信号 的连续小波变换
[0066] 对获取的标定信号 判断幅值，以最大幅值前0.03s作为开始端，取连续0.3s长信号数据 ，然后对 进行连续小波变换,获得处理后信号 。连续小波变换的尺度分别取2、8、14、20。
[0067] 步骤(4)计算到达时间
[0068] 对步骤(3)中处理得到的信号 用希尔伯特变换法进行分析，确定信号的起振点时刻，第一通道、第二通道、第三通道的起振点时刻分别为 、 、 ;
[0069] 步骤(5)建立定位数据库
[0070] 根据步骤(4)中得到的不同通道的起振点时刻 、 、 ,计算第二通道、第三通道与第一通道的到达时间差 、 ；计算得出同一网格内不同敲击点的多次敲击结果的时间差平均值 、 为；然后以网格编号，时间差平均值 和 为查询表的属性，建立定位查询表grid_tbl。
[0071] 步骤(6)现场数据采集
[0072] 通过数据采集卡、电荷放大器和加速度传感器来实现现场数据信号采集。现场数据信号 包括松动件跌落时的冲击信号和环境背景噪声。
[0073] 步骤(7)现场数据信号 的傅里叶变换
[0074] 对获取的现场数据信号 判断幅值，以最大幅值前0.03s作为开始端，取0.3s长信号数据 ，并对 进行傅里叶变换；根据信号频谱图确定 的主频率 ，然后根据公式： ，求得小波尺度 ,式中 为采样周期； 为给定小波的中心频率，该中心频率为归一化后的数值，例如db2小波，其中心频率为0.6667。
[0075] 步骤(8)对 进行连续小波变换
[0076] 将步骤(3)中连续小波变换尺度2、8、14、20与步骤(7)中求得的小波尺度 作比较，取最接近 的作为新连续小波变换尺度M，然后以尺度M对步骤(7)中得到的信号进行连续小波变换处理得到信号 。
[0077] 步骤(9)希尔伯特变换
[0078] 对步骤(8)中处理得到的信号 用希尔伯特变换法进行分析，确定该信号的起振点时刻，第一通道、第二通道、第三通道的起振点时刻分别为 、 、 。
[0079] 步骤(10)查表定位
[0080] 根据步骤(9)中确定的各通道的起振时刻，求得第二通道、第三通道与第一通道的信号到达时间差 、 ；然后根据步骤(8)中确定的小波尺度M，找到M尺度下的定位查询表grid_tbl；用到达时间差 、 与定位查询表grid_tbl中平均值的 、 分别作比较 ，求其时间距离 ，时间距离 中的最小值所对应的网格编号 即为定位网格N。
[0081] 所述的步骤(3)和步骤(8)中的连续小波变换的公式为
[0082]
[0083] 其中， 为步骤(3)中的 ，为步骤(8)中的 ； 为小波函数，此处取小波函数为db2； 为 的复共轭；为尺度因子；为时间平移因子。
即为经小波尺度 连续变换后的信号。
[0084] 所述的步骤(4)和步骤(9)中采用的希尔伯特变换到达时差定位法，是在时域通过求得信号的希尔伯特包络线，根据包络线的峰值来确定信号起振点时刻；包括以下步骤：
[0085] (1)求输入信号 的希尔伯特变换 ，根据定义：
[0086] ，其中， 即为 经希尔伯特变换后的信号。
[0087] (2)以信号 为实部，其希尔伯特变换 为虚部，构成一个新的函数如式：
[0088] 其中， 是新函数的幅值函数； 为相位函数。则 即为信号 的包络线函数。
[0089] (3)对获取的包络线函数 判断幅值，取信号 幅值突然变大前的一段信号（长度为0.1s）作为背景噪声。取背景噪声的最大幅值作为阈值 。
[0090] (4)对 求峰值，比较点数一般设置为5，即取 中大于阈值 且大于与该点相邻的5个点的所有值，则此点所对应的值 即为 的一个峰值。取求得的第一个峰值所对应的时间点t即为信号的起振时刻。
[0091] 所述的步骤(10)中时间距离的计算公式为：
[0092] ，
[0093] 其中，为网格编号， 为网格划分总数； 的最小值所对应的 即为定位网格N。
[0094] 希尔伯特定位法原理图如附图1所示，当松动部件和反应堆的压力容器发生碰撞后，碰撞能量主要以弯曲波的形式进行传播，其传播速度取决于弯曲波的群速度。假设一弯曲波由n个频率很接近的简谐波组成，且其角频率分别为 ， ，…， ，主频率为 (即各频率分量都接近于 )。合成波 可以表示成： ，式中为碰撞点与传感器之间的距离， 为波数， 为振幅。根据合成波构建新函数：
[0095]
[0096] 则信号的包络线函数可表示为：
[0097]
[0098]
[0099]
[0100] 由上式知当 取得最大值时，即时， 取最大值。根据波传播理论，弯曲波的群速度为波的角频率对波数的导数，即： ，由于各个角频率很接近，所以可得： ，因此包络线上取得最大值时
的时间点t可以表示为： ，所以冲击信号包络线的极值点所对应的时间点即为碰撞产生的弯曲波以速度 传播距离 的时刻。如附图1所示，因此信号包络线的第一峰值即为信号到达时间。
实施例
[0101] 下面结合试验，说明本发明，本发明的效果可以通过平板试验的实验结果分析说明：
[0102] 1、试验条件
[0103] 如附图2所示，本试验的试验平台由测试对象、加速度传感器、电荷放大器、数据采集卡和计算机构成。测试对象包括：金属球、卧式锅炉和钢板及其支撑。钢板尺寸为200cm*150cm*2cm。为了尽量减少环境噪声的影响，在钢板的四个边沿下均加了缓冲隔离。
缓冲隔离由钢板和橡胶板构成，缓冲隔离中钢板尺寸为20cm*20cm*1.2cm,橡胶板尺寸为
20cm*20cm*2cm。每个缓冲隔离由3块钢板和3块橡胶板组成，由底层开始分别为钢板、橡胶板、钢板、橡胶板、钢板、橡胶板，总厚度约10cm。金属球重量分别为20g、44g、100g、175g、
225g、360g、510g、640g、880g、1.12kg、1.4kg、1.8kg、2.52kg、4.1kg、9kg、11.2kg。
[0104] 2、数据采集
[0105] 数据采集包括两部分：力锤敲击信号采集和背景噪声采集。
[0106] (1)冲击信号采集：采用三个加速度传感器，加速度传感器在钢板上呈三角形排列，可参看如附图3。冲击物为力锤和不同质量金属球。标定信号为力锤敲击每个网格中心点多次冲击信号。力锤敲击信号为不同质量金属球，分别在高度10cm、15cm、20cm处敲击钢板不同点的冲击信号。
[0107] (2)背景噪声采集：采用三个加速度传感器，其中一个加速度传感器位于锅炉顶部进水口处，一个位于锅炉出水口，一个位于锅炉顶部左上方。背景噪声为锅炉从冷态到热态过程中的背景噪声，根据锅炉运行状态（锅炉未运行前、小火、大火、水泵运转、放气等）来采集数据，锅炉运行状态每改变一次采集一组数据。
[0108] 3、冲击试验结果及分析
[0109] 力锤敲击信号与背景噪声分别按信噪比为-5dB、0dB、5bB、10dB进行叠加，用于试验结分析。信噪比定义为：
[0110]
[0111] 其中 表示噪声幅值取绝对值后求平均， 表示信号幅值取绝对值后求最大值。
[0112] 实验结果对本发明提出的定位方法进行了验证。实验结果如附图4和附图5所示，从附图4中可看出金属球的冲击高度对冲击信号的频率成分影响并不大，这正好与赫兹碰撞理论相符。从附图5中可以看出信号经连续小波变换后，在保持信号能量不变的前提下，其信号能量主频率部分信号被放大，信号频率成分变得集中。从而减弱了低频噪声及其他频率成分波的干扰。
[0113] 本发明用如下方法来计算矢量误差和相对误差。
[0114] 矢量误差： ，其中 , 分别为松动件跌落位置所在网格与定位位置所在网格之间的水平间隔网格数和垂直网格间隔数。
[0115] 相对误差： ，其中 为矢量误差的模， 为网格的面积，为以三个定位传感器为顶点构成的三角型的面积。
[0116] 试验结果对本发明提出的松动件定位方法进行了验证，并且还将本发明提出的方法与其他方法进行了比较分析。
[0117] 通过试验数据分析，试验定位分析结果如表1~4所示。
[0118] 表1 小波尺度2下分析定位结果
[0119]
[0120] 表2 小波尺度8下分析定位结果(1)
[0121]
[0122] 表3 小波尺度8下分析定位结果(2)
[0123]
[0124] 表4 小波尺度20下分析定位结果
[0125]
[0126] 由表1~4可以看出,所有钢球跌落位置的定位平均相对误差最小为2.60%，最大为13.79%。根据矢量误差可以看出，定位网格与松动件实际跌落网格的间隔在1个网格左右。说明该方法对不同质量松动件都有较好的定位效果，大质量的松动件定位效果最好。
[0127] 以175g不含噪声和880g加噪声（信噪比为5dB）钢球为例，分析不同定位方法的估计效果。如表5和表6所示。信号经相同尺度 连续分解后，再经不同时间差估算方法确定时间差，最后都由查表法定位。
[0128] 表5 175g-尺度8
[0129]
[0130] 表6 880g-尺度8(5dB)
[0131]
[0132] 由表5和表6可以看出，本发明提出的方法与其他两种方法相比定位精度要高，且抗干扰能力较强。
[0133] 本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围不应当被视为仅限于实施例所所陈述的具体形式，本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。
