Method, device, equipment, medium and product for calculating river split ratio

Abstract

The invention relates to the technical field of hydrology and waterpower, and discloses a method, a device, equipment, a medium and a product for calculating a river diversion ratio, and the method comprises the steps: obtaining water and sediment, terrain, terrain data corresponding to different working conditions and flood flow of a target river reach; calculating water and sediment, topographic data and flood flow under a fixed bed to obtain first hydrological information of a river reach in a future time period; under the moving bed, calculating daily water and sediment, topographic data corresponding to each working condition scheme and flood flow, and obtaining second hydrological information of the river reach in a future time period when each flood flow occurs under different working conditions; according to the first hydrological information and the second hydrological information, the river channel section is divided into sub-sections, the flow is calculated, and a first split ratio and a second split ratio are determined according to the sub-section flow ratio; and determining a target river reach regulation and control scheme according to the first split ratio and the different second split ratio. The split ratio of a future river channel under the fixed bed and movable bed working conditions is compared.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of hydrology and hydraulics, and particularly relates to a method and device for calculating river diversion ratio, equipment, medium and product. BACKGROUND

[0002] For the actual river diversion ratio measurement, the measured flow velocity and water depth can be obtained by using the hydrological measurement method, the current flow velocity flow area is obtained by combining the river bed surface elevation, the flow of different branches on the same section of the river is calculated, and thus the diversion ratio is obtained. At present, the actual river diversion ratio measurement and calculation method is relatively mature, but for braided rivers, the river bed is constantly changing, and under such dynamic conditions, the diversion ratio is also a continuously changing quantity. For braided rivers that are still in the planning and regulation stage, when there are multiple regulation schemes to be compared and selected, the diversion ratio after the implementation of the scheme cannot be predicted by field measurement. SUMMARY

[0003] The present application provides a method, device, equipment and medium for calculating river diversion ratio to solve the problem that the diversion ratio after the implementation of the scheme cannot be predicted by field measurement.

[0004] In a first aspect, the present application provides a method for calculating river diversion ratio, comprising: obtaining long series of daily water and sediment data of a target river section, initial terrain data, modified terrain data corresponding to different working condition schemes, and flow when different frequency floods occur; Under the condition of fixed bed, the long series of daily water and sediment data, the initial terrain data and the flow when different frequency floods occur are calculated by using a two-dimensional mathematical model of water and sediment to obtain first hydrological information of the target river section in the future period; Under the condition of dynamic bed, the long series of daily water and sediment data, the modified terrain data corresponding to each working condition scheme and the flow when different frequency floods occur are calculated by using a two-dimensional mathematical model of water and sediment to obtain second hydrological information of the target river section in the future period when each frequency flood occurs after the implementation of different working condition schemes on the target river section; For each frequency flood, each section in the target river is divided into a plurality of sub-sections and the flow of each sub-section is calculated according to the first hydrological information, and the first diversion ratio of the section is determined according to the flow ratio of each sub-section in the section. For each frequency flood, each section in the target river is divided into a plurality of sub-sections and the flow of each sub-section is calculated according to the second hydrological information, and the second diversion ratio of the section is determined according to the flow ratio of each sub-section in the section when different working condition schemes are implemented on the target river section. The regulation scheme of the target river section is determined according to the first diversion ratio corresponding to each frequency flood and the second diversion ratio corresponding to different working condition schemes.

[0005] According to the historical data of the river channel and the conditions of the fixed bed and the moving bed, the influence of the fixed bed and the moving bed on the future river channel diversion is determined, and then the conditions of different working conditions are introduced, so that the influence of different working conditions on the future river channel diversion under the moving bed condition can be determined. According to the prediction results under different conditions, the diversion ratio under each condition can be calculated, which provides data support for the selection of subsequent regulation schemes.

[0006] In an optional implementation, under the moving bed condition, the two-dimensional water and sediment mathematical model is used to calculate the long series of daily water and sediment data, the modified terrain data corresponding to each working condition scheme, and the flow when different frequency floods occur, to obtain the second hydrological information of the target river section in the future period when each frequency flood occurs under the implementation of different working condition schemes, including: Under the moving bed condition, the two-dimensional water and sediment mathematical model is used to calculate the long series of daily water and sediment data and the modified terrain data corresponding to each working condition scheme, to obtain a plurality of third hydrological information of the target river section in the future period, and one working condition scheme corresponds to one set of third hydrological information. The third hydrological information includes the elevations of different positions of the target river section; The evolution terrain data of different working condition schemes is extracted according to the elevations of different positions in the third hydrological information; The evolution terrain data of different working condition schemes and the flow when different frequency floods occur are input into the two-dimensional water and sediment mathematical model to obtain the second hydrological information of the target river section in the future period when each frequency flood occurs under the implementation of different working condition schemes under the moving bed condition.

[0007] By using the two-dimensional water and sediment model, the long series of daily water and sediment data and the modified terrain data corresponding to each working condition scheme are introduced to obtain the predicted river channel terrain data in the future period as the third hydrological information, and then the flood flow data is introduced to obtain a plurality of second hydrological information under each working condition. The plurality of second hydrological information can be compared with the first hydrological information to obtain the influence of a working condition on the future river bed diversion ratio under the conditions of the fixed bed and the moving bed, and the influence of different working conditions on the future river bed diversion ratio under the moving bed condition can also be determined by comparing the plurality of second hydrological information.

[0008] In an optional implementation, for each frequency flood, each section in the target river channel is divided into a plurality of sub-sections according to the first hydrological information, and the flow of each sub-section is calculated, including: A plurality of section lines are determined according to a pre-established section line file, and each section line corresponds to an endpoint coordinate; The tree data structure algorithm is used to determine the river boundary point closest to the endpoint coordinate in the first hydrological information to obtain the section corresponding to each section line; The section is divided into a plurality of sub-sections.

[0009] A scheme for determining a section by an intersection of a section line and a river boundary is provided, which can ensure that the section can contain the flow of all tributaries of the study river reach. A scheme for manually dividing the section into sub-sections is also provided, which is used for flexible calculation of the diversion ratio.

[0010] In an alternative embodiment, the first hydrological information includes water depth, horizontal flow velocity and vertical flow velocity at different positions in the target river reach; the step of calculating the flow of each sub-section according to the first hydrological information comprises: dividing the section into a plurality of small sections; interpolating the first hydrological information to obtain water depth, horizontal flow velocity and vertical flow velocity of each small section; calculating flow velocity projection of each small section according to the normal vector, horizontal flow velocity and vertical flow velocity of each small section; calculating the flow of each small section according to the flow velocity projection, length and water depth of each small section; adding the flow of all small sections in the same sub-section to obtain the flow of each sub-section.

[0011] In the process of calculating the flow of the river, the water flow velocity in the river is not uniform, and the method of dividing the section into a plurality of small sections for separate calculation and accumulation can improve the calculation accuracy.

[0012] In an alternative embodiment, when different working condition schemes are implemented for the target river reach, for each frequency flood, each section in the target river reach is divided into a plurality of sub-sections according to the second hydrological information and the flow of each sub-section is calculated, comprising: determining a plurality of section lines according to a pre-established section line file, and end point coordinates corresponding to each section line; determining the river boundary points closest to the end point coordinates in the second hydrological information by using a tree data structure algorithm to obtain the sections corresponding to each section line; dividing the sections into a plurality of sub-sections.

[0013] A scheme for determining a section by an intersection of a section line and a river boundary is provided, which can ensure that the section can contain the flow of all tributaries of the study river reach. A scheme for manually dividing the section into sub-sections is also provided, which is used for flexible calculation of the diversion ratio.

[0014] In an alternative embodiment, the second hydrological information includes water depth, horizontal flow velocity and vertical flow velocity at different positions in the target river reach; the step of calculating the flow of each sub-section according to the second hydrological information comprises: dividing the section into a plurality of small sections; interpolating the second hydrological information to obtain water depth, horizontal flow velocity and vertical flow velocity of each small section; calculating flow velocity projection of each small section according to the normal vector, horizontal flow velocity and vertical flow velocity of each small section; The flow of each section is calculated according to the flow velocity projection, length and water depth of each section; The flow of each sub-section is obtained by adding the flow of all sections in the same sub-section.

[0015] In the river flow calculation process, the water flow velocity distribution in the river is uneven, and the calculation accuracy can be improved by dividing the section into multiple sections, respectively calculating and accumulating.

[0016] In a second aspect, the present application provides a device for calculating river diversion ratio, comprising: A data acquisition module for acquiring data required for calculation in the river to be measured; A data calculation module for calculating the river diversion ratio according to the data required for calculation in the river to be measured.

[0017] In a third aspect, the present application provides an electronic device, comprising a memory and a processor, the memory and the processor are communicatively connected, the memory stores computer instructions, and the processor executes the computer instructions to perform the method for calculating river diversion ratio of the first aspect or any of its corresponding embodiments.

[0018] In a fourth aspect, the present application provides a computer readable storage medium, which stores computer instructions for making a computer execute the method for calculating river diversion ratio of the first aspect or any of its corresponding embodiments.

[0019] In a fifth aspect, the present application provides a computer program product comprising computer instructions for making a computer execute the method for calculating river diversion ratio of the first aspect or any of its corresponding embodiments. BRIEF DESCRIPTION OF DRAWINGS

[0020] In order to more clearly illustrate the specific embodiments of the present application or the technical solutions in the prior art, the drawings needed in the specific embodiments or prior art description will be briefly introduced below. Obviously, the drawings in the following description are some embodiments of the present application, and those skilled in the art can obtain other drawings according to these drawings without creative labor.

[0021] Figure 1 is a schematic diagram of an application scenario according to an embodiment of the present application; Figure 2 is a first flowchart of a method for calculating river diversion ratio according to an embodiment of the present application; Figure 3 is a schematic diagram of a braided river calculated by the present application; Figure 4is a second flowchart of a method for calculating a riverway split ratio according to an embodiment of the present application; Figure 5 is a schematic diagram of the content contained in the third hydrological information; Figure 6 is a third flowchart of a method for calculating a riverway split ratio according to an embodiment of the present application; Figure 7 is a schematic diagram of the content contained in the cross-section line file; Figure 8 is a schematic diagram of a sub-cross-section division method; Figure 9 is a fourth flowchart of a method for calculating a riverway split ratio according to an embodiment of the present application; Figure 10 is a fifth flowchart of a method for calculating a riverway split ratio according to an embodiment of the present application; Figure 11 is a sixth flowchart of a method for calculating a riverway split ratio according to an embodiment of the present application; Figure 12 is a structural block diagram of a device for calculating a riverway split ratio according to an embodiment of the present application; Figure 13 is a hardware structure schematic diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0022] To make the objectives, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described below in connection with the drawings in the embodiments of the present application. Obviously, the described embodiments are some but not all of the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work fall within the scope of protection of the present application.

[0023] It can be understood that, before using the technical solutions disclosed in the embodiments of the present application, the type, use range, use scenario and the like of the personal information involved in the present application should be informed to the user and the authorization of the user should be obtained through appropriate means according to relevant laws and regulations.

[0024] The terms "first", "second" are only used for descriptive purposes and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined with "first", "second" can explicitly or implicitly include one or more of the features. In the description of the present application, the meaning of "multiple" is two or more, unless otherwise specifically limited.

[0025] As an optional application scenario of the embodiments of the present application, as shown inFigure 1 As shown, the system for calculating the river diversion ratio may include at least one terminal device and at least one server. Figure 1 The system is illustrated in the example, which includes a computer 101, a mobile terminal 102, and a server 103, and the terminal devices such as the computer 101 and the mobile terminal 102 are connected to the server 103 through a network 110.

[0026] Specifically, the terminal device can be a smartphone, tablet, laptop, PDA, desktop computer, game console, smart TV, smart wearable device, in-vehicle terminal, VR (Virtual Reality) device, AR (Augmented Reality) device, etc. Server 103 can be a standalone physical server, a server cluster, a distributed system, or a cloud server providing cloud services. Network 110 can be a wired or wireless network, examples of which include, but are not limited to, the Internet, corporate intranet, local area network, wide area network, mobile communication network, and combinations thereof.

