Hydraulic analysis method and hydraulic analysis apparatus

Abstract

This invention provides a hydraulic analysis method that can accurately assess the condition of a pipeline network. [Solution] The method comprises an analysis condition setting step, an analysis execution step, and an analysis result evaluation step. The analysis condition setting step includes an inflow / outflow intersection definition step that defines the inflow and outflow intersections of water to the pipeline network, and a daily average inflow (Q) calculated from the total inflow (Q) obtained by adding up the daily unit time inflow (Qt) measured at the inflow intersections and corresponding to the day to be analyzed. AVR A hydraulic analysis method comprising: a daily average water extraction amount setting step of setting the daily average water extraction amount (W1) at the extraction intersection; a grouping step of dividing the extraction intersections into groups that share a common water extraction trend; and a group-specific daily average water extraction amount setting step of setting the value obtained by allocating the daily average water extraction amount (W1) to each group based on the ratio of daily average water extraction amounts (S1) for each group obtained in advance, as the group-specific daily average water extraction amount (W2).

Description

[Technical Field]

[0001] The present invention relates to a hydraulic analysis method and hydraulic analysis apparatus for a pipeline network composed of multiple pipelines and multiple intersections where each pipeline intersects. [Background technology]

[0002] In recent years, IoT technology has attracted attention for its ability to remotely monitor hydraulic and water quality information (water pressure, flow rate, flow direction, residual chlorine concentration, etc.) in pipeline networks by placing water pressure sensors and flow velocity sensors at multiple locations within the network.

[0003] However, due to cost constraints and limitations on the placement of sensors, there are limits to expanding the number of measurement points. Therefore, it is desirable to comprehensively understand the state of the entire pipeline network by utilizing multi-point measurement data obtained through IoT technology and the results of hydraulic analysis of the pipeline network.

[0004] It is desirable to use computer-based hydraulic analysis to understand the condition of the pipeline network that constitutes the water supply system, to respond quickly to abnormalities, and to formulate an efficient pipeline renewal plan.

[0005] For example, Patent Document 1 proposes a pipeline network analysis method and a pipeline network analysis apparatus that perform hydraulic analysis on a water distribution pipeline network in which a water reservoir and individual water demand points are connected by multiple pipelines to calculate pipeline network characteristic values ​​at any given water demand point. [Prior art documents] [Patent Documents]

[0006] [Patent Document 1] Japanese Patent Publication No. 2022-091533 [Overview of the project] [Problems that the invention aims to solve]

[0007] The conventional hydraulic analysis techniques described above are methods for calculating water pressure and flow rate under certain hydraulic conditions, and have been used for evaluation in the planning and design of pipelines. The usual method of use involves using the average water volume based on actual water usage observed from metered water readings to determine the water volume taken from any extraction point in the pipeline network. This average volume is then multiplied by a time coefficient (the ratio of average water volume to the maximum water volume per hour) to calculate the water volume at maximum dynamic pressure, and confirming that there are no areas with low water pressure or high flow velocity.

[0008] However, in order to apply hydraulic analysis techniques to pipeline networks for the maintenance and management of equipment, it is necessary to accurately calculate the water pressure and flow rate at each intersection where water is drawn under certain hydraulic conditions, such as a desired day or time of day. To do this, it is necessary to calculate the amount of water drawn at each intersection that reflects actual usage.

[0009] The objective of the present invention is to provide a hydraulic analysis method and hydraulic analysis apparatus that can accurately grasp the state of a pipeline network. [Means for solving the problem]

[0010] To achieve the above objective, the first characteristic configuration of the hydraulic analysis method according to the present invention is a hydraulic analysis method for a pipeline network composed of multiple pipelines and multiple intersections where each pipeline intersects, comprising: an analysis condition setting step of setting analysis conditions for the pipeline network; an analysis execution step of performing a hydraulic analysis based on the analysis conditions and outputting analysis results; and an analysis result evaluation step of comparing and evaluating the analysis results with measured values, wherein the analysis condition setting step comprises an inflow / outflow intersection definition step of defining inflow intersections where water flows into the pipeline network and outflow intersections where water is taken out; and a daily average inflow (Q) obtained from the total inflow (Q) which is measured at the inflow intersections and obtained by adding up the daily inflow (Qt) corresponding to the day to be analyzed. AVRThe method includes: setting the average daily water volume (W1) at the extraction intersection; a grouping step that divides the extraction intersection into groups with common water extraction trends; and setting the average daily water volume (W2) by group by distributing the average daily water volume (W1) to each group based on the ratio (S1) of the average daily water volume (S1) for each group obtained in advance.

[0011] In the analysis condition setting step, the time-averaged inflow volume measured at the inflow intersection is set as the daily average outflow volume at the outflow intersection. Additionally, the outflow intersections are divided into groups with common water outflow trends, and the daily average outflow volume (W1) is allocated as the group-specific daily average outflow volume (W2) based on the ratio of the daily average outflow volumes for each group, thereby enabling analysis that reflects the actual outflow volume.

[0012] The second characteristic configuration is that, in addition to the first characteristic configuration described above, the analysis condition setting step includes a group-specific water volume correction coefficient setting step which obtains a group-specific water volume correction coefficient (S2) for each unit time from the time variation pattern of the water volume extracted by each group that is similar to the analysis target day, among the daily time variation patterns of the water volume extracted by each group that have been measured in advance for each analysis condition including season, day of the week, and event; and a group-specific water volume calculation step which calculates the group-specific water volume extracted by each unit time (W3) from the group-specific daily average water volume extracted (W2) and the group-specific water volume correction coefficient (S).

