JP5273531B2 - Inflow position estimation device, inflow position estimation method, program - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To accurately perform the inference of the inflow position of unknown water in a shunt type culvert. <P>SOLUTION: A temporary inflow position of the unknown water is set in the water supply system in the culvert, and the flow rate at the terminal point in the water supply system is calculated after an elapse of a predetermined time since the unknown water flows in, on the assumption that the unknown water flows in at the temporary inflow position. When a difference is greater than a predetermined threshold value, the difference between the calculated flow rate and the actually measured flow rate value in the terminal point at the water supply system after the elapse of the predetermined time since the unknown water flows in, and the temporary inflow position is replaced and set as a new temporary inflow position, on the determination that the processing is executed repeatedly. The flow rate is calculated based on the newly set temporary inflow position, and when the difference is smaller than a predetermined threshold value, the temporary inflow position set at this point is the inferred position, on the determination that the repetitive processing is completed. <P>COPYRIGHT: (C)2010,JPO&amp;INPIT

Description

  The present invention relates to an estimation device, an estimation method, and a program for estimating an inflow position of unknown water in a diversion pipe tube.

  2. Description of the Related Art Conventionally, a separate sewer system has been realized in which sewage such as domestic wastewater and rainwater that has fallen on a road or site are flowed through different pipes. In the sewerage sewer system, the sewage pipes and the rainwater pipes are configured independently. However, if damaged parts occur due to the aging of the pipes, rainwater and groundwater (hereinafter referred to as these) Together, it may be called “unknown water”).

  Such unknown water causes adverse effects such as an increase in operational costs at the sewage treatment plant, a reduction in purification capacity, and a decrease in sewage flow capacity due to lime components adhering to the pipe wall. Therefore, if a phenomenon that indicates the possibility of unknown water flowing in, such as an increase in influent water at a sewage treatment plant, is identified, the inflow location of unknown water is specified to avoid the negative effects of unknown water. It is necessary to take measures such as repair.

  However, conventionally, in order to identify the inflow position of unknown water, it was necessary to conduct a wide area basic survey and a detailed survey, which required a large amount of money, a large amount of human resources, and a large amount of survey time. Therefore, the request | requirement with respect to the technique for simply narrowing down the inflow position of unknown water is increasing.

In response to such a request, the rainfall data showing the time series change of the rainfall in each area to be estimated, and the unknown water quantity data showing the time series change of the unknown water quantity at the base point located downstream from these districts, Based on this, a technique has been proposed in which a correlation value between rainfall data and an unknown water amount in each district is calculated, and this correlation value is output as an unknown water generation distribution in each district (see Patent Document 1).
JP 2005-23763 A

  However, in the conventional technology, the location of unknown water inflow was estimated on a regional basis, so if an area where many pipes are concentrated is obtained as an estimation result, the unknown water flow is actually There has been a problem that much time and cost are required for the processing (detailed investigation) for identifying a specific portion of the pipe tub where the inflow of water has occurred.

  In view of the above circumstances, an object of the present invention is to provide an inflow position estimation device, an inflow position estimation method, and a program capable of estimating the inflow position of unknown water in a shunting pipe with higher accuracy. It is said.

