CN103137224A - 基于小波能量谱的核电站松动部件质量估计方法 - Google Patents

基于小波能量谱的核电站松动部件质量估计方法 Download PDF

Info

Publication number
CN103137224A
CN103137224A CN2013100313800A CN201310031380A CN103137224A CN 103137224 A CN103137224 A CN 103137224A CN 2013100313800 A CN2013100313800 A CN 2013100313800A CN 201310031380 A CN201310031380 A CN 201310031380A CN 103137224 A CN103137224 A CN 103137224A
Authority
CN
China
Prior art keywords
signal
quality
wavelet
yardstick
impact
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
CN2013100313800A
Other languages
English (en)
Other versions
CN103137224B (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 CN201310031380.0A priority Critical patent/CN103137224B/zh
Publication of CN103137224A publication Critical patent/CN103137224A/zh
Application granted granted Critical
Publication of CN103137224B publication Critical patent/CN103137224B/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.模拟核反应堆中松动件跌落时的冲击信号
Figure 2013100313800100004DEST_PATH_IMAGE001
并采集；3.对获取的冲击信号
Figure 431125DEST_PATH_IMAGE001
判断幅值；4.根据计算的信号的小波能量谱找出该图的尺度峰值
Figure DEST_PATH_IMAGE003
；5.用线性插值法建立尺度峰值函数
Figure 705297DEST_PATH_IMAGE004
；6.采集现场信号
Figure DEST_PATH_IMAGE005
；7.对获取的现场信号
Figure 437235DEST_PATH_IMAGE005
判断幅值；8.计算信号
Figure 518324DEST_PATH_IMAGE006
的小波能量谱，根据小波能量谱图找出该图尺度峰值
Figure DEST_PATH_IMAGE007
；9.将尺度峰值
Figure 258746DEST_PATH_IMAGE007
带入5建立的尺度峰值函数
Figure 21166DEST_PATH_IMAGE004
，求得松动件的质量估计值，然后对多个通道中的松动件质量估计值求平均，得到最终的松动件质量结果。本发明具有质量估计误差小、一致性好和抗干扰能力强的优点。

Description

基于小波能量谱的核电站松动部件质量估计方法
技术领域
本发明属于核工程技术领域，尤其涉及一种基于小波能量谱的核电站松动部件质量估计方法，用于估计核电站松动部件质量的质量。
技术背景
核电站一回路中松动部件的出现，就意味着反应堆系统安全性能和稳定性能的削弱。松动部件检测系统（LPMS）是核电站一回路最近本的安全检测工具之一。质量估计作为LPMS中的重要组成部分，准确的估算出松动部件的质量能够为判别松动部件类型提供重要依据，为故障诊断提供便利。对于提高核电站系统运行的安全性和稳定性起到重要作用。
现有的松动部件质量估计方法相关文献有：
[1] MAYO W . loose part mass and energy estimation [J] .Progress in Nuclear Energy, 1999, 34(3): 263-282.
一种基于碰撞信号高低频段能量比的松动件质量估计方法。该方法质量估计误差较大，一致性差等问题，尤其是抗干扰能力差。
[2]模式识别技术用于核电站松动件质量估计模型研究[J]，核动力工程，2001, 22（5）: 465-470
一种利用模式识别理论和人工神经网络理论实现对核电站松动件质量进行估计的数学模型。利用模式识别理论的特征提取技术，实现对松动件碰撞信号的特征压缩，从而形成比较少的特征空间维数，经过神经网络的学习过程，实现对核电站松动件质量大小的估计。该方法能够得到较高的质量估计精度，但需要大量的试验数据对其进行训练，且一致性较差。 
[3] 曹衍龙，程实，杨将新，郑华文，何元峰，基于小波包的质量估计方法[J]，机械工程学报，2010 .46(22):1-5
.     一种基于小波包的松动件质量估计方法。该方法通过小波包分解获取各频段小波包系数，计算各频段小波包系数的能量潜，并进行归化，得到特征值矢量；以此特征值矢量作为输入参数，以松动件质量作为输出参数，应用支持向量机((Support vector machine,  SVM)实现松动件质量估计。
发明内容
本发明的目的是为了克服现有质量估计方法存在的误差大、一致性差问题，提出了一种基于小波能量谱的核电站松动部件质量估计方法。本发明具有估计误差小、一致性好和抗干扰能力强等特点，具有较好的工程实际应用价值。
本发明解决其技术问题所采用的技术方案包括以下步骤：
步骤(1)在核反应堆一回路安装多个加速度传感器，以获取核反应堆中的环境噪和松动件跌落时产生的混杂有环境噪声的冲击信号。
步骤(2)用不同质量钢球敲击钢板来模拟核反应堆中松动件跌落时的冲击信号；通过数据采集卡、电荷放大器和加速度传感器获取钢球敲击钢板的冲击信号                                               
Figure 2013100313800100002DEST_PATH_IMAGE002
，并将这些冲击信号组成质量估计信号库；
所述的冲击信号
Figure 991345DEST_PATH_IMAGE002
为已知质量钢球敲击钢板产生的冲击信号；
步骤(3)对获取的冲击信号判断幅值，以最大幅值前0.03s为开始端，取0.3s长信号数据，对信号数据
Figure 829168DEST_PATH_IMAGE004
进行连续小波变换处理得到信号
Figure 2013100313800100002DEST_PATH_IMAGE006
。
步骤(4)计算信号
Figure 896481DEST_PATH_IMAGE006
的小波能量谱，根据小波能量谱图找出该图能量峰值所对应的尺度，即尺度峰值
Figure 2013100313800100002DEST_PATH_IMAGE008
,然后统计不同质量钢球的冲击信号所对应的尺度峰值的均值
Figure 2013100313800100002DEST_PATH_IMAGE010
。
步骤(5)用线性插值法建立不同质量钢球与尺度峰值均值之间的尺度峰值函数；
步骤(6)通过数据采集卡、电荷放大器和加速度传感器来实现现场信号采集，现场信号
Figure 730094DEST_PATH_IMAGE014
包括松动件跌落时的冲击信号和环境背景噪声；
步骤(7)对获取的现场信号
Figure 552557DEST_PATH_IMAGE014
判断幅值，以最大幅值出现前0.03s开始，取0.3s长信号数据
Figure 2013100313800100002DEST_PATH_IMAGE016
，对信号
Figure 491694DEST_PATH_IMAGE016
进行连续小波变换处理得到信号。
步骤(8)计算信号
Figure 831540DEST_PATH_IMAGE018
的小波能量谱，根据小波能量谱图找出该图尺度峰值
Figure 2013100313800100002DEST_PATH_IMAGE020
。 
步骤(9)将尺度峰值带入步骤(5)建立尺度峰值函数
Figure 165886DEST_PATH_IMAGE012
，求得松动件的质量估计值，然后对多个通道中的松动件质量估计值求平均，得到最终的松动件质量结果。
所述的步骤(3)和步骤(7)中连续小波变换的公式为：
Figure 2013100313800100002DEST_PATH_IMAGE024
,
其中，为步骤(3)中的
Figure 933597DEST_PATH_IMAGE004
，为步骤(7)中的
Figure 670609DEST_PATH_IMAGE016
；
Figure 2013100313800100002DEST_PATH_IMAGE028
为小波函数，此处取小波函数为db2；为的复共轭；
Figure 2013100313800100002DEST_PATH_IMAGE032
为尺度因子；
Figure 2013100313800100002DEST_PATH_IMAGE034
为时间平移因子。
Figure 2013100313800100002DEST_PATH_IMAGE036
即为经小波尺度
Figure 108992DEST_PATH_IMAGE032
连续变换后的信号。
所述的步骤(4)和步骤(8)中小波能量谱的公式为：
Figure 2013100313800100002DEST_PATH_IMAGE038
，
式中，
Figure 2013100313800100002DEST_PATH_IMAGE040
为步骤(4)中的，为步骤(8)中的
Figure 515756DEST_PATH_IMAGE018
；
Figure 2013100313800100002DEST_PATH_IMAGE042
称为在尺度上的小波能量谱。则被分析的时间信号
Figure 2013100313800100002DEST_PATH_IMAGE044
的能量分布为：
Figure 2013100313800100002DEST_PATH_IMAGE048
式中，
Figure 2013100313800100002DEST_PATH_IMAGE050
反映了信号经小波变换后能量随尺度的变化情况。
所述的步骤(5)中的线性插值法为：设质量数组为:
Figure 2013100313800100002DEST_PATH_IMAGE052
,与设质量数组为
Figure 2013100313800100002DEST_PATH_IMAGE054
对应的小波尺度均值数组为:
Figure 2013100313800100002DEST_PATH_IMAGE056
，变量
Figure 120492DEST_PATH_IMAGE020
为输入小波尺度，变量
Figure 52676DEST_PATH_IMAGE022
为质量估计值。则可以建立尺度峰值函数：
Figure 2013100313800100002DEST_PATH_IMAGE058
 ，(
Figure 2013100313800100002DEST_PATH_IMAGE060
)。
本发明有益效果如下：
先通过对质量估计信号库中已知质量松动件敲击钢板产生的冲击信号进行小波连续变换，然后计算信号小波能量谱，根据信号能量尺度图找到小波能量峰值处所对应的小波尺度，统计多次同质量钢球冲击信号的尺度峰值，并求平均值。根据求得的各不同质量钢球的尺度峰值均值与质量，并通过线性插值法建立碰撞质量尺度峰值函数。最后通过质量尺度峰值函数来实现核电站松动件的质量估计。小波能量谱方法能很好地反映出冲击质量的大小变化，且各通道特征具有较好的一致性。该方法能以能量图的形式将冲击质量表示出来，直观明了，从能量图上就能大致判别冲击质量。抗干扰能力较强。在5dB噪声干扰情况下，仍然保持冲击信号的主要特征不变；在低于0dB时，能量图主峰主要表现为噪声信号能量特性，而冲击信号能量峰值变为次峰。若噪声是已知的，则仍能实现质量估计。
本发明具有质量估计误差小、一致性好和抗干扰能力强的优点，具有较好的工程实际应用价值。
附图说明
图1 本发明同一质量冲击信号三通道小波能量谱图；
图2 本发明不同质量冲击信号小波能量谱图；
图3本发明传感器安装位置图；
图4 本发明中1.8kg信号与1.8kg-10dB信号时域图；
图5本发明中1.8kg信号与1.8kg-10dB信号小波能量谱图；
图6本发明中冲击质量与尺度峰值均值关系；
图7 本发明的流程图。
具体实施方式
下面结合附图和实施例对 本发明作进一步说明。
如图7所示，基于小波能量谱的核电站松动部件质量估计方法，包括以下步骤：
步骤(1)在核反应堆一回路安装多个加速度传感器，以获取核反应堆中的环境噪和松动件跌落时产生的混杂有环境噪声的冲击信号。
步骤(2)用不同质量钢球敲击钢板来模拟核反应堆中松动件跌落时的冲击信号；通过数据采集卡、电荷放大器和加速度传感器获取钢球敲击钢板的冲击信号，并将这些冲击信号组成质量估计信号库；
所述的冲击信号
Figure 211234DEST_PATH_IMAGE002
为已知质量钢球敲击钢板产生的冲击信号；
步骤(3)对获取的冲击信号
Figure 161872DEST_PATH_IMAGE002
判断幅值，以最大幅值前0.03s为开始端，取0.3s长信号数据
Figure 716482DEST_PATH_IMAGE004
，对信号数据
Figure 272228DEST_PATH_IMAGE004
进行连续小波变换处理得到信号
Figure 214776DEST_PATH_IMAGE006
。
步骤(4)计算信号
Figure 590394DEST_PATH_IMAGE006
的小波能量谱，根据小波能量谱图找出该图能量峰值所对应的尺度，即尺度峰值
Figure 683115DEST_PATH_IMAGE008
,然后统计不同质量钢球的冲击信号所对应的尺度峰值的均值
Figure 155684DEST_PATH_IMAGE010
。
步骤(5)用线性插值法建立不同质量钢球与尺度峰值均值之间的尺度峰值函数
Figure 472396DEST_PATH_IMAGE012
；
步骤(6)通过数据采集卡、电荷放大器和加速度传感器来实现现场信号
Figure 397627DEST_PATH_IMAGE014
采集，现场信号
Figure 296968DEST_PATH_IMAGE014
包括松动件跌落时的冲击信号和环境背景噪声；
步骤(7)对获取的现场信号
Figure 358465DEST_PATH_IMAGE014
判断幅值，以最大幅值出现前0.03s开始，取0.3s长信号数据
Figure 580499DEST_PATH_IMAGE016
，对信号
Figure 196288DEST_PATH_IMAGE016
进行连续小波变换处理得到信号。
步骤(8)计算信号
Figure 77974DEST_PATH_IMAGE018
的小波能量谱，根据小波能量谱图找出该图尺度峰值
Figure 267647DEST_PATH_IMAGE020
。 
步骤(9)将尺度峰值
Figure 370732DEST_PATH_IMAGE020
带入步骤(5)建立尺度峰值函数
Figure 874526DEST_PATH_IMAGE012
，求得松动件的质量估计值，然后对多个通道中的松动件质量估计值求平均，得到最终的松动件质量结果。
所述的步骤(3)和步骤(7)中连续小波变换的公式为：
Figure 208872DEST_PATH_IMAGE024
,
其中，为步骤(3)中的
Figure 169055DEST_PATH_IMAGE004
，为步骤(7)中的
Figure 728825DEST_PATH_IMAGE016
；
Figure 525879DEST_PATH_IMAGE028
为小波函数，此处取小波函数为db2；
Figure 72398DEST_PATH_IMAGE030
为的复共轭；
Figure 459834DEST_PATH_IMAGE032
为尺度因子；
Figure 427790DEST_PATH_IMAGE034
为时间平移因子。即为经小波尺度
Figure 907630DEST_PATH_IMAGE032
连续变换后的信号。
所述的步骤(4)和步骤(8)中小波能量谱的公式为：
Figure 976080DEST_PATH_IMAGE038
，
式中，
Figure 724725DEST_PATH_IMAGE040
为步骤(4)中的
Figure 308153DEST_PATH_IMAGE006
，为步骤(8)中的；
Figure 746404DEST_PATH_IMAGE042
称为在尺度
Figure 719478DEST_PATH_IMAGE032
上的小波能量谱。则被分析的时间信号
Figure 790202DEST_PATH_IMAGE044
的能量分布为：
Figure 781292DEST_PATH_IMAGE046
Figure 621072DEST_PATH_IMAGE048
式中，反映了信号经小波变换后能量随尺度的变化情况。
所述的步骤(5)中的线性插值法为：设质量数组为:
Figure 597435DEST_PATH_IMAGE052
,与其对应的小波尺度均值为:
Figure 126637DEST_PATH_IMAGE056
，变量
Figure 820923DEST_PATH_IMAGE020
为输入小波尺度，变量
Figure 675747DEST_PATH_IMAGE022
为质量估计值。则可以建立分段函数：
Figure 726059DEST_PATH_IMAGE058
 ，(
Figure 274852DEST_PATH_IMAGE060
)。
如图1所示，从同一质量钢球冲击信号三通道的小波能量谱图中不难看出，三通道传感器接收到的能量分布具有一致性，且能量峰值具有较一致的尺度，这与Hertz理论指出的碰撞主频率是相吻合的。
如图2所示，从不同质量钢球的冲击信号归一化后的小波能量谱图中可看出，随着冲击质量的增大，小波能量谱向右移动。这些明显的特征变化为质量估计提供了重要依据，即可以根据尺度峰值随着冲击物质量的变化而发生转移的情况来估计松动部件的质量。
本发明先通过对质量估计信号库中已知质量松动件敲击钢板产生的冲击信号进行小波连续变换，然后计算信号小波能量谱，根据信号能量尺度图找到小波能量峰值处所对应的小波尺度，统计多次同质量钢球冲击信号的尺度峰值，并求平均值。根据求得的各不同质量钢球的尺度峰值均值与质量，并通过线性插值法建立碰撞质量尺度峰值函数。最后通过质量尺度峰值函数来实现核电站松动件的质量估计。小波能量谱方法能很好地反映出冲击质量的大小变化，且各通道特征具有较好的一致性。该方法能以能量图的形式将冲击质量表示出来，直观明了，从能量图上就能大致判别冲击质量。抗干扰能力较强。在5dB噪声干扰情况下，仍然保持冲击信号的主要特征不变；在低于0dB时，能量图主峰主要表现为噪声信号能量特性，而冲击信号能量峰值变为次峰。若噪声是已知的，则仍能实现质量估计。
实施例
下面结合实施例对本发明作进一步说明：
本发明的效果可以通过平板试验的实验结果分析说明：
1.试验条件
本试验的试验平台由测试对象、传感器、电荷放大器、数据采集卡和计算机构成。测试对象包括：钢球、卧式锅炉和钢板及其支撑。钢板尺寸为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_IMAGE062
               
其中，
Figure DEST_PATH_IMAGE064
表示噪声幅值取绝对值后求平均，表示信号幅值取绝对值后求最大值。
4.仿真结果
实验结果对本发明提出的质量估计方法进行了验证。实验结果如附图4和附图5所示，附图5为1.8kg钢球无噪信号和加10dB噪声信号的小波能量图，从图中可看出含噪信号的小波能量主峰值表现为噪声信号特征，冲击信号特征表现为次峰值。所以只要背景噪声的频率结构是已知的，就可以根据噪声信号特征，结合信噪比情况对估值结果进行适当的修正，从而提高质量估计的准确度。
本发明用如下方法来计算绝对误差和相对误差：
绝对误差：
Figure DEST_PATH_IMAGE068
，其中
Figure DEST_PATH_IMAGE070
为松动件质量真实值与估计值之间的差值。
相对误差：
Figure DEST_PATH_IMAGE072
，其中
Figure DEST_PATH_IMAGE074
为松动件质量真实值。
试验结果对本发明提出的松动件质量估计方法进行了验证，通过试验数据分析，建立尺度峰值函数，通过尺度峰值函数来实现松松部件的质量估计。附图6为松动件质量与尺度之间的对应关系，表1为各冲击质量对应的尺度均值，从附图6和表1中不难看出，小质量的尺度峰值变化相对明显，但随着冲击质量的增大，峰值变化逐渐变得不明显，这是因为冲击质量大的信号，其能量主要集中在低频范围内，在频谱上很难区分。
表1 不同尺度峰值所对应的质量估计值
Figure DEST_PATH_IMAGE076
表2 质量估计误差
Figure DEST_PATH_IMAGE078
表2为质量估计的误差，从表中可以看出，质量估计误差最大为19.4%，最小为9.1%，质量估计误差基本在15%左右。所以本发明所提出的方法具有较高的准确率。且抗干扰能力较强，在背景噪声频谱结构已知的情况下，在信噪比为0dB时仍能较准确的进行松动件质量估计。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围不应当被视为仅限于实施例所所陈述的具体形式，本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。

Claims (1)

1. 基于小波能量谱的核电站松动部件质量估计方法，其特征在于包括以下步骤：
步骤(1)在核反应堆一回路安装多个加速度传感器，以获取核反应堆中的环境噪和松动件跌落时产生的混杂有环境噪声的冲击信号；
步骤(2)用不同质量钢球敲击钢板来模拟核反应堆中松动件跌落时的冲击信号；通过数据采集卡、电荷放大器和加速度传感器获取钢球敲击钢板的冲击信号                                               
Figure 2013100313800100001DEST_PATH_IMAGE002
，并将这些冲击信号组成质量估计信号库；
所述的冲击信号
Figure 741876DEST_PATH_IMAGE002
为已知质量钢球敲击钢板产生的冲击信号；
步骤(3)对获取的冲击信号
Figure 8910DEST_PATH_IMAGE002
判断幅值，以最大幅值前0.03s为开始端，取0.3s长信号数据
Figure 2013100313800100001DEST_PATH_IMAGE004
，对信号数据
Figure 817597DEST_PATH_IMAGE004
进行连续小波变换处理得到信号
Figure 2013100313800100001DEST_PATH_IMAGE006
；
步骤(4)计算信号的小波能量谱，根据小波能量谱图找出该图能量峰值所对应的尺度，即尺度峰值
Figure 2013100313800100001DEST_PATH_IMAGE008
,然后统计不同质量钢球的冲击信号所对应的尺度峰值的均值；
步骤(5)用线性插值法建立不同质量钢球与尺度峰值均值之间的尺度峰值函数
Figure DEST_PATH_IMAGE012
；
步骤(6)通过数据采集卡、电荷放大器和加速度传感器来实现现场信号
Figure DEST_PATH_IMAGE014
采集，现场信号
Figure 171185DEST_PATH_IMAGE014
包括松动件跌落时的冲击信号和环境背景噪声；
步骤(7)对获取的现场信号判断幅值，以最大幅值出现前0.03s开始，取0.3s长信号数据
Figure DEST_PATH_IMAGE016
，对信号
Figure 865788DEST_PATH_IMAGE016
进行连续小波变换处理得到信号
Figure DEST_PATH_IMAGE018
；
步骤(8)计算信号
Figure 326857DEST_PATH_IMAGE018
的小波能量谱，根据小波能量谱图找出该图尺度峰值
Figure DEST_PATH_IMAGE020
； 
步骤(9)将尺度峰值带入步骤(5)建立尺度峰值函数
Figure 608113DEST_PATH_IMAGE012
，求得松动件的质量估计值
Figure DEST_PATH_IMAGE022
，然后对多个通道中的松动件质量估计值求平均，得到最终的松动件质量结果；
所述的步骤(3)和步骤(7)中连续小波变换的公式为：
Figure DEST_PATH_IMAGE024
,
其中，
Figure DEST_PATH_IMAGE026
为步骤(3)中的
Figure 958936DEST_PATH_IMAGE004
，为步骤(7)中的
Figure 274511DEST_PATH_IMAGE016
；为小波函数，此处取小波函数为db2；
Figure DEST_PATH_IMAGE030
为
Figure 750622DEST_PATH_IMAGE028
的复共轭；
Figure DEST_PATH_IMAGE032
为尺度因子；为时间平移因子，即为经小波尺度
Figure 18225DEST_PATH_IMAGE032
连续变换后的信号；
所述的步骤(4)和步骤(8)中小波能量谱的公式为：
，
式中，
Figure DEST_PATH_IMAGE040
为步骤(4)中的
Figure 34723DEST_PATH_IMAGE006
，为步骤(8)中的；
Figure DEST_PATH_IMAGE042
称为在尺度
Figure 914134DEST_PATH_IMAGE032
上的小波能量谱，则被分析的时间信号
Figure DEST_PATH_IMAGE044
的能量分布为：
Figure DEST_PATH_IMAGE046
Figure DEST_PATH_IMAGE048
式中，
Figure DEST_PATH_IMAGE050
反映了信号经小波变换后能量随尺度的变化情况；
所述的步骤(5)中的线性插值法为：设质量数组为:,与设质量数组为
Figure DEST_PATH_IMAGE054
对应的小波尺度均值数组为:
Figure DEST_PATH_IMAGE056
，变量
Figure 815837DEST_PATH_IMAGE020
为输入小波尺度，变量
Figure 636026DEST_PATH_IMAGE022
为质量估计值，则可以建立尺度峰值函数：
Figure DEST_PATH_IMAGE058
 ，(
Figure DEST_PATH_IMAGE060
)。
CN201310031380.0A 2013-01-25 2013-01-25 基于小波能量谱的核电站松动部件质量估计方法 Active CN103137224B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310031380.0A CN103137224B (zh) 2013-01-25 2013-01-25 基于小波能量谱的核电站松动部件质量估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310031380.0A CN103137224B (zh) 2013-01-25 2013-01-25 基于小波能量谱的核电站松动部件质量估计方法

Publications (2)

Publication Number Publication Date
CN103137224A true CN103137224A (zh) 2013-06-05
CN103137224B CN103137224B (zh) 2015-07-15

Family

ID=48496926

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310031380.0A Active CN103137224B (zh) 2013-01-25 2013-01-25 基于小波能量谱的核电站松动部件质量估计方法

Country Status (1)

Country Link
CN (1) CN103137224B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104408303A (zh) * 2014-11-20 2015-03-11 西安交通大学 一种基于数据匹配的lpms质量估计方法
CN104965026A (zh) * 2015-05-25 2015-10-07 北京理工大学 一种使用声发射信号能量特征值预测脆性材料临界失稳的方法
CN105045965A (zh) * 2015-06-26 2015-11-11 西安交通大学 基于仿真数据库的压力容器lpms质量估计方法
CN105225433A (zh) * 2015-08-22 2016-01-06 上海核工程研究设计院 一种核电站松动部件监测报警系统及方法
CN105427904A (zh) * 2015-11-13 2016-03-23 陕西卫峰核电子有限公司 一种基于核电站松动部件定位方法
CN106409351A (zh) * 2016-09-13 2017-02-15 中广核工程有限公司 一种核电厂松动部件振动信号的处理方法及装置
CN106874896A (zh) * 2017-03-31 2017-06-20 阳江核电有限公司 一种核电站一回路部件松脱诊断系统信号特征识别的辅助学习方法及系统
CN107194135A (zh) * 2017-07-21 2017-09-22 中国工程物理研究院核物理与化学研究所 基于改进小波算法的反应堆中子噪声频谱计算方法
CN107544337A (zh) * 2017-09-19 2018-01-05 中国核动力研究设计院 一种松脱部件监测系统触发数据智能分类的方法
CN108595793A (zh) * 2018-04-10 2018-09-28 中南大学 用于铝电解槽过程控制的正常化槽电压数据时频分析方法及设备
CN110532836A (zh) * 2018-05-25 2019-12-03 中广核工程有限公司 核电站信号识别方法、装置、计算机设备及存储介质
CN112613705A (zh) * 2020-12-14 2021-04-06 中广核研究院有限公司 部件质量获取方法、装置、计算机设备和存储介质
CN112949069A (zh) * 2021-03-05 2021-06-11 北京化工大学 一种基于emd的松脱件质量估计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5268948A (en) * 1992-06-25 1993-12-07 Duke Power Company Locking assembly for nuclear fuel bundles
JPH08334442A (ja) * 1995-06-05 1996-12-17 Mitsubishi Heavy Ind Ltd 衝撃検知方法及び装置
JP2001264151A (ja) * 2000-03-23 2001-09-26 Mitsubishi Heavy Ind Ltd 衝撃物の質量検知方法及び装置
CN101718862A (zh) * 2009-11-13 2010-06-02 浙江大学 基于ar模型和小波变换的核电站松动件定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5268948A (en) * 1992-06-25 1993-12-07 Duke Power Company Locking assembly for nuclear fuel bundles
JPH08334442A (ja) * 1995-06-05 1996-12-17 Mitsubishi Heavy Ind Ltd 衝撃検知方法及び装置
JP2001264151A (ja) * 2000-03-23 2001-09-26 Mitsubishi Heavy Ind Ltd 衝撃物の質量検知方法及び装置
CN101718862A (zh) * 2009-11-13 2010-06-02 浙江大学 基于ar模型和小波变换的核电站松动件定位方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
方力先等: "小波变换在松动件检测系统报警中的应用研究", 《原子能科学技术》, vol. 38, no. 5, 30 September 2004 (2004-09-30), pages 432 - 435 *
曹衍龙等: "基于小波包的松动件质量估计方法", 《机械工程学报》, vol. 46, no. 22, 30 November 2010 (2010-11-30), pages 1 - 5 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104408303B (zh) * 2014-11-20 2017-09-12 西安交通大学 一种基于数据匹配的lpms质量估计方法
CN104408303A (zh) * 2014-11-20 2015-03-11 西安交通大学 一种基于数据匹配的lpms质量估计方法
CN104965026A (zh) * 2015-05-25 2015-10-07 北京理工大学 一种使用声发射信号能量特征值预测脆性材料临界失稳的方法
CN105045965A (zh) * 2015-06-26 2015-11-11 西安交通大学 基于仿真数据库的压力容器lpms质量估计方法
CN105045965B (zh) * 2015-06-26 2018-08-28 西安交通大学 基于仿真数据库的压力容器lpms质量估计方法
CN105225433A (zh) * 2015-08-22 2016-01-06 上海核工程研究设计院 一种核电站松动部件监测报警系统及方法
CN105427904A (zh) * 2015-11-13 2016-03-23 陕西卫峰核电子有限公司 一种基于核电站松动部件定位方法
CN106409351A (zh) * 2016-09-13 2017-02-15 中广核工程有限公司 一种核电厂松动部件振动信号的处理方法及装置
CN106874896A (zh) * 2017-03-31 2017-06-20 阳江核电有限公司 一种核电站一回路部件松脱诊断系统信号特征识别的辅助学习方法及系统
CN106874896B (zh) * 2017-03-31 2021-03-02 阳江核电有限公司 一种核电站一回路部件松脱诊断系统信号特征识别的辅助学习方法及系统
CN107194135A (zh) * 2017-07-21 2017-09-22 中国工程物理研究院核物理与化学研究所 基于改进小波算法的反应堆中子噪声频谱计算方法
CN107194135B (zh) * 2017-07-21 2020-06-23 中国工程物理研究院核物理与化学研究所 基于改进小波算法的反应堆中子噪声频谱计算方法
CN107544337B (zh) * 2017-09-19 2020-10-30 中国核动力研究设计院 一种松脱部件监测系统触发数据智能分类的方法
CN107544337A (zh) * 2017-09-19 2018-01-05 中国核动力研究设计院 一种松脱部件监测系统触发数据智能分类的方法
CN108595793A (zh) * 2018-04-10 2018-09-28 中南大学 用于铝电解槽过程控制的正常化槽电压数据时频分析方法及设备
CN108595793B (zh) * 2018-04-10 2022-03-08 中南大学 用于铝电解的正常化槽电压数据时频分析方法及设备
CN110532836A (zh) * 2018-05-25 2019-12-03 中广核工程有限公司 核电站信号识别方法、装置、计算机设备及存储介质
CN110532836B (zh) * 2018-05-25 2022-03-22 中广核工程有限公司 核电站信号识别方法、装置、计算机设备及存储介质
CN112613705A (zh) * 2020-12-14 2021-04-06 中广核研究院有限公司 部件质量获取方法、装置、计算机设备和存储介质
CN112949069A (zh) * 2021-03-05 2021-06-11 北京化工大学 一种基于emd的松脱件质量估计方法
CN112949069B (zh) * 2021-03-05 2023-10-20 北京化工大学 一种基于emd的松脱件质量估计方法

Also Published As

Publication number Publication date
CN103137224B (zh) 2015-07-15

Similar Documents

Publication Publication Date Title
CN103137224A (zh) 基于小波能量谱的核电站松动部件质量估计方法
CN103137225B (zh) 基于小波变换和希尔伯特变换的核电站松动部件定位方法
CN104820993B (zh) 一种联合粒子滤波和跟踪置前检测的水下弱目标跟踪方法
CN107461611B (zh) 基于小波和emd重构相结合的泄漏检测方法及泄漏检测装置
CN103412287B (zh) 基于lvd的线性调频信号参数估计方法
CN109143195B (zh) 一种基于全kl散度的雷达目标检测方法
CN104502899A (zh) 一种自适应的恒虚警率目标检测方法
CN105137437A (zh) 一种基于空域相位方差加权的目标检测方法
CN105044686B (zh) 一种雷达密集假目标干扰抑制方法
CN105702305B (zh) 核电站松动部件多目标三维动态跟踪方法
CN105425223B (zh) 广义帕累托杂波下稀疏距离扩展雷达目标的检测方法
CN105785339A (zh) 非均匀杂波环境下机载雷达杂波协方差矩阵的估计方法
CN109946763A (zh) 一种基于波群识别的畸形波短期预报方法
CN104391279A (zh) 基于电离层传播特性的相径扰动抑制方法
CN104568113B (zh) 一种基于模型的海洋声传播调查爆炸波自动截取方法
CN101425855B (zh) 宽带信道参数抽取方法、宽带信道仿真方法及其装置
CN107390194A (zh) 一种基于全布雷格曼散度的雷达目标检测方法
CN110376575A (zh) 一种基于阻尼参数匹配随机共振的低频线谱检测方法
CN104215939A (zh) 一种融合广义对称结构信息的知识辅助空时自适应处理方法
CN101426213B (zh) 宽带信道仿真方法及其装置
CN103487807B (zh) 一种isar最优成像时间选择方法
CN105093189B (zh) 基于gcv的机载雷达目标检测方法
CN104535977A (zh) 一种基于gsm信号的雷达目标探测方法
CN102830394A (zh) 基于多谱线积累的弱目标探测方法
Wang et al. A novel indoor ranging method using weighted altofrequent RSSI measurements

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 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.

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