JP3471449B2 - Pipe network loop selection method, initial value determination method, and pipe network calculation method - Google Patents

Description

Detailed Description of the Invention

[0001]

BACKGROUND OF THE INVENTION The present invention relates to the calculation of pipe networks.
The present invention relates to a method for selecting the loop, a method for setting an initial value, and a method for calculating a pipe network.

[0002]

2. Description of the Related Art Conventionally, in water supply business, gas business, etc.
When installing a new pipe or considering stable supply, the flow rate of the pipe network was calculated. Also, in the electric power industry, the electric wire network current was similarly calculated.

Details of such pipe network calculation are described in, for example, "Design Method and Calculation Examples of Water Supply and Distribution Pipelines, Hyundai Science and Engineering Publishing, Nobuo Matsuda, February 10, 1972, p. 81". In such a pipe network calculation, the calculator specifies multiple loops for the target pipe network, sets the initial flow rate of each pipe, and calculates for each loop by the calculation method described in the previous work. To do. This calculation is easy to program and is calculated by a calculator.

However, since the selection of the loop is done manually, it takes a lot of time to extract and set the loop from the complicated pipe network, and there are problems such as wasteful calculation due to an error.

Since the setting of the initial flow rate is also based on human intuition, there is a problem that when the initial setting value is bad, the number of times of calculation is large and the pipe network calculation takes a long time.

An invention for solving this problem and automating it is being filed by the present inventor in Japanese Patent Application No. 6-203008.

[0007]

However, in both the flow rate method and the energy method used in pipe network analysis, calculation is performed based on the relation table of pipes and nodes, and these pipes and nodes are manually input. When I did it, I was excluding isolated pipes. However, when pipe network analysis is performed by specifying the range of the target area using the mapping system, an isolated pipe may occur due to the specified range, and if there is an isolated pipe, the solution becomes an indeterminate problem. It will be impossible. Therefore, the range and the like must be entered with care so that an isolated pipe is not generated, and there is a problem that the operator's burden is heavy.

The present invention has been made in view of the above problems, and an object of the present invention is to provide a pipe network capable of automatically extracting an isolated pipe, setting a loop of the pipe network, and setting an initial value. It is to provide a calculation method.

[0009]

According to a first aspect of the present invention, there is provided a computer for loop selection of a pipe network in which at least one source node and a plurality of nodes are connected by a pipe having resistance. A loop selection method using (a) extracting and eliminating isolated pipes from the pipe network, and (b) based on the total resistance of the path from the supply source node to each node, An upstream node and a downstream node of the pipe, a step of determining the upstream pipe of each node, (c) a step of using a pipe not designated as an upstream pipe as a loop source pipe in the step, and (d) Ignoring the loop source pipe of the pipe network,
Determining a node order indicating how much downstream each node is with respect to the source node, based on an upstream node and a downstream node of each node; A loop source pipe as an origin, and a step of selecting a loop based on the node order of the upstream side node and the node order of the downstream side node of each pipe, wherein the step (a) includes flags of all the nodes, An initial value is given to the flag of the pipe of, the supply source node, the flag of the pipe connected to the supply source node, and the flag of the other node of the pipe are changed to a value other than the initial value. , For the node whose flag has changed, changing the flag of the pipe connected to the node and the flag of the other node of the pipe is repeated, and the flag of the pipe is initialized. A loop Selection of the tube network and extracting the pipe remains as an isolated pipe.

A second aspect of the present invention is a method for determining an initial value which uses a computer to determine an initial value of a pipe network in which at least one source node and a plurality of nodes are connected by a pipe having resistance. a) extracting and eliminating isolated pipes from the pipe network; and (b) based on the total resistance of the path from the source node to each node, the upstream node and the downstream node of each pipe, and In the step of determining the upstream pipe of the node, and (c) in the above step,
(D) ignoring the loop source pipe of the pipe network, and making each node downstream of the source node, a pipe not designated as an upstream pipe among all pipes; Determining the node order indicating whether the node is a side node based on the upstream side node and the downstream side node of each node, and (e) the node order of the upstream side node of each pipe with the loop source pipe as an origin. , Selecting a loop based on the node order of downstream nodes,
(F) The flow rate of the loop source pipe is set to “0”, the pipe network has a tree structure, the flow rate of the pipe supplied to the node is determined from the demand amount of the terminal node, and the demand amount of each node is added to the upstream side. Determining the flow rate of the pipe on the side and setting the flow rate to the initial value of the pipe flow rate, wherein the step (a) gives the initial values to the flags of all the nodes and the flags of all the pipes, and the source node And, changing the flag of the pipe connected to the supply source node and the flag of the other node of the pipe to a value other than the initial value, of the other node,
For a node whose flag has changed, changing the flag of the pipe connected to the node and the flag of the other node of the pipe is repeated, and the pipe with the initial value remains unchanged. It is a method for determining the initial value of the pipe network, which is characterized by extracting as.

A third aspect of the present invention is a calculation method for using a computer to calculate a pipe network in which at least one source node and a plurality of nodes are connected by a pipe having resistance, and (a) the pipe Extracting and eliminating isolated pipes from the network, and (b) based on the total resistance of the path from the supply source node to each node, the upstream and downstream nodes of each pipe, and the upstream side of each node. Determining a pipe, (c) in the above process, using a pipe not designated as an upstream pipe among all pipes as a loop source pipe, and (d) ignoring the loop source pipe of the pipe network, (E) determining a node order indicating which node is the most downstream node with respect to the source node, based on the upstream node and the downstream node of each node; The step of selecting a loop based on the node order of the upstream side node and the node order of the downstream side node of each pipe with the loop source pipe as the origin, and (f) setting the flow rate of the loop source pipe to "0" A tree structure, determining the flow rate of the pipe supplied to the node from the demand amount of the end node, determining the flow amount of the upstream pipe while adding the demand amount of each node, and the initial value of the pipe flow rate, (G) performing a pipe network calculation using the initial value obtained in the step (f) in each loop selected in the step (e), and all the steps (a) are Give the initial values to the flags of all nodes and the flags of all the pipes, and change the source node, the flag of the pipe connected to the source node, and the flag of the other node of the pipe to values other than the initial values. , Above Out of the node,
For a node whose flag has changed, changing the flag of the pipe connected to the node and the flag of the other node of the pipe is repeated, and the pipe with the initial value remains unchanged. It is a calculation method of a pipe network characterized by extracting as.

[0012]

In the first aspect of the present invention, the isolated pipe is extracted and eliminated from the pipe network, the pipe not designated as the upstream pipe for each node is used as the loop source pipe, and the loop source pipe is removed from the pipe network. In, the node order from the supply source is determined, and the loop is selected based on this node order. Selection of such a loop can be performed using a computer.

[0013]

In the second invention, in the pipe network calculation of the first invention, the flow rate of a specific pipe is set to "0", the pipe network is made into a tree structure, and the demand amount of a terminal node is supplied to that node. The flow rate of the pipe is determined, the flow rate of the upstream pipe is determined while adding the demand amount of each node, and is set as the initial value of the pipe flow rate. Here, this processing can be performed by a computer.

The third invention is the same as the second invention, except that step (d)
A combination of the steps of performing the pipe network calculation using the initial value obtained in the step (e) in each loop selected in step (e) can be performed by a computer.

Here, the pipe network is a gas pipe network, a water supply pipe network, or an electric power network. In step (a), when the isolated pipe is extracted, the pressure of the node at the end of the isolated pipe is set to “0”, the pressure of the node connected to the isolated pipe is obtained, and the isolated pipe is isolated according to the positive / negative of the pressure. It is also possible to visually display the pipes.

[0016]

[Examples] Japanese Patent Application No. 6-20300
8 together, the embodiment of the present invention will be described in detail with reference to the drawings, taking into consideration to deepen the understanding of the entire system. FIG. 1 is a flowchart of a method for calculating a pipe network according to an embodiment of the present invention. As a target pipe network, for example, a pipe network as shown in FIG. 5 is taken.

A description will be given below with reference to these figures. Step 100 First, step 100 will be described.

In step 100, the number of nodes NN indicating the number of junction points of the pipe network, the fluid demand q (NN) allocated to each node, the number of pipes PN of the pipe network, the pipe resistance RES (PN) of each pipe, Of the nodes at both ends of the pipe, one node number NPPN1 (PN) arbitrarily determined, the other node number NPPN2 (PN), the fluid supply source node number NF, and the fluid supply source node number NFIX (NF) are set as known numbers. The storage unit to be stored is defined.

NN, PN, N in parentheses in the alphanumeric characters
F indicates that there are NN pieces, PN pieces, and NF pieces of respective data. When there is a demand for fluid in the middle of the pipe, the fluid demand q (N
N) and are added. The pipe resistance R
ES (PN) is fluid pressure P, flow rate Q, specific gravity P,
Assuming that the pipe inner diameter is D, the extension is L, the parameters are M and N, and the functions are F and G, P ^m = F (Q) · G (P, L,
D) When expressed by Q ^N , it is a value of G (P, L, D) that can be calculated regardless of the flow rate Q.

In the case of the pipe network shown in FIG. 5, NN = 15, PN
= 16, NF = 2, and a node number is assigned to each node. Furthermore, each pipe is also numbered, and the resistance of each pipe is calculated in advance. In the case shown in FIG. 5, node numbers from "1" to "15" are assigned,
A pipe number from "1" to "16" is assigned to each pipe.

FIG. 2 is a diagram showing one node number NPPN1 (PN) and the other node number NPPN2 (PN) which are arbitrarily defined among the nodes at both ends of the pipe.

In FIG. 2, A indicates a pipe number, there are nodes at both ends of the pipe, and the number in the node is the node number. C indicates the resistance of the pipe. Here, if one node is “9” for the pipe 13, the other node is “2”, and NPPN1 (13) = 9, NP
PN2 (13) = 2.

In FIG. 5, since node "1" and node "4" are fluid supply source nodes, the number of fluid supply source nodes is NF = 2, and fluid supply source node number NFLX.
(1) = 1 and NFLX (2) = 4.

FIG. 3 defines the demand quantity, and shows that the demand quantity of the node "2" is q (2). This demand amount will be described with reference to FIG. Step 200 Next, step 200 will be described. Step 20
0 defines a storage unit for storing unknowns.

That is, in step 200, the upstream pipe number IPIPE (N
N), the total resistance SD (NN) from the fluid source, the node order LEVEL (NN) from the fluid source, and the upstream node number NPPN1 (N) given to the pipe.
P), the downstream node number NPPN2 (NP) and the initial flow rate QN (NP), and the storage unit for storing the values of the loop number LN and the loop source pipe LOOP (LN) are defined.

As shown in FIG. 4, the upstream pipe IPIP
E shows the pipe of the route with the least resistance from the fluid source. For example, pipe numbers "1", "2", "1"
When the pipe “1” is selected from the pipes in “3” and is set as the upstream pipe, IPIPE (2) = 1.

The total resistance SD (NN) from the fluid source is
For each node, it is the smallest of the sum of the pipe resistances RES of the plurality of paths from the fluid supply source to this node.

Node order LEVEL from fluid source
(NN) indicates how downstream each node is from the fluid supply source.

In the upstream node number NPPN1 (NP) and the downstream node number NPPN2 (NP), the node having the smaller total resistance SD from the fluid supply source is connected to the upstream pipe NPP.
N1 and the larger one is the downstream pipe NPPN2.
For example, in FIG. 2, the node "9" is replaced by the node "5".
If SD is small, NPPN1 (13) = 9 and N
PPN2 (13) = 2.

The initial flow rate QN (NP) is an initial value of the flow rate when the pipe network calculation is performed.

The loop number LN indicates the number of loops defined by this embodiment.

The loop source pipe LOOP (LN) is a pipe number determined by this embodiment, and one loop source pipe exists in each of the above-mentioned loops. Step 250 Next, step 250 will be described. Step 25
0 is a process for extracting and eliminating isolated pipes.

The process of extracting the isolated pipe in the process of step 250 will be briefly described with reference to FIGS. 10, 11 and 12. FIG. 10 is a flowchart when the processing of step 250 is processed by computer.

First, node flags N _F (NN) are provided for all the nodes from “1” to “NN”, and pipes are provided for all the pipes from “1” to “PN”. A flag P _F (PN) is provided. 11 is an explanatory diagram showing the change of the pipe flag P _F in step 250 for the example of FIG. 5, and FIG. 12 is an explanatory diagram showing the change of the node flag N _F in step 250 for the example of FIG. In the example of FIG. 5, 15 node flags N _F (NN) and 16 pipe flags P _F (PN) are used.
However, a relational table as shown in "Pipe number", "Node number at both ends of pipe" in Fig. 11, and "Node number" and "Pipe number connected to node" in Fig. 12 is provided. A state in which the pipe flag P _F of FIG. 11 and the node flag N _{F of} FIG. 12 are created for each process of step 250 will be described.

As an initial value, all node flags N _F
(NN) to "0" and all pipe flags P _F (P
N) is set to "0" (step 251). In the example of FIG. 5, the stage of step 251 includes the “initial value” column of the “pipe flag P _F ” in FIG. 11 and the “node flag N in FIG.
_It is shown in the "Initial value" column of " _F ".

Next, the node flag N of the water supply source node
_F (NN) is set to "2" (step 252). That is, in the example of FIG. 5, node “1” and node “4”
Is a fluid supply source node, the node flag N
_F (1) and node flag N _F (4) are set to “2”. In the example of FIG. 5, the result is shown in the “water supply source” column of the “node flag N _F ” of FIG. 12.

Next, the pipe flag P _F of the pipe connected to the node whose node flag N _F has changed from “0” to “2” is set to “2” (step 253). In the example of FIG. 5, as shown in the “first time” column of the “pipe flag P _F ” in FIG. 11, the pipes “1”, “4”, “12”, and , The pipe flags P _F of the pipes “3” and “6” connected to the node “4” are set to “2”.

For the pipe whose pipe flag P _F has changed from "0" to "2" in step 253, the node number NPPN1 (PN) whose node flag N _F is "2"
Find the other node number NPPN2 (PN) for, if the node flag N _F of the node is "0" to "2" (step 254). In the example of FIG.
For the pipes "1", "3", "4", "6", "12", the node numbers NPPN1 (1) = 1, NPPN1
(3) = 4, ..., the other node number NPPN2
(1) = 2, NPPN2 (3) = 3, ... Are searched, and the result as shown in the “first time” column of the “node flag N _F ” in FIG. 12 is obtained.

Next, at step 254, the node flag N _F
It is determined whether there is a node whose value changes from "0" to "2" (step 255). If it is determined in step 255 that there is a node in which the node flag N _F changes from “0” to “2”, the processing from step 253 is repeated. In the example of FIG. 5, the second processing is performed, and in step 253, the “pipe flag P _F ” of FIG. 11 is executed.
As shown in the "Second" column of the
At 54, the result as shown in the “second time” column of the “node flag N _F ” in FIG. 12 is obtained. Second step 255
At the stage, it is determined that there is a node whose node flag N _F changes from “0” to “2”, and therefore the processing from step 253 to step 254 of the third time is performed, and the “pipe flag P of FIG. _The result as shown in the "3rd" column of " _F " is obtained (step 253). In the third step 254, the node numbers NPPN1 (8) = 6 and NPPN2 (8) = 7 at both ends of the pipes “8” and “10” whose pipe flag P _F has changed from “0” to “2”. , NPP
Since N1 (10) = 6 and NPPN2 (10) = 10 are all already “2”, the node flag N _F is compared with the “second time” of the “node flag N _F ” in FIG. And there is no change at all.

When it is determined in step 255 that there is no node whose node flag N _F changes from "0" to "2", all the processes from step 251 to step 255 are completed, and at this stage, the pipe Flag P _F is “0”
Is determined to be an isolated pipe. In the example of FIG. 5, it can be seen that the pipes “14”, “15”, and “16” are isolated pipes, as shown in the “third” column of the “pipe flag P _F ” in FIG. 11. By the above processing, the isolated pipes “14”, “15” and “16” existing in FIG. 5 are extracted and excluded from the calculation targets of the pipe network. Therefore, as shown in FIG. 6, the pipe network calculation is performed for the pipes "1" to "13" that do not include isolated pipes. At this time, for example, the number of nodes NN is “10” and the number of pipes PN is “13”.

As described above, as a result of step 250, in the following processing, isolated pipes can be excluded and pipe network calculation can be performed for a pipe network that does not include an isolated pipe in the set area.

Step 300 Next, step 300 will be described. Step 30
0 calculates the total resistance SD from the source of each node,
The upstream node number NPPN1, the downstream node number NPPN2 of each pipe, and the upstream pipe number IP of each node
This is a process of determining the IPE.

The process of step 300 will be briefly described with reference to FIGS. 13 and 14. FIG. 13A shows a change in the node relation data, and FIG. 13B shows a change in the pipe relation data.

In FIG. 13A, an upstream pipe number IPIPE is provisionally set for each node. in this case,
Since the nodes “1” and “4” are supply source nodes, the upstream pipe number is not set. Then, as an initial value of the total resistance SD for each node except the supply source node,
For example, "100" is given. Next, the total upstream resistance SD is calculated based on the upstream pipe number of each node and the resistance of each pipe.

For example, in FIG. 6, the upstream pipe of the node "5" is "4", the upstream pipe of the node "6" is "7", and the resistors RES of the pipes "4" and "7" are RES.
Are "20" and "15", respectively, so the upstream total resistance SD of the node "6" is "35 (= 20 + 15)".
Becomes In this way, the upstream total resistance SD of the nodes other than each supply source node is calculated.

Next, as shown in FIG.
Considering the node set in (a) and its upstream pipe number, the upstream node number NPPN1 and the downstream node number NPPN2 are determined for each pipe. For example, in FIG.
From (a), the upstream pipe of the node "6" is pipe "7".
Therefore, in FIG. 13B, the upstream node number of the pipe “7” is “5” and the downstream node number is “6”. If the upstream node number and the downstream node number cannot be determined due to the relationship between the node number and the upstream pipe number in FIG. 13A, the upstream node number and the downstream node number are provisionally determined. In this way, the upstream node number NPPN1 and the downstream node number NPPN2 are calculated for each pipe.

Next, the upstream total resistances SD (1) and SD (2) of the upstream node and the downstream node calculated in FIG. 13 (b) for each pipe are shown in FIG. 13 (a). Obtained and compared, if SD (1) ≦ SD (2), the upstream and downstream node numbers are not changed, and S
If D (1)> SD (2), the node numbers on the upstream side and the downstream side are exchanged. For example, the pipe number is "2"
Since the upstream side node and the downstream side node of the pipe of are the node "2" and the node "3", respectively.
From (a), the upstream total resistance SD (1) = 1 of the node “2” = 1
0, the upstream total resistance SD (2) = 21 of the node “3” is read, and SD (1) ≦ SD (2), so that for the pipe “2”, the upstream node number and the downstream node number are exchanged. Does not.

On the other hand, the upstream node and the downstream node of the pipe with the pipe number "3" are the node "3" and the node "4", respectively. The upstream total resistance SD (1) = 21 and the upstream total resistance SD (2) = 0 of the node “4” are read, and SD
Since (1)> SD (2), the upstream node number and the downstream node number are exchanged for the pipe "3".

As shown in FIG. 13B, the upstream side node and the downstream side node are exchanged for the pipes having the pipe numbers "3", "9" and "11".

Corresponding to this replacement, as shown in FIG. 13 (a), the total resistance SD from the supply source node is calculated for each node. The upstream pipe of the small path is the upstream pipe, and the total resistance of the path is the upstream total resistance. FIG.
In the case of (a), the upstream pipe numbers of the nodes "3", "7", and "10" are changed to "3", "9", and "11".

Next, as shown in FIG. 14, the upstream pipe number IPI for each node is calculated using the node relation data and the pipe relation data in the same procedure as described above.
Upstream node number NPPN for PE and each pipe
1 and the downstream node number NPPN2.

FIG. 18 is a flow chart when the processing shown in step 300 is processed by computer.

The value of the total resistance SD from the fluid supply source is temporarily set to a large numerical value that cannot actually occur, for example, "100" (step 301). Next, "0.0" is set to the total resistance SD of the fluid supply source node (step 30).
2). As a result, column 1 of the total resistance on the upstream side in FIG.
This means that the leftmost numerical value of 03 has been set. Next, the switches ISW1 and ISW2 are set to "0" (step 30
3) The value of the pararouter I for picking up the pipes from "1" in ascending order is set (step 304). The total resistance SD from the fluid supply sources of the upstream node and the downstream node of the I-th parameter is compared, and the upstream node NPP is compared.
A node number having a small SD value and a node number having a large SD value are set in N1 and the downstream node NPPN2, respectively (steps 305, 306, 307).

In the pipe number "1" shown in the pipe number column of FIG. 13B, NPPN1 (1) and NPPN
Since the node numbers of 2 (2) are "1" and "2", respectively, and the total resistances SD (1) and SD (2) from the supply source of each node are "0.0" and "100", respectively, SD (1)
Since ≦ SD (2) is satisfied, NPPN1 (1) and NPPN2 (2) are not replaced in step 306, and the values are “1” and “2”, respectively.

When the total resistance SD from the fluid supply source of the downstream node number is "100", the value of ISW is set to 1 (steps 308 and 309). For example, S of the node number “2” shown in the column 103 of the total resistance on the upstream side of FIG.
Since D (2) is “100”, in this case ISW
Value becomes "1". A value SDNEW obtained by adding the total resistance SD from the fluid supply source of the upstream node NPPN1 and the resistance RES of the pipe to the total resistance SD from the fluid supply source of the downstream node NPPN2, and the value of ISW2 when SDNEW is small. Is set to "1" and the value of SDNEW is set to the total resistance SD from the fluid supply source of the downstream side node NPPN2.
And the pipe number I is determined as the value of the upstream pipe IPIPE of the downstream node NPPN2. For example, in this step, the pipe number in FIG.
In this case, the downstream node number NPPN2 (2) is “2”, and the total resistance SD (2) from the fluid supply source with the node number “2” in FIG. 13A is changed from “100” to “10”. And the upstream pipe IPIPE (2) becomes 1.

Next, by adding 1 to the value of the parameter I and setting it to I, preparation for data processing is performed for the next pipe. By comparing the value of parameter I with the number of pipes NN in step 313, step 314
Then, when data processing is not performed for all the pipes, the process proceeds to step 305 to perform data processing for the next pipe, but when data processing is performed for all the pipes, the process proceeds to step 315 and both ISW1 and ISW2 are set to "0". To see if When ISW1 is "1", the total resistance SD to the fluid supply source in the node is "100", that is, the data processing is not completed because the temporary value remains, and ISW2 is "1". At this time, since the value of the total resistance SD from the fluid supply source has changed, it is not known whether or not the data processing is completed. Therefore, the process proceeds to step 303, and the data processing is indicated again.

As described above, as shown in FIG. 13B, the upstream node number NPPN1 and the downstream node number NPPN2 are exchanged with the pipe numbers "3", "9", and "11". Node numbers “2”, “3”,
The total resistance SD from the upstream IPIPE and the fluid supply source of "5" to "10" is set or changed. Particularly, in the node number "3" of FIG. 13A, the data once set or changed in the data processing process of the pipe number "2" is changed again in the data processing process of the pipe number "3". Similar re-modifications are also found at node numbers "7" and "10". In the first data processing process (FIG. 13), step 309 and step 31
Since both ISW1 and ISW2 have been changed to 1 in 2, the process proceeds to the second data processing step (FIG. 14). In the second data processing process, the upstream node number NPPN1 and the downstream node number NPPN2 change in the processing process of the pipe number "8", and steps 310 to 31 are performed.
By the data processing of 3, the upstream pipe number IPIPE of the node number “6” changes from “7” to “8”, and the total resistance SD from the fluid supply source also changes from “35” to “23”.

In FIG. 6, the upstream pipe IPIPE of the node "6" changes from the pipe number "7" to "8". In the pipe number “10” in FIG. 14, the upstream node number NPPN1 and the downstream node number NPPN2 are exchanged, but in steps 310 and 311 S
Since DNEW ≧ SD (NH), the node relation data is not moved. In the second data processing process, step 31
I pass through 2 only once, so ISW2 = 1
Has become. Therefore, the third data processing process is started according to the determination in step 315. In the third data processing process, there is no change in data and the process of step 312 is not performed, so both ISW1 and ISW2 become "0", so the data processing process ends in step 316.

As described above, according to step 300,
Even if the node number and the pipe number are freely set in the initial stage, the upstream node number NPPN1 and the downstream node number N
It can be seen that PPN2 and upstream pipe number are determined. Step 400 Next, step 400 will be described. Step 400 is
This is a process of calculating the number of loops LN and assigning a pipe that is not designated as an upstream pipe number as a loop source pipe.

FIG. 15 is a diagram showing the processing at this time. FIG. 15A shows node relation data,
The upstream pipe number IPIPE calculated for each node in step 300, the total upstream resistance SD, and the node rank level from the supply source are shown, and "99" is given as an initial value of this node rank. Figure 15 (b)
Indicates pipe-related data, and indicates a pipe resistance RES, an upstream node number NPPN1 and a downstream node number NPPN2 for each pipe. FIG. 15C shows the relationship between the loop number and the loop source pipe number LOOP.

The process in step 400 is as follows when it is simplified. The upstream pipe number I in FIG.
The pipe specified in PIPE (here, "1",
"3", "4", "8", "9", "6", "12",
The sign of the pipe resistance RES of “11”) is made negative as shown in FIG. 15B, and the pipes (“2” and “5” here, which are not specified in the upstream pipe number IPIPE of FIG. 15A). , "7", "10", "13") are loop source pipes, and the loop source pipe number LOOP is associated with the loop number as shown in FIG.

FIG. 19 is a flow chart when the processing of step 400 is processed by computer.

The parameter N is set to "1", and preparations are made for data processing in ascending order of the node numbers "1" to "NN" (step 401). Next, the node whose total resistance SD from the fluid supply source is “0” is not processed,
For a node for which the total resistance SD is not "0", the resistance value RES of the upstream pipe of the node is made negative and LEVEL of the node is set to "99" (steps 402 and 403). The numerical value "99" means a large numerical value that is not possible in practice, and need not be "99". Step 5 for the meaning of LEVEL
00 will be described later. The provisional numerical value “99” is set in LEVEL in step 400 because the preprocessing of the next step 500 is executed in step 400 together with other processing. Although the pipe resistance is replaced with a negative number in step 403, the pipe number specified for the upstream pipe is marked so that it can be identified. Instead of using a negative number, a flag is set for the total number of pipes PN. You may prepare.

Next, "1" is added to the parameter N to prepare for the data processing of the next node number (step 40).
4). The process returns to step 402 until the data processing of all the nodes is performed in step 405, and when the data processing of all the nodes is completed, the process proceeds to step 406. As a result of the data processing in steps 401 to 405, FIG.
LEVEL in (a) and pipe resistance R in FIG. 15 (b)
Data changes like ES. At steps 406 to 410, the values of the pipe resistance RES of the pipe numbers "1" to "NP" are sequentially inspected, the number of positive pipes is counted, and the pipe number is set to the loop source pipe number LOOP.

Here, the loop source pipe will be described.
When performing the pipe network calculation, for example, when there is only one fluid supply source and no loop, the flow rate of the fluid flowing through each pipe is determined by the addition of only the fluid demand of the node, not by the pressure. However, if there are multiple sources or there are loops, it is necessary to determine the flow rate at several points. The loop source pipe is a pipe at some of these places. The flow rate corresponds to the unknown flow rate to be obtained, and is the main existence in the calculation of the pipe network flow rate. The number of unknown flow rates is NF + LL−1, where NF is the fluid supply source location and LL is the number of loops. Loop source pipe number LN defined in step 200
Means the number of unknown flow rates. Step 406
Through the data processing of step 410, the loop source pipe number LOO is associated with the loop number as shown in FIG.
P is required. Step 500 Next, step 500 will be described. Step 50
0 is a process of determining the node rank LEVEL from the supply source using the upstream node number NPPN1 and the downstream node number NPPN2 of each pipe.

The process of step 500 will be briefly described with reference to FIG.

This process is a process for calculating the node rank indicating the rank of each node from the supply source node in the pipe network in which the loop source pipe determined in step 400 is removed.

As shown in FIG. 7, in the target network, the pipe designated as the loop source pipe in step 400 is indicated by a broken line. The numbers in the triangles in the figure indicate the order of nodes from the supply source. Source node "1",
The node order of "4" is "1". The node order of the nodes “2”, “5”, and “9” on the downstream side of the supply source node “1” is “2”.

The node order of the downstream nodes "3" and "8" of the supply source node "4" is "2". Furthermore, the node order of the nodes “7” and “10” on the downstream side of the node “8” is “3”. Further, the node order of the node "6" on the downstream side of the node "7" is "4".

In this way, the node rank of each node is determined.

FIGS. 16 and 17 show data changes of the node relation data and the pipe relation data in the processing of step 500.

In FIG. 16A, the node "1",
Since "4" is a source node, the node order LEVEL from the source is changed from "99" to "1".

Next, the nodes "2", "3", "5",
Since "8" and "9" are on the downstream side of the nodes "1" and "4", the node order LEVEL can be changed from "99" to "2".

Similarly in FIG. 17, the node rank from the supply source is calculated for each node.

FIG. 20 is a flow chart when the processing of step 500 is processed by computer.

At steps 501 to 504, the node order LEVEL of the fluid supply source node is set to "1". In step 505, the pipe number is sequentially from "1" to "N".
The parameter I which is designated in ascending order up to "P" is set to "1" and the ISW is set to "0". Step 51
When the ISW is changed to "1" in 3, it is possible that the data processing has not been completed, so the processing returns to step 505 and data processing is performed again. In steps 506 to 511, the pipe resistance RES
Since the loop source pipe exceeds "0", it jumps to step 511 without performing any data processing. When the pipe resistance RES is negative, the node rank LEV from the fluid supply source of the downstream node NPPN2 (2) of the pipe number I
EL is "99" and the upstream node NPP of the pipe
Only when the node rank LEVEL from the fluid source of N1 (2) is not "99"
While changing the value of W to "1", the pipe resistance RES
Is negative, the value is returned to positive, and the value obtained by adding "1" to the node rank LEVEL from the fluid supply source of the upstream node NPPN1 (2) is added to the downstream node NPP.
Set to node order LEVEL from N2 (1) fluid source. This is repeated by the number NP of pipes in steps 511 to 512. In step 513,
If the data has changed in step 510, ie, I
When SW = 1, there is a possibility that the node order LEVEL from the fluid supply source of all the nodes has not been determined, so the process returns to step 505 and the calculation is repeated.

As shown in FIG. 16 described above, except that the node level LEVEL from the supply source of the node number "6" is "99" in the first data processing, the supply from the supply sources of other nodes is changed. The node order is set. As shown in FIG. 17, in the second day processing, the node rank LEVEL from the supply source of the node number “6” is set to “4” in step 510 and the ISW is also changed to “1”. It Therefore, the process proceeds to the third data processing by step 513. Since the step 510 is not performed in the third data processing, the ISW is held at "0", and the judgment of the step 513 leads to the end of the step 514. Step 600 Next, step 600 will be described. Step 60
0 is a process of selecting a loop with the loop source pipe LOOP as the origin and the upstream node number NPPN1, the downstream node number NPPN2, and the node rank level as indexes.

This process will be briefly described with reference to FIG.

By the processing of step 400, the loop source pipe is calculated corresponding to the loop number. In this case, 5 loops have been calculated. 8 (a)-
(E) shows the five loops. For example, looking at the loop source pipe "2", the node ranks of the nodes "2" and "3" at both ends are "2".
Therefore, the nodes “1” and “4” connected to the node “2” and the node “3” and having the node order “1” are calculated. Since the nodes “1” and “4” are supply source nodes, the loop is determined to be determined at this time. That is,
As shown in FIG. 8A, nodes “1”, “2”,
A loop connecting "3" and "4" is selected.
In the present application, a loop includes not only a closed curve but also a straight line.

As for the loop source pipe “5”, the nodes “3” and “7” at both ends of the loop source pipe are “2” and “3”, respectively. The node "8" and the node "7" are connected to each other, and the node "3" and the node "8" to the node "4" both having the node order "2" are connected (FIG. 8B).

As for the loop source pipe “7”, since the node ranks of the nodes “5” and “6” at both ends thereof are “2” and “4”, respectively, the node “6” having a large node rank is tracked. Then, the nodes “7” and “8” on the upstream side of the node “6” are connected to the node “6”, and the node “5” and the node “8” both having the node rank “2” are connected to the node “1”. ] And the node “4” (FIG. 8C).

Looking at the loop source pipe 10, the nodes "6" and "10" at both ends of the loop source pipe are "4" and "3", respectively. The node "7" is combined with the node "6", and the nodes "7" and "10" upstream of the node "7" and the node "10" both having the node rank "3" are connected to the nodes "7" and "10". Combine. In this case, since the node “7” and the node “10” are connected to the same node “8”, the loop is determined here (FIG. 8 (d)).

The loop source pipe 13 is similarly processed (FIG. 8 (e)). In this way, five loops are selected.

FIG. 21 is a flow chart when the processing of step 600 is processed by computer.

In step 601, the storage units of the downstream loop RL and the upstream loop RH are defined. For the loop source pipe number LOOP, "1" is set to the parameter I for processing data in ascending order from "1" to "LN". Steps 603 to 605 determine the node rank levels LL and LH from the fluid supply sources of the nodes downstream and upstream of the loop source pipe LOOP (I). Through steps 605 to 609, the upstream pipe IPIPE, the upstream node NPPN1 of the pipe, and the node level LEVEL from the fluid supply source of the node are obtained from the node having the higher node order. In this way, the nodes connected in order from the loop on the node side having a higher node order LEVEL from the fluid supply source are traced, and the node order LEVEL from the fluid supply source on the other loop side is matched on the way. From that time, we will keep track of the nodes that contact each other. When the downstream loop and the upstream loop reach the source node in step 610, the loop is completed, and in step 612, the loop processing for the next loop source pipe number is prepared.
This case is an example of FIGS. 8A, 8B, 8C, 8D, and 8E. In step 607, when the node number of the downstream side loop and the number of the upstream side node coincide with each other on the way, it means that the loop is completed there. Therefore, preparations are made to enter the loop processing of the next loop source pipe number. This case is shown in FIG. Step 700 Next, step 700 will be described. Step 70
0 is a process in which the flow rate of the loop source pipe is set to "0" to give an initial value with the pipe network as a tree structure.

The process of step 700 will be briefly described with reference to FIG. In FIG. 9, the flow rate of the loop source pipe indicated by the broken line is “0”. The demand amount q (NN) is given to each node.

In this case, ignoring the loop source pipe, the nodes "2", "3", "5", "6", "9",
"10" is treated as an end node, and the flow rate of the pipe connected to it is the demand amount of that node. For example, regarding the node "7", the flow rate of the downstream pipe is q (6), and the demand amount of the node "7" is q.
Since it is (7), the flow rate of the pipe on the upstream side of the node “7” is q (6) + q (7). Similarly, the flow rate on the upstream side of the node “8” is q (6) + q (7) + q (8).
It becomes + q (10). In this way, the initial flow rate of the pipe is set.

FIG. 22 is a flow chart when the processing of step 700 is processed by computer.

A storage unit QN (NP) for storing the value of the pipe initial flow rate in step 701 and a flag storage unit QN (NP) for not double counting the fluid demand of each node.
To prepare. Instead of preparing the flag storage unit QN (NP), there is also a method of saving the storage capacity by making the numerical value of the demand amount of the node once added negative. Step 7
When the flow rate is added from the upstream node of the loop source pipe at 02, FLAG is set to "0", and when the flow rate is added from the downstream node, FLAG prepared for changing FLAG to "1" is set to "0". Set to. As a result of steps 704 to 713, the loop source pipe LO
From the upstream node NPPN1 of the OP, the upstream pipe IPIPE of that node in sequence, and then the upstream node N of the pipe
Like PPN1, while tracing, the fluid demand amount of each node is not double counted, and the fluid demand amount is added up to the fluid supply source node. Means for avoiding double counting are steps 708 and 70.
By adding the fluid demand amount of each node once by 9 and changing the flag QF to 1 at the same time, the fluid demand amount of that node is not added thereafter. When the addition of the flow rate from the upstream node of the loop source pipe is completed, step 7
Since the FLAG is “1” in 12, the process returns to step 704, and in step 705, the FLAG is changed to “0” and the downstream node number of the loop source pipe is set in the node K of the flow rate addition. The procedure of the flow rate addition in the subsequent steps 707 to 711 is the same as when the upstream node number of the loop source pipe is set. Step 7
13 and step 714, the flow rate is added repeatedly by the number of loop source pipes. Step 800 Next, step 800 will be described. Step 80
At 0, the flow rate of the loop source pipe is sequentially changed so as to approach the solution by the Hardy-Cross method or the like, and it is repeated until the difference from the previous calculated value becomes a certain threshold value or less.

That is, the initial flow rate shown in FIG. 9 is given to the five loops shown in FIG. 8, and the flow rate of the loop source pipe is changed by the Hardy-Cross method or the like so as to gradually approach the solution. Pipe network calculation is repeated until the difference from the calculated value falls below a certain threshold. Regarding this calculation, for example, Maeage “Design method and calculation example of water supply distribution pipeline”
A known document such as the above is applied, and the details thereof are omitted.
Thus, according to this embodiment, it is possible to automatically set the loop necessary for calculating the flow rate of a pipe network having an arbitrary number of fluid supply source nodes, demand nodes and pipes without excluding isolated pipes in advance. Since it can be performed automatically, labor saving can be achieved.

Further, in setting the initial flow rate required for calculating the flow rate of the pipe network and having a large effect on the calculation time, the pipe flow rate which seems to have the smallest flow rate in the loop is set to 0, that is, Since the initial flow rate is set as the tree structure, the calculation time can be shortened.

The present invention can be variously modified within the scope of its technical idea. For example, step 10
In the case of 0, the digital mapping data may be processed and edited to prepare, or if the number of nodes is about 50, it may be created by hand. Also, step 6
It is also possible to incorporate 00 into step 800 and set the flow rate of the loop source pipe by the Hardy-Cross method or the like while searching for a loop.

In step 250, when the isolated pipe is extracted, as shown in FIG. 23, if the pressure of the nodes (eg, node 11 and node 13) at the end of the isolated pipe is set to "0", the node 12, Since gas or the like is consumed in Nos. 14 and 15, the pressure thereof becomes negative, and the pressure of the node (for example, node 2, node 3 etc.) of the pipe which is not isolated is positive. Only the color can be displayed by changing the color.

Although step 250 is performed centering on the node, it may be processed centering on the pipe.

As explained in step 400, since the flow rate of the loop source pipe corresponds to the unknown quantity in the pipe network flow rate calculation, the closer the value is to the flow rate obtained by the process of step 800, the shorter the calculation time. . Therefore, it is necessary to take out the pipe with the smallest flow rate from the loop. Therefore, as shown in step 700, the total resistance SD from the fluid supply source of each node is
Is calculated and it is highly possible that the pipe attached to the node with the highest resistance will be the pipe with the minimum flow rate, and the flow rate of the pipe is set to “0”.

Step 700 is an independent invention,
The tree structure is assumed assuming that the pipe with the lowest flow rate is cut out from the loop. For example, the engineer may select the loop source pipe in consideration of the demand, the pipe extension, the diameter, the pressure of the fluid supply source, etc., and set the flow rate of the loop source pipe to "0" to set the initial flow rate. .
Further, the present invention can be applied to calculation of water pipe network, calculation of gas pipe network, and calculation of electric wiring network current. In the electric wire network current calculation, the pipe may be considered as the wiring and the flow rate may be considered as the current.

[0097]

As described above in detail, according to the present invention, it is possible to provide a pipe network calculation method capable of automatically setting the loop of the pipe network and setting the initial value.

[Brief description of drawings]

FIG. 1 is a flowchart showing processing steps for selecting a loop, setting an initial value, and calculating a network according to an embodiment of the present invention.

FIG. 2 is a diagram showing the definition of a pipe

[Figure 3] Diagram showing the definition of demand

FIG. 4 is a diagram showing a node definition

FIG. 5 is a diagram showing a target pipe network in this embodiment.

FIG. 6 shows the pipe network after eliminating isolated pipes.

FIG. 7 is a diagram showing a loop source pipe and node order.

FIG. 8 is a diagram showing selected loops.

FIG. 9 is a diagram showing an example of calculating an initial value of a flow rate.

FIG. 10 is a flowchart in the case where the process of step 250 is processed by a computer.

FIG. 11 is a pipe flag P _F in step 250.
Explanatory diagram showing changes in

FIG. 12: Node flag N _F in step 250
Explanatory diagram showing changes in

FIG. 13 is a diagram showing the processing of step 300.

FIG. 14 is a diagram showing a process when the number of repetitions of step 300 is 2;

FIG. 15 is a diagram showing the processing of step 400.

FIG. 16 is a diagram showing the processing of step 500.

FIG. 17 is a diagram showing a process in Step 500 where the number of repetitions is two.

FIG. 18 is a flowchart when the processing of step 300 is computer-processed.

FIG. 19 is a flowchart when the processing of step 400 is processed by a computer.

FIG. 20 is a flowchart of a case where the processing of step 500 is processed by a computer.

FIG. 21 is a flowchart when the processing of step 600 is processed by a computer.

FIG. 22 is a flowchart of a case where the process of step 700 is processed by a computer.

FIG. 23 is a diagram showing a pressure source set at an end of an isolated pipe and a calculated pressure value.

[Explanation of symbols]

A ......... pipe number B ...... node C ......... pipe resistance

Continuation of front page (58) Fields surveyed (Int.Cl. ⁷ , DB name) G06F 17/50

Claims (9)

(57) [Claims] 1.At least one source node and multiple nodes
Of a network of pipes connected to each other by a pipe with resistance
Loop selection method that uses computer to select loop
And (A) A step of extracting and eliminating isolated pipes from the pipe network
When, (B) Total resistance of the route from the source node to each node
The upstream node and downstream node of each pipe based on
And the step of determining the upstream pipe of each node, (C) In the above process, it is not designated as an upstream pipe
A step of using the pipe as a loop source pipe, (D) Ignoring the loop source pipe of the pipe network,
How far downstream the node is from the source node
The node order indicating whether the node is a node
Determining based on the node and the downstream node, (E) Upstream of each pipe with the loop source pipe as the origin
Based on the node order of the side nodes and the node order of the downstream nodes
The process of selecting the loop Equipped with,The step (a) includes Initialize flags for all nodes and flags for all pipes
Give a value, Of the source node and the pipe connected to the source node
Initialize the flag and the flag of the other node of the pipe
Change to a value other than the value, Of the other nodes, the node whose flag has changed
Against the flag of the pipe connected to the node
Changing the flag of the other node of the pipe
repetition, An isolated pie for a pipe whose flag remains at its initial value
Loop selection method for pipe networks characterized by extracting as loops
Law. 2.At least one source node and multiple nodes
Of a network of pipes connected to each other by a pipe with resistance
Initial value determination method for determining initial value using computer
And (A) A step of extracting and eliminating isolated pipes from the pipe network
When, (B) Total resistance of the route from the source node to each node
The upstream node and downstream node of each pipe based on
And the step of determining the upstream pipe of each node, (C) In the above process, an upstream pipe out of all pipes
The process of using a pipe not specified as a loop source pipe, (D) Ignoring the loop source pipe of the pipe network,
How far downstream the node is from the source node
The node order indicating whether the node is a node
Determining based on the node and the downstream node, (E) Upstream of each pipe with the loop source pipe as the origin
Based on the node order of the side nodes and the node order of the downstream nodes
The process of selecting the loop (F) Pipe network with the flow rate of the loop source pipe set to “0”
Is a tree structure, and the node
Demand of each node by determining the flow rate of the pipe supplied to the node
While determining the flow rate of the upstream pipe while adding
The process of setting the initial value, Equipped with,The step (a) includes Initialize flags for all nodes and flags for all pipes
Give a value, Of the source node and the pipe connected to the source node
Initialize the flag and the flag of the other node of the pipe
Change to a value other than the value, Of the other nodes, the node whose flag has changed
Against the flag of the pipe connected to the node
Changing the flag of the other node of the pipe
repetition, An isolated pie for a pipe whose flag remains at its initial value
Method for determining initial value of pipe network characterized by extracting as a group
Law. 3.At least one source node and multiple nodes
Of a network of pipes connected to each other by a pipe with resistance
A method of calculating using a computer, (A) A step of extracting and eliminating isolated pipes from the pipe network
When, (B) Total resistance of the route from the source node to each node
The upstream node and downstream node of each pipe based on
And the step of determining the upstream pipe of each node, (C) In the above process, an upstream pipe out of all pipes
The process of using a pipe not specified as a loop source pipe, (D) Ignoring the loop source pipe of the pipe network,
How far downstream the node is from the source node
The node order indicating whether the node is a node
Determining based on the node and the downstream node, (E) Upstream of each pipe with the loop source pipe as the origin
Based on the node order of the side nodes and the node order of the downstream nodes
The process of selecting the loop (F) Pipe network with the flow rate of the loop source pipe set to “0”
Is a tree structure, and the node
Demand of each node by determining the flow rate of the pipe supplied to the node
While determining the flow rate of the upstream pipe while adding
The process of setting the initial value, (G) Before each loop selected in step (e)
Pipe network calculation is performed using the initial values obtained in step (f).
The steps to perform, Equipped with,The step (a) includes Initialize flags for all nodes and flags for all pipes
Give a value, Of the source node and the pipe connected to the source node
Initialize the flag and the flag of the other node of the pipe
Change to a value other than the value, Of the other nodes, the node whose flag has changed
Against the flag of the pipe connected to the node
Changing the flag of the other node of the pipe
repetition, An isolated pie for a pipe whose flag remains at its initial value
A method for calculating a pipe network, which is characterized by extracting as a group. 4. The loop selection method for a pipe network according to claim 1, wherein the pipe network is a gas pipe network, a water supply pipe network, or an electric power network. 5. The method for determining an initial value of a pipe network according to claim 2, wherein the pipe network is a gas pipe network, a water supply pipe network, or an electric power network. 6. The method of calculating a pipe network according to claim 3, wherein the pipe network is a gas pipe network, a water supply pipe network, or an electric power network. 7. In the step (a), when the isolated pipe is extracted, the pressure of the node at the end of the isolated pipe is set to "0".
The method according to claim 1, wherein the pressure of a node connected to the isolated pipe is obtained, and the isolated pipe is visually distinguished and displayed according to whether the pressure is positive or negative. 8. The pressure of the node at the end of the isolated pipe is set to "0" when the isolated pipe is extracted in the step (a).
3. The method for determining an initial value of a pipe network according to claim 2, wherein the pressure of a node connected to the isolated pipe is obtained, and the isolated pipe is visually distinguished and displayed according to the positive / negative of the pressure. 9. In the step (a), when the isolated pipe is extracted, the pressure of the node at the end of the isolated pipe is set to "0".
4. The method according to claim 3, wherein the pressure of the node connected to the isolated pipe is obtained, and the isolated pipe is visually distinguished and displayed according to the positive / negative of the pressure.