CN105702305B - 核电站松动部件多目标三维动态跟踪方法 - Google Patents
核电站松动部件多目标三维动态跟踪方法 Download PDFInfo
- Publication number
- CN105702305B CN105702305B CN201610075672.8A CN201610075672A CN105702305B CN 105702305 B CN105702305 B CN 105702305B CN 201610075672 A CN201610075672 A CN 201610075672A CN 105702305 B CN105702305 B CN 105702305B
- Authority
- CN
- China
- Prior art keywords
- signal
- loose part
- scale
- positioning
- loose
- 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.)
- Expired - Fee Related
Links
- 230000001133 acceleration Effects 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims description 45
- 229910000831 Steel Inorganic materials 0.000 claims description 14
- 238000001228 spectrum Methods 0.000 claims description 14
- 239000010959 steel Substances 0.000 claims description 14
- 238000005516 engineering process Methods 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 abstract description 6
- 238000003745 diagnosis Methods 0.000 abstract description 5
- 230000015556 catabolic process Effects 0.000 abstract 1
- 238000007794 visualization technique Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 23
- 238000012360 testing method Methods 0.000 description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000003909 pattern recognition Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 239000002826 coolant Substances 0.000 description 1
- 239000007799 cork Substances 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 238000009863 impact test Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21C—NUCLEAR REACTORS
- G21C17/00—Monitoring; Testing ; Maintaining
- G21C17/02—Devices or arrangements for monitoring coolant or moderator
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Plasma & Fusion (AREA)
- General Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种核电站松动部件多目标三维动态跟踪方法。本发明包括以下步骤：步骤(1)数据库的建立；步骤(2)三维虚拟建模；步骤(3)松动部件现场信号获取；步骤(4)松动部件定位；步骤(5)松动部件质量估计；步骤(6)动态跟踪显示。本发明通过加速度传感器周期地获取到松动部件的信号，对松动部件进行定位和质量估计后，利用三维可视化技术，实现动态跟踪松动部件，为专家诊断和故障维修提供依据。本发明抗噪能力强，定位精度高，质量估计误差小，并能同时跟踪多个松动部件目标。
Description
技术领域
本发明属于核工程技术领域，具体涉及松动部件定位、松动部件质量估计和松动部件三维动态跟踪。该方法不仅能动态跟踪松动部件状态，而且能跟踪多个不同质量的松动部件。
技术背景
为了监测核电站轻水反应堆冷却剂系统中可能存在脱落或者松动的部件，在核电站一回路安装松动部件监测系统(Loose Parts Monitoring System，LPMS)是提高核电站运行安全性，稳定性行之有效的方式。但是，传统LPMS检测到脱落的松动部件后，核电站并非能够马上停堆检修。因此需要动态跟踪松动部件，形象的显示核电站内部工作过程和松动部件跌落轨迹，有利于在停堆检修时快速找到松动部件，使维修人员的辐射量减到最小，同时为专家诊断提供依据。
三维动态显示是二维平面的扩展，目前该技术应用于数字电视、远程教育、远程工业控制、三维视频会议系统、虚拟现实等领域。在LPMS中的报警、松动部件定位以及松动部件质量估计研究依然停留在二维空间，虽然可以比较方便的检测出核电站可能存在的问题，但是缺少直观方便性，为专家诊断和维修带来困难。
现有的松动部件定位相关文献有：
[1]Ziola S M,Source location in thin plates using cross-correlation[J],Journal of the Acoustical Society of America,l991,90:2551-2556.
一种基于信号相关的到达时差定位法。定位效果虽然准确，但抗噪能力较差。
[2]Figedy S,Oksa G.Modern methods of signal processing in the loosepart monitoring system[J]，Progress in Nuclear Energy,2005,46(3-4):253-267.
一种基于小波去噪的时差估计方法，通过小波去噪，即去除噪声的影响，然后再估计时差。但该方法还是以信号的过零点作为信号到达时间，没有考虑弯曲波的复杂传播模式，所以在实际定位中当传播距离较远时仍然存在较大的误差。
[3]林京,罗志,曾亮,黄利平,高飞,华佳东.基于希尔伯特变换和数据筛选的核电站松动件定位方法:陕西,CN104376881A[P].2015-02-25.
一种基于希尔伯特变换和数据筛选的核电站松动件定位方法，以Butterworth滤波器松动部件冲击信号的噪声后，进行希尔伯特变换得到冲击信号起振时刻，并对网格中心点的名义传播速度进行数据筛选，计算方差，方差中最小值的位置即为定位结果。该方法定位精度高，搜寻速度快，但对松动部件没有动态跟踪定位，缺少直观方便性，为专家诊断和维修带来困难。
现有的松动部件质量估计方法相关文献有：
[1]MAYO W.loose part mass and energy estimation[J].Progress inNuclear Energy,1999,34(3):263-282.
一种基于碰撞信号高低频段能量比的松动件质量估计方法。该方法质量估计误差较大，一致性差等问题，尤其是抗干扰能力差。
[2]模式识别技术用于核电站松动件质量估计模型研究[J]，核动力工程，2001,22(5):465-470.
一种利用模式识别理论和人工神经网络理论实现对核电站松动件质量进行估计的数学模型。利用模式识别理论的特征提取技术，实现对松动件碰撞信号的特征压缩，从而形成比较少的特征空间维数，经过神经网络的学习过程，实现对核电站松动件质量大小的估计。该方法能够得到较高的质量估计精度，但需要大量的试验数据对其进行训练，且一致性较差。
[3]曹衍龙，程实，杨将新，郑华文，何元峰，基于小波包的质量估计方法[J]，机械工程学报，2010.46(22):1-5
一种基于小波包的松动件质量估计方法。该方法通过小波包分解获取各频段小波包系数，计算各频段小波包系数的能量潜，并进行归化，得到特征值矢量；以此特征值矢量作为输入参数，以松动件质量作为输出参数，应用支持向量机((Support vector machine,SVM)实现松动件质量估计。
[4]方力先,季田田,曾复,张卫,万鹏飞,谢永诚,王赤虎,张海丰,邓晶晶.基于小波能量谱的核电站松动部件质量估计方法:浙江,CN103137224A[P].2013-06-05
一种基于小波能量谱的核电站松动部件质量估计方法。该方法通过建立已知质量的松动部件的建立尺度函数，来求解松动部件的质量。该方法误差小、一致性好和抗干扰能力强的优点。
发明内容
本发明为了实现多目标三维动态跟踪松动部件，提出一种核电站松动部件多目标三维动态跟踪方法。本发明采用的技术方案包括以下步骤：
步骤(1)数据库的建立。
1-1由于松动部件主要发生在压力容器顶部或底部，在采集数据之前，根据反应堆压力容器顶部和底部的曲面划分多个均等的网格小分块(常用网格面积约50*50cm2)，并对网格分块进行编号，然后在压力容器顶部和底部各安装三个呈等边三角形分布的加速度传感器。
1-2通过数据采集卡、电荷放大器和加速度传感器获取力锤敲击压力容器某一固定点的冲击信号C(ti)，用于测定传播速度。具体的：加速度传感器获取力锤敲击信号后，经电荷放大器进入数据采集卡。对冲击信号C(ti)进行基于小波变换的希尔伯特变换，得到冲击信号C(ti)到达某个加速度传感器的传播时间t，测量固定点到该加速度传感器的弧长分别为L，弧长除以传播时间，得到传播速度V。
1-3.测量各网格分块中心点到三个加速度传感器的弧长Si(i＝1,2,3)，然后将弧长Si除以传播速度V，得到敲击点到每个加速度传感器的时间Ti(i＝1,2,3)，然后获取任意一个加速度传感器与另两个加速度传感器时间差然后以网格编号，时间差和为查询表的属性，建立定位数据库。
Ti＝Si/V；
步骤(2)三维虚拟建模。根据核反应堆一回路系统中设备的实际图纸，利用三维建模软件建立核反应堆一回路系统中各个子模块的3D模型，并利用纹理贴图技术、特征简化技术和模型轻量化技术尽可能地减小文件、提高系统的运算速度。该3D模型可以透视，旋转，从不同角度观察内部部件。
步骤(3)松动部件现场信号获取。通过数据采集卡、电荷放大器和加速度传感器来实现现场信号x(ti)的采集，现场信号x(ti)包括松动件跌落时的冲击信号和环境背景噪声。
步骤(4)松动部件定位。
4-1对现场信号x(ti)进行频域分析确定信号主频率fm，用连续小波变换计算得到尺度因子a。选数值2、8、14、20中最接近尺度因子a，作为最终连续小波变换尺度对现场信号x(ti)进行处理，然后通过Hilbert包络线法求得信号的到达时间差t21、t31。
4-2将t21、t31与数据库中对应尺度下的定位查询表中的时间差作比较，时间差距离最小值所对应的网格即为定位网格。时间差距离的计算公式为：
式中：i为网格编号；n为网格划分总数；dn×1的最小值所对应的i即为定位网格。
步骤(5)松动部件质量估计
5-1对现场信号x(ti)出现最大幅值出现前0.03s开始，取0.3s长信号数据g(ti)，对信号g(ti)进行连续小波变换处理得到信号m(ti)后，计算信号m(ti)的小波能量谱，根据小波能量谱图找出该图尺度峰值xa。
5-2将尺度峰值xa带入已经建立好的尺度峰值函数F(xa)，求得松动部件的质量估计值y，然后对多个通道中的松动部件质量估计值求平均，得到最终的松动部件质量。
步骤(6)动态跟踪显示。LPMS以时间T(该值可调，一般采用1h)为周期，周期地采集核电站一回路中的松动部件的现场信号。在Ti(i＝0,1,2,3……)时刻采集松动部件的现场信号，根据步骤(4)定位该松动部件的位置，根据步骤(5)估计该松动部件的质量，同时在压力容器曲面的定位位置处显示该松动部件。
在Ti+1(Ti+1＝Ti+T，i＝0,1,2,3……)时刻，再次采集松动部件的现场信号，根据步骤(4)定位该松动部件的位置，根据步骤(5)估计该松动部件的质量，若该质量与先前任意周期内松动部件质量估计值不同，则为新增定位目标，同时在压力容器曲面的定位位置处显示此新增松动部件；若该质量与先前任意周期内松动部件质量估计值相同，则在压力容器曲面的定位位置处显示该松动部件，并标记该定位目标的移动轨迹。如此周期循环，动态的构成松动部件跌落轨迹。
所述的步骤(4)连续小波变换的公式为：
其中，f(ti)为步骤(3)中的x(ti)，ψa,b(ti)为小波函数，此处取小波函数为db2；为ψa,b(ti)的复共轭；a为尺度因子；b为时间平移因子。WTf(a,b)即为经尺度因子a连续变换后的信号。
数字角频率与模拟角频率的关系：
式中：ω为数字角频率，以2π为周期；Ω为模拟角频率；f为模拟频率；fs为信号的采样频率。设ψ(t)的频域中心为ω0，则ψa,b(t)的中心频率(数字角频率)，根据数字角频率与模拟角频率的关系式，可以得到数字中心频率ωa,b＝Ωa,b·Ts，小波变换时尺度和频率之间的关系：
式中：Ts为采样周期；Fc为给定小波的中心频率，为归一化后的数值；Fa为伪频率。
所述的步骤(4)通过Hilbert包络线求时间差：
设给定信号a(t)，其希尔伯特变换的定义：
以信号a(t)为实部，其希尔伯特变换为虚部，构成一个新的函数：
式中：是新函数的幅值函数；为相位函数。
则|z(t)|即为信号a(t)的包络线函数，反映了信号能量和时间的关系。该方法通过提取信号包络线，认为包络线的首个峰值所对应的时间点即为振动信号到达测量点的时刻。
所述的步骤(5)中小波能量谱的公式为：
式中，z(ti)为步骤(5)中的m(ti)；E(a)称为在尺度a上的小波能量谱。则被分析的时域信号F(t)的能量分布为：
式中，EF反映了信号经小波变换后能量随尺度的变化情况。
步骤(5)中的尺度峰值函数F(xa)的建立是通过线性插值法建立不同质量钢球与尺度峰值均值之间的尺度峰值函数F(xa)。线性插值法：设质量数组L＝[l1,l2,l3,...ln],与质量数组为L对应的小波尺度均值数组P＝[p1,p2,p3,...pn]，变量xa为输入小波尺度，变量y为质量估计值。则可以建立尺度峰值函数F(xa)：
本发明有益效果如下：
本发明通过加速度传感器周期地获取到松动部件的信号，对松动部件进行定位和质量估计后，利用三维可视化技术，实现动态跟踪松动部件，为专家诊断和故障维修提供依据。本发明抗噪能力强，定位精度高，质量估计误差小，并能同时跟踪多个松动部件目标。
附图说明
图1本发明流程图；
图2本发明加速度传感器布局图；
图3本发明冲击信号希尔伯特包络线图；
图4本发明同一质量冲击信号三通道小波能量谱图；
图5本发明不同质量冲击信号小波能量谱图；
图6本发明动态跟踪定位结果显示图；
图7本发明中冲击质量与尺度峰值均值关系。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
如图1所示，核电站松动部件多目标动态跟踪方法，包括以下步骤：
步骤(1)数据库的建立。
1-1由于松动部件主要发生在压力容器顶部或底部，在采集数据之前，根据反应堆压力容器顶部和底部的曲面划分多个均等的网格小分块(常用网格面积约50*50cm2)，并对网格分块进行编号，然后在压力容器顶部和底部各安装三个呈等边三角形分布的加速度传感器，传感器的安装方式如图2所示。
1-2通过数据采集卡、电荷放大器和加速度传感器获取力锤敲击压力容器某一固定点的冲击信号C(ti)，用于测定传播速度。具体的：加速度传感器获取力锤敲击信号后，经电荷放大器进入数据采集卡。对冲击信号C(ti)进行基于小波变换的希尔伯特变换，得到冲击信号C(ti)到达某个加速度传感器的传播时间t，测量固定点到该加速度传感器的弧长分别为L，弧长除以传播时间，得到传播速度V。
1-3.测量各网格分块中心点到三个加速度传感器的弧长Si(i＝1,2,3)，然后将弧长Si除以传播速度V，得到敲击点到每个加速度传感器的时间Ti(i＝1,2,3)，然后获取任意一个加速度传感器与另两个加速度传感器时间差然后以网格编号，时间差和为查询表的属性，建立定位数据库。
Ti＝Si/V；
步骤(2)三维虚拟建模。根据核反应堆一回路系统中设备的实际图纸，利用三维建模软件建立核反应堆一回路系统中各个子模块的3D模型，并利用纹理贴图技术、特征简化技术和模型轻量化技术尽可能地减小文件、提高系统的运算速度。该3D模型可以透视，旋转，从不同角度观察内部部件。
步骤(3)松动部件现场信号获取。通过数据采集卡、电荷放大器和加速度传感器来实现现场信号x(ti)采集，现场信号x(ti)包括松动件跌落时的冲击信号和环境背景。
步骤(4)松动部件定位。
4-1对现场信号x(ti)进行频域分析确定信号主频率fm，用连续小波变换计算得到尺度因子a。选数值2、8、14、20中最接近尺度因子a，作为最终连续小波变换尺度对现场信号x(ti)进行处理。然后通过Hilbert包络线法求得信号的到达时间差t21、t31。
4-2将t21、t31与数据库中对应尺度下的定位查询表中的时间差作比较，时间差距离最小值所对应的网格即为定位网格。时间差距离的计算公式为：
式中：i为网格编号；n为网格划分总数；dn×1的最小值所对应的i即为定位网格。
步骤(5)松动部件质量估计
5-1对得到的现场信号x(ti)出现最大幅值出现前0.03s开始，取0.3s长信号数据g(ti)，对信号g(ti)进行连续小波变换处理得到信号m(ti)后，计算信号m(ti)的小波能量谱，根据小波能量谱图找出该图尺度峰值xa。
5-2将尺度峰值xa带入已经建立好的尺度峰值函数F(xa)，求得松动部件的质量估计值y，然后对多个通道中的松动部件质量估计值求平均，得到最终的松动部件质量。
步骤(6)动态跟踪显示。LPMS以时间T(该值可调，一般采用1h)为周期，周期地采集核电站一回路中的松动部件的冲击信号。在Ti(i＝0,1,2,3……)时刻采集松动部件的冲击信号，根据步骤(4)定位该松动部件的位置，根据步骤(5)估计该松动部件的质量，同时在压力容器曲面的定位位置处显示该松动部件。
在Ti+1(Ti+1＝Ti+T，i＝0,1,2,3……)时刻，再次采集松动部件的冲击信号，根据步骤(4)定位该松动部件的位置，根据步骤(5)估计该松动部件的质量，若该质量与先前任意周期内松动部件质量估计值不同，则为新增定位目标，同时在压力容器曲面的定位位置处显示此新增松动部件；若该质量与先前任意周期内松动部件质量估计值相同，则在压力容器曲面的定位位置处显示该松动部件，并标记该定位目标的移动轨迹。如此周期循环，动态的构成松动部件跌落轨迹。
如图6所示，时间周期设置为1h，初始松动部件定位结果为27，同样在模型中标记该位置；在下一周期，采集的松动部件冲击信号定位结果为7号网格，同样在模型中标记该位置。标记的位置构成松动部件运动轨迹。
所述的步骤(4)连续小波变换的公式为：
其中，f(ti)为步骤(3)中的x(ti)，；ψa,b(ti)为小波函数，此处取小波函数为db2；为ψa,b(ti)的复共轭；a为尺度因子；b为时间平移因子。WTf(a,b)即为经尺度因子a连续变换后的信号。
数字角频率与模拟角频率的关系：
式中：ω为数字角频率，以2π为周期；Ω为模拟角频率；f为模拟频率；fs为信号的采样频率。设ψ(t)的频域中心为ω0，则ψa,b(t)的中心频率(数字角频率)，根据数字角频率与模拟角频率的关系式，可以得到数字中心频率ωa,b＝Ωa,b·Ts，小波变换时尺度和频率之间的关系：
式中：Ts为采样周期；Fc为给定小波的中心频率，为归一化后的数值；Fa为伪频率。
所述的步骤(4)通过Hilbert包络线求时间差：
设给定信号a(t)，其希尔伯特变换的定义：
以信号a(t)为实部，其希尔伯特变换为虚部，构成一个新的函数：
式中：是新函数的幅值函数；为相位函数。
则|z(t)|即为信号a(t)的包络线函数，反映了信号能量和时间的关系。该方法通过提取信号包络线，如图3所示，从图分析得出，认为包络线的首个峰值所对应的时间点即为振动信号到达测量点的时刻。
所述的步骤(5)中小波能量谱的公式为：
式中，z(ti)为步骤(5)中的m(ti)；E(a)称为在尺度a上的小波能量谱。则被分析的时域信号F(t)的能量分布为：
式中，EF反映了信号经小波变换后能量随尺度的变化情况。
步骤(5)中的尺度峰值函数F(xa)的建立是通过线性插值法建立不同质量钢球与尺度峰值均值之间的尺度峰值函数F(xa)。线性插值法：设质量数组L＝[l1,l2,l3,...ln]，与质量数组为L对应的小波尺度均值数组P＝[p1,p2,p3,...pn]，其中pi(1≤i≤n)为对应于质量li的尺度峰值均值，变量xa为输入小波尺度，变量y为质量估计值。则可以建立尺度峰值函数F(xa)：
如图4所示，从同一质量钢球冲击信号三通道的小波能量谱图中不难看出，三通道传感器接收到的能量分布具有一致性，且能量峰值具有较一致的尺度，这与Hertz理论指出的碰撞主频率是相吻合的。
如图5所示，从不同质量钢球的冲击信号归一化后的小波能量谱图中可看出，随着冲击质量的增大，小波能量谱向右移动。这些明显的特征变化为质量估计提供了重要依据，即可以根据尺度峰值随着冲击物质量的变化而发生转移的情况来估计松动部件的质量。
实施例
下面结合实施例对本发明作进一步说明：
本发明的效果可以通过曲面实验的实验结果分析说明：
1.试验条件
本试验的试验平台由测试对象、传感器、电荷放大器、数据采集卡和计算机构成。测试对象包括：钢球、锅炉封顶及其支撑。锅炉封顶尺寸为R＝50cm，每个网格面积约为20*20cm2。为了尽量减少环境噪声的影响，在锅炉封顶的下面加了缓冲隔离，缓冲隔离由橡胶垫以及软木构成。钢球重量分别为44g、175g、225g、360g、510g、640g、880g、1120g。
2.数据采集
数据采集包括两部分：冲击信号采集和背景噪声采集。
(1)冲击信号采集：采用三个加速度传感器，加速度传感器在锅炉上呈等边三角形分布排列，可参看附图2。冲击物为不同质量钢球。冲击信号为不同质量钢球，分别在高度10cm、15cm、20cm处敲击钢板不同点的冲击信号。
(2)背景噪声采集：采用三个加速度传感器，其中一个加速度传感器位于锅炉顶部进水口处，一个位于锅炉出水口，一个位于锅炉顶部左上方。背景噪声为锅炉从冷态到热态过程中的背景噪声，根据锅炉运行状态(锅炉未运行前、小火、大火、水泵运转、放气等)来采集数据，锅炉运行状态每改变一次采集一组数据。
3.冲击试验结果及分析
冲击信号与背景噪声分别按信噪比为0dB、5bB进行叠加，用于试验结果分析。信噪比定义为：
其中，表示噪声幅值取绝对值后求平均，表示信号幅值取绝对值后求最大值。
4.仿真结果
4.1松动部件定位的仿真结果
误差计算公式为Er＝||ve||,ve＝(x,y)；其中分别跌落位置网格与查表法定位网格之间的水平网格数和垂直网格数，以网格为单位。对不同质量钢球冲击信号进行定位分析，得到结果如表1、表2所示。加入0db和5db背景噪声信号定位分析如表3所示。
表1 定位误差表1
表2 定位误差表2
注：图中误差表示为误差偏离方向+偏离距离，例“右上1”表示定位网格相较于实际网格偏离了1格的距离，偏离方向为右上方。
表3 背景噪声分析表
由定位误差表分析可得，力锤和不同质量钢球的定位误差都控制在1个网格范围左右。由于连续小波变换使频率成分更集中，定位精度有很大提高。不同质量加入相同背景噪声，定位效果稳定，小质量钢球的冲击信号加入0db信噪比背景噪声，定位仍然较准确，同一质量钢球冲击信号加入5db信噪比噪声的情况下仍能较准确定位，抗噪能力强。
4.2松动部件质量估计仿真结果
用如下方法来计算绝对误差和相对误差：
绝对误差：E＝Δm，其中Δm为松动件质量真实值与估计值之间的差值。
相对误差：其中m为松动件质量真实值。
试验结果对本发明提出的松动部件质量估计方法进行了验证，通过试验数据分析，建立尺度峰值函数，通过尺度峰值函数来实现松松部件的质量估计。附图7为松动件质量与尺度之间的对应关系，表4为各冲击质量对应的尺度均值，从附图7和表4中不难看出，小质量的尺度峰值变化相对明显，但随着冲击质量的增大，峰值变化逐渐变得不明显，这是因为冲击质量大的信号，其能量主要集中在低频范围内，在频谱上很难区分。
表4 不同尺度峰值所对应的质量估计值
| 尺度峰值 | 质量估值/g | 尺度峰值 | 质量估值/g | 尺度峰值 | 质量估值/g |
| 4 | 20.2739 | 11 | 514.6258 | 18 | 1938.8 |
| 4.5 | 32.2727 | 11.5 | 637.2932 | 18.5 | 2078 |
| 5 | 44.4041 | 12 | 714.2722 | 19 | 2217.1 |
| 5.5 | 62.2652 | 12.5 | 790.2203 | 19.5 | 2356.3 |
| 6 | 80.1262 | 13 | 866.1683 | 20 | 2496.5 |
| 6.5 | 97.9872 | 13.5 | 948.214 | 20.5 | 2829.1 |
| 7 | 119.8749 | 14 | 1031.6 | 21 | 3204.4 |
| 7.5 | 142.2741 | 14.5 | 1115 | 21.5 | 3579.7 |
| 8 | 164.6732 | 15 | 1208.7 | 22 | 3955 |
| 8.5 | 227.6932 | 15.5 | 1303 | 22.5 | 5264.4 |
| 9 | 270.5928 | 16 | 1397.4 | 23 | 6701.1 |
| 9.5 | 313.4926 | 16.5 | 1530.7 | 23.5 | 8137.9 |
| 10 | 356.3921 | 17 | 1665.2 | 24 | 9676.9 |
| 10.5 | 433.1475 | 17.5 | 1799.7 | 24.5 | 1105.3 |
表5 质量估计误差
| 质量/g | 估计误差 | 质量/g | 估计误差 | 质量/g | 估计误差 | 质量/g | 估计误差 |
| 20 | 15％ | 225 | 8％ | 880 | 11％ | 2520 | 19.4％ |
| 44 | 9.1％ | 360 | 10％ | 1120 | 13％ | 40000 | 15％ |
| 100 | 11％ | 510 | 15％ | 1400 | 13.5％ | 9000 | 17％ |
| 175 | 9％ | 640 | 12.5％ | 1800 | 10.2％ | 11200 | 18.4％ |
表5为质量估计的误差，从表中可以看出，质量估计误差最大为19.4％，最小为9.1％，质量估计误差基本在15％左右。所以本发明所提出的方法具有较高的准确率。且抗干扰能力较强，在背景噪声频谱结构已知的情况下，在信噪比为0dB时仍能较准确的进行松动件质量估计。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围不应当被视为仅限于实施例所所陈述的具体形式，本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。
Claims (4)
1.核电站松动部件多目标三维动态跟踪方法，其特征在于包括以下步骤：
步骤(1)数据库的建立；
1-1.由于松动部件主要发生在压力容器顶部或底部，在采集数据之前，根据反应堆压力容器顶部和底部的曲面划分多个均等的网格小分块，并对网格分块进行编号，然后在压力容器顶部和底部各安装三个呈等边三角形分布的加速度传感器；
1-2.通过数据采集卡、电荷放大器和加速度传感器获取力锤敲击压力容器某一固定点的冲击信号C(ti)，用于测定传播速度；具体的：加速度传感器获取力锤敲击信号后，经电荷放大器进入数据采集卡；对冲击信号C(ti)进行基于小波变换的希尔伯特变换，得到冲击信号C(ti)到达某个加速度传感器的传播时间t，测量固定点到该加速度传感器的弧长为L，弧长除以传播时间，得到传播速度V；
1-3.测量各网格分块中心点到三个加速度传感器的弧长Si，然后将弧长Si除以传播速度V，得到敲击点到每个加速度传感器的时间Ti，然后获取任意一个加速度传感器与另两个加速度传感器时间差然后以网格编号，时间差和为查询表的属性，建立定位数据库；其中i＝1,2,3；
Ti＝Si/V；
步骤(2)三维虚拟建模；根据核反应堆一回路系统中设备的实际图纸，利用三维建模软件建立核反应堆一回路系统中各个子模块的3D模型，并利用纹理贴图技术、特征简化技术和模型轻量化技术减小文件、提高系统的运算速度；该3D模型能够透视和旋转，从不同角度观察内部部件；
步骤(3)松动部件现场信号获取；通过数据采集卡、电荷放大器和加速度传感器来实现现场信号x(ti)的采集，现场信号x(ti)包括松动件跌落时的冲击信号和环境背景噪声；
步骤(4)松动部件定位；
4-1.对现场信号x(ti)进行频域分析确定信号主频率fm，用连续小波变换计算得到尺度因子a；选数值2、8、14、20中最接近尺度因子a的数值作为最终连续小波变换尺度对现场信号x(ti)进行处理，然后通过Hilbert包络线法求得信号的到达时间差t21、t31；
4-2.将t21、t31与数据库中对应尺度下的定位查询表中的时间差作比较，时间差距离最小值所对应的网格即为定位网格；时间差距离的计算公式为：
式中：i为网格编号；n为网格划分总数；dn×1的最小值所对应的i即为定位网格；
步骤(5)松动部件质量估计
5-1.对现场信号x(ti)最大幅值出现前0.03s开始，取0.3s长信号数据g(ti)，对信号数据g(ti)进行连续小波变换处理得到信号m(ti)后，计算信号m(ti)的小波能量谱，根据小波能量谱图找出该图尺度峰值xa；
5-2.将尺度峰值xa带入已经建立好的尺度峰值函数F(xa)，求得松动部件的质量估计值y，然后对多个通道中的松动部件质量估计值求平均，得到最终的松动部件质量；
步骤(6)动态跟踪显示；松动部件监测系统以时间T为周期，周期地采集核电站一回路中的松动部件的现场信号；在Ti时刻采集松动部件的现场信号，其中i＝0,1,2,3……；根据步骤(4)定位该松动部件的位置，根据步骤(5)估计该松动部件的质量，同时在压力容器曲面的定位位置处显示该松动部件；
在Ti+1时刻，再次采集松动部件的现场信号，其中Ti+1＝Ti+T，i＝0,1,2,3……；根据步骤(4)定位该松动部件的位置，根据步骤(5)估计该松动部件的质量，若该质量与先前任意周期内松动部件质量估计值不同，则为新增定位目标，同时在压力容器曲面的定位位置处显示此新增松动部件；若该质量与先前任意周期内松动部件质量估计值相同，则在压力容器曲面的定位位置处显示该松动部件，并标记该定位目标的移动轨迹；如此周期循环，动态的构成松动部件跌落轨迹。
2.根据权利要求1所述的核电站松动部件多目标三维动态跟踪方法，其特征在于所述的步骤(4)连续小波变换的公式为：
其中，f(t)为步骤(3)中的x(ti)，ψa,b(t)为小波函数，此处取小波函数为db2；为ψa,b(t)的复共轭；a为尺度因子；b为时间平移因子；WTf(a,b)即为经尺度因子a连续变换后的信号；
数字角频率与模拟角频率的关系：
式中：ω为数字角频率，以2π为周期；Ω为模拟角频率；f为模拟频率；fs为信号的采样频率；设ψ(t)的频域中心为ω0，则ψa,b(t)的中心频率，即数字角频率根据数字角频率与模拟角频率的关系式，得到数字中心频率ωa,b＝Ωa,b·Ts，小波变换时尺度和频率之间的关系：
式中：Ts为采样周期；Fc为给定小波的中心频率，为归一化后的数值；Fa为伪频率。
3.根据权利要求1所述的核电站松动部件多目标三维动态跟踪方法，其特征在于所述的步骤(4)通过Hilbert包络线求时间差：
设给定信号a(t)，其希尔伯特变换的定义：
以信号a(t)为实部，其希尔伯特变换为虚部，构成一个新的函数：
式中：是新函数的幅值函数；为相位函数；
则|z(t)|即为信号a(t)的包络线函数，反映了信号能量和时间的关系；该方法通过提取信号包络线，认为包络线的首个峰值所对应的时间点即为振动信号到达测量点的时刻。
4.根据权利要求2所述的核电站松动部件多目标三维动态跟踪方法，其特征在于所述的步骤(5)中小波能量谱的公式为：
式中，z(ti)为步骤(5)中的m(ti)；E(a)称为在尺度因子a上的小波能量谱；则经过连续小波变换处理后的时域信号F(t)的能量分布为：
式中，EF反映了信号经小波变换后能量随尺度的变化情况；
步骤(5)中的尺度峰值函数F(xa)的建立是通过线性插值法建立不同质量钢球与尺度峰值均值之间的尺度峰值函数F(xa)；线性插值法：设质量数组L＝[l1,l2,l3,...ln],与质量数组为L对应的小波尺度均值数组P＝[p1,p2,p3,...pn]，变量xa为输入小波尺度，变量y为质量估计值；则可以建立尺度峰值函数F(xa)：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610075672.8A CN105702305B (zh) | 2016-02-03 | 2016-02-03 | 核电站松动部件多目标三维动态跟踪方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610075672.8A CN105702305B (zh) | 2016-02-03 | 2016-02-03 | 核电站松动部件多目标三维动态跟踪方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105702305A CN105702305A (zh) | 2016-06-22 |
| CN105702305B true CN105702305B (zh) | 2019-01-15 |
Family
ID=56230236
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610075672.8A Expired - Fee Related CN105702305B (zh) | 2016-02-03 | 2016-02-03 | 核电站松动部件多目标三维动态跟踪方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105702305B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107194135B (zh) * | 2017-07-21 | 2020-06-23 | 中国工程物理研究院核物理与化学研究所 | 基于改进小波算法的反应堆中子噪声频谱计算方法 |
| CN107748049B (zh) * | 2017-09-04 | 2019-09-20 | 西安交通大学 | 基于椭球面近似最短路径的核电站松动件定位方法 |
| CN107633127A (zh) * | 2017-09-15 | 2018-01-26 | 西北工业大学 | 一种基于多尺度动力学分离的连接结构松动诊断方法 |
| CN107544337B (zh) * | 2017-09-19 | 2020-10-30 | 中国核动力研究设计院 | 一种松脱部件监测系统触发数据智能分类的方法 |
| CN109003686A (zh) * | 2018-05-25 | 2018-12-14 | 中广核工程有限公司 | 核电站松脱部件的定位方法、装置、计算机设备及存储介质 |
| CN110532836B (zh) * | 2018-05-25 | 2022-03-22 | 中广核工程有限公司 | 核电站信号识别方法、装置、计算机设备及存储介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4641526A (en) * | 1983-06-06 | 1987-02-10 | Hitachi, Ltd. | Method and apparatus for estimating sound source position |
| CN103137225A (zh) * | 2013-01-25 | 2013-06-05 | 杭州电子科技大学 | 基于小波变换和希尔伯特变换的核电站松动部件定位方法 |
| CN104376881A (zh) * | 2014-11-20 | 2015-02-25 | 西安交通大学 | 基于希尔伯特变换和数据筛选的核电站松动件定位方法 |
| JP6388410B2 (ja) * | 2016-08-31 | 2018-09-12 | キヤノン株式会社 | 通信装置およびその制御方法、並びにプログラム |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6388410A (ja) * | 1986-10-01 | 1988-04-19 | Babcock Hitachi Kk | 原子炉流体路で発生したル−スパ−ツの監視装置 |
-
2016
- 2016-02-03 CN CN201610075672.8A patent/CN105702305B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4641526A (en) * | 1983-06-06 | 1987-02-10 | Hitachi, Ltd. | Method and apparatus for estimating sound source position |
| CN103137225A (zh) * | 2013-01-25 | 2013-06-05 | 杭州电子科技大学 | 基于小波变换和希尔伯特变换的核电站松动部件定位方法 |
| CN104376881A (zh) * | 2014-11-20 | 2015-02-25 | 西安交通大学 | 基于希尔伯特变换和数据筛选的核电站松动件定位方法 |
| JP6388410B2 (ja) * | 2016-08-31 | 2018-09-12 | キヤノン株式会社 | 通信装置およびその制御方法、並びにプログラム |
Non-Patent Citations (1)
| Title |
|---|
| 小波变换在核电站一回路松动件定位的应用研究;方力先 等;《核科学与工程》;20010331;第21卷(第1期);50-53页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105702305A (zh) | 2016-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105702305B (zh) | 核电站松动部件多目标三维动态跟踪方法 | |
| CN103137224B (zh) | 基于小波能量谱的核电站松动部件质量估计方法 | |
| CN108444653B (zh) | 一种面阵型超声波气体泄漏点定位检测装置及方法 | |
| CN103137225B (zh) | 基于小波变换和希尔伯特变换的核电站松动部件定位方法 | |
| CN109557538A (zh) | 基于海用相参雷达测量海浪参数的方法 | |
| CN104376881B (zh) | 基于希尔伯特变换和数据筛选的核电站松动件定位方法 | |
| CN105589066B (zh) | 一种利用垂直矢量阵估计水下匀速运动航行器参数的方法 | |
| CN104408303B (zh) | 一种基于数据匹配的lpms质量估计方法 | |
| CN108646221A (zh) | 一种基于声源阵列的空间麦克风定位方法 | |
| JPH1151300A (ja) | リーク監視装置 | |
| CN103430030A (zh) | 用于通过远程测量径向速度确定流体移动的方法和装置 | |
| CN104507053A (zh) | 差分气压测高辅助 wlan 指纹定位中的楼层判别方法 | |
| Li et al. | A novel ultrasonic array signal processing scheme for wind measurement | |
| CN114187330A (zh) | 一种基于光流法的结构微幅振动工作模态分析方法 | |
| CN112433208A (zh) | 一种双站圆周探地雷达快速时域成像方法及系统 | |
| Comesana et al. | An introduction to virtual phased arrays for beamforming applications | |
| CN107884743A (zh) | 适用于任意结构声音阵列的波达方向智能估计方法 | |
| CN106225916B (zh) | 固定式海洋平台安全可靠性的定量、在线的检测方法 | |
| CN105045965B (zh) | 基于仿真数据库的压力容器lpms质量估计方法 | |
| CN102073037A (zh) | 基于自适应阈值选取技术的迭代海流反演方法 | |
| CN109085248B (zh) | 承压管道壁冲击源的定位方法、装置及系统 | |
| CN111965628B (zh) | 一种垂直出水航行体瞬时波浪参数估计方法 | |
| Li et al. | A new iterative near-field coherent subspace method for rub-impact fault localization using AE technique | |
| Zhe et al. | A Continuous Short RMS-Based Grid Time Difference Localization Method for Reactor Loose Parts | |
| CN113639934A (zh) | 气体泄漏位置三维定位方法、定位系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190115 Termination date: 20200203 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |