Water balance calculation method, device and system for water diversion project

Abstract

The present invention discloses a method, device and system for calculating water balance of a water diversion project. The method comprises generating a directed graph of a canal system based on the relative positions of canal sections, water diversions and gates in the water diversion project; generating a generalized building model based on the building types of canal sections, water diversions and gates in the water diversion project; generating a canal system topology model based on the directed graph of the canal system and the generalized building model; traversing the canal system topology model according to the water balance theory and the changes in the inflow and outflow of the canal section by using a path search algorithm to generate the paths between each terminal water diversion and the water source; when the flow of one or more water diversions changes at the same time, according to the micro-tube flow bundle theory and the paths related to the water diversions with changed flow, a bottom-up flow superposition method is used to obtain the flow change of the relevant canal section. The present invention can quickly calculate the flow change of each canal section, improve the dispatching accuracy and the work efficiency of personnel.

Description

Technical Field

[0001] The present invention belongs to the field of water conservancy and water affairs informatization and intelligence, and specifically relates to a method, device and system for calculating water balance of a water diversion project. Background Art

[0002] As a major strategic project related to the national economy and people's livelihood, large-scale water resource allocation projects have the characteristics of long water transmission lines, large flow rates, and complex operations. Scientific and reasonable dispatching and control strategies are an important prerequisite for ensuring the safe, stable, and efficient operation of water diversion projects. Water balance is the basic principle of water resource dispatching work. Achieving water balance regulation in water transmission lines is the basis for optimizing the operation of water diversion work. Therefore, ensuring the flow balance of each channel section in the water diversion project is of great significance for the scientific and efficient operation of large-scale water resource allocation projects.

[0003] The commonly used water balance control technology is the canal control theory. By combining the two theories of automatic control and optimal scheduling, the gate opening can be adjusted in real time. However, when multiple diversion outlets change flow at the same time during the water diversion process, the flow changes in each channel section are complicated, making it difficult to achieve flow balance. In addition, for gates with poor reliability, the supply is still determined by demand during scheduling, that is, the upstream supplies as much water as the diversion outlet needs, resulting in a huge workload and easy errors, which cannot meet the real-time optimal scheduling needs of large-scale water diversion projects. Therefore, the design of a reasonable water balance control technology is an urgent problem to be solved. Summary of the invention

[0004] In view of the above problems, the present invention proposes a method, device and system for calculating water balance of water diversion projects, which can quickly calculate the flow changes of each channel section and improve the scheduling accuracy and personnel work efficiency.

[0005] In order to achieve the above technical objectives and the above technical effects, the present invention is implemented through the following technical solutions:

[0006] In a first aspect, the present invention provides a method for calculating water balance of a water diversion project, comprising:

[0007] Generate a directed graph of the canal system based on the relative positions of canal sections, diversion outlets and gates in the water diversion project;

[0008] Generate a generalized building model based on the building types of canal sections, diversion ports and gates in the water diversion project;

[0009] Generate a canal topology model based on the canal directed graph and the building generalization model;

[0010] According to the water balance theory and combined with the changes in the inflow and outflow of the canal section, the path search algorithm is used to traverse the canal topology model to generate the path between each terminal water diversion outlet and the water source;

[0011] When the flow of one or more diversion outlets changes at the same time, the flow change of the relevant channel section is obtained by using the bottom-up flow superposition method based on the micro-tube flow bundle theory and the paths related to each diversion outlet where the flow changes.

[0012] Optionally, the vertices in the canal system directed graph represent water diversion outlets or gates, and the edges in the canal system directed graph represent canal sections.

[0013] Optionally, the method for generating the canal topology model includes:

[0014] Based on the canal directed graph and the generalized building model, a universal one-dimensional hydrodynamic model is used to generate a canal topology model.

[0015] Optionally, the method for generating the path between each terminal water diversion outlet and the water source includes:

[0016] Repeat the traversal steps until all vertices in the canal topology model are traversed;

[0017] The traversal step includes: starting from the vertex representing the terminal water diversion outlet, visiting the next vertex along a branch path in the canal topology model, and marking the visited vertices until the out-degree of the vertex representing the terminal water diversion outlet is 0, and then stepping back step by step to visit the next unmarked vertex.

[0018] Optionally, the vertex representing the terminal water diversion outlet is denoted by V g , the vertex representing the water source is denoted by V s , find vertex V according to the path search algorithm g With V s The path between

[0019] For any vertex V in the canal topology model i , its adjacent vertex set is recorded as:

[0020] U_adjoin(Vi)={Vx|d_out(Vi, Vx)=1}

[0021] Where U_adjoin(Vi) represents the set of adjacent points of any vertex Vi; Vx represents the vertex adjacent to any vertex Vi; d_out represents the adjacent vertex determination condition, if the value is 1, it means that the two points are adjacent;

[0022] The visited state F_visited(Vx) of vertex Vx is recorded as:

[0023]

[0024] Optionally, the method for acquiring the flow change condition of the relevant channel section includes:

[0025] The canal system is currently in a balanced state, and the flow change at the terminal diversion outlet is 0;

[0026] Input the flow changes of the terminal diversion outlets into the canal topology model;

[0027] Traverse each water diversion in the canal topology model in turn to determine whether the flow of each water diversion has changed;

[0028] If the flow rate of one or more diversion outlets changes, the paths from each diversion outlet where the flow rate has changed to the water source are searched upwards respectively, and the flow rate change value is superimposed on each channel section. Finally, the flow superposition value of each relevant channel section is used as the flow change of each channel section.

[0029] Optionally, the expression for the flow change of each channel section is:

[0030]

[0031] In the formula, Q <i，j> Indicates channel section<i,j> The sum of the changes in flow, Q k It indicates that the diversion outlet where the flow rate has changed is superimposed on the diversion outlet flow of the channel section.

[0032] Optionally, the water balance calculation method for the water diversion project further includes:

[0033] The flow changes of relevant channel sections are displayed using a front-end web page based on HTML5 technology.

[0034] In a second aspect, the present invention provides a water balance calculation device for a water diversion project, comprising:

[0035] The canal system directed graph generation module is used to generate a canal system directed graph based on the relative positions of canal sections, diversion outlets and gates in the water diversion project;

[0036] The building generalization model generation module is used to generate the building generalization model based on the building types of the canal section, water diversion outlet and gate in the water diversion project;

[0037] A canal topology model generation module, used to generate a canal topology model based on the canal directed graph and the building generalization model;

[0038] A path generation module is used to traverse the canal topology model based on water balance theory and changes in the inflow and outflow of the canal section, and to generate paths between each terminal water diversion outlet and the water source using a path search algorithm;

[0039] The balance calculation module is used to obtain the flow change of the relevant channel section by using the bottom-up flow superposition method based on the micro-tube flow theory and the paths related to the water outlets where the flow changes when the flow of one or more water outlets changes at the same time.

[0040] In a third aspect, the present invention provides a water balance system for a water diversion project, comprising a storage medium and a processor;

[0041] The storage medium is used to store instructions;

[0042] The processor is configured to operate according to the instructions to execute the method according to any one of the first aspects.

[0043] Compared with the prior art, the present invention has the following beneficial effects:

[0044] The present invention is different from the traditional canal control theory. It adopts the graph theory technology in the network topology structure, abstracts the water diversion outlets and gates as "points", abstracts the canal sections as "edges", and expresses their logical relationship in the computer. The water balance calculation of each canal section is realized through the path search algorithm. It can not only ensure the balance of scheduling flow when the flow of each canal section changes complexly, but also greatly reduce the workload of stations along the line and improve the accuracy of scheduling.

[0045] Furthermore, the present invention also outputs the flow change of the channel section through a visualized front-end web page to provide a reference basis for scheduling decisions. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] In order to make the content of the present invention more clearly understood, the present invention is further described in detail below based on specific embodiments and in conjunction with the accompanying drawings, wherein:

[0047] Figure 1 A flow chart of water balance calculation based on graph theory according to an embodiment of the present invention;

[0048] Figure 2 A simple directed diagram of graph theory of an embodiment of the present invention;

[0049] Figure 3 A schematic diagram of the decomposition and cross-section of a micro-tube flow beam according to an embodiment of the present invention;