[0013] From the daily time variation patterns of the water volume extracted by each group, which have been measured in advance, a group-specific water volume extraction correction coefficient (S) for each unit time that matches the season, day of the week, event, etc. to be analyzed can be obtained. Then, from the group-specific daily average water volume extracted (W2) and the group-specific water volume extraction correction coefficient (S), a group-specific water volume extracted per unit time that is closer to the actual situation (W3) can be obtained.

[0014] The third characteristic configuration, in addition to the second characteristic configuration described above, includes a water extraction amount correction coefficient calculation step for calculating a unit time water extraction amount correction coefficient (m) that corrects the total value of the group-by-group water extraction amounts (W3) per unit time so as to match the unit time inflow water amount.

[0015] The total value of the group-by-group water extraction amounts (W3) per unit time and the unit time inflow water amount can be made to match by the unit time water extraction amount correction coefficient (m).

[0016] The fourth characteristic configuration, in addition to the third characteristic configuration described above, is that the analysis condition setting step includes a water extraction amount correction coefficient calculation step for setting the product of the group-by-group water extraction amount correction coefficient (S) and the unit time water extraction amount correction coefficient (m) as the water extraction amount correction coefficient (K).

[0017] The fifth characteristic configuration, in addition to the fourth characteristic configuration described above, is that the analysis execution step repeatedly executes a process of outputting the water pressure, flow rate, or flow direction, which is the analysis value at the extraction intersection, by performing the hydraulic analysis under the analysis conditions in which the water extraction amount correction coefficient (K) is adjusted.

[0018] The sixth characteristic configuration, in addition to the fifth characteristic configuration described above, is that the analysis result evaluation step includes a deviation degree evaluation step for evaluating the deviation degree of the analysis value obtained by the analysis execution step with respect to the measured value, and an analysis condition determination step for determining the optimal analysis condition from among the plurality of analysis conditions based on the evaluation result of the deviation degree.

[0019] By repeating the execution of the hydraulic analysis based on different values of the group correction coefficient and determining the optimal analysis condition with the smallest deviation degree with respect to the measured value of the water pressure, flow rate, or flow direction, an accurate analysis result for the entire pipeline network can be obtained.

[0020] The first characteristic configuration of the hydraulic analysis device according to the present invention is a hydraulic analysis device for a pipeline network composed of multiple pipelines and multiple intersections where each pipeline intersects, comprising: an analysis condition setting unit for setting analysis conditions for the pipeline network; an analysis execution unit for performing hydraulic analysis based on the analysis conditions and outputting analysis results; and an analysis result evaluation unit for comparing and evaluating the analysis results with measured values, wherein the analysis condition setting unit comprises an inflow / outflow intersection definition processing unit for defining inflow intersections where water flows into the pipeline network and outflow intersections where water is taken out, and a daily average inflow (Q) obtained from the total inflow (Q) which is measured at the inflow intersections and obtained by adding up the daily inflow (Qt) corresponding to the day to be analyzed. AVR The system includes: a daily average water extraction volume setting processing unit that sets the daily average water extraction volume (W1) at the extraction intersection as the daily average water extraction volume (W1); a grouping processing unit that divides the extraction intersection into groups with common water extraction trends; and a group-specific daily average water extraction volume setting processing unit that sets the value obtained by allocating the daily average water extraction volume (W1) to each group as the group-specific daily average water extraction volume (W2) based on the ratio (S1) of the daily average water extraction volume for each group that has been acquired in advance.

[0021] The second characteristic configuration is that, in addition to the first characteristic configuration described above, the analysis condition setting unit includes a group-specific water volume correction coefficient acquisition processing unit that acquires a group-specific water volume correction coefficient (S2) for each unit of time from the time variation pattern of the water volume extracted by each group that matches the day of analysis, among the time variation patterns of the water volume extracted by each group that have been measured in advance for each analysis condition including season, day of the week, and event; and a group-specific water volume calculation processing unit that calculates the group-specific water volume extracted by each unit of time (W3) from the group-specific daily average water volume extracted (W2) and the group-specific water volume correction coefficient (S).

[0022] The third characteristic configuration is that, in addition to the second characteristic configuration described above, the analysis condition setting unit includes a unit time water extraction correction processing unit that calculates and corrects a unit time water extraction correction coefficient (m) that corrects the sum of the group-specific water extraction amounts (W3) per unit time to match the unit time inflow amount (Qt).

[0023] The fourth characteristic configuration is that, in addition to the third characteristic configuration described above, the analysis condition setting unit includes a water volume correction coefficient calculation processing unit that sets the water volume correction coefficient (K) as the product of the group-specific water volume correction coefficient (S) and the water volume correction coefficient (m) per unit time.

[0024] The fifth characteristic configuration is that, in addition to the fourth characteristic configuration described above, the analysis execution processing unit performs the hydraulic analysis multiple times by adjusting the extraction water volume correction coefficient (K) under analysis conditions, thereby outputting the analysis values ​​of water pressure, flow rate, or flow direction at the extraction intersection.

[0025] The fifth characteristic configuration is that, in addition to the fourth characteristic configuration described above, the analysis result evaluation unit includes a deviation degree evaluation processing unit that evaluates the degree of deviation of the analysis value obtained by the analysis execution processing unit from the measured value, and an analysis condition determination processing unit that determines the optimal analysis condition from among the plurality of analysis conditions based on the deviation degree evaluation result. [Effects of the Invention]

[0026] As described above, the present invention provides a hydraulic analysis method and hydraulic analysis apparatus that can accurately grasp the state of a pipeline network. [Brief explanation of the drawing]

