CN109388904B - Method and system for calculating leakage amount based on DMA partition flow - Google Patents
Method and system for calculating leakage amount based on DMA partition flow Download PDFInfo
- Publication number
- CN109388904B CN109388904B CN201811267944.XA CN201811267944A CN109388904B CN 109388904 B CN109388904 B CN 109388904B CN 201811267944 A CN201811267944 A CN 201811267944A CN 109388904 B CN109388904 B CN 109388904B
- Authority
- CN
- China
- Prior art keywords
- flow
- data
- calculating
- flow data
- module
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Databases & Information Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Software Systems (AREA)
- Probability & Statistics with Applications (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Examining Or Testing Airtightness (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measuring Volume Flow (AREA)
Abstract
The invention discloses a method and a system for calculating leakage based on DMA partition flow, which comprises the following steps: installing a flow meter at an inflow point of the DMA partition, and setting flow meter data uploading interval time; uploading the sampling data of the flowmeter to an SCADA system in real time, calling the sampling data in the SCADA system by the water supply management system, and storing the sampling data in a database of the water supply management system; and acquiring first flow data and second flow data according to the sampling data, determining a minimum night flow time period, and judging whether the second flow data conforms to normal distribution or not, thereby determining the leakage loss of the DMA partition. By calculating the leakage quantity after the second flow data is subjected to the normality judgment, the accuracy of a leakage quantity calculation result can be improved, and the calculation efficiency can be improved.
Description
Technical Field
The invention relates to the technical field of leakage calculation of a water supply network, in particular to a method and a system for calculating leakage based on DMA partition flow.
Background
With the rapid development of urbanization, the importance of the water supply system as a part of municipal utilities is further highlighted, and the problem of water supply leakage is troubling water supply companies, water administration departments and the like. According to' 2016 Town statistics yearbook of urban and rural statistics, it is shown that the water leakage of 2016 urban public water supply in China is 78 billionth of cubic meters, and a large amount of water resources are lost. Because the leakage has various forms, the leakage can be caused by various reasons, such as the breakage of a pipeline main body, the breakage of a pipeline joint, the water leakage of pipeline accessories (a brake, a valve, a fire hydrant and the like), and the calculation of the actual leakage is extremely complicated.
By analyzing the nighttime flow data of the measured metering cell (DMA) of the water supply pipe network, an approximate real leakage numerical value can be obtained in time. The leakage amount calculation method provided by the prior art has the following defects: firstly, the determination of the night minimum flow time period is not scientific enough, the time period generally adopted is 2; secondly, abnormal data points are not removed, abnormal data points may be generated due to meter problems or network problems in the data transmission process in the data acquisition process, the abnormal points need to be eliminated in the calculation process, and the calculation result can be more accurate and is closer to the actual leakage value.
Disclosure of Invention
In view of this, the present invention provides a method for calculating a leakage amount based on DMA partition flow, which includes: installing a flow meter at an inflow point of the DMA partition, and setting flow meter data uploading interval time; uploading the sampling data of the flowmeter to an SCADA system in real time, calling the sampling data in the SCADA system by the water supply management system, and storing the sampling data in a database of the water supply management system; acquiring all sampling data within one month, and calculating first flow data of the time period by taking 30 minutes as a reference; calculating a first average value u1 and a first standard deviation sigma 1 of all first flow data in the same time period within one month according to the first flow data, and removing the first flow data except for [ u1-3 sigma 1, u1+3 sigma 1] to obtain second flow data; according to the second flow data, calculating a minimum flow time period of 2 continuous hours in one month, and determining a minimum night flow time period of the DMA partition according to the minimum flow time period; acquiring all second flow data in the time period within one month according to the minimum night flow time period, setting an alpha value, calling a kstest method in MATLAB to perform normality judgment, and judging whether the second flow data accords with normal distribution or not according to a judgment result; if the second flow data accords with the normal distribution, calculating a second average value u and a second standard deviation sigma of the second flow data; and calculating u-2 sigma as the leakage loss amount of the DMA partition according to the second average value u and the second standard deviation sigma.
Further, if the second flow data does not conform to the normal distribution, the alpha value is adjusted, and the kstest method in MATLAB is called again to perform normality judgment.
Further, if the second stream data does not conform to the normal distribution after the alpha value is adjusted, the sampled data is re-acquired.
Further, the leakage amount of all the partitions is ranked and displayed by adopting icons.
Further, the alpha value is 0.05 with an interval of less than 10 minutes.
The invention also provides a system for calculating the leakage based on the DMA partition flow, which comprises the following steps: the system comprises a DMA system, an SCADA system and a water supply management system, wherein the DMA system and the water supply management system are respectively and electrically connected with the SCADA system; the DMA system divides the area to be tested into a plurality of DMA subareas; the flow-in point of the DMA partition is connected with a flow meter, and the flow meter uploads sampling data to an SCADA system in real time according to preset interval time; the water supply management system comprises a database, a first calculation module, a data rejection module, a second calculation module, a judgment module, a third calculation module and an output module, wherein the first calculation module, the data rejection module and the second calculation module are all electrically connected with the database, the data rejection module is electrically connected with the first calculation module and the second calculation module, the judgment module is electrically connected with the second calculation module and the third calculation module, and the third calculation module is electrically connected with the output module; the water supply management system calls the sampling data in the SCADA system and stores the sampling data into a database; the first calculation module is used for acquiring all sampling data within one month from the database and calculating first flow data of the time period by taking 30 minutes as a reference; the data removing module is used for calculating a first average value u1 and a first standard deviation sigma 1 of all first flow data in the same time period within one month according to the first flow data, removing the first flow data except for [ u1-3 sigma 1, u1+3 sigma 1], obtaining second flow data and storing the second flow data in a database; the second calculation module is used for calculating a minimum flow time period of 2 continuous hours in one month according to the second flow data, determining a minimum night flow time period of the DMA partition according to the minimum flow time period, and acquiring all the second flow data in the time period in one month according to the minimum night flow time period; the judging module is used for setting an alpha value, calling a kstest method in MATLAB to judge the normality, and judging whether the second flow data accords with normal distribution or not according to a judging result; if the second flow data accords with the normal distribution, the third calculation module calculates a second average value u and a second standard deviation sigma of the second flow data; and the output module is used for calculating u-2 sigma as the leakage loss of the DMA partition according to the second average value u and the second standard deviation sigma and outputting the leakage loss.
Further, if the judging module judges that the second flow data does not conform to the normal distribution, the judging module adjusts an alpha value and recalls a kstest method in the MATLAB to carry out normality judgment.
Furthermore, if the second flow data does not conform to the normal distribution after the determining module adjusts the alpha value, the first calculating module retrieves the sampling data.
Further, the output module ranks the leakage amount of all the partitions and displays the leakage amount by adopting a chart.
Further, the alpha value set by the judging module is 0.05, and the preset interval time of the flowmeter is less than 10 minutes.
Compared with the prior art, the method and the system for calculating the leakage quantity based on the DMA partition flow at least realize the following beneficial effects:
1. the leakage quantity calculation method provided by the invention can judge whether the acquired flow data conforms to normal distribution by using MATLAB, the calculation result conforming to normal distribution is closer to the real leakage data, and the accuracy of the calculation result can be improved.
2. The leakage quantity calculation method provided by the invention can select the night minimum flow time period according to the actual situation, the time period determination mode is more scientific, and the leakage data difference caused by water consumption of a user at night can be avoided to the greatest extent;
3. the leakage quantity calculation method provided by the invention can eliminate abnormal data in the flow data, can avoid result deviation caused by the abnormal data, and enables the calculation result of the leakage quantity to be more accurate;
4. the leakage quantity calculation method provided by the invention has the advantages of simpler check mode, simple operation and convenient realization by adopting a computer, can save labor, improve efficiency and reduce cost, can reduce errors possibly caused by manual operation, and further enhances the accuracy of a calculation result.
Of course, it is not necessary for any product in which the present invention is practiced to achieve all of the above-described technical effects simultaneously.
Other features of the present invention and advantages thereof will become apparent from the following detailed description of exemplary embodiments thereof, which proceeds with reference to the accompanying drawings.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description, serve to explain the principles of the invention.
FIG. 1 is a flow chart of a method for calculating a leakage amount based on DMA partition flow provided by the present invention;
FIG. 2 is a flow chart of another method provided by the present invention for calculating the amount of leakage based on DMA partition traffic;
fig. 3 is a system structural diagram for calculating a leakage amount based on DMA partition traffic according to an embodiment of the present invention.
Detailed Description
Various exemplary embodiments of the present invention will now be described in detail with reference to the accompanying drawings. It should be noted that: the relative arrangement of the components and steps, the numerical expressions and numerical values set forth in these embodiments do not limit the scope of the present invention unless it is specifically stated otherwise.
The following description of at least one exemplary embodiment is merely illustrative in nature and is in no way intended to limit the invention, its application, or uses.
Techniques, methods, and apparatus known to those of ordinary skill in the relevant art may not be discussed in detail but are intended to be part of the specification where appropriate.
In all examples shown and discussed herein, any particular value should be construed as merely illustrative, and not limiting. Thus, other examples of the exemplary embodiments may have different values.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, further discussion thereof is not required in subsequent figures.
Example one
In order to calculate the leakage more accurately, please refer to fig. 1, fig. 1 is a flowchart of a method for calculating the leakage based on the DMA partition traffic according to the present invention. As shown in FIG. 1, the present invention provides a method for calculating a leakage amount based on DMA partition traffic, comprising:
s1: installing a flow meter at an inflow point of the DMA partition, and setting flow meter data uploading interval time;
s2: uploading the sampling data of the flowmeter to an SCADA system in real time, calling the sampling data in the SCADA system by the water supply management system, and storing the sampling data in a database of the water supply management system;
s3: acquiring all sampling data within one month, and calculating first flow data of the time period by taking 30 minutes as a reference;
s4: calculating a first average value u1 and a first standard deviation sigma 1 of all first flow data in the same time period within one month according to the first flow data, and removing the first flow data except for [ u1-3 sigma 1, u1+3 sigma 1] to obtain second flow data;
s5: according to the second flow data, calculating a minimum flow time period of 2 continuous hours in one month, and determining a minimum night flow time period of the DMA partition according to the minimum flow time period;
s6: acquiring all second flow data in the time period within one month according to the minimum night flow time period, setting an alpha value, calling a kstest method in MATLAB to perform normality judgment, and judging whether the second flow data accords with normal distribution or not according to a judgment result;
s7: if the second flow data accords with the normal distribution, calculating a second average value u and a second standard deviation sigma of the second flow data;
s8: and calculating u-2 sigma as the leakage loss amount of the DMA partition according to the second average value u and the second standard deviation sigma.
Specifically, in the method for calculating the leakage amount based on the DMA partition flow provided by this embodiment, the sampling data of the flow meter is uploaded to the SCADA system in real time, the water supply management system calls the sampling data from the SCADA system, and calculates a first average value u1 and a first standard deviation σ 1 of all first flow data in one half hour in one month based on half-hour, removes abnormal data according to the first average value u1 and the first standard deviation σ 1 to obtain second flow data, determines a minimum night flow time period according to the second flow data, and performs normality determination on the second flow data in 2 hours by using MATLAB, and calculates a second average value u and a second standard deviation σ if the second flow data are in accordance with a normal distribution, and calculates the leakage amount of the DMA partition according to the second average value u and the second standard deviation σ.
The method for calculating the leakage amount based on the DMA partition flow provided by the embodiment at least has the following beneficial effects:
1. the leakage amount calculation method provided by this embodiment can determine whether the acquired flow data conforms to normal distribution by using MATLAB, and the calculation result conforming to normal distribution is closer to the real leakage data, so that the accuracy of the calculation result can be improved.
2. The leakage amount calculation method provided by the embodiment can select the night minimum flow time period according to the actual situation, the time period determination mode is more scientific, and the leakage data difference caused by the water consumption of the user at night can be avoided to the greatest extent;
3. the leakage quantity calculation method provided by the embodiment can eliminate abnormal data in the flow data, and can avoid result deviation caused by the abnormal data, so that the calculation result of the leakage quantity is more accurate;
4. the leakage quantity calculation method provided by the embodiment has the advantages of simpler check mode, simplicity in operation and convenience in realization by adopting a computer, can save labor, improve efficiency and reduce cost, can reduce errors possibly caused by manual operation, and further enhances the accuracy of a calculation result.
Example two
Referring to FIG. 2, FIG. 2 is a flow chart of another method for calculating the amount of leakage based on DMA partition flow according to the present invention. As shown in FIG. 2, the present invention provides a method for calculating a leakage amount based on DMA partition traffic, comprising:
s1: installing a flow meter at an inflow point of the DMA partition, and setting the interval time for uploading flow meter data, wherein the interval time is less than 10 minutes;
s2: uploading the sampling data of the flowmeter to an SCADA system in real time, calling the sampling data in the SCADA system by the water supply management system, and storing the sampling data in a database of the water supply management system;
s3: acquiring all sampling data within one month, and calculating first flow data of the time period by taking 30 minutes as a reference;
s4: calculating a first average value u1 and a first standard deviation sigma 1 of all first flow data in the same time period within one month according to the first flow data, and removing the first flow data except for [ u1-3 sigma 1, u1+3 sigma 1] to obtain second flow data;
s5: according to the second flow data, calculating a minimum flow time period of 2 continuous hours in one month, and determining a minimum night flow time period of the DMA partition according to the minimum flow time period;
s6: acquiring all second flow data in the time period within one month according to the minimum night flow time period, setting an alpha value to be 0.05, and calling a kstest method in MATLAB to perform normality judgment;
s7: judging whether the second flow data accords with normal distribution or not according to the judgment result;
s81: if the second flow data does not conform to the normal distribution, adjusting an alpha value, and recalling a kstest method in MATLAB to perform normality judgment;
s82: if the second flow data does not conform to the normal distribution after the alpha value is adjusted, returning to the step 3 to obtain the data again;
s91: if the second flow data accords with the normal distribution, calculating a second average value u and a second standard deviation sigma of the second flow data;
s92: calculating u-2 sigma as the leakage loss of the DMA partition according to the second average value u and the second standard deviation sigma;
s10: and ranking the leakage amount of all the partitions, and displaying by adopting a chart.
Specifically, in the method for calculating the leakage loss based on the DMA partition flow provided in this embodiment, the sampling data of the flow meter is uploaded to the SCADA system in real time, the water supply management system calls the sampling data from the SCADA system, and calculates a first average value u1 and a first standard deviation σ 1 of all first flow data within one half hour in one month on the basis of half an hour, removes abnormal data according to the first average value u1 and the first standard deviation σ 1 to obtain second flow data, determines a minimum night flow time period according to the second flow data, and performs normality determination on the second flow data for 2 hours by using MATLAB, if the second flow data do not conform to the normal distribution, adjusts an alpha value (a parameter value of normality check in MATLAB) to re-determine whether the second flow data conform to the normal distribution, and if the second flow data do not conform to the normal distribution, discards the sampling data and re-acquires the sampling data; and if the second flow data are in accordance with the normal distribution, calculating a second average value u and a second standard deviation sigma, calculating the leakage loss of the DMA partition according to the second average value u and the second standard deviation sigma, ranking the leakage loss of all the partitions, and displaying the ranked leakage loss in a chart form.
The method for calculating the leakage amount based on the DMA partition flow provided by the embodiment at least has the following beneficial effects:
1. the leakage amount calculation method provided by this embodiment can determine whether the acquired flow data conforms to normal distribution by using MATLAB, and the calculation result conforming to normal distribution is closer to the true leakage data, so that the accuracy of the calculation result can be improved.
2. The leakage amount calculation method provided by the embodiment can select the night minimum flow time period according to the actual situation, the time period determination mode is more scientific, and the leakage data difference caused by water consumption of the user at night can be avoided to the greatest extent;
3. the leakage quantity calculation method provided by the embodiment can eliminate abnormal data in the flow data, and can avoid result deviation caused by the abnormal data, so that the calculation result of the leakage quantity is more accurate;
4. the method for calculating the leakage loss is simpler in inspection mode, simple to operate and convenient to realize by adopting a computer, can save labor, improve efficiency and reduce cost, can reduce errors possibly caused by manual operation, and further enhances the accuracy of a calculation result;
5. according to the leakage amount calculation method provided by the embodiment, the leakage amount is displayed by using the chart, and the display mode is more intuitive.
EXAMPLE III
Referring to fig. 3, fig. 3 is a system structure diagram for calculating a leakage amount based on DMA partition flow according to an embodiment of the present invention. As shown in fig. 3, the present embodiment provides a system for calculating a leakage amount based on DMA partition traffic, including: the system comprises a DMA system 1, an SCADA system 2 and a water supply management system 3, wherein the DMA system 1 and the water supply management system 3 are respectively electrically connected with the SCADA system 2.
The DMA system 1 divides the area to be measured into a plurality of DMA subareas; and a flow-in point of the DMA partition is connected with a flowmeter, the flowmeter uploads sampling data to an SCADA system in real time according to preset interval time, and the preset interval time is less than 10 minutes.
The water supply management system 3 comprises a database 30, a first calculation module 31, a second calculation module 32, a data removing module 33, a judgment module 34, a third calculation module 35 and an output module 36, wherein the first calculation module 31, the data removing module 33 and the second calculation module 32 are all electrically connected with the database 30, the data removing module 33 is electrically connected with the first calculation module 31 and the second calculation module 32, the judgment module 34 is electrically connected with the second calculation module 32 and the third calculation module 35, and the third calculation module 35 is electrically connected with the output module 36.
The water supply management system 3 calls the sampled data in the SCADA system 2 and stores the sampled data in the database 30.
The first calculation module 32 is configured to obtain all the sampled data in one month from the database 30, and calculate the first flow data for the time period based on 30 minutes.
The data removing module 33 is configured to calculate a first average value u1 and a first standard deviation σ 1 of all first flow data in the same time period within one month according to the first flow data, remove the first flow data other than [ u1-3 σ 1, u1+3 σ 1], obtain a second flow data, and store the second flow data in the database 30.
The second calculating module 32 is configured to calculate a minimum flow time period of 2 consecutive hours in one month according to the second flow data, determine a minimum night flow time period of the DMA partition according to the minimum flow time period, and obtain all the second flow data in the time period in one month according to the minimum night flow time period.
The judging module 34 is configured to set an alpha value, where the alpha value is 0.05, call a kstest method in the MATLAB to perform normality judgment, judge whether the second flow data conforms to normal distribution according to a judgment result, if the second flow data does not conform to the normal distribution, the judging module 34 adjusts the alpha value, and re-call the kstest method in the MATLAB to perform normality judgment, and if the second flow data does not conform to the normal distribution after the judging module 34 adjusts the alpha value, the first calculating module 31 re-acquires the sampling data; if the second flow data conforms to the normal distribution, the third calculating module 35 calculates a second average u and a second standard deviation σ of the second flow data;
and the output module 36 is configured to calculate u-2 σ as the leakage amount of the DMA partition according to the second average value u and the second standard deviation σ, and output the leakage amount. The output module 36 ranks the leakage amount of all the partitions and displays the leakage amount by using a chart.
Specifically, the system for calculating the leakage amount based on the DMA partition flow provided by this embodiment includes: the system comprises a DMA system 1, a SCADA system 2 and a water supply management system 3, wherein the water supply management system 3 comprises a database 30, a first calculation module 31, a second calculation module 32, a data elimination module 33, a judgment module 34, a third calculation module 35 and an output module 36. The flow meter uploads the sampling data to the SCADA system in real time according to preset interval time, and the water supply management system calls the sampling data from the SCADA system and stores the sampling data in the database 30. The first calculation module 31 obtains all sampling data in one month from the database 30, calculates the first flow data of the time period on the basis of half an hour, the data removing module 33 calculates a first average value u1 and a first standard deviation sigma 1 of all the first flow data in the half hour in one month, removes abnormal data according to the first average value u1 and the first standard deviation sigma 1 to obtain second flow data, and stores the second flow data in the database 30, then the second calculation module 32 determines a minimum night flow time period according to the second flow data, the determination module 34 uses MATLAB to perform normality determination on the second flow data in the 2 hours, if the second flow data do not accord with the normal distribution, an alpha value (a parameter value of the normality check in the MATLAB) is adjusted to determine whether the second flow data accord with the normal distribution again, if the second flow data do not accord with the normal distribution, the batch of sampling data is abandoned, and the first calculation module 31 obtains the sampling data again; if the second flow data conforms to the normal distribution, the third calculation module 35 calculates a second average value u and a second standard deviation σ, and the output module 36 calculates and outputs the leakage loss amount of the DMA partition according to the second average value u and the second standard deviation σ. The output module 36 ranks the leakage amount of all the partitions and displays the leakage amount in a chart form.
The system for calculating the leakage amount based on the DMA partition flow provided by the embodiment at least has the following beneficial effects:
1. the leakage calculation system provided by this embodiment can determine whether the acquired flow data conforms to normal distribution by using MATLAB, and the calculation result conforming to normal distribution is closer to the real leakage data, so that the accuracy of the calculation result can be improved.
2. The leakage amount calculation system provided by the embodiment can select the night minimum flow time period according to the actual situation, the time period determination mode is more scientific, and the leakage data difference caused by water consumption of the user at night can be avoided to the greatest extent;
3. the leakage calculation system provided by the embodiment can eliminate abnormal data in the flow data, and can avoid result deviation caused by the abnormal data, so that the calculation result of the leakage is more accurate;
4. the detection mode of the leakage calculation system provided by the embodiment is simpler, the operation is simple, the computer implementation is convenient, the labor can be saved, the efficiency is improved, the cost is reduced, the error possibly caused by manual operation can be reduced, and the accuracy of the calculation result is further enhanced;
5. the leakage calculation system provided by the invention utilizes the chart to display the leakage, and the display mode is more visual.
By the embodiment, the method and the system for calculating the leakage quantity based on the DMA partition flow, provided by the invention, at least have the following beneficial effects:
1. the leakage quantity calculation method provided by the invention can judge whether the acquired flow data conforms to normal distribution by using MATLAB, the calculation result conforming to normal distribution is closer to the real leakage data, and the accuracy of the calculation result can be improved.
2. The leakage quantity calculation method provided by the invention can select the night minimum flow time period according to the actual situation, the time period determination mode is more scientific, and the leakage data difference caused by water consumption of a user at night can be avoided to the greatest extent;
3. the leakage quantity calculation method provided by the invention can eliminate abnormal data in the flow data, can avoid result deviation caused by the abnormal data, and enables the calculation result of the leakage quantity to be more accurate;
4. the leakage calculation method provided by the invention has the advantages of simpler check mode, simple operation and convenience for realization by a computer, can save labor, improve efficiency and reduce cost, can reduce errors possibly caused by manual operation, and further enhances the accuracy of a calculation result.
Although some specific embodiments of the present invention have been described in detail by way of examples, it should be understood by those skilled in the art that the above examples are for illustrative purposes only and are not intended to limit the scope of the present invention. It will be appreciated by those skilled in the art that modifications may be made to the above embodiments without departing from the scope and spirit of the invention. The scope of the invention is defined by the appended claims.
Claims (10)
1. A method for calculating leakage based on DMA partition flow, comprising:
installing a flow meter at an inflow point of the DMA partition, and setting flow meter data uploading interval time;
uploading the sampling data of the flowmeter to an SCADA system in real time, calling the sampling data in the SCADA system by a water supply management system, and storing the sampling data into a database of the water supply management system;
acquiring all sampling data within one month, and calculating first flow data of one month by taking 30 minutes as a reference;
calculating a first average value u1 and a first standard deviation sigma 1 of all first flow data in the same time period within one month according to the first flow data, and removing the first flow data except for [ u1-3 sigma 1, u1+3 sigma 1] to obtain second flow data;
according to the second flow data, calculating a minimum flow time period of 2 continuous hours in one month, and determining a minimum night flow time period of the DMA partition according to the minimum flow time period;
acquiring all second flow data in the minimum night flow time period within one month according to the minimum night flow time period, setting an alpha value, calling a kstest method in MATLAB to perform normality judgment, and judging whether the second flow data accords with normal distribution or not according to a judgment result;
calculating a second mean u and a second standard deviation σ of the second flow data if the second flow data conforms to the normal distribution;
and calculating u-2 sigma as the leakage amount of the DMA partition according to the second average value u and the second standard deviation sigma.
2. The method of calculating the amount of leakage based on DMA partition traffic of claim 1,
and if the second flow data does not conform to the normal distribution, adjusting the alpha value, and recalling a kstest method in MATLAB to perform normality judgment.
3. The method of calculating the amount of leakage based on DMA partition traffic of claim 2,
and if the second flow data still does not conform to the normal distribution after the alpha value is adjusted, the sampling data is obtained again.
4. The method of calculating the amount of leakage based on DMA partition traffic of claim 1,
and ranking the leakage amount of all the partitions, and displaying by adopting a chart.
5. The method of calculating the amount of leakage based on DMA partition traffic of claim 1,
the alpha value is 0.05 and the interval time is less than 10 minutes.
6. A system for calculating a leakage amount based on DMA partition traffic, comprising: the system comprises a DMA system, an SCADA system and a water supply management system, wherein the DMA system and the water supply management system are respectively and electrically connected with the SCADA system;
the DMA system divides a region to be detected into a plurality of DMA subareas; the flow-in point of the DMA subarea is connected with a flowmeter, and the flowmeter uploads sampling data to an SCADA system in real time according to preset interval time;
the water supply management system comprises a database, a first calculation module, a data elimination module, a second calculation module, a judgment module, a third calculation module and an output module, wherein the first calculation module, the data elimination module and the second calculation module are all electrically connected with the database, the data elimination module is electrically connected with the first calculation module and the second calculation module, the judgment module is electrically connected with the second calculation module and the third calculation module, and the third calculation module is electrically connected with the output module;
the water supply management system calls the sampling data in the SCADA system and stores the sampling data into the database;
the first calculation module is used for acquiring all sampling data in one month from the database and calculating first flow data of one month time by taking 30 minutes as a reference;
the data removing module is used for calculating a first average value u1 and a first standard deviation sigma 1 of all first flow data in the same time period within one month according to the first flow data, removing the first flow data except for [ u1-3 sigma 1, u1+3 sigma 1], obtaining second flow data and storing the second flow data into the database;
the second calculation module is used for calculating a minimum flow time period of 2 continuous hours in one month according to the second flow data, determining a minimum night flow time period of the DMA partition according to the minimum flow time period, and acquiring all the second flow data in the minimum night flow time period in one month according to the minimum night flow time period;
the judging module is used for setting an alpha value, calling a kstest method in MATLAB to judge the normality, and judging whether the second flow data accords with normal distribution or not according to a judging result;
the third calculation module calculates a second average value u and a second standard deviation σ of the second flow data if the second flow data conforms to the normal distribution;
and the output module is used for calculating u-2 sigma as the leakage loss of the DMA partition according to the second average value u and the second standard deviation sigma and outputting the leakage loss.
7. The system for calculating the amount of leakage based on DMA partition traffic of claim 6,
and if the judging module judges that the second flow data does not conform to the normal distribution, the judging module adjusts the alpha value and recalls a kstest method in MATLAB to carry out normality judgment.
8. The system for calculating the amount of leakage based on DMA partition traffic of claim 7,
and if the second flow data still does not conform to the normal distribution after the alpha value is adjusted by the judging module, the first calculating module acquires the sampling data again.
9. The system for calculating the amount of leakage based on DMA partition traffic of claim 6,
and the output module ranks the leakage quantities of all the partitions and displays the leakage quantities by adopting a chart.
10. The system for calculating the amount of leakage based on DMA partition traffic of claim 6,
the alpha value set by the determination module is 0.05, and the interval time preset by the flowmeter is less than 10 minutes.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811267944.XA CN109388904B (en) | 2018-10-29 | 2018-10-29 | Method and system for calculating leakage amount based on DMA partition flow |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811267944.XA CN109388904B (en) | 2018-10-29 | 2018-10-29 | Method and system for calculating leakage amount based on DMA partition flow |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109388904A CN109388904A (en) | 2019-02-26 |
| CN109388904B true CN109388904B (en) | 2023-04-07 |
Family
ID=65427145
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811267944.XA Active CN109388904B (en) | 2018-10-29 | 2018-10-29 | Method and system for calculating leakage amount based on DMA partition flow |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109388904B (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110185092B (en) * | 2019-04-19 | 2020-07-07 | 浙江大学 | Dynamic DMA partition-based leakage monitoring method for town water supply system |
| CN110569248B (en) * | 2019-07-31 | 2021-08-03 | 杭州电子科技大学 | Improved SPC (selective pressure control) residential water supply leakage monitoring and early warning method |
| CN113405026B (en) * | 2020-03-16 | 2022-11-11 | 中国港湾工程有限责任公司 | Newly-added leakage judgment and leakage rate estimation method suitable for water supply pipe network in underdeveloped area |
| CN111982211B (en) * | 2020-07-28 | 2022-09-16 | 杭州电子科技大学 | Method for detecting water supply leakage amount of community |
| CN112559969B (en) * | 2020-12-10 | 2022-04-29 | 北部湾大学 | Small leakage detection method based on accumulation sum algorithm |
| CN112487364B (en) * | 2020-12-10 | 2022-03-11 | 北部湾大学 | Small leakage detection method based on exponential weighted moving average algorithm |
| CN113864664B (en) * | 2021-09-29 | 2023-08-15 | 广东粤海水务投资有限公司 | Pipe network leakage early warning method and system based on flow distribution probability calculation |
| CN114135794B (en) * | 2021-11-22 | 2023-11-24 | 杭州数梦工场科技有限公司 | Method and device for detecting leakage of water network |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103953098A (en) * | 2014-04-22 | 2014-07-30 | 中国科学院生态环境研究中心 | Pipe network loss monitoring method |
| CN106677263A (en) * | 2017-02-16 | 2017-05-17 | 中国科学院生态环境研究中心 | Pipe network real leakage loss determination method and system |
| CN107612749A (en) * | 2017-10-17 | 2018-01-19 | 广东青藤环境科技有限公司 | A kind of DMA zone meterings and leakage loss analyzing and positioning monitoring big data platform and method |
| CN107767012A (en) * | 2017-08-25 | 2018-03-06 | 浙江钛合仪器有限公司 | A kind of water supply network water leakage management system and its method for building up and system application |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2554289B (en) * | 2015-05-13 | 2020-09-23 | Nec Corp | Water-leak state estimation system, method, and recording medium |
-
2018
- 2018-10-29 CN CN201811267944.XA patent/CN109388904B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103953098A (en) * | 2014-04-22 | 2014-07-30 | 中国科学院生态环境研究中心 | Pipe network loss monitoring method |
| CN106677263A (en) * | 2017-02-16 | 2017-05-17 | 中国科学院生态环境研究中心 | Pipe network real leakage loss determination method and system |
| CN107767012A (en) * | 2017-08-25 | 2018-03-06 | 浙江钛合仪器有限公司 | A kind of water supply network water leakage management system and its method for building up and system application |
| CN107612749A (en) * | 2017-10-17 | 2018-01-19 | 广东青藤环境科技有限公司 | A kind of DMA zone meterings and leakage loss analyzing and positioning monitoring big data platform and method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109388904A (en) | 2019-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109388904B (en) | Method and system for calculating leakage amount based on DMA partition flow | |
| CN108984873B (en) | Water supply network real-time leakage detection method, device, system and storage medium | |
| CN112035784B (en) | Method and system for determining loss of low-voltage transformer area based on electricity consumption acquisition data | |
| CN110647447B (en) | Abnormal instance detection method, device, equipment and medium for distributed system | |
| CN112597263B (en) | Pipe network detection data abnormity judgment method and system | |
| CN107741577A (en) | A kind of Source of Gateway Meter degree of accuracy on-line monitoring and analysis method and system | |
| CN111998918A (en) | Error correction method, error correction device and flow sensing system | |
| CN104534285A (en) | Energy consumption anomaly monitoring method and device | |
| CN115841398B (en) | Water supply network leakage analysis method and system based on Internet of things | |
| CN201373758Y (en) | Natural gas flow monitoring and checking device through computer measurement | |
| CN112985713A (en) | Pipe network leakage monitoring method and system based on edge calculation | |
| US9205293B2 (en) | Fluid flow testing system | |
| CN108801320B (en) | Diagnosis method of natural gas measurement system | |
| CN202886434U (en) | Civil metering ammeter electricity theft and leakage monitoring system | |
| CN108062616A (en) | A kind of multiplexing electric abnormality detection model with constraints | |
| CN116186950A (en) | Processing method, system, device and storage medium for gas pipe network simulation | |
| CN203746069U (en) | Low-voltage line loss analyzing system in public transformer power supply area | |
| CN113033842A (en) | Control management method and control management system for leakage of water supply network system | |
| CN113432666B (en) | Agricultural underground water exploitation amount measuring method based on dynamic calculation of electric-water conversion coefficient | |
| CN114331128A (en) | Intelligent water affair management method, system, equipment and medium based on Internet of things | |
| CN201867649U (en) | Testing system for testing collection precision of instrument control system | |
| CN112377821A (en) | Pipeline leakage troubleshooting method and system based on intelligent water meter platform big data | |
| CN213337805U (en) | Digital electric energy meter with error test function | |
| CN113986897A (en) | Multi-source data fusion method and device based on hydrological robot | |
| CN212779731U (en) | Pressure transmitter verification and calibration device and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |