WO2016174958A1 - 漏水発生位置推定装置、システムおよび方法 - Google Patents

漏水発生位置推定装置、システムおよび方法 Download PDF

Info

Publication number
WO2016174958A1
WO2016174958A1 PCT/JP2016/058918 JP2016058918W WO2016174958A1 WO 2016174958 A1 WO2016174958 A1 WO 2016174958A1 JP 2016058918 W JP2016058918 W JP 2016058918W WO 2016174958 A1 WO2016174958 A1 WO 2016174958A1
Authority
WO
WIPO (PCT)
Prior art keywords
occurrence
value
water
position estimation
sensor
Prior art date
Application number
PCT/JP2016/058918
Other languages
English (en)
French (fr)
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 GB1716409.6A priority Critical patent/GB2555535B/en
Publication of WO2016174958A1 publication Critical patent/WO2016174958A1/ja

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E03—WATER SUPPLY; SEWERAGE
    • E03B—INSTALLATIONS OR METHODS FOR OBTAINING, COLLECTING, OR DISTRIBUTING WATER
    • E03B7/00—Water main or service pipe systems
    • E—FIXED CONSTRUCTIONS
    • E03—WATER SUPPLY; SEWERAGE
    • E03B—INSTALLATIONS OR METHODS FOR OBTAINING, COLLECTING, OR DISTRIBUTING WATER
    • E03B7/00—Water main or service pipe systems
    • E03B7/07—Arrangement of devices, e.g. filters, flow controls, measuring devices, siphons or valves, in the pipe systems
    • E03B7/071—Arrangement of safety devices in domestic pipe systems, e.g. devices for automatic shut-off
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/26—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors
    • G01M3/28—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/26—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors
    • G01M3/28—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds
    • G01M3/2807—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds for pipes
    • G01M3/2815—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds for pipes using pressure measurements
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/15—Leakage reduction or detection in water storage or distribution