[0027] [Figure 1] Functional block diagram of the hydraulic analysis device. [Figure 2] (a) is a diagram illustrating the pipeline network, (b) is a diagram illustrating the flow equation, and (c) is a diagram illustrating the closed pipeline equation. [Figure 3] A flowchart showing the overall process of hydraulic analysis. [Figure 4] A flowchart illustrating the processes for the analysis execution step and the analysis result evaluation step. [Figure 5] Flowchart showing the solution search process [Figure 6](a) is an explanatory diagram of the average daily water volume W1, and (b) is an explanatory diagram of the average daily water volume W2 by group. [Figure 7] (b) is an explanatory diagram of the pipeline network under analysis, where intersections are grouped. (b) is an explanatory diagram of the relationship between the measured daily average water volume extracted from each group's extraction intersection and the group's daily average water volume W2. [Figure 8] Diagram illustrating the daily water volume W3 by group. [Figure 9] (a) is an explanatory diagram of past daily water volume data that serves as the basis for calculating the group-specific water volume correction coefficient S2, and (b) is an explanatory diagram of the group-specific water volume correction coefficient S2 when the group is a school. [Figure 10] (a) is a table explaining the relationship between the total amount of water extracted per unit time by group W3 and the amount of water flowing into the pipeline network, and (b) is a water flow diagram showing the relationship between the total amount of water extracted per unit time by group W3 before correction and the amount of water flowing into the pipeline network per unit time Qt. [Figure 11] This diagram illustrates the relationship between the total amount of water extracted per unit time by group (W3) after correction, and the amount of water flowing into the pipeline network per unit time (Qt). [Figure 12] (a) is an explanatory diagram of the correction coefficient, and (b) is a table explaining the relationship between the total value of the group-specific water extraction volume W3 per unit time after correction and the water inflow volume into the pipeline network. [Figure 13] Diagram illustrating the analysis results when the water extraction volume correction coefficient K is changed. [Figure 14] (a) is an explanatory diagram of the scoring table used to evaluate the degree of deviation, and (b) is an explanatory diagram of the placement of sensors installed in the pipeline network. [Figure 15] (a) and (b) are diagrams illustrating the algorithm for deriving the optimal correction coefficient K for each group. [Figure 16] Diagram illustrating the algorithm for deriving the optimal correction coefficient K for each group. [Modes for carrying out the invention]

[0028] Below, an example of the hydraulic analysis method and hydraulic analysis apparatus according to the present invention will be described based on the drawings. [Configuration of the hydraulic analysis device] The hydraulic analysis device 1 consists of a general-purpose personal computer with a motherboard equipped with a CPU, a memory board equipped with semiconductor memory, and other components, to which storage devices such as hard disks and SSDs, a touch-panel LCD display, and input / output devices such as a keyboard and mouse are connected.

[0029] The storage device has an OS program installed to manage the system, and a hydraulic analysis program is further installed as an application program executed by the CPU under the management of the OS program. The hydraulic analysis program is stored and distributed on recording media consisting of non-volatile optical discs such as CD-ROMs and DVD-ROMs, or non-volatile semiconductor memory such as USB memory, or it is provided via a network from a cloud server and installed on a personal computer to function as the hydraulic analysis device 1.

[0030] Figure 1 shows the functional blocks of the hydraulic analysis device 1, which are realized by the CPU and the hydraulic analysis program executed on the CPU. The hydraulic analysis device 1 consists of a calculation processing unit 2, a data storage unit 3, a display unit 4 using a liquid crystal display, and an input unit 5 equipped with a keyboard and mouse. It is configured to perform hydraulic analysis on a water distribution pipeline network that connects water reservoirs and individual water demand points with multiple pipelines and to calculate pipeline network characteristic values ​​at any given water demand point. Pipeline network characteristic values ​​refer to characteristic values ​​that allow us to understand the movement and state of water in the pipeline network, and include values ​​such as water pressure, flow direction, flow velocity, water volume, and residual chlorine concentration.

[0031] The calculation processing unit 2 comprises an analysis job management unit 20, an analysis condition setting unit 21, an analysis execution unit 22, and an analysis result evaluation unit 23. The analysis condition setting unit 21 is a functional block that sets analysis conditions for the pipeline network, the analysis execution unit 22 is a functional block that performs hydraulic analysis based on the analysis conditions and outputs the analysis results, and the analysis result evaluation unit 23 is a functional block that compares and evaluates the analysis results with measured values. The analysis job management unit 20 is a functional block that controls the operation of each functional block, the analysis condition setting unit 21, the analysis execution unit 22, and the analysis result evaluation unit 23, and oversees the hydraulic analysis process.

[0032] The analysis condition setting unit 21 includes an inflow / outflow intersection definition processing unit 21A, a daily average water extraction volume setting processing unit 21B, a grouping processing unit 21C, a group-specific daily average water extraction volume setting processing unit 21D, a group-specific water extraction volume correction coefficient acquisition processing unit 21E, a group-specific water extraction volume calculation processing unit 21F, a unit time water extraction volume correction processing unit 21G, and a water extraction volume correction coefficient calculation processing unit 21H, among others.

[0033] The inflow / outflow intersection definition processing unit 21A is a functional block that defines the inflow intersections where water flows into the pipeline network and the outflow intersections where water is taken out. The daily average outflow water volume setting processing unit 21B sets the unit time inflow water volume Qt(m³) measured at the inflow intersection. 3 This is a functional block that sets the average daily inflow (=Q / 24), which is calculated from the total inflow Q (=ΣQt) obtained by adding up the daily amounts of ), to the average daily outflow W1 at the outflow intersection.

[0034] The grouping processing unit 21C is a functional block that divides multiple extraction intersections into groups that share a common water extraction trend. The group-specific daily average water extraction volume setting processing unit 21D is a functional block that sets the group-specific daily average water extraction volume W2 (=W1 × S1) by allocating the daily average water extraction volume W1 to each group based on the previously acquired ratio S1 of daily average water extraction volumes for each group.

[0035] The group-specific water volume extraction correction coefficient acquisition processing unit 21E is a functional block that acquires the group-specific water volume extraction correction coefficient S2 for each unit of time from the time variation patterns of the water volume extraction for each group that are measured in advance for at least one set of analysis conditions including season, day of the week, and event, and that match the target analysis condition. The group-specific water volume extraction calculation processing unit 21F is a functional block that calculates the group-specific water volume extraction W3 (=W2 × S2) for each unit of time from the group-specific daily average water volume extraction W2 and the group-specific water volume extraction correction coefficient S2.

[0036] The unit time water volume extraction processing unit 21G is a functional block that calculates and applies a unit time water volume extraction coefficient m to correct the sum of the group-specific water volume extraction W3 for each unit time so that it matches the unit time inflow water volume Qt. The water volume extraction coefficient calculation processing unit 21H is a functional block that sets the product S2·m of the group-specific water volume extraction coefficient S2 and the unit time water volume extraction coefficient m as the water volume extraction coefficient K.

[0037] The analysis execution unit 22 outputs the analyzed values ​​of water pressure, flow rate, or flow direction at the extraction intersection by performing a hydraulic analysis under analysis conditions in which the extraction water volume correction coefficient (K) has been adjusted, and by performing a process multiple times to output the analyzed values ​​of water pressure, flow rate, or flow direction at the extraction intersection. The analysis result evaluation unit 23 includes a deviation degree evaluation processing unit 23A that evaluates the degree of deviation of the analyzed values ​​obtained by the analysis execution unit 22 from the measured values, and an analysis condition determination processing unit 23B that determines the optimal analysis conditions from among multiple analysis conditions based on the deviation degree evaluation results.

[0038] The data storage unit 3 includes a pipeline network diagram storage area 30 in which a pipeline network diagram representing the pipeline network connecting from the water storage tank to each water demand point is stored, a pipeline information storage area 31 in which component information (pipe type, diameter, pipeline extension, etc.) and laying information (laying date, constructor, etc.) of pipelines and other components constituting the pipeline network diagram are stored, a group information storage area 32 in which conditions for grouping and dividing each pipeline and intersection constituting the pipeline network diagram and group information of the pipelines grouped and divided according to the conditions are stored, an analysis result storage area 33 in which the analysis results by the analysis execution unit 22 are stored, and a measured data storage area 34 in which inflow data flowing into the inflow intersection from the water storage tank and measured data such as water pressure and flow rate at the extraction intersection which is a water demand point are stored, etc.

[0039] In Fig. 2(a), a water pipeline network 6 is illustrated. The water pipeline network 6 is composed of a plurality of pipelines 8 and a plurality of intersections 9 where each pipeline 8 intersects. Water supplied from the upstream water storage tank 7 flows into the inflow intersection 9A, and water is taken out from a plurality of extraction intersections 9B through the pipeline 8. In the water pipeline network 6, the connection point between pipelines which is a water demand point is referred to as a node (synonymous with intersection).

[0040] When calculating the water head of each node by a water analysis using the node head method in the analysis execution unit 22, the Hazen-Williams formula H = 10.666×(L×Q 1.85 ) / (C 1.85 ×d 4.87 ) and the flow equation which is the flow continuity condition equation at the node (water demand point 9B) illustrated in Fig. 2(b) Σ(±Q ij ) = P i and the closed pipeline equation shown in Fig. 2(c) Σ(±H ks ) - δE k = 0 are obtained as simultaneous solutions.

[0041] Here, H is the frictional head loss of the pipe (m), L is the pipeline extension (m), Q is the flow rate (m 3 / s), d is the actual inner diameter of the pipe (m), C is the velocity coefficient, and Q iP is the flow rate in each pipe connected to the node of interest. i is the water supply from the node. Furthermore, the closed pipe equation states that the water in the pipeline network must be such that the total energy loss is minimized; that is, in a pipeline network with J pipelines, ΣQ j H j The flow is minimized as j=1~J. The flow equation is used as a constraint: ΣQ j H j →By solving for min, the closed-pipe equation can be obtained.

[0042] The value of the flow velocity coefficient C remains constant regardless of the time of day or the amount of water, and varies depending on the roughness of the inner wall of the pipe. For example, in the case of cast iron pipes, it is 130-140 for newer pipes, and decreases to 60-70 for older pipes with rust deposits on the inner wall.

[0043] By reading necessary information such as pipeline length, pipe diameter, and height of each node including the reservoir from the pipeline information storage area 31, and setting the amount of water flowing into the inflow intersection and the amount of water taken out from the water demand point (outlet intersection), the flow direction of each pipeline is determined, and the pipeline 8 necessary for distributing water from the reservoir 7 to the node (water demand point 9B) is identified.

[0044] The group information storage area 32 stores attributes for dividing multiple water extraction intersections included in the pipeline network into groups with common water extraction trends. Figure 9(a) illustrates the trend of daily water demand. For example, attribute information such as schools where water demand occurs mainly during the daytime on weekdays, commercial areas where demand occurs from daytime to nighttime on both weekdays and weekends, industrial areas where demand occurs mainly on weekdays, water tank type residences where demand occurs intermittently, and direct pressure type residences where demand peaks in the morning and evening are stored along with measured values ​​of the amount of water extracted for each season, including weekdays, weekends, and days when events such as baseball games are held.

