基于小波变换融合盲源分离算法的泄漏声波特征提取方法
【专利摘要】本发明公开了一种基于小波变换融合盲源分离算法的泄漏声波特征提取方法，包括以下步骤：在被测管道上设置传感器，通过传感器对泄漏点进行信号采集，获取泄漏声波采集信号；利用小波变换对泄漏声波采集信号进行预处理，获取观测信号，并对观测信号采用盲源分离算法进行处理，获取目标信号；对步骤二中的目标信号进行评价，并对观测信号组成进行优选。本发明的有益效果是：本发明提供的基于小波变换融合盲源分离算法的泄漏声波特征提取方法通过泄漏时刻采样点偏差和幅值损失两个评价参数对目标处理信号进行评价，该方法能够对泄漏时刻进行准确定位，同时对微弱信号的泄漏幅值的补偿作用明显。
【专利说明】
基于小波变换融合盲源分离算法的泄漏声波特征提取方法
技术领域
[0001]本发明涉及管道检测技术领域，尤其是一种基于小波变换融合盲源分离算法的泄漏声波特征提取方法。
【背景技术】
[0002]目前可以应用于油气管道的泄漏监测方法有许多种，其中，声波法与传统的质量平衡法、负压波法、瞬态模型法等相比具有诸多优点:灵敏度高、定位精度高、误报率低、检测时间短、适应性强;测量的是管线流体中的微弱动态压力变化量，与管线运行压力的绝对值无关;响应频率更宽，检测范围更宽等。
[0003]输气管道发生泄漏时产生声波信号，随着传播距离的增加泄漏信号产生衰减，波形特征被噪声覆盖，为提取有效的泄漏特征，国内外学者进行了大量的研究，根据调研结果，现阶段国内外涉及输气管道泄漏声波特征提取方法的专利主要有:
[0004]美国专利US6389881公开了一种利用音波技术进行管道泄漏检测的技术，该技术利用传感器采集管内动态压力，采用模式匹配滤波技术对信号进行滤波处理，排除噪声，降低干扰，提高了定位精度；
[0005]中国专利200710177617.0公开了一种基于压力和声波信息融合的泄漏检测方法，该方法分别采集管道上下游压力和声波信号(0.2-20HZ内)，经过数据滤波、特征级融合和决策级融合三个层次的处理获得最终检测结果，并利用基于相关分析、小波分析等融合的定位方法进行泄漏定位，提高了泄漏检测的准确性和定位精度。
[0006]中国专利201510020155.6公开了一种基于声波幅值的油气管道泄漏定位方法，该方法采用经过小波分析处理后得到低频段声波幅值来进行泄漏检测和定位，提出了一种不考虑声速及时间差的泄漏定位方法。
[0007]中国专利CN104614069A公开了一种基于联合近似对角化盲源分离算法的电力设备故障音检测方法，步骤包括:采用麦克风阵列;采用基于联合近似对角化盲源分离算法针对采用麦克风阵列采集的声音信号分离各个独立声源信号；提取独立声源信号的Mel频率倒谱系数MFCC作为声音特征参数，通过模式匹配算法识别声音信号，将待测试声音模板与所有的参考样本模板进行匹配后，匹配距离最小的参考样本模板就是电力设备工作音识别的结果。
[0008]现有的专利只是小波变换或者盲源分离算法单一处理方法的应用，对两种方法的融合技术没有描述，具体表现为:
[0009](I)小波变换能够提取低频段的信号特征，是应用最为普遍的信号处理方法，但同时小波变换在信号提取时也存在较为明显的缺陷，在实际应用中，低频段信号特征的获取需要对原始信号进行深层分解，在泄漏时刻的定位以及泄漏幅值的获得上容易出现较大偏差，容易造成时间差的计算误差，使得定位误差较大;泄漏幅值损失容易造成泄漏波形的失真，容易造成泄漏的漏判和误判。
[0010](2)为解决这一问题，采用盲源分离算法对信号进行处理，经研究发现盲源分离能够准确定位泄漏时刻，且在泄漏幅值方面不但没有损失，反而有所补偿，尤其是在信号较为微弱时补偿更为明显，但在应用时，盲源分离同样存在较为明显的缺陷:盲源分离得到的目标信号顺序、种类不能确定。
【发明内容】
[0011]本发明的目的是为克服上述现有技术的不足，提供一种基于小波变换融合盲源分离算法的泄漏声波特征提取方法。
[0012]为实现上述目的，本发明采用下述技术方案:
[0013]基于小波特征近似信号融合盲源分离算法的泄漏声波特征提取方法，包括以下步骤:
[0014]步骤一:在被测管道上设置传感器，通过传感器对泄漏点进行信号采集，获取泄漏声波采集信号；
[0015]步骤二:利用小波变换对泄漏声波采集信号进行预处理，获取多个细节信号，将泄漏声波采集信号和多个细节信号作为观测信号，并对观测信号采用盲源分离算法进行处理，获取目标信号；
[0016]步骤三，对步骤二中的目标信号进行评价，并对观测信号组成进行优选。
[0017]优选的，所述步骤一中，声波传感器采用动态压力传感器。
[0018]优选的，所述步骤二中，小波变换采用的小波基为sym8，分解层数根据传感器采集的泄漏声波采集信号中含有的信号种类决定，所述信号种类包括泄漏声波信号，背景噪声以及流动噪声。
[0019]优选的，所述背景噪声包括动力设备的运转，管道外部环境的噪声以及硬件设备、电路产生的噪声;所述流动噪声包括流体流动产生的瑞流噪声。
[0020]优选的，所述步骤二中，观测信号获取方法具体步骤如下:
[0021]步骤S201:将泄漏声波采集信号作为原始信号，确定小波分解层数N，N大于等于2，将原始信号作为待分解信号，进行小波分解，分解首次分别获取第一层细节信号和第一层近似信号；
[0022]步骤S202:将第一层近似信号作为待分解信号，对待分解信号进行小波分解，分别获取待分解信号对应的第二层细节信号和第二层近似信号；
[0023]步骤S203:将第N-1层近似信号作为待分解信号，重复执行步骤S202，直至达到分解层数N，第N-1层近似信号小波分解对应第N层细节信号和第N层近似信号；
[0024]步骤S204:选取第一层至第N层对应的各层细节信号以及第N层近似信号作为观测信号。
[0025]进一步优选的，所述根据信号种类确认分解层数的方法为:分解层数等于传感器采集的泄漏声波信号中含有的信号种类减I。
[0026]优选的，所述步骤二中，采用盲源分离算法进行处理获取目标信号的数目的方式有两种:一是目标信号的总数等于观测信号的总数，即当观测信号有m个，则目标信号也有m个;二是直接定义目标信号数目为一个，即当观测信号有m个，目标信号有且只有一个。
[0027]优选的，所述步骤三中，利用泄漏时刻采样点偏差和幅值损失作为评价参数。
[0028]所述泄漏时刻采样点偏差指目标信号的泄漏时刻采样点与原始信号的泄漏时刻采样点之差。所述泄漏时刻采样点偏差越小代表泄漏时刻定位更准确。
[0029]所述幅值损失指目标信号的泄漏幅值与原始信号的泄漏幅值的差值与原始信号泄漏幅值的绝对值之比。幅值损失为负值，且该值绝对值越大，代表幅值补偿越显著。
[0030]本发明的有益效果是:
[0031]1.本发明提供的基于小波变换融合盲源分离算法的泄漏声波特征提取方法通过泄漏时刻采样点偏差和幅值损失两个评价参数对目标处理信号进行评价，该方法能够对泄漏时刻进行准确定位，同时对微弱信号的泄漏幅值的补偿作用明显；
[0032]2.本发明解决了现阶段小波变换在泄漏时刻定位以及泄漏幅值偏差误差较大，盲源分离目标信号顺序、种类不能确定的问题，提高了声波法泄漏检测与定位技术的适用性；
[0033]3.本发明方法简单，操作方便，对提取油气管道声波法泄漏检测与定位方法中的泄漏声波特征适用性强。
【附图说明】
[0034]图1是本发明实施例提供的基于小波变换融合盲源分离算法的泄漏声波特征提取方法的原理图；
[0035]图2是本发明实施例提供的基于小波变换融合盲源分离算法的泄漏声波特征提取方法实施流程示意图；
[0036]图3是本发明实施例提供的基于小波变换融合盲源分离算法的泄漏声波特征提取方法处理前的原始信号示意图；
[0037]图4a是本发明实施例提供的基于小波变换融合盲源分离算法的泄漏声波特征提取方法处理后得到的第一个目标信号示意图；
[0038]图4b是本发明实施例提供的基于小波变换融合盲源分离算法的泄漏声波特征提取方法处理后得到的第二个目标信号示意图；
[0039]图4c是本发明实施例提供的基于小波变换融合盲源分离算法的泄漏声波特征提取方法处理后得到的第三个目标信号示意图；
[0040]图5是本发明实施例提供的基于小波变换融合盲源分离算法的泄漏声波特征提取方法处理后得到的一个目标信号示意图。
【具体实施方式】
[0041 ]下面结合附图和实施例对本发明进一步说明。
[0042]如图1所示，基于小波变换融合盲源分离算法的泄漏声波特征提取方法的流程图参考图1，本发明根据
【发明内容】
给出的技术方案采用下述实验参数对本发明进行实验验证:
[0043]实验参数如下:原始信号为2MPa下0.6mm泄漏孔径时距离泄漏点109m的传感器采集的信号，参考图3，分解层数N为2，采样频率f为3000Hz，小波变换分解函数为sym8或者db40
[0044]实施例1:该实施例中，采用的盲源分离算法进行处理获取目标信号的数目的方式为目标信号的总数等于观测信号。原始信号为A0，小波变换分解后得到的第一层观测信号为D1，第二层观测信号为D2和A2，即用于步骤二中盲源分离算法的观测信号为A2，D2，D1，再分别对A2，D2，D1进行盲源分离，得到目标信号。
[0045]如图3、图4a、图4b、图4c所不，图4表不米用本发明提供的方法获取的二个目标f目号。
[0046]实施例2:该实施例中，采用盲源分离算法进行处理获取目标信号的数目的方式为定义目标信号有且只有一个。图5表示采用本发明提供的方法获取的一个目标信号。
[0047]参考图4a、图4b和图4c，根据实施例1的实验附图可以看出，原始信号幅值为-5.06159kPa，原始信号泄漏时刻对应采样点为46441;由图4a可以看出，泄漏声波信号幅值为-9.53160kPa，泄漏声波信号泄漏时刻对应采样点为46441。经过发明得到的3个信号中与原始信号接近的目标信号为泄漏声波信号，其余2个分别为背景噪声和流动噪声，且泄漏声波信号、背景噪声和流动噪声的排序依次为1、2、3，所以，本发明将泄漏声波信号中含有的信号种类进行分类，并且对目标信号的顺序进行了确认。该目标信号数目的确定方式为优选方式，因为本发明考虑到实验过程中不仅要从泄漏声波采集信号中获取泄漏声波信号，还要对从泄漏声波采集信号中获取的背景噪声和流动噪声进行进一步的研究，为此，优选上述目标信号数目确认方式。
[0048]参考图5，根据实施例1的实验附图可以看出，泄漏声波信号幅值为-9.53160kPa，泄漏声波信号泄漏时刻对应的采样点为46441。根据上述数据可以看出，本发明提供的方法得到的目标信号的泄漏时刻采样点与原始信号的幅值偏差为0，且本发明得到的目标信号的泄漏声波幅值大于原始信号的幅值，因此本发明提供的方法对泄漏时刻定位较为准确，且幅值损失为-88.31%，即泄漏幅值补偿作用明显。若本发明仅仅从泄漏声波采集信号中获取泄漏声波信号，可采用实施例2采用的目标信号数目确定方式，目标信号定义为有且只有一个，所以不用考虑目标信号的排序、种类区分的问题。
[0049]综上所述，本发明提供的方法通过泄漏时刻采样点偏差和幅值损失两个评价参数对目标信号进行评价，该方法能够对泄漏声波的泄漏时刻进行准确定位，同时对微弱信号的泄漏幅值的补偿作用明显，因此，有效降低了现阶段小波变化在泄漏时刻误差较大的问题。
[0050]因此，本发明盲源分离后可以将目标信号进行有效的分类，进而提高了声波法泄漏检测与定位的实用性。
[0051]上述虽然结合附图对本发明的【具体实施方式】进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
【主权项】
1.基于小波变换融合盲源分离算法的泄漏声波特征提取方法，其特征是，包括以下步骤: 步骤一:在被测管道上设置传感器，通过传感器对泄漏点进行信号采集，获取泄漏声波米集?目号； 步骤二:利用小波变换对泄漏声波采集信号进行预处理，获取多个细节信号，将泄漏声波采集信号和多个细节信号作为观测信号，并对观测信号采用盲源分离算法进行处理，获取目标信号； 步骤三，对步骤二中的目标信号进行评价，并对观测信号组成进行优选。2.如权利要求1所述的基于小波变换融合盲源分离算法的泄漏声波特征提取方法，其特征是，所述步骤一中，声波传感器采用动态压力传感器。3.如权利要求1所述的基于小波变换融合盲源分离算法的泄漏声波特征提取方法，其特征是，所述步骤二中，小波变换采用的小波基为sym8，分解层数由传感器采集的泄漏声波采集信号中含有的信号种类决定，所述信号种类包括泄漏声波信号，背景噪声以及流动噪声。4.如权利要求3所述的基于小波变换融合盲源分离算法的泄漏声波特征提取方法，所述背景噪声包括动力设备的运转，管道外部环境的噪声以及硬件设备、电路产生的噪声;所述流动噪声包括流体流动产生的瑞流噪声。5.如权利要求3所述的基于小波变换融合盲源分离算法的泄漏声波特征提取方法，所述步骤二中，观测信号获取方法具体步骤如下: 步骤S201:将泄漏声波采集信号作为原始信号，确定小波分解层数N，N大于等于2，将原始信号作为待分解信号，进行小波分解，分解首次分别获取第一层细节信号和第一层近似信号； 步骤S202:将第一层近似信号作为待分解信号，对待分解信号进行小波分解，分别获取待分解信号对应的第二层细节信号和第二层近似信号； 步骤S203:将第N-1层近似信号作为待分解信号，重复执行步骤S202，直至达到分解层数N，第N-1层近似信号小波分解对应第N层细节信号和第N层近似信号； 步骤S204:选取第一层至第N层对应的各层细节信号以及第N层近似信号作为观测信号。6.如权利要求5所述的基于小波变换融合盲源分离算法的泄漏声波特征提取方法，所述根据信号种类确认分解层数的方法为:分解层数等于泄漏声波信号种类减I。7.如权利要求5所述的基于小波变换融合盲源分离算法的泄漏声波特征提取方法，所述目标信号数量定义方式为目标信号的总数等于观测信号的总数。8.如权利要求5所述的基于小波变换融合盲源分离算法的泄漏声波特征提取方法，所述目标信号数量定义方式为目标信号数目为一个。9.如权利要求1所述的基于小波变换融合盲源分离算法的泄漏声波特征提取方法，所述步骤三中，利用泄漏时刻采样点偏差和幅值损失作为评价参数。10.如权利要求8所述的基于小波变换融合盲源分离算法的泄漏声波特征提取方法，所述泄漏时刻采样点偏差指目标信号的泄漏时刻采样点与原始信号的泄漏时刻采样点之差；所述幅值损失指目标信号的泄漏幅值与原始信号的泄漏幅值的差值与原始信号泄漏幅值 的绝对值之比。
【文档编号】G06F17/14GK105927861SQ201610246962
【公开日】2016年9月7日
【申请日】2016年4月20日
【发明人】刘翠伟, 石海信, 梁金禄, 方丽萍, 李玉星, 张玉乾, 韩金珂, 梁杰
【申请人】中国石油大学(华东), 钦州学院