基于希尔伯特变换和数据筛选的核电站松动件定位方法技术

基于希尔伯特变换和数据筛选的核电站松动件定位方法，在被撞击物上划分网格；安装加速度传感器；计算每个网格中心点到两两加速度传感器的距离差；采集各个加速度传感器接收到的撞击信号；采用移动平均法消除信号的低频噪声，采用Butterworth滤波器将环境背景噪声滤除；对消噪后的冲击信号进行希尔伯特变换，得到冲击信号的起振点时刻；根据网格中心点到各传感器的传播距离差和各通道之间的到达时间差计算每个网格中心点的名义传播速度；对所有网格中心点的名义传播速度进行数据筛选，计算其方差；找所有方差中最小值，并记录其中心点的坐标，将定位结果显示出来。本发明专利技术具有定位精度高，搜寻速度快的优点。

【技术实现步骤摘要】

本专利技术涉及核电站
，具体涉及基于希尔伯特变换和数据筛选的核电站松动件定位方法。技术背景核电站中存在大量的螺钉、螺帽等连接件，由于高速水流的不断冲击，出现腐蚀、折旧而松动甚至跌落，还有在系统测试、补给燃料、大修阶段从外界进入系统的金属碎片，这都会使得系统运行的稳定性和可靠性降低，甚至影响到整个核电站的安全。松动件定位作为松动件检测系统的重要组成部分，准确地定位松动件有利于在停堆检修时快速找到松动件，并进行相应的修补，尽量减少维修人员暴露在核辐射下的时间，保证维修人员的安全，对于核电站的稳定性和安全性有很大的帮助。现有的松动件定位相关文献有：G..Por,J.Kiss,I.Sorosanszky,G..Szappanos，Development of a false alarm free advanced loose parts monitoring system(ALPS)[J]，Progress in Nuclear Energy，2003,43(1-4):243-251.一种基于信号SPRT(Sequence Probability Ratio Test，顺序概率比测试)的时差估计方法，先用噪声的AR模型对信号进行白化，然后通过计算白化信号的SPRT来估算信号到达时差。Yong Beum Kim,Seon Jae Kim,Hae Dong Chung,Yong Won Park,<br>Jin Ho Park,A Study on Technique to Estimate Impact Location of Loose Part Using Wigner-Ville Distribution[J],Progress in Nuclear Energy,2003,43(1-4):261-266.一种基于魏格纳-威利分布的松动件冲击位置的估计方法研究，提出了对信号进行魏格纳-威利变换，获取时频图，进而获取松动件跌落位置。该方法准确性高，但计算量很大。每个信号的时频域线都是不一样的，而且需要手工来绘制，不利于实现自动化。S.Figedy,G..Oksa,Modern methods of signal processing in the loose part monitoring system[J],Progress in Nuclear Energy,2005,46(3-4):253-267.一种基于小波去噪的时差估计方法，通过小波去噪，即去除噪声的影响，然后再估计时差。该方法由于去除了噪声的影响，所以在低信噪比的情况下仍然有较好的估算效果，但是由于该方法还是以信号的过零点作为信号到达时间，没有考虑弯曲波的复杂传播模式，所以在实际定位时当传播距离较远时仍然存在较大误差。Seong-Nam Jeong,Kyoung-Hang Woo,Eoun-Taeg Hwang,Won-HoChoi,A Study on the Estimation Method of Impact Position Using the Frequency Analysis[J],Strategic Technology.The 1st International Forum on,2006:392-395.一种基于频率分析的冲击位置估计方法研究，提出了一种基于频散的松动件定位方法。该方法准确性高，稳定性也比较高，但前期需要建立起比较完善的数据库，工作量大。
技术实现思路
为了达到上述目的，本专利技术采取的技术方案为：基于希尔伯特变换和数据筛选的核电站松动件定位方法，包括以下步骤：1)根据定位精度要求，在被撞击物上划分网格，然后对划分好的网格编号1～N；2)在被撞击物上按等边三角形布置安装3个加速度传感器，以获取松动件跌落时产生的撞击信号f(t)；3)根据被撞击物的几何形状计算每个网格中心点到各两两加速度传感器的距离差di,j；4)通过数据采集卡同步采集各个加速度传感器接收到的撞击信号f(t)，撞击信号f(t)包括松动件跌落时的冲击信号s(t)和环境背景噪声n(t)；5)采用移动平均法消除撞击信号f(t)的低频噪声，再采用8阶Butterworth滤波器将环境背景噪声n(t)滤除，得到消除了噪声干扰的冲击信号s(t)；6)对步骤5)中处理得到的冲击信号s(t)用希尔伯特变换进行分析，得到冲击信号s(t)的起振点时刻；由于无法获知冲击信号s(t)跌落的时刻，只能从3个加速度传感器的冲击信号s(t)中获得波传播的时间差，每两个通道之间的到达时间差为ti,j；7)根据网格中心点到各加速度传感器的传播距离差di,j和各通道之间的到达时间差ti,j计算每个网格中心点的名义传播速度vi,j；8)计算弯曲波在结构中传播的速度上下限vmax和vmin；9)对所有网格中心点的名义传播速度vi,j进行数据筛选：判断名义传播速度vi,j是否在弯曲波传播的速度区间[vmax,vmin]内，若网格中心点的名义传播速度vi,j均在速度区间[vmax,vmin]内，则计算该网格中心点的速度的方差D(v)，否则给方差D(v)赋以一个大值；10)搜索所有网格中心点的方差D(v)中最小值，记录其中心点的坐标；11)定位结果显示。所述的步骤3)中，对于平面，根据公式(1)计算网格中心到各个加速度传感器的距离差：di,j=(x-xi)2+(y-yi)2-(x-xj)2+(y-yi)2---(1)]]>其中(x,y)为网格中心点的坐标，(xi,yi)为加速度传感器i的坐标，(xj,yj)为加速度传感器j的坐标，i＝1,2,3；j＝1,2,3；i≠j；对于半球面，根据公式(2)计算网格中心到各个加速度传感器的距离差：其中r为球的半径，为网格中心点的球坐标，为加速度传感器i的球坐标，为加速度传感器j的球坐标，i＝1,2,3；j＝1,2,3；i≠j。所述的步骤5)中，根据公式(3)消除撞击信号f(t)低频噪声：f′(t)=(f(t)+f(t-1)+...f(t-n+1))n---(3)]]>其中f'(t)为移动平均后的撞击信号，n为移动平均的时期个数，n设置为7。所述的步骤6)中，根据包络线的峰值来确定冲击信号s(t)起振点时刻，具体包括以下步骤：6.1)求冲击信号s(t)的希尔伯特变换根据定义：s^(t)=∫-∞∞s(t)1&pi本文档来自技高网...