[1] In order to solve the above-described problem, an inflow position estimation device according to one aspect of the present invention includes an assumed inflow position setting unit that sets an assumed inflow position of unknown water in a pipe water supply system, and the assumed inflow position. When it is assumed that unknown water has flowed in, after a predetermined time has passed since the unknown water flowed in, a flow rate calculation unit for calculating a flow rate at the terminal point of the water transmission system, and the unknown water When the difference between the actually measured flow rate measured at the terminal point of the water supply system after the predetermined time has passed and the calculated flow rate is greater than a preset threshold value If the difference between the actual flow rate measurement value and the calculated flow rate is smaller than a preset threshold value , it is determined that the assumed inflow position setting unit and the flow rate calculation unit repeat processing. Position setting And an end determination unit that determines the assumed inflow position set at this time as an estimated position, and the end inflow position setting unit includes the end determination unit. When it is determined that the process is to be repeated, the assumed inflow position is replaced with a new assumed inflow position, and the flow rate calculation unit is configured to newly set the assumed inflow when the end determination unit determines to repeat the process. The flow rate is calculated based on the position.

  [2] In addition, the assumed inflow position setting unit in the inflow position estimation apparatus according to one aspect of the present invention, when repeating the process, if the actual flow rate value is larger than the calculated flow rate, If a point on the downstream side of the assumed inflow position set in the process is replaced with a new assumed inflow position, and the measured flow rate is smaller than the calculated flow rate, the immediately preceding process is performed. The point upstream of the assumed inflow position set in step 1 may be set as a new assumed inflow position.

  [3] Further, the assumed inflow position setting unit in the inflow position estimation apparatus according to one aspect of the present invention firstly sets the start end of the main pipe in the water transmission system as the start point, the end point of the main pipe as the end point, and the start point and the end point Is set as the assumed inflow position, and when the process is repeated, if the measured flow rate is larger than the calculated flow rate, the assumed inflow position set in the immediately preceding process is set. If the intermediate point between the start point and the end point is replaced with the new assumed inflow position without changing the start point and the end point, the measured flow rate is smaller than the calculated flow rate. The assumed inflow position set in the process is set as the end point, the start point is left as it is, and an intermediate point between the start point and the end point is replaced and set as a new assumed inflow position. Good.

  [4] In addition, the assumed inflow position setting unit in the inflow position estimation apparatus according to an aspect of the present invention, when repeating the process, starts from the start end of the other pipe having the merge position on the downstream side, A confluence position may be set as the end point, and an intermediate point between the start point and the end point may be set as the new assumed inflow position.

  One aspect of the present invention may be specified as an inflow position estimation method performed by the above-described inflow position estimation apparatus. One embodiment of the present invention may be specified as a computer program for causing a computer to execute the inflow position estimation method.

  According to the present invention, it is possible to estimate the inflow position of unknown water in a diverter type pipe rod with higher accuracy.

  FIG. 1 is a functional block diagram illustrating a functional configuration of the inflow position estimation apparatus 1. As shown in the figure, the inflow position estimation device 1 includes a water supply system data storage unit 11, an assumed inflow position setting unit 12, an initial value storage unit 13, a flow rate calculation unit 14, and an end determination unit 15.

  The water supply system data storage unit 11 stores data (water supply system data) related to a pipe network (water supply system) that is an estimation target of the unknown water inflow position. The water supply system data includes, for example, information on the length of the pipe rod, the pipe diameter of the pipe rod, the gradient of the pipe rod, and the joining position of the main pipe and the sub pipe.

  The assumed inflow position setting unit 12 sets a position (assumed inflow position) on which the unknown water is assumed to flow on the water supply system based on the water supply system data. In addition, the assumed inflow position setting unit 12 ends the process in the entire inflow position estimating apparatus 1 when there is no pipe to set a new assumed inflow position.

