CN111365623B - Negative pressure wave inflection point identification method based on linear fitting - Google Patents

Negative pressure wave inflection point identification method based on linear fitting 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
Chinese (zh)
Other versions
CN111365623A (en
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/en
Publication of CN111365623A publication Critical patent/CN111365623A/en
Application granted granted Critical
Publication of CN111365623B publication Critical patent/CN111365623B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17DPIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00Protection or supervision of installations
    • F17D5/005Protection or supervision of installations of gas pipelines, e.g. alarm
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17DPIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00Protection or supervision of installations
    • F17D5/02Preventing, monitoring, or locating loss
    • F17D5/06Preventing, monitoring, or locating loss using electric or acoustic means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L11/00Measuring 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/02Measuring 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/025Measuring 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00Investigating fluid-tightness of structures
    • G01M3/02Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/26Investigating 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/28Investigating 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/2807Investigating 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/2815Investigating 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

The invention belongs to the technical field of pipeline leakage detection, and relates to a negative pressure wave inflection point identification method based on linear fitting. The method comprises the steps of respectively selecting appropriate data to perform linear fitting on a stationary section of a negative pressure wave denoising signal before leakage of a pipeline and a descending section of the negative pressure wave denoising signal after leakage of the pipeline; the inflection point of the negative pressure wave is rapidly and accurately obtained by solving the intersection point of the two straight lines. The method improves the accuracy and speed of identifying the inflection point of the negative pressure wave on the premise of keeping the signal characteristics of the original negative pressure wave as much as possible, has simple principle, can be directly embedded into most program systems, and has low cost and wide application range. The real-time online detection of the negative pressure wave inflection point and other pipeline leakage detection and positioning methods related to negative pressure wave inflection point identification have positive significance in the aspects of leakage point positioning accuracy and positioning speed.

Description

Negative pressure wave inflection point identification method based on linear fitting
Technical Field
The invention belongs to the technical field of pipeline leakage detection, and relates to a negative pressure wave inflection point identification method based on linear fitting.
Background
The pipeline is just like a life line and bears the important task of safely transporting petroleum, natural gas, water resources and the like and relating to national civilian resources. However, as the service life of the pipeline increases, pipeline leakage often occurs due to corrosion, aging, weld defects, geological condition changes, third party damage, and the like. Especially for long-distance high-pressure pipelines, once the pipelines leak, the leakage amount is large, the repair difficulty is high, the destructiveness and the harmfulness of the pipelines are far larger than those of low-level pipelines, and serious leakage accidents are caused. Therefore, the leakage phenomenon of the high-pressure pipeline is detected in time, the leakage position of the pipeline is accurately positioned, and the method has practical significance for safe operation of the pipeline, guarantee of life and property safety of people and reduction of environmental pollution.
At present, a plurality of methods for researching leakage and positioning of the pipeline exist, and the main methods comprise a pressure gradient method, a sound wave method, a vibration signal method and a negative pressure wave method. Among these methods, the negative pressure wave method has become one of the main means for detecting the leakage of the pipeline at present because of its characteristics of simple principle and sensor arrangement and relatively accurate measurement. The pipeline leakage detection and positioning method based on the negative pressure wave method mainly comprises the step of obtaining the arrival time difference of the negative pressure waves acquired by the sensor to position a leakage point, so that the arrival time of the negative pressure waves acquired by the sensor is required to be as accurate as possible.
The time of arrival of the negative pressure wave is determined by the time of occurrence of the inflection point of the negative pressure wave, however, the current methods for identifying the inflection point of the negative pressure wave are few and mainly focus on a threshold method and a singular point identification method. The threshold value method is too conservative, so that the negative pressure wave inflection point identification is delayed frequently, and the positioning error is large. The singular point identification method is used for identifying the singular point of the negative pressure wave signal instead of identifying the inflection point so as to determine the arrival time of the negative pressure wave. However, it requires multi-scale wavelet transformation of the negative pressure wave signal, which causes singular point deviation as the scale increases, thereby causing measurement errors; in addition, the workload is large, the calculation time is long, and the position of the pipeline leakage occurring cannot be located at the first time. Therefore, the method for accurately and rapidly identifying the inflection point of the negative pressure wave has important practical significance.
Disclosure of Invention
In order to solve the problems, the invention provides a negative pressure wave inflection point identification method based on linear fitting. The negative pressure wave signals are basically kept stable before leakage by analyzing a large amount of denoised pipeline leakage negative pressure wave signals; after the leakage occurs, the negative pressure wave signal is linearly decreased. It is like a sudden turning of a horizontal straight line, a new straight line is formed, and the intersection point of the two straight lines is the inflection point of the negative pressure wave. By solving the inflection point, the time of arrival of the negative pressure wave can be obtained more accurately.
The technical scheme of the invention is as follows:
a negative pressure wave inflection point identification method based on linear fitting comprises the following specific steps:
step 1: and de-noising the negative pressure wave signal generated by the pipeline leakage and acquired by the sensor.
Step 2: selecting P time points and corresponding data values of the denoising signal stationary section obtained by the processing in the step 1 for linear fitting to obtain a straight line L1=a1t+b1
And step 3: and averaging the data values corresponding to the P selected time points to set the data values as a leakage threshold S.
And 4, step 4: and comparing the data values of all points in all the data of the negative pressure wave signal with S in sequence from small to large according to the corresponding time. When the data values corresponding to the N continuous time points are lower than the set S, the pipeline leakage is considered to occur; extracting the 1 st time point of the N time points, and performing linear fitting on the next M time points and corresponding data values thereof to obtain a straight line L2=a2t+b2
And 5: solving two straight lines L1And L2The intersection point of the two points is the inflection point of the negative pressure wave, so that the time t (b) when the sensor detects the arrival of the negative pressure wave can be obtained2-b1)/(a1-a2)。
The invention has the advantages that: according to the negative pressure wave inflection point identification method based on linear fitting, on the premise that denoising processing is carried out only once and original negative pressure wave signals are not greatly changed, original signal characteristics are kept as far as possible, and the negative pressure wave inflection point identification accuracy is improved. The method has simple principle, can be directly embedded into most program systems, and has the advantages of high detection speed, low cost and wide application range. The real-time online detection of the negative pressure wave inflection point and other pipeline leakage detection and positioning methods related to negative pressure wave inflection point identification have positive significance in the aspects of leakage point positioning accuracy and positioning speed.
Drawings
Fig. 1 is a schematic diagram of a negative pressure wave inflection point identification method based on linear fitting.
Fig. 2 is a partially enlarged view of the inflection point position of the negative pressure wave in fig. 1.
In the figure: 1 denoising a signal; 2 straight line L1=a1t+b1(ii) a 3 straightLine L2=a2t+b2(ii) a 4 negative pressure wave inflection point.
Detailed Description
The following detailed description of the invention refers to the accompanying drawings.
The schematic diagram of the negative pressure wave inflection point identification method based on linear fitting is shown in fig. 1.
The specific process is as follows:
step 1: two fiber grating pressure sensors are respectively arranged at two ends of one gas transmission steel pipeline, and the fiber grating sensors are selected because of high sensitivity, so that the method has higher measurement precision. The data acquisition equipment is a fiber grating demodulator, and the sampling frequency is 1000 Hz. And (3) carrying out denoising processing on a negative pressure wave signal generated by the pipeline leakage and acquired by one of the sensors to obtain a denoising signal 1 shown in fig. 1. The wavelet transformation is used for the denoising processing of the negative pressure wave signal, the denoising method can obviously reduce the time-course noise of the negative pressure wave signal, the details of the original signal are still reserved, and the solving precision of the method is further improved.
Step 2: as can be seen from the denoised signal 1 in FIG. 1, the stationary part of the denoised signal 1 is about 6000 points, and the inflection point of the negative pressure wave appears between about 6000 and 8000 points. And selecting a numerical value close to the inflection point of the predicted negative pressure wave to perform linear fitting, so that the inflection point identification precision is higher. Thus, according to the method, the value of P is chosen to be 1000, specifically from the 5000 th point to the 6000 th point of the de-noised signal 1. Performing linear fitting on the 1000 time points and the corresponding data values to obtain a straight line L1=a1t+b 1 2。
And step 3: averaging the data values corresponding to the 1000 time points selected in step 2, and setting the average as a leakage threshold
Figure BDA0002404041740000041
And 4, step 4: for each point data value in all the data of the negative pressure wave signal, the data values are sequentially matched with the data values of all the points in the sequence from small to large according to the corresponding time
Figure BDA0002404041740000043
A comparison is made. When the data values corresponding to N continuous time points are lower than the set data values
Figure BDA0002404041740000042
The pipeline leakage is considered to occur; extracting the 1 st time point of the N time points, and performing linear fitting on the next M time points and corresponding data values thereof to obtain a straight line L2=a2t+b 2 3。
If the value N is set to be too small, the method is too sensitive, and pipeline vibration or fluid disturbance in the pipeline can be identified as leakage to generate false alarm; if the setting is too large, the calculation amount is increased, and the time for identifying the inflection point of the negative pressure wave is prolonged. The value of M is not suitable to be set to be too large, because the pressure is suddenly increased due to the upstream flow compensation after the negative pressure wave passes through, the measured value is increased to some extent, and the leakage section straight line L is formed2=a2t+b 23, certain fluctuation occurs, and the accuracy of linear fitting is influenced.
As can be seen from the de-noising signal 1 in fig. 1, the descending section of the de-noising signal 1 of the negative pressure wave has about 5000 points, and according to the method, the value of N is set to 1000; as can be seen from the denoised signal 1 in fig. 2, a certain rising fluctuation occurs in the measured value between 6900 th point and 7000 th point of the descent segment, so the value M should not be too large, and is set to 100. If in different pipeline models and practical engineering application, appropriate N and M values are set according to the characteristics of the pipeline leakage signal so as to ensure the measurement accuracy of the algorithm.
And 5: solving two straight lines L1And L2The intersection point of (a) is the inflection point 4 of the negative pressure wave, so that the time t (b) when the sensor detects the arrival of the negative pressure wave can be obtained2-b1)/(a1-a2). According to the method, the arrival time of the negative pressure wave is calculated to be 6.843 s. And (5) repeating the steps 1 to 5, and obtaining that the inflection point arrival time of the negative pressure wave acquired by the other sensor is 6.729 s. The arrangement distance of the two sensors is 48m, the leakage position of the pipeline is positioned by using a traditional negative pressure wave method, and the positioning error is only 0.5 m. To verify the reliability of the method, a total of 40 sets of tests were carried outThe average positioning error is 1.1m, which shows that the method can effectively improve the identification precision of the inflection point of the negative pressure wave.