[0045] [Structure of the hydraulic analysis method] The procedure for the hydraulic analysis method performed using the hydraulic analysis device 1 described above will now be explained. The hydraulic analysis method according to the present invention is a hydraulic analysis method that performs a hydraulic analysis on a water distribution pipeline network 6, which connects a water reservoir 7 and individual water demand points, which are extraction intersections 9B, with multiple pipelines 8, in order to calculate pipeline network characteristic values ​​at any water demand point. When the hydraulic analysis program is started and the analysis job management unit 20 is started up, a series of analysis processes are executed.

[0046] First, an input guidance screen is displayed via the analysis condition setting unit 21 to set the necessary analysis conditions. The operator enters the necessary information according to the instructions on the input guidance screen to set the analysis conditions. Then, the analysis execution unit 22 is started and the analysis is performed. The analysis result evaluation unit 23 evaluates the analysis results, and based on the results, some of the analysis conditions are changed, and the analysis is performed again by the analysis execution unit 22. The analysis and evaluation are repeated multiple times until the desired analysis result is obtained.

[0047] The analysis job management unit 20 performs the following steps: an analysis condition setting step, which is executed by the analysis condition setting unit 21 to set analysis conditions for the pipeline network; an analysis execution step, which is executed by the analysis execution unit 22 to perform hydraulic analysis based on the analysis conditions and output analysis results; and an analysis result evaluation step, which is executed by the analysis result evaluation unit 23 to compare and evaluate the analysis results with measured values. These steps are described in detail below.

[0048] As shown in Figure 3, the analysis condition setting step consists of an inflow / outflow intersection definition step (SA1), a daily average water extraction rate setting step (SA2), a grouping step (SA3), and a group-specific daily average water extraction rate setting step (SA4).

[0049] In the inflow / outflow intersection definition step (SA1), inflow intersections where water flows from reservoirs into the pipeline network and outflow intersections where water is extracted at the point of demand are defined. In this example, it is assumed that a portion of the entire pipeline network is pre-selected as blocks, and that water is supplied to the block pipeline network from a single inflow intersection and extracted from multiple outflow intersections.

[0050] In the daily average water extraction setting step (SA2), as shown in Figure 6(a), the daily average water extraction (=Q / 24) obtained from the total water extraction Q, which is the sum of the hourly water extraction Qt measured at each of the 24 hours of the day at the inflow intersection, is set as the daily average water extraction (W1) at the extraction intersection. The daily average water extraction (W1) is the sum of the daily average water extraction amounts from each extraction intersection of the block.

[0051] In the grouping step (SA3), each extraction intersection is grouped into groups that share a common water extraction trend, specifically, groups that share a common pattern of water extraction volume as explained in Figure 9(a). Figure 7(a) shows an example of a block pipeline network where each group is distinguished by a different color (shades of color on the diagram) so that the grouped intersections are the same color, and Figure 7(b) shows examples of the names of the groups to be divided.

[0052] In the group-specific daily average water extraction setting step (SA4), as shown in Figure 6(b), the daily average water extraction W1 is allocated (prorated) to each group based on the ratio S1 of the previously acquired daily average water extraction amounts for each group (existing information), and this value is set as the group-specific daily average water extraction amount W2.

[0053] Specifically, the average daily water volume for each group is read from the past outflow volume for each group (exemplified in Figure 9(a)) that matches the season, day of the week, etc., stored in the actual data storage area 34 and is the subject of analysis. Based on the ratio S1 of the average daily water volume (existing information) calculated based on the ratio of the average daily water volume (existing information) as shown in Figure 7(b), the average daily water volume (W1) is apportioned and set as the average daily water volume for each group. For example, the total value of the known average daily water volume (236.3 m³) 3 Maintain the ratio S1 of the daily average water volume (12.3, 37.9, ...) of each group to the total daily average water volume (124.4 m³) 3 The average daily water volume (6.6, 19.89, ...) of each group is allocated proportionally to the total.

[0054] Next, the analysis condition setting unit 21 performs the group-specific water extraction volume correction coefficient setting step (SA5) and the group-specific water extraction volume calculation step (SA6). In the group-specific water volume correction coefficient setting step (SA5), a group-specific water volume correction coefficient S2 for each unit time, corresponding to 24 hours in a day, is obtained for each group from the daily time variation pattern of the water volume extracted for each group, which has been measured in advance (see Figure 9(a)). Here, the 24-hour average of the water volume correction coefficient S2 is set to 1. Figure 9(b) shows the water volume correction coefficients for the group of schools.

[0055] In the group-specific water extraction calculation step (SA6), the group-specific water extraction amount W3 (=W2 × S2) per unit time is calculated from the group-specific daily average water extraction amount W2 and the group-specific water extraction amount correction coefficient S2. Figure 8 is set up so that the group-specific water extraction amount W3 obtained by multiplying the group-specific daily average water extraction amount by the group-specific water extraction amount correction coefficient S2 at each time is shown as a bar graph of varying shades.

[0056] When comparing the sum of the hourly group-specific water extraction volume W3 shown in Figure 8, which was set in this manner, with the unit time inflow volume Qt shown in Figure 6(a), as shown in Figures 10(a) and (b), the summation average of the daily average inflow volume W1 and the group-specific daily average water extraction volume W2 both match at 124.4, but the sum of the unit time inflow volume Qt and the group-specific daily average water extraction volume W2 for each hour do not match.

[0057] Therefore, the analysis condition setting unit 21 calculates the unit time water volume correction step (SA7) and the water volume correction coefficient calculation step (SA8). In the unit time water extraction rate correction step (SA7), a unit time water extraction rate correction coefficient m is calculated to correct the sum of the group-specific water extraction rates W3 per unit time so that it matches the unit time inflow rate, and the sum of the group-specific water extraction rates W3 per unit time is corrected. In the group correction coefficient calculation step (SA8), the product of the group-specific water extraction rate correction coefficient S2 and the unit time water extraction rate correction coefficient m is calculated as the water extraction rate correction coefficient K.