[0050] Figure 4 A flow chart of a canal water diversion depth-first search algorithm according to an embodiment of the present invention;

[0051] Figure 5 A flow chart of a canal water balance analysis algorithm according to an embodiment of the present invention;

[0052] FIG. 6 is a flow diagram of a channel system when the water diversion outlet changes in different situations according to an embodiment of the present invention. DETAILED DESCRIPTION

[0053] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the scope of protection of the present invention.

[0054] The application principle of the present invention is described in detail below in conjunction with the accompanying drawings.

[0055] Example 1

[0056] Graph theory, as an important branch of discrete mathematics, originated from the theory written by mathematician Euler to solve the "Seven Bridges Problem" of Königsberg. Graph theory is a mathematical model that describes the connection between things by forming a topological graph with points and edges. It is an important means to solve network search and map path planning. In the field of water conservancy, graph structure can be used to assist hydrological analysis. The specific meaning of graph theory is as follows:

[0057] As an abstract model for expressing the relationship between things, graph theory can be used to perform complex analysis and calculation after abstract modeling of things. Figure 2 As shown in Figure 1, a graph consists of two elements: vertices (V) and edges (E). Edges depend on vertices to exist. Vertices are represented by V. i Indicates, (i=0,1,2......), while using E<Vi,Vj> It represents (i=0,1,2......,j=0,1,2......). Since water diversion projects have similar point and line characteristics, it is reasonable to use graph theory to analyze them.

[0058] The micro-tube flow theory refers to: using the Lagrangian viewpoint as the path of water flow in the canal system, decomposing the water body into tiny independent particles, and assuming that the fixed particles follow an independent path from the starting point to the end point, that is, the path of the particles from the water source to the water diversion is an independent micro-tube flow, which simplifies the calculation process. For the shared canal section water volume, it is simplified to the superposition of micro-tube flow beams of independent particle movement from the water source to the water diversion. The flow change of each end water inlet will independently cause the flow change of the common canal section. Figure 3 As shown, the specific meaning is: the flow reaching V4, V5, and V6 is all from the canal section<V1,V2> Provided, but in the channel section<V1,V2> Internally, the flow bundles are considered to exist independently and the three are unrelated to each other. Figure 3 In the middle, line 1 represents the flow to V4, line 2 represents the flow to V5, and line 3 represents the flow to V6.

[0059] like Figure 4As shown in the figure, the path search algorithm represents: according to the principle of maximum path depth, the first path starts from a vertex and visits the unvisited vertices in turn until the last unvisited vertex, at which point the first path ends; if there are still unvisited vertices in the graph at this time, start from the first unvisited vertex and repeat the above steps to obtain a new path; follow this method until all vertices have been visited.

[0060] Based on the above analysis, a method for calculating water balance of a water diversion project is provided in an embodiment of the present invention. Figure 1 As shown, the following steps are included:

[0061] (1) Generate a directed graph of the canal system based on the relative positions of the canal sections, diversion outlets, and gates in the water diversion project;

[0062] (2) Generate a generalized building model based on the building types of canal sections, diversion outlets, and gates in the water diversion project;

[0063] (3) generating a canal topology model based on the canal directed graph and the building generalization model;

[0064] (4) Based on the water balance theory and the changes in the inflow and outflow of the canal section, a path search algorithm is used to traverse the canal topology model to generate paths between each terminal water diversion outlet and the water source; the path search algorithm has the characteristics of stacking and popping operations, taking up less space, and running slowly;

[0065] (5) When the flow rates of one or more diversion outlets change simultaneously, the flow rate changes of the relevant channel sections are obtained by using the bottom-up flow superposition method based on the micro-tube flow bundle theory and the paths related to each diversion outlet where the flow rate changes.

[0066] In a specific implementation of an embodiment of the present invention, the method for generating the canal directed graph is: abstracting the water diversion and gate as "points", abstracting the canal section as "edges", and establishing the canal directed graph, that is, the vertices in the canal directed graph represent the water diversion or gates, and the edges in the canal directed graph represent the canal sections.

[0067] In a specific implementation of the embodiment of the present invention, the method for generating the generalized building model is: by merging parallel canals, deleting water channels and water purification stations, etc., the hydraulic structures are appropriately simplified to establish the generalized building model.