【技术保护点】
基于希尔伯特变换和数据筛选的核电站松动件定位方法，其特征在于，包括以下步骤：1)根据定位精度要求，在被撞击物上划分网格，然后对划分好的网格编号1～N；2)在被撞击物上按等边三角形布置安装3个加速度传感器，以获取松动件跌落时产生的撞击信号f(t)；3)根据被撞击物的几何形状计算每个网格中心点到各两两加速度传感器的距离差di,j；4)通过数据采集卡同步采集各个加速度传感器接收到的撞击信号f(t)，撞击信号f(t)包括松动件跌落时的冲击信号s(t)和环境背景噪声n(t)；5)采用移动平均法消除撞击信号f(t)的低频噪声，再采用8阶Butterworth滤波器将环境背景噪声n(t)滤除，得到消除了噪声干扰的冲击信号s(t)；6)对步骤5)中处理得到的冲击信号s(t)用希尔伯特变换进行分析，得到冲击信号s(t)的起振点时刻；由于无法获知冲击信号s(t)跌落的时刻，只能从3个加速度传感器的冲击信号s(t)中获得波传播的时间差，每两个通道之间的到达时间差为ti,j；7)根据网格中心点到各加速度传感器的传播距离差di,j和各通道之间的到达时间差ti,j计算每个网格中心点的名义传播速度vi,j；8)计算弯曲波在结构中传播的速度上下限vmax和vmin；9)对所有网格中心点的名义传播速度vi,j进行数据筛选：判断名义传播速度vi,j是否在弯曲波传播的速度区间[vmax,vmin]内，若网格中心点的名义传播速度vi,j均在速度区间[vmax,vmin]内，则计算该网格中心点的速度的方差D(v)，否则给方差D(v)赋以一个大值；10)搜索所有网格中心点的方差D(v)中最小值，记录其中心点的坐标；11)定位结果显示。...

【技术特征摘要】
1.基于希尔伯特变换和数据筛选的核电站松动件定位方法，其
特征在于，包括以下步骤：
1)根据定位精度要求，在被撞击物上划分网格，然后对划分好
的网格编号1～N；
2)在被撞击物上按等边三角形布置安装3个加速度传感器，以
获取松动件跌落时产生的撞击信号f(t)；
3)根据被撞击物的几何形状计算每个网格中心点到各两两加速
度传感器的距离差di,j；
4)通过数据采集卡同步采集各个加速度传感器接收到的撞击信
号f(t)，撞击信号f(t)包括松动件跌落时的冲击信号s(t)和环境背景噪
声n(t)；
5)采用移动平均法消除撞击信号f(t)的低频噪声，再采用8阶
Butterworth滤波器将环境背景噪声n(t)滤除，得到消除了噪声干扰的
冲击信号s(t)；
6)对步骤5)中处理得到的冲击信号s(t)用希尔伯特变换进行分
析，得到冲击信号s(t)的起振点时刻；由于无法获知冲击信号s(t)跌落
的时刻，只能从3个加速度传感器的冲击信号s(t)中获得波传播的时
间差，每两个通道之间的到达时间差为ti,j；
7)根据网格中心点到各加速度传感器的传播距离差di,j和各通道
之间的到达时间差ti,j计算每个网格中心点的名义传播速度vi,j；
8)计算弯曲波在结构中传播的速度上下限vmax和vmin；
9)对所有网格中心点的名义传播速度vi,j进行数据筛选：判断名
义传播速度vi,j是否在弯曲波传播的速度区间[vmax,vmin]内，若网格中
心点的名义传播速度vi,j均在速度区间[vmax,vmin]内，则计算该网格中
心点的速度的方差D(v)，否则给方差D(v)赋以一个大值；
10)搜索所有网格中心点的方差D(v)中最小值，记录其中心点的
坐标；
11)定位结果显示。
2.根据权利要求1所述的定位方法，其特征在于：所述的步骤3)
中，对于平面，根据公式(1)计算网格中心到各个加速度传感器的
距离差：
di,j=(x-xi)2+(y-yi)2-(x-xj)2+(y-yj)2---(1)]]>其中(x,y)为网格中心点的坐标，(xi,yi)为加速度传感器i的坐标，
(xj,yj)为加速度传感器j的坐标，i＝1,2,3；j＝1,2,3；i≠j；
对于半球面，根据公式(2)计算网格中心到各个加速度传感器
的距离差：
其中r为球的半径，为网格中心点的球坐标，为加速
度传感器i的球坐标，为加速度传感器j的球坐标，
i＝1,2,3；j＝1,2...

【专利技术属性】
技术研发人员：林京，罗志，曾亮，黄利平，高飞，华佳东，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：陕西;61