Claims (1)

1. A negative pressure wave inflection point identification method based on linear fitting is characterized by comprising the following steps:
step 1: carrying out denoising processing on a negative pressure wave signal generated by pipeline leakage and acquired by a sensor;
step 2: selecting P time points and corresponding data values of the denoising signal stationary section obtained by the processing in the step 1 for linear fitting to obtain a straight line L1=a1t+b1
And step 3: averaging the data values corresponding to the P selected time points, and setting the data values as a leakage threshold S;
and 4, step 4: sequentially comparing data values of all points in all data of the negative pressure wave signal with S according to the sequence of the corresponding time from small to large; when the data values corresponding to the N continuous time points are lower than the set S, the pipeline leakage is considered to occur; extracting the 1 st time point of the N time points, and performing linear fitting on the next M time points and corresponding data values thereof to obtain a straight line L2=a2t+b2
And 5: solving two straight lines L1And L2The intersection point of the two points is the inflection point of the negative pressure wave, and the time t (b) when the sensor detects the arrival of the negative pressure wave is obtained2-b1)/(a1-a2)。
CN202010155928.2A 2020-03-09 2020-03-09 Negative pressure wave inflection point identification method based on linear fitting Active CN111365623B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010155928.2A CN111365623B (en) 2020-03-09 2020-03-09 Negative pressure wave inflection point identification method based on linear fitting

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010155928.2A CN111365623B (en) 2020-03-09 2020-03-09 Negative pressure wave inflection point identification method based on linear fitting