The initial value storage unit 13 stores various initial values used in the calculation in the flow rate calculation unit 14.
The flow rate calculation unit 14 is a water supply system after a predetermined time has elapsed since the unknown water flowed in, assuming that unknown water is flowing in at the assumed flow-in position set by the assumed flow-in position setting unit 12. The flow rate at the end point of is calculated. The end of the water supply system is, for example, a pump station or a sewage treatment plant.

  The end determination unit 15 determines whether or not to end the process for the pipe rod currently being processed, based on the calculation result in the flow rate calculation unit 14 and the actual flow rate value input from the outside. The actual flow rate measurement value is an actual measurement value of the flow rate at the terminal point of the water supply system after a predetermined time has elapsed since the unknown water flowed in. When the end determination unit 15 determines not to end the process, the assumed inflow position setting unit 12 sets a new assumed inflow position for the pipe so that the flow rate calculation unit 14 determines the flow rate based on the assumed inflow position. calculate. On the other hand, if the end determination unit 15 determines to end the process, the end determination unit 15 adds the current assumed inflow position to the inflow position estimation result.

  Next, the process of the flow rate calculation unit 14 will be described. When the inside of the pipe tub is an open channel not filled with water, assuming that the flow of water in the pipe tub is a one-dimensional flow having a constant density, the fluid mass conservation equation is expressed as shown in Equation 1.

  In Equation 1, x is the distance from the starting end of the pipe rod to be treated, A is the cross-sectional area of the flowing water, V is the average flow velocity, and q is the inflow rate per unit volume from the outside into the pipe rod. Specifically, the unit volume is calculated by subtracting the inflow rate at the end position when no unknown water inflow occurs (for example, in fine weather) from the inflow rate at the end position when inflow of unknown water (for example, in rainy weather). The value of q is obtained by converting it into the per-flow rate. In addition, since the inflow flow rate when the inflow of unknown water does not occur has different patterns depending on the day, this value may be generalized. For example, it may be generalized for each day of the week by calculating an average value using a plurality of inflows for the same day of the week. Similarly, it may be generalized for each month or season by calculating an average value.

  Next, assuming that the flow of water in the pipe rod is a one-dimensional flow having a constant density, the fluid momentum conservation equation is expressed as shown in Equation 2.

In Equation 2, g represents gravitational acceleration, f ′ represents a frictional resistance coefficient, h represents the height from the bottom of the pipe rod to the water surface, and R represents the diameter depth. The diameter depth is a value obtained by dividing the cross-sectional area of flowing water by the length of the bottom surface in contact with water, which is the wet side length. By converting the partial differential equations of Equations 1 and 2 into ordinary differential equations by the characteristic curve method, the basic equations in the case of C ⁺ (Equations 3 and 4) and the basic equations in the case of C ⁻ (Equation 5 And Equation 6) can be obtained.

In Equations 3 to 6, y is the water depth, S is the gradient, c is the velocity of propagation of the gravitational wave, and _Ay is a partial differential of the flowing water cross-sectional area with the water depth. By discretizing Equations 3 to 6 of the basic equations by the difference method, Equations 7 to 10 can be obtained.

V _i and y _i in Equation 7 represent the average flow velocity V and the water depth y at point i after Δt seconds, respectively. At this time, each value of V _Ri , V _Si , y _Ri , and y _Si in Equations 7 to 10 can be calculated by a specified time interval method. FIG. 2 is a conceptual diagram showing the concept of the specified time interval method. Each value mentioned above can be represented like the following formula | equation 11-Formula 14 by the regular time interval method.

In Equations 11 to 14, each value of V _i , V _i , y _i , and y _i at t0 is required. These values are stored in advance in the initial value storage unit 13 as initial values. Yes. These values are values calculated on the basis of the flow rate data of pipes measured in a natural flow state where unknown water does not flow. Each value is calculated in advance by solving the following equations 15 to 19 for V and y, and is stored in the initial value storage unit 13. Also, the initial value storage unit 13 stores other parameter values in advance.

  In Equations 15 to 19, Q is the flow rate, A is the flowing water cross-sectional area, V is the flow velocity, y is the water depth, d is the pipe radius, n is the Manning's roughness coefficient, R is the depth, and the flowing water cross-sectional area is wet. A value divided by the length (the length of water in contact with the wall surface in the cross section), I represents the hydrodynamic gradient. Moreover, ψ represents an angle between a line segment extending from the center of the circle to the contact point between the water surface and the wall surface of the pipe and a line segment passing through the center of the circle and perpendicular to the water surface. FIG. 3 is a cross-sectional view of the tube, showing the values of d, ψ, A, and y. FIG. 3A shows the case where y is d or less, and FIG. 3B shows the case where y is larger than d.

  By repeating the processing using Expression 7 to Expression 14, the flow rate at each point when the inflow time has elapsed from the time point t0 can be calculated. Therefore, the flow rate calculation unit 14 calculates the flow rate at the end point when the inflow time has elapsed from the time point t0 by repeatedly executing the processes using Equations 7 to 14.

  FIG. 4 is a schematic diagram showing an outline of the water supply system. In FIG. 4, the pump station is the end, the band extending leftward from the pump station represents the main pipe, and the scale displayed below the band representing the main pipe represents the distance from the start end (left end) of the main pipe. Further, a band extending upward from a point 1000 m from the start end of the main pipe represents the sub pipe, and a scale displayed on the left of the band representing the sub pipe represents a distance from the start end (upper end) of the sub pipe. In addition, a band extending leftward from a point 170 m from the start end of the secondary pipe represents a secondary pipe (hereinafter referred to as “secondary secondary pipe”) when the secondary pipe is used as a main pipe. The displayed scale represents the distance from the start end (left end) of the sub-sub pipe. The water flowing through the secondary pipe flows from the start to the secondary pipe, the water flowing through the secondary pipe flows from the start to the main pipe, and the water flowing through the main pipe flows from the start to the pumping station (terminal). ).