[0027] This invention provides a method for calculating river diversion ratios. By introducing different operating conditions to calculate historical river data and topographic data, the method aims to predict future river diversion under different operating conditions and calculate the diversion ratio.

[0028] According to an embodiment of the present invention, a method for calculating the diversion ratio of a river channel is provided. It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions. Furthermore, although a logical order is shown in the flowchart, in some cases, the steps shown or described may be executed in a different order than that shown here.

[0029] This embodiment provides a method for calculating river diversion ratios, which can be used on mobile terminals such as mobile phones and tablets. Figure 2 This is a flowchart of a method for calculating the river diversion ratio according to an embodiment of the present invention, as shown below. Figure 2 As shown, the process includes the following steps: Step S201: Obtain long-term daily water and sediment data, initial topographic data, modified topographic data corresponding to different working conditions, and flow rates during floods of different frequencies for the target river section.

[0030] Since the method provided in this embodiment of the invention is used to calculate the diversion ratio of a river channel, the target river segment obtained should contain multiple tributaries. For example, the target river segment can be as follows: Figure 3 The braided river shown.

[0031] The long series of daily water and sediment data refers to the quantitative data of the water body and the sediment contained therein in the long time range, and the initial terrain data refers to the elevation data of the riverbed terrain to be calculated. Different working conditions refer to the case of changing the structure of the riverbed by artificial intervention, and the modified terrain data refers to the quantitative data of the change of the riverbed elevation data by artificial intervention under different working conditions. As an example, different working conditions can refer to the process of artificial dredging, dam building, etc., and the modified terrain data refers to the quantitative data of the change of the riverbed elevation in the corresponding dredging and dam building process. In the specific implementation process, the modified riverbed terrain under some working conditions will not change with the water and sediment scouring, and the modified terrain data is also a fixed value. The flow of different frequency floods refers to the flow size data of different floods.

[0032] Step S202, under the fixed bed condition, the long series of daily water and sediment data, the initial terrain data and the flow of different frequency floods are calculated by using the two-dimensional mathematical model of water and sediment to obtain the first hydrological information of the target river section in the future period.

[0033] The fixed bed state refers to the case where the riverbed is made of rigid materials that cannot be scoured and moved during the river simulation calculation. In this state, the water flow can scour the sediment on the riverbed, but the shape of the riverbed remains fixed. The first hydrological information refers to the prediction data obtained by using the two-dimensional mathematical model of water and sediment to predict the future river conditions under the fixed bed state according to the long series of daily water and sediment data, the initial terrain data and the flow of different frequency floods in the past period. The first hydrological information can be used as comparison information for introducing the prediction results of the moving bed state in the subsequent process.

[0034] Step S203, under the moving bed condition, the long series of daily water and sediment data, the modified terrain data corresponding to each working condition scheme and the flow of different frequency floods are calculated by using the two-dimensional mathematical model of water and sediment to obtain the second hydrological information of the target river section in the future period when each frequency flood occurs under the implementation of different working condition schemes.

[0035] The moving bed state refers to the case where the riverbed is made of materials that can be scoured and moved during the river simulation calculation. In this state, the shape of the riverbed will change constantly with the action of the water flow. The second hydrological information refers to the multiple prediction results obtained by predicting each working condition under the condition of different working conditions, combined with the long series of daily water and sediment data in the past period, the modified terrain data corresponding to each working condition scheme and the flow of different frequency floods. The second hydrological information can be compared with the first hydrological information to determine the influence of the fixed bed state and the moving bed state on the prediction results of the future diversion, and the second hydrological information can also be compared with each other to determine the influence of different working conditions on the prediction results of the future diversion under the moving bed state.

[0036] In step S204, for each frequency flood, each section in the target river is divided into multiple sub-sections according to the first hydrological information, and the flow of each sub-section is calculated, and the first flow distribution ratio of the section is determined according to the flow ratio of each sub-section in the section.

[0037] The frequency flood is used to describe the size of the flood flow, the section refers to a vertical cutting surface set by a person on the river, and the flow distribution is analyzed by measuring the flow on the cutting surface. The sub-section refers to multiple small sections divided by a person on the section, and the flood flow of each sub-section is calculated, and the flood flow of each sub-section is compared with each other to obtain the flow ratio data of each sub-section corresponding to each flow under the first hydrological information, as the first flow distribution ratio.

[0038] In step S205, for each frequency flood, each section in the target river is divided into multiple sub-sections according to the second hydrological information, and the flow of each sub-section is calculated, and the second flow distribution ratio of the section is determined according to the flow ratio of each sub-section in the section.

[0039] On the basis of the above steps, the conditions of different working conditions are introduced, the flow of each sub-section under each working condition is calculated, and the flood flow of each sub-section is compared with each other to obtain the flow ratio data of each sub-section corresponding to each flow under the second hydrological information, as the second flow distribution ratio.

[0040] In step S206, the control scheme of the target river section is determined according to the first flow distribution ratio corresponding to each frequency flood and the second flow distribution ratio corresponding to different working conditions.

[0041] The first flow distribution ratio corresponding to each frequency flood and the second flow distribution ratio under different working conditions are compared, the influence of the movable bed state on the fixed bed state of the river flow distribution can be judged, and the influence of different working conditions on the river flow distribution can be judged. According to the prediction result of the influence of different working conditions on the flow distribution, the control scheme of the target river section can be determined.

[0042] The method for calculating the river flow distribution ratio provided in this embodiment can judge the influence of the fixed bed and the movable bed on the future river flow distribution according to the historical data of the river and by introducing the conditions of the fixed bed and the movable bed, and then the influence of different working conditions on the future river flow distribution can be judged under the condition of the movable bed. According to the prediction result under different conditions, the flow distribution ratio under each condition can be calculated, which provides data support for the selection of subsequent control schemes.

[0043] In this embodiment, a method for calculating the river flow distribution ratio is provided, Figure 4is a flow chart of a method for calculating a river diversion ratio according to an embodiment of the present application, as shown in Figure 4 The flow chart comprises the following steps: In step S301, long series of daily water and sediment data of a target river section, initial terrain data, modified terrain data corresponding to different working condition schemes, and flow when floods of different frequencies occur are obtained. For details, refer to step S201 of the embodiment shown in Figure 2 The step S201 of the embodiment shown in is not repeated here.