[0068] In a specific implementation of the embodiment of the present invention, the method for generating the canal topology model includes:

[0069] Based on the canal directed graph and the generalized building model, a general one-dimensional hydrodynamic model software (such as SWMM) is used for modeling to generate a canal topology model, which can simulate, predict and manage hydraulic processes, and provide a visual interface to manage the topological relationships of canal sections and hydraulic structures.

[0070] In a specific implementation of the embodiment of the present invention, the method for generating the path between each terminal water diversion outlet and the water source includes:

[0071] Repeat the traversal steps until all vertices in the canal topology model are traversed;

[0072] The traversal step includes: starting from the vertex representing the terminal water diversion, visiting the next vertex along a branch path in the canal topology model, marking the visited vertices, until the out-degree of the vertex representing the terminal water diversion is 0, and then stepping back to visit the next unmarked vertex. Specifically:

[0073] The vertex representing the terminal water diversion outlet is denoted as V g , the vertex representing the water source is denoted by V s , find vertex V according to the path search algorithm g With V s The path between

[0074] For any vertex V in the canal topology model i , its adjacent vertex set is recorded as:

[0075] U_adjoin(Vi)={Vx|d_out(Vi, Vx)=1}

[0076] Where U_adjoin(Vi) represents the set of adjacent points of any vertex Vi; Vx represents the vertex adjacent to any vertex Vi; d_out represents the adjacent vertex determination condition, if the value is 1, it means that the two points are adjacent;

[0077] The visited state F_visited(Vx) of vertex Vx is recorded as:

[0078]

[0079] In a specific implementation of the embodiment of the present invention, the direction of the line segment in the canal system directed graph is defined according to the flow direction of the water flow, such as Figure 5 As shown, the method for obtaining the flow change of the relevant channel section includes:

[0080] The canal system is currently in a balanced state (i.e. the water volume deviation of all canal sections is set to 0), and the flow change of the terminal diversion outlet is 0;

[0081] Input the flow changes of the terminal diversion outlets into the canal topology model;

[0082] Traverse each water diversion in the canal topology model in turn to determine whether the flow of each water diversion has changed;

[0083] If the flow rate of one or more diversion outlets changes, the paths from each diversion outlet with changed flow rate to the water source are searched upwards, and the flow rate change value is superimposed on each channel section. Finally, the flow superposition value of each relevant channel section is used as the flow rate change of each channel section. Specifically, the expression of the flow rate change of each channel section is:

[0084]

[0085] In the formula, Q <i，j> Indicates channel section<i,j> The sum of the changes in flow, Q k It indicates that the diversion outlet where the flow rate has changed is superimposed on the diversion outlet flow of the channel section.

[0086] In a specific implementation of the embodiment of the present invention, the water balance calculation method of the water diversion project further includes:

[0087] The flow changes of relevant channel sections are displayed using a front-end web page based on HTML5 technology.

[0088] In the specific implementation process, the present invention adopts a design mode of front-end and back-end separation. The front-end is responsible for interactive display, and the back-end is responsible for the calculation of the flow change of the relevant channel section (i.e., water balance calculation). The front-end is developed using HTML5 (HyperTextMarkup Language 5), and the back-end uses graph theory to develop the flow change analysis program of the channel section. The front-end and back-end interact through the agreed interface. The back-end program reads the *.inp file (i.e., the channel topology model) after the modeling is completed, and can directly identify the water diversion, gates, and channel sections.

[0089] The water balance calculation method for a water diversion project in an embodiment of the present invention is described in detail below in conjunction with a specific implementation method.

[0090] As shown in Figure 6, a diagram of the flow change in a channel section of a project using the water balance calculation method for a water diversion project proposed in an embodiment of the present invention is shown. In the figure, changes in the flow demand at the water diversion outlet are represented by circles; dotted pipes indicate that the flow has not changed, and solid pipes indicate that the flow has changed. A positive value indicates an increase in flow, and a negative value indicates a decrease in flow.

[0091] Figure 6(a) is a generalized graph of a certain project area. The arrows indicate the flow direction. At this time, there is no flow adjustment plan for each diversion. Figure 6(b) shows that when the diversion N-13 changes, the flow of N-13 needs to increase by 10m 3 / s, the flow rates of N-18, N-16, N-14, N-10, and N-9 upstream need to be increased by 10m 3 / s, the flow adjustment of each channel section is shown in the figure; Figure 6(c) shows that when the diversion outlets N-13 and N-19 change, the flow of N-13 needs to increase by 10m 3 / s, N-19 traffic needs to be increased by 15m 3 / s, at this time, the flow of the upstream channel section is superimposed, and the flow of N-18, N-16, N-14, N-10, and N-9 upstream needs to increase by 25m 3 / s; Similarly, Figure 6(d) shows that on the basis of Figure 6(c), the flow demand of water diversion port N-23 is reduced by 20m 3 / s, the flow adjustment of each channel section; Figure 6(e) shows the flow demand of the water diversion outlet N-27 increased by 25m on the basis of Figure 6(d). 3 / s, the flow adjustment of each channel section; Figure 6(f) shows the flow demand of the water diversion outlet N-4 reduced by 30m on the basis of Figure 6(e). 3 / s, the flow adjustment situation of each upstream channel section.

[0092] Example 2

[0093] Based on the same inventive concept as that of Example 1, an embodiment of the present invention provides a water balance calculation device for a water diversion project, comprising:

[0094] The canal system directed graph generation module is used to generate a canal system directed graph based on the relative positions of canal sections, diversion outlets and gates in the water diversion project;

[0095] The building generalization model generation module is used to generate the building generalization model based on the building types of the canal section, water diversion outlet and gate in the water diversion project;

[0096] A canal topology model generation module, used to generate a canal topology model based on the canal directed graph and the building generalization model;

[0097] A path generation module is used to traverse the canal topology model based on water balance theory and changes in the inflow and outflow of the canal section, and to generate paths between each terminal water diversion outlet and the water source using a path search algorithm;

[0098] The balance calculation module is used to obtain the flow change of the relevant channel section by using the bottom-up flow superposition method based on the micro-tube flow theory and the paths related to the water outlets where the flow changes when the flow of one or more water outlets changes at the same time.

[0099] In a specific implementation of the embodiment of the present invention, the water balance calculation device for the water diversion project further includes:

[0100] The display module is used to display the flow changes of relevant channel sections using the front-end web page based on HTML5 technology.

[0101] In a specific implementation of an embodiment of the present invention, the method for generating the canal directed graph is: abstracting the water diversion and gate as "points", abstracting the canal section as "edges", and establishing the canal directed graph, that is, the vertices in the canal directed graph represent the water diversion or gates, and the edges in the canal directed graph represent the canal sections.

[0102] In a specific implementation of the embodiment of the present invention, the method for generating the generalized building model is: by merging parallel canals, deleting water channels and water purification stations, etc., the hydraulic structures are appropriately simplified to establish the generalized building model.

[0103] In a specific implementation of the embodiment of the present invention, the method for generating the canal topology model includes:

[0104] Based on the canal directed graph and the generalized building model, a general one-dimensional hydrodynamic model software (such as SWMM) is used to perform modeling to generate a canal topology model.

[0105] In a specific implementation of the embodiment of the present invention, the method for generating the path between each terminal water diversion outlet and the water source includes:

[0106] Repeat the traversal steps until all vertices in the canal topology model are traversed;

[0107] The traversal step includes: starting from the vertex representing the terminal water diversion, visiting the next vertex along a branch path in the canal topology model, marking the visited vertices, until the out-degree of the vertex representing the terminal water diversion is 0, and then stepping back to visit the next unmarked vertex. Specifically:

[0108] The vertex representing the terminal water diversion outlet is denoted as V g , the vertex representing the water source is denoted by V s , find vertex V according to the path search algorithm g With V s The path between

[0109] For any vertex V in the canal topology model i , its adjacent vertex set is recorded as:

[0110] U_adjoin(Vi)={Vx|d_out(Vi, Vx)=1}

[0111] Where U_adjoin(Vi) represents the set of adjacent points of any vertex Vi; Vx represents the vertex adjacent to any vertex Vi; d_out represents the adjacent vertex determination condition, if the value is 1, it means that the two points are adjacent;

