一种地震预警信息处理方法

技术领域

本发明属于地震预警信息处理领域，具体涉及判断地震预警信号是否正常播放的方法。

背景技术

地震预警，是指在地震发生后，利用地震波传播速度小于电波传播速度的特点，提前对地震波尚未到达的地方进行预警。地震预警信息的播报对于减少公众生命财产因地震所带来的损失具有重大意义。

公众接收来自于地震预警播报终端的信息，因此地震预警信息在终端是否正常播报直接影响地震放在工作。但是目前研究或专利技术中尚缺乏地震预警信息在终端是否正常播报的反馈与判断机制。

发明内容

本发明是鉴于上述的不足而完成的，其目的在于提供一种地震预警信息处理方法，反馈与判断地震预警信息在终端是否正常播报。

为了实现上述目的，本发明采用了以下技术方案：

首先通过拾音器收集预警播报终端的地震预警信息。

对收集的地震预警信号S(t)进行第一次分数阶傅里叶变换(FRFT)。在信号的时频域内进行坐标旋转，设旋转角度为α时，所采集信号的最大值带宽范围相较于其他旋转角度而言最小，则根据

其中，j为虚数单位。

将信号变换到分数阶域后，找出信号变换后的最大值位置和带宽范围，在以最优阶次p变换得到的分数阶域内，设横轴为X

在分数阶域内，对采集的地震预警信号进行分解，获得多重信号分量，分解方法采用经验模态(EMD)分解；

计算各信号分量与原采集信号的能量比值并降序排列，信号能量采用积分的形式获得，设原信号为S(X

按排列顺序叠加计算能量比值之和，当能量占比叠加值首次超过预设值时，记录参与叠加的n重信号分量。

对能量占比最高的第一重信号分量以外的n-1重信号分量分别进行FRFT变换，变换过程的最优阶次计算方法与第一次分数阶傅里叶变换时的方法一致，各分量经过变换后，找出其最大值位置和带宽范围。

将记录的n重信号分量的最大值位置和带宽范围与预存的正常地震预警信号进行对比，对比差值小于等于预设阈值时判断为地震预警信号正常播放，对比差值大于预设阈值时判断为地震预警信号未正常播放。将判断结果反馈地震监测与预警系统。

附图说明

通过结合附图和下面实施例的具体描述，本发明的实施过程和优点将更加直观和易于理解，其中：

图1为一种地震预警信息处理方法流程图；

图2为地震预警信号原始时频图；

图3为分数阶域内EMD分解分量图；

图4为处理完成后的地震预警信号时频图。

具体实施方式

下面结合附图对本发明的实施过程做进一步说明。

本发明的一种地震预警信息处理方法流程如图1所示，主要步骤包括：

首先通过拾音器收集预警播报终端的地震预警信息。

对采集的地震预警信号S(t)进行第一次分数阶傅里叶变换(FRFT)。在信号的时频域内进行坐标旋转，设旋转角度为α时，所采集信号的最大值带宽范围相较于其他旋转角度而言最小，则根据

其中，j为虚数单位。

下一步将信号变换到分数阶域后，找出信号变换后的最大值位置和带宽范围，在以最优阶次p变换得到的分数阶域内，设横轴为X

随后在分数阶域内，对采集的地震预警信号进行分解，获得多重信号分量，分解方法采用经验模态(EMD)分解；

计算各信号分量与原采集信号的能量比值并降序排列，信号能量采用积分的形式获得，设原信号为S(X

按排列顺序叠加计算能量比值之和，当能量占比叠加值首次超过预设值时，记录参与叠加的n重信号分量。

对能量占比最高的第一重信号分量以外的n-1重信号分量分别进行FRFT变换，变换过程的最优阶次计算方法与第一次分数阶傅里叶变换时的方法一致，各分量经过变换后，找出其最大值位置和带宽范围。

将记录的n重信号分量的最大值位置和带宽范围与预存的正常地震预警信号进行对比，对比差值小于等于预设阈值时判断为地震预警信号正常播放，对比差值大于预设阈值时判断为地震预警信号未正常播放。将判断结果反馈地震监测与预警系统。

实施例1：

通过拾音器收集预警播报终端的地震预警信息，其在时频域内的分布如图2 所示。信号中含有较多的噪声信号，不利于直接进行分析。这里首先对采集的地震预警信号S(t)进行第一次分数阶傅里叶变换(FRFT)，通过在信号的时频域内进行坐标旋转，确定旋转角度为13°时，所采集信号的最大值带宽范围相较于其他旋转角度而言最小，则根据

其中，j为虚数单位。

将信号变换到分数阶域后，找出信号变换后的最大值位置和带宽范围，在以最优阶次p变换得到的分数阶域内，设横轴为X

随后在分数阶域内，对采集的地震预警信号进行经验模态(EMD)分解，获得多重信号分量，结果如图3所示，其中横坐标进行了归一化处理。

计算各信号分量与原采集信号的能量比值并降序排列，得到的结果如表1 所示。能量占比为1.62％的剩余信号由于占比较低，未在图3和表1中列出。

表1信号分量的能量比值

按排列顺序叠加计算能量比值之和，预设能量占比叠加值大于等于80％时，记录参与叠加的信号分量，从表中可以看出前两重信号的能量占比叠加值已超过预设值。

对第二重信号分量信号分量进行FRFT变换，变换过程的最优阶次计算方法与第一次分数阶傅里叶变换时的方法一致，旋转角度为19°时，所采集信号的最大值带宽范围相较于其他旋转角度而言最小，则根据

将记录的第一重和第二重信号分量进行其最大值带宽范围的滤波，保留各最大值所在带宽的信号，并对各滤波后的分量进行分数阶傅里叶逆变换，变换到时域内并叠加，以此获得处理后的地震预警信号，如图4所示。

将记录的第一重和第二重信号分量的最大值位置和带宽范围与预存的正常地震预警信号进行对比，对比结果如表2所示，结果的差值均在10％的范围内，判定地震预警信号正常播放，将判断结果反馈地震监测与预警系统。

表2地震预警信号是否正常播放对比表

尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。