[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN111365623B - 一种基于线性拟合的负压波拐点识别方法 - Google Patents

一种基于线性拟合的负压波拐点识别方法 Download PDF

Info

Publication number
CN111365623B
CN111365623B CN202010155928.2A CN202010155928A CN111365623B CN 111365623 B CN111365623 B CN 111365623B CN 202010155928 A CN202010155928 A CN 202010155928A CN 111365623 B CN111365623 B CN 111365623B
Authority
CN
China
Prior art keywords
negative pressure
pressure wave
inflection point
leakage
pipeline
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
CN202010155928.2A
Other languages
English (en)
Other versions
CN111365623A (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN202010155928.2A priority Critical patent/CN111365623B/zh
Publication of CN111365623A publication Critical patent/CN111365623A/zh
Application granted granted Critical
Publication of CN111365623B publication Critical patent/CN111365623B/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/005—Protection or supervision of installations of gas pipelines, e.g. alarm
    • 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
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L11/00—Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00
    • G01L11/02—Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00 by optical means
    • G01L11/025—Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00 by optical means using a pressure-sensitive optical fibre
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/26—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors
    • G01M3/28—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds
    • G01M3/2807—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds for pipes
    • G01M3/2815—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds for pipes using pressure measurements

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Examining Or Testing Airtightness (AREA)

Abstract

本发明属于管道泄漏检测技术领域，涉及到一种基于线性拟合的负压波拐点识别方法。该方法通过对管道发生泄漏前，负压波去噪信号的平稳段，以及管道发生泄漏后，负压波去噪信号的下降段，分别选取适当的数据进行线性拟合；通过求解两条直线的交点，以迅速且精确地获得负压波的拐点。本发明在尽可能保留了原始负压波信号特征的前提下，提高了负压波拐点识别的精度和速度，并且该方法原理简单，可直接嵌入绝大部分程序系统，成本低廉，适用范围广。对负压波拐点的实时在线检测，以及其他涉及到负压波拐点识别的管道泄漏检测与定位方法，在泄漏点定位精度以及定位速度方面，都将产生积极的意义。

Description

一种基于线性拟合的负压波拐点识别方法
技术领域
本发明属于管道泄漏检测技术领域，涉及到一种基于线性拟合的负压波拐点识别方法。
背景技术
管道就像是一条生命线，承担着安全运输石油、天然气、水资源等，关系国计民生资源的重大任务。但随着管道服役时间的增长，由于腐蚀、老化、焊缝缺陷、地质条件变化以及第三方破坏等原因，管道泄漏现象时常发生。特别对于长输高压管道来说，管道一旦发生泄漏，泄漏量大，修复难度高，其破坏性，危害性将远大于低级别的管道，导致严重的泄漏事故。因此，及时检测到高压管道的泄漏现象并准确地定位管道泄漏位置，对管道的安全运营以及保证人民生命财产安全，减少环境污染方面都具有切实意义。
目前研究管道泄漏与定位的方法有很多，主要方法包括压力梯度法、声波法、振动信号法和负压波法。在这些方法中，负压波法因其具有原理和传感器布置简单，测量相对准确的特点，已经成为目前管道泄漏检测的主要手段之一。基于负压波法的管道泄漏检测与定位方法，主要是获得传感器采集到负压波到达的时间差进行泄漏点的定位，这就要求传感器采集到负压波到达的时间要尽量精确。
负压波达到的时间是通过负压波拐点发生的时间确定的，然而，目前关于负压波拐点识别的方法还较少，主要集中在阈值法和奇异点识别方法。阈值法过于保守，经常导致负压波拐点识别延后，定位误差较大。奇异点识别法是通过对负压波信号奇异点的识别代替对拐点的识别，来确定负压波到达的时间。然而，它需要对负压波信号进行多尺度的小波变换，随着尺度的增加，这会导致奇异点偏离，从而导致测量误差；此外，工作量大，计算时间长，往往不能第一时间定位管道泄漏发生位置。因此，研发一种可以精确、迅速地识别负压波拐点的方法具有重要的现实意义。
发明内容
为解决上述问题，本发明提供一种基于线性拟合的负压波拐点识别方法。通过分析大量的去噪后的管道泄漏负压波信号发现，在泄漏之前，负压波信号是基本保持稳定的；发生泄漏之后，负压波信号呈线性降低。这就好比一条水平直线突然发生转折，形成一条新的直线，而两条直线的交点即为负压波的拐点。通过对拐点的求解，便可更加准确地获得负压波的达到时间。
本发明的技术方案是：
一种基于线性拟合的负压波拐点识别方法，具体步骤如下：
步骤1：对传感器采集到的管道泄漏所产生的负压波信号进行去噪处理。
步骤2：对步骤1处理得到的去噪信号平稳段，选取P个时刻点及其对应的数据值进行线性拟合，得到直线L1＝a1t+b1。
步骤3：对所选取的P个时刻点所对应的数据值求平均值，设为泄漏阈值S。
步骤4：对负压波信号全部数据中各点数据值，按其对应的时刻由小到大的顺序，依次与S进行比较。当有连续N个时刻点对应的数据值低于所设的S，则认为管道泄漏发生；提取这N个时刻点的第1个时刻点，并对其后M个时刻点及其对应的数据值进行线性拟合，得到直线L2＝a2t+b2。
步骤5：求解两条直线L1和L2的交点即为负压波拐点，便可获得传感器检测到负压波到达的时间t＝(b2-b1)/(a1-a2)。
本发明的效果和益处：本发明提出的一种基于线性拟合的负压波拐点识别方法，在仅进行一次去噪处理而未大幅改变原始负压波信号的前提下，尽可能的保留了原始信号特征，并提高了负压波拐点识别的精度。该方法原理简单，可直接嵌入绝大部分程序系统，检测速度快，成本低廉，适用范围广。对负压波拐点的实时在线检测，以及其他涉及到负压波拐点识别的管道泄漏检测与定位方法，在泄漏点定位精度以及定位速度方面，都将产生积极的意义。
附图说明
图1是一种基于线性拟合的负压波拐点识别方法示意图。
图2是图1中负压波拐点位置的局部放大图。
图中：1去噪信号；2直线L1＝a1t+b1；3直线L2＝a2t+b2；4负压波拐点。
具体实施方式
以下结合技术方案和附图详细叙述本发明的具体实施方式。
本发明提出的一种基于线性拟合的负压波拐点识别方法示意图如图1所示。
具体过程是：
步骤1：在一根输气钢管道两端各安装两个光纤光栅类压力传感器，之所以选择光纤光栅类传感器，是因为其灵敏度高，可使本方法测量精度更高。数据采集设备为光纤光栅解调仪，采样频率为1000Hz。对其中一个传感器采集到的管道泄漏所产生的负压波信号进行去噪处理，获得图1所示的去噪信号1。负压波信号去噪处理使用小波变换，该去噪方法可使负压波信号时程噪声明显减小，并且依然保留原始信号的细节，将进一步提高本方法的求解精度。
步骤2：由图1去噪信号1可以看出，去噪信号1的平稳段大约有6000个点，且负压波拐点出现的位置大约在第6000点到第8000点之间。选择离预计负压波拐点较近的数值，进行线性拟合，可使拐点识别精度更高。因此，根据本方法，P值选取为1000，具体是该去噪信号1的第5000点到第6000点。对这1000个时刻点及其对应的数据值进行线性拟合，得到直线L1＝a1t+b1 2。
步骤3：对步骤2中所选取的1000个时刻点所对应的数据值求平均值，设为泄漏阈值
Figure BDA0002404041740000041
步骤4：对负压波信号全部数据中各点数据值，按其对应的时刻由小到大的顺序，依次与
Figure BDA0002404041740000043
进行比较。当有连续N个时刻点对应的数据值低于所设的
Figure BDA0002404041740000042
则认为管道泄漏发生；提取这N个时刻点的第1个时刻点，并对其后M个时刻点及其对应的数据值进行线性拟合，得到直线L2＝a2t+b2 3。
N值若设置过小，会导致本发明过于灵敏，会将管道振动或管内流体扰动识别为泄漏，产生误报；若设定过大，也会增加计算量，延长负压波拐点识别的时间。M值不宜设置过大，因为负压波通过后，上游流量补入导致压力骤升，会使测量值有所升高，使得泄漏段直线L2＝a2t+b2 3出现一定波动，影响其线性拟合的精度。
由图1去噪信号1可以看出，负压波去噪信号1的下降段大约有5000个点，根据本方法，N值设为1000；由图2去噪信号1可以看出，在下降段第6900点到第7000点之间测量值出现一定回升波动，因此M值不宜过大，设为100即可。若在不同管道模型以及实际工程应用中，也应先根据管道泄漏信号的特征，设定合适的N和M值，以保证该算法的测量精度。
步骤5：求解两条直线L1和L2的交点即为负压波拐点4，便可获得传感器检测到负压波到达的时间t＝(b2-b1)/(a1-a2)。按本方法，计算得到负压波到达时间为6.843s。重复步骤1至步骤5，获得另外一个传感器所采集得到的负压波拐点到达时间为6.729s。两个传感器布置间距为48m，利用传统负压波法对管道泄漏位置进行定位，定位误差仅为0.5m。为验证本方法的可靠性，共进行40组测试，平均定位误差为1.1m，表明本方法可有效提高负压波拐点的识别精度。

Claims (1)

1.一种基于线性拟合的负压波拐点识别方法，其特征在于，步骤如下：
步骤1：对传感器采集到的管道泄漏所产生的负压波信号进行去噪处理；
步骤2：对步骤1处理得到的去噪信号平稳段，选取P个时刻点及其对应的数据值进行线性拟合，得到直线L1＝a1t+b1；
步骤3：对所选取的P个时刻点所对应的数据值求平均值，设为泄漏阈值S；
步骤4：对负压波信号全部数据中各点数据值，按其对应的时刻由小到大的顺序，依次与S进行比较；当有连续N个时刻点对应的数据值低于所设的S，则认为管道泄漏发生；提取这N个时刻点的第1个时刻点，并对其后M个时刻点及其对应的数据值进行线性拟合，得到直线L2＝a2t+b2；
步骤5：求解两条直线L1和L2的交点即为负压波拐点，便获得传感器检测到负压波到达的时间t＝(b2-b1)/(a1-a2)。
CN202010155928.2A 2020-03-09 2020-03-09 一种基于线性拟合的负压波拐点识别方法 Active CN111365623B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010155928.2A CN111365623B (zh) 2020-03-09 2020-03-09 一种基于线性拟合的负压波拐点识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010155928.2A CN111365623B (zh) 2020-03-09 2020-03-09 一种基于线性拟合的负压波拐点识别方法

Publications (2)

Publication Number Publication Date
CN111365623A CN111365623A (zh) 2020-07-03
CN111365623B true CN111365623B (zh) 2021-01-19

Family

ID=71204334

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010155928.2A Active CN111365623B (zh) 2020-03-09 2020-03-09 一种基于线性拟合的负压波拐点识别方法

Country Status (1)

Country Link
CN (1) CN111365623B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113188055B (zh) * 2021-05-31 2022-11-25 东北大学 基于负压波衰减驱动的管道泄漏自适应动态补偿定位方法
CN113639208B (zh) * 2021-07-16 2022-12-06 汕头大学 一种基于极限逼近的负压波信号拐点定位方法及系统
CN113933002B (zh) * 2021-08-12 2024-07-09 吉林大学 一种长距离大型输水压力管道水压试验漏点识别的方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105987284A (zh) * 2015-01-27 2016-10-05 中国石油天然气股份有限公司 管道泄漏监测方法及装置
CN106481989A (zh) * 2016-12-01 2017-03-08 山东省科学院激光研究所 管道泄漏检测方法及装置
CN107084313A (zh) * 2017-05-16 2017-08-22 云南大红山管道有限公司 矿浆管道泄漏定位报警系统及方法
CN107620868A (zh) * 2017-11-08 2018-01-23 山东省科学院激光研究所 管道泄漏检测方法及装置
CN108019622A (zh) * 2018-02-05 2018-05-11 吉林大学 一种基于压力差的管道泄露定位的计算方法
CN108443714A (zh) * 2018-04-09 2018-08-24 山东省科学院激光研究所 确定管道泄漏的方法及装置
CN109114424A (zh) * 2017-06-26 2019-01-01 中国石油天然气股份有限公司 管道泄漏点位置的确定方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102853261A (zh) * 2011-06-27 2013-01-02 国际商业机器公司 确定输送管道中的流体泄漏量的方法和装置
CN103822097B (zh) * 2012-11-16 2016-11-16 国际商业机器公司 在流体输送管道中估计负压波的波速的方法及装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105987284A (zh) * 2015-01-27 2016-10-05 中国石油天然气股份有限公司 管道泄漏监测方法及装置
CN106481989A (zh) * 2016-12-01 2017-03-08 山东省科学院激光研究所 管道泄漏检测方法及装置
CN107084313A (zh) * 2017-05-16 2017-08-22 云南大红山管道有限公司 矿浆管道泄漏定位报警系统及方法
CN109114424A (zh) * 2017-06-26 2019-01-01 中国石油天然气股份有限公司 管道泄漏点位置的确定方法及装置
CN107620868A (zh) * 2017-11-08 2018-01-23 山东省科学院激光研究所 管道泄漏检测方法及装置
CN108019622A (zh) * 2018-02-05 2018-05-11 吉林大学 一种基于压力差的管道泄露定位的计算方法
CN108443714A (zh) * 2018-04-09 2018-08-24 山东省科学院激光研究所 确定管道泄漏的方法及装置

Also Published As

Publication number Publication date
CN111365623A (zh) 2020-07-03

Similar Documents

Publication Publication Date Title
CN111365623B (zh) 一种基于线性拟合的负压波拐点识别方法
Lang et al. A multiple leaks’ localization method in a pipeline based on change in the sound velocity
CN108050396B (zh) 一种流体管道泄漏源监测定位系统及方法
CN104180166B (zh) 一种基于管道压力数据的管道泄漏检测方法
CN104063588B (zh) 基于多源数据融合的管道腐蚀缺陷尺寸的预测方法
CN107620868B (zh) 管道泄漏检测方法及装置
CN109469834B (zh) 一种液体管道泄漏检测方法、装置及系统
CN103775832B (zh) 基于瞬变流反问题方法的输油管道漏失检测的装置
WO2015072130A1 (ja) 漏洩判定システムおよび漏洩判定方法
WO2021022315A1 (en) Method and system to monitor pipeline condition
WO2023088471A1 (zh) 基于多传感信息融合的管道缺陷检测及定位方法
WO2024148884A1 (zh) 一种多传感器参数相互融合的清管器实时定位方法及装置
CN114993229A (zh) 一种基于多传感器数据融合的椭圆度在线检测方法
CN111828845A (zh) 基于人工智能的管道自动检漏方法
CN103853144A (zh) 基于采油生产数据的现场传感器故障检测方法
CN114088308B (zh) 一种基于低反射啁啾光栅阵列的输运管道拾振检漏方法
CN1584531A (zh) 基于压力信号检测油气管道泄漏的方法
Qu et al. Outlier Detection and Forecasting for Bridge Health Monitoring Based on Time Series Intervention Analysis.
CN1955636B (zh) 检查涡轮机组部件的圆柱部和梢部之间的连接区域的方法
CN117168713B (zh) 一种氦氢质谱检漏仪信号快速处理系统及方法
CN107330264A (zh) 一种桥梁监测数据可靠性的验证方法
CN108303173B (zh) 一种分布式光纤传感管道扰动事件检测方法
CN115062891A (zh) 一种长距离有压输水管道渗漏快速定位方法
CN112413413B (zh) 结合深度学习和多次测量技术的管道泄漏监测及定位方法
CN114330442A (zh) 一种基于k-近邻法的管道应变特征分类计算方法及系统

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