[0044] In step S302, under the condition of a fixed bed, the long series of daily water and sediment data, the initial terrain data, and the flow when floods of different frequencies occur are calculated by using a two-dimensional mathematical model of water and sediment, to obtain first hydrological information of the target river section in a future period. For details, refer to step S202 of the embodiment shown in Figure 2 The step S202 of the embodiment shown in is not repeated here.

[0045] In step S303, under the condition of a moving bed, the long series of daily water and sediment data, the modified terrain data corresponding to each working condition scheme, and the flow when floods of different frequencies occur are calculated by using a two-dimensional mathematical model of water and sediment, to obtain second hydrological information of the target river section in a future period when floods of each frequency occur, while different working condition schemes are implemented on the target river section.

[0046] Specifically, the above step S303 comprises: In step S3031, under the condition of a moving bed, the long series of daily water and sediment data and the modified terrain data corresponding to each working condition scheme are calculated by using a two-dimensional mathematical model of water and sediment, to obtain a plurality of sets of third hydrological information of the target river section in a future period, one set of third hydrological information corresponding to one working condition scheme, and the third hydrological information containing elevations of different positions of the target river section.

[0047] The third hydrological information refers to a plurality of prediction results obtained by predicting future river diversion under different working conditions by introducing the long series of daily water and sediment data and the local terrain data file modified from the initial terrain data by using different working condition schemes. The local terrain data file includes elevation information of different positions in the modified part of the terrain. The moving bed condition can include a plurality of working conditions, and each working condition corresponds to one set of third hydrological information after modifying the initial terrain. For example, water and sediment data within 5 years or 10 years can be collected as the long series of daily water and sediment data, and the data includes daily water and sediment data within 5 years or 10 years. The third hydrological information can be Y1, Y2, and Y3, corresponding to working condition 1, working condition 2, and working condition 3, respectively.

[0048] In step S3032, the evolution terrain data of different working condition schemes is extracted according to the elevations of different positions in the third hydrological information.

[0049] The evolution topography data refers to the topography elevation data of the riverbed after the initial topography is modified and the riverbed is eroded for a period of time in the moving bed state.

[0050] The third hydrological information obtained by the two-dimensional mathematical model of water and sediment through calculation on the long series of daily water and sediment data and the modified topography data corresponding to each working condition scheme can include the horizontal coordinates, the vertical coordinates, the horizontal flow velocity, the vertical flow velocity, the water level, the elevation, the water depth, the flow velocity, and the scouring and silting thickness of multiple coordinate points in the target river section. When the third hydrological information is used to extract the evolution topography data, the evolution topography data can be obtained according to the horizontal coordinates, the vertical coordinates, and the elevation of the multiple coordinate points.

[0051] Exemplarily, the third hydrological information can be presented in the form as shown in Figure 5 In the file as shown in Figure 5 The first three lines are file descriptions, and "x", "y", "u", "v", "z", "zb", "h", "vel", and "dz" in the second line Variables represent the names of each column of data starting from the fourth line, where x represents the horizontal coordinate in the horizontal direction, y represents the vertical coordinate in the vertical direction, u represents the horizontal flow velocity, v represents the vertical flow velocity, z represents the water level, zb represents the elevation data, h represents the water depth, vel represents the flow velocity, and dz represents the scouring and silting thickness.

[0052] In step S3033, the evolution topography data of different working condition schemes and the flow rates when different frequency floods occur are input into the two-dimensional mathematical model of water and sediment to obtain the second hydrological information of the target river section in the future period when each frequency flood occurs under the moving bed condition when the different working condition schemes are implemented on the target river section.

[0053] Exemplarily, the flow rates when different frequency floods occur can be Q1, Q2, and Q3, representing the flow rate data when three different frequency floods occur. Each evolution topography data needs to be calculated in combination with Q1, Q2, and Q3, and therefore, in this example, each evolution topography data can obtain three second hydrological information. The same as the third hydrological information, the first hydrological information and the second hydrological information can include the horizontal coordinates, the vertical coordinates, the horizontal flow velocity, the vertical flow velocity, the water level, the elevation, the water depth, the flow velocity, and the scouring and silting thickness of the river branches.

[0054] In step S304, for each frequency flood, each section in the target river is divided into multiple sub-sections and the flow rates of the sub-sections are calculated according to the first hydrological information, and the first split ratio of the section is determined according to the flow rate ratio of the sub-sections in the section. For details, refer to step S204 of the embodiment as shown in Figure 2

[0055] ​Step S305, for each frequency flood, the second hydrological information is used to divide each section in the target river into a plurality of sub-sections and calculate the flow of each sub-section, and the second split ratio of the section is determined according to the flow ratio of each sub-section in the section. For details, please refer to Figure 2 Step S205 of the embodiment shown will not be repeated here.

[0056] Step S306, according to the first split ratio corresponding to each frequency flood and the second split ratio corresponding to different working condition schemes, the regulation scheme of the target river section is determined. For details, please refer to Figure 2 Step S206 of the embodiment shown will not be repeated here.

[0057] The method for calculating the river split ratio provided in this embodiment uses a two-dimensional water and sediment model, introduces long series of daily water and sediment data and modified terrain data corresponding to each working condition scheme, obtains predicted river terrain data in the future period under different working conditions, and then introduces flood flow data to obtain a plurality of sets of second hydrological information under each working condition. The plurality of sets of second hydrological information can be compared with the first hydrological information to obtain the influence of a working condition on the future river bed split ratio under fixed bed conditions and dynamic bed conditions, or the plurality of sets of second hydrological information can be compared with each other to determine the influence of different working conditions on the future river bed split ratio under dynamic bed conditions.

[0058] In this embodiment, a method for calculating a river split ratio is provided, Figure 6 is a flow chart of a method for calculating a river split ratio according to an embodiment of the present application, as shown in the figure, the flow includes the following steps: Figure 6 Step S401, long series of daily water and sediment data, initial terrain data, modified terrain data corresponding to different working condition schemes, and flow when different frequency floods occur of the target river section are obtained. For details, please refer to Figure 2 Step S201 of the embodiment shown will not be repeated here.

[0059] Step S402, under fixed bed conditions, the long series of daily water and sediment data, the initial terrain data, and the flow when different frequency floods occur are calculated by using a two-dimensional water and sediment mathematical model to obtain the first hydrological information of the target river section in the future period. For details, please refer to Figure 2 Step S202 of the embodiment shown will not be repeated here.

[0060] Step S403, under dynamic bed conditions, the long series of daily water and sediment data, the modified terrain data corresponding to each working condition scheme, and the flow when different frequency floods occur are calculated by using a two-dimensional water and sediment mathematical model to obtain the second hydrological information of the target river section when each frequency flood occurs in the future period when different working condition schemes are implemented on the target river section. For details, please refer to Figure 2 ​Step S203 of the illustrated embodiment will not be described here again.

[0061] Step S404, for each frequency flood, each section in the target river is divided into a plurality of sub-sections according to the first hydrological information and the flow of each sub-section is calculated, including: Step S4041, a plurality of section lines are determined according to the pre-established section line file, and the end point coordinates corresponding to each section line; The section line file records a plurality of section line information, wherein the section line information includes the end point coordinates of the two ends of the section line, and the connecting line of the end point coordinates of the section line is the section line. The end point coordinates of the section line should be set outside the area of the study river section to ensure that the section line can completely cover the study river section.

[0062] Exemplarily, the section line file is as shown in Figure 7 Each row is a section line, the first column indicates the section number, the second and third columns are the x and y coordinates of one end of the section line, and the fourth and fifth columns are the x and y coordinates of the other end of the section line.

[0063] Step S4042, the tree data structure algorithm is used to determine the river boundary point closest to the end point coordinates in the first hydrological information to obtain the section corresponding to each section line. Wherein, determining the river boundary point closest to the end point coordinates means determining the two intersection points of the section line and the river boundary, and taking the two intersection points as the section end points of the section to which the section line belongs.

[0064] Step S4043, the section is divided into a plurality of sub-sections.

[0065] In an optional embodiment, the section can be divided into a plurality of sub-sections by artificial division, each sub-section serving as a calculation unit and representing a tributary. The subsequent calculation of the flow ratio is also calculated in units of sub-sections.

[0066] Exemplarily, in the example as shown in Figure 8 The section is divided into 5 sub-sections.

[0067] Step S405, for each frequency flood, each section in the target river is divided into a plurality of sub-sections according to the second hydrological information and the flow of each sub-section is calculated, and the second flow ratio of the section is determined according to the flow ratio of each sub-section in the section. For details, please refer to step S205 of the illustrated embodiment. Figure 2 Step S205 of the illustrated embodiment will not be described here again.

[0068] Step S406, according to the first flow ratio corresponding to each frequency flood and the second flow ratio corresponding to different working condition schemes, a control scheme of the target river section is determined. For details, please refer toFigure 2 The step S206 of the embodiment is not described herein again.

[0069] The method for calculating the river diversion ratio provided in the embodiment provides a scheme for determining a cross section by a cross point of a cross section line and a river boundary, which can ensure that the cross section can contain the flow of all tributaries of the study river section, and gives a scheme for manually dividing the cross section into sub-cross sections, which is used for flexible calculation of the diversion ratio.

[0070] In the embodiment, a method for calculating a river diversion ratio is provided, Figure 9 is a flow chart of the method for calculating the river diversion ratio according to the embodiment of the application, as shown in the figure, the flow includes the following steps: Figure 9 The flow includes the following steps: Step S501, obtaining long series daily water and sediment data of a target river section, initial terrain data, modified terrain data corresponding to different working condition schemes, and flow when different frequency floods occur. For details, please refer to Figure 2 The step S201 of the embodiment is not described herein again.

[0071] Step S502, under the condition of fixed bed, using a water and sediment two-dimensional mathematical model to calculate the long series daily water and sediment data, the initial terrain data, and the flow when different frequency floods occur, to obtain first hydrological information of the target river section in a future period. For details, please refer to Figure 2 The step S202 of the embodiment is not described herein again.

[0072] Step S503, under the condition of moving bed, using a water and sediment two-dimensional mathematical model to calculate the long series daily water and sediment data, the modified terrain data corresponding to each working condition scheme, and the flow when different frequency floods occur, to obtain second hydrological information of the target river section in the future period when each frequency flood occurs when the target river section is implemented with different working condition schemes. For details, please refer to Figure 2 The step S203 of the embodiment is not described herein again.

[0073] Step S504, the first hydrological information contains water depth, horizontal flow velocity and vertical flow velocity at different positions in the target river section; the step of calculating the flow of each sub-cross section according to the first hydrological information includes: Step S5041, dividing the cross section into a plurality of small sections; Exemplarily, the cross section can be divided into 500 small sections, and the number of small sections is not limited herein. The more the number of small sections is divided, the higher the calculation accuracy is.

[0074] Step S5042, performing interpolation operation on the first hydrological information to obtain water depth, horizontal flow velocity and vertical flow velocity of each small section; The water depth of each section refers to the vertical distance from the bottom of the water body to the water surface in the section. The horizontal flow velocity refers to the speed and direction of the water flow particles moving in the horizontal direction. The vertical flow velocity refers to the speed of the water flow particles moving in the vertical direction.

[0075] In step S5043, the flow velocity projection of each section is calculated according to the normal vector, the horizontal flow velocity and the vertical flow velocity of each section.

[0076] The normal vector of each section refers to the unit normal vector of the cross-sectional area of each section.

[0077] Exemplarily, the flow velocity projection of each section can be calculated by uv = u * fx + v * fy, where h is the water depth of each section, u is the horizontal flow velocity, v is the vertical flow velocity, f is the normal vector, uv is the flow velocity projection result, fx is the projection of the normal vector in the x direction, and fy is the projection of the normal vector in the y direction.

[0078] In step S5044, the flow of each section is calculated according to the flow velocity projection, the length and the water depth of each section.

[0079] Exemplarily, the flow of each section can be calculated by q = uv * L * h. Where q is the total flow of the section, uv is the flow velocity projection result, L is the length of each section, and h is the water depth of the section.

[0080] In step S5045, the flows of all sections in the same sub-section are added to obtain the flow of each sub-section.

[0081] This process first determines the sub-section to which each section belongs, and then adds the flows of all sections in the same sub-section. The sum of the flows of all sections in the same sub-section is the flow of the sub-section.

[0082] In step S505, for each frequency flood, each section in the target river is divided into multiple sub-sections according to the second hydrological information, and the flow of each sub-section is calculated, and the second split ratio of the section is determined according to the flow ratio of each sub-section in the section. For details, please refer to step S205 of the embodiment shown in Figure 2 The step S205 of the embodiment shown in is not repeated here.

[0083] In step S506, the control scheme of the target river section is determined according to the first split ratio corresponding to each frequency flood and the second split ratio corresponding to different working conditions. For details, please refer to step S206 of the embodiment shown in Figure 2 The step S206 of the embodiment shown in is not repeated here.

[0084] The method for calculating the river split ratio provided in this embodiment can improve the calculation accuracy by dividing the section into multiple sections, calculating and accumulating the flow of each section, and the flow velocity distribution in the river is not uniform in the calculation process of the river flow.

[0085] A method for calculating a river diversion ratio is provided in the embodiments of the present application, Figure 10 is a flow chart of a method for calculating a river diversion ratio according to an embodiment of the present application, as shown in the figure, the flow includes the following steps: Figure 10 In step S601, long series of daily water and sediment data of a target river section, initial terrain data, modified terrain data corresponding to different working condition schemes, and flow when floods of different frequencies occur are obtained. For details, please refer to step S201 of the embodiment shown in Figure 2

[0086] In step S602, under the condition of fixed bed, the long series of daily water and sediment data, the initial terrain data, and the flow when floods of different frequencies occur are calculated by using a two-dimensional mathematical model of water and sediment, to obtain first hydrological information of the target river section in a future period. For details, please refer to step S202 of the embodiment shown in Figure 2

[0087] In step S603, under the condition of movable bed, the long series of daily water and sediment data, the modified terrain data corresponding to each working condition scheme, and the flow when floods of different frequencies occur are calculated by using a two-dimensional mathematical model of water and sediment, to obtain second hydrological information of the target river section in a future period when floods of each frequency occur while implementing different working condition schemes on the target river section. For details, please refer to step S203 of the embodiment shown in Figure 2

[0088] In step S604, for each frequency flood, each cross section in the target river is divided into multiple sub-cross sections according to the first hydrological information, and the flow of each sub-cross section is calculated, and the first diversion ratio of the cross section is determined according to the flow ratio of each sub-cross section in the cross section. For details, please refer to step S204 of the embodiment shown in Figure 2

[0089] In step S605, for each frequency flood, each cross section in the target river is divided into multiple sub-cross sections according to the second hydrological information, and the flow of each sub-cross section is calculated, including: In step S6051, a plurality of cross section lines and the endpoint coordinates corresponding to each cross section line are determined according to a pre-established cross section line file.

[0090] The cross section line file records a plurality of cross section line information, wherein the cross section line information includes endpoint coordinates of both ends of the cross section line, and the connecting line of the endpoint coordinates of the cross section line is the cross section line. The endpoint coordinates of the cross section line should be set outside the area of the study river section, so as to ensure that the cross section line can completely cover the study river section.

[0091] ​​​​​Step S6052, determining the river boundary point closest to the endpoint coordinate in the second hydrological information by using the tree data structure algorithm, to obtain the cross section corresponding to each cross section line.

[0092] The determination of the river boundary point closest to the endpoint coordinate refers to determining two intersection points of the cross section line and the river boundary, and taking the two intersection points as the cross section endpoints of the cross section to which the cross section line belongs.

[0093] Step S6053, dividing the cross section into a plurality of sub-cross sections.

[0094] In an optional embodiment, the cross section can be divided into a plurality of sub-cross sections in a manually divided manner, each sub-cross section representing a tributary as a calculation unit, and the subsequent calculation of the flow distribution ratio is also performed in units of the sub-cross sections.

[0095] Step S606, determining the regulation scheme of the target river section according to the first flow distribution ratio corresponding to each frequency flood and the second flow distribution ratio corresponding to different working condition schemes. For details, please refer to Figure 2 Step S206 of the embodiment shown in FIG. 2 will not be repeated here.

[0096] The method for calculating the river flow distribution ratio provided in this embodiment provides a scheme for determining the cross section through the intersection points of the cross section line and the river boundary, which can ensure that the cross section can contain the flow of all tributaries of the study river section. The scheme for manually dividing the cross section into sub-cross sections is provided for flexible calculation of the flow distribution ratio.

[0097] In this embodiment, a method for calculating the river flow distribution ratio is provided, Figure 11 is a flowchart of a method for calculating the river flow distribution ratio according to an embodiment of the present application, as shown in FIG. 1, which includes the following steps: Figure 11 Step S701, obtaining the long series daily water and sediment data of the target river section, the initial terrain data, the modified terrain data corresponding to different working condition schemes, and the flow when different frequency floods occur. For details, please refer to Figure 2 Step S201 of the embodiment shown in FIG. 2 will not be repeated here.

[0098] Step S702, under the condition of fixed bed, calculating the long series daily water and sediment data, the initial terrain data, and the flow when different frequency floods occur by using the two-dimensional mathematical model of water and sediment, to obtain the first hydrological information of the target river section in the future period. For details, please refer to Figure 2 Step S202 of the embodiment shown in FIG. 2 will not be repeated here.

[0099] ​Step S703, under the moving bed condition, the long series of daily water and sediment data, the modified terrain data corresponding to each working condition scheme, and the flow when different frequency floods occur are calculated by using the two-dimensional mathematical model of water and sediment, to obtain the second hydrological information of the target river section in the future period when each frequency flood occurs when the different working condition schemes are implemented on the target river section. For details, please refer to Figure 2 Step S203 of the embodiment shown will not be repeated here.

