CN109344784A - A kind of leakage orienting method merging hydrophone and Wavelet Denoising Method - Google Patents

A kind of leakage orienting method merging hydrophone and Wavelet Denoising Method Download PDF

Info

Publication number
CN109344784A
CN109344784A CN201811186445.8A CN201811186445A CN109344784A CN 109344784 A CN109344784 A CN 109344784A CN 201811186445 A CN201811186445 A CN 201811186445A CN 109344784 A CN109344784 A CN 109344784A
Authority
CN
China
Prior art keywords
hydrophone
delay value
wavelet
time delay
leakage
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.)
Pending
Application number
CN201811186445.8A
Other languages
Chinese (zh)
Inventor
张楠
李舒
郭欣伟
江红
张攀
殷会娟
吴凯
李凌琪
韩金旭
孙维营
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yellow River Institute of Hydraulic Research
Original Assignee
Yellow River Institute of Hydraulic Research
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 Yellow River Institute of Hydraulic Research filed Critical Yellow River Institute of Hydraulic Research
Priority to CN201811186445.8A priority Critical patent/CN109344784A/en
Publication of CN109344784A publication Critical patent/CN109344784A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02Preprocessing
    • G06F2218/04Denoising
    • G06F2218/06Denoising by applying a scale-space analysis, e.g. using wavelet analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/22Source localisation; Inverse modelling

Abstract

The invention belongs to aqueduct leakage orienting technical fields, in particular to a kind of leakage orienting method for merging hydrophone and Wavelet Denoising Method, one hydrophone is respectively installed included in the two ends of aqueduct, the information of acoustic wave generated in mine water transmission process is obtained using hydrophone, obtains untreated time series original acoustic wave information source;Wavelet coefficient noise reduction process model is established, noise-removed filtering is carried out to the original non-stationary signal source of acquisition;Setup delay value model determines time delay value to the signal after denoising;Location model is established, using the time delay value of above-mentioned acquisition, positions breakthrough.The present invention improves the accuracy and science of leakage orienting, improves the utilization rate of mine water resource.

Description