[0112] The visited state F_visited(Vx) of vertex Vx is recorded as:

[0113]

[0114] In a specific implementation of the embodiment of the present invention, the method for obtaining the flow change condition of the relevant channel section includes:

[0115] The canal system is currently in a balanced state, and the flow change at the terminal diversion outlet is 0;

[0116] Input the flow changes of the terminal diversion outlets into the canal topology model;

[0117] Traverse each water diversion in the canal topology model in turn to determine whether the flow of each water diversion has changed;

[0118] If the flow rate of one or more diversion outlets changes, the paths from each diversion outlet with changed flow rate to the water source are searched upwards, and the flow rate change value is superimposed on each channel section. Finally, the flow superposition value of each relevant channel section is used as the flow rate change of each channel section. Specifically, the expression of the flow rate change of each channel section is:

[0119]

[0120] In the formula, Q <i，j> Indicates channel section<i,j> The sum of the changes in flow, Q k It indicates that the diversion outlet where the flow rate has changed is superimposed on the diversion outlet flow of the channel section.

[0121] Example 3

[0122] Based on the same inventive concept as that of Example 1, a water balance system for a water diversion project is provided in an embodiment of the present invention, comprising a storage medium and a processor;

[0123] The storage medium is used to store instructions;

[0124] The processor is configured to operate according to the instruction to execute the method according to any one of the embodiments 1

[0125] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0126] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0127] These computer program instructions may also be stored in a computer readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0128] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0129] The embodiments of the present invention are described above in conjunction with the accompanying drawings, but the present invention is not limited to the above-mentioned specific implementation methods. The above-mentioned specific implementation methods are merely illustrative and not restrictive. Under the enlightenment of the present invention, ordinary technicians in this field can also make many forms without departing from the scope of protection of the purpose of the present invention and the claims, which all fall within the protection of the present invention.

[0130] The above shows and describes the basic principles and main features of the present invention and the advantages of the present invention. It should be understood by those skilled in the art that the present invention is not limited to the above embodiments. The above embodiments and descriptions are only for explaining the principles of the present invention. Without departing from the spirit and scope of the present invention, the present invention may have various changes and improvements, which fall within the scope of the present invention to be protected. The scope of protection of the present invention is defined by the attached claims and their equivalents.

Claims

1. A method for calculating water balance of a water diversion project, characterized in that: include: Generate a directed graph of the canal system based on the relative positions of canal sections, diversion outlets and gates in the water diversion project; Generate a generalized building model based on the building types of canal sections, diversion ports and gates in the water diversion project; Generate a canal topology model based on the canal directed graph and the building generalization model; According to the water balance theory and combined with the changes in the inflow and outflow of the canal section, the path search algorithm is used to traverse the canal topology model to generate the path between each terminal water diversion outlet and the water source; When the flow rate of one or more diversion outlets changes at the same time, the flow rate change of the relevant channel section is obtained by using the bottom-up flow superposition method based on the micro-tube flow bundle theory and the paths related to each diversion outlet where the flow rate changes; The method for generating the path between each terminal water diversion outlet and the water source includes: Repeat the traversal steps until all vertices in the canal topology model are traversed; The traversal step includes: starting from the vertex representing the terminal water diversion outlet, visiting the next vertex along a branch path in the canal topology model, and marking the visited vertices until the out-degree of the vertex representing the terminal water diversion outlet is 0, and then step-by-step back to visit the next unmarked vertex; The vertex representing the terminal water diversion outlet is denoted as V g , the vertex representing the water source is denoted by V s , find vertex V according to the path search algorithm g With V s The path between For any vertex V in the canal topology model i , its adjacent vertex set is recorded as: U_adjoin(Vi)={Vx|d_out(Vi, Vx)=1} Where U_adjoin(Vi) represents the set of adjacent points of any vertex Vi; Vx represents the vertex adjacent to any vertex Vi; d_out represents the adjacent vertex determination condition, if the value is 1, it means that the two points are adjacent; The visited state F_visited(Vx) of vertex Vx is recorded as: The method for obtaining the flow change of the relevant channel section includes: The canal system is currently in a balanced state, and the flow change at the terminal diversion outlet is 0; Input the flow changes of the terminal diversion outlets into the canal topology model; Traverse each water diversion in the canal topology model in turn to determine whether the flow of each water diversion has changed; If the flow rate of one or more diversion outlets changes, the paths from each diversion outlet with changed flow rate to the water source are searched upwards, and the flow rate change value is superimposed on each channel section. Finally, the flow superposition value of each relevant channel section is used as the flow rate change of each channel section. The expression of flow change in each channel section is: In the formula, Q ＜i，j＞ represents the sum of the changes in flow in channel section <i, j>, Q k It indicates that the diversion outlet where the flow rate has changed is superimposed on the diversion outlet flow of the channel section.

