CN103137225A - 基于小波变换和希尔伯特变换的核电站松动部件定位方法 - Google Patents

基于小波变换和希尔伯特变换的核电站松动部件定位方法 Download PDF

Info

Publication number
CN103137225A
CN103137225A CN2013100314165A CN201310031416A CN103137225A CN 103137225 A CN103137225 A CN 103137225A CN 2013100314165 A CN2013100314165 A CN 2013100314165A CN 201310031416 A CN201310031416 A CN 201310031416A CN 103137225 A CN103137225 A CN 103137225A
Authority
CN
China
Prior art keywords
signal
grid
transform
wavelet
channel
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
Application number
CN2013100314165A
Other languages
English (en)
Other versions
CN103137225B (zh
Inventor
方力先
季田田
曾复
张卫
万鹏飞
谢永诚
王赤虎
张海丰
邓晶晶
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Dianzi University
Shanghai Nuclear Engineering Research and Design Institute Co Ltd
Original Assignee
Hangzhou Dianzi University
Shanghai Nuclear Engineering Research and Design Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hangzhou Dianzi University, Shanghai Nuclear Engineering Research and Design Institute Co Ltd filed Critical Hangzhou Dianzi University
Priority to CN201310031416.5A priority Critical patent/CN103137225B/zh
Publication of CN103137225A publication Critical patent/CN103137225A/zh
Application granted granted Critical
Publication of CN103137225B publication Critical patent/CN103137225B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种基于小波变换和希尔伯特变换的核电站松动部件定位方法。本发明包括如下步骤：步骤(1)网格划分，对划分好的网格编号；步骤(2)采集标定信号采集；步骤(3)标定信号
Figure 2013100314165100004DEST_PATH_IMAGE002
的连续小波变换；步骤(4)计算到达时间；步骤(5)建立定位数据库；步骤(6)采集现场数据；步骤(7)现场数据信号
Figure 2013100314165100004DEST_PATH_IMAGE004
的傅里叶变换；步骤(8)对进行连续小波变换；步骤(9)希尔伯特变换；步骤(10)查表定位。本发明具有实时性好、定位精度高、抗干扰能力强及计算快速方便的优点。

Description

基于小波变换和希尔伯特变换的核电站松动部件定位方法
技术领域
本发明属于核工程技术领域，尤其涉及一种基于小波变换和希尔伯特变换的核电站松动部件定位方法，用于估计核电站松动部件跌落位置。
技术背景
核电站中由于组件老化、构件疲劳而造成金属部件的松动或脱落。以及在主系统施工、加料、维护期间因疏忽无意中遗留下的金属部件,都会加大系统中其他部分的组件损伤和材料磨损，从而削弱反应堆系统运行的安全性能和稳定性能。此外,松动部件会增加控制杆的潜在干扰并增加主系统中放射性残渣的积累量，对核电厂工作人员的安全也造成威胁。松动部件定位作为核电站松动部件监测系统的主要组成部分，准确的定位松动部件有利于在停堆检修时快速找到松动部件，并进行相应的处理，尽量减少维修人员暴露在核辐射下的时间，保证维修人员的安全，是评价一回路系统内部件机械安全性和提高核反应堆安全性的可靠手段。
现有的松动件定位方法相关文献有：
[1] J. F. Doyle, An experimental method for determining the location and time of initiation of an unknown dispersing pulse[J]，Experimental Mechanics, 1987, 27: 229-233.
    一种基于高阶统计量的时差估计方法，利用高阶统计量对高斯噪声的不敏感性这一优点，以二阶或更高阶的相关函数来估计被噪声污染的信号间的时间延迟。与传统相关函数相比，则能够去除高斯噪声的影响，但当噪声为有色时误差还是较大。
[2] 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来估算信号到达时差。
[3] 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.
一种基于小波去噪的时差估计方法，通过小波去噪，即去除噪声的影响，然后再估计时差。该方法由于去除了噪声的影响，所以在低信噪比的情况下仍然有较好的估算效果，但是由于该方法还是以信号的过零点作为信号到达时间，没有考虑弯曲波的复杂传播模式，所以在实际定位时当传播距离较远时仍然存在较大误差。
[4] Gee Yong Park . Se Woo Cheon ,Cheol Kwon Lee, et al. An Estimation Method for Impact Location of Loose Parts [J] .Progress in Nuclear Energy, 2005, 48(2006): 360-370.
一种松动件冲击位置的估计方法，提出了一种基于信号包络线的确定冲击信号到达传感器的起始时间的方法，但是现场实验数据的分析结果很不理想。
发明内容
本发明的目的是为了克服现有定位技术在实际应用中定位精度差、抗干扰能力弱以及计算速度慢等缺陷，提出了一种基于小波变换和希尔伯特变换的核电站松动部件定位方法，该方法抗干扰能力强，定位精度高、可通过计算机实现自动快速定位。
本发明解决其技术问题所采用的技术方案包括以下步骤：
步骤(1)网格划分
    根据定位精度要求，在被撞击物上划分网格，定位精度要求越高，网格划分越小;网格长宽为 
Figure 2013100314165100002DEST_PATH_IMAGE002
，然后对划分好的网格编号。
步骤(2)标定信号采集
通过数据采集卡、电荷放大器和加速度传感器来实现标定信号的采集，将三个加速度传感器安装在被撞击物上，其中加速度传感器采集力锤敲击信号后，经电荷放大器进入数据采集卡；每个通道的标定信号包括力锤敲击信号和本底噪声；
所述的力锤敲击信号为力锤敲击不同网格内不同点所得的信号；
步骤(3)标定信号的连续小波变换
对获取的标定信号
Figure 490730DEST_PATH_IMAGE004
判断幅值，以最大幅值前0.03s作为开始端，取连续0.3s长信号数据
Figure 2013100314165100002DEST_PATH_IMAGE006
，然后对
Figure 598363DEST_PATH_IMAGE006
进行连续小波变换,获得处理后信号。连续小波变换的尺度分别取2、8、14、20。
步骤(4)计算到达时间
对步骤(3)中处理得到的信号用希尔伯特变换法进行分析，确定信号
Figure 106147DEST_PATH_IMAGE008
的起振点时刻，第一通道、第二通道、第三通道的起振点时刻分别为
Figure 2013100314165100002DEST_PATH_IMAGE010
、
Figure 2013100314165100002DEST_PATH_IMAGE012
、
Figure 2013100314165100002DEST_PATH_IMAGE014
;
步骤(5)建立定位数据库
根据步骤(4)中得到的不同通道的起振点时刻
Figure 115560DEST_PATH_IMAGE010
、
Figure 382594DEST_PATH_IMAGE012
、
Figure 407925DEST_PATH_IMAGE014
,计算第二通道、第三通道与第一通道的到达时间差
Figure 2013100314165100002DEST_PATH_IMAGE016
、
Figure 2013100314165100002DEST_PATH_IMAGE018
；计算得出同一网格内不同敲击点的多次敲击结果的时间差平均值
Figure 2013100314165100002DEST_PATH_IMAGE020
、
Figure 2013100314165100002DEST_PATH_IMAGE022
为；然后以网格编号，时间差平均值
Figure 263755DEST_PATH_IMAGE020
和
Figure 555321DEST_PATH_IMAGE022
为查询表的属性，建立定位查询表grid_tbl。
步骤(6)现场数据采集
通过数据采集卡、电荷放大器和加速度传感器来实现现场数据信号采集。现场数据信号
Figure 2013100314165100002DEST_PATH_IMAGE024
包括松动件跌落时的冲击信号和环境背景噪声。
步骤(7)现场数据信号
Figure 371967DEST_PATH_IMAGE024
的傅里叶变换
对获取的现场数据信号
Figure 436875DEST_PATH_IMAGE024
判断幅值，以最大幅值前0.03s作为开始端，取0.3s长信号数据
Figure 2013100314165100002DEST_PATH_IMAGE026
，并对
Figure 756998DEST_PATH_IMAGE026
进行傅里叶变换；根据信号频谱图确定
Figure 186842DEST_PATH_IMAGE026
的主频率
Figure 2013100314165100002DEST_PATH_IMAGE028
，然后根据公式：
Figure 2013100314165100002DEST_PATH_IMAGE030
，求得小波尺度
Figure 2013100314165100002DEST_PATH_IMAGE032
,式中
Figure 2013100314165100002DEST_PATH_IMAGE034
为采样周期；
Figure 2013100314165100002DEST_PATH_IMAGE036
为给定小波的中心频率，该中心频率为归一化后的数值，例如db2小波，其中心频率为0.6667。
步骤(8)对进行连续小波变换
将步骤(3)中连续小波变换尺度2、8、14、20与步骤(7)中求得的小波尺度
Figure 32754DEST_PATH_IMAGE032
作比较，取最接近
Figure 738542DEST_PATH_IMAGE032
的作为新连续小波变换尺度M，然后以尺度M对步骤(7)中得到的信号进行连续小波变换处理得到信号。
步骤(9)希尔伯特变换
对步骤(8)中处理得到的信号
Figure 2013100314165100002DEST_PATH_IMAGE040
用希尔伯特变换法进行分析，确定该信号的起振点时刻，第一通道、第二通道、第三通道的起振点时刻分别为
Figure 2013100314165100002DEST_PATH_IMAGE042
、
Figure 2013100314165100002DEST_PATH_IMAGE044
、
Figure 2013100314165100002DEST_PATH_IMAGE046
。
步骤(10)查表定位
根据步骤(9)中确定的各通道的起振时刻，求得第二通道、第三通道与第一通道的信号到达时间差
Figure 2013100314165100002DEST_PATH_IMAGE048
、
Figure 2013100314165100002DEST_PATH_IMAGE050
；然后根据步骤(8)中确定的小波尺度M，找到M尺度下的定位查询表grid_tbl；用到达时间差、
Figure 940646DEST_PATH_IMAGE050
与定位查询表grid_tbl中平均值的
Figure 317662DEST_PATH_IMAGE020
、
Figure 815640DEST_PATH_IMAGE022
分别作比较 ，求其时间距离 ，时间距离
Figure 977500DEST_PATH_IMAGE052
中的最小值所对应的网格编号即为定位网格N。
    所述的步骤(3)和步骤(8)中的连续小波变换的公式为
Figure 2013100314165100002DEST_PATH_IMAGE056
其中，
Figure 2013100314165100002DEST_PATH_IMAGE058
为步骤(3)中的
Figure 557167DEST_PATH_IMAGE008
，为步骤(8)中的；
Figure 2013100314165100002DEST_PATH_IMAGE060
为小波函数，此处取小波函数为db2；为
Figure 651210DEST_PATH_IMAGE060
的复共轭；
Figure 2013100314165100002DEST_PATH_IMAGE064
为尺度因子；
Figure 2013100314165100002DEST_PATH_IMAGE066
为时间平移因子。
Figure 2013100314165100002DEST_PATH_IMAGE068
即为经小波尺度
Figure 610070DEST_PATH_IMAGE032
连续变换后的信号。
所述的步骤(4)和步骤(9)中采用的希尔伯特变换到达时差定位法，是在时域通过求得信号的希尔伯特包络线，根据包络线的峰值来确定信号起振点时刻；包括以下步骤：
(1)求输入信号的希尔伯特变换
Figure 2013100314165100002DEST_PATH_IMAGE072
，根据定义：
       
Figure 2013100314165100002DEST_PATH_IMAGE074
 ，其中，
Figure 733491DEST_PATH_IMAGE072
即为
Figure 154108DEST_PATH_IMAGE070
经希尔伯特变换后的信号。
(2)以信号
Figure 423415DEST_PATH_IMAGE070
为实部，其希尔伯特变换
Figure 864761DEST_PATH_IMAGE072
为虚部，构成一个新的函数如式： 
Figure 2013100314165100002DEST_PATH_IMAGE076
 
其中， 
Figure 2013100314165100002DEST_PATH_IMAGE078
是新函数的幅值函数；
Figure 2013100314165100002DEST_PATH_IMAGE080
为相位函数。则即为信号
Figure 150511DEST_PATH_IMAGE070
的包络线函数。
(3)对获取的包络线函数
Figure 702715DEST_PATH_IMAGE082
判断幅值，取信号
Figure 560950DEST_PATH_IMAGE082
幅值突然变大前的一段信号（长度为0.1s）作为背景噪声。取背景噪声的最大幅值作为阈值
Figure 2013100314165100002DEST_PATH_IMAGE084
。
(4)对
Figure 665872DEST_PATH_IMAGE082
求峰值，比较点数一般设置为5，即取中大于阈值且大于与该点相邻的5个点的所有值，则此点所对应的值
Figure 2013100314165100002DEST_PATH_IMAGE086
即为
Figure 147035DEST_PATH_IMAGE082
的一个峰值。取求得的第一个峰值所对应的时间点t即为信号的起振时刻。
所述的步骤(10)中时间距离的计算公式为：
Figure 2013100314165100002DEST_PATH_IMAGE088
，
Figure 2013100314165100002DEST_PATH_IMAGE090
其中，
Figure 962807DEST_PATH_IMAGE054
为网格编号，
Figure 2013100314165100002DEST_PATH_IMAGE092
为网格划分总数；
Figure 2013100314165100002DEST_PATH_IMAGE093
的最小值所对应的
Figure 252842DEST_PATH_IMAGE054
即为定位网格N。
本发明的有益效果如下：
先根据采集信号的频谱结构确定信号主频率
Figure 818953DEST_PATH_IMAGE028
，根据频率与小波尺度
Figure 917359DEST_PATH_IMAGE032
的关系公式计算得到小波尺度
Figure 104365DEST_PATH_IMAGE032
。然后用尺度为
Figure 491484DEST_PATH_IMAGE032
的小波函数对信号进行连续小波变换。然后用希尔伯特变换法对连续小波变换后的信号进行分析计算，求得信号起振时刻，最后用查表法实现定位。小波变化具有时-频局部化特性。经连续小波变换后的信号可以实现特定时-频成分信号的突显，使信号所包含的能量主频率较集中，且小波变换是等能量变换。这就减小了波的频散效应对波速的影响。因为基于希尔伯特变换的时差估计法采用的是能量最大的弯曲波成分的速度，而经小波连续变换后的信号正好满足该条件，所以该方法对冲击信号起振点的确定较准确。又因为基于希尔伯特变换的时差估计法中认为信号希尔伯特包络线的第一个峰值为信号到达时刻，而不是信号过零点，所以该方法抗干扰能力强。最后通过查表法实现定位，而不像三角形定位法需要求解双曲线方程，从而不会受到多解和无解问题的困扰，而且计算速度快，计算方便，虽然在系统运行前需要做大量的前期工作，但定位精度能够得到保障。
     本发明具有实时性好、定位精度高、抗干扰能力强及计算快速方便的优点。
附图说明
图1为希尔伯特定位法原理图；
图2为平板试验示意图；
图3为敲击点位置及传感器安装位置图；
图4 为510g钢球频谱图；
图5为1.12kg钢球小波变换前后频谱比较图；
图6 880g无噪信号与880g加噪信号（信噪比为5dB）；
图7为本发明流程图。
具体实施方式
    下面结合附图和实施例对本发明作进一步的说明。
如图7所示，基于小波变换和希尔伯特变换的核电站松动部件定位方法，包括以下步骤：
步骤(1)网格划分
    根据定位精度要求，在被撞击物上划分网格，定位精度要求越高，网格划分越小;网格长宽为
Figure DEST_PATH_IMAGE094
，然后对划分好的网格编号。
步骤(2)标定信号采集
通过数据采集卡、电荷放大器和加速度传感器来实现标定信号的采集，将三个加速度传感器安装在被撞击物上，其中加速度传感器采集力锤敲击信号后，经电荷放大器进入数据采集卡；每个通道的标定信号
Figure DEST_PATH_IMAGE095
包括力锤敲击信号和本底噪声；
所述的力锤敲击信号为力锤敲击不同网格内不同点所得的信号；
步骤(3)标定信号
Figure 720340DEST_PATH_IMAGE095
的连续小波变换
对获取的标定信号
Figure 610935DEST_PATH_IMAGE095
判断幅值，以最大幅值前0.03s作为开始端，取连续0.3s长信号数据
Figure DEST_PATH_IMAGE096
，然后对进行连续小波变换,获得处理后信号
Figure DEST_PATH_IMAGE097
。连续小波变换的尺度分别取2、8、14、20。
步骤(4)计算到达时间
对步骤(3)中处理得到的信号
Figure 705241DEST_PATH_IMAGE097
用希尔伯特变换法进行分析，确定信号
Figure 206630DEST_PATH_IMAGE097
的起振点时刻，第一通道、第二通道、第三通道的起振点时刻分别为
Figure 217311DEST_PATH_IMAGE010
、
Figure 982005DEST_PATH_IMAGE012
、
Figure 78137DEST_PATH_IMAGE014
;
步骤(5)建立定位数据库
根据步骤(4)中得到的不同通道的起振点时刻
Figure 622031DEST_PATH_IMAGE010
、
Figure 549536DEST_PATH_IMAGE012
、
Figure 688393DEST_PATH_IMAGE014
,计算第二通道、第三通道与第一通道的到达时间差
Figure 68559DEST_PATH_IMAGE016
、
Figure 646171DEST_PATH_IMAGE018
；计算得出同一网格内不同敲击点的多次敲击结果的时间差平均值
Figure 631445DEST_PATH_IMAGE020
、 为；然后以网格编号，时间差平均值
Figure 310130DEST_PATH_IMAGE020
和
Figure 425853DEST_PATH_IMAGE022
为查询表的属性，建立定位查询表grid_tbl。
步骤(6)现场数据采集
通过数据采集卡、电荷放大器和加速度传感器来实现现场数据信号采集。现场数据信号
Figure DEST_PATH_IMAGE098
包括松动件跌落时的冲击信号和环境背景噪声。
步骤(7)现场数据信号
Figure 593530DEST_PATH_IMAGE098
的傅里叶变换
对获取的现场数据信号判断幅值，以最大幅值前0.03s作为开始端，取0.3s长信号数据
Figure DEST_PATH_IMAGE099
，并对进行傅里叶变换；根据信号频谱图确定
Figure 846897DEST_PATH_IMAGE099
的主频率
Figure 541184DEST_PATH_IMAGE028
，然后根据公式：
Figure DEST_PATH_IMAGE100
，求得小波尺度,式中为采样周期；
Figure 492325DEST_PATH_IMAGE036
为给定小波的中心频率，该中心频率为归一化后的数值，例如db2小波，其中心频率为0.6667。
步骤(8)对
Figure 369014DEST_PATH_IMAGE099
进行连续小波变换
将步骤(3)中连续小波变换尺度2、8、14、20与步骤(7)中求得的小波尺度
Figure 457056DEST_PATH_IMAGE032
作比较，取最接近
Figure 22292DEST_PATH_IMAGE032
的作为新连续小波变换尺度M，然后以尺度M对步骤(7)中得到的信号
Figure DEST_PATH_IMAGE101
进行连续小波变换处理得到信号。
步骤(9)希尔伯特变换
对步骤(8)中处理得到的信号
Figure 283509DEST_PATH_IMAGE102
用希尔伯特变换法进行分析，确定该信号的起振点时刻，第一通道、第二通道、第三通道的起振点时刻分别为
Figure 14704DEST_PATH_IMAGE042
、
Figure 70385DEST_PATH_IMAGE044
、
Figure 824714DEST_PATH_IMAGE046
。
步骤(10)查表定位
根据步骤(9)中确定的各通道的起振时刻，求得第二通道、第三通道与第一通道的信号到达时间差、
Figure 929121DEST_PATH_IMAGE050
；然后根据步骤(8)中确定的小波尺度M，找到M尺度下的定位查询表grid_tbl；用到达时间差、
Figure 194066DEST_PATH_IMAGE050
与定位查询表grid_tbl中平均值的
Figure 734769DEST_PATH_IMAGE020
、
Figure 706136DEST_PATH_IMAGE022
分别作比较 ，求其时间距离
Figure 339505DEST_PATH_IMAGE093
，时间距离中的最小值所对应的网格编号
Figure 209558DEST_PATH_IMAGE054
即为定位网格N。
    所述的步骤(3)和步骤(8)中的连续小波变换的公式为
Figure DEST_PATH_IMAGE103
其中，
Figure DEST_PATH_IMAGE104
为步骤(3)中的
Figure 35431DEST_PATH_IMAGE097
，为步骤(8)中的
Figure 338237DEST_PATH_IMAGE099
；
Figure DEST_PATH_IMAGE105
为小波函数，此处取小波函数为db2；
Figure DEST_PATH_IMAGE106
为的复共轭；
Figure DEST_PATH_IMAGE064A
为尺度因子；
Figure DEST_PATH_IMAGE066A
为时间平移因子。
Figure 748544DEST_PATH_IMAGE068
即为经小波尺度
Figure 101028DEST_PATH_IMAGE032
连续变换后的信号。
所述的步骤(4)和步骤(9)中采用的希尔伯特变换到达时差定位法，是在时域通过求得信号的希尔伯特包络线，根据包络线的峰值来确定信号起振点时刻；包括以下步骤：
(1)求输入信号
Figure 105893DEST_PATH_IMAGE070
的希尔伯特变换
Figure 871724DEST_PATH_IMAGE072
，根据定义：
        ，其中，
Figure 561648DEST_PATH_IMAGE072
即为
Figure 409519DEST_PATH_IMAGE070
经希尔伯特变换后的信号。
(2)以信号为实部，其希尔伯特变换为虚部，构成一个新的函数如式： 
Figure 780741DEST_PATH_IMAGE076
 
其中， 是新函数的幅值函数；
Figure 274356DEST_PATH_IMAGE080
为相位函数。则
Figure 364672DEST_PATH_IMAGE082
即为信号
Figure 844457DEST_PATH_IMAGE070
的包络线函数。
(3)对获取的包络线函数
Figure 299709DEST_PATH_IMAGE082
判断幅值，取信号
Figure 261849DEST_PATH_IMAGE082
幅值突然变大前的一段信号（长度为0.1s）作为背景噪声。取背景噪声的最大幅值作为阈值。
(4)对
Figure 191945DEST_PATH_IMAGE082
求峰值，比较点数一般设置为5，即取中大于阈值且大于与该点相邻的5个点的所有值，则此点所对应的值
Figure 106177DEST_PATH_IMAGE086
即为
Figure 88783DEST_PATH_IMAGE082
的一个峰值。取求得的第一个峰值所对应的时间点t即为信号的起振时刻。
所述的步骤(10)中时间距离的计算公式为：
Figure 885838DEST_PATH_IMAGE088
，
Figure 556991DEST_PATH_IMAGE090
其中，
Figure 730483DEST_PATH_IMAGE054
为网格编号，
Figure 537902DEST_PATH_IMAGE092
为网格划分总数；
Figure 505858DEST_PATH_IMAGE093
的最小值所对应的
Figure 929886DEST_PATH_IMAGE054
即为定位网格N。
希尔伯特定位法原理图如附图1所示，当松动部件和反应堆的压力容器发生碰撞后，碰撞能量主要以弯曲波的形式进行传播，其传播速度取决于弯曲波的群速度。假设一弯曲波由n个频率很接近的简谐波组成，且其角频率分别为，
Figure DEST_PATH_IMAGE110
，…，
Figure DEST_PATH_IMAGE112
，主频率为 (即各频率分量都接近于
Figure 736430DEST_PATH_IMAGE114
)。合成波
Figure DEST_PATH_IMAGE116
可以表示成： 
Figure DEST_PATH_IMAGE118
，式中
Figure DEST_PATH_IMAGE120
为碰撞点与传感器之间的距离，
Figure DEST_PATH_IMAGE122
为波数，为振幅。根据合成波构建新函数：  
Figure DEST_PATH_IMAGE126
 
Figure DEST_PATH_IMAGE128
                   则信号的包络线函数可表示为：
    
             
Figure DEST_PATH_IMAGE134
 
由上式知当
Figure DEST_PATH_IMAGE136
取得最大值时，即
Figure DEST_PATH_IMAGE138
时，
Figure DEST_PATH_IMAGE140
取最大值。根据波传播理论，弯曲波的群速度为波的角频率对波数的导数，即： 
Figure DEST_PATH_IMAGE142
，由于各个角频率很接近，所以可得： 
Figure DEST_PATH_IMAGE144
，因此包络线上取得最大值时的时间点t可以表示为： 
Figure DEST_PATH_IMAGE146
，所以冲击信号包络线的极值点所对应的时间点即为碰撞产生的弯曲波以速度
Figure DEST_PATH_IMAGE148
传播距离
Figure DEST_PATH_IMAGE150
的时刻。如附图1所示，因此信号包络线的第一峰值即为信号到达时间。
实施例
下面结合试验，说明本发明，本发明的效果可以通过平板试验的实验结果分析说明：
1、试验条件
如附图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。
2、数据采集
数据采集包括两部分：力锤敲击信号采集和背景噪声采集。
(1)冲击信号采集：采用三个加速度传感器，加速度传感器在钢板上呈三角形排列，可参看如附图3。冲击物为力锤和不同质量金属球。标定信号为力锤敲击每个网格中心点多次冲击信号。力锤敲击信号为不同质量金属球，分别在高度10cm、15cm、20cm处敲击钢板不同点的冲击信号。
(2)背景噪声采集：采用三个加速度传感器，其中一个加速度传感器位于锅炉顶部进水口处，一个位于锅炉出水口，一个位于锅炉顶部左上方。背景噪声为锅炉从冷态到热态过程中的背景噪声，根据锅炉运行状态（锅炉未运行前、小火、大火、水泵运转、放气等）来采集数据，锅炉运行状态每改变一次采集一组数据。
3、冲击试验结果及分析
力锤敲击信号与背景噪声分别按信噪比为-5dB、0dB、5bB、10dB进行叠加，用于试验结分析。信噪比定义为：
Figure DEST_PATH_IMAGE152
               
其中
Figure DEST_PATH_IMAGE154
表示噪声幅值取绝对值后求平均，
Figure DEST_PATH_IMAGE156
表示信号幅值取绝对值后求最大值。
实验结果对本发明提出的定位方法进行了验证。实验结果如附图4和附图5所示，从附图4中可看出金属球的冲击高度对冲击信号的频率成分影响并不大，这正好与赫兹碰撞理论相符。从附图5中可以看出信号经连续小波变换后，在保持信号能量不变的前提下，其信号能量主频率部分信号被放大，信号频率成分变得集中。从而减弱了低频噪声及其他频率成分波的干扰。
本发明用如下方法来计算矢量误差和相对误差。
矢量误差：
Figure DEST_PATH_IMAGE158
，其中,
Figure DEST_PATH_IMAGE160
分别为松动件跌落位置所在网格与定位位置所在网格之间的水平间隔网格数和垂直网格间隔数。
相对误差：
Figure DEST_PATH_IMAGE162
，其中为矢量误差的模，
Figure 636746DEST_PATH_IMAGE094
为网格的面积，为以三个定位传感器为顶点构成的三角型的面积。
试验结果对本发明提出的松动件定位方法进行了验证，并且还将本发明提出的方法与其他方法进行了比较分析。
通过试验数据分析，试验定位分析结果如表1~4所示。
 
表1  小波尺度2下分析定位结果
Figure DEST_PATH_IMAGE166
表2 小波尺度8下分析定位结果(1)
Figure DEST_PATH_IMAGE168
表3  小波尺度8下分析定位结果(2)
表4 小波尺度20下分析定位结果
Figure DEST_PATH_IMAGE172
由表1~4可以看出,所有钢球跌落位置的定位平均相对误差最小为2.60%，最大为13.79%。根据矢量误差可以看出，定位网格与松动件实际跌落网格的间隔在1个网格左右。说明该方法对不同质量松动件都有较好的定位效果，大质量的松动件定位效果最好。
以175g不含噪声和880g加噪声（信噪比为5dB）钢球为例，分析不同定位方法的估计效果。如表5和表6所示。信号经相同尺度
Figure 158306DEST_PATH_IMAGE032
连续分解后，再经不同时间差估算方法确定时间差，最后都由查表法定位。
表5 175g-尺度8
Figure DEST_PATH_IMAGE174
表6  880g-尺度8(5dB)
Figure DEST_PATH_IMAGE176
由表5和表6可以看出，本发明提出的方法与其他两种方法相比定位精度要高，且抗干扰能力较强。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围不应当被视为仅限于实施例所所陈述的具体形式，本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。

Claims (1)

1.基于小波变换和希尔伯特变换的核电站松动部件定位方法，其特征在于包括以下步骤：
步骤(1)网格划分
    根据定位精度要求，在被撞击物上划分网格，定位精度要求越高，网格划分越小;网格长宽为 
Figure 2013100314165100001DEST_PATH_IMAGE002
，然后对划分好的网格编号；
步骤(2)标定信号采集
通过数据采集卡、电荷放大器和加速度传感器来实现标定信号的采集，将三个加速度传感器安装在被撞击物上，其中加速度传感器采集力锤敲击信号后，经电荷放大器进入数据采集卡；每个通道的标定信号包括力锤敲击信号和本底噪声；
所述的力锤敲击信号为力锤敲击不同网格内不同点所得的信号；
步骤(3)标定信号
Figure 599832DEST_PATH_IMAGE004
的连续小波变换
对获取的标定信号
Figure 237487DEST_PATH_IMAGE004
判断幅值，以最大幅值前0.03s作为开始端，取连续0.3s长信号数据，然后对
Figure 668731DEST_PATH_IMAGE006
进行连续小波变换,获得处理后信号
Figure 2013100314165100001DEST_PATH_IMAGE008
，连续小波变换的尺度分别取2、8、14、20；
步骤(4)计算到达时间
对步骤(3)中处理得到的信号
Figure 50033DEST_PATH_IMAGE008
用希尔伯特变换法进行分析，确定信号
Figure 889813DEST_PATH_IMAGE008
的起振点时刻，第一通道、第二通道、第三通道的起振点时刻分别为
Figure 2013100314165100001DEST_PATH_IMAGE010
、
Figure 2013100314165100001DEST_PATH_IMAGE012
、
Figure 2013100314165100001DEST_PATH_IMAGE014
；
步骤(5)建立定位数据库
根据步骤(4)中得到的不同通道的起振点时刻、、
Figure 8095DEST_PATH_IMAGE014
,计算第二通道、第三通道与第一通道的到达时间差
Figure 2013100314165100001DEST_PATH_IMAGE016
、
Figure 2013100314165100001DEST_PATH_IMAGE018
；计算得出同一网格内不同敲击点的多次敲击结果的时间差平均值
Figure 2013100314165100001DEST_PATH_IMAGE020
、
Figure 2013100314165100001DEST_PATH_IMAGE022
为；然后以网格编号，时间差平均值
Figure 390797DEST_PATH_IMAGE020
和
Figure 42358DEST_PATH_IMAGE022
为查询表的属性，建立定位查询表grid_tbl；
步骤(6)现场数据采集
通过数据采集卡、电荷放大器和加速度传感器来实现现场数据信号采集，现场数据信号
Figure 2013100314165100001DEST_PATH_IMAGE024
包括松动件跌落时的冲击信号和环境背景噪声；
步骤(7)现场数据信号
Figure 212308DEST_PATH_IMAGE024
的傅里叶变换
对获取的现场数据信号判断幅值，以最大幅值前0.03s作为开始端，取0.3s长信号数据
Figure 2013100314165100001DEST_PATH_IMAGE026
，并对
Figure 982742DEST_PATH_IMAGE026
进行傅里叶变换；根据信号频谱图确定
Figure 133100DEST_PATH_IMAGE026
的主频率
Figure 2013100314165100001DEST_PATH_IMAGE028
，然后根据公式：，求得小波尺度
Figure 2013100314165100001DEST_PATH_IMAGE032
,式中
Figure 2013100314165100001DEST_PATH_IMAGE034
为采样周期；
Figure 2013100314165100001DEST_PATH_IMAGE036
为给定小波的中心频率，该中心频率为归一化后的数值；
步骤(8)对
Figure 885287DEST_PATH_IMAGE026
进行连续小波变换
将步骤(3)中连续小波变换尺度2、8、14、20与步骤(7)中求得的小波尺度
Figure 146504DEST_PATH_IMAGE032
作比较，取最接近
Figure 815383DEST_PATH_IMAGE032
的作为新连续小波变换尺度M，然后以尺度M对步骤(7)中得到的信号
Figure 2013100314165100001DEST_PATH_IMAGE038
进行连续小波变换处理得到信号
Figure 2013100314165100001DEST_PATH_IMAGE040
；
步骤(9)希尔伯特变换
对步骤(8)中处理得到的信号
Figure 527573DEST_PATH_IMAGE040
用希尔伯特变换法进行分析，确定该信号的起振点时刻，第一通道、第二通道、第三通道的起振点时刻分别为
Figure 2013100314165100001DEST_PATH_IMAGE042
、
Figure 2013100314165100001DEST_PATH_IMAGE044
、
Figure 2013100314165100001DEST_PATH_IMAGE046
；
步骤(10)查表定位
根据步骤(9)中确定的各通道的起振时刻，求得第二通道、第三通道与第一通道的信号到达时间差
Figure 2013100314165100001DEST_PATH_IMAGE048
、
Figure 2013100314165100001DEST_PATH_IMAGE050
；然后根据步骤(8)中确定的小波尺度M，找到M尺度下的定位查询表grid_tbl；用到达时间差
Figure 203274DEST_PATH_IMAGE048
、
Figure 504067DEST_PATH_IMAGE050
与定位查询表grid_tbl中平均值的
Figure 89769DEST_PATH_IMAGE020
、
Figure 254034DEST_PATH_IMAGE022
分别作比较 ，求其时间距离
Figure 2013100314165100001DEST_PATH_IMAGE052
，时间距离中的最小值所对应的网格编号
Figure 2013100314165100001DEST_PATH_IMAGE054
即为定位网格N；
    所述的步骤(3)和步骤(8)中的连续小波变换的公式为
Figure 2013100314165100001DEST_PATH_IMAGE056
其中，
Figure 2013100314165100001DEST_PATH_IMAGE058
为步骤(3)中的
Figure 49745DEST_PATH_IMAGE008
，为步骤(8)中的
Figure 427636DEST_PATH_IMAGE026
；为小波函数，此处取小波函数为db2；
Figure 2013100314165100001DEST_PATH_IMAGE062
为
Figure 153016DEST_PATH_IMAGE060
的复共轭；
Figure 2013100314165100001DEST_PATH_IMAGE064
为尺度因子；为时间平移因子，即为经小波尺度
Figure 304774DEST_PATH_IMAGE032
连续变换后的信号；
所述的步骤(4)和步骤(9)中采用的希尔伯特变换到达时差定位法，是在时域通过求得信号的希尔伯特包络线，根据包络线的峰值来确定信号起振点时刻，包括以下步骤：
(1)求输入信号
Figure 2013100314165100001DEST_PATH_IMAGE070
的希尔伯特变换
Figure 2013100314165100001DEST_PATH_IMAGE072
，根据定义：
       
Figure 2013100314165100001DEST_PATH_IMAGE074
 ，其中，即为
Figure 838490DEST_PATH_IMAGE070
经希尔伯特变换后的信号；
(2)以信号为实部，其希尔伯特变换
Figure 357513DEST_PATH_IMAGE072
为虚部，构成一个新的函数如式： 
Figure 2013100314165100001DEST_PATH_IMAGE076
 
其中， 
Figure 2013100314165100001DEST_PATH_IMAGE078
是新函数的幅值函数；为相位函数，则
Figure DEST_PATH_IMAGE082
即为信号的包络线函数；
(3)对获取的包络线函数
Figure 608814DEST_PATH_IMAGE082
判断幅值，取信号
Figure 285783DEST_PATH_IMAGE082
幅值突然变大前0.1s长度信号作为背景噪声，取背景噪声的最大幅值作为阈值
Figure DEST_PATH_IMAGE084
；
(4)对
Figure 520455DEST_PATH_IMAGE082
求峰值，比较点数一般设置为5，即取
Figure 33083DEST_PATH_IMAGE082
中大于阈值
Figure 240073DEST_PATH_IMAGE084
且大于与该点相邻的5个点的所有值，则此点所对应的值
Figure DEST_PATH_IMAGE086
即为的一个峰值，取求得的第一个峰值所对应的时间点t即为信号的起振时刻；
所述的步骤(10)中时间距离的计算公式为：
Figure DEST_PATH_IMAGE088
，
Figure DEST_PATH_IMAGE090
其中，
Figure 763907DEST_PATH_IMAGE054
为网格编号，
Figure DEST_PATH_IMAGE092
为网格划分总数；的最小值所对应的
Figure 50531DEST_PATH_IMAGE054
即为定位网格N。
CN201310031416.5A 2013-01-25 2013-01-25 基于小波变换和希尔伯特变换的核电站松动部件定位方法 Active CN103137225B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310031416.5A CN103137225B (zh) 2013-01-25 2013-01-25 基于小波变换和希尔伯特变换的核电站松动部件定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310031416.5A CN103137225B (zh) 2013-01-25 2013-01-25 基于小波变换和希尔伯特变换的核电站松动部件定位方法

Publications (2)

Publication Number Publication Date
CN103137225A true CN103137225A (zh) 2013-06-05
CN103137225B CN103137225B (zh) 2015-09-30

Family

ID=48496927

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310031416.5A Active CN103137225B (zh) 2013-01-25 2013-01-25 基于小波变换和希尔伯特变换的核电站松动部件定位方法

Country Status (1)

Country Link
CN (1) CN103137225B (zh)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103712758A (zh) * 2013-12-12 2014-04-09 西安交通大学 利用响应信息判定粘弹夹层结构松动程度的方法
CN104111108A (zh) * 2014-07-16 2014-10-22 云南电力试验研究院（集团）有限公司电力研究院 一种旋转机械扭振冲击信号特征提取方法
CN104408303A (zh) * 2014-11-20 2015-03-11 西安交通大学 一种基于数据匹配的lpms质量估计方法
CN105387883A (zh) * 2015-10-28 2016-03-09 国家电网公司 一种电站电力部件松动定位系统及定位方法
CN105427904A (zh) * 2015-11-13 2016-03-23 陕西卫峰核电子有限公司 一种基于核电站松动部件定位方法
CN105702305A (zh) * 2016-02-03 2016-06-22 杭州电子科技大学 核电站松动部件多目标三维动态跟踪方法
CN105865654A (zh) * 2016-03-23 2016-08-17 东南大学 一种声波测温信号的选取方法及锅炉测温方法
CN106409351A (zh) * 2016-09-13 2017-02-15 中广核工程有限公司 一种核电厂松动部件振动信号的处理方法及装置
CN107544337A (zh) * 2017-09-19 2018-01-05 中国核动力研究设计院 一种松脱部件监测系统触发数据智能分类的方法
CN108921082A (zh) * 2018-06-28 2018-11-30 哈尔滨工程大学 基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法
CN109003686A (zh) * 2018-05-25 2018-12-14 中广核工程有限公司 核电站松脱部件的定位方法、装置、计算机设备及存储介质
CN109581516A (zh) * 2018-11-29 2019-04-05 桂林理工大学 曲波域统计量自适应阈值探地雷达数据去噪方法及系统
CN109975136A (zh) * 2019-04-09 2019-07-05 西南交通大学 一种基于小波包分析的钢框架结构损伤识别方法
CN110658053A (zh) * 2019-08-29 2020-01-07 中国空间技术研究院 一种基于小波变换的卫星组件冲击试验条件制定系统及方法
CN111827370A (zh) * 2019-04-17 2020-10-27 福建农林大学 基于小波系数相位角变化的桩基损伤位置判别方法
CN112612004A (zh) * 2020-12-14 2021-04-06 中广核研究院有限公司 部件定位方法、装置、计算机设备和存储介质
CN113688351A (zh) * 2021-08-06 2021-11-23 北京每日优鲜电子商务有限公司 用于检测物品重量的方法、装置、电子设备和可读介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4955269A (en) * 1988-02-04 1990-09-11 Westinghouse Electric Corp. Turbine blade fatigue monitor
CN102520075A (zh) * 2011-12-30 2012-06-27 哈尔滨商业大学 基于谐波小波包的直升机部件声发射信号传播特性的测试方法
CN102539150A (zh) * 2012-01-17 2012-07-04 电子科技大学 基于连续小波变换的旋转机械部件的自适应故障诊断方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4955269A (en) * 1988-02-04 1990-09-11 Westinghouse Electric Corp. Turbine blade fatigue monitor
CN102520075A (zh) * 2011-12-30 2012-06-27 哈尔滨商业大学 基于谐波小波包的直升机部件声发射信号传播特性的测试方法
CN102539150A (zh) * 2012-01-17 2012-07-04 电子科技大学 基于连续小波变换的旋转机械部件的自适应故障诊断方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GEE YONG PARK ET AL.: "An estimation method for impact location of loose parts", 《PROGRESS IN NUCLAER ENERGY》 *
方力先 陈仲仪: "小波变换在核电站一回路松动件定位的应用研究", 《核科学与工程》 *
杨将新 等: "基于希尔伯特变换的核电站松动部件定位方法", 《机械工程学报》 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103712758B (zh) * 2013-12-12 2016-03-30 西安交通大学 利用响应信息判定粘弹夹层结构松动程度的方法
CN103712758A (zh) * 2013-12-12 2014-04-09 西安交通大学 利用响应信息判定粘弹夹层结构松动程度的方法
CN104111108A (zh) * 2014-07-16 2014-10-22 云南电力试验研究院（集团）有限公司电力研究院 一种旋转机械扭振冲击信号特征提取方法
CN104408303A (zh) * 2014-11-20 2015-03-11 西安交通大学 一种基于数据匹配的lpms质量估计方法
CN104408303B (zh) * 2014-11-20 2017-09-12 西安交通大学 一种基于数据匹配的lpms质量估计方法
CN105387883A (zh) * 2015-10-28 2016-03-09 国家电网公司 一种电站电力部件松动定位系统及定位方法
CN105427904A (zh) * 2015-11-13 2016-03-23 陕西卫峰核电子有限公司 一种基于核电站松动部件定位方法
CN105702305A (zh) * 2016-02-03 2016-06-22 杭州电子科技大学 核电站松动部件多目标三维动态跟踪方法
CN105702305B (zh) * 2016-02-03 2019-01-15 杭州电子科技大学 核电站松动部件多目标三维动态跟踪方法
CN105865654A (zh) * 2016-03-23 2016-08-17 东南大学 一种声波测温信号的选取方法及锅炉测温方法
CN105865654B (zh) * 2016-03-23 2018-07-27 东南大学 一种声波测温信号的选取方法及锅炉测温方法
CN106409351A (zh) * 2016-09-13 2017-02-15 中广核工程有限公司 一种核电厂松动部件振动信号的处理方法及装置
CN107544337A (zh) * 2017-09-19 2018-01-05 中国核动力研究设计院 一种松脱部件监测系统触发数据智能分类的方法
CN109003686A (zh) * 2018-05-25 2018-12-14 中广核工程有限公司 核电站松脱部件的定位方法、装置、计算机设备及存储介质
CN108921082A (zh) * 2018-06-28 2018-11-30 哈尔滨工程大学 基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法
CN108921082B (zh) * 2018-06-28 2021-11-12 哈尔滨工程大学 基于经验小波和希尔伯特变换的柴油机燃烧共振频率提取方法
CN109581516A (zh) * 2018-11-29 2019-04-05 桂林理工大学 曲波域统计量自适应阈值探地雷达数据去噪方法及系统
CN109581516B (zh) * 2018-11-29 2020-12-25 桂林理工大学 曲波域统计量自适应阈值探地雷达数据去噪方法及系统
CN109975136A (zh) * 2019-04-09 2019-07-05 西南交通大学 一种基于小波包分析的钢框架结构损伤识别方法
CN109975136B (zh) * 2019-04-09 2021-06-25 西南交通大学 一种基于小波包分析的钢框架结构损伤识别方法
CN111827370A (zh) * 2019-04-17 2020-10-27 福建农林大学 基于小波系数相位角变化的桩基损伤位置判别方法
CN110658053A (zh) * 2019-08-29 2020-01-07 中国空间技术研究院 一种基于小波变换的卫星组件冲击试验条件制定系统及方法
CN112612004A (zh) * 2020-12-14 2021-04-06 中广核研究院有限公司 部件定位方法、装置、计算机设备和存储介质
CN112612004B (zh) * 2020-12-14 2024-04-09 中广核研究院有限公司 部件定位方法、装置、计算机设备和存储介质
CN113688351A (zh) * 2021-08-06 2021-11-23 北京每日优鲜电子商务有限公司 用于检测物品重量的方法、装置、电子设备和可读介质

Also Published As

Publication number Publication date
CN103137225B (zh) 2015-09-30

Similar Documents

Publication Publication Date Title
CN103137225B (zh) 基于小波变换和希尔伯特变换的核电站松动部件定位方法
CN103137224B (zh) 基于小波能量谱的核电站松动部件质量估计方法
CN104376881B (zh) 基于希尔伯特变换和数据筛选的核电站松动件定位方法
CN109812390A (zh) 一种风力发电机组的叶片净空监测方法
CN105702305B (zh) 核电站松动部件多目标三维动态跟踪方法
Liu et al. A new method for arrival time determination of impact signal based on HHT and AIC
CN101586997A (zh) 一种拉索振动基频的计算方法
CN108627741A (zh) 一种基于故障指示器的行波—阻抗法双端带支路配电网故障定位方法
CN109946763A (zh) 一种基于波群识别的畸形波短期预报方法
CN103149552B (zh) 一种多普勒天气雷达径向速度场的退模糊方法
CN103792513B (zh) 一种雷声定位系统及方法
CN107607092B (zh) 一种基于无人机遥感的波浪与浮体运动预报系统
CN102628955A (zh) 一种地震预警震级获取方法
WO2022179247A1 (zh) 基于全索系索力测试和误差自适应分析的损伤判定方法
CN110134976A (zh) 一种机载激光测深信号提取方法及系统
CN104217112A (zh) 一种基于多类型信号的电力系统低频振荡分析方法
CN105116396A (zh) 一种连续波雷达多普勒回波检测方法
CN102879714A (zh) 变压器局部放电检测定位方法
CN110007299A (zh) 一种基于混合坐标伪谱技术的微弱目标检测跟踪方法
CN105676018A (zh) 一种基于稳态学习机制的电网低频振荡模糊综合预警方法
CN113970739A (zh) 一种空管一次雷达自适应风电场杂波识别和抑制方法
CN107843406B (zh) 基于脉动压力相关函数的空腔模态波运动特征确定方法
CN102546052B (zh) 一种自适应门限控制信号参数统计测量方法
CN104215939A (zh) 一种融合广义对称结构信息的知识辅助空时自适应处理方法
CN113344275A (zh) 一种基于lstm模型的浮式平台波浪爬升在线预报方法

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
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 310018 No. 2 street, Xiasha Higher Education Zone, Hangzhou, Zhejiang

Patentee after: HANGZHOU DIANZI University

Patentee after: SHANGHAI NUCLEAR ENGINEERING RESEARCH & DESIGN INSTITUTE Co.,Ltd.

Address before: 310018 No. 2 street, Xiasha Higher Education Zone, Hangzhou, Zhejiang

Patentee before: HANGZHOU DIANZI University

Patentee before: SHANGHAI NUCLEAR ENGINEERING RESEARCH & DESIGN INSTITUTE

Address after: 310018 No. 2 street, Xiasha Higher Education Zone, Hangzhou, Zhejiang

Patentee after: HANGZHOU DIANZI University

Patentee after: Shanghai Nuclear Engineering Research and Design Institute Co.,Ltd.

Address before: 310018 No. 2 street, Xiasha Higher Education Zone, Hangzhou, Zhejiang

Patentee before: HANGZHOU DIANZI University

Patentee before: SHANGHAI NUCLEAR ENGINEERING RESEARCH & DESIGN INSTITUTE Co.,Ltd.