Definitions

  • the present invention relates to an apparatus, a system, and a method for estimating a leakage occurrence position in a distribution pipe network.
  • Patent document 1 states that “a calculation method for determining the geographical location of anomalies, which is a test designed to statistically determine one probable geographical location of anomalies within a district or zone. Multiple executions, where each test execution produces a single result, and the calculation method combines the results of multiple tests to generate a score for the likely location where the anomaly was determined. '' is doing.
  • the present invention provides a leak occurrence position estimation system and method that can estimate the leak occurrence position even when there is an influence of a change in the water distribution state due to other factors including a change in demand.
  • the disclosed water leakage occurrence position estimation device includes a measurement value collection unit that collects first measurement values of a plurality of sensors arranged in a distribution pipe network, and a prediction in a predetermined water distribution state based on the first measurement values of the sensors.
  • the indicator value calculation unit that calculates the index value from the predicted measurement value in the same distribution state as the distribution state that measured the measured value of the water, and water leakage based on the arrangement of multiple sensors and the comparison of the index values between the multiple sensors It has a position estimation part which estimates the area which satisfies the positional relationship with the generation position as the position of the occurrence of water leakage.
  • the leak occurrence position can be estimated even if there is a change in the water distribution state due to other factors.
  • FIG. 1 is a configuration diagram of a leak occurrence position estimation system 100 that estimates a leak occurrence position based on a measurement value of a pressure sensor in a distribution pipe network.
  • the leak occurrence position estimation system 100 includes a leak occurrence position estimation device 101, a leak occurrence event detection device 102, and a measurement value collection device 103.
  • the measurement value collection device 103 collects sensor measurement values from the sensor 191 that measures the state of the target distribution pipe network, and transmits the collected sensor measurement values to the water leakage occurrence position estimation device 101 and the water leakage occurrence event detection device 102.
  • the water leakage occurrence event detection device 102 receives the sensor measurement value from the measurement value collection device 103, detects an abnormality including the occurrence of water leakage based on the received sensor measurement value, and leaks the detected abnormality occurrence information as event data. It transmits to the generation
  • the event data transmitted by the water leakage occurrence event detection device 102 includes information on the type of event, the DMA (District Metered Area) where the abnormality occurred, and the time of occurrence of the abnormality.
  • the event type includes at least the occurrence of water leakage.
  • the water leakage occurrence event detection device 102 uses an arbitrary technique such as statistical machine learning for the abnormality occurrence detection processing.
  • the water leakage occurrence position estimation device 101 includes an event data reception unit 111, a measurement value prediction unit 112, an index value calculation unit 113, a position estimation unit 114, a measurement value collection unit 131, a position display unit 132, and measurement units.
  • Each storage unit includes a value storage unit 121, a sensor arrangement storage unit 122, and a position storage unit 123.
  • the event data reception unit 111 receives event data including a water leakage occurrence event from the water leakage occurrence event detection device 102, and outputs the received event data to the measurement value prediction unit 112. When the event data is water leakage occurrence event data, the event data reception unit 111 outputs the water leakage occurrence event data to the index value calculation unit 113.
  • the measurement value prediction unit 112 reads the measurement value data from the measurement value storage unit 121, inputs event data from the event data reception unit 111, and obtains sensor measurement values in the same water distribution state as the water distribution state from which the measurement value data was acquired.
  • the predicted sensor measurement value (hereinafter, predicted measurement value) is output to the index value calculation unit 113. Details of the prediction process and the water distribution state will be described later with reference to FIG.
  • the index value calculation unit 113 reads the measurement value from the measurement value storage unit 121, inputs the predicted measurement value from the measurement value prediction unit 112, and inputs the water leakage occurrence event data from the event data reception unit 111.
  • the index value calculation unit 113 calculates the water leakage occurrence event data from the measurement values read from the measurement value storage unit 121 and the predicted measurement values input from the measurement value prediction unit 112 of each sensor in the water leakage occurrence DMA included in the event data.
  • the index value of the positional relationship between the water leakage position corresponding to and the position of the sensor in the pipe network is calculated, and the calculated index value is output to the position estimation unit 114. Details of the index value calculation process will be described later with reference to FIG.
  • the position estimation unit 114 inputs the index value of each sensor from the index value calculation unit 113, reads the sensor arrangement information from the sensor arrangement storage unit 122, and the position of the index value relative to the water leakage occurrence position based on the magnitude comparison between the sensors An area that satisfies the relationship is estimated as the position of occurrence of water leakage, and the estimated position of water leakage is stored in the position storage unit 123. Details of the position estimation processing will be described later with reference to FIG.
  • the measurement value storage unit 121 receives the measurement values of a plurality of sensors installed in the distribution pipe network to be managed from the measurement value collection unit 131, and stores the received measurement values.
  • the stored measurement value is read by the index value calculation unit 113 and the measurement value prediction unit 112.
  • the sensor arrangement storage unit 122 stores sensor arrangement information including the position of the sensor installed in the management target water distribution pipe network.
  • the stored sensor arrangement information is read by the position estimation unit 114.
  • the sensor arrangement information includes, in addition to the position of the sensor, for example, the type of sensor such as a pressure sensor, a flow rate sensor, and a listening level sensor, and the installation altitude of the sensor.
  • the position storage unit 123 inputs a water leak occurrence position that is an estimation result of the water leak occurrence position from the position estimation unit 114, and stores the input water leak occurrence position.
  • the stored water leak occurrence position is read by the position display unit 132.
  • the measurement value collection unit 131 receives the measurement values of a plurality of sensors installed in the distribution pipe network to be managed from the measurement value collection device 103 and stores them in the measurement value storage unit 121.
  • the position display unit 132 reads out the leakage occurrence position that is the estimation result of the leakage occurrence position from the position storage unit 123 and presents it to the operator of the leakage occurrence position estimation device 101. For example, the position display unit 132 displays the water leakage occurrence position on the display window to the operator. It is only necessary that the operator can recognize the water leakage occurrence position, and it may be displayed on a smart device such as a smartphone or tablet held by the operator. In addition, when the water leak occurrence position satisfies a preset condition, for example, when a specific area is estimated as the water leak occurrence position, the position display unit 132 performs a push-type notification by e-mail, an alarm, or the like. Also good.
  • FIG. 2 is a hardware configuration diagram of the water leak occurrence position estimation apparatus 101.
  • the leakage occurrence position estimation device 101 is dedicated hardware or a computer to which the CPU 201, the memory 202, the communication control unit 203, the input unit 204, the display unit 205, and the peripheral device IF unit 206 are connected by a bus 210.
  • the CPU 201 executes a program on memory 202.
  • the memory 202 temporarily stores programs, tables, and the like.
  • the communication control unit 203 is connected to the network 220.
  • the input unit 204 is a keyboard, a mouse, or the like.
  • the display unit 205 is the display illustrated in FIG.
  • the peripheral device IF unit 207 is an interface with a printer or the like.
  • the water leak occurrence position estimation device 101 in FIG. 1 is realized by the CPU 201 executing the program of each processing unit.
  • the water leak occurrence position estimation device 101 and the water leak occurrence event detection device 102 may be realized as different programs on one dedicated hardware or computer.
  • FIG. 3 is a configuration example of one DMA in the water distribution pipe network to be managed by the water leak occurrence position estimation apparatus 101.
  • DMA is a part of the water distribution network, and there are few adjacent pipes and pipes with inflow and outflow of water. In many cases, the number of pipes is limited to one. It is measured.
  • the water distribution network is composed of a number of DMAs.
  • FIG. 3 is an example of DMA, and includes pipes such as a water pipe 351.
  • An area 340 indicated by a broken line is a DMA in which the inflow / outflow pipe is limited to one inflow pipe from the reservoir 301 and the flow rate sensor 310 is installed in the inflow pipe, and therefore the area 340 is referred to as a DMA 340.
  • the DMA 340 is provided with a plurality of pressure sensors 320-325.
  • the sensor 191 that collects the measurement value by the measurement value collection device 103 is the flow sensor 310 and the pressure sensor 320-325.
  • the pressure sensor 320 measures the pressure at the point where it flows into the DMA 340.
  • the pressure at the inflow point is referred to as inflow pressure.
  • FIG. 4 is a diagram showing a time-series change (trend) of the measurement value of the sensor accompanying the occurrence of water leakage.
  • FIG. 4 shows a distribution flow rate trend 411, an inflow pressure trend 412, and a pressure 1 trend 413 before and after the occurrence of water leakage.
  • the distribution flow rate trend 411 is a time-series change in the measurement value of the flow rate sensor 310.
  • the inflow pressure trend 412 is a time-series change in the measurement value of the pressure sensor 320.
  • the pressure 1 trend 413 is a time series change of one measurement value of the pressure sensors 321 to 325.
  • the measurement value collection unit 131 collects the measurement values of each sensor at a predetermined cycle, for example, a 1-5 minute cycle, and stores it in the measurement value storage unit 121.
  • the time 461 is the water leakage occurrence time.
  • the water leakage occurrence event detection device 102 detects the occurrence of water leakage and transmits the water leakage occurrence event data to the event data receiving unit 111.
  • the leak occurrence position estimation device 101 starts the execution of the index value calculation unit 113 in response to the input of the leak occurrence event data to the index value calculation unit 113 after the time 462 when the event data reception unit 111 receives the leak occurrence event data. Then, the process of estimating the location of water leakage is started. In response to reception of the leakage occurrence event data by the event data receiving unit 111, it may be said that the leakage occurrence position estimation device 101 starts the estimation process of the leakage occurrence position. Thereafter, the leakage occurrence position estimation device 101 repeats the leakage occurrence position estimation process periodically, for example, at times 463 to 465.
  • FIG. 5 is a diagram illustrating a prediction example of the measurement value before the occurrence of water leakage by the measurement value prediction unit 112.
  • the vertical axis is a measured value of the water distribution flow rate 411 by the flow rate sensor 310.
  • the horizontal axis represents the pressure loss, that is, the difference between the measured value of the inflow pressure 412 by the pressure sensor 320 and the measured value 413 of the pressure sensor.
  • the measurement value prediction unit 112 generates a prediction curve (broken line 511 in FIG. 5) for predicting the measurement value for each sensor and for each time zone between event occurrence times included in each event data.
  • the types of events include, in addition to the occurrence of water leakage, valve opening, valve closing, and water leakage isolation (repair).
  • the actual value (white circle mark) 501 of the relationship between the water distribution flow rate 411 and the pressure loss is the water leakage occurrence time included in a certain water leakage occurrence event data input by the measurement value prediction unit 112 and the event (leakage leakage) Event other than the occurrence event) Indicates the actual value of the time zone between the occurrence time.
  • the actual value (hatched shaded circle) 502 indicates the actual value at the time after the water leakage occurrence time.
  • the actual value 501 is the actual value before the occurrence of water leakage
  • the actual value 502 is the actual value after the occurrence of water leakage.
  • the measured value prediction unit 112 estimates a prediction curve 511 indicating a correlation between the actual values using a plurality of actual values 501 in the time zone before the water leakage occurrence time.
  • the measurement value prediction unit 112 also estimates a prediction curve error range 521-522 indicating an error level of the correlation.
  • a statistical method such as a least square method or a statistical machine learning method is used.
  • the measurement value prediction unit 112 calculates a predicted measurement value in the same water distribution state as the input measurement value.
  • the distribution state is a distribution flow rate including at least the boundary head and the amount of leakage. Since the head of a certain point is the sum of the elevation of the distribution pipe at that point and the value obtained by converting the pressure value in this distribution pipe into the elevation, the boundary head corresponds to the inflow pressure 412 here.
  • the measured value predicting unit 112 displays the actual value 502 on the prediction curve 511.
  • the pressure loss value at the same water distribution flow rate is calculated (the pressure loss corresponding to the water distribution flow rate of the actual value 502 is obtained on the prediction curve 511). Further, a value obtained by subtracting the calculated pressure loss value from the inflow pressure 412 corresponding to the input measurement value is calculated, and this value is set as a predicted measurement value. Further, the measurement value prediction unit 112 similarly calculates the error range of the predicted measurement value using the prediction curve error range.
  • the distribution flow rate of the actual value 502 includes an increase in the flow rate due to water leakage.
  • shaft and a horizontal axis is an example, and if a water distribution flow rate including the boundary head and the amount of water leaks can be made the same, it can use arbitrary prediction models.
  • the measurement value prediction unit 112 generates a prediction curve for each time period between each event occurrence time other than the occurrence of water leakage, so that the water leakage occurrence position estimation device 101 has a plurality of events other than the occurrence of water leakage. Even if water leakage occurs, the location of water leakage can be estimated. For example, if water leakage occurs several hours after opening or closing of a valve due to construction, a more probable location of water leakage can be estimated by generating a prediction curve from actual values after opening or closing of the valve. .
  • FIG. 6 is a diagram showing a geographical change in the index value under the condition that the distribution flow rate including the boundary head and the amount of leakage is the same in the distribution state.
  • FIG. 6 illustrates a one-dimensional example for convenience.
  • the horizontal axis of the graph is the distance from the distribution reservoir 301.
  • the left end is the distribution reservoir 301, and the distance from the distribution reservoir 301 is farther to the right.
  • the outflows 611 and 612 indicate the distribution of outflows that combined demand and water leakage after the occurrence of water leakage and before the occurrence of water leakage. Since the water distribution flow rate including the water leakage amount is the same, the outflow amount at other points decreases due to the increase in the outflow amount at the water leakage occurrence position 650.
  • Water heads 621 and 622 indicate the distribution of water heads after the occurrence of water leakage and before the occurrence of water leakage.
  • the measurement value prediction unit 112 calculates the value at the pressure sensor of the water head 622 before the occurrence of water leakage.
  • a value (difference) obtained by subtracting the water head 622 before the occurrence of water leakage from the water head 621 after the occurrence of water leakage is the water head difference 631. Since the altitude does not change before and after the occurrence of water leakage, the water head difference 631 is equal to the decrease value of the pressure value due to the occurrence of water leakage.
  • the location of water leakage is not a point where the head difference is negative, for example, downstream of the point 664. Further, as the points 661 to 664 are closer to the water leakage occurrence position, the water head difference takes a larger value. The relationship between the distance from the location of the water leakage and the head difference is not affected by the decrease in demand downstream of the water leakage due to the pressure drop caused by water leakage.
  • the index value calculation unit 113 calculates, as an estimated value of the water head difference, a decrease value from the predicted pressure value of the pressure measurement value at the pressure sensor after the water leakage occurrence event as an index value for each pressure sensor in the pipe network. .
  • the index value calculation unit 113 subtracts the predicted pressure value of the pressure sensor input from the measurement value prediction unit 112 from the measurement value of the pressure sensor read from the measurement value storage unit 121 for each pressure sensor in the pipe network, The difference is output to the position estimation unit 114 as an index value.
  • the index value calculation unit 113 receives the error range of the predicted pressure value from the measurement value prediction unit 112 together with the calculation of the index value, and estimates the position where the range obtained by inverting the positive / negative direction of the error range as the error range of the index value Output to the unit 114.
  • the index value calculation unit 113 calculates the index value for each of a plurality of times after the water leakage occurrence time recorded in the received water leakage event data. Then, the index values at a plurality of times may be output to the position estimating unit 114 as time series. As an example of a plurality of times, it is possible to select a time at a predetermined interval (predetermined period) after the water leakage occurrence time 462, as in the times 462 to 465 in FIG.
  • the position estimation unit 114 is an area where the index value is not downstream of the negative pressure sensor.
  • the area near the largest pressure sensor (greater than the index value of other pressure sensors) is estimated as the position of the occurrence of water leakage.
  • the position estimation unit 114 uses a fixed upstream / downstream relationship stored in the sensor arrangement storage unit 122 as the upstream / downstream relationship in the pipe network. For example, the area where the sensor arrangement storage unit 122 further divides the DMA 340 and the upstream / downstream relationship between each area and the sensor are stored in the sensor arrangement storage unit 122, and the position estimation unit 114 is the upstream / downstream between the area and the sensor. Using the relationship, it is possible to limit the area that is estimated as the water leakage occurrence position.
  • the position estimation unit 114 may estimate the water leakage occurrence position by using a size comparison considering the error range using the error range of the index value input from the index value calculation unit 113. Between the pressure sensor of the index value for the determination of the pressure sensor with the negative index value, that is, the comparison of the index value with 0 (determination of positive / negative of the index value) and the identification of the pressure sensor with the maximum index value In the size comparison, the position estimation unit 114 uses, for example, a size comparison considering the following error range.
  • the size comparison between the pressure sensors by the position estimation unit 114 will be described taking the comparison of index values of the two pressure sensors A and B as an example.
  • the position estimating unit 114 determines that the index value of the pressure sensor A is greater than the index value of the pressure sensor B when the error upper limit value of the index value of the pressure sensor B is smaller than the error lower limit value of the index value of the pressure sensor A. judge.
  • the position estimation unit 114 estimates the area where the areas estimated from the index values at the plurality of times overlap as the water leakage occurrence position. Good.
  • the position estimation unit 114 can estimate the leak occurrence position from the index value at each time, but by narrowing down to the area estimated as the leak occurrence position at a plurality of times, the position estimation unit 114 is more accurate based on the transition of the index value over a long period of time. The location of water leakage can be estimated.
  • the index value calculation unit 113 and the position estimation unit 114 compare the water head difference on the condition that the water flow rate including the water leakage amount is the same, and thus the influence of the change in the demand amount is affected. Even if it exists, the relationship between a water leak generation position and a sensor position can be specified.
  • Fig. 7 shows an example of the screen display of the location of water leakage.
  • the position display unit 132 displays the water leakage occurrence position on the screen.
  • a position display window 701 displayed on the display or the like by the position display unit 132 includes a DMA selection box 702, a position display panel 703, and an index value display panel 704.
  • the position display unit 132 displays the leak occurrence position related to the DMA included in the leak occurrence event data received by the leak occurrence position estimation apparatus 101. However, the position display unit 132 changes the DMA to be displayed on the position display panel 703 when the operator operates the DMA selection box 702.
  • the position display unit 132 displays the water leak occurrence position estimated by the position estimation unit 114 on the position display panel 703. In FIG. 7, the water leakage occurrence position is displayed as a shaded portion of the region 761.
  • the position display unit 132 displays the index value of each sensor calculated by the index value calculation unit 113 on the index value display panel 704. In FIG. 7, the sensor ID, the index value, and the sensor type are displayed.
  • the leakage occurrence position estimation device 101 can estimate the leakage occurrence position and present the estimated leakage occurrence position to the operator even if there is an influence of a state change due to other factors including a change in demand.
  • a water leak occurrence position estimation device for estimating the upstream / downstream relationship between adjacent sensors and areas from the pressure sensor measurement value, sensor installation altitude value, and pipe network connection relationship between adjacent sensors will be described. Since most of the configuration and processing of the water leakage occurrence position estimation device are the same as those in the first embodiment, the differences will be mainly described.
  • FIG. 8 is a configuration example of one DMA having a plurality of inflow paths in the water distribution pipe network that is a target of water leakage occurrence position estimation. Compared with the DMA of FIG. 3, a reservoir 802, a flow sensor 811 and a pressure sensor 829 are added.
  • the flow direction in the pipe network changes depending on the distribution state of the boundary head and the distribution flow rate.
  • the tendency for the upstream / downstream relationship between regions to change will intensify.
  • the position estimation unit 114 of the present embodiment additionally inputs the sensor installation altitude value and the pipe network connection relationship between the adjacent sensors from the sensor arrangement storage unit 122, and the upstream / downstream relationship between the adjacent sensor and the region, It is estimated from the pressure sensor measurement value, the sensor installation altitude value, and the pipe network connection relationship between adjacent sensors.
  • the upstream / downstream relationship between the sensors can be estimated by calculating the head value of the sensor position from the pressure sensor measurement value.
  • the position estimation unit 114 adds the sensor installation altitude value to the value obtained by converting the pressure sensor measurement value into the altitude value, and calculates the head value of each sensor. Subsequently, the position estimation unit 114 determines which one of the adjacent sensors connected by the pipe network is upstream or downstream from the comparison of the head values between the sensors. From the upstream / downstream relationship of each sensor and the pipe connection relationship, the upstream / downstream relationship is also estimated for the area near the sensor.
  • the position estimation unit 114 estimates a water leak occurrence position based on the estimated upstream / downstream relationship. Since this estimation process is the same as that of Example 1, description is abbreviate
  • the leak occurrence position estimation device 101 can estimate the leak occurrence position with higher accuracy even in a DMA having a plurality of inflow paths.
  • a water leakage occurrence position estimation device for managing a water distribution pipe network in which an acoustic level sensor is installed in addition to a pressure sensor will be described. Since most of the configuration and processing of the water leakage occurrence position estimation device are the same as those in the first embodiment, the differences will be mainly described.
  • FIG. 9 is a configuration example of one DMA in which an acoustic level sensor is installed in a distribution pipe network that is a target of water leakage occurrence position estimation. Compared with the DMA of FIG. 3, acoustic level sensors 941-944 are added.
  • the measurement value prediction unit 112 calculates, for each acoustic level sensor, measurement value data (acoustic level) in the same time zone on a different day from the measurement value data after the occurrence of water leakage among the measurement value data before the occurrence of water leakage.
  • the measured value is output to the index value calculation unit 113.
  • the index value calculation unit 113 outputs a value obtained by subtracting the input predicted measurement value from the input measurement value for each acoustic level sensor to the position estimation unit 114 as an index value.
  • the position estimation unit 114 estimates the area near the acoustic level sensor having the maximum index value as the position of the water leakage occurrence within the range of the water leakage occurrence position estimated from the pressure sensor information. And stored in the position storage unit 123.
  • the water leak occurrence position estimation device 101 can estimate the water leak occurrence position with higher accuracy by using the measurement value of the acoustic level sensor in addition to the pressure sensor.
  • the embodiment of the present invention is not limited to the above-described examples, and includes various modifications.
  • the above-described embodiments have been described in detail for easy understanding of the present invention, and are not necessarily limited to those having all the configurations described.
  • a part of the configuration of one embodiment can be replaced with the configuration of another embodiment, and the configuration of another embodiment can be added to the configuration of one embodiment.
  • Each of the above-described configurations, functions, processing units, processing means, and the like may be realized by hardware by designing a part or all of them with, for example, an integrated circuit.
  • Each of the above-described configurations, functions, and the like may be realized by software by interpreting and executing a program that realizes each function by the processor.
  • Information such as programs, tables, and files that realize each function can be stored in a memory, a hard disk, a recording device such as an SSD (Solid State Drive), or a recording medium such as an IC card, an SD card, or a DVD.
  • SSD Solid State Drive
  • 100 Water leakage occurrence position estimation system
  • 101 Water leakage occurrence position estimation apparatus
  • 102 Water leakage occurrence event detection apparatus
  • 103 Measurement value collection apparatus
  • 111 Event data reception section
  • 112 Measurement value prediction section
  • 113 Index value calculation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Hydrology & Water Resources (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Examining Or Testing Airtightness (AREA)

Abstract

開示する漏水発生位置推定装置は、配水管網に配置された複数のセンサの第１の計測値を収集する計測値収集部、センサの第１の計測値に基づいて、所定の配水状態における予測計測値を予測する計測値予測部、漏水発生をイベントとするイベントデータの受信に応答して、漏水発生の後のセンサの第２の計測値と、計測値予測部によって予測された、第２の計測値を計測した配水状態と同じ配水状態における予測計測値とから指標値を計算する指標値計算部、および、複数のセンサの配置と、複数のセンサ間の指標値の大小比較に基づく漏水発生位置との位置関係を満足する地域を漏水発生の位置として推定する位置推定部を有する。

Description

漏水発生位置推定装置、システムおよび方法
　本発明は、配水管網中の漏水発生位置の推定装置、システムおよび方法に関する。
　特許文献１は、「異常の地理的位置を決定する計算方法で、前記計算方法は、地区またはゾーン内での異常の確からしい地理的位置を一つ統計的に決定するよう設計されたテストを複数実行し、なお、各テストの実行は結果一つを生成し、また、前記計算方法は、複数のテストの結果を統合して異常の決定された確からしい位置に関するスコアを生成する」と記載している。
ＵＳ　Ｐａｔｅｎｔ　８，５８３，３８６
　特許文献１の方法で漏水発生を異常としてその位置を決定する場合、漏水発生とほぼ同時に他要因による配水状態の変化があると、漏水発生による配水状態の変化に加えて他要因による配水状態の変化も圧力計などのメータデータに影響する。したがって、特許文献１の方法では、他要因による配水状態の変化の影響を受けて漏水発生位置として誤った位置を提示してしまう。他要因による配水状態の変化の例としては、漏水発生の影響による需要量の減少がある。
　そこで、本発明では、需要量変化を含む他要因による配水状態の変化の影響があっても、漏水発生位置を推定できる漏水発生位置推定装置システムおよび方法を提供する。
　開示する漏水発生位置推定装置は、配水管網に配置された複数のセンサの第１の計測値を収集する計測値収集部、センサの第１の計測値に基づいて、所定の配水状態における予測計測値を予測する計測値予測部、漏水発生をイベントとするイベントデータの受信に応答して、漏水発生の後のセンサの第２の計測値と、計測値予測部によって予測された、第２の計測値を計測した配水状態と同じ配水状態における予測計測値とから指標値を計算する指標値計算部、および、複数のセンサの配置と、複数のセンサ間の指標値の大小比較に基づく漏水発生位置との位置関係を満足する地域を漏水発生の位置として推定する位置推定部を有する。
　開示する漏水発生位置推定装置によれば、他要因による配水状態の変化があっても漏水発生位置を推定できる。
漏水発生位置推定システムの構成図である。 漏水発生位置推定装置のハードウェア構成図である。 ＤＭＡの構成例である。 漏水発生にともなうセンサの計測値の時系列変化を示す図である。 漏水発生前の計測値の予測例を示す図である。 境界水頭および漏水量を含む配水流量の同じ条件における、指標値の地理的変化を示す図である。 漏水発生位置の画面表示の例である。 複数の流入経路を有するＤＭＡの構成例である。 音響レベルセンサが設置されたＤＭＡの構成例である。
　以下、実施例を図面を用いて説明する。
　図１は、配水管網中の圧力センサの計測値にもとづいて漏水発生位置を推定する漏水発生位置推定システム１００の構成図である。漏水発生位置推定システム１００は、漏水発生位置推定装置１０１、漏水発生イベント検知装置１０２、および計測値収集装置１０３を含む。
　計測値収集装置１０３は、対象の配水管網の状態を計測するセンサ１９１からセンサ計測値を収集し、収集したセンサ計測値を漏水発生位置推定装置１０１および漏水発生イベント検知装置１０２に送信する。
　漏水発生イベント検知装置１０２は、計測値収集装置１０３からセンサ計測値を受信し、受信したセンサ計測値に基づいて漏水発生を含む異常発生を検知し、検知した異常発生の情報をイベントデータとして漏水発生位置推定装置１０１に送信する。漏水発生イベント検知装置１０２が送信するイベントデータは、イベントの種類、異常が発生したＤＭＡ(District Metered Area)および異常発生の時刻に関する情報を含む。また、イベントの種類には、少なくとも漏水発生を含む。漏水発生イベント検知装置１０２は、異常発生検知処理に、統計的機械学習などの任意の技術を利用する。
　漏水発生位置推定装置１０１は、イベントデータ受信部１１１、計測値予測部１１２、指標値計算部１１３、位置推定部１１４、計測値収集部１３１、および位置表示部１３２の各処理部、並びに、計測値記憶部１２１、センサ配置記憶部１２２、および位置記憶部１２３の各記憶部、を有する。
　イベントデータ受信部１１１は、漏水発生イベント検知装置１０２より、漏水発生イベントを含むイベントデータを受信し、受信したイベントデータを計測値予測部１１２に出力する。イベントデータ受信部１１１は、イベントデータが漏水発生イベントデータである場合は、漏水発生イベントデータを指標値計算部１１３に出力する。
　計測値予測部１１２は、計測値記憶部１２１から計測値データを読み出し、イベントデータ受信部１１１からのイベントデータを入力し、計測値データを取得した配水状態と同じの配水状態におけるセンサ計測値を予測し、予測したセンサ計測値（以下、予測計測値）を指標値計算部１１３に出力する。予測処理、および、配水状態の詳細は、図５の説明で後述する。
　指標値計算部１１３は、計測値記憶部１２１から計測値を読み出し、計測値予測部１１２からの予測計測値を入力し、イベントデータ受信部１１１からの漏水発生イベントデータを入力する。指標値計算部１１３は、イベントデータに含まれる漏水発生ＤＭＡ内の各センサの、計測値記憶部１２１から読み出した計測値と計測値予測部１１２から入力した予測計測値とから、漏水発生イベントデータに対応する漏水位置と管網中のセンサの位置との位置関係の指標値を計算し、求めた指標値を位置推定部１１４に出力する。指標値の計算処理の詳細は、図６の説明で後述する。
　位置推定部１１４は、指標値計算部１１３からの各センサの指標値を入力し、センサ配置記憶部１２２からセンサ配置情報を読み出し、指標値のセンサ間の大小比較に基づく漏水発生位置との位置関係を満足する地域を漏水発生の位置として推定し、推定した漏水発生位置を位置記憶部１２３に格納する。位置推定処理の詳細は、図６の説明で後述する。
　計測値記憶部１２１は、計測値収集部１３１から管理対象の配水管網に設置された複数のセンサの計測値を受信し、受信した計測値を格納する。格納した計測値は、指標値計算部１１３および計測値予測部１１２により読み出される。
　センサ配置記憶部１２２は、管理対象の配水管網に設置されたセンサの位置を含むセンサ配置情報を格納する。格納したセンサ配置情報は、位置推定部１１４により読み出される。センサ配置情報は、センサの位置に加えて、例えば圧力センサ、流量センサ、音聴レベルセンサのようなセンサの種別と、センサの設置標高等を含む。
　位置記憶部１２３は、位置推定部１１４から漏水発生位置の推定結果である漏水発生位置を入力し、入力した漏水発生位置を格納する。格納した漏水発生位置は、位置表示部１３２により読み出される。
　計測値収集部１３１は、計測値収集装置１０３から管理対象の配水管網に設置された複数のセンサの計測値を受信し、計測値記憶部１２１に格納する。
　位置表示部１３２は、位置記憶部１２３から漏水発生位置の推定結果である漏水発生位置を読み出し、漏水発生位置推定装置１０１の操作者に提示する。例えば、位置表示部１３２は、ディスプレイのウィンドウに漏水発生位置を操作者に表示する。操作者が漏水発生位置を認識できればよく、操作者の持つスマートフォンやタブレット等のスマートデバイス上に表示してもよい。また、位置表示部１３２は、漏水発生位置が予め設定された条件を満たす場合、例えば、特定の地域が漏水発生位置と推定された場合には、メールやアラーム等でプッシュ型の通知を行ってもよい。
　図２は、漏水発生位置推定装置１０１のハードウェア構成図である。漏水発生位置推定装置１０１は、ＣＰＵ２０１、メモリ２０２、通信制御部２０３、入力部２０４、表示部２０５、および周辺機器ＩＦ部２０６が、バス２１０により接続される専用ハードウェアまたはコンピュータである。
　ＣＰＵ２０１は、メモリ２０２上のプログラムを実行する。メモリ２０２は、プログラム、テーブル等を一時記憶する。通信制御部２０３は、ネットワーク２２０と接続されている。入力部２０４は、キーボード、マウス等である。表示部２０５は、図１の説明のディスプレイである。周辺機器ＩＦ部２０７は、プリンタ等とのインタフェースである。
　図１と図２との対比から明らかなように、図１の漏水発生位置推定装置１０１は、ＣＰＵ２０１が各処理部のプログラムを実行することで実現している。漏水発生位置推定装置１０１と漏水発生イベント検知装置１０２を、一つの専用ハードウェアやコンピュータ上に異なるプログラムとして実現してもよい。
　図３は、漏水発生位置推定装置１０１の管理対象となる配水管網のうち１つのＤＭＡの構成例である。一般的に、ＤＭＡは、配水管網の一部分で、隣接する管網と水の流入出のある管が少数、多くの場合は１つに限定されており、また全ての流入出管で流量が計測される。配水管網は多数のＤＭＡから構成される。
　図３は、ＤＭＡの一例であり、配水管３５１等の管を含む。破線で示す区域３４０は、流入出管が配水池３０１からの１つの流入管に限定され、その流入管には流量センサ３１０が設置されたＤＭＡであるので、区域３４０をＤＭＡ３４０と呼ぶ。ＤＭＡ３４０には、複数の圧力センサ３２０－３２５が設置されている。
　計測値収集装置１０３が計測値を収集するセンサ１９１は、ＤＭＡ３４０の例では、流量センサ３１０および圧力センサ３２０－３２５である。圧力センサ３２０は、ＤＭＡ３４０に流入する地点での圧力を計測する。流入地点の圧力を、以下、流入圧力と呼ぶ。圧力センサ３２１－３２５のように、ＤＭＡの流入地点ではなく、ＤＭＡの内部に設置されているセンサを、管網中のセンサを呼ぶ。
　図４は、漏水発生にともなうセンサの計測値の時系列変化（トレンド）を示す図である。図４を参照して、計測値記憶部１２１に格納されるセンサの計測値の例と、漏水発生位置推定処理の開始について説明する。図４は、漏水発生の前後における配水流量トレンド４１１、流入圧力トレンド４１２、圧力１トレンド４１３を示している。配水流量トレンド４１１は、流量センサ３１０の計測値の時系列変化である。流入圧力トレンド４１２は、圧力センサ３２０の計測値の時系列変化である。圧力１トレンド４１３は、圧力センサ３２１－３２５の一つの計測値の時系列変化である。
　計測値収集部１３１は、図４のように、所定の周期、例えば１－５分周期で、各センサの計測値を収集し、計測値記憶部１２１に格納する。図４の例では、時刻４６１が漏水発生時刻である。時刻４６２に漏水発生イベント検知装置１０２が漏水発生を検知し、漏水発生イベントデータをイベントデータ受信部１１１に送信する。
　漏水発生位置推定装置１０１は、イベントデータ受信部１１１が漏水発生イベントデータを受信した時刻４６２以降に、指標値計算部１１３に漏水発生イベントデータの入力に応答した指標値計算部１１３の実行開始をもって、漏水発生位置の推定処理を始める。漏水発生イベントデータのイベントデータ受信部１１１による受信に応答して、漏水発生位置推定装置１０１は漏水発生位置の推定処理を始めると言ってもよい。漏水発生位置推定装置１０１は、以降、例えば時刻４６３－４６５のように定期的に、漏水発生位置の推定処理を繰り返す。
　図５は、計測値予測部１１２による、漏水発生前の計測値の予測例を示す図である。縦軸は、流量センサ３１０による配水流量４１１の計測値である。また、横軸は圧力損失、すなわち、圧力センサ３２０による流入圧力４１２の計測値と、圧力センサの計測値４１３との差である。
　計測値予測部１１２は、センサごとに、また各イベントデータに含まれるイベント発生時刻の間の時間帯ごとに、計測値を予測する予測カーブ（図５の破線５１１）を生成する。イベントの種別は、漏水発生に加えて、バルブ開放、バルブ閉止、漏水の隔離（修理）などがある。
　図５において、配水流量４１１と圧力損失との関係の実績値（白丸印）５０１は、計測値予測部１１２が入力したある漏水発生イベントデータに含まれる漏水発生時刻と、その直前のイベント（漏水発生イベント以外のイベント）発生時刻との間の時間帯の実績値を表す。また、実績値（斜線網掛けの丸印）５０２は、漏水発生時刻後の時刻の実績値を示す。換言すると、実績値５０１は漏水発生より前の実績値であり、実績値５０２は、漏水発生後の実績値である。
　計測値予測部１１２は、漏水発生時刻の前の時間帯の複数の実績値５０１を用いて、実績値間の相関関係を示す予測カーブ５１１を推定する。また、計測値予測部１１２は、予測カーブ５１１に加えて、相関関係の誤差レベルを示す予測カーブ誤差範囲５２１－５２２も推定する。予測カーブ５１１および予測カーブ誤差範囲の推定処理には、例えば最小二乗法などの統計的手法や統計的機械学習の手法を利用する。
　続いて、計測値予測部１１２は、入力した計測値と同じ配水状態における予測計測値を計算する。配水状態は、少なくとも境界水頭、および漏水量を含む配水流量である。ある地点の水頭は、その地点の配水管の標高と、この配水管内の圧力値を標高に換算した値の和であるため、ここでは、境界水頭は流入圧力４１２に相当する。
　漏水発生後の実績値５０２を、入力した計測値から計算された点（配水流量と圧力損失の関係を示す点）とすると、計測値予測部１１２は、予測カーブ５１１上で、実績値５０２と同じ配水流量における圧力損失値を算出する（予測カーブ５１１上で、実績値５０２の配水流量に対応する圧力損失を求める。）。さらに、入力した計測値に対応する流入圧力４１２から算出した圧力損失値を引いた値を計算し、この値を予測計測値とする。また、計測値予測部１１２は、予測カーブ誤差範囲を用いて、予測計測値の誤差範囲を同様に計算する。なお、実績値５０２の配水流量は、漏水による流量の増加分を含んでいる。また、上記の縦軸・横軸の選び方は例であり、境界水頭および漏水量を含む配水流量を同じにして予測ができれば、任意の予測モデルを利用することができる。
　なお、計測値予測部１１２が漏水発生以外の各イベント発生時刻の間の時間帯ごとに予測カーブを生成することで、漏水発生位置推定装置１０１は、複数の漏水発生以外のイベントが生じた後に漏水が発生したとしても、漏水発生位置を推定できる。例えば、工事に伴うバルブの開放や閉止があった数時間後に、漏水が発生した場合、バルブの開放や閉止の後の実績値から予測カーブを生成することでより確からしい漏水発生位置を推定できる。
　図６は、配水状態である、境界水頭および漏水量を含む配水流量が同じ条件における、指標値の地理的変化を示す図である。現実の配水管網は、二次元的に広がっているが、図６では便宜上、一次元的な例で説明する。グラフ横軸は、配水池３０１からの距離であり、左端が配水池３０１で、右にいくほど配水池３０１からの距離が遠い地点であることを示す。
　流出量６１１および６１２は、漏水発生後および漏水発生前の、需要と漏水を合わせた流出量の分布を示す。漏水量を含む配水流量が同じとしているため、漏水発生位置６５０での流出量増加により、それ以外の地点での流出量は減少する。
　水頭６２１および６２２は、漏水発生後および漏水発生前の、水頭の分布を示す。計測値予測部１１２は、漏水発生前の水頭６２２の圧力センサでの値を計算している。
　漏水発生後の水頭６２１から、漏水発生前の水頭６２２を引いた値（差）が、水頭差６３１である。漏水発生前後で、標高は変化しないため、水頭差６３１は、圧力値の漏水発生による減少値と等しい。
　漏水による流出量の変化により、漏水発生位置は、水頭差が負の地点、例えば地点６６４の下流側ではない。また、地点６６１－６６４のように、漏水発生位置に近いほど、水頭差が大きい値をとる。漏水発生位置との距離と水頭差の関係は、漏水を原因とする圧力低下により、漏水下流で需要が減少することの影響を受けない。
　指標値計算部１１３は、水頭差の推定値として、管網中の各圧力センサについて、漏水発生イベント後の圧力センサでの圧力計測値の、予測圧力値からの減少値を指標値として計算する。
　すなわち、指標値計算部１１３は、管網中の各圧力センサについて、計測値記憶部１２１から読み出した圧力センサの計測値から、計測値予測部１１２から入力した圧力センサの予測圧力値を引き、その差を指標値として位置推定部１１４に出力する。
　また、指標値計算部１１３は、指標値の計算と共に、計測値予測部１１２より予測圧力値の誤差範囲を入力し、誤差範囲の正負方向を反転した範囲を、指標値の誤差範囲として位置推定部１１４に出力する。
　また、指標値計算部１１３は、最新の圧力計測値にもとづいて指標値を計算するのに加えて、受信した漏水イベントデータに記録された漏水発生時刻以降の複数の時刻についてそれぞれ指標値を計算し、複数の時刻の指標値を時系列として位置推定部１１４に出力してもよい。複数の時刻の例としては、図4の時刻４６２－４６５の各時刻のように、漏水発生時刻４６２以降の所定間隔（所定周期）での時刻を選ぶことができる。
　位置推定部１１４は、複数のセンサの配置と、複数のセンサの指標値の、センサ間の大小比較とに基づいて、指標値が負の圧力センサの下流ではない地域で、指標値がＤＭＡ内で最大（他の圧力センサの指標値より大）の圧力センサの近隣である地域を漏水発生の位置として推定する。
　位置推定部１１４は、管網中の上下流関係として、センサ配置記憶部１２２に格納された固定的な上下流関係を利用する。例えば、センサ配置記憶部１２２がＤＭＡ３４０を更に細かく分割したエリアと、各エリアとセンサ間の上下流関係をセンサ配置記憶部１２２に格納しておき、位置推定部１１４はエリアとセンサ間の上下流関係を用いて、漏水発生位置と推定するエリアを限定することができる。
　位置推定部１１４は、指標値計算部１１３から入力した指標値の誤差範囲を用いて、誤差範囲を考慮した大小比較により漏水発生位置を推定してもよい。指標値が負の圧力センサの判定、すなわち、指標値の０との大小比較（指標値の正負の判定）と、指標値が最大の圧力センサの特定のための、指標値の圧力センサ間での大小比較において、位置推定部１１４は、例えば下記の誤差範囲を考慮した大小比較を用いる。
　位置推定部１１４による圧力センサ間の大小比較を、２つの圧力センサＡおよびＢの指標値の比較を例にとって説明する。位置推定部１１４は、圧力センサＡの指標値の誤差下限値よりも、圧力センサＢの指標値の誤差上限値が小さいときに、圧力センサＡの指標値が圧力センサＢの指標値より大きいと判定する。
　位置推定部１１４は、指標値計算部１１３から複数の時刻の指標値を時系列で入力する場合、複数の時刻の指標値からそれぞれ推定される地域の重なる地域を漏水発生位置として推定してもよい。
　位置推定部１１４は、各時刻の指標値からでも漏水発生位置を推定できるが、複数の時刻で漏水発生位置として推定される地域に絞り込むことで、長期間の指標値の推移に基づいたより正確な漏水発生位置を推定できる。
　漏水発生位置推定装置１０１は、指標値計算部１１３および位置推定部１１４が、漏水量をふくむ配水流量が同じであることを条件にして水頭差を比較することで、需要量の変化の影響があっても漏水発生位置とセンサ位置との関係を特定できる。
　図７は、漏水発生位置の画面表示の例である。位置表示部１３２は、漏水発生位置を画面表示する。位置表示部１３２がディスプレイ等に表示する位置表示ウィンドウ７０１は、ＤＭＡ選択ボックス７０２、位置表示パネル７０３、指標値表示パネル７０４を有する。
　位置表示部１３２は、漏水発生位置推定装置１０１が受信した漏水発生イベントデータに含まれるＤＭＡに関する漏水発生位置を表示する。ただし、位置表示部１３２は、操作者がＤＭＡ選択ボックス７０２を操作することによって、位置表示パネル７０３に表示する対象のＤＭＡを変更する。
　位置表示部１３２は、位置表示パネル７０３に、位置推定部１１４が推定した漏水発生位置を表示する。図７では、領域７６１の網掛け部として漏水発生位置を表示している。位置表示部１３２は、指標値表示パネル７０４に、指標値計算部１１３が計算した各センサの指標値を表示する。図７では、センサＩＤ、指標値、センサ種別を表示している。
　以上のように、漏水発生位置推定装置１０１は、需要量変化を含む他要因による状態変化の影響があっても、漏水発生位置を推定し、推定した漏水発生位置を操作者に提示できる。
　本実施例では、圧力センサ計測値、センサ設置標高値、および隣接センサ間の管網接続関係から、隣接センサおよび地域間の上下流関係を推定する漏水発生位置推定装置について説明する。漏水発生位置推定装置の構成と処理の大部分は実施例１と同様であるので、相違点を中心に説明する。
　図８は、漏水発生位置推定の対象となる配水管網のうち複数の流入経路を有する１つのＤＭＡの構成例である。図３のＤＭＡと比較して、配水池８０２、流量センサ８１１、圧力センサ８２９が追加されている。
　図８のように、複数の流入経路（配水池）を有するＤＭＡ３４０の場合、漏水発生に加えて、境界水頭や配水流量の配水状態に依存して管網中の流向が変化し、隣接センサおよび地域間の上下流関係が変化する傾向が強まる。
　そこで、本実施例の位置推定部１１４は、センサ設置標高値、および隣接センサ間の管網接続関係をセンサ配置記憶部１２２から追加して入力し、隣接センサおよび地域間の上下流関係を、圧力センサ計測値、センサ設置標高値、および隣接センサ間の管網接続関係から推定する。
　図６の例のように、上流から下流へ向けて水頭は減少するので、圧力センサ計測値からセンサ位置の水頭値を計算することで、センサ間の上下流関係を推定できる。
　位置推定部１１４は、圧力センサ計測値を標高値に換算した値に、センサ設置標高値を加算し、各センサの水頭値を計算する。続いて位置推定部１１４は、センサ間の水頭値の大小比較から、管網で接続された隣接センサとのいずれが上流または下流かを判定する。各センサの上下流関係と、管路接続関係から、センサ付近の地域についても上下流関係を推定する。
　位置推定部１１４は、推定した上下流関係にもとづいて、漏水発生位置を推定する。この推定処理は、実施例１と同様のため、説明は省略する。
　上下流関係をセンサ計測値から動的に推定することで、漏水発生位置推定装置１０１は、複数の流入経路を有するＤＭＡなどであっても、より精度よく漏水発生位置を推定できる。
　本実施例では、圧力センサに加えて、音響レベルセンサが設置された配水管網を管理対象とする漏水発生位置推定装置について説明する。漏水発生位置推定装置の構成と処理の大部分は実施例１と同様であるので、相違点を中心に説明する。
　図９は、漏水発生位置推定の対象となる配水管網のうち、音響レベルセンサの設置された１つのＤＭＡの構成例である。図３のＤＭＡと比較して、音響レベルセンサ９４１－９４４が追加されている。
　計測値予測部１１２は、各音響レベルセンサについて、漏水発生前の計測値データのうち、漏水発生後の計測値データと異なる日の同じ時間帯の計測値データ（音響レベル）を計算し、予測計測値として指標値計算部１１３に出力する。
　指標値計算部１１３は、各音響レベルセンサについて、入力した計測値から、入力した予測計測値を引いた値を、指標値として位置推定部１１４に出力する。
　位置推定部１１４は、実施例１のように、圧力センサの情報から推定した漏水発生位置の範囲内で、さらに、指標値が最大の音響レベルセンサの近隣である地域を漏水発生の位置として推定し、位置記憶部１２３に格納する。
　漏水発生位置推定装置１０１は、圧力センサに加えて音響レベルセンサの計測値を活用することで、より精度の高い漏水発生位置を推定できる。
　なお、本発明の実施形態は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。また、上記の各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路で設計する等によりハードウェアで実現してもよい。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによりソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリや、ハードディスク、ＳＳＤ（Ｓｏｌｉｄ　Ｓｔａｔｅ　Ｄｒｉｖｅ）等の記録装置、または、ＩＣカード、ＳＤカード、ＤＶＤ等の記録媒体に置くことができる。
　１００：漏水発生位置推定システム、１０１：漏水発生位置推定装置、１０２：漏水発生イベント検知装置、１０３：計測値収集装置、１１１：イベントデータ受信部、１１２：計測値予測部、１１３：指標値計算部、１１４：位置推定部、１３１：計測値収集部、１３２：位置表示部。

Claims (15)

  1. 　配水管網に配置された複数のセンサの第１の計測値を収集する計測値収集部、
    　前記センサの前記第１の計測値に基づいて、所定の配水状態における予測計測値を予測する計測値予測部、
    　漏水発生をイベントとするイベントデータの受信に応答して、前記漏水発生の後の前記センサの第２の計測値と、前記計測値予測部によって予測された、前記第２の計測値を計測した配水状態と同じ配水状態における前記予測計測値とから指標値を計算する指標値計算部、および
    　複数の前記センサの配置と、複数の前記センサ間の前記指標値の大小比較に基づく漏水発生位置との位置関係を満足する地域を前記漏水発生の位置として推定する位置推定部を有することを特徴とする漏水発生位置推定装置。
  2. 　請求項１に記載の漏水発生位置推定装置であって、
    　前記配水状態は、少なくとも境界水頭および配水流量で表すことを特徴とする漏水発生位置推定装置。
  3. 　請求項２に記載の漏水発生位置推定装置であって、
    　複数の前記センサは複数の圧力センサであり、
    　前記指標値計算部は、前記漏水発生の後の前記各圧力センサの圧力計測値の、前記予測計測値である予測圧力値からの減少値を、前記各圧力センサの前記指標値とし、
    　前記位置推定部は、前記配水管網の、前記指標値が負である、複数の前記圧力センサのうちの第１の圧力センサの下流ではない地域で、前記指標値が、複数の前記センサのうちの他の圧力センサの前記指標値より大である第２の圧力センサの近隣である地域を前記漏水発生の位置として推定することを特徴とする漏水発生位置推定装置。
  4. 　請求項２に記載の漏水発生位置推定装置であって、
    　前記指標値計算部は、さらに前記指標値の誤差上限値と誤差下限値を計算し、
    　前記位置推定部は、複数の前記センサの第１のセンサと第２のセンサについて、前記第１のセンサの前記指標値の前記誤差下限値よりも、前記第２のセンサの前記指標値の前記誤差上限値が小さいときに、前記第１のセンサの前記指標値が前記第２のセンサの前記指標値より大きいとする判定に基づいて前記漏水発生の位置を推定することを特徴とする漏水発生位置推定装置。
  5. 　請求項２に記載の漏水発生位置推定装置であって、
    　前記指標値計算部は、前記漏水イベントデータに含まれる前記漏水発生の時刻以降の複数の時刻において前記指標値を計算し、
    　前記位置推定部が、複数の前記時刻の前記指標値のそれぞれから推定される複数の前記地域の重なる地域を前記漏水発生の位置として推定することを特徴とする漏水発生位置推定装置。
  6. 　請求項３に記載の漏水発生位置推定装置であって、
    　前記位置推定部は、前記圧力センサ計測値、前記圧力センサの設置標高値、および隣接する前記圧力センサとの間の管網接続関係から、隣接する前記圧力センサおよび前記地域間の上下流関係を推定し、推定した前記上下流関係に基づいて前記漏水発生の位置を推定することを特徴とする漏水発生位置推定装置。
  7. 　請求項３に記載の漏水発生位置推定装置であって、
    　音響レベルセンサが前記配水管網に配置され、
    　前記指標値計算部は、前記漏水発生の後の前記音響レベルセンサの音響レベル計測値の、前記予測計測値である予測音響レベル値からの増加値を、前記音響レベルセンサの前記指標値とし、
    　前記位置推定部は、前記漏水発生の位置として、前記圧力センサの前記指標値に基づいて推定した前記地域の中で、前記音響レベルセンサの前記指標値が最大の音響レベルセンサの近隣である地域を前記漏水発生の位置として推定することを特徴とする漏水発生位置推定装置。
  8. 　配水管網の地区内における漏水発生の位置を推定する漏水発生位置推定システムであって、
    　前記配水管網に設置された複数のセンサの計測値を収集する計測値収集装置、
    　前記計測値収集装置から受信した前記計測値から漏水発生をイベントとして検知する漏水発生イベント検知装置、並びに、
    　前記配水管網に配置された前記センサの第１の計測値を収集する計測値収集部、
    　前記センサの前記第１の計測値に基づいて、所定の配水状態における予測計測値を予測する計測値予測部、
    　前記漏水発生のイベントのイベントデータの受信に応答して、前記漏水発生の後の前記センサの第２の計測値と、前記計測値予測部によって予測された、前記第２の計測値を計測した配水状態と同じ配水状態における前記予測計測値とから指標値を計算する指標値計算部、および
    　複数の前記センサの配置と、複数の前記センサ間の前記指標値の大小比較に基づく漏水発生位置との位置関係を満足する地域を前記漏水発生の位置として推定する位置推定部を含む漏水発生位置推定装置を有する漏水発生位置推定システム。
  9. 　配水道網の地区内における漏水発生の位置を推定する漏水発生位置推定装置における漏水発生位置推定方法であって、前記漏水発生位置推定装置は、
    　前記配水管網に設置された複数のセンサの第１の計測値を収集し、
    　前記センサの前記第１の計測値に基づいて、所定の配水状態における予測計測値を予測し、
    　漏水発生をイベントとするイベントデータの受信に応答して、前記漏水発生の後の前記センサの第２の計測値と、前記第２の計測値を計測した配水状態と同じ配水状態において予測された前記予測計測値とから指標値を計算し、
    　複数の前記センサの配置と、複数の前記センサ間の前記指標値の大小比較に基づく漏水発生位置との位置関係を満足する地域を前記漏水発生の位置として推定することを特徴とする漏水発生位置推定方法。
  10. 　請求項９に記載の漏水発生位置推定方法であって、
    　前記配水状態は、少なくとも境界水頭および配水流量で表すことを特徴とする漏水発生位置推定方法。
  11. 　請求項１０に記載の漏水発生位置推定方法であって、
    　複数の前記センサは複数の圧力センサであり、
    　前記漏水発生位置推定装置は、
    　前記漏水発生の後の前記各圧力センサの圧力計測値の、前記予測計測値である予測圧力値からの減少値を、前記各圧力センサの前記指標値とし、
    　前記配水管網の、前記指標値が負である、複数の前記圧力センサのうちの第１の圧力センサの下流ではない地域で、前記指標値が、複数の前記センサのうちの他の圧力センサの前記指標値より大である第２の圧力センサの近隣である地域を前記漏水発生の位置として推定することを特徴とする漏水発生位置推定方法。
  12. 　請求項１０に記載の漏水発生位置推定方法であって、
    　前記漏水発生位置推定装置は、
    　さらに前記指標値の誤差上限値と誤差下限値を計算し、
    　複数の前記センサの第１のセンサと第２のセンサについて、前記第１のセンサの前記指標値の前記誤差下限値よりも、前記第２のセンサの前記指標値の前記誤差上限値が小さいときに、前記第１のセンサの前記指標値が前記第２のセンサの前記指標値より大きいとする判定に基づいて前記漏水発生の位置を推定することを特徴とする漏水発生位置推定方法。
  13. 　請求項１０に記載の漏水発生位置推定方法であって、
    　前記漏水発生位置推定装置は、
    　前記漏水イベントデータに含まれる前記漏水発生の時刻以降の複数の時刻において前記指標値を計算し、
    　複数の前記時刻の前記指標値のそれぞれから推定される複数の前記地域の重なる地域を前記漏水発生の位置として推定することを特徴とする漏水発生位置推定方法。
  14. 　請求項１１に記載の漏水発生位置推定方法であって、
    　前記漏水発生位置推定装置は、
    　前記圧力センサ計測値、前記圧力センサの設置標高値、および隣接する前記圧力センサとの間の管網接続関係から、隣接する前記圧力センサおよび前記地域間の上下流関係を推定し、推定した前記上下流関係に基づいて前記漏水発生の位置を推定することを特徴とする漏水発生位置推定方法。
  15. 　請求項１１に記載の漏水発生位置推定方法であって、
    　音響レベルセンサが前記配水管網に配置され、
    　前記漏水発生位置推定装置は、
    　前記漏水発生の後の前記音響レベルセンサの音響レベル計測値の、前記予測計測値である予測音響レベル値からの増加値を、前記音響レベルセンサの前記指標値とし、
    　前記漏水発生の位置として、前記圧力センサの前記指標値に基づいて推定した前記地域の中で、前記音響レベルセンサの前記指標値が最大の音響レベルセンサの近隣である地域を前記漏水発生の位置として推定することを特徴とする漏水発生位置推定方法。
PCT/JP2016/058918 2015-04-28 2016-03-22 漏水発生位置推定装置、システムおよび方法 WO2016174958A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
GB1716409.6A GB2555535B (en) 2015-04-28 2016-03-22 Water leak occurence position estimation device, system, and method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2015-090907 2015-04-28
JP2015090907A JP6471035B2 (ja) 2015-04-28 2015-04-28 漏水発生位置推定装置、システムおよび方法

Publications (1)

Publication Number Publication Date
WO2016174958A1 true WO2016174958A1 (ja) 2016-11-03

Family

ID=57198302

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2016/058918 WO2016174958A1 (ja) 2015-04-28 2016-03-22 漏水発生位置推定装置、システムおよび方法

Country Status (3)

Country Link
JP (1) JP6471035B2 (ja)
GB (1) GB2555535B (ja)
WO (1) WO2016174958A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2572841A (en) * 2018-02-23 2019-10-16 Hitachi Ltd Water pressure gauge arrangement support system and method

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2567180A (en) * 2017-10-05 2019-04-10 Homeserve Plc Leak detection method and apparatus
JP7183463B1 (ja) * 2022-06-28 2022-12-05 株式会社ニュージェック 漏水検知システムおよび漏水検知方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009192329A (ja) * 2008-02-13 2009-08-27 Toshiba Corp 配水管路の漏水診断装置及び漏水診断方法
JP2010048058A (ja) * 2008-08-25 2010-03-04 Toshiba Corp 漏水節点推定装置
JP2014145603A (ja) * 2013-01-28 2014-08-14 Hitachi Ltd 漏水推定装置及びシステム及び方法
JP2014222424A (ja) * 2013-05-14 2014-11-27 株式会社日立製作所 配水制御装置および方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009192329A (ja) * 2008-02-13 2009-08-27 Toshiba Corp 配水管路の漏水診断装置及び漏水診断方法
JP2010048058A (ja) * 2008-08-25 2010-03-04 Toshiba Corp 漏水節点推定装置
JP2014145603A (ja) * 2013-01-28 2014-08-14 Hitachi Ltd 漏水推定装置及びシステム及び方法
JP2014222424A (ja) * 2013-05-14 2014-11-27 株式会社日立製作所 配水制御装置および方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2572841A (en) * 2018-02-23 2019-10-16 Hitachi Ltd Water pressure gauge arrangement support system and method
GB2572841B (en) * 2018-02-23 2020-10-21 Hitachi Ltd Water pressure gauge arrangement support system and method

Also Published As

Publication number Publication date
JP2016205065A (ja) 2016-12-08
GB201716409D0 (en) 2017-11-22
GB2555535B (en) 2020-10-07
JP6471035B2 (ja) 2019-02-13
GB2555535A (en) 2018-05-02

Similar Documents

Publication Publication Date Title
JP6394726B2 (ja) 運用管理装置、運用管理方法、及びプログラム
US10401250B2 (en) Leakage detection and leakage location in supply networks
JP5901140B2 (ja) システムの高い可用性のためにセンサデータを補間する方法、コンピュータプログラム、システム。
JP6304767B2 (ja) センサ監視装置、センサ監視方法、及びセンサ監視プログラム
KR102118818B1 (ko) 밸브 진단 방법 및 밸브 진단 장치
JP6318053B2 (ja) 漏水分布推定装置
JP5164954B2 (ja) 機器診断方法及び機器診断装置
JP6625851B2 (ja) 漏水診断装置、漏水診断方法及びコンピュータプログラム
JP2013073414A (ja) プラントのセンサ診断装置およびセンサ診断方法
EP2613207A2 (en) Adaptive trend-change detection and function fitting system and method
JP6471035B2 (ja) 漏水発生位置推定装置、システムおよび方法
EP3594653B1 (en) Diagnosis cost output device, diagnosis cost output method, and diagnosis cost output
US20220412845A1 (en) Anomaly detection and failure prediction for predictive monitoring of industrial equipment and industrial measurement equipment
JP5253218B2 (ja) 機器診断方法
JP5792667B2 (ja) センサ診断装置及びセンサ診断方法
US20230280240A1 (en) Abnormality diagnosis device and abnormality diagnosis method
JP2017002554A (ja) 管路異常検知装置および方法
JP7079897B2 (ja) 資産としての流量計
JP7095743B2 (ja) 流体漏洩診断装置、流体漏洩診断システム、流体漏洩診断方法、及び、流体漏洩診断プログラム
JP2020129158A (ja) 異常検出装置、シミュレータ、プラント監視システム、異常検出方法及びプログラム
JP2019144044A (ja) ガス器具判別装置
JPWO2019031371A1 (ja) 状態分析装置、状態分析方法及びプログラム
JP6800065B2 (ja) 管理システム、情報処理方法、プログラム、通信装置
WO2015087378A1 (ja) 水道施設網の異常検知装置及び異常検知方法
JP6856124B2 (ja) 分析装置、分析方法及びコンピュータプログラム

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16786235

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 201716409

Country of ref document: GB

Kind code of ref document: A

Free format text: PCT FILING DATE = 20160322

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 16786235

Country of ref document: EP

Kind code of ref document: A1