双谱分析技术在管道监测系统应用研究.docVIP

双谱分析技术在管道监测系统应用研究 　　【摘 要】本文利用双谱分析技术，分析基于光时域反射探测技术的管道监测系统中的常态信号、铲车激励信号和人激励行为的信号；通过双谱分析得到不同激励行为的分析图谱，基于双谱分析的图谱，设计针对不同激励行为的识别算法流程，实现了对不同激励行为的高准确识别率和极低的误识别率，其识别准确性99.5%，误识别率为0%，具有较大的实际应用价值。 　　【关键词】光时域反射探测技术；管道监测；谱分析；识别率 　　1 概述 　　20世纪以来国内经济高速发展，管道施工建设爆发式增长，但地下管道安全管理缺乏规范性，导致违章占压和新增设施与燃气管线安全间距不足形成安全隐患；油气管道与密闭空间、人员密集区、市政和民用管道、公路、铁路、河流等各类目标交叉穿越，安全距离不符合规范要求；同时前期建设的管网老化、腐蚀严重、存在很大安全隐患[1、7～9]。随着管道安全事故的频发和影响恶劣，基于光时域反射探测技术的管道安全监测系统越来越受到广泛关注，而目前基于光时域反射探测技术存在信号弱，噪声非线性等问题。专利（201110214678.6）主要涉及了差分法和多小波综合判断进行定位识别和检测；专利（201110214679.0）中涉及短时傅里叶变换的方法对事件进行监测；以上专利中主要提出了光时域反射探测技术中的线性噪声的影响，可以对事件进行定位，但未涉及事件行为类型的判断，而实际工程中，系统信号呈现多元化，复杂化，并存在相互交叉耦合等问题。本文利用双谱分析技术分析基于光时域反射探测技术的管道监测系统中的常态信号、铲车激励信号和人激励行为的信号；通过双谱分析得到不同激励行为的分析图谱，基于双谱分析的图谱，设计针对不同激励行为的识别算法流程，实现了对不同激励行为的高准确识别率和极低的误识别率。 　　2 谱分析原理 　　2.1 高阶谱定义 　　实际中高斯过程的信号向量服从高斯分布，相应的，其高阶累积量为零，而高阶矩却不一定。因此，信号分析时常常使用的是高阶累积量及其谱，高阶累积量谱又简称为高阶谱。下面给出高阶谱的定义。 　　假设随机过程x（t）的k阶累积量函数 绝对可和，即 　　或者 　　， 则k阶累积量谱 存在且连续，定义为k阶累积量的k-1维FFT： 　　其中 。 　　实际应用计算中，高阶累计量及其谱的值都需要估计，随着k值得增大，阶次增加，估计计算量也将大大增加，因此一般只考虑 ，即功率谱、双谱和三谱这三种情况。 　　2.2 双谱分析 　　在高阶谱中，双谱的阶数最低，处理方法最简单，但包含高阶谱的所有特性。双谱能够用来恢复、提取信号的相位信息，还能够系统的非线性特性。此外，高斯过程的三阶累积量为零，而非高斯过程则未必。因此，双谱可用来抑制高斯有色噪声的干扰，还有助于判断信号过程的类型。双谱从更高阶的概率结构上来表征随机信号，能够弥补二阶统计知识在这些方面的不足[2～6]。 　　在上面定义了双谱，这节就离散的确定性信号和含噪声的有限能量信号给出双谱的定义。 　　假设 为一离散的确定性信号，其双谱定义为 　　假设 为含噪声的有限能量信号，其双谱定义为 　　其中 为信号频率 的分量。 　　定义式中的双谱 是一个复数，它同时包括幅值 和相位 。 　　3 应用分析 　　3.1 双谱分析方法 　　在双谱分析方法中可分为参数法和非参数法，参数法要基于一个模型（AR模型、MA模型、ARMA模型），其基本思想是将待估计的数据序列通过选定的模型，然后估计出模型的参数，最后利用估计的参数进行双谱估计。 　　其中非参数法是用离散傅立叶变换（DFT）实现的，可分为直接方法（基于观测数据的傅立叶变换）和间接方法（基于观测数据的自相关函数的傅立叶变换）两种。 　　利用双谱分析DVS系统的信号的应用中主要利用间接法计算双谱，具体计算过程如下： 　　（1）将待估计的信号x（n），n=0，1，…，N-1，分成K段，每段的数据长度为M。数据段之间存在重叠部分，因此K≥N/M。第i段数据表示为xi（n），其中i=0，1，…，K-1，n=0，1，…，M-1； 　　（2）每段数据取平均值，并用该段的实际数据减去均值 　　（3）估计每段数据的3阶累计量值； 　　其中s1=max（0，-m，-n），s2=min（M-1，M-1-m，M-1-n）。 　　（4）计算所有数据段的三阶累积量的均值，得到整个信号的三阶累积量估计值； 　　（5）计算双谱估计值 　　式中LM-1， 是二维的窗函数。 　　3.2 信号双谱分析 　　依据以上介绍的双谱估计方法，分别对基于光时域反射探测技术的管道检测系统的常态、人激励、铲车激励等几种常见的第三方事件行为的双谱分析，具体分析结果如下： 　　注：以下图中x轴和y轴的表示频率，z轴