CN108662442A - The localization method and device of pipe leakage - Google Patents
The localization method and device of pipe leakage Download PDFInfo
- Publication number
- CN108662442A CN108662442A CN201810478059.XA CN201810478059A CN108662442A CN 108662442 A CN108662442 A CN 108662442A CN 201810478059 A CN201810478059 A CN 201810478059A CN 108662442 A CN108662442 A CN 108662442A
- Authority
- CN
- China
- Prior art keywords
- velocity
- leakage
- pipeline
- information
- signal
- 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
Links
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
The invention discloses a kind of localization method of pipe leakage and devices, it is related to pipeline inspection technology field, main purpose is existing unintelligible for Small leak and the pressure signal failing edge inflection point simmered, it is larger to cause to fail to report higher and position error, and the problem of installation pressure sensor can bring danger is punched using pipeline.Including：Obtain the ultrasonic velocity sensor for being separately mounted to pipeline first and last end, flow velocity signal, the velocity of sound signal of flow sensor acquisition；The leakage transient information at first and last end is calculated according to the flow velocity signal, the velocity of sound signal according to the leakage transient information propagation attenuation computation model pre-established；Head end information fusion signal and final word fusion signal after being merged into row information according to the leakage transient information to the first and last end determine delay time, and calculate distance of the leakage point apart from pipeline head end ultrasonic sensor using preset leak point positioning formula.
Description
Technical field
The present invention relates to a kind of pipeline inspection technology fields, more particularly to the localization method and dress of a kind of pipe leakage
It sets.
Background technology
Pipeline transportation is one of the important way of liquid resource long distance transmission, during pipeline transportation, the leakage of pipeline
Situation, which then becomes, to be influenced to ensure an important factor for liquid resource stablizes transmission.Wherein, in order to when pipe leakage occurs, can and
When carry out remedy, the important topic of pipe leakage being positioned to for research.
Although currently, existing negative pressure wave method carry out pipe leakage positioning can detect rapidly to be more than 3% it is normal defeated
The leakage of amount, and leak point positioning is carried out using its time difference according to the inflection point time of first and last side pressure force signal failing edge, but
It is that unintelligible for Small leak and the pressure signal failing edge inflection point that simmers, it is larger to cause to fail to report higher and position error, and makes
Danger can be brought by punching installation pressure sensor with pipeline, influence the locating effect of pipe leakage.
Invention content
In view of this, the present invention provides a kind of localization method and device of pipe leakage, main purpose be it is existing for
Small leak and the pressure signal failing edge inflection point simmered are unintelligible, and it is larger to cause to fail to report higher and position error, and uses pipeline
The problem of punching installation pressure sensor can bring danger, influence the locating effect of pipe leakage.
According to the present invention on one side, a kind of localization method of pipe leakage is provided, including：
Obtain be separately mounted to pipeline first and last end ultrasonic velocity sensor, flow sensor acquisition flow velocity signal,
Velocity of sound signal, the flow velocity signal include the flow velocity signal of pipeline head end and the flow velocity signal of pipe end, the velocity of sound signal
The velocity of sound signal of velocity of sound signal and pipe end including pipeline head end；
Believed according to the flow velocity signal, the velocity of sound according to the leakage transient information propagation attenuation computation model pre-established
The leakage transient information at first and last end number is calculated, the leakage transient information propagation attenuation computation model is the movement according to pipeline
The transient information computation model at the pipeline first and last end that equation and continuity equation combination ultrasonic velocity attenuation characteristic are determined；
Head end information fusion signal after being merged into row information according to the leakage transient information to the first and last end and end
Client information fusion signal determines delay time, and calculates leakage point apart from pipeline head end using preset leak point positioning formula
The distance of ultrasonic sensor.
Further, ultrasonic velocity sensor, the flow sensor that the acquisition is separately mounted to pipeline first and last end are adopted
Before the flow velocity signal of collection, velocity of sound signal, the method further includes：
The Transient Flow in Pipeline equation that liquid is established according to the equation of motion, continuity equation isWherein, v is flow velocity, and ρ is fluid density, and D is pipe
Road internal diameter, alpha levels axle clamp angle, λ are darcy hydraulic simulation experiment, and x is the axial distance along pipeline, and t is time, awFor suction wave
Velocity of wave, g are acceleration of gravity, acFor the velocity of sound in liquid, a0For the velocity of sound in liquid under one atmospheric pressure of room temperature, K is proportionality coefficient；
When ignoring elevation variation, then the node i at the j moment established according to gsin α=0 is along C+And C-Finite difference side
Journey is respectively Wherein, aci,jIndicate the ultrasonic acoustic at i-node j moment
Speed；vi,jIndicate the flow velocity at i-node j moment；
When Sudden Leakage point is the i-node j moment, and assume that j-1 moment pipeline fluids are in steady state condition, in j+1 moment i+2
The parameter of node remains unchanged, and solving to obtain leakage transient information propagation attenuation computation model according to the finite difference equations is
Wherein, the boundary condition of leakage point meetsFmFor the information of pipe end
Transition, Δ FmFor the increment of the information transition of pipe end, acmAnd vmThe respectively velocity of sound and flow velocity of pipe end, LxFor leakage
For point to the distance of pipeline first section, Δ x is run of designing, Δ x=awΔ t, FsFor the information transition of pipeline head end, Δ FsFor pipe
The increment of the information transition of road head end, acsAnd vsThe respectively velocity of sound and flow velocity of pipeline head end, LxFor leakage point to pipeline first section
Distance, Fu-For the information transition before leakage point.
Further, the head end information after the basis merges the leakage transient information at the first and last end into row information is melted
Before conjunction signal and final word fusion signal determine delay time, the method further includes：
Using the correspondence when pipeline leaks between the head end velocity of sound, flow and the end velocity of sound, flow to described
The leakage transient information at first and last end is merged into row information, obtains the leakage transient information F of first and endsAnd FmRespectively Fs=acs-
Ka0ρawvs, Fm=acm+Ka0ρawvm。
Further, the head end information after the basis merges the leakage transient information at the first and last end into row information is melted
It closes signal and final word fusion signal determines delay time, and leakage point is calculated using preset leak point positioning formula
Distance apart from pipeline head end ultrasonic sensor includes：
The leakage transient information F of head end after being merged respectively to information using wavelet analysis methodsAnd the leakage wink of end
Become information FmMutation time point is extracted, the mutation time is clicked through according to cross-correlation analysis method and pre-set time window
The rough time to peak search of row determined for the first delay time, and according to cross-correlation analysis hair and pre-set time window pair
The mutation is to determine for the second delay time by carrying out the search of accurate time to peak, will first delay time and described the
The sum of two delay times are determined as delay time；
Utilize preset leak point positioning formulaLeakage point is calculated apart from pipeline head end supersonic sensing
The distance of device, wherein the tDFor delay time, L is duct length.
According to the present invention on one side, a kind of positioning device of pipe leakage is provided, including：
Acquiring unit is separately mounted to ultrasonic velocity sensor, the flow sensor at pipeline first and last end and adopts for obtaining
Flow velocity signal, the velocity of sound signal of collection, the flow velocity signal include the flow velocity signal of pipeline head end and the flow velocity signal of pipe end,
The velocity of sound signal includes the velocity of sound signal of pipeline head end and the velocity of sound signal of pipe end；
First computing unit, for according to the leakage transient information propagation attenuation computation model pre-established according to the stream
Fast signal, the velocity of sound signal calculate the leakage transient information at first and last end, and the leakage transient information propagation attenuation calculates mould
Type is the pipeline first and last end determined according to the equation of motion and continuity equation combination ultrasonic velocity attenuation characteristic of pipeline
Transient information computation model；
Second computing unit is believed for the head end after being merged into row information according to the leakage transient information to the first and last end
Breath fusion signal and final word fusion signal determine delay time, and are calculated and let out using preset leak point positioning formula
Distance of the leak source apart from pipeline head end ultrasonic sensor.
Further, described device further includes：
First establishing unit, the Transient Flow in Pipeline equation for establishing liquid according to the equation of motion, continuity equation areWherein, v is flow velocity, and ρ is fluid density, and D is pipe
Road internal diameter, alpha levels axle clamp angle, λ are darcy hydraulic simulation experiment, and x is the axial distance along pipeline, and t is time, awFor suction wave
Velocity of wave, g are acceleration of gravity, acFor the velocity of sound in liquid, a0For the velocity of sound in liquid under one atmospheric pressure of room temperature, K is proportionality coefficient；
Second establishes unit, for when ignoring elevation variation, then the node i at the j moment established according to gsin α=0 to be along C+And C-Finite difference equations be respectively Wherein, aci,jIndicate the ultrasonic acoustic at i-node j moment
Speed；vi,jIndicate the flow velocity at i-node j moment；
Unit is solved, for being the i-node j moment when Sudden Leakage point, and assumes that j-1 moment pipeline fluids are in stable state
Operating mode is remained unchanged in the parameter of j+1 moment i+2 nodes, is solved to obtain leakage transient information according to the finite difference equations
Propagation attenuation computation model is Wherein, the boundary condition of leakage point meetsFmFor the information transition of pipe end, Δ FmFor the information of pipe end
The increment of transition, acmAnd vmThe respectively velocity of sound and flow velocity of pipe end, LxFor leakage point to the distance of pipeline first section, Δ x is
Run of designing, Δ x=awΔ t, FsFor the information transition of pipeline head end, Δ FsFor the increment of the information transition of pipeline head end, acs
And vsThe respectively velocity of sound and flow velocity of pipeline head end, LxFor leakage point to the distance of pipeline first section, Fu-For the information before leakage point
Transition.
Further, described device further includes：
Integrated unit, for utilizing pair when pipeline leaks between the head end velocity of sound, flow and the end velocity of sound, flow
It should be related to and the leakage transient information at the first and last end is merged into row information, obtain the leakage transient information F of first and endsAnd FmPoint
It Wei not Fs=acs-Ka0ρawvs, Fm=acm+Ka0ρawvm。
Further, second computing unit is specifically used for the head after being merged respectively to information using wavelet analysis method
The leakage transient information F at endsAnd the leakage transient information F of endmExtract mutation time point, according to cross-correlation analysis method and
Pre-set time window searches for the rough time to peak of mutation time point progress and determined for the first delay time, and according to
Cross-correlation analysis is sent out and pre-set time window is to determine second by the accurate time to peak search of progress to the mutation
The sum of first delay time and second delay time are determined as delay time by delay time；Utilize preset leakage
Point location formulaCalculate distance of the leakage point apart from pipeline head end ultrasonic sensor, wherein the tD
For delay time, L is duct length.
According to another aspect of the invention, a kind of storage medium is provided, at least one is stored in the storage medium can
It executes instruction, the executable instruction makes processor execute the corresponding operation of localization method such as above-mentioned pipe leakage.
In accordance with a further aspect of the present invention, a kind of terminal is provided, including：Processor, memory, communication interface and communication
Bus, the processor, the memory and the communication interface complete mutual communication by the communication bus；
For the memory for storing an at least executable instruction, it is above-mentioned that the executable instruction makes the processor execute
The corresponding operation of localization method of pipe leakage.
By above-mentioned technical proposal, technical solution provided in an embodiment of the present invention at least has following advantages：
The present invention provides a kind of localization method of pipe leakage and devices, obtain be separately mounted to pipeline first and last end first
Ultrasonic velocity sensor, flow sensor acquisition flow velocity signal, velocity of sound signal, the flow velocity signal includes pipeline head end
Flow velocity signal and pipe end flow velocity signal, the velocity of sound signal includes the velocity of sound signal and pipe end of pipeline head end
Velocity of sound signal；According to the leakage transient information propagation attenuation computation model pre-established according to the flow velocity signal, the velocity of sound
Signal calculates the leakage transient information at first and last end, and the leakage transient information propagation attenuation computation model is the fortune according to pipeline
The transient information at the pipeline first and last end that dynamic equation and continuity equation combination ultrasonic velocity attenuation characteristic are determined calculates mould
Type；Head end information fusion signal after being merged into row information according to the leakage transient information to the first and last end and final word
Fusion signal determines delay time, and calculates leakage point apart from pipeline head end ultrasonic wave using preset leak point positioning formula
The distance of sensor.With it is existing unintelligible for Small leak and the pressure signal failing edge inflection point simmered, cause to fail to report it is higher and
Position error is larger, and can bring danger using pipeline punching installation pressure sensor, influences the positioning effect of pipe leakage
Fruit is compared, and the embodiment of the present invention establishes velocity of sound generation and propagation attenuation computation model by the mechanism of transmission of ultrasonic velocity,
And ultrasonic velocity and flow signal being subjected to signal fused so that information fusion signal can more burst signal knee,
The time difference that pipeline first and last client information merges signal can be more accurately extracted by improving cross-correlation analysis method, so as to carry
The positioning accuracy of high Small leak fast and effeciently positions Small leak point, improves the efficiency of pipe leakage positioning.
Above description is only the general introduction of technical solution of the present invention, in order to better understand the technical means of the present invention,
And can be implemented in accordance with the contents of the specification, and in order to allow above and other objects of the present invention, feature and advantage can
It is clearer and more comprehensible, below the special specific implementation mode for lifting the present invention.
Description of the drawings
By reading the detailed description of hereafter preferred embodiment, various other advantages and benefit are common for this field
Technical staff will become clear.Attached drawing only for the purpose of illustrating preferred embodiments, and is not considered as to the present invention
Limitation.And throughout the drawings, the same reference numbers will be used to refer to the same parts.In the accompanying drawings：
Fig. 1 shows a kind of localization method flow chart of pipe leakage provided in an embodiment of the present invention；
Fig. 2 shows the generation of the velocity of sound at a kind of leakage point provided in an embodiment of the present invention and propagate schematic diagram；
Fig. 3 shows a kind of characteristic curve network schematic diagram provided in an embodiment of the present invention；
Fig. 4 shows the schematic diagram of pipeline first and last segment signal catastrophe point time difference provided in an embodiment of the present invention a kind of；
Fig. 5 shows a kind of pipeline first and last end provided in an embodiment of the present invention flow changing curve schematic diagram；
Fig. 6 shows a kind of pipeline first and last end ultrasonic velocity change curve schematic diagram provided in an embodiment of the present invention；
Fig. 7 shows a kind of pipeline first and last end pressure change curve schematic diagram provided in an embodiment of the present invention；
Fig. 8 shows a kind of pipeline head end information fusion signal schematic representation provided in an embodiment of the present invention；
Fig. 9 shows a kind of pipe end information fusion signal schematic representation provided in an embodiment of the present invention；
Figure 10 shows a kind of first time cross-correlation analysis schematic diagram provided in an embodiment of the present invention；
Figure 11 shows a kind of second of cross-correlation analysis schematic diagram provided in an embodiment of the present invention；
Figure 12 shows a kind of cross-correlation analysis schematic diagram of pipeline first and last end signal provided in an embodiment of the present invention；
Figure 13 shows under a kind of different leakage rates provided in an embodiment of the present invention that pipeline head end information merges signal contrast
Curve synoptic diagram；
Figure 14 shows under a kind of different leakage rates provided in an embodiment of the present invention that pipe end information merges signal contrast
Curve synoptic diagram；
Figure 15 shows a kind of positioning device structure schematic diagram of pipe leakage provided in an embodiment of the present invention；
Figure 16 shows the positioning device structure schematic diagram of another pipe leakage provided in an embodiment of the present invention；
Figure 17 shows a kind of terminal structure schematic diagrames provided in an embodiment of the present invention.
Specific implementation mode
The exemplary embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although showing the disclosure in attached drawing
Exemplary embodiment, it being understood, however, that may be realized in various forms the disclosure without should be by embodiments set forth here
It is limited.On the contrary, these embodiments are provided to facilitate a more thoroughly understanding of the present invention, and can be by the scope of the present disclosure
Completely it is communicated to those skilled in the art.
An embodiment of the present invention provides a kind of localization methods of pipe leakage, as shown in Figure 1, the method includes：
101, the ultrasonic velocity sensor for being separately mounted to pipeline first and last end, the flow velocity letter of flow sensor acquisition are obtained
Number, velocity of sound signal.
Wherein, the flow velocity signal includes the flow velocity signal of pipeline head end and the flow velocity signal of pipe end, the velocity of sound
Signal includes the velocity of sound signal of pipeline head end and the velocity of sound signal of pipe end.It is small in order to accurately determine in the embodiment of the present invention
The leakage point of leakage rate is separately installed with ultrasonic velocity sensor, flow sensing in head end and the end of pipeline to be detected
Device, to acquire flow velocity signal, the velocity of sound signal of pipeline head end and end respectively.
It should be noted that in the present embodiment, the phenomenon that leakage point is burst, in order at the time of accurate description leakage point and
Node situation, if i-node j moment Sudden Leakages, and assume that j-1 moment pipeline fluids are in steady state condition, the velocity of sound at leakage point
Generation and propagate schematic diagram, as shown in Figure 2.
In addition, the data of acquisition are by from gathered data NI data collecting systems DAQ-9184 on, the selection of I/O fasteners
It is the NI-9208 Dynamic Data Acquiring fasteners of NI companies, 16 Channel Synchronous acquire, and sample rate is up to 500S/s, and resolution ratio is
24bit.Data sampling rate is 100HZ, and all experiments are in Intel Pentium processor (2.90GHz) and 6GB RAM
PC machine on use 2015 software of MATLAB R2014a and Labview.
102, according to the leakage transient information propagation attenuation computation model pre-established according to the flow velocity signal, the sound
Fast signal calculates the leakage transient information at first and last end.
Wherein, the leakage transient information propagation attenuation computation model is according to the equation of motion of pipeline and continuity side
The transient information computation model at the pipeline first and last end that journey combination ultrasonic velocity attenuation characteristic is determined, the leakage at shown first and last end
Transient information is generated by the transition of the velocity of sound and the transition of flow, small due to occurring in pipeline in the embodiment of the present invention
When leakage rate, the head end velocity of sound is reduced, and flow increases and the end velocity of sound is reduced, and flow is reduced, therefore can be by measuring pipeline
The leakage transient information at first and last end monitors whether pipeline leaks, and accurately embodies leakage information, small to be accurately positioned
The anchor point for leaking and simmering.
103, according to the first and last end leakage transient information into row information merge after head end information fusion signal with
And final word fusion signal determines delay time, and leakage point is calculated apart from pipeline using preset leak point positioning formula
The distance of head end ultrasonic sensor.
For the embodiment of the present invention, due to including velocity of sound signal and flow signal in the transient information of first and last end, need
The leakage transient information at first and last end is merged into row information, signal and end are then merged by the head end information after fusion
Information fusion signal determines that delay time, the delay time are the sum of the time catastrophe point of first terminal fusion signal, so
Distance of the leakage point apart from pipeline head end ultrasonic sensor is calculated according to preset leak point positioning formula afterwards.
For the embodiment of the present invention, in order to pre-establish out leakage transient information propagation attenuation computation model, to obtain
The data direct modeling operation of acquisition is got, it is described to obtain the ultrasonic velocity sensor for being separately mounted to pipeline first and last end, stream
The flow velocity signal of quantity sensor acquisition, before velocity of sound signal, the method further includes：
The Transient Flow in Pipeline equation that liquid is established according to the equation of motion, continuity equation isWherein, v is flow velocity, and ρ is fluid density, and D is pipe
Road internal diameter, alpha levels axle clamp angle, λ are darcy hydraulic simulation experiment, and x is the axial distance along pipeline, and t is time, awFor suction wave
Velocity of wave, g are acceleration of gravity, acFor the velocity of sound in liquid, a0For the velocity of sound in liquid under one atmospheric pressure of room temperature, K is proportionality coefficient；
When ignoring elevation variation, then the node i at the j moment established according to gsin α=0 is along C+And C-Finite difference equations be respectively Wherein, aci,jIndicate the ultrasonic acoustic at i-node j moment
Speed；vi,jIndicate the flow velocity at i-node j moment；When Sudden Leakage point is the i-node j moment, and assume at j-1 moment pipeline fluids
It in steady state condition, is remained unchanged in the parameter of j+1 moment i+2 nodes, is solved to obtain leakage wink according to the finite difference equations
Become information propagation attenuation computation model into Wherein, the boundary condition of leakage point meetsFmFor the information transition of pipe end, Δ FmFor the information of pipe end
The increment of transition, acmAnd vmThe respectively velocity of sound and flow velocity of pipe end, LxFor leakage point to the distance of pipeline first section, Δ x is
Run of designing, Δ x=awΔ t, FsFor the information transition of pipeline head end, Δ FsFor the increment of the information transition of pipeline head end, acs
And vsThe respectively velocity of sound and flow velocity of pipeline head end, LxFor leakage point to the distance of pipeline first section, Fu_For the information before leakage point
Transition.
Setting is needed to manage in order to ensure the foundation that the equation of motion, contiguity equation can be stablized the embodiment of the present invention
Road caliber, tubing, wall thickness and tube fluid property remain unchanged within the scope of pipe range.Join for the equation of motion and continuity equation
It is solved after vertical, is one group of hyperboloid partial differential equation to composition, is solved by characteristics method.
In addition, by being the positioning to the leakage point of oil pipeline in this present embodiment, the working fluid in pipeline is hydro carbons
The oil of substance, acoustic characteristic meet the rule of Kneser liquid, i.e., at a certain temperature, the velocity of sound is linear with increasing for pressure
Ground increases, to existWherein, p is pressure, acFor the velocity of sound in liquid, a0For in liquid under one atmospheric pressure of room temperature
The velocity of sound, K are proportionality coefficient, can be inferred to the equation of motion and continuity equationWhen for using characteristics method, need to ignore
Elevation changes, and the C established+Characteristic curve and its characteristic equation are And C-It is special
Sign line and its characteristic equation are Then utilize finite difference calculus to above-mentioned C+、C-
Characteristic equation carry out numerical solution, characteristic curve network schematic diagram, as shown in Figure 3, wherein pipeline is divided into N section, and every section is managed
A length of Δ x, along duct length x-axis marked as 1,2 ..., N+1 indicates the node between pipeline section respectively, along time t axis decile
It is Δ t at time horizon.
In the embodiment of the present invention, according to the flow velocity before pipe leakage, the flow velocity that leakage counts out money, the flow velocity after leakage, leakage
Flow velocity after point, the velocity of sound before leakage point, the velocity of sound after leakage point, and the boundary condition for limiting leakage point needs to meet：To obtain the sonic velocity change of leakage pointObtain the velocity of sound and leakage rate, the density of medium, suction wave wave that leakage generates
Speed is related with the acoustic velocity value of normal atmosphere pressure.
Further solve, along C+(i+1, j+1) is arrived from point (i+, j), obtaining equation isBoth sides take increment to obtain simultaneously,And it enablesFu+For the information transition after leakage point.Then due to leakage
Transient information is to be transmitted from leakage point to pipe ends, therefore remained unchanged in the parameter of j+1 moment i+2 nodes, obtains Δ
aci+1,j+1=Ka0ρawΔvi+1,j+1, wherein the Parameters variation of j+2 moment i+2 nodes can be expressed asFurther obtainIn conjunction with Fu+Definition,
It obtainsSimilarly, it can be expressed as in the Parameters variation of j+3 moment i+3 nodesTo push over out
AndWherein, F is definedu-It is for the information transition before leakage point
Further, in the embodiment of the present invention, in order to avoid at leakage point due to acoustic velocity value, flow speed value change and cause
The unstability for leaking transient information needs to merge the leakage transient information at first and last end into row information, and the basis is to described
Head end information after the leakage transient information at first and last end is merged into row information merges signal and final word fusion signal determination
Go out before delay time, the method further includes：Utilize the head end velocity of sound, flow and the end velocity of sound, stream when pipeline leaks
Correspondence between amount merges the leakage transient information at the first and last end into row information, obtains the leakage transition of first and end
Information FsAnd FmRespectively Fs=acs-Ka0ρawvs, Fm=acm+Ka0ρawvm。
For the embodiment of the present invention, since acoustic velocity value at leakage point and flow speed value all change, and with certain velocity of wave
Signal propagation is carried out to pipeline first and last end, leakage transient information can carry out signal by ultrasonic sensor and flowmeter and adopt
Collection, therefore, leakage transient information include not only velocity of sound signal, but also include flow signal, are declined according to leakage transient information propagation
Subtract two formula in computation model, can infer when pipeline leaks, the head end velocity of sound is reduced, and flow increases, and end
The velocity of sound is held to reduce, flow is reduced, therefore, the leakage transient information F at first and last endsAnd FmIt is all corresponding to reduce.In addition, going out by inference
Fs=acs-Ka0ρawvs、Fm=acm+Ka0ρawvmMonitor whether pipeline goes out to measure the leakage transient information at pipeline first and last end
Now leak.
In the embodiment of the present invention, in order to calculating the further of distance of the leakage point apart from pipeline head end ultrasonic sensor
Refinement, the head end information after the basis merges the leakage transient information at the first and last end into row information merge signal and end
Client information fusion signal determines delay time, and calculates leakage point apart from pipeline head end using preset leak point positioning formula
The distance of ultrasonic sensor includes：The leakage transient information F of head end after being merged respectively to information using wavelet analysis methodsWith
And the leakage transient information F of endmMutation time point is extracted, according to cross-correlation analysis method and pre-set time window to institute
It states the rough time to peak search of mutation time point progress and determined for the first delay time, and according to cross-correlation analysis hair and in advance
The time window of setting is to determine for the second delay time by the accurate time to peak search of progress to the mutation, and described first is prolonged
Slow the sum of time and second delay time are determined as delay time；Utilize preset leak point positioning formula
Calculate distance of the leakage point apart from pipeline head end ultrasonic sensor, wherein the tDFor delay time, L is duct length.
For the embodiment of the present invention, in order to extract the peak point of time catastrophe point, and rationally choosing by cross-correlation function
Position and the length for selecting two signal segments as cross-correlation, to ensure that cross-correlation function has apparent peak point, first with
Wavelet analysis method extracts the catastrophe point (time) that pipeline first and last end measures fusion signal, and time window is arranged according to the position of catastrophe point
Signal length is merged to extract pipeline first and last end.Start timing when the catastrophe point that first merges signal reaches certain threshold value,
Stop timing when the catastrophe point that second merges signal reaches certain threshold value, then records interval time.
In the embodiment of the present invention, it is assumed that the time span for acquiring signal is 0~t0, the time that wavelet packet analysis measures is t1,
So measure in wavelet packet and signal acquisition length intercepted on time basis, retain head end fusion signal corresponds to the time be 0
~(t0-t1) signal length；Terminal fusion signal corresponds to t1~t0Signal length.Corresponding time span signal will be retained to make
For the signal data of new time series.The schematic diagram of pipeline first and last segment signal catastrophe point time difference, as shown in Figure 4.For the first time
Rough time to peak is carried out using cross-correlation analysis to search for, it is 0~(t respectively to acquire first terminal fusion signal time0-t1) and
t1~t0, it is spaced gathered data and carries out cross-correlation analysis, estimates the delay time t of signalD1.Second using correlation analysis into
The accurate time to peak search of row, in the delay time t of the first correlation analysis estimationD1It is front and back that certain sampling period is respectively taken to carry out the
Secondary correlation analysis retrieves maximum value tD2.The delay time of estimation after searching for twice is tD=tD1+tD2, tD1For for the first time
Rough search cross-correlation peak value, tD2For second of coarse-fine fine searching cross-correlation peak value.
Then when pipe leakage occurs, the cross-correlation function that signal is merged by pipeline head and end information is analyzed, and is obtained
The peak value R of signal is merged to two client informationssm(τ) obtains delay time tD, wherein leakage transient information is gone by wavelet packet
After making an uproar, the signal at first and last end is FsAnd FmCross-correlation function beT is the leakage velocity of sound
The propagation period between two sonacs.Finally utilize preset leak point positioning formulaIt calculates and lets out
Distance of the leak source apart from pipeline head end ultrasonic sensor.
For the embodiment of the present invention, specific application scenarios can be as follows, but not limited to this, including：Select pipeline
Length is 2000 meters, and internal diameter of the pipeline is DN 50, and the first and last end of pipeline is separately installed with ultrasonic velocity sensor, flow sensing
Device, temperature sensor, and the flow for setting current fluid conveying is installed different leakage apertures and is let out to simulate difference as 12m3/h
Leakage degree, according to live physical condition, the lift of pipe outlet pump is 120 meters of H2O water columns.In addition the density of selection fluid is
1000kg/m3, ultrasonic velocity 1520m/s, suction wave velocity of wave is 1120m/s, and wherein the value of K is 2.13 × 10-6.When away from
Leak valve at 400 meters from pipeline head end is opened, and Small leak amount is 2.5L/min at this time, in order to by flow, the velocity of sound, pressure etc. no
Same parameter is compared, the Daubechies 4 of flow signal, ultrasonic velocity signal, pressure signal through wavelet packet of acquisition
(db4) it decomposes and carries out signal de-noising, it has been found by contrast that when Small leak occurs, corresponding discharge, velocity of sound signal and pressure signal difference
As shown in Fig. 5, Fig. 6, Fig. 7, then by leaking the pipe obtained after transient information propagation attenuation computation model and signal fused
Road head end information fusion signal, pipe end information fusion signal difference are as shown in Figure 8, Figure 9, wherein there are one fusion signals
Obviously inflection point.Then information fusion signal is subjected to cross-correlation analysis.First, the information fusion sigtnal interval 1 is adopted
Sampling point is sampled, and is then carried out first time cross-correlation analysis, is obtained correlation analysis peak value, i.e. delay time tD1, such as Figure 10 institutes
The first time cross-correlation analysis shown, then, in time to peak tD1Around, i.e., 50 sampled points are respectively taken in prime information fusion signal,
Cross-correlation analysis is carried out again, obtains time to peak tD2, second of cross-correlation analysis as shown in figure 11 finally obtain delay time
tDFor -1.22s, merge signal improves cross-correlation function, the classical cross-correlation function that merges signal, ultrasonic velocity signal
The classical cross-correlation function of classical cross-correlation function and pipeline first and last side pressure force signal is as shown in figure 12.Since Small leak amount is
2.5L/min, at 400 meters of pipeline head end, the first and last end signal catastrophe point cross-correlation time difference should be -1.2s. can be with from Figure 12
Find out, the cross-correlation analysis that improves that information merges signal calculates the time difference of acquisition as -1.22s, and information merges the classics of signal mutually
The time difference that correlation analysis calculates acquisition is -0.8s, and the time difference that the classical cross-correlation analysis calculating of velocity of sound signal obtains is -0.52s,
And the peak point that the classical cross-correlation analysis of pressure calculates acquisition is more, which is peak point caused by leakage for bad determination, from
And it is the pipeline first and last client information under 2.5L/min, 1.5L/min and 1L/min according to leakage rate shown in Figure 13, Figure 14, is given
Signal contrast curve is merged, with the reduction of leakage rate, the amplitude of information fusion waveform is gradually reduced.When leakage rate reaches 1L/
When min, the more difficult identification of waveform inflection point is determined position error and is accordingly increased so causing to leak, and passes through and improves cross-correlation analysis estimation
The time difference, and leakage point position is calculated using preset leak point positioning formula, it obtains being happened at apart from pipeline head end when leak position
600m, 800m and 1400m, the position location error directly obtained compared to existing classical cross-correlation analysis are obviously reduced.
The present invention provides a kind of localization methods of pipe leakage, with existing under Small leak and the pressure signal simmered
Drop is unintelligible along inflection point, and it is larger to cause to fail to report higher and position error, and punches installation pressure sensor using pipeline and can bring
Danger, the locating effect for influencing pipe leakage are compared, and the embodiment of the present invention is established by the mechanism of transmission of ultrasonic velocity
The velocity of sound generates and propagation attenuation computation model, and ultrasonic velocity and flow signal carried out signal fused so that information
Fusion signal can more burst signal knee, and pipeline first and last end can be more accurately extracted by improving cross-correlation analysis method
Information merges the time difference of signal, so as to improve the positioning accuracy of Small leak, fast and effeciently positions Small leak point, improves pipe
The efficiency of road leakage positioning.
Further, as the realization to method shown in above-mentioned Fig. 1, an embodiment of the present invention provides a kind of pipe leakages
Positioning device, as shown in figure 15, which includes：Acquiring unit 21, the first computing unit 22, the second computing unit 23.
Acquiring unit 21, for obtaining the ultrasonic velocity sensor, the flow sensor that are separately mounted to pipeline first and last end
Flow velocity signal, the velocity of sound signal of acquisition, the flow velocity signal include the flow velocity signal of pipeline head end and the flow velocity letter of pipe end
Number, the velocity of sound signal includes the velocity of sound signal of pipeline head end and the velocity of sound signal of pipe end；
First computing unit 22, for according to the leakage transient information propagation attenuation computation model pre-established according to
Flow velocity signal, the velocity of sound signal calculate the leakage transient information at first and last end, and the leakage transient information propagation attenuation calculates
Model is the pipeline first and last determined according to the equation of motion and continuity equation combination ultrasonic velocity attenuation characteristic of pipeline
The transient information computation model at end；
Second computing unit 23, for the head end after being merged into row information according to the leakage transient information to the first and last end
Information merges signal and final word fusion signal determines delay time, and is calculated using preset leak point positioning formula
Distance of the leakage point apart from pipeline head end ultrasonic sensor.
Further, as shown in figure 16, described device further includes：First establishing unit 24, second establish unit 25, solve
Unit 26.
First establishing unit 24, the Transient Flow in Pipeline equation for establishing liquid according to the equation of motion, continuity equation
ForWherein, v is flow velocity, and ρ is fluid density, and D is
Internal diameter of the pipeline, alpha levels axle clamp angle, λ are darcy hydraulic simulation experiment, and x is the axial distance along pipeline, and t is time, awFor negative pressure
Wave velocity of wave, g are acceleration of gravity, acFor the velocity of sound in liquid, a0For the velocity of sound in liquid under one atmospheric pressure of room temperature, K is ratio system
Number；
Second establishes unit 25, ignores elevation variation, the then node i at the j moment established according to gsin α=0 for working as
Along C+And C-Finite difference equations be respectively Wherein, aci,jIndicate the ultrasonic acoustic at i-node j moment
Speed；vi,jIndicate the flow velocity at i-node j moment；
Unit 26 is solved, for being the i-node j moment when Sudden Leakage point, and assumes that j-1 moment pipeline fluids are in steady
State operating mode is remained unchanged in the parameter of j+1 moment i+2 nodes, is solved to obtain leakage transition letter according to the finite difference equations
Ceasing propagation attenuation computation model is Wherein, the boundary condition of leakage point meetsFmFor the information transition of pipe end, Δ FmFor the information of pipe end
The increment of transition, acmAnd vmThe respectively velocity of sound and flow velocity of pipe end, LxFor leakage point to the distance of pipeline first section, Δ x is
Run of designing, Δ x=awΔ t, FsFor the information transition of pipeline head end, Δ FsFor the increment of the information transition of pipeline head end, acs
And vsThe respectively velocity of sound and flow velocity of pipeline head end, LxFor leakage point to the distance of pipeline first section, Fu-For the information before leakage point
Transition.
Further, as shown in figure 16, described device further includes：
Integrated unit 27, for utilizing when pipeline leaks between the head end velocity of sound, flow and the end velocity of sound, flow
Correspondence merges the leakage transient information at the first and last end into row information, obtains the leakage transient information F of first and endsAnd Fm
Respectively Fs=acs-Ka0ρawvs, Fm=acm+Ka0ρawvm。
Second computing unit 23 is specifically used for the leakage of the head end after being merged respectively to information using wavelet analysis method
Transient information FsAnd the leakage transient information F of endmMutation time point is extracted, according to cross-correlation analysis method and is pre-set
Time window the search of rough time to peak carried out to the mutation time point determined for the first delay time, and according to cross-correlation point
Analysis hair and pre-set time window are to determine for the second delay time by the accurate time to peak search of progress to the mutation,
The sum of first delay time and second delay time are determined as delay time；Utilize preset leak point positioning formulaCalculate distance of the leakage point apart from pipeline head end ultrasonic sensor, wherein the tDFor delay time,
L is duct length.
The present invention provides the propagation machines that a kind of positioning device of pipe leakage, the embodiment of the present invention pass through ultrasonic velocity
Reason establishes velocity of sound generation and propagation attenuation computation model, and ultrasonic velocity and flow signal has been carried out signal fused, makes
Signal knee can more be burst by obtaining information fusion signal, and pipeline can be more accurately extracted by improving cross-correlation analysis method
First and last client information merges the time difference of signal, so as to improve the positioning accuracy of Small leak, fast and effeciently positions Small leak point,
Improve the efficiency of pipe leakage positioning.
A kind of storage medium is provided according to an embodiment of the present invention, and it is executable that the storage medium is stored at least one
Instruction, the computer executable instructions can perform the localization method of the pipe leakage in above-mentioned any means embodiment.
Figure 17 shows a kind of structural schematic diagram of the terminal provided according to an embodiment of the present invention, the present invention is specifically real
Example is applied not limit the specific implementation of terminal.
As shown in figure 17, which may include：Processor (processor) 502, communication interface
(Communications Interface) 504, memory (memory) 506 and communication bus 508.
Wherein：Processor 502, communication interface 504 and memory 506 complete mutual lead to by communication bus 508
Letter.
Communication interface 504, for being communicated with the network element of miscellaneous equipment such as client or other servers etc..
Processor 502, for executing program 510, in the localization method embodiment that can specifically execute above-mentioned pipe leakage
Correlation step.
Specifically, program 510 may include program code, which includes computer-managed instruction.
Processor 502 may be central processor CPU or specific integrated circuit ASIC (Application
Specific Integrated Circuit), or be arranged to implement the integrated electricity of one or more of the embodiment of the present invention
Road.The one or more processors that terminal includes can be same type of processor, such as one or more CPU；Can also be
Different types of processor, such as one or more CPU and one or more ASIC.
Memory 506, for storing program 510.Memory 506 may include high-speed RAM memory, it is also possible to further include
Nonvolatile memory (non-volatile memory), for example, at least a magnetic disk storage.
Program 510 specifically can be used for so that processor 502 executes following operation：
Obtain be separately mounted to pipeline first and last end ultrasonic velocity sensor, flow sensor acquisition flow velocity signal,
Velocity of sound signal, the flow velocity signal include the flow velocity signal of pipeline head end and the flow velocity signal of pipe end, the velocity of sound signal
The velocity of sound signal of velocity of sound signal and pipe end including pipeline head end；
Believed according to the flow velocity signal, the velocity of sound according to the leakage transient information propagation attenuation computation model pre-established
The leakage transient information at first and last end number is calculated, the leakage transient information propagation attenuation computation model is the movement according to pipeline
The transient information computation model at the pipeline first and last end that equation and continuity equation combination ultrasonic velocity attenuation characteristic are determined；
Head end information fusion signal after being merged into row information according to the leakage transient information to the first and last end and end
Client information fusion signal determines delay time, and calculates leakage point apart from pipeline head end using preset leak point positioning formula
The distance of ultrasonic sensor.
Obviously, those skilled in the art should be understood that each module of the above invention or each step can be with general
Computing device realize that they can be concentrated on a single computing device, or be distributed in multiple computing devices and formed
Network on, optionally, they can be realized with the program code that computing device can perform, it is thus possible to which they are stored
It is performed by computing device in the storage device, and in some cases, it can be with different from shown in sequence execution herein
The step of going out or describing, either they are fabricated to each integrated circuit modules or by them multiple modules or
Step is fabricated to single integrated circuit module to realize.In this way, the present invention is not limited to any specific hardware and softwares to combine.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, any made by repair
Change, equivalent replacement, improvement etc., should all include within protection scope of the present invention.
Claims (10)
1. a kind of localization method of pipe leakage, which is characterized in that including：
Obtain the ultrasonic velocity sensor for being separately mounted to pipeline first and last end, the flow velocity signal of flow sensor acquisition, the velocity of sound
Signal, the flow velocity signal include the flow velocity signal of pipeline head end and the flow velocity signal of pipe end, and the velocity of sound signal includes
The velocity of sound signal of pipeline head end and the velocity of sound signal of pipe end；
According to the leakage transient information propagation attenuation computation model pre-established according to the flow velocity signal, the velocity of sound signal meter
The leakage transient information at first and last end is calculated, the leakage transient information propagation attenuation computation model is the equation of motion according to pipeline
And the transient information computation model at pipeline first and last end that continuity equation combination ultrasonic velocity attenuation characteristic is determined；
Head end information fusion signal and end after being merged into row information according to the leakage transient information to the first and last end are believed
Breath fusion signal determines delay time, and calculates leakage point apart from pipeline head end ultrasound using preset leak point positioning formula
The distance of wave sensor.
2. according to the method described in claim 1, it is characterized in that, described obtain the ultrasonic wave for being separately mounted to pipeline first and last end
Before sonic sensor, the flow velocity signal of flow sensor acquisition, velocity of sound signal, the method further includes：
The Transient Flow in Pipeline equation that liquid is established according to the equation of motion, continuity equation isWherein, v is flow velocity, and ρ is fluid density, and D is pipe
Road internal diameter, alpha levels axle clamp angle, λ are darcy hydraulic simulation experiment, and x is the axial distance along pipeline, and t is time, awFor suction wave
Velocity of wave, g are acceleration of gravity, acFor the velocity of sound in liquid, a0For the velocity of sound in liquid under one atmospheric pressure of room temperature, K is proportionality coefficient；
When ignoring elevation variation, then the node i at the j moment established according to gsin α=0 is along C+And C-Finite difference equations be respectively
Wherein, aci,jIndicate the ultrasonic velocity at i-node j moment；vi,jIndicate the flow velocity at i-node j moment；
When Sudden Leakage point is the i-node j moment, and assume that j-1 moment pipeline fluids are in steady state condition, in j+1 moment i+2
The parameter of node remains unchanged, and solving to obtain leakage transient information propagation attenuation computation model according to the finite difference equations is
Wherein, the boundary condition of leakage point meetsFmFor the information of pipe end
Transition, Δ FmFor the increment of the information transition of pipe end, acmAnd vmThe respectively velocity of sound and flow velocity of pipe end, LxFor leakage
For point to the distance of pipeline first section, Δ x is run of designing, Δ x=awΔ t, FsFor the information transition of pipeline head end, Δ FsFor pipe
The increment of the information transition of road head end, acsAnd vsThe respectively velocity of sound and flow velocity of pipeline head end, LxFor leakage point to pipeline first section
Distance, Fu-For the information transition before leakage point.
3. according to the method described in claim 1, it is characterized in that, the basis to the leakage transient information at the first and last end into
Before head end information fusion signal and final word fusion signal after row information fusion determine delay time, the method
Further include：
Using the correspondence when pipeline leaks between the head end velocity of sound, flow and the end velocity of sound, flow to the first and last
The leakage transient information at end is merged into row information, obtains the leakage transient information F of first and endsAnd FmRespectively Fs=acs-Ka0ρ
awvs, Fm=acm+Ka0ρawvm。
4. according to the method described in claim 3, it is characterized in that, the basis to the leakage transient information at the first and last end into
Head end information fusion signal and final word fusion signal after row information fusion determine delay time, and are let out using preset
Leak source ranging formula calculates distance of the leakage point apart from pipeline head end ultrasonic sensor：
The leakage transient information F of head end after being merged respectively to information using wavelet analysis methodsAnd the leakage transient information of end
FmMutation time point is extracted, the mutation time point is carried out according to cross-correlation analysis method and pre-set time window rough
Time to peak search determined for the first delay time, and according to cross-correlation analysis hair and pre-set time window to described prominent
Change is to determine for the second delay time by the accurate time to peak search of progress, and first delay time and described second are postponed
The sum of time is determined as delay time；
Utilize preset leak point positioning formulaLeakage point is calculated apart from pipeline head end ultrasonic sensor
Distance, wherein the tDFor delay time, L is duct length.
5. a kind of positioning device of pipe leakage, which is characterized in that including：
Acquiring unit, for obtaining the ultrasonic velocity sensor for being separately mounted to pipeline first and last end, flow sensor acquires
Flow velocity signal, velocity of sound signal, the flow velocity signal include the flow velocity signal of pipeline head end and the flow velocity signal of pipe end, described
Velocity of sound signal includes the velocity of sound signal of pipeline head end and the velocity of sound signal of pipe end；
First computing unit, for being believed according to the flow velocity according to the leakage transient information propagation attenuation computation model pre-established
Number, the velocity of sound signal calculate the leakage transient information at first and last end, the leakage transient information propagation attenuation computation model is
According to the wink at the pipeline first and last end that the equation of motion of pipeline and continuity equation combination ultrasonic velocity attenuation characteristic are determined
Become information computation model；
Second computing unit melts for the head end information after being merged into row information according to the leakage transient information to the first and last end
It closes signal and final word fusion signal determines delay time, and leakage point is calculated using preset leak point positioning formula
Distance apart from pipeline head end ultrasonic sensor.
6. device according to claim 5, which is characterized in that described device further includes：
First establishing unit, the Transient Flow in Pipeline equation for establishing liquid according to the equation of motion, continuity equation areWherein, v is flow velocity, and ρ is fluid density, and D is pipe
Road internal diameter, alpha levels axle clamp angle, λ are darcy hydraulic simulation experiment, and x is the axial distance along pipeline, and t is time, awFor suction wave
Velocity of wave, g are acceleration of gravity, acFor the velocity of sound in liquid, a0For the velocity of sound in liquid under one atmospheric pressure of room temperature, K is proportionality coefficient；
Second establishes unit, for when ignoring elevation variation, then the node i at the j moment established according to gsin α=0 to be along C+And C-
Finite difference equations be respectively Wherein, aci,jIndicate the ultrasonic acoustic at i-node j moment
Speed；vi,jIndicate the flow velocity at i-node j moment；
Unit is solved, for being the i-node j moment when Sudden Leakage point, and assumes that j-1 moment pipeline fluids are in stable state work
Condition is remained unchanged in the parameter of j+1 moment i+2 nodes, is solved to obtain leakage transient information biography according to the finite difference equations
Broadcasting decay calculation model is Wherein, the boundary condition of leakage point meetsFmFor the information transition of pipe end, Δ FmFor the information of pipe end
The increment of transition, acmAnd vmThe respectively velocity of sound and flow velocity of pipe end, LxFor leakage point to the distance of pipeline first section, Δ x is
Run of designing, Δ x=awΔ t, FsFor the information transition of pipeline head end, Δ FsFor the increment of the information transition of pipeline head end, acs
And vsThe respectively velocity of sound and flow velocity of pipeline head end, LxFor leakage point to the distance of pipeline first section, Fu-For the information before leakage point
Transition.
7. device according to claim 6, which is characterized in that described device further includes：
Integrated unit, for being closed using the head end velocity of sound, flow are corresponding between the end velocity of sound, flow when pipeline leaks
It is to be merged into row information to the leakage transient information at the first and last end, obtains the leakage transient information F of first and endsAnd FmRespectively
Fs=acs-Ka0ρawvs, Fm=acm+Ka0ρawvm。
8. device according to claim 7, which is characterized in that
Second computing unit, the leakage transition for being specifically used for the head end after being merged respectively to information using wavelet analysis method are believed
Cease FsAnd the leakage transient information F of endmMutation time point is extracted, according to cross-correlation analysis method and pre-set time
Window carries out the search of rough time to peak to the mutation time point and determined for the first delay time, and according to cross-correlation analysis send out with
And pre-set time window is to determine for the second delay time by the accurate time to peak search of progress to the mutation, it will be described
The sum of first delay time and second delay time are determined as delay time；Utilize preset leak point positioning formulaCalculate distance of the leakage point apart from pipeline head end ultrasonic sensor, wherein the tDFor delay time,
L is duct length.
9. a kind of storage medium, it is stored with an at least executable instruction in the storage medium, the executable instruction makes processing
Device executes the corresponding operation of localization method of the pipe leakage as described in any one of claim 1-4.
10. a kind of terminal, including：Processor, memory, communication interface and communication bus, the processor, the memory and
The communication interface completes mutual communication by the communication bus；
The memory makes the processor execute as right is wanted for storing an at least executable instruction, the executable instruction
Ask the corresponding operation of the localization method of the pipe leakage described in any one of 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810478059.XA CN108662442B (en) | 2018-08-23 | 2018-08-23 | The localization method and device of pipe leakage |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810478059.XA CN108662442B (en) | 2018-08-23 | 2018-08-23 | The localization method and device of pipe leakage |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108662442A true CN108662442A (en) | 2018-10-16 |
| CN108662442B CN108662442B (en) | 2019-08-06 |
Family
ID=63776711
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810478059.XA Active CN108662442B (en) | 2018-08-23 | 2018-08-23 | The localization method and device of pipe leakage |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108662442B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109340585A (en) * | 2018-10-26 | 2019-02-15 | 北京化工大学 | A kind of pipeline leakage testing and localization method |
| CN109469834A (en) * | 2018-12-07 | 2019-03-15 | 中国石油大学(北京) | A kind of fluid pipeline leakage detection method, apparatus and system |
| CN110131594A (en) * | 2019-05-20 | 2019-08-16 | 辽宁石油化工大学 | A kind of leakage locating method based on ultrasonic velocity variation |
| CN110159937A (en) * | 2019-07-02 | 2019-08-23 | 北京中计新业科技发展有限公司 | A kind of pipeline leakage detection method of naval vessel pipe network |
| CN110513603A (en) * | 2019-08-13 | 2019-11-29 | 常州大学 | A kind of nonmetal pipeline leakage locating method based on inverse Transient Analysis |
| CN112413413A (en) * | 2020-11-20 | 2021-02-26 | 成都艾斯皮埃尔科技有限公司 | Pipeline leakage monitoring and positioning method combining deep learning and multiple measurement technology |
| CN114508704A (en) * | 2020-11-16 | 2022-05-17 | 中国石油天然气股份有限公司 | Pipeline leakage detection method and device and storage medium |
| CN117762063A (en) * | 2023-12-04 | 2024-03-26 | 安徽省江北华衍水务有限公司 | Waterproof safety supervision control system for water works |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003214600A (en) * | 2002-01-17 | 2003-07-30 | Jfe Engineering Kk | Fluid passing method to porous body and its device |
| CN101413628A (en) * | 2008-11-26 | 2009-04-22 | 哈尔滨工业大学 | Method for performing gas pipeline leakage position by using instant change on-line diagnosis coupling excitation frequency response |
| CN101960303A (en) * | 2008-11-11 | 2011-01-26 | 韩国电力系统开发设计有限公司 | Apparatus for measuring fluid leakage from a valve using ultrasonic wave, sound, and temperature variations, and method for measuring fluid leakage using same |
| CN104747912A (en) * | 2015-04-23 | 2015-07-01 | 重庆邮电大学 | Fluid conveying pipe leakage acoustic emission time-frequency positioning method |
| CN105987284A (en) * | 2015-01-27 | 2016-10-05 | 中国石油天然气股份有限公司 | Pipeline leakage monitoring method and device |
| CN106934511A (en) * | 2015-12-30 | 2017-07-07 | 海能达通信股份有限公司 | Method for scheduling task and system, server, user terminal |
| CN108036201A (en) * | 2017-12-21 | 2018-05-15 | 廊坊市蓝德采油技术开发有限公司 | A kind of Leak Detection in Oil Pipeline Using method based on negative pressure wave method and traffic trends method |
-
2018
- 2018-08-23 CN CN201810478059.XA patent/CN108662442B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003214600A (en) * | 2002-01-17 | 2003-07-30 | Jfe Engineering Kk | Fluid passing method to porous body and its device |
| CN101960303A (en) * | 2008-11-11 | 2011-01-26 | 韩国电力系统开发设计有限公司 | Apparatus for measuring fluid leakage from a valve using ultrasonic wave, sound, and temperature variations, and method for measuring fluid leakage using same |
| CN101413628A (en) * | 2008-11-26 | 2009-04-22 | 哈尔滨工业大学 | Method for performing gas pipeline leakage position by using instant change on-line diagnosis coupling excitation frequency response |
| CN105987284A (en) * | 2015-01-27 | 2016-10-05 | 中国石油天然气股份有限公司 | Pipeline leakage monitoring method and device |
| CN104747912A (en) * | 2015-04-23 | 2015-07-01 | 重庆邮电大学 | Fluid conveying pipe leakage acoustic emission time-frequency positioning method |
| CN106934511A (en) * | 2015-12-30 | 2017-07-07 | 海能达通信股份有限公司 | Method for scheduling task and system, server, user terminal |
| CN108036201A (en) * | 2017-12-21 | 2018-05-15 | 廊坊市蓝德采油技术开发有限公司 | A kind of Leak Detection in Oil Pipeline Using method based on negative pressure wave method and traffic trends method |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109340585A (en) * | 2018-10-26 | 2019-02-15 | 北京化工大学 | A kind of pipeline leakage testing and localization method |
| CN109469834A (en) * | 2018-12-07 | 2019-03-15 | 中国石油大学(北京) | A kind of fluid pipeline leakage detection method, apparatus and system |
| CN110131594A (en) * | 2019-05-20 | 2019-08-16 | 辽宁石油化工大学 | A kind of leakage locating method based on ultrasonic velocity variation |
| CN110159937A (en) * | 2019-07-02 | 2019-08-23 | 北京中计新业科技发展有限公司 | A kind of pipeline leakage detection method of naval vessel pipe network |
| CN110513603A (en) * | 2019-08-13 | 2019-11-29 | 常州大学 | A kind of nonmetal pipeline leakage locating method based on inverse Transient Analysis |
| CN110513603B (en) * | 2019-08-13 | 2021-09-28 | 常州大学 | Non-metal pipeline leakage positioning method based on inverse transient analysis method |
| CN114508704A (en) * | 2020-11-16 | 2022-05-17 | 中国石油天然气股份有限公司 | Pipeline leakage detection method and device and storage medium |
| CN114508704B (en) * | 2020-11-16 | 2024-04-30 | 中国石油天然气股份有限公司 | Pipeline leakage detection method and device and storage medium |
| CN112413413A (en) * | 2020-11-20 | 2021-02-26 | 成都艾斯皮埃尔科技有限公司 | Pipeline leakage monitoring and positioning method combining deep learning and multiple measurement technology |
| CN112413413B (en) * | 2020-11-20 | 2022-03-01 | 成都艾斯皮埃尔科技有限公司 | Pipeline leakage monitoring and positioning method combining deep learning and multiple measurement technology |
| CN117762063A (en) * | 2023-12-04 | 2024-03-26 | 安徽省江北华衍水务有限公司 | Waterproof safety supervision control system for water works |
| CN117762063B (en) * | 2023-12-04 | 2024-06-04 | 安徽省江北华衍水务有限公司 | Waterproof safety supervision control system for water works |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108662442B (en) | 2019-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108662442B (en) | The localization method and device of pipe leakage | |
| Covas et al. | Case studies of leak detection and location in water pipe systems by inverse transient analysis | |
| Li et al. | A review of methods for burst/leakage detection and location in water distribution systems | |
| Covas et al. | Standing wave difference method for leak detection in pipeline systems | |
| CN103939749B (en) | Flow circuits based on big data leakage intelligent adaptive monitoring system and method | |
| Kalpakli Vester et al. | POD analysis of the turbulent flow downstream a mild and sharp bend | |
| CN1184931A (en) | Method and apparatus for detecting and positioning leakage of fluid transferring pipeline | |
| Tuck et al. | Analysis of transient signals in simple pipeline systems with an extended blockage | |
| CN112711844A (en) | Pipeline leakage positioning, leakage amount early warning and automatic processing method and system | |
| Nerella et al. | Fluid transients and wave propagation in pressurized conduits due to valve closure | |
| CN111853552A (en) | Heat supply pipeline leakage online detection positioning method | |
| CN107463737A (en) | A kind of computational methods and device of fluid pipeline leakage rate | |
| Azpiroz et al. | CFD modelling of bypass pigs with a deflector disk | |
| CN108197070A (en) | Natural gas not exclusively blocks pipeline method for numerical simulation | |
| CN101308364B (en) | Water supply pipe network affair model modelling approach based on transient flow analysis | |
| CN107014451A (en) | The method of ultrasonic flow sensor coefficient is speculated based on generalized regression nerve networks | |
| CN109737317A (en) | The infrasound positioning system and method for fluid line leakage | |
| Abdelhafidh et al. | Hybrid mechanism for remote water pipeline monitoring system | |
| CN111062127A (en) | Pipeline leakage point detection method and device, storage medium and terminal | |
| CN109340583A (en) | Heating network leakage monitoring system and method | |
| Covas et al. | An assessment of the application of inverse transient analysis for leak detection: Part II—Collection and application of experimental data | |
| BR112020012925B1 (en) | METHOD AND SYSTEM FOR DETECTING AN IRREGULARITY WITHIN A FLUID CHANNEL, AND COMPUTER READABLE STORAGE MEDIUM | |
| Shehadeh et al. | Modelling the effect of incompressible leakage patterns on rupture area in pipeline | |
| Abdullahi et al. | Simulation and detection of blockage in a pipe under mean fluid flow using acoustic wave propagation technique | |
| CN113063100A (en) | Pipeline leakage detection method and device |
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 |