Nothing Special   »   [go: up one dir, main page]

CN108332063B - 一种基于互相关的管道泄漏定位方法 - Google Patents

一种基于互相关的管道泄漏定位方法 Download PDF

Info

Publication number
CN108332063B
CN108332063B CN201810083681.0A CN201810083681A CN108332063B CN 108332063 B CN108332063 B CN 108332063B CN 201810083681 A CN201810083681 A CN 201810083681A CN 108332063 B CN108332063 B CN 108332063B
Authority
CN
China
Prior art keywords
cross
sensors
leakage
pipeline
correlation
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
Application number
CN201810083681.0A
Other languages
English (en)
Other versions
CN108332063A (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.)
Institute of Acoustics CAS
Original Assignee
Institute of Acoustics CAS
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 Institute of Acoustics CAS filed Critical Institute of Acoustics CAS
Priority to CN201810083681.0A priority Critical patent/CN108332063B/zh
Publication of CN108332063A publication Critical patent/CN108332063A/zh
Application granted granted Critical
Publication of CN108332063B publication Critical patent/CN108332063B/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

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

本发明公开了一种基于互相关的管道泄漏定位方法，基于下述检测装置实现：在待检测管道上安装两个传感器，传感器为加速度传感器或者水听器；每个传感器各连接数据采集仪，用于同步采集管道泄漏时引起的流体中传播的泄漏声信号；所述方法包括：获取两个传感器拾取的泄漏声时域信号，计算互功率谱密度函数；对互功率谱密度函数加权后进行傅里叶逆变换得到互相关函数，获取互相关函数峰值对应的时间延迟估计，从而通过定位公式反推出泄漏点相对第一传感器的距离。本发明的方法无需对传感器信号进行滤波处理；能够有效消除反射影响，从而更准确实现复杂工况下管道泄漏定位。

Description

一种基于互相关的管道泄漏定位方法
技术领域
本发明涉及管道输送领域，特别是涉及到复杂工况下城市供水管道的泄漏定位方法，具体涉及一种基于互相关的管道泄漏定位方法。
背景技术
随着城市化进程的加速和人口的递增，水资源的需求量不断增加。我国城市供水管网规划建设、运营管理中，普遍存在地下基础设施腐蚀、老化，管线信息缺乏，过度的交通负载和人为破坏，监测及运营管理手段落后等问题，导致管道泄漏频发。管道泄漏通常占水损失总量的70％以上，造成直接和间接的能源、经济损失和安全隐患；二则，当前综合智能城市管廊建设迫切需要管道泄漏检测与定位方法的创新和技术转化。因此，针对地下供水管网具有多分支、多节点的复杂工况，一种及时、准确的管道泄漏定位方法尤为重要。目前通常使用声学法进行城市供水管道泄漏检测，特别是采用相关检漏法实现漏点定位。相关检漏法采用互相关分析法计算泄漏声沿管道传播时到达安装在管道上的两端的传感器拾取信号之间的时延，需要已知声波传播速度以及两个测点的管道距离来确定泄漏点的位置。但该方法在实际管道漏损定位时有三个问题是难以解决的：一是当检测非金属管道、水压低、深埋管等工况时，则对漏点定位精度不高，甚至误判；二是通常需预知泄漏声和背景噪声的先验知识，选取泄漏声占主导的非频散频段(文献[1]：Y.Gao,M.J.Brennan,P.F.Joseph,J.M.Muggleton,O.Hunaidi,A model of the correlation function ofleak noise in buried plastic pipes,Journal of Sound and Vibration 277(1-2)(2004)133-148.)，先对传感器拾取信号进行滤波处理，然实际中很难获得这些先验知识；三是当泄漏声在管道的分支、节点处产生声波的反射(文献[2]：Y.Gao,M.J.Brennan,P.F.Joseph,On the effect of reflections on time delay estimation for leakdetection in buried plastic water pipes,Journal of Sound and Vibration 325(3)(2009)649-663.)，将对漏点定位不准确。
发明内容
本发明的目的在于克服目前管道泄漏定位方法存在的上述问题，提出了一种基于互相关的管道泄漏定位方法。
为了实现上述目的，本发明提出了一种基于互相关的管道泄漏定位方法，基于下述检测装置实现：在待检测管道上安装两个传感器，传感器为加速度传感器或者水听器；每个传感器各连接数据采集仪，用于同步采集管道泄漏时引起的流体中传播的泄漏声信号；所述方法包括：获取两个传感器拾取的泄漏声时域信号，计算互功率谱密度函数；对互功率谱密度函数加权后进行傅里叶逆变换得到互相关函数，获取互相关函数峰值对应的时间延迟估计，从而通过定位公式反推出泄漏点相对第一传感器的距离。
作为上述方法的一种改进，所述方法具体包括：
步骤1)获取两个传感器拾取的泄漏声压时域信号xk(t)，k＝1,2；
步骤2)在频域上建立两组时域信号之间的互功率谱密度函数：
Figure BDA0001561769380000021
其中，X1(ω)、X2(ω)分别表示x1(t)、x2(t)的傅里叶变换，*表示复共轭；
步骤3)对
Figure BDA0001561769380000022
乘以频率加权函数ω2之后进行傅里叶逆变换，得到互相关函数：
Figure BDA0001561769380000023
其中，
Figure BDA0001561769380000024
表示傅里叶逆变换；
步骤4)获取互相关函数
Figure BDA0001561769380000025
峰值对应的时间延迟估计ΔT；
步骤5)通过定位公式求出泄漏点相对第一传感器的距离d1：
d1＝(D-cΔT)/2
其中，D为两个传感器之间的距离，c为泄漏声传播速度。
本发明的优势在于：
1、本发明提出的方法直接对传感器拾取的时域信号经行傅里叶变换，通过在频域上对互谱密度函数加权后作傅里叶逆变换求取互相关结果；无需掌握泄漏声信号和背景噪声的频域信息，简单易实现，对漏损检测人员的相关专业知识要求低；
2、本发明的方法无需对传感器信号进行滤波处理；能够有效消除反射影响，从而更准确实现复杂工况下管道泄漏定位。
附图说明
图1为本发明的原理示意图；
图2为本发明的方法的流程图；
图3为n＝0时，互相关函数的示意图；
图4为n＝2时，互相关函数的示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细的说明。
如图1所示，本发明的原理为：
在时域流程图中，xk(t),(k＝1,2),为两个传感器拾取的时域信号，对其进行预滤波后对延迟信号x2(t+τ)和x1(t)相乘并作积分，求得互相关函数
Figure BDA0001561769380000031
其中
Figure BDA0001561769380000032
Figure BDA0001561769380000033
表示傅里叶逆变换；时域原理也可表示为
Figure BDA0001561769380000034
其中
Figure BDA0001561769380000035
表示卷积，RBCC(τ)为基本互相关函数。
该处理过程可视为，在频域上先求两个传感器拾取信号的互功率谱密度(CSD)，再对互谱
Figure BDA0001561769380000036
进行加权后作逆傅里叶变换(IFT)得到互相关函数
Figure BDA0001561769380000037
由于对传感器拾取的时域信号直接求n阶导数运算，
Figure BDA0001561769380000038
因此加权函数ωn，可视为当忽略其相变in时的求导运算所产生的系数，其中
Figure BDA0001561769380000039
表示傅里叶变换。
因为供水管道泄漏声信号具有低频、窄带特性，设定管道中泄漏声主要集中在f0内，特别是在小于几Hz和大于f0的频段上，传感器拾取的信号通常为背景噪声占主导。当对互谱密度函数加权后，使得几Hz内的低频信号相对较高频段上的贡献量大幅降低。因此直接消除了低频背景噪声的干扰，而不用人为进行滤波选取去除背景噪声的操作。另外，通常管道中声波反射对拾取信号在低频段的影响大，可见通过对互谱密度函数进行加权后，可同时有效去除反射影响。
如图2所示，本发明方法基于互相关的原理，提出一种用于管道泄漏的定位方法，包括以下步骤：
步骤一：在待测管道上安装两个传感器实现同步采集，两个传感器各拾取一组时域信号；
所述同步采集为数据采集仪在同一时刻采样；
所述传感器为声信号-电信号传感器；
所述传感器感知管道泄漏时引起的流体中传播的泄漏声压信号xk(t),(k＝1,2)；
步骤二：在频域上，建立两组传感器拾取信号之间的互功率谱密度函数，乘以加权函数之后，对其进行傅里叶逆变换，得到互相关函数
Figure BDA0001561769380000041
其中加权函数为ωn，n＝0,2,…,N；。特别地，当n＝0,
Figure BDA0001561769380000042
给出基本互相关函数RBCC(τ)。
所述互相关函数可以用公式表示为
Figure BDA0001561769380000043
其中
Figure BDA0001561769380000044
表示傅里叶逆变换，*表示复共轭；
步骤三：求取互相关
Figure BDA0001561769380000045
时，一般地可设置n＝1-4，可使互相关函数的峰值得到锐化，从而有利于峰值点的提取。通常当n＝2时加权函数为最优，取峰值对应的时间延迟估计进行漏点定位。
实例：
步骤一：地下供水管道发生泄漏时，沿待测管线选取两处检修井，并在管道两端各安装一个传感器，传感器类型可以是加速度传感器或者水听器。传感器连接数据采集仪，同步采集两组时域信号。
步骤二：计算两组传感器拾取的时域信号xk(t),(k＝1,2)的互功率谱密度：
Figure BDA0001561769380000046
其中*表示复共轭。对互功率谱密度加权，其中加权函数设为ωn(n＝2)，对
Figure BDA0001561769380000047
乘以频率加权函数ω2之后，进行傅里叶逆变换，得到互相关函数：
Figure BDA0001561769380000048
步骤三：实验中可以通过定位公式反演计算出实际时间延迟，时延公式可以表示为：
Figure BDA0001561769380000049
式中，D为传感器之间的距离，单位：m；d1为第一个传感器与泄漏源之间的距离，单位：m；c为泄漏声传播速度，单位：m/s。本实施例中D＝102.6m，d1＝73.5m，c＝479m/s，求出精确时间延迟估计为τ＝-0.0927s。
采用本发明方法，根据步骤二中互相关函数的峰值求得时间延迟，求得的互相关函数如图3所示，时间延迟为τ1＝-0.09s，相对误差(|τ|-|τ1|)/τ等于2.9％，能够实现泄漏定位。本发明提出的方法中当n＝0时，即为基本互相关方法，结果如图4所示，时间延迟估计为0.326s，无法进行泄漏定位。
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (1)

1.一种基于互相关的管道泄漏定位方法，基于下述检测装置实现：在待检测管道上安装两个传感器，传感器为加速度传感器或者水听器；每个传感器各连接数据采集仪，用于同步采集管道泄漏时引起的流体中传播的泄漏声信号；所述方法包括：获取两个传感器拾取的泄漏声时域信号，计算互功率谱密度函数；对互功率谱密度函数加权后进行傅里叶逆变换得到互相关函数，获取互相关函数峰值对应的时间延迟估计，从而通过定位公式反推出泄漏点相对第一传感器的距离；
所述方法具体包括：
步骤1)获取两个传感器拾取的泄漏声压时域信号xk(t)，k＝1,2；
步骤2)在频域上建立两组时域信号之间的互功率谱密度函数：
Figure FDA0002368120470000011
其中，X1(ω)、X2(ω)分别表示x1(t)、x2(t)的傅里叶变换，*表示复共轭；
步骤3)对
Figure FDA0002368120470000012
乘以频率加权函数ω2之后进行傅里叶逆变换，得到互相关函数：
Figure FDA0002368120470000013
其中，
Figure FDA0002368120470000014
表示傅里叶逆变换；
步骤4)获取互相关函数
Figure FDA0002368120470000015
峰值对应的时间延迟估计ΔT；
步骤5)通过定位公式求出泄漏点相对第一传感器的距离d1：
d1＝(D-cΔT)/2
其中，D为两个传感器之间的距离，c为泄漏声传播速度。
CN201810083681.0A 2018-01-29 2018-01-29 一种基于互相关的管道泄漏定位方法 Active CN108332063B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810083681.0A CN108332063B (zh) 2018-01-29 2018-01-29 一种基于互相关的管道泄漏定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810083681.0A CN108332063B (zh) 2018-01-29 2018-01-29 一种基于互相关的管道泄漏定位方法

Publications (2)

Publication Number Publication Date
CN108332063A CN108332063A (zh) 2018-07-27
CN108332063B true CN108332063B (zh) 2020-04-24

Family

ID=62926023

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810083681.0A Active CN108332063B (zh) 2018-01-29 2018-01-29 一种基于互相关的管道泄漏定位方法

Country Status (1)

Country Link
CN (1) CN108332063B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109595472B (zh) * 2018-11-29 2020-07-24 北京市燃气集团有限责任公司 基于频谱分析的中低压庭院管线异常流量识别方法及装置
CN109538938B (zh) * 2018-12-05 2021-03-16 浙江海洋大学 一种用于荒漠区深埋输油管道检漏的“千足虫”仿生爬行设备
CN110501417B (zh) * 2019-09-20 2020-06-30 华北电力大学 一种基于声学的管道内蒸汽湿度在线监测系统及方法
CN111024351B (zh) * 2019-12-10 2021-10-22 西安理工大学 一种基于二维波束聚焦理论的冲击监测定位方法
CN110939874B (zh) * 2019-12-16 2021-08-31 重庆邮电大学 基于四阶累积量的管道泄漏振动信号自适应时延估计方法
CN111314010B (zh) * 2020-02-20 2021-03-02 山东大学 一种多通道延时可调的fx相关器及其实现方法
CN111350948A (zh) * 2020-04-13 2020-06-30 安徽理工大学 一种基于波束形成的管道泄漏位置计算方法
CN111720755B (zh) * 2020-04-15 2022-09-09 厦门矽创微电子科技有限公司 一种家用管道探漏定位方法和系统
CN113048404B (zh) * 2021-03-12 2022-08-16 常州大学 一种城市气体管道微小泄漏诊断方法
CN113607347B (zh) * 2021-07-23 2024-04-23 山东非金属材料研究所 面向真空热压成型复合装甲材料泄漏的检测方法及装置
CN114877262A (zh) * 2022-04-20 2022-08-09 成都千嘉科技股份有限公司 基于ai芯片燃气管道泄漏的声发射检测方法
CN115143396A (zh) * 2022-06-03 2022-10-04 北京建筑大学 一种基于能量谱的供热管道泄漏检测系统与方法
CN114719198B (zh) * 2022-06-09 2022-10-14 武汉易维科技股份有限公司 一种基于渗噪积分域建模的管道渗漏智能检测方法及装置
CN115856426A (zh) * 2022-11-23 2023-03-28 吉林大学 基于粒子群自适应归一化的随机共振微弱信号检测方法
CN115789533A (zh) * 2022-11-24 2023-03-14 新兴铸管股份有限公司 燃气管道泄漏定位方法、装置、电子设备及存储介质

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102411138A (zh) * 2011-07-13 2012-04-11 北京大学 一种机器人声源定位方法
CN102305925A (zh) * 2011-07-22 2012-01-04 北京大学 一种机器人连续声源定位方法
CN102854494B (zh) * 2012-08-08 2015-09-09 Tcl集团股份有限公司 一种声源定位方法及装置
KR101525329B1 (ko) * 2013-12-30 2015-06-03 한국원자력연구원 모드분리기법을 이용한 매설배관의 누설위치 추정방법
CN104747912B (zh) * 2015-04-23 2017-04-12 重庆邮电大学 流体输送管道泄漏声发射时频定位方法
CN105811921B (zh) * 2015-08-31 2017-06-23 维沃移动通信有限公司 一种抑制工频谐波干扰的方法及滤波器

Also Published As

Publication number Publication date
CN108332063A (zh) 2018-07-27

Similar Documents

Publication Publication Date Title
CN108332063B (zh) 一种基于互相关的管道泄漏定位方法
CN104595729B (zh) 一种基于声波幅值的油气管道泄漏定位方法
Li et al. Leak location in gas pipelines using cross-time–frequency spectrum of leakage-induced acoustic vibrations
CN104595730B (zh) 一种基于声波幅值衰减模型的油气管道泄漏定位方法
CN101246467B (zh) 自适应阈值检漏结合多尺度快速时延搜索的泄漏定位方法
CN104895556B (zh) 油井动液面远程监测方法及系统
CN102563361A (zh) 基于常规数据和音波信号的输气管道泄漏检测定位装置及方法
CN102537670B (zh) 一种管道泄漏诊断方法
CN102606891A (zh) 漏水检测仪及检测系统、及检测方法
CN101551064A (zh) 供水管道泄漏检测定位的信号处理方法
AU2017393649B2 (en) Leakage positioning method based on speed difference
CN106289121B (zh) 一种变径管等效管长的计算方法
CN110469782B (zh) 一种基于自适应滤波器的管道泄漏定位装置
CN105840987A (zh) 一种基于压力波和声波的管道泄漏加权定位方法及装置
CN105953080B (zh) 基于同侧传感器布置的声波泄漏定位方法
CN101592288A (zh) 一种管道泄漏辨识方法
CN108181059B (zh) 基于小波信号的多相流管道泄漏声波信号识别方法
CN111122068A (zh) 一种基于光纤水听器的输水管道渗漏位置定位方法
CN113375875A (zh) 给水管道慢性泄漏的声学检测系统及方法
KR101381469B1 (ko) 매설배관 누설 탐지용 상호상관함수기법의 정확도 향상을 위한 기계 잡음 제거 방법
CN118604735A (zh) 一种埋地pe燃气管道定位方法、装置和存储介质
CN106678553B (zh) 一种泄漏动态压力波在管内气体中传播速度的计算方法
KR20130064403A (ko) 매설배관 누설 탐지용 상호상관함수기법의 정확도 향상을 위한 기계 잡음 제거 방법
CN117743765A (zh) 基于声波的输水管道泄漏检测定位方法
CN106195648B (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