2. A method for calculating water balance of a water diversion project according to claim 1, characterized in that: The vertices in the canal system directed graph represent water diversion outlets or gates, and the edges in the canal system directed graph represent canal sections.

3. The method for calculating water balance of a water diversion project according to claim 1, characterized in that: The method for generating the canal topology model comprises: Based on the canal directed graph and the generalized building model, a universal one-dimensional hydrodynamic model is used to generate a canal topology model.

4. The method for calculating water balance of a water diversion project according to claim 1, characterized in that: The water balance calculation method of the water diversion project also includes: The flow changes of relevant channel sections are displayed using a front-end web page based on HTML5 technology.

5. A water balance calculation device for a water diversion project, characterized in that: include: The canal system directed graph generation module is used to generate a canal system directed graph based on the relative positions of canal sections, diversion outlets and gates in the water diversion project; The building generalization model generation module is used to generate the building generalization model based on the building types of the canal section, water diversion outlet and gate in the water diversion project; A canal topology model generation module, used to generate a canal topology model based on the canal directed graph and the building generalization model; A path generation module is used to traverse the canal topology model based on water balance theory and changes in the inflow and outflow of the canal section, and to generate paths between each terminal water diversion outlet and the water source using a path search algorithm; The balance calculation module is used to obtain the flow change of the relevant channel section by using the bottom-up flow superposition method based on the micro-tube flow bundle theory and the paths related to each water outlet where the flow changes when the flow of one or more water outlets changes at the same time; The method for generating the path between each terminal water diversion outlet and the water source includes: Repeat the traversal steps until all vertices in the canal topology model are traversed; The traversal step includes: starting from the vertex representing the terminal water diversion outlet, visiting the next vertex along a branch path in the canal topology model, and marking the visited vertices until the out-degree of the vertex representing the terminal water diversion outlet is 0, and then step-by-step back to visit the next unmarked vertex; The vertex representing the terminal water diversion outlet is denoted as V g , the vertex representing the water source is denoted by V s , find vertex V according to the path search algorithm g With V s The path between For any vertex V in the canal topology model i , its adjacent vertex set is recorded as: U_adjoin(Vi)={Vx|d_out(Vi, Vx)=1} Where U_adjoin(Vi) represents the set of adjacent points of any vertex Vi; Vx represents the vertex adjacent to any vertex Vi; d_out represents the adjacent vertex determination condition, if the value is 1, it means that the two points are adjacent; The visited state F_visited(Vx) of vertex Vx is recorded as: The method for obtaining the flow change of the relevant channel section includes: The canal system is currently in a balanced state, and the flow change at the terminal diversion outlet is 0; Input the flow changes of the terminal diversion outlets into the canal topology model; Traverse each water diversion in the canal topology model in turn to determine whether the flow of each water diversion has changed; If the flow rate of one or more diversion outlets changes, the paths from each diversion outlet with changed flow rate to the water source are searched upwards, and the flow rate change value is superimposed on each channel section. Finally, the flow superposition value of each relevant channel section is used as the flow rate change of each channel section. The expression of flow change in each channel section is: In the formula, Q ＜i，j＞ represents the sum of the changes in flow in channel section <i, j>, Q k It indicates that the diversion outlet where the flow rate has changed is superimposed on the diversion outlet flow of the channel section.

6. A water balance system for a water diversion project, characterized in that: The said comprises a storage medium and a processor; The storage medium is used to store instructions; The processor is configured to operate according to the instructions to perform the method according to any one of claims 1-4.

Images