FIG. 5 is a flowchart illustrating an operation example of the inflow position estimation apparatus 1. Hereinafter, the operation of the inflow position estimation apparatus 1 will be described with reference to FIGS. 4 and 5.
First, an actual flow rate value is input to the inflow position estimation device 1 (step S101). This input is performed by an input unit (not shown), and may be input by a user operating an input device such as a keyboard, or may be input by being transmitted from another information processing device via a network. Alternatively, the value stored in the recording medium may be read and input.

  Next, the assumed inflow position setting unit 12 selects a main pipe as a pipe to be processed (step S102). Then, the assumed inflow position setting unit 12 sets an assumed inflow position for the main pipe using a bisection algorithm (step S103). Specifically, the assumed inflow position setting unit 12 first assumes the start point of the pipe to be processed (currently the main pipe) as the start point, the end point of the pipe to be processed as the end point, and assumes an intermediate point between the start point and the end point Set as inflow position. In the case of FIG. 4, a point of 750 m in the main pipe is first set as the assumed inflow position.

  The flow rate calculation unit 104 calculates the flow rate (hereinafter referred to as “flow rate estimated value”) at the end position of the main pipe after the inflow time has elapsed from the inflow start time based on the currently set assumed inflow position (step S104). . The inflow start time and the inflow time will be described later.

  Next, the end determination unit 105 subtracts the actual flow rate value from the estimated flow rate value calculated by the flow rate calculation unit 104 to calculate an estimated actual measurement difference (step S105). If the estimated actual difference is greater than or equal to a preset threshold value (NO in step S106), the end determination unit 105 determines not to end the process for the tube currently being processed.

  In this case, the assumed inflow position setting unit 12 newly sets the assumed inflow position. The assumed inflow position setting unit 12 determines whether an assumed inflow position that is newer than the assumed inflow position set immediately before should be set on the upstream side or the downstream side, based on the previous calculation result in the flow rate calculation unit 14. To do. Specifically, when the previous flow rate estimated value is larger than the actual flow rate measured value, that is, when the estimated actual difference is 0 or more (step S107—NO), the assumed inflow position setting unit 12 sets the assumed inflow position set immediately before. It is determined that a new assumed inflow position is set on the upstream side. Then, the assumed inflow position setting unit 12 sets the end point as the immediately preceding assumed inflow position and the intermediate point between the start point and the end point as the new assumed inflow position, with the start point being the same as immediately before (step S108). In the case of FIG. 4, when the immediately preceding assumed inflow position is 750 m, the end point is a point of 750 m, and the new assumed inflow position is a point of 375 m.

  On the other hand, when the previous flow rate estimated value is smaller than the actual flow rate measured value, that is, when the estimated actual measured difference is less than 0 (YES in step S107), the assumed inflow position setting unit 12 is set to be less than the assumed inflow position set immediately before. It is determined that a new assumed inflow position is set on the downstream side. Then, the assumed inflow position setting unit 12 sets the start point as the immediately preceding assumed inflow position and the intermediate point between the start point and the end point as the new assumed inflow position while the end point remains the same as immediately before (step S109). In the case of FIG. 4, when the immediately preceding assumed inflow position is 750 m, the start point is a 750 m point, and the new assumed inflow position is a 1125 m point.