[0058] As a result, as shown in Figure 11, the unit time water extraction rate correction coefficient m can be used to match the unit time water inflow rate with the total value of the group-specific water extraction rate W3 per unit time, and the setting of the analysis conditions by the analysis condition setting unit 21 is completed. Figure 12(a) shows an example of the state before and after correction of the total value of the group-specific water extraction rate W3 per unit time at 8 o'clock.

[0059] However, as shown in Figure 12(b), the average value of the group-specific water extraction volume W3 at each time point and the daily average inflow volume Q AVR Based on this, some groups will have differences in the average daily water extraction volume (W2) values ​​set for each group that do not fall within the acceptable range. Therefore, in order to bring this difference within the acceptable range, the analysis execution step is performed, and the analysis conditions are adjusted to adjust the water extraction volume correction coefficient K, and the process of outputting the water pressure, flow rate, or flow direction, which are the analysis values ​​at the extraction intersection, is performed multiple times (SA9). In the deviation evaluation step, the deviation of the analysis values ​​from the measured values ​​is evaluated, and in the analysis condition determination step, the optimal analysis conditions are determined from among the multiple analysis conditions based on the deviation evaluation results (SA10). Such processing is performed for each of the required analysis target days.

[0060] Figure 4 shows the detailed procedures for the analysis execution step and the deviation evaluation step shown in step SA9. First, the number of loops for a day, 24 times, from 0:00 to 23:00 is set (SB1), and the number of analyses to be repeated for each loop is set (SB2). Next, the analysis conditions obtained by adopting the extraction water volume correction coefficient K obtained above are set (SB3), and the analysis is executed by the analysis execution unit 22 (SB4), and analysis results such as flow velocity, water pressure, and flow direction are obtained (SB5).

[0061] The analysis result evaluation unit 23 evaluates the degree of deviation by comparing the analysis results with the measured values ​​(SB6). If the degree of deviation does not converge to the set value, SB10 (SB7, N) is performed, the process proceeds to step SB10, the loop count is updated (added by 1), the extraction water volume correction coefficient K is adjusted, and then the process from step SB3 onwards is repeated. For example, the gradient descent method can be suitably used as the algorithm for updating the extraction water volume correction coefficient K.

[0062] Figure 13 shows the analysis results for multiple patterns when the extraction water flow rate correction coefficient K is changed under the condition that the unit time inflow water flow rate Qt is satisfied.

[0063] In step SB7, if the predetermined number of loops set in step SB2 is reached, or if it is determined that the deviation has converged to the set value, a list of analysis results is generated and stored in the memory unit (SB8). The process from step SB2 to step SB8 is repeated 24 times from 0:00 to 23:00 and then terminates (SB9).

[0064] Figure 5 shows the details of the process in step SB10 of Figure 4. Based on the analysis results for 24 hours, analysis results with the smallest discrepancy (score) between measured values ​​and analyzed values ​​for each hour are selected (SC1), and it is evaluated whether the daily average water volume extracted meets the acceptable range for all groups (SC2).

[0065] If the average daily water volume extracted meets the acceptable range for all groups, the process is terminated (SC3,Y). If it does not meet the range (SC3,N), the analysis result that does not meet the range is replaced with another analysis result, and the process from step SC1 is repeated.

[0066] The formula for calculating the score used to evaluate the degree of deviation is as follows: Score = Deviation / Allowable Deviation Degree of deviation = Σ(analytical value - measured value) 2 Allowable deviation = Σ (sensor measurement accuracy) 2 A low score indicates that the value matches the actual value, while a high score indicates that the value deviates from the actual value. The above scores are calculated based on actual values ​​from water pressure and flow velocity sensors installed at multiple locations in the pipeline network. The evaluation is based on the individual scores for water pressure and flow velocity, as well as the sum of these scores to arrive at a comprehensive score. The comprehensive score may also be obtained by weighted addition, where certain values ​​are given greater weight.

[0067] Figure 14(a) shows an example of the score. The numbers indicating the measurement points correspond to the numbers of the outlet intersections in the pipeline network in Figure 14(b).

[0068] From the results of multiple analysis patterns, a combination is selected that has a low score throughout the 24-hour period and satisfies the acceptable range for the average water volume extracted. By setting such values ​​as analysis conditions, an accurate state of the entire pipeline network can be obtained.

[0069] Below, we illustrate the gradient descent method, an example of an algorithm that expresses the degree of deviation between the analysis and the measured value as a score, and derives the optimal solution for the correction coefficient K for each group that minimizes the score. Gradient descent is one of the mathematical optimization methods used to find the minimum (or maximum) value of a function, and a specific example is explained below.

[0070] As shown in Figure 15(a), in Step 1, the analysis is performed using the correction coefficient K = 0.92 shown in the circled number 1 in No. I, the score is calculated and entered in the table, and the amount of change in the correction coefficient K is set. In Step 2, as shown in Figure 15(b), the analysis is performed under the condition that the amount of change is added to only one of the groups, and a score is calculated to obtain eight different patterns. The amount of change is not particularly limited.

[0071] In Step 3, as shown in Figure 16, the score change (the difference from the score shown in the circled number 0 in No. I) is calculated, and the amount of change in the correction coefficient K is calculated according to the ratio of the change. In Step 4, the amount of change in the correction coefficient is added to the correction coefficient of No. I, and the process returns to Step 1 to analyze it as No. II. In this way, a combination with a small score that satisfies the acceptable range of average water volume is determined.