A kind of leakage orienting method merging hydrophone and Wavelet Denoising Method
Technical field
The invention belongs to aqueduct leakage orienting technical field, in particular to a kind of fusion hydrophone and Wavelet Denoising Method Leakage orienting method.
Background technique
When mine water during transportation, when aqueduct leaks, it will break the pressure balance in water-supply-pipe, manage Water in road is sprayed from breakthrough under pressure, and water and the breakthrough and the soil body of pipeline etc. of ejection touch Hit, to generate the vibration of different form and frequency, while acoustic signals will along pipeline upstream with downstream travel.Wherein, high Frequency signal is by rapid decay, and low frequency signal can propagate larger distance, the sound source of generation generally leak leaked at mouth noise, High-pressure water head strike note, dielectric friction sound.Wherein, when high pressure water flow strike note reaches ground, greater attenuation, operator are generated The reasons such as subjectivity, locating accuracy are lower;The frequency of dielectric friction sound is lower, and ambient enviroment interferes larger, such letter Number more difficult detection;Noise at breakthrough mainly passes through pipe wall and manages interior medium to be propagated along pipeline two sides.Meanwhile mine water There are kinds of risks during utilization: if northern area is larger with respect to day and night temperature, conveyance conduit is under soil hot and cold alternation It is also easy to produce piping failure, causes conveying water loss, water supply company is actually arrived according to measuring equipment verification respectively with water enterprise Factory's water, as there are pipe leakage phenomenons in transmission process, it will be related to enterprise and water supply company's reference water cost etc. and ask Topic causes quality of water conveyance to change, will increase enterprise's water treatment difficulty if pipe leakage point harmful substance enters, and increases Water treatment cost.
Domestic and foreign scholars are made that a large amount of thoroughgoing and painstaking researchs to aqueduct leakage detection method, propose many realities With effective theory and method.Mainly there are passive detection method, sound to listen detection method, correlation analysis detection method, probe method etc. in pipe.
(1) passive detection method.Mainly to find bright leakage.When there is water leakage situation generation on ground, immediately to aqueduct It carries out disease inspection and repairs.Related personnel can be arranged constantly to check discovery leakage problems, can also be leaked by user report Phenomenon.This method has the advantages that small investment, but concentrates on the inspection to leakage losses problems such as bright leakages, it is generally the case that discovery is seeped A large amount of water resource waste is had resulted in when leakage.Therefore, this method can be applied to the pipeline in underground environment, and surrounding is without river Deng.
(2) sound listens detection method.The method can carry out audition identification to leakage noise by bibcock and road surface.It further include probing Positioning.Former can be applied to be pre-positioned pipe leakage position, clue and range including leak etc.;Latter two is used for The specific location that leakage occurs is found out, referred to as leakage is accurately positioned.According to Britain's Water Research Center Research statistics, sound listens detection method The underground leak of Yue Kecha 80%, while sound listens the ability quality of the examined personnel of detection method, sound to listen the property and management of equipment Work etc. is affected.
(3) correlation analysis detection method.It may filter that interference signal by relevant treatment, therefore correlator detection leakage also can It is enough to be carried out on daytime.But the leak detection effect of Leak Noise Correlator will receive the influence of many factors, material, interface shape such as pipeline Formula, hydraulic pressure height, surrounding thickness of earth covering etc..It is more difficult complete using conventional detection method when Pipeline Crossing Program building or river bed At detection, however Leak Noise Correlator has apparent superiority in this regard, can guarantee to work normally.According to the use of Britain Experience, Leak Noise Correlator can be applied to the small-bore aqueduct of general water leakage situation, and the precision that testing result has can be extremely 90%;Meanwhile 200m should not be lower than by the distance arrangement of two sensors too far.
(4) probe method in managing.This detection method mainly develops on the basis of the technologies such as magnetic flux, ultrasound, vortex, video recording ?.Detector is placed in aqueduct and has signal device, can be moved in aqueduct with medium in pipe It is dynamic, breakthrough is found out according to the position of the electric signal and generation that leak spot noise.Simultaneously, additionally it is possible to utilize radioactive tracer Method, when conveyance conduit leaks, tracer will be diffused into the surrounding medium outside pipe, at this time can be in no tracer In the case of using specific probe incude the survivor of the position, and then find leak location.The shortcomings that such method is Detection lacks continuity, can only discontinuity detection, and block or other accident will seriously affect detection work and be normally carried out.
Summary of the invention
Aiming at the problems existing in the prior art, the present invention provides a kind of leakage orienting for merging hydrophone and Wavelet Denoising Method Method improves the accuracy and science of leakage orienting, improves the utilization rate of mine water resource.
To achieve the goals above, the following technical solution is employed by the present invention:
A kind of leakage orienting method merging hydrophone and Wavelet Denoising Method provided by the invention comprising the steps of:
Step 1, one hydrophone is respectively installed in the two ends of aqueduct, using hydrophone in mine water transmission process The information of acoustic wave of generation is obtained, and untreated time series original acoustic wave information source is obtained;
Step 2, wavelet coefficient noise reduction process model is established, noise-removed filtering is carried out to the original non-stationary signal source of acquisition;
Step 3, setup delay value model determines time delay value to the signal after denoising;
Step 4, location model is established, the time delay value obtained using step 3 positions breakthrough.
Further, the step 2 the specific implementation process is as follows:
Step 201, wavelet transformation is done to signals and associated noises, obtains one group of wavelet coefficient w, according to the linear property of wavelet transformation Matter, wavelet coefficient w are formed by stacking by signal wavelet coefficient and noise wavelet coefficients;
Step 202, threshold process is carried out to wavelet coefficient w, obtains estimation wavelet coefficient w', formula is as follows:
Wherein, sgn (w) is Semisoft threshold function table, λ1For the minimum threshold of threshold function table, λ2Most for threshold function table Big threshold value;
Step 203, wavelet reconstruction is carried out using estimation wavelet coefficient w', obtaining estimation signal is the signal x after denoising.
Further, the step 3 the specific implementation process is as follows:
Step 301, time delay value and Embedded dimensions are primarily determined according to the sigtnal interval after denoising, time delay value is greater than letter here Number received time interval, Embedded dimensions >=2;
Step 302, correlation integral is established according to C-C method, forms corresponding relationship, the Embedded dimensions of time delay value and correlation integral With the corresponding relationship of correlation integral;The correlation integral of time delay value, Embedded dimensions are as follows:
dij=| | Xi-Xj||
If x < 0, θ (x)=0;If x >=0, θ (x)=1;
Wherein, τ is time delay value, and m is Embedded dimensions, and x is the point of phase space, and r is the radius of phase space, dijFor in phase space The distance between any two points i, j, n are signal length;
Step 303, time delay value and Embedded dimensions are determined;
Correlation integral is a cumulative distribution function, indicates that distance is less than the probability of r between any two points in phase space, here The Infinite Norm of the difference of the vector of distance between points indicates, defines test statistics:
S1(m, n, r, τ)=C (m, n, r, τ)-Cm(1,n,r,τ)
The test statistics of the definition strategy average using piecemeal:
ΔS2(m, τ)=max { S2(m,rj,τ)-min{S2(m,rj,τ)}
ΔS2(m, τ) has measured S2(m,rj, τ) and for~τ to the maximum deviation of all radius r, optimal time delay value is Δ S2(m, First local minimum of τ)~τ;
In formula, ri=i*0.5 σ, σ=std (x), σ are the standard deviation of time series;
Comprehensively considerWithDefine index:
S2corThe global minimum of (τ) is Embedded dimensions m.
Further, the step 4 the specific implementation process is as follows:
t1+t2=S/V
t1-t2=τ, and t1> t2
L=V*t1
Leakage is obtained by above formula, distance of the point apart from upstream hydrophone occurs, that is, positions the calculation formula of breakthrough:
L=(S+V* τ)/2
Wherein, t1The time of the information of acoustic wave of breakthrough, t are measured for upstream hydrophone2Breakthrough is measured for downstream hydrophone Information of acoustic wave time, L is that distance of the point apart from upstream hydrophone occurs for leakage, and S is the spacing of two hydrophones, and V is letter Number along pipe transmmision speed, τ is time delay value, and k is the bulk modulus of water, and ρ is the density of water, and C is description pipeline constraint The dimensionless group influenced on acoustic wave propagation velocity is acted on, taking 1, D for water supply line is internal diameter of the pipeline, and E is pipe material Young's modulus of elasticity, e are pipe thickness.
Compared with prior art, the invention has the following advantages that
1, the invention discloses a kind of leakage orienting methods for merging hydrophone and Wavelet Denoising Method, it is intended to unconventional water source (mine water) is handled using transmission process information of acoustic wave process optimization, judges aqueduct leak location.This method is listened using water Device obtains water and conveys information of acoustic wave, at any time from the complexity in metric source, uncertainty, disordering angular quantification information of acoustic wave Between dynamic changing process uncertain factor is significantly cut down to water signal to original non-stationary information of acoustic wave noise-removed filtering It influences, setup delay value model improves time delay accuracy, improves the accuracy and science of leakage orienting, improve mine The utilization rate of well water resource.
2, the foundation of wavelet coefficient noise reduction process model provides signal at some to non-stationary signal source noise reduction process Between point on situation of change, highlight the variation characteristic of actual signal in original signal dominant frequency range, reduce objective factor pair The interference of actual signal.
3, the foundation of time delay value model, solve leakage acoustic signals dominant frequency components receive leakage acoustic propagation it is caused Time delay effect, determine breakthrough signal transmission in delay time.
4, the foundation of location model, by pipe leakage point two sides hydrophone received signal, quick, precise positioning leakage Point.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with It obtains other drawings based on these drawings.
Fig. 1 is the flow chart of a kind of leakage orienting method of of the invention fusion hydrophone and Wavelet Denoising Method;
Fig. 2 is the signal graph after original signal and denoising;
Fig. 3 is the waveform diagram of time delay value Yu correlation integral corresponding relationship;
Fig. 4 is the waveform diagram of Embedded dimensions Yu correlation integral corresponding relationship;
Fig. 5 is breakthrough positioning schematic.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
Core of the invention is to provide a kind of leakage orienting method for merging hydrophone and Wavelet Denoising Method, and it is fixed to improve leakage The accuracy of position.
Present invention will be described in further detail below with reference to the accompanying drawings and specific embodiments.
Embodiment one
As shown in Figure 1, the present embodiment provides a kind of leakage orienting methods for merging hydrophone and Wavelet Denoising Method, comprising following Step:
Hydrophone is mounted on the two ends of aqueduct by step S101, operator, defeated to mine water using hydrophone The information of acoustic wave generated during sending is obtained, and untreated time series original acoustic wave information source is obtained;
Step S102 establishes wavelet coefficient noise reduction process model, carries out denoising filter to the original non-stationary signal source of acquisition Wave;
The characteristics of small echo (Wavelet) analysis method be low frequency part frequency resolution with higher and it is lower when Between resolution ratio be suitable for detecting normal signal in high frequency section temporal resolution with higher and lower frequency resolution The transient state abnormal phenomena of middle entrainment simultaneously shows its ingredient.Wavelet Denoising Method is exactly the priori knowledge using particular problem, according to signal With the wavelet coefficient of noise there is mechanism of different nature on different scale, construct respective rule, uses other in wavelet field Mathematical method handles the wavelet coefficient of signals and associated noises.In practice, useful signal is usually expressed as low frequency signal or relatively flat Steady signal, noise signal then show as high-frequency signal.Denoising process specifically:
Step S1021 does wavelet transformation to signals and associated noises, obtains one group of wavelet coefficient w, according to the linear property of wavelet transformation Matter, wavelet coefficient w are formed by stacking by signal wavelet coefficient and noise wavelet coefficients;
Step S1022 carries out threshold process to wavelet coefficient w, obtains estimation wavelet coefficient w', so that both w' and w Difference is as small as possible, and formula is as follows:
Wherein, sgn (w) is Semisoft threshold function table, λ1For the minimum threshold of threshold function table, λ2Most for threshold function table Big threshold value;
Step S1023 carries out wavelet reconstruction using estimation wavelet coefficient w', and obtaining estimation signal is the signal after denoising X, as shown in Figure 2.
Step S103, setup delay value model determine time delay value to the signal after denoising;Specific implementation process are as follows:
Step S1031 primarily determines time delay value and Embedded dimensions according to the sigtnal interval after denoising, and time delay value is greater than here The received time interval of signal, Embedded dimensions >=2;
Step S1032 establishes correlation integral according to C-C method, forms corresponding relationship, the insertion dimension of time delay value and correlation integral Several and correlation integral corresponding relationship, as shown in Figure 3 and Figure 4;
C-C method is the effective ways of determining chaos time sequence time delay value and Embedded dimensions, and C-C method first defines association product Point, statistic is reconstructed, the appropriate value according to distance between any two points in the determining insertion dimension of BDO statistical conclusions, phase space Range realizes the statistical relationship with the time, determines the estimated value of time delay value τ and Embedded dimensions m.
The correlation integral of time delay value, Embedded dimensions are as follows:
dij=| | Xi-Xj||
If x < 0, θ (x)=0;If x >=0, θ (x)=1;
Wherein, τ is time delay value, and m is Embedded dimensions, and x is the point of phase space, and r is the radius of phase space, dijFor in phase space The distance between any two points i, j, n are signal length;
Step S1033, determines time delay value and Embedded dimensions;
Correlation integral is a cumulative distribution function, indicates that distance is less than the probability of r between any two points in phase space, here The Infinite Norm of the difference of the vector of distance between points indicates, defines test statistics:
S1(m, n, r, τ)=C (m, n, r, τ)-Cm(1,n,r,τ)
The test statistics of the definition strategy average using piecemeal:
ΔS2(m, τ)=max { S2(m,rj,τ)-min{S2(m,rj,τ)}
ΔS2(m, τ) has measured S2(m,rj, τ) and for~τ to the maximum deviation of all radius r, optimal time delay value is Δ S2(m, First local minimum of τ)~τ;
In formula, ri=i*0.5 σ, σ=std (x), σ are the standard deviation of time series;
Comprehensively considerWithDefine index:
S2corThe global minimum of (τ) is Embedded dimensions m.
Step S104, establishes location model, and the time delay value obtained using step S103 positions breakthrough.
When aqueduct leaks, it will break the pressure balance in aqueduct, fluid will be flowed out out of pipeline, And in duct wall generate rub so that generate leakage noise, while acoustic signals will along pipeline upstream with downstream travel.Its In, high-frequency signal is by rapid decay, and low frequency signal can propagate larger distance, there is such as ShiShimonoseki between each variable in pipeline System, as shown in Figure 5:
t1+t2=S/V
t1-t2=τ, and t1> t2
L=V*t1
Leakage is derived by by above formula, distance calculation formula of the point apart from upstream hydrophone, i.e. positioning breakthrough occurs Calculation formula:
L=(S+V* τ)/2
Wherein, t1The time of the information of acoustic wave of breakthrough, t are measured for upstream hydrophone2Breakthrough is measured for downstream hydrophone Information of acoustic wave time, L is that distance of the point apart from upstream hydrophone, the i.e. distance at the end distance A occur for leakage, and S is two water The spacing of device is listened, V is speed of the signal along pipe transmmision, and τ is time delay value, and k is the bulk modulus of water, and ρ is the density of water, C is the dimensionless group for describing pipeline effect of contraction and influencing on acoustic wave propagation velocity, and taking 1, D for water supply line is in pipeline Diameter, E are the Young's modulus of elasticity of pipe material, and e is pipe thickness.
L is calculated according to above-mentioned formula, according to upstream hydrophone initial set up position, pushes away L distance toward downstream direction, i.e., For breakthrough.
The foundation of three independent models is coupled into pipe leakage location model, reduces objective factor to actual signal Disturbance improves the science and accuracy of leakage point location.Meanwhile calculation method of the invention is without complicated programming, benefit With MATLAB software related command routine call correlation module, low to personnel's technical requirements, easy to accomplish, generalization is good.
Embodiment two
The present embodiment provides a specific leakage orienting test, detailed processes are as follows:
Step S201 prepares two BA-FH3200 type fibre optic hydrophones and its corollary equipment, hose, force (forcing) pump etc., In, hose length 400m, diameter DN25mm, wall thickness 2mm.
Fibre optic hydrophone is connect by step S202 with data acquisition equipment, and line delay of going forward side by side calibration, calibration, demodulation etc. are protected Equipment is demonstrate,proved to operate normally.
Aqueduct upstream is connected water pump by step S203, provides water source by a water pot, and pressurize using force (forcing) pump, pressure It is by force 0.5MPa;Aqueduct downstream is connected with collection device and (for collecting test water, reduces water resource waste), and formed It floods and flows out.A breakthrough L is arranged on aqueduct1, leakage area 1cm2
Step S204, by two hydrophone (S1、S2) external is used to be individually positioned in breakthrough L1Two sides.At three kinds Under operating condition, tested respectively.
Operating condition 1, S1In upstream side, from breakthrough L1Distance L1S1For 120m;S2In downstream side, from breakthrough L1Distance L1S2 For 90m;
Operating condition 2, S1In upstream side, from breakthrough L1Distance L1S1For 130m;S2In downstream side, from breakthrough L1Distance L1S2 For 100m;
Operating condition 3, S1In upstream side, from breakthrough L1Distance L1S1For 140m;S2In downstream side, from breakthrough L1Distance L1S2 For 110m.
Step S205, Open valve supply water, and guarantee that ambient enviroment is quiet, the initial sound letter that observation fibre optic hydrophone is collected into Number wave character, and carry out unit comparison.When signal stabilization, stores signal and be copied to Wavelet-based Denoising software journey Sequence, drop wave denoising as shown in figure 3, determining time delay value using local minimum, and calculate breakthrough position according to time delay value.Three The results are shown in Table 1 for the analysis of kind operating condition:
Operating condition L1S1Theoretical value/m L1S2Theoretical value/m L1S1Measured value/m Error
1 120 90 117.3 2.25%
2 130 100 126.7 2.53%
3 140 110 134.66 3.1%
Table 1
The leakage orienting method of fusion hydrophone and Wavelet Denoising Method through the invention, it is theoretical as can be seen from the test results The deviation of value and measured value, with the increase of distance, error increased in a certain range, but amplitude is smaller, keep In 5% error range, precise positioning may be implemented.
Finally, it should be noted that embodiment described above, only a specific embodiment of the invention, to illustrate the present invention Technical solution, rather than its limitations, scope of protection of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair It is bright to be described in detail, those skilled in the art should understand that: anyone skilled in the art In the technical scope disclosed by the present invention, it can still modify to technical solution documented by previous embodiment or can be light It is readily conceivable that variation or equivalent replacement of some of the technical features;And these modifications, variation or replacement, do not make The essence of corresponding technical solution is detached from the spirit and scope of technical solution of the embodiment of the present invention, should all cover in protection of the invention Within the scope of.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (4)

1. a kind of leakage orienting method for merging hydrophone and Wavelet Denoising Method, which is characterized in that comprise the steps of:
Step 1, one hydrophone is respectively installed in the two ends of aqueduct, using hydrophone to generating in mine water transmission process Information of acoustic wave obtained, obtain untreated time series original acoustic wave information source;
Step 2, wavelet coefficient noise reduction process model is established, noise-removed filtering is carried out to the original non-stationary signal source of acquisition;
Step 3, setup delay value model determines time delay value to the signal after denoising;
Step 4, location model is established, the time delay value obtained using step 3 positions breakthrough.
2. the leakage orienting method of fusion hydrophone and Wavelet Denoising Method according to claim 1, which is characterized in that the step Rapid 2 the specific implementation process is as follows:
Step 201, wavelet transformation is done to signals and associated noises, obtains one group of wavelet coefficient w, it is small according to the linear behavio(u)r of wavelet transformation Wave system number w is formed by stacking by signal wavelet coefficient and noise wavelet coefficients;
Step 202, threshold process is carried out to wavelet coefficient w, obtains estimation wavelet coefficient w', formula is as follows:
Wherein, sgn (w) is Semisoft threshold function table, λ1For the minimum threshold of threshold function table, λ2For the maximum threshold of threshold function table Value;
Step 203, wavelet reconstruction is carried out using estimation wavelet coefficient w', obtaining estimation signal is the signal x after denoising.
3. the leakage orienting method of fusion hydrophone and Wavelet Denoising Method according to claim 2, which is characterized in that the step Rapid 3 the specific implementation process is as follows:
Step 301, time delay value and Embedded dimensions are primarily determined according to the sigtnal interval after denoising, time delay value connects greater than signal here The time interval of receipts, Embedded dimensions >=2;
Step 302, correlation integral is established according to C-C method, forms corresponding relationship, Embedded dimensions and the pass of time delay value and correlation integral Join the corresponding relationship of integral;The correlation integral of time delay value, Embedded dimensions are as follows:
dij=| | Xi-Xj||
If x < 0, θ (x)=0;If x >=0, θ (x)=1;
Wherein, τ is time delay value, and m is Embedded dimensions, and x is the point of phase space, and r is the radius of phase space, dijIt is any in phase space The distance between two o'clock i, j, n are signal length;
Step 303, time delay value and Embedded dimensions are determined;
Correlation integral is cumulative distribution function, indicates that distance is less than the probability of r between any two points in phase space, here point with The distance between point is indicated with the Infinite Norm of the difference of vector, defines test statistics:
S1(m, n, r, τ)=C (m, n, r, τ)-Cm(1,n,r,τ)
The test statistics of the definition strategy average using piecemeal:
ΔS2(m, τ)=max { S2(m,rj,τ)-min{S2(m,rj,τ)}
ΔS2(m, τ) has measured S2(m,rj, τ) and for~τ to the maximum deviation of all radius r, optimal time delay value is Δ S2(m, τ)~τ First local minimum;
In formula, ri=i*0.5 σ, σ=std (x), σ are the standard deviation of time series;
Comprehensively considerWithDefine index:
S2corThe global minimum of (τ) is Embedded dimensions m.
4. the leakage orienting method of fusion hydrophone and Wavelet Denoising Method according to claim 3, which is characterized in that the step Rapid 4 the specific implementation process is as follows:
t1+t2=S/V
t1-t2=τ, and t1> t2
L=V*t1
Leakage is derived by by above formula, distance calculation formula of the point apart from upstream hydrophone occurs, that is, positions the meter of breakthrough Calculate formula:
L=(S+V* τ)/2
Wherein, t1The time of the information of acoustic wave of breakthrough, t are measured for upstream hydrophone2The sound of breakthrough is measured for downstream hydrophone The time of wave information, L are that distance of the point apart from upstream hydrophone occurs for leakage, and S is the spacing of two hydrophones, and V is signal edge The speed of pipe transmmision, τ are time delay value, and k is the bulk modulus of water, and ρ is the density of water, and C is description pipeline effect of contraction On the dimensionless group that acoustic wave propagation velocity influences, taking 1, D for water supply line is internal diameter of the pipeline, and E is the Young of pipe material Elasticity modulus, e are pipe thickness.
CN201811186445.8A 2018-10-12 2018-10-12 A kind of leakage orienting method merging hydrophone and Wavelet Denoising Method Pending CN109344784A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811186445.8A CN109344784A (en) 2018-10-12 2018-10-12 A kind of leakage orienting method merging hydrophone and Wavelet Denoising Method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811186445.8A CN109344784A (en) 2018-10-12 2018-10-12 A kind of leakage orienting method merging hydrophone and Wavelet Denoising Method

Publications (1)

Publication Number Publication Date
CN109344784A true CN109344784A (en) 2019-02-15

Family

ID=65309705

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811186445.8A Pending CN109344784A (en) 2018-10-12 2018-10-12 A kind of leakage orienting method merging hydrophone and Wavelet Denoising Method

Country Status (1)

Country Link
CN (1) CN109344784A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111122068A (en) * 2019-12-05 2020-05-08 郑州大学 Water pipeline leakage position positioning method based on optical fiber hydrophone
CN111899749A (en) * 2020-07-14 2020-11-06 上海建工集团股份有限公司 Noise reduction method for monitoring operation sound of concrete pumping pipeline
CN114719198A (en) * 2022-06-09 2022-07-08 武汉易维科技股份有限公司 Intelligent pipeline leakage detection method and device based on noise penetration integral domain modeling
CN116151049A (en) * 2023-04-23 2023-05-23 合肥众安睿博智能科技有限公司 Intelligent inertial navigation distance measurement management system based on ultrasonic pre-calibration data

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101196872A (en) * 2007-11-19 2008-06-11 清华大学 Leakage detecting and locating method based on pressure and sound wave information amalgamation

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101196872A (en) * 2007-11-19 2008-06-11 清华大学 Leakage detecting and locating method based on pressure and sound wave information amalgamation

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
吴华稳: "混沌时间序列分析及在铁路货运量预测中的应用研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 *
张鑫: "城市供水管网声波检漏及定位技术研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *
焦李成等: "《自然计算、机器学习与图像理解前沿》", 31 August 2008 *
王世礼等: "基于新阈值函数的小波阈值去噪算法", 《传感器与微系统》 *
相瑞等: "小波变换与偏微分方程的图像去噪算法", 《通信技术》 *
胡娟: "基于小波变换和中值滤波的图像去噪方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
褚选选: "基于噪声的长距离输水管道渗漏监测研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111122068A (en) * 2019-12-05 2020-05-08 郑州大学 Water pipeline leakage position positioning method based on optical fiber hydrophone
CN111899749A (en) * 2020-07-14 2020-11-06 上海建工集团股份有限公司 Noise reduction method for monitoring operation sound of concrete pumping pipeline
CN111899749B (en) * 2020-07-14 2023-08-29 上海建工集团股份有限公司 Noise reduction method for monitoring operation sound of concrete pumping pipeline
CN114719198A (en) * 2022-06-09 2022-07-08 武汉易维科技股份有限公司 Intelligent pipeline leakage detection method and device based on noise penetration integral domain modeling
CN116151049A (en) * 2023-04-23 2023-05-23 合肥众安睿博智能科技有限公司 Intelligent inertial navigation distance measurement management system based on ultrasonic pre-calibration data

Similar Documents

Publication Publication Date Title
CN109344784A (en) A kind of leakage orienting method merging hydrophone and Wavelet Denoising Method
Cui-wei et al. Experimental study on acoustic propagation-characteristics-based leak location method for natural gas pipelines
Khulief et al. Acoustic detection of leaks in water pipelines using measurements inside pipe
CN104595729B (en) A kind of oil and gas pipeline leakage localization method based on magnitudes of acoustic waves
CN109357171B (en) Underground pipeline leakage monitoring and positioning method and device
CN101592288B (en) Method for identifying pipeline leakage
CN100456010C (en) Method for detecting leakage of oil gas pipe based on pressure signal knee
Liu et al. Leakage monitoring research and design for natural gas pipelines based on dynamic pressure waves
CN104373821A (en) Natural gas pipeline safety monitoring device based on acoustical science active spurring
Bakhtawar et al. Review of water leak detection and localization methods through hydrophone technology
Zhang et al. Acoustic method of high-pressure natural gas pipelines leakage detection: Numerical and applications
CN103062628A (en) Urban gas buried pipeline leakage detection and locating method and applications thereof
CN106289121A (en) A kind of computational methods of reducer pipe equivalence pipe range
CN103644965B (en) A kind of sound lining acoustic impedance measurement method under the conditions of sound wave glancing incidence
Cui et al. Variable step normalized LMS adaptive filter for leak localization in water-filled plastic pipes
Zeng et al. Inverse wave reflectometry method for hydraulic transient-based pipeline condition assessment
Han et al. Localization of CO2 gas leakages through acoustic emission multi-sensor fusion based on wavelet-RBFN modeling
Liu et al. Leak location study of water pipeline based on CEEMDAN-CC at low SNR
Ozevin et al. Reliable monitoring of leak in gas pipelines using acoustic emission method
Hao et al. Leakage signal analysis of urban gas pipeline based on improved variational mode decomposition
CN112525437A (en) Underwater identification method for leakage noise of large-scale water delivery building
Wu et al. Novel signal denoising approach for acoustic leak detection
Travers Acoustic monitoring of prestressed concrete pipe
Hao et al. Research on gas pipeline multi-point leak signal processing and source locating using VMD, BSS and relative entropy
Gong et al. Distributed deterioration detection and location in single pipes using the impulse response function

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190215

RJ01 Rejection of invention patent application after publication