Then, based on the assumed inflow position newly set in step S108 or step S109, the processes after step S104 are repeatedly executed.
In the process of step S106, when the estimated actual difference is less than the threshold value (step S106-YES), the end determination unit 15 determines to end the process for the pipe being currently processed, and estimates the current assumed inflow position. It adds to a result (step S110).

  Next, the assumed inflow position setting unit 12 joins the sub pipe in the pipe that was being processed immediately before (for example, if the pipe that was being processed just before is a sub pipe, In step S111, it is determined whether or not there is a secondary pipe that is downstream of the assumed inflow position set immediately before. When such a secondary pipe exists, the assumed inflow position setting unit 12 selects the secondary pipe as a pipe to be processed next (step S102), and sets an assumed inflow position for the secondary pipe (step S102). S103). At this time, the assumed inflow position setting unit 12 uses the start end of the sub-pipe as the start point, the end point (confluence point) of the sub-pipe as the end point, and the start point and the end point similarly to the setting process performed for the main pipe first. Is set as the assumed inflow position. And the inflow position estimation apparatus 1 repeatedly performs the process after step S104.

  On the other hand, if there is no secondary pipe in the downstream side of the assumed inflow position where the merging position is set immediately before among the secondary pipes in the pipe rod that was being processed immediately before (step S111-YES), the assumed inflow The position setting unit 12 does not set the assumed inflow position, and the end determination unit 15 outputs the current estimation result (step S112).

  In the case of FIG. 4, it is assumed that the pipe tube that was being processed immediately before is the main pipe and the assumed inflow position set immediately before is downstream of 1000 m (when it is between 1000 m and 1500 m). The inflow position setting unit 12 does not perform processing on the sub pipe and the sub sub pipe. On the other hand, when the pipe rod that was being processed immediately before is the main pipe and the assumed inflow position set immediately before is upstream of 1000 m (when it is between 0 m and 1000 m), the assumed inflow position setting The unit 12 performs processing for the secondary pipe. Furthermore, when the pipe tube being processed immediately before is a secondary pipe and the assumed inflow position set immediately before is upstream from 170 m (when it is between 0 m and 170 m), the assumed inflow position The setting unit 12 performs processing on the sub-sub pipe, and when it is on the downstream side (when it is between 170 m to 750 m), the assumed inflow position setting unit 12 does not perform processing on the sub-sub pipe.

  FIG. 6 is a graph showing the calculation result of the flow rate calculation unit 14. The measured flow rate in FIG. 6 is the flow rate at the end point after a lapse of a predetermined time (hereinafter referred to as “inflow time”) from the time when a large amount of unknown water began to flow (hereinafter referred to as “inflow start time”). It is a measured value. The point in time when a large amount of unknown water starts to flow can be set as the point in time when the rainfall in the area where the pipe to be processed exceeds the threshold, based on, for example, rainfall data indicating the time series of rainfall. The length of the inflow time is longer than the time from when a large amount of unknown water starts to flow until the time when the unknown water reaches the end and the flow rate at the end begins to increase, and the flow rate at the end is constant. It is preset as a time shorter than the time up to the time. This inflow time is appropriately determined by the user based on the time-series data of the flow rate actual measurement values. For example, when a rainfall amount equal to or greater than the threshold value is observed at time T1, the actual flow rate value at the end starts to increase at time T2, and the actual flow rate value at the end becomes constant at time T3, the time T1 is the inflow start point. The value between (T2-T1) and (T3-T1) is the inflow time. The inflow time is set to, for example, 100 seconds, 400 seconds, or 10,000 seconds. In addition, the inflow of unknown water is estimated based on the rainfall data at this time, set in advance as an initial value, and stored in advance in the initial value storage unit 13.

  The graph of FIG. 6 is a graph showing the calculation result (estimated amount) of the flow rate at each point of the main pipe after the inflow time has elapsed since the start of the inflow when the inflow position is set in the secondary pipe. The axis represents the distance from the start of the main pipe, the vertical axis represents the flow rate calculation result at each point, and the rightmost end of the horizontal axis represents the end of the main pipe. In the case of FIG. 6, since the inflow position is set in the secondary pipe, the flow rate rapidly increases at a point of 1000 m where the secondary pipe joins the main pipe. Then, the flow rate decreases from the junction point to the end point. Therefore, the closer the inflow position of unknown water is to the end, the greater the end flow rate at the same point. Accordingly, when the flow rate calculation result is larger than the actual flow rate measurement value, the assumed inflow position set at that time is set closer to the end (that is, on the downstream side) than the actual inflow position. Conversely, when the flow rate calculation result is smaller than the actual flow rate measurement value, the assumed inflow position set at that time is set farther from the end (that is, upstream) than the actual inflow position. Therefore, the assumed inflow position setting unit 12 operates as described above, so that the assumed inflow position is gradually set closer to the actual inflow position. And it becomes possible to obtain a more accurate estimation result by estimating the assumption inflow position when the difference becomes less than a threshold as an inflow position of unknown water.