[0072] In the example described above, the configuration of the computer constituting the hydraulic analysis device is not specifically explained, but it may be configured as a standalone system, or it may be configured to collect values ​​from sensors that detect water pressure and flow velocity placed at key points in the pipeline network using IoT technology, store them in a database of actual measurement data, and perform the hydraulic analysis described above based on the actual measurement data stored in the database. Alternatively, the hydraulic analysis device may be configured as a cloud computer, and an API for operating the cloud computer may be ported to a terminal computer connected to the network, so that hydraulic analysis can be performed on the cloud computer via remote operation.

[0073] By employing the hydraulic analysis device and method described above, it becomes possible to easily perform remote monitoring of water distribution blocks constituting the pipeline network, optimize water distribution pressure, eliminate low flow velocity areas, maintain residual chlorine concentration, and estimate the risk of turbidity generation.

[0074] For example, by connecting the aforementioned cloud computer and terminal computer via a network, it becomes possible to utilize multi-point measurement data and pipe network analysis results based on that data, enabling high-precision remote monitoring of water pressure, flow direction, and flow velocity across the entire pipe network. Mobile computers such as smartphones can also be used as terminal computers. For example, by specifying the water distribution block to be analyzed and the analysis date on the terminal computer, appropriate analysis conditions can be automatically set on the server computer, the analysis can be executed, and the results can be viewed on the terminal computer.

[0075] Furthermore, it becomes possible to predict water pressure, flow rate, or flow direction at each intersection of the pipeline network with high accuracy and over a wide area. This allows for the estimation of locations with abnormal water pressure or flow velocity, taking into account the daily and temporal variations in water pressure, flow rate, or flow direction, enabling proper water management (management of water pressure and flow rate). For example, if the water pressure or flow rate obtained from the analysis deviates from a predetermined tolerance range, an anomaly detection program that outputs information identifying the corresponding location and abnormal values ​​can be incorporated into the hydraulic analysis device.

[0076] Furthermore, for pipelines with chronically low flow velocities, the flow velocity can be equalized by valve operation or switching of the water distribution system. Therefore, when the water pressure or flow rate obtained from the analysis consistently shows values ​​near the upper or lower limit of the allowable range, a water pressure adjustment function should be incorporated into the hydraulic analysis device to identify the flow control valves or switching valves to supply water to the relevant pipeline and output operational information. In this case, by performing the analysis in advance, it is also possible to incorporate a function that shows what value the flow velocity of the target pipeline will be adjusted to as a result of operating the flow control valves or switching valves. Conversely, it is also possible to have a function that identifies the flow control valves or switching valves necessary to adjust the flow velocity of the target pipeline to a target value and presents the operational information.

[0077] Similarly, by using hydraulic analysis to identify pipelines where low flow velocity or frequent changes in flow direction occur, it becomes possible to identify locations where stagnation may lead to low residual chlorine concentrations. This also enables the identification of flow control valves and switching valves necessary to resolve these conditions, and provides information on their operation.

[0078] Furthermore, by identifying pipelines with low daily maximum flow velocities and the potential for sediment accumulation through hydraulic analysis, it becomes possible to estimate locations where turbidity is likely to occur when flow direction changes or flow velocity increases. Therefore, it is possible to manage the analysis results over time, implement a history monitoring function that detects when flow direction changes or flow velocity increases exceeding preset values, and provide an alarm function that identifies pipelines at risk of turbidity and alerts the administrator. Each of the above functions can be realized by incorporating them into the hydraulic analysis program in advance.

[0079] If the analysis job management unit 20 is equipped with a function to automatically execute the hydraulic analysis method described above under various conditions such as weekdays, holidays, and event days, and to store the corresponding analysis conditions in the data storage unit 3, then when an analysis is needed, it will be possible to immediately perform the analysis based on the analysis conditions stored in the data storage unit 3.

[0080] The embodiments described above are merely one embodiment of the present invention, and the scope of the present invention is not limited by this description. The specific configuration of the hydraulic analysis device can be appropriately modified and designed within the scope that achieves the effects of the present invention. [Explanation of symbols]

[0081] 1:Hydraulic analysis equipment 2: Arithmetic section 3: Data storage unit 20: Analysis Job Management Department 21: Analysis Condition Setting Section 21A: Inflow / Outflow Intersection Definition Processing Unit 21B: Daily average water volume setting processing unit 21C: Grouping Processing Unit 21D: Group-specific daily average water volume setting processing unit 21E: Group-specific water volume extraction correction coefficient acquisition processing unit 21F: Group-specific water volume calculation processing unit 21G: Unit Time Extraction Water Volume Correction Coefficient Calculation Processing Unit 21H: Extraction water volume correction coefficient calculation processing unit 22: Analysis Execution Unit 23: Analysis Result Evaluation Department 23A: Deviation degree evaluation processing unit 23B: Analysis Condition Determination Processing Unit

Claims

1. A hydraulic analysis method for a pipeline network consisting of multiple pipelines and multiple intersections where each pipeline intersects, The process includes: an analysis condition setting step for setting analysis conditions for the pipeline network; an analysis execution step for performing hydraulic analysis based on the analysis conditions and outputting analysis results; and an analysis result evaluation step for comparing the analysis results with measured values ​​and evaluating them. The aforementioned step of setting analysis conditions is: The inflow / outflow intersection definition step defines inflow intersections where water flows into the pipeline network and outflow intersections where water is taken out, The average daily inflow (Q) is calculated from the total inflow (Q) obtained by adding the daily unit time inflow (Qt) corresponding to the day being analyzed, which is measured at the aforementioned inflow intersection. ＡＶＲ A daily average water volume setting step in which the daily average water volume (W1) at the aforementioned water volume intersection is set, A grouping step in which the aforementioned extraction intersections are divided into groups that share a common water extraction trend, A group-specific daily average water extraction amount setting step, in which the daily average water extraction amount (W1) is allocated to each group based on the ratio (S1) of the daily average water extraction amount for each group obtained in advance, and the group-specific daily average water extraction amount (W2) is set as the value obtained by allocating the daily average water extraction amount (W1) to each group, A hydraulic analysis method that includes [details omitted].