[0100] Step S704, for each frequency flood, each section in the target river is divided into a plurality of sub-sections according to the first hydrological information, and the flow of each sub-section is calculated, and the first flow distribution ratio of the section is determined according to the flow ratio of each sub-section in the section. For details, please refer to Figure 2 Step S204 of the embodiment shown will not be repeated here.

[0101] Step S705, the second hydrological information includes the water depth, horizontal flow velocity and vertical flow velocity at different positions in the target river section; the step of calculating the flow of each sub-section according to the second hydrological information includes: Step S7051, the section is divided into a plurality of small sections.

[0102] Exemplarily, the section can be divided into 500, 600, 700 small sections, which is not limited here. The more the number of small sections is divided, the higher the calculation accuracy is.

[0103] Step S7052, the second hydrological information is subjected to an interpolation operation to obtain the water depth, horizontal flow velocity and vertical flow velocity of each small section.

[0104] The water depth of each small section refers to the vertical distance from the bottom of the water body to the water surface in the small section. The horizontal flow velocity refers to the speed and direction of the water flow particle moving in the horizontal direction. The vertical flow velocity refers to the speed of the water flow particle moving in the vertical direction.

[0105] Step S7053, the normal vector, horizontal flow velocity and vertical flow velocity of each small section are calculated to calculate the flow velocity projection of each small section.

[0106] The normal vector of each small section refers to the unit normal vector of the flow area of each small section.

[0107] Exemplarily, the flow velocity projection of each small section can be calculated by uv = u * fx + v * fy, wherein h is the water depth of each small section, u is the horizontal flow velocity, v is the vertical flow velocity, f is the normal vector, uv is the flow velocity projection result, fx is the projection of the normal vector in the x direction, and fy is the projection of the normal vector in the y direction.

[0108] Step S7054, the flow of each small section is calculated according to the flow velocity projection, length and water depth of each small section.

[0109] Exemplarily, the flow of each section can be calculated by q=uv*L*h, where q is the total flow of the section, uv is the flow velocity projection result, L is the length of each section, and h is the water depth of the section.

[0110] In step S7055, the flows of all the sections in the same sub-section are added to obtain the flow of each sub-section.

[0111] The process first determines the sub-section to which each section belongs, and adds the flows of all the sections in the same sub-section to obtain the flow of the sub-section.

[0112] In step S706, the control scheme of the target river section is determined according to the first diversion ratio corresponding to each frequency flood and the second diversion ratio corresponding to different working condition schemes. For details, refer to Figure 2 The step S206 of the embodiment shown will not be described herein.

[0113] The method for calculating the river diversion ratio provided in the embodiment can improve the calculation accuracy by dividing the section into multiple sections, calculating and accumulating the flows of the sections respectively, when the water flow distribution in the river is uneven in the river flow calculation process.

[0114] In the embodiment, a device for calculating the river diversion ratio is also provided, which is used to implement the above embodiments and preferred embodiments, and will not be described herein. As used below, the term'module' can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, hardware or a combination of software and hardware is also possible and is contemplated.

[0115] The device for calculating the river diversion ratio provided in the embodiment includes: Figure 12 as shown, comprising: The data acquisition module 801 is configured to acquire data required for calculation in the river to be measured.

[0116] The data calculation module 802 is configured to calculate the river diversion ratio according to the data required for calculation in the river to be measured.

[0117] The device for calculating the river diversion ratio provided in the embodiment can perform the method for calculating the river diversion ratio provided in any embodiment of the present application, and has the corresponding function modules and advantages of performing the method. The further function description of each module and unit is the same as that of the corresponding embodiment, and will not be described herein.

[0118] Figure 13 The structure schematic diagram of an electronic device provided in the embodiment of the present application is shown.

[0119] The following will be specifically describedFigure 13 which shows a structural schematic diagram suitable for use to implement an electronic device in embodiments of the present application. The electronic device can include a processor (e.g., a central processing unit, a graphics processing unit, etc.) 901, which can perform various appropriate actions and processes according to programs stored in a read-only memory (ROM) 902 or programs loaded from a storage 908 into a random access memory (RAM) 903. Various programs and data required for operation of the electronic device are also stored in the RAM 903. The processor 901, the ROM 902, and the RAM 903 are connected to each other through a bus 904. An input / output (I / O) interface 905 is also connected to the bus 904.

[0120] Generally, the following devices can be connected to the I / O interface 905: an input device 906 including, for example, a touch screen, a touch pad, a keyboard, a mouse, a camera, a microphone, an accelerometer, a gyroscope, etc.; an output device 907 including, for example, a liquid crystal display (LCD), a speaker, a vibrator, etc.; a storage 908 including, for example, a magnetic tape, a hard disk, etc.; and a communication device 909. The communication device 909 can allow the electronic device to communicate wirelessly or wiredly with other devices to exchange data. Although Figure 13 The electronic device with various devices is shown, but it should be understood that all of the shown devices are not required to be implemented or possessed, and more or less devices can be alternatively implemented or possessed.

[0121] In particular, the processes described above with reference to the flowcharts can be implemented as a computer software program according to embodiments of the present application. For example, embodiments of the present application include a computer program product including a computer program carried on a non-transitory computer readable medium, the computer program containing program code for executing the method shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network through the communication device 909, or installed from the storage 908, or installed from the ROM 902. When the computer program is executed by the processor 901, the above-described functions defined in the method of diverting a flow of a river according to embodiments of the present application are performed.

[0122] Figure 13 The electronic device shown is merely an example and should not impose any limitation on the functions and the range of use of embodiments of the present application.

[0123] The embodiments of the present application further provide a computer readable storage medium, and the method according to the embodiments of the present application can be implemented in hardware, firmware, or recorded in a storage medium, or stored in a remote storage medium or a non-transitory machine readable storage medium and downloaded to a local storage medium through network, so that the method described herein can be processed by such software on a storage medium using a general purpose computer, a special purpose processor, or programmable or special hardware. The storage medium can be a magnetic disk, an optical disk, a read-only memory, a random access memory, a flash memory, a hard disk, or a solid state disk, etc. Further, the storage medium can also include a combination of the above-mentioned memories. It can be understood that the computer, the processor, the microprocessor controller, or the programmable hardware includes a storage component that can store or receive software or computer code, when the software or computer code is accessed and executed by the computer, the processor, or the hardware, a method for calculating the riverway diversion ratio shown in the above embodiments is implemented.