<Modification>
In the inflow position estimation apparatus 1 described above, since the inflow position is estimated for each pipe as the inflow position estimation result, a plurality of inflow positions may be output. For example, in the water transmission system shown in FIG. 2, there are cases where two points are obtained as an estimation result, from the junction point of the main pipe and the sub pipe to the 600 m upstream point of the main pipe and from the joint point to the 400 m point upstream of the sub pipe. is there. In this case, whether the inflow of unknown water has occurred in the vicinity of both of the two locations, even though the actual inflow of unknown water is one of these two locations. It becomes necessary to do.

  FIG. 7 is a functional block diagram illustrating a modification of the functional configuration of the inflow position estimation apparatus 1. In response to the above-described problem, the inflow position estimation device 1 may be configured to further include an inflow position estimation unit 16 as shown in FIG.

  The inflow position estimation unit 16 estimates one inflow position from among a plurality of estimation results. Specifically, the inflow position estimation unit 16 changes the flow rate at the terminal end point of the main pipe until the stable time elapses from the inflow start time (that is, the flow rate at each time point) for each inflow position included in the estimation result. Is calculated. The stabilization time is a time longer than the inflow time, and is a time required for the flow rate at the end point of the main pipe to be constant from the start of the inflow. The value of the stabilization time is appropriately determined by the user based on the time series data of the flow rate actual measurement values. For example, when a rainfall amount equal to or greater than the threshold value is observed at time T1, the actual flow rate value at the end starts to increase at time T2, and the actual flow rate value at the end becomes constant at time T3, from (T3-T1) A larger value is the stabilization time. The stabilization time is set to 200 seconds, 600 seconds, or 15000 seconds, for example.

  FIG. 8 is a graph showing the time change of the flow rate at the end point of the main pipe. The horizontal axis in FIG. 8 represents time, the left end represents the inflow start time, and the right end represents the elapsed time of the stabilization time. Moreover, the vertical axis | shaft of FIG. 8 represents the flow volume in the terminal point of the main pipe | tube in each time. The solid line (symbol a) in FIG. 8 is a graph in the case of an unknown water inflow position at a point 600 m upstream from the junction of the main pipe and the sub pipe in the water transmission system shown in FIG. Moreover, the graph of the broken line (code | symbol b) of FIG. 8 is a graph in case the inflow position of unknown water exists in the 400-m upstream point of a subpipe from the same confluence | merging point. In FIG. 8, the broken line graph appears only partially (between time t1 and time t2), but the broken line graph overlaps the solid line graph between t0 and t1 and between t2 and t3. Exist.

  Signs A and B have the same change in flow rate after time t2. Therefore, when the flow rate calculation part 14 is calculating the flow volume after t2, as above-mentioned, these two points are obtained as an estimation result. However, as is clear from FIG. 8, the graph of the symbol A and the graph of the symbol B have different slopes in a transient state (particularly, a state between time t1 and time t2). Accordingly, the inflow position estimation unit 16 selects and outputs an assumed inflow position having a value closest to the inclination in the transient state of the actually measured value among the inclinations in the transient state of the assumed inflow position of each estimation result. It is possible to narrow down the estimation result to one.

  More specifically, the flow rate position estimation unit 16 receives an input of an actual measurement value of a temporal change in flow rate at the end point of the main pipe. Further, the flow rate position estimation unit 16 calculates the temporal change in the flow rate at the terminal end point of the main pipe for all the assumed inflow positions included in the estimation result after the process of step S111. Then, the inflow position estimation unit 16 calculates the gradient in the transient state for all the assumed inflow positions included in the estimation results and the input actual measurement values, and assumes the inflow position corresponding to the inclination closest to the inclination in the actual measurement values. Is selected as the estimation result. In this case, the inflow position estimation device 1 may output only the assumed inflow position selected by the flow rate position estimation unit 16 or outputs the assumed inflow position selected by the flow rate position estimation unit 16 as the most likely candidate. Further, other estimation results may be output as other candidates.

