一种分布式光纤传感器振动信号分类方法及识别分类系统与流程

本发明属于信号检测技术领域，具体涉及一种分布式光纤传感器振动信号分类方法及识别分类系统。

背景技术：

分布式光纤传感扰动检测技术，因其抗电磁干扰能力强、灵敏度高、动态范围大、便于快速捕捉和高精度监测动态变化信息等技术优点，日益成为入侵行为安全监测的主流技术。例如，在石油管线中，利用分布式光纤扰动监测系统代替人工巡逻，及时发现偷油，漏油，破坏管线等行为，智能化管理管线，确保管线安全。在安防领域中，由于光纤的隐蔽性与抗电磁性，分布式光纤扰动监测系统十分适合实时监测各类突发事件，且监测成功率较普通光电系统高。

现有的光纤周界系统的振动信号识别方法主要有两种，一种是采用基于入侵振动信号的时域特征进行分析，如振动信号幅度或者过零率等。该方法由于进行模式识别的特征参数数量少，因此不能准确地区分各种外界振动信号，误报率高。另一种是采用信号的频谱技术，根据多维的频谱特征作为判别的依据。该方法计算量大，算法复杂，实时性较差。并且，分布式光纤传感器的信号是基于光在光纤内干涉形成，并不能很好的适合已有模型。

技术实现要素：

为解决上述问题，本发明提出一种针对光纤传感振动信号识别方法及系统，实现了对光纤传感信号的特征提取，准确识别分类，解决了光纤传感入侵信号识别分类准确率低的技术问题。

本发明采用如下技术方案，一种分布式光纤传感器振动信号分类方法，具体步骤如下：

s1：光纤电缆探测由外界振动产生的信号，使光信号携带外界振动信息，通过光电探测器得到振动信号；

s2：对振动信号进行归一化、去噪、分帧和加窗的预处理；

s3：计算振动信号的短时能量和短时过零率；

s4：依据噪声信号的短时能量、短时过零率和振动信号的短时能量、短时过零率的差异，设定短时能量、短时过零率门限；若振动信号的短时能量和短时过零率超过设定的门限，则判断是扰动事件，截取有效数据段，进入s5，否则返回s1；

s5：绘制振动信号时频域上的谱图；

s6：提取振动信号的mel频率倒谱系数；

s7：建立基于扰动事件信号的mel频率倒谱系数和时频域上谱图的深度学习识别模型；

s8：基于振动信号的mel频率倒谱系数和时频域上的谱图与步骤7的识别模型进行匹配。

优选地，s1中所述光纤传感振动信号是由外界振动产生的振动信号采集来的，采集的方法为相位解调光信号，使光信号携带外界振动信息，通过光电探测器得到振动信号。

优选地，所述s6中提取振动信号mel频率倒谱系数的步骤如下：

61)对预处理后的振动信号的每一帧做快速傅里叶变换，得到每一帧的频率谱图；

62)用梅尔滤波器每一帧的频率谱图进行功率谱滤波，并求出滤波后每帧的能量总和；

63)求出每帧能量经滤波后的对数；

64)保留滤波器的系数对每帧能量的对数做离散余弦变换。

优选地，所述s7中基于扰动事件信号的mel频率倒谱系数和时频域上的谱图训练深度学习识别模型的步骤如下：

71)对光纤传感系统实际采集到的扰动事件的信号进行分类，并标签和编码。

72)对每个扰动事件信号分别提取mel频率倒谱系数并绘制时频域的谱图。

73)固定每个扰动事件信号时频域谱图的尺寸得到一个三维数组，并在第三维上扩展mel频率倒谱系数的参数，将mel频率倒谱系数的参数从二维拉伸成一维，根据固定尺度第二维的维度铺满整个第三维，得到多个n*n*4的三维特征数组。

74)将提取好的特征数组按扰动事件类别进行分类，并标签和编码制成一个训练集。

75)构建对光纤传感振动信号识别系统的网络识别模型，使用训练集对光纤传感振动信号识别系统的网络模型进行训练，使之满足每个扰动事件信号与预设类别相符。

优选地，所述s8基于振动信号的mel频率倒谱系数和时频域上的谱图与识别模型进行匹配。

81)对传入的扰动事件信号分别提取mel频率倒谱系数并绘制时频域的谱图。

82)固定每个扰动事件信号时频域谱图的尺寸得到一个三维数组，并在第三维上扩展mel频率倒谱系数的参数，将mel频率倒谱系数的参数从二维拉伸成一维，根据固定尺度第二维的维度铺满整个第三维，得到多个n*n*4的三维特征数组。

83)将提取好的特征数组使用预先训练的网络进行分类。

一种针对光纤传感振动信号识别分类的系统，系统包括：光纤信号探测模块和光纤信号识别模块，其中，

所述光纤信号探测系统包括光纤光缆，光电探测器和数据采集卡，光纤光缆携带的振动信息通过光电探测器转换为电信号，经数据采集卡传回主机。

所述光纤信号识别系统包括信号处理模块，特征提取模块和信号分类识别模块，电信号传入信号处理模块对原始信号做预处理，预处理后的信号的短时能量和短时过零率超过阈值的话会传入特征提取模块，提取传入信号的mel频率倒谱系数和时频域上的谱图，并将两类特征融合，最后将融合的特征传入信号分类识别模块。

