CN109668058B - 基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法 - Google Patents
基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法 Download PDFInfo
- Publication number
- CN109668058B CN109668058B CN201811585943.XA CN201811585943A CN109668058B CN 109668058 B CN109668058 B CN 109668058B CN 201811585943 A CN201811585943 A CN 201811585943A CN 109668058 B CN109668058 B CN 109668058B
- Authority
- CN
- China
- Prior art keywords
- linear prediction
- pipeline
- signal
- leakage
- water supply
- 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.)
- Active
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
本发明提供基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法,属于供水管网泄漏检测与定位技术领域。本发明首先采集管道未通水时的环境背景噪声信号,而后在同一环境背景下分别采集管道正常时的声音信号、管道在漏损时的声音信号;对采集的信号分别计算其lyapunov指数、短时过零率、线性预测倒谱系数LPCC,并建立B‑P神经网络;在待测管道上采集声音信号,分别计算其lyapunov指数、短时过零率、线性预测倒谱系数特征值并输入所建立的B‑P神经网络,进行漏损辨识。本发明解决了现有供水管道检漏技术依靠人的经验辨识、漏损辨识精度不高的问题。本发明可用于供水管道漏损精确辨识。
Description
技术领域
本发明涉及供水管道漏损辨识方法,属于供水管网泄漏检测与定位技术领域。
背景技术
水是生命之源、发展之本,水资源是关系国计民生的一项重要资源;我国城市的供水管网漏失率较高,据统计,我国600多个城市供水管网的平均漏损率(管网漏损率是管网漏水量与供水总量之比)超过15%,最高达70%以上,全国城市供水年漏损量近100亿立方米;而乡镇供水管道的漏失情况更为严重。
目前主流的主动漏失监测方法有基于流量监测的方法和基于声波振动信号监测的方法两种。
基于流量监测的漏失监测方法:
将检漏目标区域与外界连接的其中一个阀门上安装流量计,并关闭除此之外的其它所有阀门,则在深夜用水量最小时,就可以通过流量计来判断该区域内有无漏失。如果有,则通过关闭阀门进一步缩小区域范围,直到缩小至管段范围。如果现场条件不允许只开启一个阀门,则可以留一个进水口和一个出水口,并在这两处安装流量计。通过两个流量计及小区内用户水表的读数来判断该区域内有无漏失。
基于流量监测的方法的缺点则在于,通过逐步缩小检测区域的方法,可以最终找到漏失点的位置及漏失量的大小,但是这种方法,工作强度非常大,经常需要在深夜进行,而且在一个区域内频繁地开闭一些阀门,会导致管道内水流频繁发生变化,使管垢脱落,带来水质恶化风险。
基于声波振动信号监测的方法:
利用声波振动信号进行供水管网泄漏检测与定位是目前自来水公司采用的主要手段。管网漏失存在时,由于水流与管壁的摩擦,会产生声波振动信号,基于声波振动信号监测的方法即是通过对这些声学信号的监测进而发现管网漏失。这种方法又称为听音检漏法,是主动检漏法的一种,是指通过工人或仪器对管道进行“听诊”的方式而发现漏失的方法。这种方法最初是工人通过将听音棒连到管道构筑物(闸阀、消防栓等)或管道上方地面上来听取管道中水流声音的。有经验的工人可以根据听到的声音来判断有无漏失,以及漏失的大概位置。显然这种方法对工人的能力和经验提出了很高的要求,并且效率较低。
目前供水管道漏水检测的常用方式是先通过压力流量监测或用户反应确定漏损大致区域,而后通过关闭通水的阀门,确定漏损的大致区域,而后通过听漏杆和阀栓接触,有经验的工人用人耳判断该管段是否发生漏损,最后通过电子放大听漏仪确定漏点的准确位置。听漏杆目前在检漏队伍中几乎人手一根,用以确定管段是否漏水,缺点是工作人员需要不断的比较来搜索漏损管段,为了找到一个漏点需要花费较多的人力和时间,同时定位精度不高可能导致大片开挖管道,且使用效果较大程度依赖于工人的经验;物业小区等若发生漏损,非专业人员使用听漏杆等设备,无法对漏损进行判断。
发明内容
本发明为解决现有供水管道检漏技术依靠人的经验辨识、漏损辨识精度不高的问题,提供了基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法。
本发明所述基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法,通过以下技术方案实现:
步骤一、采集管道未通水时的环境背景噪声信号,而后在同一环境背景下分别采集管道正常时的声音信号、管道在漏损时的声音信号;
步骤二、对采集的信号分别计算其lyapunov指数、短时过零率、线性预测倒谱系数LPCC;lyapunov指数就是李雅谱诺夫指数;
步骤三、以步骤二中计算得到的lyapunov指数、短时过零率、线性预测倒谱系数特征,建立B-P神经网络;
步骤四、在待测管道上采集声音信号,分别计算其lyapunov指数、短时过零率、线性预测倒谱系数特征值;
步骤五、将步骤四中得到的待测管道上的特征值输入所建立的B-P神经网络,进行漏损辨识。
本发明最为突出的特点和显著的有益效果是:
本发明所涉及的基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法,通过采集信号的lyapunov系数、线性预测倒谱系数IPCC、短时过零率特征,建立B-P神经网络;能够智能化、快速、准确的对管道是否发生漏损进行判断,避免了现有检漏工作中对工人经验的依赖性过高的状况,使得非专业认识也能进行管道漏损辨识;并且本发明方法大大提高了辨识准确率,仿真实验对100组声音信号进行辨识,漏损辨识精度高达99%,节约了无效开挖的成本。
附图说明
图1为本发明中漏损信号采集结构示意图;
图2为本发明流程图;
图3为实施例中的神经网络训练集(R=0.99996)回归曲线图;其中,Date表示数据,Fit为拟合,Y指的实际输出;T为实际目标;
图4为实施例中的神经网络验证集(R=0.99726)回归曲线图;
图5为实施例中的神经网络测试集(R=0.99989)回归曲线图;
图6为实施例中的神经网络总体(R=0.99956)回归曲线图;
图7为实施例中的神经网络性能曲线图;
图8为实施例中采用本发明方法辨识漏损的预测结果与实际漏损情况对比图;
其中,1.管道,2.阀栓(阀门或者消火栓),3.加速度传感器,4.电荷放大器,5.动态采集分析仪,6.上位机。
具体实施方式
具体实施方式一:结合图1、图2对本实施方式进行说明,本实施方式给出的基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法,具体包括以下步骤:
步骤一、采集管道未通水时的环境背景噪声信号,而后在同一环境背景下分别采集管道正常时的声音信号、管道在漏损时的声音信号。采集过程如图1所示,将加速度传感器3与水管1上的阀栓2连接,加速度传感器3将收集到的声音信号转化为电信号,然后将电信号传递至电荷放大器4进行放大,而后通过动态采集分析仪5连接到上位机6进行后续分析;
步骤二、对采集的信号分别计算其lyapunov指数、短时过零率、线性预测倒谱系数LPCC;lyapunov指数就是李雅谱诺夫指数;
步骤三、以步骤二中计算得到的lyapunov指数、短时过零率、线性预测倒谱系数特征,建立B-P神经网络;
步骤四、在待测管道上采集声音信号,分别计算其lyapunov指数、短时过零率、线性预测倒谱系数特征值;
步骤五、将步骤四中得到的待测管道上的特征值输入所建立的B-P神经网络,进行漏损辨识。
具体实施方式二:本实施方式与具体实施方式一不同的是,所述lyapunov指数的计算过程具体包括:
A1、通过自相关函数法计算信号S的时间延迟τ;
A2、通过对信号S的傅里叶变换求其平均周期T′;
A3、对信号S进行关联维数c计算,进而确定嵌入维数m;
A4、利用时间延迟τ、平均周期T′、嵌入维数m对待测信号S进行相空间重构,得到重构后信号相空间Y(ti);i=0,…,n;其中,t0表示时间序列的起始点,tn为时间序列终点;
A5、计算重构后信号相空间的起始点Y(t0)和与其最邻近点Y0(t0)的距离L0;
A6、追踪A5中两点的时间演化,直至某时刻ti,Y(ti)和与其最邻近点Y0(ti)的距离Li超过阈值ε,计算此次迭代的LMD值:
A7、在Y(ti)附近另找一个点Y1(ti),计算Y(ti)与Y1(ti)之间的距离Li′,使得Li′≤ε,并且Li与Li′线段的夹角最小;
A8、重复A6、A7,直至到达时间序列终点,迭代总次数为M;
A9、计算M次LMD值的平均值,得到Lyapunov指数为:
Lyapunov指数是衡量系统动力学特性的一个重要定量指标,它表征了系统在相空间中相邻轨道间收敛或发散的平均指数率。对于系统是否存在动力学混沌,可以从最大Lyapunov指数是否大于零非常直观的判断出来:一个正的Lyapunov指数,意味着在系统相空间中,无论初始两条轨线的间距多么小,其差别都会随着时间的演化而呈指数率的增加以致达到无法预测,这就是混沌现象。管道漏损信号存在着一定的混沌现象,不同的漏口条件及管道条件会对泄漏造成较大差别,漏损信号与噪声信号的Lyapunov指数存在着差别,漏损信号的Lyapunov指数在通常在1.55~2.20之间,而空白信号的Lyapunov指数通常在0.85~1.65之间。
其他步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式二不同的是,步骤A3中维数m>2c+1。
其他步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一不同的是,所述短时过零率的具体计算过程包括:
首先对信号进行分帧处理,保证分帧后的数组大小一致;然后利用matlab中的sign函数(符号函数)计算每帧信号短时过零率。
在离散时间语音信号情况下,如果相邻的采样具有不同的代数符号就称为发生了过零。单位时间内过零的次数就称为过零率。在语音信号识别中用以识别清浊音,通短时过零率反应了信号的连续性、及短时间内的频率特性。漏损信号与空白信号及噪声信号的短时过零率存在差别,漏损信号的短时过零率在2200~3300之间,而空白信号的短时过零率在0~4000之间。
其他步骤及参数与具体实施方式一、二或三相同。
具体实施方式五:本实施方式与具体实施方式一至四不同的是,所述线性预测倒谱系数的具体计算过程包括:
B1、将信号S进行分帧处理,形成汉明窗为256*256的矩阵P1;
B2、将分帧后的矩阵P1进行行加窗;
B3、进行快速傅里叶变换(FFT变换),用自相关法计算每帧的12阶线性预测参数LPC;
B4、求LPC的倒谱,得到线性预测倒谱系数(LPCC)。
线性预测倒谱系数是基于声道模型的重要特征参数。LPCC是丢弃了信号生成过程中的激励信息,之后用十多个倒谱系数可以代表共振峰的特性,所以可以在语音识别中取得很好的性能。漏损信号能够在管道裸露处,被人耳所识别,其线性预测倒谱系数与空白信号及噪声信号在特征值上存在着差异性。
其他步骤及参数与具体实施方式一至四相同。
实施例
采用以下实施例验证本发明的有益效果:
(1)如图1所示,将加速度传感器连接到管道暴露处的阀栓上,加速度传感器与电荷放大器连接,放大倍数调制100pc/unit;再通过动态采集分析仪连接到上位机(笔记本电脑);动态采集分析仪采用江苏TAT5912,采样频率10KHZ,分析频率3.91KHZ;分别采集了背景噪声信号50组,管道正常时的声音信号50组,管道漏损时的声音信号100组。
(2)在笔记本电脑上利用matlab对信号处理:对采集到的200组信号计算lyapunov系数、LPCC线性预测倒谱系数、短时过零率特征;
(3)将漏损时的声音信号、正常时的声音信号、背景噪声信号整理,形成神经网络输入数据组;
(4)将漏损时的声音信号、正常时的声音信号、背景噪声信号添加标签,漏损信号标签为1,噪声信号标签为0,正常运行信号标签为2;
(5)对神经网络输入参数设置:输入矩阵为200行、14列的向量矩阵,隐层节点设置为10个,截止误差(神经网络训练停止的误差)设计为0.0001,最大迭代次数设置为300次;
(6)建立神经网络:
其中,神经网络训练集(R=0.99996)回归曲线图如图3所示;神经网络验证集(R=0.99726)回归曲线图如图4所示;神经网络测试集(R=0.99989)回归曲线图如图5所示;神经网络总体(R=0.99956)回归曲线图如图6所示;
如图7所示,当迭代次数M=9次,最小均方误差为0.00056562,能够满足辨识要求,同时R值(相关度)达到0.95以上,说明信号相关性良好。
(7)将50组实际漏水信号、25组噪声信号、及25组管道正常通水信号,共计100组待测信号(样本)输入神经网络进行漏损识别:
对神经网络输出结果进行归类,当输出值在(0,0.5)之间时,使输出值=0;
当输出值在[0.5,1.5]之间时,使输出值=1;当输出值在(1.5,2.5)之间时,使输出值=2;当输出值小于等于0或大于等于2.5时,输出值为0;
测试结果如图8所示,可以看出,本发明方法在漏损辨识准确率达到99%。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (4)
1.基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法,其特征在于,所述方法具体包括以下步骤:
步骤一、采集管道未通水时的环境背景噪声信号,而后在同一环境背景下分别采集管道正常时的声音信号、管道在漏损时的声音信号;
步骤二、对采集的信号分别计算其lyapunov指数、短时过零率、线性预测倒谱系数LPCC;lyapunov指数就是李雅谱诺夫指数;
所述lyapunov指数的计算过程具体包括:
A1、通过自相关函数法计算信号S的时间延迟τ;
A2、通过对信号S的傅里叶变换求其平均周期T′;
A3、对信号S进行关联维数c计算,进而确定嵌入维数m;
A4、利用时间延迟τ、平均周期T′、嵌入维数m对待测信号S进行相空间重构,得到重构后信号相空间Y(ti);i=0,...,n;其中,t0表示时间序列的起始点,tn为时间序列终点;
A5、计算重构后信号相空间的起始点Y(t0)和与其最邻近点Y0(t0)的距离L0;
A6、追踪A5中两点的时间演化,直至某时刻ti,Y(ti)和与其最邻近点Y0(ti)的距离Li超过阈值ε,计算此次迭代的LMD值:
A7、在Y(ti)附近另找一个点Y1(ti),计算Y(ti)与Y1(ti)之间的距离Li′,使得Li′≤ε,并且Li与Li′线段的夹角最小;
A8、重复A6、A7,直至到达时间序列终点,迭代总次数为M;
A9、计算M次LMD值的平均值,得到Lyapunov指数为:
步骤三、以步骤二中计算得到的lyapunov指数、短时过零率、线性预测倒谱系数,建立B-P神经网络;
步骤四、在待测管道上采集声音信号,分别计算其lyapunov指数、短时过零率、线性预测倒谱系数;
步骤五、将步骤四中得到的待测管道上的特征值输入所建立的B-P神经网络,进行漏损辨识。
2.根据权利要求1所述基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法,其特征在于,步骤A3中维数m>2c+1。
3.根据权利要求1所述基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法,其特征在于,所述短时过零率的具体计算过程包括:
首先对信号进行分帧处理,保证分帧后的数组大小一致;然后利用sign函数计算每帧信号短时过零率。
4.根据权利要求1~3任意一项所述基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法,其特征在于,所述线性预测倒谱系数的具体计算过程包括:
B1、将信号S进行分帧处理,形成汉明窗为256*256的矩阵P1;
B2、将分帧后的矩阵P1进行加窗;
B3、进行快速傅里叶变换,用自相关法计算每帧的12阶线性预测参数LPC;
B4、求LPC的倒谱,得到线性预测倒谱系数LPCC。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811585943.XA CN109668058B (zh) | 2018-12-24 | 2018-12-24 | 基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811585943.XA CN109668058B (zh) | 2018-12-24 | 2018-12-24 | 基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109668058A CN109668058A (zh) | 2019-04-23 |
CN109668058B true CN109668058B (zh) | 2021-01-29 |
Family
ID=66147235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811585943.XA Active CN109668058B (zh) | 2018-12-24 | 2018-12-24 | 基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109668058B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110543719B (zh) * | 2019-08-30 | 2022-09-27 | 哈尔滨工业大学 | 一种基于空间计量模型的供水管道漏失预测方法 |
WO2021111602A1 (ja) * | 2019-12-05 | 2021-06-10 | 日本電信電話株式会社 | ブロー成形不良検知装置、ブロー成形不良検知方法、プログラム |
CN113780381B (zh) * | 2021-08-28 | 2022-07-01 | 特斯联科技集团有限公司 | 一种人工智能漏水检测方法及装置 |
CN113670531B (zh) * | 2021-09-13 | 2023-12-01 | 哈尔滨工业大学 | 基于相位和幅值衰减的多探头阵列检测供水管道漏损的方法及系统 |
CN115950590B (zh) * | 2023-03-15 | 2023-05-30 | 凯晟动力技术(嘉兴)有限公司 | 气体发动机泄露预警系统 |
-
2018
- 2018-12-24 CN CN201811585943.XA patent/CN109668058B/zh active Active
Non-Patent Citations (1)
Title |
---|
基于非线性理论和信息融合的说话人识别;侯丽敏;《中国博士学位论文全文数据库 信息科技辑》;20070115(第1期);第I136-10页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109668058A (zh) | 2019-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109668058B (zh) | 基于线性预测倒谱系数和lyapunov指数的供水管道漏损辨识方法 | |
Song et al. | Gas leak detection in galvanised steel pipe with internal flow noise using convolutional neural network | |
Chuang et al. | Leak detection in water distribution pipes based on CNN with mel frequency cepstral coefficients | |
Zhang et al. | Deep convolutional neural network probability imaging for plate structural health monitoring using guided waves | |
CN108520753B (zh) | 基于卷积双向长短时记忆网络的语音测谎方法 | |
Wang et al. | Study on the MFCC similarity-based voice activity detection algorithm | |
CN101592288B (zh) | 一种管道泄漏辨识方法 | |
CN104595729A (zh) | 一种基于声波幅值的油气管道泄漏定位方法 | |
CN111520615B (zh) | 基于线谱对和三次插值搜索的管网漏损识别与定位方法 | |
CN102915728B (zh) | 声音分段设备和方法以及说话者识别系统 | |
CN106090630B (zh) | 基于集成分类器的液体管道检漏方法及其系统 | |
CN101996635A (zh) | 基于重音突显度的英语发音质量评价方法 | |
Mei et al. | Adaptive signal enhancement based on improved VMD-SVD for leak location in water-supply pipeline | |
CN101452698A (zh) | 一种自动嗓音谐噪比分析方法 | |
Guo et al. | Novel leakage detection and localization method based on line spectrum pair and cubic interpolation search | |
Han et al. | Localization of CO2 gas leakages through acoustic emission multi-sensor fusion based on wavelet-RBFN modeling | |
CN105909979B (zh) | 基于小波变换融合盲源分离算法的泄漏声波特征提取方法 | |
AU2021101586A4 (en) | A System and a Method for Non-Intrusive Speech Quality and Intelligibility Evaluation Measures using FLANN Model | |
Zhengjie et al. | Pressure vessel leakage detection method based on online acoustic emission signals | |
JP2002538514A (ja) | 周波数スペクトラムにおける確率論的信頼度を用いた音声検出方法 | |
Liu et al. | Gas leak detection system in compressor stations based on a microphone array and multi-channel frequency Transformer | |
Shao et al. | Segmented Kalman filter based antistrong transient impact method for vortex flowmeter | |
Han et al. | Identification of Pipeline Leak Sizes Based on Chaos-Grey Wolf-Support Vector Machine | |
CN110487911B (zh) | 基于盲源分离的压力容器声发射信号检测的方法 | |
Chen et al. | Large-scale analysis of formant frequency estimation variability in conversational telephone speech. |
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 |