You may make it implement | achieve the function of the inflow position estimation apparatus 1 in embodiment mentioned above with a computer. In that case, a program for realizing this function may be recorded on a computer-readable recording medium, and the program recorded on the recording medium may be read into a computer system and executed. Here, the “computer system” includes an OS and hardware such as peripheral devices. The “computer-readable recording medium” refers to a storage device such as a flexible medium, a magneto-optical disk, a portable medium such as a ROM and a CD-ROM, and a hard disk incorporated in a computer system. Furthermore, the “computer-readable recording medium” dynamically holds a program for a short time like a communication line when transmitting a program via a network such as the Internet or a communication line such as a telephone line. In this case, a volatile memory inside a computer system serving as a server or a client in that case may be included and a program that holds a program for a certain time. The program may be a program for realizing a part of the functions described above, and may be a program capable of realizing the functions described above in combination with a program already recorded in a computer system.
The embodiment of the present invention has been described in detail with reference to the drawings. However, the specific configuration is not limited to this embodiment, and includes designs and the like that do not depart from the gist of the present invention.

It is a functional block diagram showing the functional structure of an inflow position estimation apparatus. It is a conceptual diagram showing the concept of a regular time interval method. It is sectional drawing of a pipe | tube, and is a figure showing the value of d, (psi), A, and y. It is the schematic showing the outline of a water supply system. It is a flowchart showing the operation example of an inflow position estimation apparatus. It is a graph showing the calculation result of a flow volume calculation part. It is a functional block diagram showing the modification of the function structure of an inflow position estimation apparatus. It is a graph showing the time change of the flow volume in the terminal point of a main pipe.

Explanation of symbols

DESCRIPTION OF SYMBOLS 1 ... Inflow position estimation apparatus, 11 ... Water supply system data storage part, 12 ... Assumed inflow position setting part, 13 ... Initial value storage part, 14 ... Flow rate calculation part, 15 ... End determination part, 16 ... Inflow position estimation part

Claims (6)