2. The aforementioned step of setting analysis conditions is: A group-specific water volume correction coefficient setting step, which obtains a group-specific water volume correction coefficient (S2) for each unit time from the time variation pattern of the water volume extracted for each group that is similar to the analysis target day, among the daily time variation patterns of the water volume extracted for each group that have been measured in advance for each analysis condition including at least season, day of the week, and event, A group-specific water volume calculation step that calculates the group-specific water volume (W3) per unit time from the group-specific average daily water volume (W2) and the group-specific water volume correction coefficient (S2), A hydraulic analysis method according to claim 1, comprising:

3. The aforementioned step of setting analysis conditions is: A unit time water extraction correction step involves calculating and correcting a unit time water extraction correction coefficient (m) that corrects the sum of the group-specific water extraction amounts (W3) per unit time so that it matches the unit time inflow amount (Qt), A hydraulic analysis method according to claim 2, including the following:

4. The aforementioned step of setting analysis conditions is: The hydraulic analysis method according to claim 3, comprising a step of calculating an extraction water volume correction coefficient, in which the product of the group-specific extraction water volume correction coefficient (S) and the unit time extraction water volume correction coefficient (m) is set as the extraction water volume correction coefficient (K).

5. The aforementioned analysis execution step is: The hydraulic analysis method according to claim 4, wherein the process of outputting water pressure, flow rate, or flow direction, which are analytical values ​​at the extraction intersection, is performed multiple times by performing the hydraulic analysis under analytical conditions in which the extraction water volume correction coefficient (K) has been adjusted.

6. The aforementioned analysis result evaluation step is: A deviation degree evaluation step is performed to evaluate the degree of deviation of the analysis value obtained in the analysis execution step with respect to the measured value, An analysis condition determination step in which the optimal analysis condition is determined from among the multiple analysis conditions based on the evaluation result of the degree of deviation, A hydraulic analysis method according to claim 5, including the following:

7. A hydraulic analysis device for a pipeline network consisting of multiple pipelines and multiple intersections where each pipeline intersects, It comprises: an analysis condition setting unit that sets analysis conditions for the pipeline network; an analysis execution unit that performs hydraulic analysis based on the analysis conditions and outputs analysis results; and an analysis result evaluation unit that compares the analysis results with measured values ​​and evaluates them. The aforementioned analysis condition setting unit is: An inflow / outflow intersection definition processing unit defines inflow intersections where water flows into the pipeline network and outflow intersections where water is taken out, The average daily inflow (Q) is calculated from the total inflow (Q) obtained by adding the daily unit time inflow (Qt) corresponding to the day being analyzed, which is measured at the aforementioned inflow intersection. ＡＶＲ A daily average water volume setting processing unit sets the daily average water volume (W1) at the aforementioned water volume intersection, A grouping processing unit that divides the aforementioned extraction intersections into groups that share a common water extraction trend, A group-specific daily average water extraction volume setting processing unit sets the value obtained by allocating the daily average water extraction volume (W1) to each group based on the ratio (S1) of the daily average water extraction volume for each group, which has been acquired in advance, as the group-specific daily average water extraction volume (W2). A hydraulic analysis device that includes [specific components / devices].

8. The aforementioned analysis condition setting unit is: A group-specific water volume correction coefficient acquisition processing unit acquires a group-specific water volume correction coefficient (S2) for each unit time from the time variation pattern of the water volume extracted for each group that matches the day of analysis, among the time variation patterns of the water volume extracted for each group that have been measured in advance for each analysis condition including at least season, day of the week, and event, and A group-specific water extraction volume calculation processing unit calculates the group-specific water extraction volume (W3) per unit time from the group-specific average daily water extraction volume (W2) and the group-specific water extraction volume correction coefficient (S), A hydraulic analyzer according to claim 7, including the following:

9. The aforementioned analysis condition setting unit is: A unit time water volume correction processing unit calculates and corrects a unit time water volume correction coefficient (m) that corrects the sum of the group-specific water volume (W3) per unit time to match the unit time inflow water volume (Qt), A hydraulic analyzer according to claim 8, including the following:

10. The aforementioned analysis condition setting unit is: The hydraulic analysis apparatus according to claim 9, which includes a water extraction volume correction coefficient calculation processing unit that sets the product of the group-specific water extraction volume correction coefficient (S) and the unit time water extraction volume correction coefficient (m) as the water extraction volume correction coefficient (K).

11. The aforementioned analysis execution processing unit, The hydraulic analysis apparatus according to claim 10, wherein the hydraulic analysis is performed multiple times under analysis conditions in which the extraction water volume correction coefficient (K) has been adjusted, thereby outputting the analysis values ​​of water pressure, flow rate, or flow direction at the extraction intersection.

12. The aforementioned analysis result evaluation unit is: A deviation degree evaluation processing unit evaluates the degree of deviation of the analysis value obtained by the analysis execution processing unit from the measured value, An analysis condition determination processing unit that determines the optimal analysis condition from among the plurality of analysis conditions based on the evaluation result of the degree of deviation, A hydraulic analyzer according to claim 11, including the following:

Images