[0124] Part of the present application can be applied as a computer program product, for example, computer program instructions, when executed by a computer, the operation of the computer can call or provide the method and / or technical solutions according to the present application. Those skilled in the art should understand that the form of computer program instructions in computer readable medium includes but is not limited to source file, executable file, installation package file, etc. Correspondingly, the way of computer program instructions executed by computer includes but is not limited to: the computer directly executes the instructions, or the computer compiles the instructions and then executes the corresponding compiled program, or the computer reads and executes the instructions, or the computer reads and installs the instructions and then executes the corresponding installed program. Here, the computer readable medium can be any available computer readable storage medium or communication medium accessible to the computer.

[0125] Although the embodiments of the present application are described in conjunction with the accompanying drawings, various modifications and changes can be made by those skilled in the art without departing from the spirit and scope of the present application, and such modifications and changes fall within the scope defined by the appended claims.

Claims

1. A method of analog computation of the split ratio in a river course, characterized by, The method comprises: obtaining long series of daily water and sediment data of a target river section, initial terrain data, modified terrain data corresponding to different working condition schemes respectively, and flow when different frequency floods occur; under the condition of fixed bed, calculating the long series of daily water and sediment data, the initial terrain data and the flow when different frequency floods occur by using a two-dimensional mathematical model of water and sediment, to obtain first hydrological information of the target river section in a future period; under the condition of moving bed, calculating the long series of daily water and sediment data, the modified terrain data corresponding to each working condition scheme respectively and the flow when different frequency floods occur by using the two-dimensional mathematical model of water and sediment, to obtain second hydrological information of the target river section in the future period when each frequency flood occurs while implementing different working condition schemes on the target river section; for each frequency flood, dividing each section in the target river channel into multiple sub-sections according to the first hydrological information and calculating the flow of each sub-section, and determining a first split ratio of the section according to the flow ratio of each sub-section in the section; for each frequency flood, dividing each section in the target river channel into multiple sub-sections according to the second hydrological information and calculating the flow of each sub-section while implementing different working condition schemes on the target river section, and determining a second split ratio of the section according to the flow ratio of each sub-section in the section; determining a regulation scheme of the target river section according to the first split ratio corresponding to each frequency flood and the second split ratio corresponding to different working condition schemes.

2. The method of claim 1, wherein, under the condition of moving bed, calculating the long series of daily water and sediment data and the modified terrain data corresponding to each working condition scheme respectively by using the two-dimensional mathematical model of water and sediment, to obtain multiple groups of third hydrological information of the target river section in the future period, one group of third hydrological information corresponding to one working condition scheme, and the third hydrological information containing elevations of different positions of the target river section; extracting evolution terrain data of different working condition schemes according to the elevations of different positions in the third hydrological information; inputting the evolution terrain data of different working condition schemes and the flow when different frequency floods occur into the two-dimensional mathematical model of water and sediment, to obtain the second hydrological information of the target river section in the future period when each frequency flood occurs while implementing different working condition schemes on the target river section under the condition of moving bed. for each frequency flood, dividing each section in the target river channel into multiple sub-sections according to the first hydrological information and calculating the flow of each sub-section, and determining a first split ratio of the section according to the flow ratio of each sub-section in the section, comprises:

3. The method of claim 1, wherein, determining multiple section lines and endpoint coordinates corresponding to each section line according to a pre-established section line file; determining the section corresponding to each section line by using a tree data structure algorithm to determine the river boundary point closest to the endpoint coordinates in the first hydrological information; dividing the section into multiple sub-sections. ​ 4. The method according to claim 1 or 3, characterized in that, The first hydrological information comprises water depth, horizontal flow velocity and vertical flow velocity at different positions in the target river section; the step of calculating the flow of each sub-section according to the first hydrological information comprises: dividing the section into a plurality of small sections; interpolating the first hydrological information to obtain water depth, horizontal flow velocity and vertical flow velocity of each small section; calculating flow velocity projection of each small section according to normal vector, horizontal flow velocity and vertical flow velocity of each small section; calculating flow of each small section according to flow velocity projection, length and water depth of each small section; adding up flow of all small sections in the same sub-section to obtain flow of each sub-section.

5. The method of claim 1, wherein, When different working condition schemes are implemented on the target river section, for each frequency flood, each section in the target river is divided into a plurality of sub-sections and flow of each sub-section is calculated according to the second hydrological information, comprising: determining a plurality of section lines and end point coordinates corresponding to each section line according to a pre-established section line file; determining river boundary points closest to the end point coordinates in the second hydrological information by using a tree data structure algorithm to obtain sections corresponding to each section line; dividing the sections into a plurality of sub-sections.

6. The method according to claim 1 or 5, characterized in that, The second hydrological information comprises water depth, horizontal flow velocity and vertical flow velocity at different positions in the target river section; the step of calculating the flow of each sub-section according to the second hydrological information comprises: dividing the section into a plurality of small sections; interpolating the second hydrological information to obtain water depth, horizontal flow velocity and vertical flow velocity of each small section; calculating flow velocity projection of each small section according to normal vector, horizontal flow velocity and vertical flow velocity of each small section; calculating flow of each small section according to flow velocity projection, length and water depth of each small section; adding up flow of all small sections in the same sub-section to obtain flow of each sub-section.

7. A device for simulating and calculating the flow split ratio within a river channel, characterized in that, The device comprises: a data acquisition module for acquiring data required for calculation in a river to be measured; a data calculation module for calculating river diversion ratio in the river to be measured according to the data required for calculation in the river to be measured.

8. An electronic device, comprising: comprise: a memory and a processor, which are communicatively connected, the memory stores computer instructions, and the processor executes the computer instructions to perform the method for simulating calculation of river diversion ratio in any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer readable storage medium stores computer instructions for causing a computer to perform the method for simulating calculation of river diversion ratio in any one of claims 1 to 6.

10. A computer program product, characterised in that, comprise computer instructions for causing a computer to perform the method for simulating calculation of river diversion ratio in any one of claims 1 to 6.

Images