一种基于EMD的松脱件质量估计方法

技术领域

本发明涉及核工程技术领域，特别涉及一种基于EMD的松脱件质量估计方法。

背景技术

松脱件监测系统(Loose Parts Monitoring System，LPMS)是保证核反应堆安全、稳定运行的系统之一，其中松脱件质量估计是LPMS中重要的一环，准确的质量估计能对松脱件类型及危害程度进行判断，为后续反应堆检修提供合理建议。现有松脱件质量估计方法如下：

(1)基于碰撞信号频率比的质量估计方法。该方法将信号频谱高低频段面积比作为特征识别不同质量冲击信号，频率分布易受传播距离、跌落高度等多种因素影响，此特征不能很好地描述信号随质量变化规律，导致质量估计误差大。

(2)基于有限元分析的质量估计方法。该方法通过有限元模拟冲击信号传播过程得到用于质量估计的信号图谱，根据信号频率及幅度即可从信号图谱中判断松脱件质量及对应碰撞速度。但由于此方法质量估计依赖于有限元模型，难以实现实时在线的松脱件质量估计。

(3)基于小波包的质量估计方法。该方法通过对信号进行小波包处理，归一化各频段小波包系数能量谱得到特征值矢量，将特征值矢量输入支持向量机进行质量估计。主要存在问题为：算法复杂，计算时间较长，实际应用价值较差。

发明内容

本发明为解决上述问题提出了一种基于EMD的松脱件质量估计方法。本方法通过EMD对信号进行预处理，提取信号主要IMF成分imf

基于EMD的松脱件质量估计方法包含以下3个步骤：(1)冲击信号采集及预处理；(2)基于信号主要IMF成分的质量估计特征提取；(3)质量估计函数模型建立及质量估计，详述如下：

步骤一：冲击信号采集与EMD处理

首先通过在反应堆压力容器底部布置3个加速度传感器实时采集不同质量松脱件冲击信号用于后续质量估计函数模型建立；记所采集冲击信号为A(t)；

基于Hertz理论，以小球冲击平板为例，小球质量M可表述为：

M(V

式中，f为信号主频率，k

由式可知，为准确估算小球质量，关键在于碰撞速度V

式中，imf

步骤二：主要IMF成分确定与特征提取

信号A(t)经EMD处理后，得到代表不同频率成分的各阶IMF分量imf

imf

式中，A

根据Hertz碰撞理论如式1所示，质量M取决于信号主频率f及碰撞速度V

M

式中，rms为信号有效值，N为采样点数，x

步骤三：质量估计函数模型建立及质量估计

采用平滑曲线拟合方法，得到由质量估计特征M

式中，a、b为待确定参数，通过平滑曲线拟合方法对数据进行拟合后即可得出a、b具体值。

令不同高度下质量估计函数为

其中，a

任选一高度下质量估计函数曲线，令最小质量冲击信号所得特征值为M

根据点(M

采集实时信号经步骤一与步骤二处理后，代入质量估计函数模型M

本发明具有以下有益效果：

(1)本发明所提质量估计特征能较好地描述信号随质量变化规律，且通过统一不同跌落高度质量估计曲线减小了跌落高度对质量估计影响，提高了质量估计精度；

(2)本发明通过计算实时信号特征并将其代入质量估计曲线即可得出质量估值，运算速度快，满足实时监测的要求；且质量估计曲线可根据后续不同质量松脱件冲击信号逐步完善，具有较好的实用性。

附图说明

图1为本方法流程图。

图2为本发明传感器布置示意图。

图3(a)(b)分别为为110g小球于74cm高度跌落所产生冲击信号波形及频谱。

图4(a)(b)分别为为110g小球于74m高度跌落所产生冲击信号的imf

图5(a)(b)(c)分别为10m、40m、74m不同跌落高度条件下，质量估计特征M

图6为不同跌落高度条件下，质量估计函数曲线。

图7为最终质量估计函数曲线。

具体实施方式

下面结合图和实施例对本发明做详细描述。

以松脱件模拟碰撞试验为实施例。使用标准椭圆封头模拟反应堆压力容器底部，封头直径为1800mm，厚度为10mm，材料为Q345R，封头下表面焊接三个刚性支座用于支撑。

如图1所示，一种基于EMD的松脱件质量估计方法，包括以下三个步骤：

步骤一：冲击信号采集与EMD处理

首先通过在试验封头底部布置4个加速度传感器实时采集不同质量松脱件冲击信号用于后续质量估计函数模型建立，记所采集冲击信号为A(t)；

基于Hertz理论，以小球冲击平板为例，小球质量M可表述为：

M(V

式中，f为信号主频率，k

由式可知，为准确估算小球质量，关键在于碰撞速度V

式中，imf

在本实施例中，传感器布置方式如图2所示。所使用小球质量分别为10g、47g、110g、175g、260g、375g，跌落高度分别为100mm、400mm和740mm。为避免偶然因素影响，每种条件下小球重复跌落10次。取110g小球74cm高度跌落条件为例，其原信号波形频谱如图3所示。

步骤二：主要IMF成分确定与特征提取

信号A(t)经EMD处理后，得到代表不同频率成分的各阶IMF分量imf

imf

式中，A

根据Hertz碰撞理论如式1所示，质量M取决于信号主频率f及碰撞速度V

M

式中，rms为信号有效值，N为采样点数，x

在本实施例中，取110g小球74cm高度跌落条件为例，imf

步骤三：质量估计函数模型建立及质量估计

采用平滑曲线拟合方法，得到由质量估计特征M

式中，a、b为待确定参数，通过平滑曲线拟合方法对数据进行拟合后即可得出a、b具体值。

令不同高度下质量估计函数为

其中，a

任选一高度下质量估计函数曲线，令最小质量冲击信号所得特征值为M

根据点(M

采集实时信号经步骤一与步骤二处理后，代入质量估计函数模型M

在本实施例中，使用matlab软件通过平滑曲线拟合方法对不同高度下数据进行拟合，得到式7函数模型中待确定参数a、b值，如表1所示，不同高度下质量估计函数模型如图6所示。

表1不同高度下质量估计函数模型待确定参数值

将各高度下质量估计函数模型统一后得到最终质量估计函数模型，其函数方程表述如下：

曲线如图7所示。

统计不同质量冲击信号所对应质量估计结果如表2所示，定义质量估计误差e如下式：

式中，m为估计质量值，m

表2不同质量冲击信号质量估计结果

以上所述为本发明的基本原理及实现形式，配合各图予以说明。本发明并不局限于实施例的限制，上述实施例和说明书所描述的为本发明原理，在不脱离本发明精神和范围前提下的变化和改进都落入要求保护的本发明范围内。