CN105927863A - DMA zone pipe network leakage online detecting and positioning system and detecting and positioning method thereof - Google Patents

DMA zone pipe network leakage online detecting and positioning system and detecting and positioning method thereof Download PDF

Info

Publication number
CN105927863A
CN105927863A CN201610302392.6A CN201610302392A CN105927863A CN 105927863 A CN105927863 A CN 105927863A CN 201610302392 A CN201610302392 A CN 201610302392A CN 105927863 A CN105927863 A CN 105927863A
Authority
CN
China
Prior art keywords
data
pressure
time
time period
flows
Prior art date
Application number
CN201610302392.6A
Other languages
Chinese (zh)
Other versions
CN105927863B (en
Inventor
王孝良
刘颖
Original Assignee
大连理工大学
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 大连理工大学 filed Critical 大连理工大学
Priority to CN201610302392.6A priority Critical patent/CN105927863B/en
Publication of CN105927863A publication Critical patent/CN105927863A/en
Application granted granted Critical
Publication of CN105927863B publication Critical patent/CN105927863B/en

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17DPIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00Protection or supervision of installations
    • F17D5/02Preventing, monitoring, or locating loss
    • F17D5/06Preventing, monitoring, or locating loss using electric or acoustic means

Abstract

The invention relates to the technical field of tap water supply pipe network leakage measurement, in particular to a DMA zone pipe network leakage online detecting and positioning system and a detecting and positioning method thereof. The system comprises a DMA zone tap water main pipeline, a first DMA zone tap water branch pipeline, a second DMA zone tap water branch pipeline,..., an N-1<th> DMA zone tap water branch pipeline, an N<th> DMA zone tap water branch pipeline, a data collecting unit, pressure detecting units, sound wave detecting units, a system server and a work station. The data collecting unit is installed on the tap water main pipeline in the position of an inlet of a DMA zone, and data are uploaded to the system server through a GPRS/3G/4G module. The pressure detecting units or the sound wave detecting units, or the pressure detecting units and the sound wave detecting units are installed on the first DMA zone tap water branch pipeline, the second DMA zone tap water branch pipeline,..., the N-1<th> DMA zone tap water branch pipeline and the N<th> DMA zone tap water branch pipeline correspondingly, and the pressure detecting units and the sound wave detecting units upload data to the system server through the GPRS/3G/4G module. The system server communicates with the work station through the industrial Ethernet. The DMA zone pipe network leakage online detecting and positioning system and the detecting and positioning method thereof are suitable for underground water supply pipe networks, installation is easy, and leakage detecting and positioning of terminal pipeline networks of all the urban DMA zones are well achieved.

Description

DMA subregion pipeline network leak on-line checking alignment system and detection and location method thereof
Technical field
The present invention relates to a kind of DMA subregion pipeline network leak on-line checking alignment system and detection and location method thereof, Belong to tap water supply pipe network to leak hunting technical field.
Background technology
It is that warp improves in urban water supply company that public supply mains are revealed " the non-income water rate " caused high The major obstacles of Ji benefit.Public supply mains leakage is mostly derived from the magnanimity of the little leak point of water supply network end Reveal accumulation.If the most accurately detecting, positioning the leak point of public supply mains, it will make to water resource Become waste greatly.
Along with the development of modern city, the optimization of public supply mains is transformed with decoupling, is carried out reasonable DMA Subregion is common recognition and the development trend of Urban water supply.Water system pipe network independent measure subregion (DMA) Reasonably optimizing construction and upgrading are to simplify pipe network topology, balance pipe network load, raising water supply security and strengthening to supply Water system monitoring and the only way which must be passed of management level.
Summary of the invention
For the problem overcoming prior art to exist, the present invention proposes a kind of DMA subregion pipeline network leak and examines online Surveying alignment system and detection and location method thereof, this system is suitable for groundwater supply pipe network environment, is easily installed, becomes This is cheap, solves leak diagnostics and the orientation problem of public supply mains DMA subregion end pipe network well, Leakage point can be detected in time, decrease a large amount of losses of water resource.
In order to realize foregoing invention purpose, the technical solution used in the present invention is: a kind of DMA subregion pipe network is let out Leakage on-line checking alignment system, including DMA subregion tap water main line, DMA subregion the 1st, 2 ... N-1, N tap water bye-pass, data acquisition unit, pressure sensing cell, sonic detection unit, system server and Work station, described data acquisition unit, including pressure transducer, flow transducer and data acquisition RTU module, Wherein pressure transducer, flow transducer are arranged on the tap water main line of DMA subregion porch, data Gather RTU module be arranged on the borehole wall of adjacent valve well, described pressure transducer, flow transducer respectively with In data acquisition RTU module, two RS-485 interfaces are connected, for collecting DMA subregion porch Pressure on tap water main line, data on flows are by being built in the GPRS/3G/4G of data acquisition RTU module Module by timing or real-time in the way of data are uploaded to system server;Described DMA subregion the 1st, 2 ... N-1, On N tap water bye-pass, it is separately installed with pressure sensing cell or sonic detection unit or pressure according to field working conditions Power detector unit and sonic detection unit, described pressure sensing cell, including being arranged on DMA subregion tap water Pressure transducer on bye-pass and the data acquisition RTU module being arranged on the borehole wall of adjacent valve well, described Sonic detection unit, including the pipeline sonic transducer being arranged on DMA subregion tap water bye-pass and be arranged on Data acquisition RTU module on the borehole wall of adjacent valve well, described pressure transducer, pipeline sonic transducer are respectively It is connected with RS-485 interface in data acquisition RTU module, for the tap water that will collect in DMA subregion Pressure on bye-pass and sonic data by be built in the GPRS/3G/4G module of data acquisition RTU module with Data are uploaded to system server by timing or real-time mode, by industry between system server and work station Ethernet carries out communication;Described work station, including data collection control unit, statistical model parameter calculation unit, Leak hunting positioning unit and DMA partition monitor of pipe network diagnosis unit, sound characteristics computational analysis unit, leakage point shows Show unit, wherein, described data collection control unit, also include gathering target configuration subelement, data acquisition Pattern configurations subelement, data transmit-receive and communications protocol resolve subelement, real-time database data management subelement, go through History bank interface subelement and located in connection data acquisition Collaborative Control subelement, based on data acquisition unit, pressure The data that detector unit and sonic detection unit gather, through statistical model parameter calculation unit and pipe network diagnosis unit Jointly completing judgement and the judgement in leakage orientation that whether pipeline network leak exists, sound characteristics computational analysis unit completes Judged that the pipeline sonic data leaking around orientation calculates.
The detection and location method of described a kind of DMA subregion pipeline network leak on-line checking alignment system, including following Step:
Step 1, set up DMA water rationing pipe network operation state rule statistical models, described statistical models Setting up premise is that pipeline is properly functioning and No leakage, according to actual condition, is divided into n time period by 24 hours, Have 1 flow monitoring point, p pressure monitoring point in DMA subregion, specifically include following sub-step:
A in () single time period, in DMA subregion entry data collecting unit, flow transducer gathers m DMA Subregion tap water main line data on flows, gathers 1 data on flows every time, and each data on flows all includes this stream The acquisition time of amount data, finds the maximum max of data on flows in this time period1And minima min1, peak-to-peak Value P1It is calculated by formula (1),
P1=max1-min1(1)
The meansigma methods of DMA subregion tap water main line data on flows in this time periodCalculated by formula (2) Obtain,
X &OverBar; 1 = 1 m &Sigma; i = 1 m X i - - - ( 2 )
In formula, the DMA subregion tap water total number of main line data on flows that m gathers in representing this time period, Xi Represent the i-th data in m DMA subregion tap water main line data on flows, single time in this time period The variance yields s of DMA subregion tap water main line data on flows in section1 2It is calculated by formula (3),
s 1 2 = 1 m - 1 &Sigma; i = 1 m ( X i - X &OverBar; 1 ) 2 - - - ( 3 )
M-discharge model when () sets up b, duplicon step (a) n time, in calculating respectively one day during n difference Between DMA subregion tap water main line data on flows peak-to-peak value, meansigma methods and variance yields, first record the 1st in section DMA subregion tap water main line data on flows peak-to-peak value, meansigma methods and variance yields in the individual time period, and with stream Amount data and the corresponding relation of this data acquisition time, with the time as transverse axis, in this time period, DMA subregion enters The data on flows of mouth data acquisition unit flow transducer is that the longitudinal axis sets up curve, i.e. sets up model Q1, then with DMA subregion tap water main line data on flows peak-to-peak value, meansigma methods and variance yields in 2nd time period, and With the corresponding relation of data on flows Yu this data acquisition time, with the time as transverse axis, in this time period, DMA divides The data on flows of district's entry data collecting unit flow transducer is that the longitudinal axis sets up curve, i.e. sets up model Q2, depend on This analogizes DMA subregion tap water main line data on flows peak-to-peak value, meansigma methods and side in the n-th time period of record Difference, and with the corresponding relation of data on flows Yu this data acquisition time, with the time as transverse axis, in this time period The data on flows of DMA subregion entry data collecting unit flow transducer is that the longitudinal axis sets up curve, i.e. sets up model Qn, according to the change dynamic corrections above-mentioned model Q of data1~Qn
C in () single time period, in pressure sensing cell, pressure transducer gathers k pressure at single pressure monitoring point Force data, gathers 1 pressure data every time, and each pressure data all includes the acquisition time of this pressure data, So within the single time period, single monitoring point gathers pressure data k altogether, finds number pressure in this time period According to maximum max2And minima min2, peak-to-peak value P2Obtained by formula (4),
P2=max2-min2 (4)
The meansigma methods of single pressure monitoring point pressure dataIt is calculated by formula (5),
X &OverBar; 2 = 1 k &Sigma; j = 1 k X j - - - ( 5 )
In formula, the k total number of pressure data that in representing this time period, single pressure monitoring point gathers, XjWhen representing this Between in section in jth data, single pressure monitoring point number pressure in the pressure data that gathers of single pressure monitoring point According to variance yields s2 2It is calculated by formula (6),
s 2 2 = 1 k - 1 &Sigma; j = 1 k ( X j - X &OverBar; 2 ) 2 - - - ( 6 )
D () duplicon step (c) p time, in calculating the single time period, the 1st pressure detecting point is to pth pressure The peak-to-peak value of monitoring point pressure data, meansigma methods and variance yields;
M-pressure model when () sets up e, duplicon step (d) n time, in calculating one day in the 1st time period Within the n-th time period the 1st~the peak-to-peak value of p pressure monitoring point pressure data, meansigma methods and variance yields, first 1st pressure monitoring point pressure data peak-to-peak value, meansigma methods and variance yields in first the 1st time period of record, and With the corresponding relation of pressure data Yu this pressure data acquisition time, with the time as transverse axis, in the 1st time period 1st pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up model P11;Record the 1st time period Interior 2nd pressure monitoring point pressure data peak-to-peak value, meansigma methods and variance yields, and with pressure data and this pressure The corresponding relation of data acquisition time, with the time as transverse axis, the 2nd pressure monitoring point number pressure in this time period Set up curve according to for the longitudinal axis, i.e. set up model P12;Pth pressure monitoring point pressure in the rest may be inferred this time period Data peak-to-peak value, meansigma methods and variance yields, and with the corresponding relation of pressure data Yu this pressure data acquisition time, With the time as transverse axis, in this time period, pth pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up Model P1p, in like manner, respectively in the 2nd time period of record the 1st~p pressure monitoring point pressure data peak-to-peak value, Meansigma methods and variance yields, with the corresponding relation of pressure data in the 2nd time period Yu this pressure data acquisition time, With the time as transverse axis, in this time period, pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up model P21~P2p, the rest may be inferred record respectively in the n-th time period the 1st~p pressure monitoring point pressure data peak-to-peak value, Meansigma methods and variance yields, with the corresponding relation of pressure data in the n-th time period Yu this pressure data acquisition time, With the time as transverse axis, in this time period, pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up model Pn1~Pnp.Change dynamic corrections above-mentioned model P according to data11~P1p、P21~P2pUntil Pn1~Pnp
F () sets up flow rate pressure model, in recording one day, DMA subregion tap water main line data on flows and p are individual Pressure monitoring point data, first, with different DMA subregion tap water main line datas on flows and the 1st pressure The corresponding relation of monitoring point pressure data, flow is transverse axis, and pressure is that the longitudinal axis sets up the 1st pressure monitoring point stream Amount-pressure model A1, then, with different DMA subregion tap water main line datas on flows and the 2nd pressure prison The corresponding relation of measuring point pressure data, flow is transverse axis, pressure be the longitudinal axis set up the 2nd pressure monitoring point flow- Pressure model A2, the rest may be inferred supervises with pth pressure with different DMA subregion tap water main line datas on flows The corresponding relation of measuring point pressure data, flow is transverse axis, pressure be the longitudinal axis set up pth pressure monitoring point flow- Pressure model Ap
Whether step 2, leak diagnostics unit judges exist leakage point and leakage orientation, according to practical situation, if Widow time length calculated by devise a stratagem, and widow time length is designated as t, and system of the present invention can set the threshold of each item data Value, specifically includes following sub-step:
A () sets current slot as the d time period, counted by formula in step 1 sub-step (a) (1) Calculate the peak-to-peak value of DMA subregion tap water main line data on flows in current the d time period, by step 1 DMA subregion tap water main line flow in formula (2) calculates current the d time period in sub-step (a) The meansigma methods of data, calculates DMA in current the d time period by formula in step 1 sub-step (a) (3) The variance yields of subregion tap water main line data on flows, record peak-to-peak value in the above current the d time period, Meansigma methods and variance yields, with the corresponding relation of data on flows Yu this data acquisition time, with the time as transverse axis, when In front the d time period, data on flows is that the longitudinal axis sets up curve, i.e. sets up model S, by peak-to-peak value in model S, Meansigma methods and variance yields and the time m-discharge model Q of foundation in step 1 sub-step (b)1~QnIn corresponding time Between the model Q of sectiondCompared with middle peak-to-peak value, meansigma methods are corresponding with variance yields, if the difference of peak-to-peak value is not less than being System threshold value, the difference of meansigma methods and the difference of variance yields all exceed system thresholds, then judge that this DMA subregion is deposited In leakage, if the difference of peak-to-peak value, meansigma methods and variance yields all exceedes system thresholds, then within this time period with The length of window t calculated flow rate statistical average of defaultVariance yields st 2If, this time period initial time For T1, the end time is T2, this time period is divided into the 1st window T1~T1+ t, the 2nd window t~T1 + 2t ..., x-th window T1+ xt~T2, calculated flow rate statistical average in y-th windowSide Difference st 2, meansigma methodsCalculated by formula (7),
X &OverBar; t = 1 e &Sigma; f = 1 e X f - - - ( 7 )
In formula, e represents the total number of data on flows in y-th window, XfRepresent data on flows in y-th window The f data in total number, variance yields st 2Calculated by formula (8),
s t 2 = 1 e - 1 &Sigma; f = 1 e ( X f - X &OverBar; t ) 2 - - - ( 8 )
Find the 1st meansigma methods of data on flows, variance yields discharge model Q m-with time corresponding in this time periodd Meansigma methods, variance yields compare the window more all exceeding system thresholds, record initial time t of this window1;Look for The meansigma methods of last 1 data on flows, variance yields discharge model Q m-with time corresponding within this time perioddFlat Average, variance yields compare the window more all exceeding system thresholds, record initial time t of this window2, the 1st Occur to last, abnormal window occurs that the time difference between abnormal window is designated as Δ t, Δ t by formula (9) Calculate,
Δ t=t2-t1 (9)
If Δ t exceedes the maximum time interval of default, then judge to there is leakage;
If b in judging this DMA subregion in () step 2 sub-step (a) there is leakage, then according to step in pipe network The peak-to-peak value of p pressure monitoring point, root in 1 sub-step (c), formula (4) calculates current the d time period The flat of p pressure monitoring point in current the d time period is calculated according to formula in step 1 sub-step (c) (5) Average and calculate p pressure prison in current the d time period according to formula in step 1 sub-step (c) (6) The variance yields of measuring point, if the peak-to-peak value of u pressure monitoring point, meansigma methods and variance yields one of which or two Or three exceed system correspondence threshold value, then find this pressure monitoring point in current the d time period maximum and Minima, according to the acquisition time that the maximum of this pressure monitoring point in current the d time period is corresponding, in step The model S that rapid 2 sub-steps (a) are set up finds the flow value B of correspondencemax, time simultaneously according to current the d Between the maximum of this pressure monitoring point in section, the model A set up in step 1 sub-step (f)uIn find work as The flow value C that in front the d time period, the maximum of this pressure monitoring point is correspondingmax;During according to current the d Between the minima of this pressure monitoring point is corresponding in section time point, in the model S that step 2 sub-step (a) is set up In find correspondence flow value Bmin, simultaneously according to the minima of this pressure monitoring point in current the d time period, The model A set up in step 1 sub-step (f)uIn find this pressure monitoring point in current the d time period Flow value C corresponding to minimamin, by BmaxWith CmaxCompare, BminWith CminCompare, if more than described Two all exceed system thresholds, then judge that this pressure monitoring point place pipeline section exists leakage, complete leakage point Orientation judges;
Step 3, according in step 2, the orientation of leakage point is judged, leak point positioning unit starting be installed on should The sonic transducer of orientation relevant position completes sonic data collection and calculates, according to located in connection survey calculation method Calculate leakage point away from the distance closing on sonic transducer.
The medicine have the advantages that a kind of DMA subregion pipeline network leak on-line checking alignment system, including DMA Subregion tap water main line, DMA subregion the 1st, 2 ... N-1, N tap water bye-pass, data acquisition unit, Pressure sensing cell, sonic detection unit, system server and work station, described data acquisition unit, including Pressure transducer, flow transducer and data acquisition RTU module, wherein pressure transducer, flow transducer peace Being contained on the tap water main line of DMA subregion porch, data acquisition RTU module is arranged on the well of valve well On wall, described pressure transducer, flow transducer connect with two RS-485 in data acquisition RTU module respectively Mouth is connected, for the pressure on the tap water main line collecting DMA subregion porch, data on flows being led to Cross the GPRS/3G/4G module being built in data acquisition RTU module by timing or real-time in the way of data are uploaded To system server;Described DMA subregion the 1st, 2 ... on N-1, N tap water bye-pass, according to field working conditions It is separately installed with pressure sensing cell or sonic detection unit or pressure sensing cell and sonic detection unit, described Pressure sensing cell, including the pressure transducer being arranged on DMA subregion tap water bye-pass and be arranged on neighbour Data acquisition RTU module on the borehole wall of nearly valve well, described sonic detection unit, including being arranged on DMA Pipeline sonic transducer on subregion tap water bye-pass and be arranged on the data acquisition on the borehole wall of adjacent valve well RTU module, described pressure transducer, pipeline sonic transducer respectively with RS-485 in data acquisition RTU module Interface is connected, and pressure and sonic data on the tap water bye-pass that will collect in DMA subregion pass through Be built in the GPRS/3G/4G module of data acquisition RTU module by timing or real-time in the way of data are uploaded to System server, carries out communication by EPA between system server and work station;With prior art phase Ratio, the present invention be suitable for groundwater supply pipe network, be easily installed, with low cost, solve each DMA in city well The end pipe network leak diagnostics of subregion and orientation problem, find leak point in time, processed in time, decrease A large amount of losses of water resource.
Accompanying drawing explanation
Fig. 1 is the population structure block diagram of present system.
Fig. 2 is the work station internal structure block diagram in present system.
Fig. 3 is mounted in the data acquisition unit structured flowchart on DMA subregion tap water main line.
Fig. 4 is mounted in the pressure sensing cell structured flowchart on DMA subregion pipe network bye-pass.
Fig. 5 is mounted in the sonic detection cellular construction block diagram on DMA subregion pipe network bye-pass.
Fig. 6 is mounted in the pressure sensing cell on DMA subregion pipe network bye-pass and sonic detection cellular construction Block diagram.
Fig. 7 is the inventive method flow chart of steps.
Detailed description of the invention
The invention will be further described below in conjunction with the accompanying drawings.
As shown in Fig. 1,2,3,4,5,6, a kind of DMA subregion pipeline network leak on-line checking alignment system, Including DMA subregion tap water main line, DMA subregion the 1st, 2 ... N-1, N tap water bye-pass, data Collecting unit, pressure sensing cell, sonic detection unit, system server and work station, described data acquisition Unit, including pressure transducer, flow transducer and data acquisition RTU module, wherein pressure transducer, stream Quantity sensor is arranged on the tap water main line of DMA subregion porch, and data acquisition RTU module is arranged on On the borehole wall of adjacent valve well, described pressure transducer, flow transducer respectively with in data acquisition RTU module Two RS-485 interfaces are connected, for by the pressure on the tap water main line collecting DMA subregion porch, Data on flows by be built in the GPRS/3G/4G module of data acquisition RTU module with timing or real-time by the way of Data are uploaded to system server;Described DMA subregion the 1st, 2 ... on N-1, N tap water bye-pass, root It is separately installed with pressure sensing cell or sonic detection unit or pressure sensing cell and sonic detection according to field working conditions Unit, described pressure sensing cell, including the pressure transducer being arranged on DMA subregion tap water bye-pass And the data acquisition RTU module being arranged on the borehole wall of adjacent valve well, described sonic detection unit, including peace It is contained in the pipeline sonic transducer on DMA subregion tap water bye-pass and is arranged on the borehole wall of adjacent valve well Data acquisition RTU module, described pressure transducer, pipeline sonic transducer respectively with data acquisition RTU module Upper RS-485 interface is connected, the pressure on the tap water bye-pass that will collect in DMA subregion and sound wave Data by be built in the GPRS/3G/4G module of data acquisition RTU module with timing or real-time by the way of by number According to being uploaded to system server, between system server and work station, carry out communication by EPA;Described Work station, including data collection control unit, statistical model parameter calculation unit, pipe network diagnosis unit, sound spy Property computational analysis unit, leakage point leak hunting positioning unit and DMA partition monitor display unit, wherein, described Data collection control unit, also includes gathering target configuration subelement, data acquisition scheme configuration subelement, number Subelement, real-time database data management subelement, history library interface subelement and phase is resolved according to transmitting-receiving and communications protocol Close position data collecting Collaborative Control subelement, based on data acquisition unit, pressure sensing cell and sonic detection The data that unit gathers, jointly completing pipeline network leak through statistical model parameter calculation unit and pipe network diagnosis unit is The judgement of no existence and the judgement in leakage orientation, sound characteristics computational analysis unit completes to be judged to leak around orientation Pipeline sonic data calculate.The detailed process of system work is: pipe network diagnosis unit and leakage point leak hunting location Unit is by collection target configuration subelement, the data acquisition scheme in DMA partition data acquisition controlling unit Configuration subelement completion system initial configuration, then, data acquisition unit, pressure sensing cell, sound wave are examined Survey the related data that collects of unit and be uploaded to system server.According to the data gathered, statistical model parameter Computing unit calculates water supply statistical model, and pipe network diagnosis unit diagnostic tube Running State, sonic detection unit divides Analysis pipe network acoustical signal characteristic, completes pipe network operation state and the analysis of gathered acoustical signal.Then, then by letting out The leak source positioning unit that leaks hunting is assigned measurement of correlation order and calculates leakage point position, it is achieved the detection to pipeline network leak With location, and by each for pipe network item data and to pipeline network leak detection with location result be shown in DMA subregion On monitoring display interface.
The detection and location method of described a kind of DMA subregion pipeline network leak on-line checking alignment system, including following Step:
Step 1, set up DMA water rationing pipe network operation state rule statistical models, described statistical models Setting up premise is that pipeline is properly functioning and No leakage, according to actual condition, is divided into n time period by 24 hours, Have 1 flow monitoring point, p pressure monitoring point in DMA subregion, specifically include following sub-step:
A in () single time period, in DMA subregion entry data collecting unit, flow transducer gathers m DMA Subregion tap water main line data on flows, gathers 1 data on flows every time, and each data on flows all includes this stream The acquisition time of amount data, finds the maximum max of data on flows in this time period1And minima min1, peak-to-peak Value P1It is calculated by formula (1),
P1=max1-min1 (1)
The meansigma methods of DMA subregion tap water main line data on flows in this time periodCalculated by formula (2) Obtain,
X &OverBar; 1 = 1 m &Sigma; i = 1 m X i - - - ( 2 )
In formula, the DMA subregion tap water total number of main line data on flows that m gathers in representing this time period, Xi Represent the i-th data in m DMA subregion tap water main line data on flows, single time in this time period The variance yields s of DMA subregion tap water main line data on flows in section1 2It is calculated by formula (3),
s 1 2 = 1 m - 1 &Sigma; i = 1 m ( X i - X &OverBar; 1 ) 2 - - - ( 3 )
M-discharge model when () sets up b, duplicon step (a) n time, in calculating respectively one day during n difference Between DMA subregion tap water main line data on flows peak-to-peak value, meansigma methods and variance yields, first record the 1st in section DMA subregion tap water main line data on flows peak-to-peak value, meansigma methods and variance yields in the individual time period, and with stream Amount data and the corresponding relation of this data acquisition time, with the time as transverse axis, in this time period, DMA subregion enters The data on flows of mouth data acquisition unit flow transducer is that the longitudinal axis sets up curve, i.e. sets up model Q1, then with DMA subregion tap water main line data on flows peak-to-peak value, meansigma methods and variance yields in 2nd time period, and With the corresponding relation of data on flows Yu this data acquisition time, with the time as transverse axis, in this time period, DMA divides The data on flows of district's entry data collecting unit flow transducer is that the longitudinal axis sets up curve, i.e. sets up model Q2, depend on This analogizes DMA subregion tap water main line data on flows peak-to-peak value, meansigma methods and side in the n-th time period of record Difference, and with the corresponding relation of data on flows Yu this data acquisition time, with the time as transverse axis, in this time period The data on flows of DMA subregion entry data collecting unit flow transducer is that the longitudinal axis sets up curve, i.e. sets up model Qn, according to the change dynamic corrections above-mentioned model Q of data1~Qn
C in () single time period, in pressure sensing cell, pressure transducer gathers k pressure at single pressure monitoring point Force data, gathers 1 pressure data every time, and each pressure data all includes the acquisition time of this pressure data, So within the single time period, single monitoring point gathers pressure data k altogether, finds number pressure in this time period According to maximum max2And minima min2, peak-to-peak value P2Obtained by formula (4),
P2=max2-min2 (4)
The meansigma methods of single pressure monitoring point pressure dataIt is calculated by formula (5),
X &OverBar; 2 = 1 k &Sigma; j = 1 k X j - - - ( 5 )
In formula, the k total number of pressure data that in representing this time period, single pressure monitoring point gathers, XjWhen representing this Between in section in jth data, single pressure monitoring point number pressure in the pressure data that gathers of single pressure monitoring point According to variance yields s2 2It is calculated by formula (6),
s 2 2 = 1 k - 1 &Sigma; j = 1 k ( X j - X &OverBar; 2 ) 2 - - - ( 6 )
D () duplicon step (c) p time, in calculating the single time period, the 1st pressure detecting point is to pth pressure The peak-to-peak value of monitoring point pressure data, meansigma methods and variance yields;
M-pressure model when () sets up e, duplicon step (d) n time, in calculating one day in the 1st time period Within the n-th time period the 1st~the peak-to-peak value of p pressure monitoring point pressure data, meansigma methods and variance yields, first 1st pressure monitoring point pressure data peak-to-peak value, meansigma methods and variance yields in first the 1st time period of record, and With the corresponding relation of pressure data Yu this pressure data acquisition time, with the time as transverse axis, in the 1st time period 1st pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up model P11;Record the 1st time period Interior 2nd pressure monitoring point pressure data peak-to-peak value, meansigma methods and variance yields, and with pressure data and this pressure The corresponding relation of data acquisition time, with the time as transverse axis, the 2nd pressure monitoring point number pressure in this time period Set up curve according to for the longitudinal axis, i.e. set up model P12;Pth pressure monitoring point pressure in the rest may be inferred this time period Data peak-to-peak value, meansigma methods and variance yields, and with the corresponding relation of pressure data Yu this pressure data acquisition time, With the time as transverse axis, in this time period, pth pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up Model P1p, in like manner, respectively in the 2nd time period of record the 1st~p pressure monitoring point pressure data peak-to-peak value, Meansigma methods and variance yields, with the corresponding relation of pressure data in the 2nd time period Yu this pressure data acquisition time, With the time as transverse axis, in this time period, pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up model P21~P2p, the rest may be inferred record respectively in the n-th time period the 1st~p pressure monitoring point pressure data peak-to-peak value, Meansigma methods and variance yields, with the corresponding relation of pressure data in the n-th time period Yu this pressure data acquisition time, With the time as transverse axis, in this time period, pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up model Pn1~Pnp.Change dynamic corrections above-mentioned model P according to data11~P1p、P21~P2pUntil Pn1~Pnp
F () sets up flow rate pressure model, in recording one day, DMA subregion tap water main line data on flows and p are individual Pressure monitoring point data, first, with different DMA subregion tap water main line datas on flows and the 1st pressure The corresponding relation of monitoring point pressure data, flow is transverse axis, and pressure is that the longitudinal axis sets up the 1st pressure monitoring point stream Amount-pressure model A1, then, with different DMA subregion tap water main line datas on flows and the 2nd pressure prison The corresponding relation of measuring point pressure data, flow is transverse axis, pressure be the longitudinal axis set up the 2nd pressure monitoring point flow- Pressure model A2, the rest may be inferred supervises with pth pressure with different DMA subregion tap water main line datas on flows The corresponding relation of measuring point pressure data, flow is transverse axis, pressure be the longitudinal axis set up pth pressure monitoring point flow- Pressure model Ap
Whether step 2, leak diagnostics unit judges exist leakage point and leakage orientation, according to practical situation, if Widow time length calculated by devise a stratagem, and widow time length is designated as t, and system of the present invention can set the threshold of each item data Value, specifically includes following sub-step:
A () sets current slot as the d time period, counted by formula in step 1 sub-step (a) (1) Calculate the peak-to-peak value of DMA subregion tap water main line data on flows in current the d time period, by step 1 DMA subregion tap water main line flow in formula (2) calculates current the d time period in sub-step (a) The meansigma methods of data, calculates DMA in current the d time period by formula in step 1 sub-step (a) (3) The variance yields of subregion tap water main line data on flows, record peak-to-peak value in the above current the d time period, Meansigma methods and variance yields, with the corresponding relation of data on flows Yu this data acquisition time, with the time as transverse axis, when In front the d time period, data on flows is that the longitudinal axis sets up curve, i.e. sets up model S, by peak-to-peak value in model S, Meansigma methods and variance yields and the time m-discharge model Q of foundation in step 1 sub-step (b)1~QnIn corresponding time Between the model Q of sectiondCompared with middle peak-to-peak value, meansigma methods are corresponding with variance yields, if the difference of peak-to-peak value is not less than being System threshold value, the difference of meansigma methods and the difference of variance yields all exceed system thresholds, then judge that this DMA subregion is deposited In leakage, if the difference of peak-to-peak value, meansigma methods and variance yields all exceedes system thresholds, then within this time period with The length of window t calculated flow rate statistical average of defaultVariance yields st 2If, this time period initial time For T1, the end time is T2, this time period is divided into the 1st window T1~T1+ t, the 2nd window t~T1 + 2t ..., x-th window T1+ xt~T2, calculated flow rate statistical average in y-th windowSide Difference st 2, meansigma methodsCalculated by formula (7),
X &OverBar; t = 1 e &Sigma; f = 1 e X f - - - ( 7 )
In formula, e represents the total number of data on flows in y-th window, XfRepresent data on flows in y-th window The f data in total number, variance yields st 2Calculated by formula (8),
s t 2 = 1 e - 1 &Sigma; f = 1 e ( X f - X &OverBar; t ) 2 - - - ( 8 )
Find the 1st meansigma methods of data on flows, variance yields discharge model Q m-with time corresponding in this time periodd Meansigma methods, variance yields compare the window more all exceeding system thresholds, record initial time t of this window1;Look for The meansigma methods of last 1 data on flows, variance yields discharge model Q m-with time corresponding within this time perioddFlat Average, variance yields compare the window more all exceeding system thresholds, record initial time t of this window2, the 1st Occur to last, abnormal window occurs that the time difference between abnormal window is designated as Δ t, Δ t by formula (9) Calculate,
Δ t=t2-t1 (9)
If Δ t exceedes the maximum time interval of default, then judge to there is leakage;
If b in judging this DMA subregion in () step 2 sub-step (a) there is leakage, then according to step in pipe network The peak-to-peak value of p pressure monitoring point, root in 1 sub-step (c), formula (4) calculates current the d time period The flat of p pressure monitoring point in current the d time period is calculated according to formula in step 1 sub-step (c) (5) Average and calculate p pressure prison in current the d time period according to formula in step 1 sub-step (c) (6) The variance yields of measuring point, if the peak-to-peak value of u pressure monitoring point, meansigma methods and variance yields one of which or two Or three exceed system correspondence threshold value, then find this pressure monitoring point in current the d time period maximum and Minima, according to the acquisition time that the maximum of this pressure monitoring point in current the d time period is corresponding, in step The model S that rapid 2 sub-steps (a) are set up finds the flow value B of correspondencemax, time simultaneously according to current the d Between the maximum of this pressure monitoring point in section, the model A set up in step 1 sub-step (f)uIn find work as The flow value C that in front the d time period, the maximum of this pressure monitoring point is correspondingmax;During according to current the d Between the minima of this pressure monitoring point is corresponding in section time point, in the model S that step 2 sub-step (a) is set up In find correspondence flow value Bmin, simultaneously according to the minima of this pressure monitoring point in current the d time period, The model A set up in step 1 sub-step (f)uIn find this pressure monitoring point in current the d time period Flow value C corresponding to minimamin, by BmaxWith CmaxCompare, BminWith CminCompare, if more than described Two all exceed system thresholds, then judge that this pressure monitoring point place pipeline section exists leakage, complete leakage point Orientation judges;
Step 3, according in step 2, the orientation of leakage point is judged, leak point positioning unit starting be installed on The sonic transducer of this orientation relevant position completes sonic data collection and calculates, according to located in connection survey calculation side Method calculates leakage point away from the distance closing on sonic transducer.

Claims (2)

1. a DMA subregion pipeline network leak on-line checking alignment system, including DMA subregion tap water Main line, DMA subregion the 1st, 2 ... N-1, N tap water bye-pass, data acquisition unit, pressure are examined Survey unit, sonic detection unit, system server and work station, it is characterised in that: described data acquisition list Unit, including pressure transducer, flow transducer and data acquisition RTU module, wherein pressure transducer, Flow transducer is arranged on the tap water main line of DMA subregion porch, and data acquisition RTU module is pacified Be contained on the borehole wall of adjacent valve well, described pressure transducer, flow transducer respectively with data acquisition RTU In module, two RS-485 interfaces are connected, for being responsible for by the tap water collecting DMA subregion porch Pressure on road, data on flows by be built in the GPRS/3G/4G module of data acquisition RTU module with Data are uploaded to system server by timing or real-time mode;Described DMA subregion the 1st, 2 ... N-1, On N tap water bye-pass, according to field working conditions be separately installed with pressure sensing cell or sonic detection unit or Pressure sensing cell and sonic detection unit, described pressure sensing cell, including being arranged on DMA subregion certainly Pressure transducer on water bye-pass and the data acquisition RTU mould being arranged on the borehole wall of adjacent valve well Block, described sonic detection unit, including the pipeline sound sensing being arranged on DMA subregion tap water bye-pass Device and the data acquisition RTU module being arranged on the borehole wall of adjacent valve well, described pressure transducer, pipe Road sonic transducer is connected with RS-485 interface in data acquisition RTU module respectively, for collecting DMA The pressure on tap water bye-pass in subregion and sonic data are by being built in data acquisition RTU module GPRS/3G/4G module by timing or real-time in the way of data are uploaded to system server, system server And carry out communication by EPA between work station;Described work station, including data collection control unit, Statistical model parameter calculation unit, pipe network diagnosis unit, sound characteristics computational analysis unit, leakage point leak hunting fixed Bit location and DMA partition monitor display unit, wherein, described data collection control unit, also include adopting It is single that collection target configuration subelement, data acquisition scheme configuration subelement, data transmit-receive and communications protocol resolve son Unit, real-time database data management subelement, history library interface subelement and located in connection data acquisition Collaborative Control Subelement, the data gathered based on data acquisition unit, pressure sensing cell and sonic detection unit, through system Meter model parameter calculation unit and pipe network diagnosis unit jointly complete judgement that whether pipeline network leak exist and let out The judgement in leakage orientation, sound characteristics computational analysis unit completes to be judged to leak the pipeline sonic data around orientation Calculate.
A kind of inspection of DMA subregion pipeline network leak on-line checking alignment system Measure method for position, it is characterised in that comprise the following steps:
Step 1, set up DMA water rationing pipe network operation state rule statistical models, described statistics mould Premise set up by type is that pipeline is properly functioning and No leakage, according to actual condition, is divided into n time by 24 hours Section, has 1 flow monitoring point, p pressure monitoring point, specifically includes following sub-step in DMA subregion Rapid:
A in () single time period, in DMA subregion entry data collecting unit, flow transducer gathers m time DMA subregion tap water main line data on flows, gathers 1 data on flows every time, and each data on flows is all wrapped Containing the acquisition time of this data on flows, find the maximum max of data on flows in this time period1And minima min1, peak-to-peak value P1It is calculated by formula (1),
P1=max1-min1 (1)
The meansigma methods of DMA subregion tap water main line data on flows in this time periodCounted by formula (2) Obtain,
X &OverBar; 1 = 1 m &Sigma; i = 1 m X i - - - ( 2 )
In formula, the DMA subregion tap water total number of main line data on flows that m gathers in representing this time period, XiRepresent the i-th data in m DMA subregion tap water main line data on flows in this time period, single The variance yields s of DMA subregion tap water main line data on flows in the individual time period1 2Calculated by formula (3) Obtain,
s 1 2 = 1 m - 1 &Sigma; i = 1 m ( X i - X &OverBar; 1 ) 2 - - - ( 3 )
M-discharge model when () sets up b, duplicon step (a) n time, n difference in calculating respectively a day In time period, DMA subregion tap water main line data on flows peak-to-peak value, meansigma methods and variance yields, first remember Record DMA subregion tap water main line data on flows peak-to-peak value in the 1st time period, meansigma methods and variance yields, And with the corresponding relation of data on flows Yu this data acquisition time, with the time as transverse axis, DMA in this time period The data on flows of subregion entry data collecting unit flow transducer is that the longitudinal axis sets up curve, i.e. sets up model Q1, then with DMA subregion tap water main line data on flows peak-to-peak value, meansigma methods in the 2nd time period And variance yields, and with the corresponding relation of data on flows Yu this data acquisition time, with the time as transverse axis, this time Between in section the data on flows of DMA subregion entry data collecting unit flow transducer be that the longitudinal axis sets up curve, I.e. set up model Q2, the rest may be inferred records DMA subregion tap water main line flow number in the n-th time period According to peak-to-peak value, meansigma methods and variance yields, and with the corresponding relation of data on flows Yu this data acquisition time, with Time is transverse axis, the data on flows of DMA subregion entry data collecting unit flow transducer in this time period Set up curve for the longitudinal axis, i.e. set up model Qn, according to the change dynamic corrections above-mentioned model Q of data1~Qn
C in () single time period, in pressure sensing cell, pressure transducer gathers k time at single pressure monitoring point Pressure data, gathers 1 pressure data every time, and each pressure data all includes the collection of this pressure data Time, so within the single time period, single monitoring point gathers pressure data k altogether, finds this time period The maximum max of interior pressure data2And minima min2, peak-to-peak value P2Obtained by formula (4),
P2=max2-min2 (4)
The meansigma methods of single pressure monitoring point pressure dataIt is calculated by formula (5),
X &OverBar; 2 = 1 k &Sigma; j = 1 k X j - - - ( 5 )
In formula, the k total number of pressure data that in representing this time period, single pressure monitoring point gathers, XjRepresent Jth data in the pressure data that in this time period, single pressure monitoring point gathers, single pressure monitoring point The variance yields s of pressure data2 2It is calculated by formula (6),
s 2 2 = 1 k - 1 &Sigma; j = 1 k ( X j - X &OverBar; 2 ) 2 - - - ( 6 )
D () duplicon step (c) p time, in calculating the single time period, the 1st pressure detecting point is pressed to pth Peak-to-peak value, meansigma methods and the variance yields of power monitoring point pressure data;
M-pressure model when () sets up e, duplicon step (d) n time, the 1st time in calculating one day To the in the n-th time period the 1st~the peak-to-peak value of p pressure monitoring point pressure data, meansigma methods and variance in section 1st pressure monitoring point pressure data peak-to-peak value, meansigma methods and side in value, first the 1st time period of record Difference, and with the corresponding relation of pressure data Yu this pressure data acquisition time, with the time as transverse axis, the 1st In the individual time period, the 1st pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up model P11;Note Record in the 1st time period the 2nd pressure monitoring point pressure data peak-to-peak value, meansigma methods and variance yields, and with The corresponding relation of pressure data and this pressure data acquisition time, with the time as transverse axis, in this time period the 2nd Individual pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up model P12;The rest may be inferred this time period Interior pth pressure monitoring point pressure data peak-to-peak value, meansigma methods and variance yields, and with pressure data and this pressure The corresponding relation of force data acquisition time, with the time as transverse axis, in this time period, pth pressure monitoring presses Force data is that the longitudinal axis sets up curve, i.e. sets up model P1p, in like manner, respectively in the 2nd time period of record the 1~p pressure monitoring point pressure data peak-to-peak value, meansigma methods and variance yields, with pressure in the 2nd time period Data and the corresponding relation of this pressure data acquisition time, with the time as transverse axis, pressure monitoring in this time period Point pressure data are that the longitudinal axis sets up curve, i.e. set up model P21~P2p, the rest may be inferred when recording n-th respectively Between in section the 1st~p pressure monitoring point pressure data peak-to-peak value, meansigma methods and variance yields, with the n-th time Pressure data and the corresponding relation of this pressure data acquisition time in section, with the time as transverse axis, in this time period Pressure monitoring point pressure data is that the longitudinal axis sets up curve, i.e. sets up model Pn1~Pnp.Change according to data Dynamic corrections above-mentioned model P11~P1p、P21~P2pUntil Pn1~Pnp
F () sets up flow rate pressure model, in recording one day DMA subregion tap water main line data on flows and P pressure monitoring point data, first, with different DMA subregion tap water main line datas on flows and the 1st The corresponding relation of individual pressure monitoring point pressure data, flow is transverse axis, and pressure is that the longitudinal axis sets up the 1st pressure Monitoring point flow rate pressure model A1, then, with different DMA subregion tap water main line datas on flows with The corresponding relation of the 2nd pressure monitoring point pressure data, flow is transverse axis, and pressure is that the longitudinal axis sets up the 2nd Pressure monitoring point flow rate pressure model A2, the rest may be inferred with different DMA subregion tap water main line flows Data and the corresponding relation of pth pressure monitoring point pressure data, flow is transverse axis, and pressure is that the longitudinal axis is set up Pth pressure monitoring point flow rate pressure model Ap
Whether step 2, leak diagnostics unit judges exist leakage point and leakage orientation, according to practical situation, Setup algorithm widow time length, widow time length is designated as t, and system of the present invention can set each item number According to threshold value, specifically include following sub-step:
A () sets current slot as the d time period, by formula in step 1 sub-step (a) (1) Calculate the peak-to-peak value of DMA subregion tap water main line data on flows in current the d time period, by step DMA subregion tap water supervisor in formula (2) calculates current the d time period in rapid 1 sub-step (a) The meansigma methods of road data on flows, when calculating current the d by formula in step 1 sub-step (a) (3) Between the variance yields of DMA subregion tap water main line data on flows in section, record the above current d Peak-to-peak value, meansigma methods and variance yields in time period, with the corresponding relation of data on flows Yu this data acquisition time, With the time as transverse axis, in current the d time period, data on flows is that the longitudinal axis sets up curve, i.e. sets up model S, By peak-to-peak value, meansigma methods and variance yields and the time m-flow of foundation in step 1 sub-step (b) in model S Model Q1~QnIn the model Q of corresponding time perioddCompared with middle peak-to-peak value, meansigma methods are corresponding with variance yields, If the difference of peak-to-peak value all exceedes system threshold not less than system thresholds, the difference of meansigma methods and the difference of variance yields Value, then judge that this DMA subregion exists leakage, if the difference of peak-to-peak value, meansigma methods and variance yields all exceedes System thresholds, then with the length of window t calculated flow rate statistical average of default within this time period Variance yields st 2If this time period initial time is T1, the end time is T2, this time period is divided into the 1st Individual window T1~T1+ t, the 2nd window t~T1+ 2t ..., x-th window T1+ xt~T2, at y Calculated flow rate statistical average in individual windowVariance yields st 2, meansigma methodsCalculated by formula (7),
X &OverBar; t = 1 e &Sigma; f = 1 e X f - - - ( 7 )
In formula, e represents the total number of data on flows in y-th window, XfRepresent flow number in y-th window According to the f data in total number, variance yields st 2Calculated by formula (8),
s t 2 = 1 e - 1 &Sigma; f = 1 e ( X f - X &OverBar; t ) 2 - - - ( 8 )
Find the 1st meansigma methods of data on flows, variance yields discharge model m-with time corresponding in this time period QdMeansigma methods, variance yields compare the window more all exceeding system thresholds, record initial time t of this window1; The meansigma methods of last 1 data on flows, variance yields discharge model Q m-with time corresponding in finding this time periodd Meansigma methods, variance yields compare the window more all exceeding system thresholds, record initial time t of this window2, 1st occurs to last, abnormal window occurs that the time difference between abnormal window is designated as Δ t, Δ t and leads to Cross formula (9) to calculate,
Δ t=t2-t1 (9)
If Δ t exceedes the maximum time interval of default, then judge to there is leakage;
If b in judging this DMA subregion in () step 2 sub-step (a) there is leakage, then according to step in pipe network The peak-to-peak of p pressure monitoring point in formula (4) calculates current the d time period in rapid 1 sub-step (c) Value, calculate p pressure prison in current the d time period according to formula in step 1 sub-step (c) (5) The meansigma methods of measuring point and calculating in current the d time period according to formula in step 1 sub-step (c) (6) The variance yields of p pressure monitoring point, if the peak-to-peak value of u pressure monitoring point, meansigma methods and variance yields its In one or two or three exceed system correspondence threshold value, then find this pressure prison in current the d time period The maximum of measuring point and minima, corresponding according to the maximum of this pressure monitoring point in current the d time period Acquisition time, find in the model S that step 2 sub-step (a) is set up correspondence flow value Bmax, Simultaneously according to the maximum of this pressure monitoring point in current the d time period, in step 1 sub-step (f) The model A of middle foundationuIn find the stream that the maximum of this pressure monitoring point is corresponding in current the d time period Value Cmax;According to the time point that the minima of this pressure monitoring point in current the d time period is corresponding, The model S that step 2 sub-step (a) is set up finds the flow value B of correspondencemin, simultaneously according to current d The minima of this pressure monitoring point in the individual time period, the model A set up in step 1 sub-step (f)uIn Find the flow value C that the minima of this pressure monitoring point interior of current the d time period is correspondingmin, by Bmax With CmaxCompare, BminWith CminCompare, if described above two all exceed system thresholds, then judge this pressure There is leakage in place, power monitoring point pipeline section, completes the orientation to leakage point and judge;
Step 3, according in step 2, the orientation of leakage point is judged, leak point positioning unit starting be installed on The sonic transducer of this orientation relevant position completes sonic data collection and calculates, according to located in connection survey calculation Method calculates leakage point away from the distance closing on sonic transducer.
CN201610302392.6A 2016-05-07 2016-05-07 DMA subregions pipeline network leak on-line checking alignment system and its detection localization method CN105927863B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610302392.6A CN105927863B (en) 2016-05-07 2016-05-07 DMA subregions pipeline network leak on-line checking alignment system and its detection localization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610302392.6A CN105927863B (en) 2016-05-07 2016-05-07 DMA subregions pipeline network leak on-line checking alignment system and its detection localization method

Publications (2)

Publication Number Publication Date
CN105927863A true CN105927863A (en) 2016-09-07
CN105927863B CN105927863B (en) 2018-06-29

Family

ID=56835484

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610302392.6A CN105927863B (en) 2016-05-07 2016-05-07 DMA subregions pipeline network leak on-line checking alignment system and its detection localization method

Country Status (1)

Country Link
CN (1) CN105927863B (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106895268A (en) * 2017-03-28 2017-06-27 智润科技有限公司 Set up the method and monitoring system of tree structure leakage loss condition monitoring pipe network
CN107355688A (en) * 2017-07-14 2017-11-17 水联网技术服务中心(北京)有限公司 A kind of LeakView urban water supplies pipe network model Control management system
CN107766468A (en) * 2017-09-30 2018-03-06 上海邦芯物联网科技有限公司 A kind of Intelligent pipe net monitoring system based on big data
CN108916663A (en) * 2018-09-08 2018-11-30 北京逸智联科技有限公司 A kind of municipal network water supply monitoring system and monitoring method
CN108980635A (en) * 2017-06-05 2018-12-11 北京菁华昱创节能设备有限公司 A kind of water system protective device of intelligent and high-efficiency
CN109555979A (en) * 2018-12-10 2019-04-02 清华大学 A kind of water supply network leakage monitoring method
CN109668054A (en) * 2017-10-13 2019-04-23 中国石油天然气股份有限公司 A kind of pipe cleaner positioning system and method
CN110185092A (en) * 2019-04-19 2019-08-30 浙江大学 A kind of urban water supply system leakage monitoring method based on dynamic DMA subregion
CN110836784A (en) * 2018-08-16 2020-02-25 广州极飞科技有限公司 Spraying system, fault detection method based on spraying system and fault type determination method
CN111006137A (en) * 2019-12-18 2020-04-14 北京无线电计量测试研究所 Water supply pipeline leakage monitoring and leakage positioning method and system
CN111271608A (en) * 2020-03-05 2020-06-12 北京中竞国际能源科技有限公司 Leakage management system and method for compressed air system
CN111609324A (en) * 2020-05-29 2020-09-01 北京化工大学 Pipeline leakage detection method and device

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202082629U (en) * 2011-04-06 2011-12-21 黄定军 Natural gas pipeline leakage monitoring system
WO2012059390A1 (en) * 2010-11-05 2012-05-10 Siemens Aktiengesellschaft Leakage detection and leakage location in supply networks
CN103939749A (en) * 2014-04-28 2014-07-23 东北大学 Oil pipeline network leakage intelligent self-adaptation monitoring system and method based on big data
CN205090197U (en) * 2015-07-09 2016-03-16 南京声宏毅霆网络科技有限公司 Leak detection system and monitoring facilities of pipeline

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012059390A1 (en) * 2010-11-05 2012-05-10 Siemens Aktiengesellschaft Leakage detection and leakage location in supply networks
CN202082629U (en) * 2011-04-06 2011-12-21 黄定军 Natural gas pipeline leakage monitoring system
CN103939749A (en) * 2014-04-28 2014-07-23 东北大学 Oil pipeline network leakage intelligent self-adaptation monitoring system and method based on big data
CN205090197U (en) * 2015-07-09 2016-03-16 南京声宏毅霆网络科技有限公司 Leak detection system and monitoring facilities of pipeline

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨进等: "自适应供水管网泄漏检测定位仪器系统", 《仪表技术与传感器》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106895268A (en) * 2017-03-28 2017-06-27 智润科技有限公司 Set up the method and monitoring system of tree structure leakage loss condition monitoring pipe network
CN108980635A (en) * 2017-06-05 2018-12-11 北京菁华昱创节能设备有限公司 A kind of water system protective device of intelligent and high-efficiency
CN107355688A (en) * 2017-07-14 2017-11-17 水联网技术服务中心(北京)有限公司 A kind of LeakView urban water supplies pipe network model Control management system
CN107766468A (en) * 2017-09-30 2018-03-06 上海邦芯物联网科技有限公司 A kind of Intelligent pipe net monitoring system based on big data
CN109668054A (en) * 2017-10-13 2019-04-23 中国石油天然气股份有限公司 A kind of pipe cleaner positioning system and method
CN109668054B (en) * 2017-10-13 2020-11-03 中国石油天然气股份有限公司 Pipe cleaner positioning system and method
CN110836784A (en) * 2018-08-16 2020-02-25 广州极飞科技有限公司 Spraying system, fault detection method based on spraying system and fault type determination method
CN108916663A (en) * 2018-09-08 2018-11-30 北京逸智联科技有限公司 A kind of municipal network water supply monitoring system and monitoring method
CN109555979A (en) * 2018-12-10 2019-04-02 清华大学 A kind of water supply network leakage monitoring method
CN110185092A (en) * 2019-04-19 2019-08-30 浙江大学 A kind of urban water supply system leakage monitoring method based on dynamic DMA subregion
CN111006137A (en) * 2019-12-18 2020-04-14 北京无线电计量测试研究所 Water supply pipeline leakage monitoring and leakage positioning method and system
CN111271608A (en) * 2020-03-05 2020-06-12 北京中竞国际能源科技有限公司 Leakage management system and method for compressed air system
CN111609324A (en) * 2020-05-29 2020-09-01 北京化工大学 Pipeline leakage detection method and device

Also Published As

Publication number Publication date
CN105927863B (en) 2018-06-29

Similar Documents

Publication Publication Date Title
US9410833B1 (en) Methods and apparatus for fluid flow measurement
US10775213B2 (en) Devices and system for channeling and automatic monitoring of fluid flow in fluid distribution systems
TWI502299B (en) System and method for monitoring resources in a water utility network
CN104236626B (en) The integrated on-line monitoring system of drainage pipeline liquid level and flow
RU2386992C2 (en) Operating parametre transducer incorporating diagnostics capabilities
CN101761780B (en) Gas pipeline leakage detecting and positioning device and method thereof
US6789411B2 (en) Detecting leaks in pipes
US10387590B2 (en) Techniques for iterative reduction of uncertainty in water distribution networks
CN104595729B (en) A kind of oil and gas pipeline leakage localization method based on magnitudes of acoustic waves
US20140152465A1 (en) System and method for identifying likely geographical locations of anomalies in a water utility network
CA2440397C (en) Sewer flow monitoring method and system
JP4612696B2 (en) Water leakage diagnosis device and water leakage diagnosis method for water distribution pipe
US8072340B2 (en) Water leakage monitoring system
US8766806B2 (en) Method and apparatus for real time enhancing of the operation of a fluid transport pipeline
Whittle et al. WATERWISE@ SG: A testbed for continuous monitoring of the water distribution system in singapore
CN101506629B (en) Flow measurement diagnostics
TWI549072B (en) Method, medium, apparatus, and system for sensing events affecting liquid flow in a liquid distribution system
CN103221726B (en) Method and apparatus for controlling gas flow via a gas shut-off valve assembly
AU2014324846B2 (en) Non-intrusive sensor system
Loureiro et al. Water distribution systems flow monitoring and anomalous event detection: A practical approach
CN105864642B (en) A kind of pipeline leakage testing device and detection method
CN102235575B (en) Data processing method and system for checking pipeline leakage
JP5094374B2 (en) Flow measurement device and gas supply system
CN101684894B (en) Method and device for monitoring pipeline leakage
CN203147291U (en) System capable of monitoring pipeline leakage by means of infrasonic waves, flow balance and negative pressure waves

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
GR01 Patent grant