Publications (2)

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

Family

ID=71204334

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010155928.2A Active CN111365623B (en) 2020-03-09 2020-03-09 Negative pressure wave inflection point identification method based on linear fitting

Country Status (1)

Country Link
CN (1) CN111365623B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113188055B (en) * 2021-05-31 2022-11-25 东北大学 Pipeline leakage self-adaptive dynamic compensation positioning method based on negative pressure wave attenuation driving
CN113639208B (en) * 2021-07-16 2022-12-06 汕头大学 Negative pressure wave signal inflection point positioning method and system based on extreme approximation
CN113933002A (en) * 2021-08-12 2022-01-14 吉林大学 Method for identifying leakage point of long-distance large-scale water delivery pressure pipeline hydrostatic test

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105987284A (en) * 2015-01-27 2016-10-05 中国石油天然气股份有限公司 Pipeline leakage monitoring method and device
CN106481989A (en) * 2016-12-01 2017-03-08 山东省科学院激光研究所 pipeline leakage detection method and device
CN107084313A (en) * 2017-05-16 2017-08-22 云南大红山管道有限公司 Ore slurry pipeline leaks positioning alarm system and method
CN107620868A (en) * 2017-11-08 2018-01-23 山东省科学院激光研究所 Pipeline leakage detection method and device
CN108019622A (en) * 2018-02-05 2018-05-11 吉林大学 A kind of computational methods of the pipeline leakage positioning based on pressure differential
CN108443714A (en) * 2018-04-09 2018-08-24 山东省科学院激光研究所 Determine the method and device of pipe leakage
CN109114424A (en) * 2017-06-26 2019-01-01 中国石油天然气股份有限公司 Pipe leakage point method for determining position and device

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102853261A (en) * 2011-06-27 2013-01-02 国际商业机器公司 Method and device for determining leakage amount of fluid in conveying pipeline
CN103822097B (en) * 2012-11-16 2016-11-16 国际商业机器公司 The method and device of the velocity of wave of suction wave is estimated in fluid-transporting tubing

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105987284A (en) * 2015-01-27 2016-10-05 中国石油天然气股份有限公司 Pipeline leakage monitoring method and device
CN106481989A (en) * 2016-12-01 2017-03-08 山东省科学院激光研究所 pipeline leakage detection method and device
CN107084313A (en) * 2017-05-16 2017-08-22 云南大红山管道有限公司 Ore slurry pipeline leaks positioning alarm system and method
CN109114424A (en) * 2017-06-26 2019-01-01 中国石油天然气股份有限公司 Pipe leakage point method for determining position and device
CN107620868A (en) * 2017-11-08 2018-01-23 山东省科学院激光研究所 Pipeline leakage detection method and device
CN108019622A (en) * 2018-02-05 2018-05-11 吉林大学 A kind of computational methods of the pipeline leakage positioning based on pressure differential
CN108443714A (en) * 2018-04-09 2018-08-24 山东省科学院激光研究所 Determine the method and device of pipe leakage