An assumed inflow position setting unit for setting an assumed inflow position of unknown water in the pipe water supply system;
When it is assumed that unknown water is flowing in at the assumed inflow position, a flow rate calculation unit that calculates a flow rate at a terminal point of the water supply system after a predetermined time has passed since the unknown water flowed in; ,
The difference between the actual measured flow rate at the end point of the water transmission system after the lapse of the predetermined time from the inflow of the unknown water and the calculated flow rate is a preset threshold value. If it is greater than the threshold value , it is determined that the assumed inflow position setting unit and the flow rate calculation unit repeat the process, and the difference between the actual flow rate value and the calculated flow rate is smaller than a preset threshold value. the assumption inflow position setting unit and determines to end the iterative process of the flow rate calculation unit includes a termination determination unit to estimate position assumptions inflow position set at this time, the,
The assumed inflow position setting unit is configured to replace the assumed inflow position with a new assumed inflow position when the end determination unit determines to repeat the process,
The flow rate calculation unit is an inflow position estimation device that calculates the flow rate based on a newly set assumed inflow position when the end determination unit determines that the process is to be repeated. The assumed inflow position setting unit, when repeating the process, if the measured flow rate is larger than the calculated flow rate, a point downstream of the assumed inflow position set in the immediately preceding process. Is replaced with a new assumed inflow position, and if the measured flow rate is smaller than the calculated flow rate, a point upstream of the assumed inflow position set in the immediately preceding process is newly set. Set as a hypothetical inflow position,
The inflow position estimation apparatus according to claim 1.   The assumed inflow position setting unit first sets a start point of the main pipe in the water supply system as a start point, an end point of the main pipe as an end point, sets an intermediate point between the start point and the end point as the assumed inflow position, and repeats the process. In the case where the actual flow rate value is larger than the calculated flow rate, the assumed inflow position set in the immediately preceding process is set as the start point, the end point is left as it is and the intermediate point between the start point and the end point. When the point is replaced with a new assumed inflow position and the measured flow rate is smaller than the calculated flow rate, the assumed inflow position set in the immediately preceding process is set as the end point and the start point. The inflow position estimation apparatus according to claim 2, wherein the intermediate point between the start point and the end point is set as a new assumed inflow position without changing. When the process is repeated, the assumed inflow position setting unit sets the start point of the other pipe having the merge position on the downstream side as the start point, the merge position as the end point, and an intermediate point between the start point and the end point. Set as the new assumed inflow position,
The inflow position estimation apparatus according to claim 3. An assumed inflow position setting step in which the inflow position estimating device sets an assumed inflow position of unknown water in the water supply system of the pipe;
When the inflow position estimation device assumes that unknown water is flowing in at the assumed inflow position, the flow rate at the end point of the water supply system after a predetermined time has passed since the unknown water flowed in A flow rate calculating step for calculating
The difference between the flow rate actually measured at the terminal point of the water supply system after the predetermined time has elapsed after the unknown water has flowed in and the calculated flow rate. Is larger than a preset threshold value, it is determined that the processes of the assumed inflow position setting step and the flow rate calculation step are repeatedly executed, and a difference between the actually measured flow rate value and the calculated flow rate is preset. If the less than the threshold value, the assumption inflow position setting step and determined to end the iterative process of the flow rate calculating step, organic and termination determination step for the estimated position assumptions inflow position set at this point, the And
When the inflow position estimation device determines to repeat the process in the end determination step, the assumed inflow position setting step replaces and sets the assumed inflow position as a new assumed inflow position,
The inflow position estimation device calculates the flow rate based on a newly set assumed inflow position in the flow rate calculation step when determining repetition of processing in the end determination step.
Inflow position estimation method. An assumed inflow position setting step for setting an assumed inflow position of unknown water in the pipe water supply system;
A flow rate calculating step for calculating a flow rate at a terminal point of the water transmission system after a predetermined time has elapsed since the unknown water flowed in, assuming that unknown water is flowing in at the assumed inflow position; ,
The difference between the actual measured flow rate at the end point of the water transmission system after the lapse of the predetermined time from the inflow of the unknown water and the calculated flow rate is a preset threshold value. In the case where the difference between the actual flow rate measurement value and the calculated flow rate is smaller than a preset threshold value. , A program for causing a computer to execute an end determination step in which it is determined that the repeated process of the assumed inflow position setting step and the flow rate calculation step is to be ended and the assumed inflow position set at this time is an estimated position. Yes,
If it is determined that the process is repeated in the end determination step, the assumed inflow position setting step is performed by replacing the assumed inflow position with a new assumed inflow position,
When it is determined that the process is repeated in the end determination step, the flow rate is calculated based on a newly set assumed inflow position in the flow rate calculation step.
program.

Images