发明所达到的有益效果：本发明是一种针对光纤传感振动信号识别方法及系统，实现了对入侵信号行为种类的判断并做出及时的预警，解决了光纤传感入侵信号识别分类准确率低的技术问题。本发明参考了声音信号在低频域和光纤振动信号波形振动相似的特性，利用已经成熟的mel频率倒谱系数特征提取法，并佐之时频域上的谱图，将两个特征共同作为卷积神经网络的特征参数进行分类识别。实现了对入侵信号行为种类识别精度的提升，在安防领域具有较强的普适性。

附图说明

图1是本发明实施例的一种针对光纤传感振动信号识别方法的流程示意图；

图2是提取mel频率倒谱系数的流程示意图；

图3是本发明实施例的一种针对光纤传感振动信号识别分类的系统的模块示意图。

具体实施方式

下面根据附图并结合实施例对本发明的技术方案作进一步阐述。

图1是本发明实施例的一种针对光纤传感振动信号识别方法的流程示意图，一种分布式光纤传感器振动信号分类方法，包括以下步骤：

s1：光纤光缆探测由外界振动产生的信号，设备采样率和参数的设置以光缆铺设的具体位置来调整，采集到的信号有非入侵信号和入侵信号，还有环境噪声的干扰。

s2:对采集到的信号做预处理，根据实际的采样率和采集到的离散信号的点数确定合适的分帧窗函数，由于信号频率集中在低频域，要预加重信号的高频部分，再对其做小波去噪处理。

s3：对分帧后的振动信号分别计算短时能量和短时过零率，公式如下：

短时能量：

其中m表示窗的起点，n表示帧的长度，w(m)是窗函数，x(n)是一帧传感系统输出信号。

短时过零率：

其中，sgn[]为符号函数，

x(n)是一帧传感系统输出信号。

s4：依据噪声信号的短时能量、短时过零率和振动信号的短时能量、短时过零率的差异，设定短时能量、短时过零率门限；若振动信号的短时能量和短时过零率超过设定的门限，则判断是扰动事件，截取有效数据段，进入s5，否则返回s1；

s5：利用傅里叶变换绘制待处理信号时频域上的谱图。

s6：提取振动信号的mel频率倒谱系数，如图2所示：

提取振动信号的mel频率倒谱系数的步骤如下：

61)对预处理后的振动信号的每一帧做快速傅里叶变换，得到每一帧信号的频率谱图；

62)采用13个三角滤波器组成滤波器组，每个滤波器都对每帧的频率谱图进行滤波，梅尔尺度的转换公式为m(f)＝1125ln(1+f/700.0)，再把每帧的每个滤波器滤波后的能量加起来。

63)基于人的听觉原理，使用倒数均值法求出每帧的每个滤波器能量的对数。

64)滤波器组是有重叠的，每帧使用滤波器所得的能量是有相关的，使用离散余弦变换对能量做去相关处理。

s7：基于扰动事件信号的mel频率倒谱系数和时频域上的谱图训练深度学习识别模型；

71)根据铺设实际情况采集到的各类信号进行标签分类，截取信号短时能量和短时过零率都超出阈值的部分，利用短时傅里叶变换对已分类的信号批处理绘制时频域上的谱图，将其保存为jpg格式，并依次分类和编号；

72)再对分类后的信号固定帧数提取mel频率倒谱系数，每帧通过13个三角滤波器，固定提取保留有效信息的50帧的参数，将其保存为txt或者npy格式；

73)将信号的谱图和mel频率倒谱系数一一对应，谱图固定尺寸大小为896*1184*3，并在第三维上扩展mel频率倒谱系数的参数，将二维的参数拉伸成一维数组并铺满整个第三维，每帧信号提取了一个896*1184*4的数组，制成一个训练集。用pytorch的框架进行模型训练，光纤传感振动信号识别系统的网络识别模型共有七层，第一、第三、第五层为卷积层，第二和第四层为降采样层，第六层全连接层，最后一层为输出层。其中，卷积层通过多个滤波器提取三维特征数组的多个特征图谱，降采样层是利用了图像特征的局部相关性原理，减小了数据量的计算，提取了特征图谱的多个子图。全连接层计算输入向量和权重向量之间的点积，再加上偏置。利用softmax判断该训练集属于哪个类别。

s8：基于振动信号的mel频率倒谱系数和时频域上的谱图与识别模型进行匹配。

81)对传入的扰动事件信号分别提取固定帧数的mel频率倒谱系数并绘制时频域的谱图。

82)固定每个扰动事件信号时频域谱图的尺寸得到一个三维数组，并在第三维上扩展mel频率倒谱系数的参数，将mel频率倒谱系数的参数从二维拉伸成一维，根据固定尺度第二维的维度铺满整个第三维，得到n*n*4的三维特征数组。

83)使用预先训练的网络对提取好的特征数组进行分类。

一种针对光纤传感振动信号识别分类系统，系统包括：光纤信号探测模块和光纤信号识别模块如图3所示，其中，

所述光纤信号探测系统包括光纤光缆，光纤传感模块，光电探测器和数据采集卡，光纤携带的振动信息通过光电探测器转换为电信号，经数据采集卡传回主机。

所述光纤信号识别系统包括信号处理模块，特征提取模块和信号分类识别模块，经光纤传感系统采集到的光纤传感振动信号传入信号处理模块对原始信号做预处理，预处理后的信号的短时能量和短时过零率超过阈值的话会传入特征提取模块，提取传入信号的频率倒谱系数和时频域上的谱图，并将两类特征融合，最后将融合的特征信息送入信号分类识别模块。

以上所述仅为本发明的优选实施方式，用以说明本发明的技术方案，而非对其限制；应当指出：对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。