Also Published As

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

Similar Documents

Publication Publication Date Title
CN111365623B (en) Negative pressure wave inflection point identification method based on linear fitting
CN104180166B (en) A kind of pipeline leakage detection method based on pipeline pressure data
Lang et al. A multiple leaks’ localization method in a pipeline based on change in the sound velocity
Aliev Noise control of the Beginning and Development Dynamics of Accidents
CN105840987B (en) It is a kind of that localization method and device are weighted based on the pipe leakage of pressure wave and sound wave
Hou et al. An improved negative pressure wave method for natural gas pipeline leak location using FBG based strain sensor and wavelet transform
CN109469834B (en) Liquid pipeline leakage detection method, device and system
CN103775832B (en) Based on the device that the petroleum pipeline leakage of transient flow Inverse Problem Method detects
US20220275610A1 (en) Method and system to monitor pipeline condition
WO2023088471A1 (en) Pipeline defect detection method based on multi-sensing information fusion and pipeline defect positioning method based on multi-sensing information fusion
CN114993229A (en) Multi-sensor data fusion-based ovality online detection method
CN1261744C (en) Method for detecting leakage of oil-gas pipeline based on pressure signals
CN115876288A (en) Electronic instrument fault analysis method and system based on big data
CN113626970A (en) Method and system for evaluating corrosion residual life of common pipe gallery pipeline
CN111828845A (en) Automatic pipeline leakage detection method based on artificial intelligence
CN1955636B (en) Testing process for the coupling zone between the cylindrical part and the clearance of a turbomachine part
Qu et al. Outlier Detection and Forecasting for Bridge Health Monitoring Based on Time Series Intervention Analysis.
CN115062891A (en) Method for quickly positioning leakage of long-distance pressure water delivery pipeline
CN202868318U (en) Pipe network tiny leak defect detecting system based on simulation technology
Zhang et al. Designing a reliable leak detection system for West Products Pipeline
CN108303173B (en) Distributed optical fiber sensing pipeline disturbance event detection method
CN112413413B (en) Pipeline leakage monitoring and positioning method combining deep learning and multiple measurement technology
CN114330442A (en) Pipeline strain characteristic classification calculation method and system based on K-nearest neighbor method
CN115935241B (en) Multi-parameter mutually-fused pipe cleaner real-time positioning method and device
CN112178360A (en) System and method for rapidly positioning infrasonic pipeline robot

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