CN111734961A - 一种天然气管道泄漏检测方法 - Google Patents

一种天然气管道泄漏检测方法 Download PDF

Info

Publication number
CN111734961A
CN111734961A CN202010588128.XA CN202010588128A CN111734961A CN 111734961 A CN111734961 A CN 111734961A CN 202010588128 A CN202010588128 A CN 202010588128A CN 111734961 A CN111734961 A CN 111734961A
Authority
CN
China
Prior art keywords
algorithm
gull
optimal
sound wave
imf
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010588128.XA
Other languages
English (en)
Other versions
CN111734961B (zh
Inventor
路敬祎
岳继康
王冬梅
董宏丽
刘庆强
刘远红
刘均
张彦生
宋金波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeast Petroleum University
Original Assignee
Northeast Petroleum University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northeast Petroleum University filed Critical Northeast Petroleum University
Priority to CN202010588128.XA priority Critical patent/CN111734961B/zh
Publication of CN111734961A publication Critical patent/CN111734961A/zh
Application granted granted Critical
Publication of CN111734961B publication Critical patent/CN111734961B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17D—PIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00—Protection or supervision of installations
    • F17D5/02—Preventing, monitoring, or locating loss
    • F17D5/06—Preventing, monitoring, or locating loss using electric or acoustic means
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Mechanical Engineering (AREA)
  • Examining Or Testing Airtightness (AREA)

Abstract

一种天然气管道泄漏检测方法，涉及管道泄漏检测技术领域，其解决了管道泄漏检测系统误报率高的问题。所述方法包含以下步骤：利用声学传感器采集声波信号；利用海鸥算法优化变分模态分解算法；利用基于优化的变分模态分解阈值去噪方法对所述的声波信号进行降噪预处理，获取去噪的声波信号；分别提取声波信号的时频特征构建高维特征向量矩阵，利用局部线性嵌入算法对所述的高维特征向量矩阵进行降维，提取有利于分类的敏感特征向量；利用海鸥优化算法搜索最佳的惩罚因子和核函数使最小二乘支持向量机的性能最优，将所述的敏感特征向量作为最小二乘支持向量机的训练样本，根据采集的声波信号作为测试样本，来检测天然气管道是否发生泄漏。

Description

一种天然气管道泄漏检测方法
技术领域
本发明涉及管道泄漏检测技术领域，具体而言，涉及一种天然气管道泄漏检测方法。
背景技术
管道运输是一种主要的运输方式，具有管理方便、运输成本低、稳定、安全等特点，是输送石油、天然气等流体资源的必然选择。随着天然气田的不断开发和天然气管网的不断扩大，管道腐蚀、老化、偷气现象日益严重。这些现象对输气管道的安全运行构成了潜在的严重威胁。因此，必须采取泄漏控制措施来检测天然气管道系统中的泄漏。气体管道泄漏可以通过多种检测方法进行评估，如质量平衡、负压波、瞬态模型、分布式光纤、声波法。声波法具有灵敏度高、定位精度高、误报率低、检测时间短、适应性强等优点。因此，我们使用这种检测方法来获取声波信号。变分模态分解算法是2014年Dragomiretskiy等人提出的一种自适应信号分解方法，然而，变分模态分解算法的参数的设置不准确会严重影响信号分解的结果。声波法采集的管道泄漏信号容易受到环境噪声影响，常常使得管道泄漏检测系统产生误报，漏报等问题。特征提取构造的特征向量中存在着变化相同或相似的冗余特征，如果直接将高维特征向量作为最小二乘支持向量机的输入，会加大训练难度，导致较高的误报率。选取不同的惩罚因子和核函数会得到不同的最小二乘支持向量机，如果以经验选取惩罚因子和核函数，可能会降低检测准确率。
文献号为CN111022941A的现有技术提供了一种基于LSTM循环神经网络的天然气管道泄漏检测方法，其包括以下步骤：步骤1：采集管道泄漏音频数据；步骤2：对所述音频数据采用改进的梅尔频率倒谱系数MFCC方法进行特征提取；步骤3：构建一个基于长短期记忆LSTM循环神经网络的天然气管道泄漏检测模型；步骤4：对所述管道微泄漏检测模型进行训练及参数优化，并验证模型的鲁棒性。该文献在管道泄漏检测模型的参数优化过程中，基于LSTM循环神经网络的天然气管道泄漏检测方法能够对管道状态进行检测，但该现有方法对变化相同或相似的冗余特征的剔除没有考虑，训练难度较大，影响训练时间。
发明内容
本发明提供了一种天然气管道泄漏检测方法，以解决现有管道泄漏检测系统产生误报率高、漏报等问题。
本发明为解决上述技术问题提供了如下技术方案：
利用声学传感器采集不同种类的声波信号，利用海鸥算法优化变分模态分解算法，使变分模态分解算法根据声波信号的特点自适应的选择最佳的分解层数和罚参数；利用基于优化的变分模态分解阈值去噪方法对所述的声波信号进行降噪预处理，获取去噪的声波信号；提取声波信号的时频特征构建高维特征向量矩阵，通过局部线性嵌入算法对所述的高维特征向量矩阵进行降维，提取有利于分类的敏感特征向量，降低向量机的训练难度；将所述降维后的敏感特征向量作为最小二乘支持向量机的训练样本，利用海鸥优化算法搜索最小二乘支持向量机的最优参数，使最小二乘向量机的性能最优，根据采集的声波信号作为测试样本，来检测天然气管道是否发生泄漏。
进一步地，利用海鸥算法优化变分模态分解算法，以排列熵的均值与方差的比值作为粒子群算法的适应度函数，搜索最优的信号分解层数和惩罚参数。海鸥算法优化变分模态分解算法的步骤如下：
步骤1：输入声波信号，设置变分模态分解算法需要优化的参数范围，初始化海鸥算法各项参数，包括种群规模N、最大迭代次数L以及位置X等。
步骤2：计算适应度值，每个海鸥的位置作为运行一次变分模态分解算法的参数，并求得该位置所有本征模态函数的适应度值。
步骤3：更新海鸥最佳位置和适应度，如果当前海鸥的适应度值小于上次迭代的最佳适应度值，则用当前海鸥的适应度值代替上次迭代的最佳适应度值，用当前海鸥的位置代替上次迭代的最佳位置，否则最佳适应度值和最佳位置不变。
步骤4：重复步骤2～3，进行循环迭代，直到达到最大迭代次数后输出最优的分解层数和惩罚参数。
进一步地，适应度值的计算公式如下：
Figure BDA0002555416560000021
Figure BDA0002555416560000022
其中，HP(IMFk)为本征模态函数(IMF)的排列熵；mean(HP(IMFk))为所有排列熵的均值；var(HP(IMFk))为所有排列熵的方差；对IMFk进行相空间重构排序，再映射为符号序列，每种符号序列出现的概率即为P(.)，N为P(.)的个数。
优选地，基于变分模态分解的阈值去噪方法，其步骤如下：
步骤1：首先利用优化后的变分模态分解算法对所述的声波信号进行分解，得到若干本征模态函数(IMF)；
步骤2：计算每个本征模态函数的方差密度和平均周期之积；
Figure BDA0002555416560000031
Figure BDA0002555416560000032
Figure BDA0002555416560000033
其中，E(.)为数学期望，M为IMFn长度，
Figure BDA0002555416560000034
为IMFn的平均周期，K为IMFs的个数。
步骤3：将满足以下条件的第一个IMF视为跳变点：
Figure BDA0002555416560000035
步骤4：计算所有IMF在跳变点之前的平均值：
Figure BDA0002555416560000036
步骤5：当满足以下条件时，确定非有效IMF：
Figure BDA0002555416560000037
对非有效IMF进行软阈值处理：
Figure BDA0002555416560000038
其中,j表示非有效IMF的个数，M为信号的长度。
步骤6：对非有效IMF进行阈值处理后，对所有IMF求和，得到去噪的声波信号
优选地，所述的时频特征值包含：均值、绝对均值、标准方差、均方根值、方根幅值、绝对最大值、最大值、最小值、峰峰值、偏度、峭度、峰值因子、裕度因子、波形因子、脉冲因子。由上述时频域特征值组成高维特征向量，通过局部线性嵌入算法对所述的高维特征向量进行降维，提取有利于分类的敏感特征向量。
进一步地，利用海鸥优化算法搜索最小二乘支持向量机的最优参数，所述需要优化的参数包含惩罚因子和核函数，将所述的敏感特征向量作为优化后最小二乘支持向量机的训练样本。其优化步骤如下：
步骤1：初始化海鸥种群，最小二乘支持向量机的惩罚因子和核函数参数作为海鸥的位置X，种群规模设置为N，最大迭代次数为L。
步骤2：计算适应度值，对训练集进行交叉验证意义下的准确率作为海鸥优化算法的适应度函数值。
步骤3：更新海鸥最佳位置和适应度，如果当前海鸥的适应度值小于上次迭代的最佳适应度值，则用当前海鸥的适应度值代替上次迭代的最佳适应度值，用当前海鸥的位置代替上次迭代的最佳位置，否则最佳适应度值和最佳位置不变。
步骤4：海鸥算法完成搜索的终止条件为循环达到最大迭代次数，设置的迭代次数为L。
步骤5：当所述的适应度函数完成收敛，即对训练集进行交叉验证意义下的准确率达到最大时，输出最优的惩罚因子和核函数。
本发明具有以下有益效果：利用海鸥算法优化变分模态分解算法，解决了变分模态算法参数设置不准确导致信号分解效果不佳的问题。利用基于变分模态分解的阈值去噪方法对所述的声波信号进行降噪预处理，有效解决了环境噪声对管道泄漏检测造成干扰的问题。通过局部线性嵌入算法对高维特征向量进行降维，提取有利于分类的敏感特征向量，剔除了变化相同或相似的冗余特征，降低了最小二乘支持向量机的训练难度，减少了训练时间。通过海鸥优化算法搜索最小二乘支持向量机最优的惩罚因子和核函数参数，提高了最小二乘支持向量机的预测准确率。该检测方法简单、准确率高、检测成本低。
附图说明
为了更为清楚的说明本发明的技术方案，下面对本发明的所需要的附图进行说明。
图1为本发明所述的天然气管道泄漏检测方法的流程框图；
图2为本发明中海鸥算法优化变分模态分解算法(步骤S100的子步骤)的流程框图；
图3为本发明中管道声波信号去噪(步骤S110的子步骤)流程框图；
图4为本发明中提取管道声波信号敏感特征的(步骤S120的子步骤)流程框图；
图5为本发明中海鸥优化最小二乘支持向量机的(步骤S130的子步骤)流程框图；
图6为本发明中检测管道工况的(步骤S140的子步骤)流程框图；
图7为本发明实施例提供的管道声波信号去噪效果图；图中：(a)为泄漏信号去噪效果图，(b)为敲击信号去噪效果图。
图8为本发明实施例提供的时频域特征向量变化图；
图9为本发明实施例提供的局部线性嵌入降维后提取的敏感特征向量图；
图10为本发明实施例提供的海鸥算法优化最小二乘支持向量机的适应度收敛曲线及最佳惩罚因子和核函数；
图11为本发明实施例提供的优化后最小二乘支持向量机的预测分类与实际类别对比图。
具体实施方式
下面将结合附图对本发明的具体实施方式进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为一种天然气管道泄漏检测方法流程实施例，用于对天然气管道的工况进行检测，所述方法包含以下步骤：
步骤S100:利用声学传感器采集管道的声波信号,利用海鸥算法优化变分模态分解算法；
本实施例中，采集的声波信号包括但不局限于正常信号、敲击信号、泄漏信号。利用海鸥算法优化变分模态分解算法，以排列熵的均值与方差的比值作为粒子群算法的适应度函数，搜索最优的信号分解层数和惩罚参数。
步骤S110:利用基于变分模态分解的阈值去噪方法对所述的声波信号进行降噪预处理，获取去噪的声波信号；
本实施例中，通过对所述的声波信号进行降噪预处理，可以降低管道内其他噪声对有效声波信号的干扰，最终获得去噪的声波信号。
步骤S120:通过局部线性嵌入算法提取声波信号的敏感特征向量；
本实施例中，所述的正常信号、敲击信号、泄漏信号时频曲线之间存在明显差异，但众多特征向量存在变化相同或相似的冗余特征，因此局部线性嵌入算法提取声波信号的敏感特征向量，可以降低向量机的训练难度，提高检测准确率。
步骤S130:利用海鸥算法优化最小二乘支持向量机；
本实施例中，利用海鸥优化算法迭代搜索最优的惩罚因子和核函数，训练的过程就是寻找最佳参数组合使向量机的性能最优。
步骤S140:检测天然气管道的运行工况；
本实施例中，将采集的声波信号通过处理后作为测试集，输入到优化后的最小二乘支持向量机里面，判断天然气管道的运行工况。
请参照图2，本实施例中，进一步的，所述步骤S100可以包括如下子步骤：
步骤200：输入声波信号，设置参数范围，初始化海鸥算法；
本实施例中，输入声学传感器采集的声波信号，设置变分模态分解算法需要优化的参数范围，初始化海鸥算法各项参数，包括种群规模N、最大迭代次数L以及位置xi等。
步骤210：循环海鸥位置，运行变分模态分解算法，计算适应度值；
本实施例中，每个海鸥的位置作为运行一次变分模态分解算法的参数，并求得该位置所有本征模态函数的适应度值。
所述适应度值的计算公式如下：
Figure BDA0002555416560000061
Figure BDA0002555416560000062
其中，HP(IMFk)为本征模态函数(IMF)的排列熵；mean(HP(IMFk))为所有排列熵的均值；var(HP(IMFk))为排列熵的方差；对IMFk进行相空间重构排序，再映射为符号序列，每种符号序列出现的概率即为P(.),N为P(.)的个数。
步骤220：更新海鸥最佳位置和适应度；
本实施例中，更新海鸥最佳位置和适应度，如果当前海鸥的适应度值小于上次迭代的最佳适应度值，则用当前海鸥的适应度值代替上次迭代的最佳适应度值，用当前海鸥的位置代替上次迭代的最佳位置，否则最佳适应度值和最佳位置不变。
步骤230：输出最优的分解层数和罚参数；
本实施例中，重复步骤2～3，进行循环迭代，直到达到最大迭代次数后输出最优的分解层数和罚参数。
请参照图3，本实施例中，进一步的，所述步骤S110可以包括如下子步骤：
步骤S300:对所述的声波信号进行变分模态分解，获取分解后的本征模态函数(IMF)；
本实施例中，本征模态函数从低频到高频依次排列，其中，包含有效本征模态函数和非有效本征模态函数。
步骤S310:计算每个本征模态函数的方差密度和平均周期之积；
本实施例中，本征模态函数的方差密度和平均周期之积为常数，
Figure BDA0002555416560000071
Figure BDA0002555416560000072
Figure BDA0002555416560000073
其中，E(.)为数学期望，M为IMFn长度，
Figure BDA0002555416560000074
为IMFn的平均周期，K为IMFs的个数。
步骤S320:确定跳变点；
本实施例中，将满足以下条件的第一个IMF视为跳变点：
Figure BDA0002555416560000075
步骤S330:计算所有IMF在跳变点之前的平均值；
本实施例中，平均值的计算公式如下：
Figure BDA0002555416560000076
步骤S340:确定非有效IMF，对非有效IMF进行软阈值处理；
本实施例中，当满足以下条件时，确定非有效IMF：
Figure BDA0002555416560000077
对非有效IMF进行软阈值处理：
Figure BDA0002555416560000078
其中，j表示非有效IMF的个数，M为信号的长度。
步骤S350:对非有效IMF进行阈值处理后，对所有IMF求和，得到去噪的声波信号。经过上述步骤得到去噪后的声波信号如图6所示。
请参照图4，本实施例中，进一步的，所述步骤S120可以包括如下子步骤：
步骤S400:提取声波信号的时频域特征值，构建高维特征向量；
本实施例中，所述的时频域特征包含均值、绝对均值、标准方差、均方根值、方根幅值、绝对最大值、最大值、最小值、峰峰值、偏度、峭度、峰值因子、裕度因子、波形因子、脉冲因子。需要说明的是，所述特征值还可以包含其他与泄漏检测有关的指标，这里不做特别限定，由这些特征值构建高维特征向量，所述的时频域特征值的变化图如图7所示。
步骤S410:通过局部线性嵌入算法对所述高维特征向量进行降维，提取有利于分类的敏感特征向量；
本实施例中，所述的局部线性嵌入算法是一种非线性降维方法，能将高维特征模型映射到低维空间，剔除变化特点相同或相似的冗余的特征参数，提取对工况识别有利的敏感特征。所述提取的敏感特征值如图8所示，其中红色“+”表示敲击信号的特征值，蓝色“Δ”表示泄漏信号特征值，黑色“◇”表示正常信号特征值，由图可知，同一管道工况的数据聚集在一起，而且管道工况特征值分类明显，有利于后续的分类。
请参照图5，本实施例中，进一步的，所述步骤S130可以包括如下子步骤：
步骤S500:初始化海鸥种群；
本实施例中，最小二乘支持向量机的惩罚因子和核函数参数作为海鸥的位置，随机初始化海鸥的位置X，种群规模设置为N，最大迭代次数为L。
步骤S510:计算适应度值；
本实施例中，对训练集进行交叉验证意义下的准确率作为海鸥优化算法的适应度函数值。
步骤S520:更新海鸥最佳位置和适应度；
本实施例中，如果当前海鸥的适应度值大于上次迭代的最佳适应度值，则用当前海鸥的适应度值代替上次迭代的最佳适应度值，用当前海鸥的位置代替上次迭代的最佳位置，否则最佳适应度值和最佳位置不变。
步骤S530:是否满足终止条件；
本实施例中，海鸥算法完成搜索的终止条件为循环达到最大迭代次数，设置的迭代次数为L。
步骤S540:输出最优的惩罚因子和核函数；
本实施例中，当所述的适应度函数完成收敛，即对训练集进行交叉验证意义下的准确率达到最大时，输出最优的惩罚因子和核函数。图9为海鸥算法优化最小二乘支持向量机的适应度收敛曲线及最佳惩罚因子和核函数，由图可以看出，所述适应度函数在69代完成收敛，所述适应度函数达到最高为93.3333％，惩罚因子为26.6648，核函数为0.16308。
请参照图6，本实施例中，进一步的，所述步骤S140可以包括如下子步骤：
步骤S600:利用最佳参数训练最小二乘支持向量机；
本实施例中，利用步骤S440得到的最佳惩罚因子和核函数来训练最小二乘支持向量机，得到所述最小二乘支持向量机最优分类面函数。
步骤S610:采集的声波信号作为测试集，判断管道的运行工况；
本实施例中，天然气管道采集的声波信号作为测试集，输入到训练好的最小二乘支持向量机里面，如检测结果为正常时，结果输出为0；当检测为敲击时，结果输出为1；当检测结果为泄漏时，结果输出为2。图10为优化后最小二乘支持向量机的预测分类与实际类别对比图，由图可知，本发明对天然气管道的运行工况具有较高的检测准确率。
综上所述，首先利用声波传感器采集不同种类的声波信号，利用海鸥算法优化变分模态分解，使变分模态分解算法根据信号的特点自适应的选的最佳的分解层数和罚参数，解决了变分模态分解参数设置不准确导致的信号分解效果不佳的问题；利用基于变分模态分解算法的阈值去噪方法对声波信号进行去噪，有效解决了环境噪声对管道泄漏检测造成干扰的问题；通过局部线性嵌入算法对高维特征向量进行降维，提取有利于分类的敏感特征，剔除了变化相同或相似的冗余特征，降低了最小二乘支持向量机的训练难度，减少了训练时间。通过海鸥优化算法搜索最小二乘支持向量机最优的惩罚因子和核函数参数，解决了人为经验选取参数的弊端，提高了最小二乘支持向量机的预测准确率。该检测方法简单、准确率高、检测成本低。
以上所述仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均包含在本发明的保护范围内。

Claims (6)

1.一种天然气管道泄漏检测方法，其特征在于，所述方法包括：
利用声学传感器采集声波信号，利用海鸥算法优化变分模态分解算法，使变分模态分解算法根据声波信号的特点自适应的选择最佳的分解层数和罚参数；
利用基于变分模态分解的阈值去噪方法对所述的声波信号进行降噪预处理，获取去噪的声波信号；
提取声波信号的时频特征构建高维特征向量矩阵，利用局部线性嵌入算法对所述的高维特征向量矩阵进行降维，提取有利于分类的敏感特征向量，降低向量机的训练难度；
将所述降维后的敏感特征向量作为最小二乘支持向量机的训练样本，利用海鸥优化算法搜索最小二乘支持向量机的最优参数，使最小二乘向量机的性能最优，根据采集的声波信号作为测试样本，来检测天然气管道是否发生泄漏。
2.根据权利要求1所述的天然气管道泄漏检测方法，其特征在于，变分模态分解算法的参数的设置不准确会影响信号分解的结果，因此，利用海鸥算法优化变分模态分解算法，以排列熵的均值与方差的比值作为粒子群算法的适应度函数，搜索最优的信号分解层数和惩罚参数；海鸥算法优化变分模态分解算法的步骤如下：
步骤1：输入声波信号，设置变分模态分解算法需要优化的参数范围，初始化海鸥算法各项参数，包括种群规模N、最大迭代次数L以及位置X；
步骤2：计算适应度值，每个海鸥的位置作为运行一次变分模态分解算法的参数，并求得该位置所有本征模态函数的适应度值；
步骤3：更新海鸥最佳位置和适应度，如果当前海鸥的适应度值小于上次迭代的最佳适应度值，则用当前海鸥的适应度值代替上次迭代的最佳适应度值，用当前海鸥的位置代替上次迭代的最佳位置，否则最佳适应度值和最佳位置不变；
步骤4：重复步骤2～3，进行循环迭代，直到达到最大迭代次数后输出最优的分解层数和罚参数。
3.根据权利要求2所述的天然气管道泄漏检测方法，其特征在于，适应度函数的计算公式如下：
Figure FDA0002555416550000011
Figure FDA0002555416550000012
其中，HP(IMFk)为本征模态函数(IMF)的排列熵；mean(HP(IMFk))为所有排列熵的均值；var(HP(IMFk))为排列熵的方差；对IMFk进行相空间重构排序，再映射为符号序列，每种符号序列出现的概率即为P(.),N为P(.)的个数。
4.根据权利要求1所述的天然气管道泄漏检测方法，其特征在于，基于变分模态分解的阈值去噪方法，其步骤如下：
步骤1：首先利用优化后的变分模态分解算法对所述的声波信号进行分解，得到若干本征模态函数(IMF)；
步骤2：计算每个本征模态函数的方差密度和平均周期之积；
Figure FDA0002555416550000021
Figure FDA0002555416550000022
Figure FDA0002555416550000023
其中，E(.)为数学期望，M为IMFn长度，
Figure FDA0002555416550000024
为IMFn的平均周期，K为IMFs的个数。
步骤3：将满足以下条件的第一个IMF视为跳变点：
Figure FDA0002555416550000025
步骤4：计算所有IMF在跳变点之前的平均值：
Figure FDA0002555416550000026
步骤5：当满足以下条件时，视为非有效IMF：
Figure FDA0002555416550000027
对非有效IMF进行软阈值处理：
Figure FDA0002555416550000028
其中，j表示非有效IMF的个数，M为信号的长度。
步骤6：对非有效IMF进行阈值处理后，对所有IMF求和，得到去噪的声波信号。
5.根据权利要求1所述的天然气管道泄漏检测方法，其特征在于，利用海鸥优化算法搜索最小二乘支持向量机的最优参数，将所述的敏感特征向量作为最小二乘支持向量机的训练样本，所述需要优化的参数包含惩罚因子和核函数；优化步骤如下：
步骤1：初始化海鸥种群，最小二乘支持向量机的惩罚因子和核函数参数作为海鸥的位置X，种群规模设置为N，最大迭代次数为L；
步骤2：计算适应度值，对训练集进行交叉验证意义下的准确率作为海鸥优化算法的适应度函数值；
步骤3：更新海鸥最佳位置和适应度，如果当前海鸥的适应度值大于上次迭代的最佳适应度值，则用当前海鸥的适应度值代替上次迭代的最佳适应度值，用当前海鸥的位置代替上次迭代的最佳位置，否则最佳适应度值和最佳位置不变；
步骤4：海鸥算法完成搜索的终止条件为循环达到最大迭代次数，设置的迭代次数为L。
步骤5：当所述的适应度函数完成收敛，即对训练集进行交叉验证意义下的准确率达到最大时，输出最优的惩罚因子和核函数。
6.根据权利要求1所述的天然气管道泄漏检测方法，其特征在于，所述的时频特征包含：均值、绝对均值、标准方差、均方根值、方根幅值、绝对最大值、最大值、最小值、峰峰值、偏度、峭度、峰值因子、裕度因子、波形因子、脉冲因子；由上述时频域特征值组成高维特征向量，通过局部线性嵌入算法对所述的高维特征向量进行降维。
CN202010588128.XA 2020-06-24 2020-06-24 一种天然气管道泄漏检测方法 Active CN111734961B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010588128.XA CN111734961B (zh) 2020-06-24 2020-06-24 一种天然气管道泄漏检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010588128.XA CN111734961B (zh) 2020-06-24 2020-06-24 一种天然气管道泄漏检测方法

Publications (2)

Publication Number Publication Date
CN111734961A true CN111734961A (zh) 2020-10-02
CN111734961B CN111734961B (zh) 2021-09-17

Family

ID=72650983

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010588128.XA Active CN111734961B (zh) 2020-06-24 2020-06-24 一种天然气管道泄漏检测方法

Country Status (1)

Country Link
CN (1) CN111734961B (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363094A (zh) * 2020-11-17 2021-02-12 西北工业大学 一种基于海鸥算法的三分量磁传感器误差校正方法
CN112858461A (zh) * 2021-02-01 2021-05-28 南昌航空大学 一种基于优化变分模态分解的巴克豪森信号消噪方法
CN113110024A (zh) * 2021-04-08 2021-07-13 扬州大学 一种基于椭圆焦半径改进型海鸥算法的风力机叶片振动控制器设计方法
CN113364315A (zh) * 2021-05-08 2021-09-07 南昌大学 一种基于海鸥优化算法的三相逆变器控制方法
CN113393909A (zh) * 2021-06-28 2021-09-14 广西民族大学 化工动态优化问题混合海鸥优化方法、系统、计算机设备
CN113432805A (zh) * 2021-07-13 2021-09-24 南昌航空大学 一种供水管道泄漏的无线监测装置及方法
CN113932158A (zh) * 2021-10-25 2022-01-14 东北石油大学 一种长输管道泄漏信号的采集传输装置
CN114263855A (zh) * 2021-11-19 2022-04-01 合肥工业大学 一种天然气运输管道泄露的预测方法及其应用
CN114638170A (zh) * 2022-03-31 2022-06-17 河北工业大学 基于数据挖掘技术的电能计量装置故障诊断方法及系统
CN115204316A (zh) * 2022-09-07 2022-10-18 中国人民解放军陆军装甲兵学院 一种基于支持向量机的车载网络故障诊断方法
CN115234850A (zh) * 2022-07-20 2022-10-25 西安科技大学 一种基于声波振动的光纤传感器气体泄漏检测方法
CN115640511A (zh) * 2022-12-24 2023-01-24 武汉商启网络信息有限公司 一种基于支持向量机的管道泄露检测方法及装置
CN115935241A (zh) * 2023-01-09 2023-04-07 国机传感科技有限公司 一种多参数相互融合的清管器实时定位方法及装置
CN117093899A (zh) * 2023-08-22 2023-11-21 山东大学 基于相异差和双流扩维图的燃气管网泄漏检测方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008078323A1 (en) * 2006-12-27 2008-07-03 Eltav Wireless Monitoring Ltd. Device and system for monitoring valves
WO2009135639A1 (en) * 2008-05-05 2009-11-12 Aktiebolaget Skf Determining hydrostatic leakage
CN108870090A (zh) * 2018-06-22 2018-11-23 大连理工大学 基于最小二乘支持向量机信息融合的管道泄漏检测方法
CN109654384A (zh) * 2019-01-29 2019-04-19 南京工业大学 基于pso-vmd算法的管道泄漏检测装置及检测方法
CN109681789A (zh) * 2019-03-01 2019-04-26 重庆邮电大学 基于变分模态分解的流体管道泄漏声振动相干检测定位方法
CN111024433A (zh) * 2019-12-30 2020-04-17 辽宁大学 一种改进灰狼算法优化支持向量机的工业装备健康状态检测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008078323A1 (en) * 2006-12-27 2008-07-03 Eltav Wireless Monitoring Ltd. Device and system for monitoring valves
WO2009135639A1 (en) * 2008-05-05 2009-11-12 Aktiebolaget Skf Determining hydrostatic leakage
CN108870090A (zh) * 2018-06-22 2018-11-23 大连理工大学 基于最小二乘支持向量机信息融合的管道泄漏检测方法
CN109654384A (zh) * 2019-01-29 2019-04-19 南京工业大学 基于pso-vmd算法的管道泄漏检测装置及检测方法
CN109681789A (zh) * 2019-03-01 2019-04-26 重庆邮电大学 基于变分模态分解的流体管道泄漏声振动相干检测定位方法
CN111024433A (zh) * 2019-12-30 2020-04-17 辽宁大学 一种改进灰狼算法优化支持向量机的工业装备健康状态检测方法

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363094A (zh) * 2020-11-17 2021-02-12 西北工业大学 一种基于海鸥算法的三分量磁传感器误差校正方法
CN112858461A (zh) * 2021-02-01 2021-05-28 南昌航空大学 一种基于优化变分模态分解的巴克豪森信号消噪方法
CN113110024A (zh) * 2021-04-08 2021-07-13 扬州大学 一种基于椭圆焦半径改进型海鸥算法的风力机叶片振动控制器设计方法
CN113364315A (zh) * 2021-05-08 2021-09-07 南昌大学 一种基于海鸥优化算法的三相逆变器控制方法
CN113393909A (zh) * 2021-06-28 2021-09-14 广西民族大学 化工动态优化问题混合海鸥优化方法、系统、计算机设备
CN113432805A (zh) * 2021-07-13 2021-09-24 南昌航空大学 一种供水管道泄漏的无线监测装置及方法
CN113932158A (zh) * 2021-10-25 2022-01-14 东北石油大学 一种长输管道泄漏信号的采集传输装置
CN114263855B (zh) * 2021-11-19 2024-04-26 合肥工业大学 一种天然气运输管道泄漏的预测方法及其应用
CN114263855A (zh) * 2021-11-19 2022-04-01 合肥工业大学 一种天然气运输管道泄露的预测方法及其应用
CN114638170A (zh) * 2022-03-31 2022-06-17 河北工业大学 基于数据挖掘技术的电能计量装置故障诊断方法及系统
CN115234850A (zh) * 2022-07-20 2022-10-25 西安科技大学 一种基于声波振动的光纤传感器气体泄漏检测方法
CN115204316A (zh) * 2022-09-07 2022-10-18 中国人民解放军陆军装甲兵学院 一种基于支持向量机的车载网络故障诊断方法
CN115640511A (zh) * 2022-12-24 2023-01-24 武汉商启网络信息有限公司 一种基于支持向量机的管道泄露检测方法及装置
CN115935241A (zh) * 2023-01-09 2023-04-07 国机传感科技有限公司 一种多参数相互融合的清管器实时定位方法及装置
CN115935241B (zh) * 2023-01-09 2023-05-12 国机传感科技有限公司 一种多参数相互融合的清管器实时定位方法及装置
CN117093899A (zh) * 2023-08-22 2023-11-21 山东大学 基于相异差和双流扩维图的燃气管网泄漏检测方法及系统

Also Published As

Publication number Publication date
CN111734961B (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
CN111734961B (zh) 一种天然气管道泄漏检测方法
CN107590516B (zh) 基于光纤传感数据挖掘的输气管道泄漏检测识别方法
CN111853555B (zh) 一种基于动态过程的供水管网暗漏识别方法
CN101592288B (zh) 一种管道泄漏辨识方法
CN115169479A (zh) 污水处理过程远程监控方法、系统及存储介质
CN113904786B (zh) 一种基于线路拓扑分析和潮流特性的虚假数据注入攻击辨识方法
CN113670616B (zh) 一种轴承性能退化状态检测方法及系统
CN111222583B (zh) 一种基于对抗训练与关键路径提取的图像隐写分析方法
CN113343123B (zh) 一种生成对抗多关系图网络的训练方法和检测方法
CN113343402B (zh) 基于多层卷积稀疏编码的管道腐蚀等级评估方法
CN114662327A (zh) 一种用于auv舱内泄露监测的辅助方法
CN114547608A (zh) 一种基于降噪自编码核密度估计的网络安全态势评估方法
CN117092980B (zh) 一种基于大数据的电气故障检测控制系统
CN116680639A (zh) 一种基于深度学习的深海潜水器传感器数据的异常检测方法
CN113987910A (zh) 一种耦合神经网络与动态时间规划的居民负荷辨识方法及装置
CN116629686A (zh) 一种评估企业能耗数据的方法及装置
CN116502091A (zh) 一种基于lstm和注意力机制的网络入侵检测方法
CN113407425B (zh) 基于BiGAN与OTSU的内部用户行为检测方法
CN114553468A (zh) 一种基于特征交叉与集成学习的三级网络入侵检测方法
CN114186583A (zh) 一种储油罐罐壁腐蚀检测异常信号恢复方法及系统
CN116701948B (zh) 管道故障诊断方法及系统、存储介质和管道故障诊断设备
CN114118308B (zh) 基于约束能量最小化变分自编码的高光谱目标检测方法
CN116428529A (zh) 一种基于AlexNet卷积神经网络的供暖管道泄漏检测方法
CN117031194B (zh) 用于配电网的超声波隐患检测方法及系统
CN114579981B (zh) 一种跨域漏洞检测方法、系统、存储介质和电子设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant