| 主权项 |
核电站松动部件多目标三维动态跟踪方法，其特征在于包括以下步骤：步骤(1)数据库的建立；1‑1由于松动部件主要发生在压力容器顶部或底部，在采集数据之前，根据反应堆压力容器顶部和底部的曲面划分多个均等的网格小分块，并对网格分块进行编号，然后在压力容器顶部和底部各安装三个呈等边三角形分布的加速度传感器；1‑2通过数据采集卡、电荷放大器和加速度传感器获取力锤敲击压力容器某一固定点的冲击信号C(t<sub>i</sub>)，用于测定传播速度；具体的：加速度传感器获取力锤敲击信号后，经电荷放大器进入数据采集卡；对冲击信号C(t<sub>i</sub>)进行基于小波变换的希尔伯特变换，得到冲击信号C(t<sub>i</sub>)到达某个加速度传感器的传播时间t，测量固定点到该加速度传感器的弧长分别为L，弧长除以传播时间，得到传播速度V；1‑3.测量各网格分块中心点到三个加速度传感器的弧长S<sub>i</sub>，然后将弧长S<sub>i</sub>除以传播速度V，得到敲击点到每个加速度传感器的时间T<sub>i</sub>，然后获取任意一个加速度传感器与另两个加速度传感器时间差<img file="FDA0000921329120000011.GIF" wi="195" he="91" />然后以网格编号，时间差<img file="FDA0000921329120000012.GIF" wi="58" he="83" />和<img file="FDA0000921329120000013.GIF" wi="66" he="88" />为查询表的属性，建立定位数据库；其中i＝1,2,3；T<sub>i</sub>＝S<sub>i</sub>/V；<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><msub><mi>T</mi><mn>21</mn></msub><mo>‾</mo></mover><mo>=</mo><msub><mi>T</mi><mn>2</mn></msub><mo>-</mo><msub><mi>T</mi><mn>1</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000921329120000014.GIF" wi="299" he="83" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><msub><mi>T</mi><mn>31</mn></msub><mo>‾</mo></mover><mo>=</mo><msub><mi>T</mi><mn>3</mn></msub><mo>-</mo><msub><mi>T</mi><mn>1</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000921329120000015.GIF" wi="299" he="83" /></maths>步骤(2)三维虚拟建模；根据核反应堆一回路系统中设备的实际图纸，利用三维建模软件建立核反应堆一回路系统中各个子模块的3D模型，并利用纹理贴图技术、特征简化技术和模型轻量化技术尽可能地减小文件、提高系统的运算速度；该3D模型可以透视，旋转，从不同角度观察内部部件；步骤(3)松动部件现场信号获取；通过数据采集卡、电荷放大器和加速度传感器来实现现场信号x(t<sub>i</sub>)的采集，现场信号x(t<sub>i</sub>)包括松动件跌落时的冲击信号和环境背景噪声；步骤(4)松动部件定位；4‑1对现场信号x(t<sub>i</sub>)进行频域分析确定信号主频率f<sub>m</sub>，用连续小波变换计算得到尺度因子a；选数值2、8、14、20中最接近尺度因子a，作为最终连续小波变换尺度对现场信号x(t<sub>i</sub>)进行处理，然后通过Hilbert包络线法求得信号的到达时间差t<sub>21</sub>、t<sub>31</sub>；4‑2将t<sub>21</sub>、t<sub>31</sub>与数据库中对应尺度下的定位查询表中的时间差作比较，时间差距离最小值所对应的网格即为定位网格；时间差距离的计算公式为：<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>d</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><mover><mrow><msub><mi>T</mi><mn>21</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><mo>‾</mo></mover><mo>-</mo><msub><mi>t</mi><mn>21</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mover><mrow><msub><mi>T</mi><mn>31</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><mo>‾</mo></mover><mo>-</mo><msub><mi>t</mi><mn>31</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>,</mo><mrow><mo>(</mo><mn>1</mn><mo>≤</mo><mi>i</mi><mo>≤</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000921329120000021.GIF" wi="963" he="99" /></maths>式中：i为网格编号；n为网格划分总数；d<sub>n×1</sub>的最小值所对应的i即为定位网格；步骤(5)松动部件质量估计5‑1对现场信号x(t<sub>i</sub>)出现最大幅值出现前0.03s开始，取0.3s长信号数据g(t<sub>i</sub>)，对信号g(t<sub>i</sub>)进行连续小波变换处理得到信号m(t<sub>i</sub>)后，计算信号m(t<sub>i</sub>)的小波能量谱，根据小波能量谱图找出该图尺度峰值x<sub>a</sub>；5‑2将尺度峰值x<sub>a</sub>带入已经建立好的尺度峰值函数F(x<sub>a</sub>)，求得松动部件的质量估计值y，然后对多个通道中的松动部件质量估计值求平均，得到最终的松动部件质量；步骤(6)动态跟踪显示；LPMS以时间T(该值可调，一般采用1h)为周期，周期地采集核电站一回路中的松动部件的现场信号；在T<sub>i</sub>时刻采集松动部件的现场信号，其中i＝0,1,2,3……；根据步骤(4)定位该松动部件的位置，根据步骤(5)估计该松动部件的质量，同时在压力容器曲面的定位位置处显示该松动部件；在T<sub>i+1</sub>时刻，再次采集松动部件的现场信号，其中T<sub>i+1</sub>＝T<sub>i</sub>+T，i＝0,1,2,3……；根据步骤(4)定位该松动部件的位置，根据步骤(5)估计该松动部件的质量，若该质量与先前任意周期内松动部件质量估计值不同，则为新增定位目标，同时在压力容器曲面的定位位置处显示此新增松动部件；若该质量与先前任意周期内松动部件质量估计值相同，则在压力容器曲面的定位位置处显示该松动部件，并标记该定位目标的移动轨迹；如此周期循环，动态的构成松动部件跌落轨迹。 |
