River channel erosion and deposition feature judgment method and device, computer equipment and storage medium

Abstract

The invention relates to the technical field of reservoir regulation and control, and discloses a river channel erosion and deposition characteristic judgment method and device, computer equipment and a storage medium. Flow topographic information of a target river channel is acquired, and a hydrodynamic model is established based on the flow topographic information and a preset initial roughness factor set; respectively inputting the respectively obtained water flow information of the first time period and the second time period into a hydrodynamic model corresponding to the target river channel to obtain a test result of the first water level changing along with time and a test result of the second water level changing along with time; adjusting the initial roughness factor set based on the first water level change-over-time test result to obtain a first roughness factor set, and adjusting the initial roughness factor set based on the second water level change-over-time test result to obtain a second roughness factor set; based on the first roughness factor set and the second roughness factor set, it is determined that the target river channel is in the scouring or silting state, the scouring and silting state of the target river channel can be accurately determined, and meanwhile resource consumption and time cost are remarkably reduced.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of reservoir regulation, and particularly relates to a river channel erosion and deposition feature judgment method and device, computer equipment and a storage medium. BACKGROUND

[0002] The judgment method of river channel sediment deposition generally includes a measurement method and a mathematical model calculation method. The sediment deposition can cause the water level to rise, and in particular, in the case where the water volume is unchanged, whether the deposition occurs can be determined by water level monitoring, cross section measurement or riverbed measurement, or the mathematical model calculation method can simulate the sediment transport and deposition process by establishing a water and sediment movement model. The accuracy of the measurement method is often related to the sampling point position, and the accuracy deviation can be caused by the dependence on the sampling point. The mathematical model simulates the sediment transport and deposition process, and needs to combine a large amount of historical data and hydrological information, and the calculation resources and time cost are more. SUMMARY

[0003] Therefore, the present application provides a river channel erosion and deposition feature judgment method and device, computer equipment and a storage medium to solve the problem of accuracy deviation of the existing deposition judgment method and more calculation resources and time cost.

[0004] In a first aspect, the present application provides a river channel erosion and deposition feature judgment method, which includes: obtaining the flow topography information of a target river channel, and establishing a fluid dynamics model corresponding to the target river channel based on the flow topography information and a preset initial roughness set corresponding to the target river channel; obtaining the water flow information of the target river channel in a first time period and a second time period respectively, wherein there is an abnormal change in water flow between the first time period and the second time period, and the water flow information at least includes the actual water level change result of the target river channel in the first time period and the second time period respectively; inputting the water flow information in the first time period into the fluid dynamics model corresponding to the target river channel to obtain a first water level change test result, and inputting the water flow information in the second time period into the fluid dynamics model corresponding to the target river channel to obtain a second water level change test result; adjusting the initial roughness set based on the first water level change test result and the actual water level change result of the target river channel in the first time period to obtain a first roughness set, and adjusting the initial roughness set based on the second water level change test result and the actual water level change result of the target river channel in the second time period to obtain a second roughness set; and determining whether the target river channel is in an erosion or deposition state based on the first roughness set and the second roughness set.

[0005] In an optional implementation, the fluid dynamic model comprises a steady flow dynamic model and a non-steady flow dynamic model, and the establishing of the fluid dynamic model corresponding to the target river channel based on the flow topography information and preset initial roughness comprises: dividing the target river channel into a plurality of river sections based on all the observation stations along the target river channel as section nodes; setting the initial roughness of each river section based on the flow topography characteristics corresponding to the river section; and establishing the steady flow dynamic model and the non-steady flow dynamic model corresponding to the target river channel based on the flow topography information of each river section and the corresponding initial roughness.

[0006] In an optional implementation, the flow information comprises at least the inflow of the target river channel and the water level before the dam, and the inputting of the flow information of the first time period into the fluid dynamic model corresponding to the target river channel to obtain the first water level change test result comprises: taking the inflow of the target river channel in the first time period as the upper boundary condition and the water level before the dam of the target river channel in the first time period as the lower boundary condition, and calculating the first water level value set of each position along the target river channel by using the steady flow dynamic model; and inputting the first water level value set of each position into the non-steady flow dynamic model based on the inflow of the target river channel in the first time period as the upper boundary condition and the water level before the dam of the target river channel in the first time period as the lower boundary condition, to obtain the first water level change test result.

[0007] In an optional implementation, the flow information comprises at least the inflow of the target river channel and the water level before the dam, and the inputting of the flow information of the second time period into the fluid dynamic model corresponding to the target river channel to obtain the second water level change test result comprises: taking the inflow of the target river channel in the second time period as the upper boundary condition and the water level before the dam of the target river channel in the second time period as the lower boundary condition, and calculating the second water level value set of each position along the target river channel by using the steady flow dynamic model; and inputting the second water level value set of each position into the non-steady flow dynamic model based on the inflow of the target river channel in the second time period as the upper boundary condition and the water level before the dam of the target river channel in the second time period as the lower boundary condition, to obtain the second water level change test result.

[0008] In an optional embodiment, the initial roughness set is composed of initial roughnesses set for each river section, and the initial roughness set is adjusted based on the first water level variation with time test result and the actual water level variation with time result of the target river channel in the first time period to obtain a first roughness set, or the initial roughness set is adjusted based on the second water level variation with time test result and the actual water level variation with time result of the target river channel in the second time period to obtain a second roughness set, including: based on the target condition that the first water level variation with time test result is consistent with the actual water level variation with time result of the target river channel in the first time period, the initial roughness set for each river section is adjusted, and the adjusted roughness corresponding to each river section is combined into the first roughness set, or based on the target condition that the second water level variation with time test result is consistent with the actual water level variation with time result of the target river channel in the second time period, the initial roughness set for each river section is adjusted, and the adjusted roughness corresponding to each river section is combined into the second roughness set.

[0009] In an optional embodiment, determining whether the target river channel is in a scouring or silting state based on the first roughness set and the second roughness set includes: extracting the first roughness and the second roughness corresponding to the current river section from the first roughness set and the second roughness set respectively; determining whether the river section is in a scouring or silting state based on the first roughness and the second roughness; repeating the steps of extracting the first roughness and the second roughness corresponding to the current river section from the first roughness set and the second roughness set respectively until the scouring or silting state corresponding to all river sections in the target river channel is obtained.

[0010] In an optional embodiment, based on the first roughness set and the second roughness set, it is determined whether the target river channel is in a scouring or silting state by the following formula:

[0011]

[0012] Where r represents the rate of change. When r is positive, the target river channel is determined to be in a siltation state. When r is negative, the target river channel is determined to be in a scouring state. K1 represents the second roughness set; K0 represents the first roughness set.

[0013] In a second aspect, the present application provides a river erosion and deposition feature determination device, which comprises: a dynamic model establishing module, configured to acquire flow topographic information of a target river, and establish a fluid dynamic model corresponding to the target river based on the flow topographic information and an initial roughness set corresponding to the target river; a water flow information acquiring module, configured to acquire water flow information of the target river in a first time period and a second time period respectively, wherein there is an abnormal change in water flow between the first time period and the second time period, and the water flow information at least comprises actual water level change results of the target river in the first time period and the second time period respectively; a water level change test module, configured to input the water flow information in the first time period into the fluid dynamic model corresponding to the target river to obtain a first water level change test result, and input the water flow information in the second time period into the fluid dynamic model corresponding to the target river to obtain a second water level change test result; a roughness adjusting module, configured to adjust the initial roughness set based on the first water level change test result and the actual water level change result of the target river in the first time period to obtain a first roughness set, and adjust the initial roughness set based on the second water level change test result and the actual water level change result of the target river in the second time period to obtain a second roughness set; and a river state determining module, configured to determine whether the target river is in an erosion or deposition state based on the first roughness set and the second roughness set.

[0014] In a third aspect, the present application provides a computer device, which comprises a memory and a processor, the memory and the processor are communicatively connected with each other, the memory stores computer instructions, and the processor executes the computer instructions to perform the river erosion and deposition feature determination method of the first aspect or any of the corresponding embodiments thereof.

[0015] In a fourth aspect, the present application provides a computer readable storage medium, which stores computer instructions, and the computer instructions are used to make a computer perform the river erosion and deposition feature determination method of the first aspect or any of the corresponding embodiments thereof.

[0016] The present application has the following technical effects:

[0017] The river channel erosion and deposition feature judgment method provided in the embodiment, obtains flow topographic information of a target river channel, and establishes a fluid dynamics model corresponding to the target river channel based on the flow topographic information and an initial roughness set corresponding to the target river channel; obtains flow information of the target river channel in a first time period and a second time period respectively; inputs the flow information in the first time period into the fluid dynamics model corresponding to the target river channel to obtain a first water level change with time test result, and inputs the flow information in the second time period into the fluid dynamics model corresponding to the target river channel to obtain a second water level change with time test result; adjusts the initial roughness set based on the first water level change with time test result and an actual water level change with time result of the target river channel in the first time period to obtain a first roughness set, and adjusts the initial roughness set based on the second water level change with time test result and the actual water level change with time result of the target river channel in the second time period to obtain a second roughness set; and determines that the target river channel is in an erosion or deposition state based on the first roughness set and the second roughness set, which not only can accurately determine the erosion and deposition state of the target river channel, but also significantly reduces resource consumption and time cost. BRIEF DESCRIPTION OF DRAWINGS

[0018] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, hereinafter, a brief introduction will be given to the drawings needed to be used in the description of the embodiments or the prior art. Obviously, the drawings described hereinafter are some embodiments of the present application, and other drawings can also be obtained by those skilled in the art without any creative effort on the basis of these drawings.

[0019] Figure 1 is a flowchart of a river channel erosion and deposition feature judgment method according to an embodiment of the present application;

[0020] Figure 2 is an example diagram of roughness and river channel features according to an embodiment of the present application;

[0021] Figure 3a is an example diagram of comparison of flow rates in June and September according to an embodiment of the present application;

[0022] Figure 3b is an example diagram of comparison of reservoir water levels in June and September according to an embodiment of the present application;

[0023] Figure 4 is a flowchart of another river channel erosion and deposition feature judgment method according to an embodiment of the present application;

[0024] Figure 5 is a structural block diagram of a river channel erosion and deposition feature judgment device according to an embodiment of the present application;

[0025] Figure 6Schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0026] To make the purpose, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without making creative efforts shall fall within the scope of protection of the present invention.

[0027] According to an embodiment of the present invention, an embodiment of a method for determining river channel scouring and deposition characteristics is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0028] In this embodiment, a method for determining the characteristics of river scouring and silting is provided, which can be used in the above-mentioned computer device. Figure 1 FIG. 1 is a flow chart of a method for determining river scouring and silting characteristics according to an embodiment of the present invention. Figure 1 As shown, the process includes the following steps:

[0029] Step S101 : acquiring flow topography information of a target river channel, and establishing a fluid dynamics model corresponding to the target river channel based on the flow topography information and a preset initial roughness set corresponding to the target river channel.

[0030] The embodiment of the present invention can collect flow topographic information of the target river channel, wherein the flow topographic information can include data such as the cross-sectional morphology of the river channel (such as riverbed height, bank slope, river channel width), water level of each observation station, measured inflow flow, outflow flow and topographic cross-sectional data. For example, for a large river channel reservoir, the entire reservoir area is about 572 kilometers, and the roughness along the entire length is unknown. Roughness is a parameter that characterizes the roughness of the riverbed surface and affects the water flow resistance. The corresponding initial roughness set can be set for the target river channel according to factors such as riverbed material (rock, sediment, vegetation), river channel morphology (straight or curved), etc. The specific determination of the initial roughness based on river channel characteristics can refer to the following. Figure 2 The reference roughness setting shown in FIG. 4 is shown, wherein the initial roughness set can be a single-element set containing only one value, that is, the initial roughness can be set for the target river channel as a whole based on the riverbed state of the entire river channel as the initial roughness set. The target river channel can also be divided into multiple river channel areas, and the corresponding initial roughness can be set for each river channel area based on the riverbed state and other information of each river channel area to form an initial roughness set without limitation.

[0031] The embodiment of the present invention is based on the basic equations of fluid mechanics and simplifies the river channel into a one-dimensional, two-dimensional or three-dimensional mathematical model. Taking the one-dimensional hydrodynamic model as an example, a large river-type reservoir with tributaries flowing into the reservoir area is established. When establishing a one-dimensional fluid dynamic model of the reservoir, the flow movement of the main and tributary rivers is taken into account at the same time. The main and tributary rivers of the reservoir are regarded as single river channels, and the confluence points of the rivers become bifurcations. Then the one-dimensional hydrodynamic model may include three parts: the single river channel flow motion equation, the bifurcations connection equation and the boundary conditions.

[0032] Among them, the water flow continuity equation used in the fluid dynamics model to describe the water flow motion is:

[0033]

[0034] The equation of motion for water flow is:

[0035]

[0036] Among them, the subscript i is the section number; Q is the flow rate; Z is the water level; A is the water area; q L is the lateral inflow per unit length of the river section; t is time; x is the coordinate along the flow path; K is the roughness.

[0037] The flow rate in and out of each branch point must be balanced with the actual increase or decrease rate of water volume in the branch point, that is:

[0038]

[0039] Among them, Ω is the water storage capacity of the branch point. If the point is generalized as a geometric point, then Ω=0.

[0040] Dynamic connection condition: If the branch point can be generalized as a geometric point, the water flow in and out of each branch channel is gentle, and there is no sudden change in water level, then the water level of each branch channel section should be equal, that is:

[0041]

[0042] Boundary conditions: In the calculation, boundary conditions are not given for a single river channel separately. Instead, boundary conditions are given for the main and tributary rivers of the reservoir included in the calculation scope as a whole. The flow process is given for the inlet of each main and tributary river, and the water level process or flow process or water level-flow relationship is given for the model outlet.

[0043] Numerical method: Use the four-point implicit difference format to discretize. The four-point implicit format is also known as the Preissmann format, which is a forward difference format that averages (or weighted averages) the four adjacent points. The derivative of t is taken as the average of the forward time difference quotients on the adjacent nodes, and the derivative of x is taken as the average or weighted average of the forward spatial difference quotients of the adjacent two layers. For the M points in the grid, let

[0044]

[0045] wherein f is any function, f i n is f(x i , t n ); theta is a weighting factor having a value that is a positive number less than or equal to 1.

[0046] In step S102, the water flow information of the target river channel in the first time period and the second time period is obtained respectively, wherein there is an abnormal change of water flow between the first time period and the second time period, and the water flow information at least includes the actual water level change result of the target river channel in the first time period and the second time period respectively.

[0047] The embodiment of the present application considers that the sediment mainly concentrates in the flood season, so if it is necessary to evaluate the erosion and deposition of the target river channel, the state before and after the arrival of the sediment before and after the flood can be collected, not only the water level and flow close, but also the sediment concentration close and less, in order to reduce the interference caused by the difference of the flow capacity of the river channel boundary at different elevations, such as Figure 3a and Figure 3b As shown in the formula, the selected reservoir water level change process level2 in September is close to the selected reservoir water level change process level1 in June, and the water level difference is at most ±1m, in order to reduce the influence of the inflow flow difference on the flow capacity of the river channel boundary, the selected inflow flow Q2 in September should be close to the selected inflow flow Q2 change process in June, and the difference between the inflow flows is at most 1000m 3 / s, through comparison, the two time periods of June 3-6 and September 3-6 can meet the requirements of reservoir water level and inflow flow at the same time; therefore, the water flow information of the target river channel in June 3-6 (the first time period) and September 3-6 (the second time period) can be obtained respectively, wherein the water flow can include but is not limited to inflow flow, water level before the dam and actual water level change result with time.

[0048] In step S103, the water flow information in the first time period is input into the fluid dynamics model corresponding to the target river channel to obtain the first water level change test result, and the water flow information in the second time period is input into the fluid dynamics model corresponding to the target river channel to obtain the second water level change test result.

[0049] The embodiment of the present application can input the water flow in the first time period into the fluid dynamics model corresponding to the target river channel, so as to obtain the first water level change test result, and can input the water flow in the second time period into the fluid dynamics model corresponding to the target river channel, so as to obtain the second water level change test result.

[0050] Step S104, based on the first water level change test result and the actual water level change result of the target river channel in the first time period, adjusting the initial roughness set to obtain the first roughness set, and based on the second water level change test result and the actual water level change result of the target river channel in the second time period, adjusting the initial roughness set to obtain the second roughness set.

[0051] The embodiment of the present application can compare the first water level change test result with the actual water level change result of the target river channel in the first time period, adjust the initial roughness set to obtain the first roughness set, such as roughness, large flow, low water level, and the need to reduce roughness, only as an example, the second water level change test result can also be compared with the actual water level change result of the target river channel in the second time period, and the initial roughness set is adjusted to obtain the second roughness set.

[0052] Step S105, based on the first roughness set and the second roughness set, determining whether the target river channel is in a scouring or silting state.

[0053] The embodiment of the present application generally shows the form of scouring state that the riverbed sediment is carried away by the water flow, avoiding the decrease of roughness, corresponding to the decrease of roughness; the form of silting state is that the sediment deposition causes the riverbed to be lifted, the surface roughness increases, and the roughness value increases, so the change of the first roughness set and the second roughness set can be compared to determine whether the roughness increases or decreases, and then the target river channel is determined to be in a scouring state or a silting state, only as an example.

[0054] The river channel scouring and silting feature judgment method provided by the embodiment acquires the flow topographic information of the target river channel, and establishes a fluid dynamic model corresponding to the target river channel based on the flow topographic information and the initial roughness set corresponding to the target river channel; the flow information of the target river channel in the first time period and the second time period is acquired respectively; the flow information in the first time period is input into the fluid dynamic model corresponding to the target river channel to obtain the first water level change test result, and the flow information in the second time period is input into the fluid dynamic model corresponding to the target river channel to obtain the second water level change test result; based on the first water level change test result and the actual water level change result of the target river channel in the first time period, the initial roughness set is adjusted to obtain the first roughness set, and based on the second water level change test result and the actual water level change result of the target river channel in the second time period, the initial roughness set is adjusted to obtain the second roughness set; based on the first roughness set and the second roughness set, it is determined whether the target river channel is in a scouring or silting state, which not only can accurately determine the scouring and silting state of the target river channel, but also significantly reduces the resource consumption and time cost.

[0055] The embodiment provides a river channel scouring and silting feature judgment method, which can be used for the computer device, Figure 4 is a flow chart of the river channel scouring and silting feature judgment method according to the embodiment of the application, as shown in the figure, and the flow chart comprises the following steps: Figure 4

[0056] In step S401, the flow topography information of the target river channel is acquired, and a fluid dynamic model corresponding to the target river channel is established based on the flow topography information and an initial roughness set corresponding to the target river channel.

[0057] Specifically, the fluid dynamic model comprises a steady flow dynamic model and an unsteady flow dynamic model, and the step S401 comprises the following steps.

[0058] In step S4011, the target river channel is divided into a plurality of river sections based on all along-observation stations of the target river channel as section nodes.

[0059] In step S4012, the initial roughness of each river section is set based on the flow topography features of the river section.

[0060] In step S4013, the steady flow dynamic model and the unsteady flow dynamic model corresponding to the target river channel are established based on the flow topography information of each river section and the corresponding initial roughness.

[0061] According to the embodiment of the application, the target river channel is divided into a plurality of river sections from the upstream to the downstream according to the distribution of the along-observation stations of the target river channel and the change of the hydraulic characteristics, that is, the along-observation stations are taken as the section nodes, as shown in Table 1, so that the roughness of each river section can be adjusted subsequently, and the initial roughness of each river section can be set based on the flow topography features (roughness river bed composition and bed surface features, plane shape and water flow state and bank wall features) of each river section, as shown in Table 2. Figure 2 Then, the steady flow dynamic model and the unsteady flow dynamic model corresponding to the target river channel can be established based on the flow topography information of each river section and the corresponding initial roughness, and the river flow characteristics can be described from the steady state and the dynamic state respectively.

[0062] Table 1

[0063]

[0064]

[0065] ​The present invention divides the river channel by using observation stations along the river as nodes, and sets an initial roughness for each river channel to dynamically adapt to the characteristics of the river section, which can match different river channel terrain conditions. A constant flow dynamics model is established to finely analyze the steady-state water flow and provide initial steady-state conditions for the non-constant flow dynamics model. The non-constant flow dynamics model dynamically simulates the river channel scouring and silting process. The combination of the constant flow dynamics model and the non-constant flow dynamics model can accurately predict the river channel scouring and silting status.

[0066] Step S402: Obtain the flow information of the target river in the first time period and the second time period respectively. Figure 1 Step S102 of the illustrated embodiment will not be described in detail here.

[0067] There is an abnormal change in water flow between the first time period and the second time period, and the water flow information at least includes the actual water level change results of the target river channel in the first time period and the second time period respectively.

[0068] Step S403: input the water flow information of the first time period into the fluid dynamics model corresponding to the target river channel to obtain the first water level change test result over time, and input the water flow information of the second time period into the fluid dynamics model corresponding to the target river channel to obtain the second water level change test result over time.

[0069] The water flow information includes at least the incoming flow of the target river channel and the water level in front of the dam.

[0070] Specifically, the inflow flow of the target river channel in the first time period is used as the upper boundary condition, and the water level in front of the dam of the target river channel in the first time period is used as the lower boundary condition, and the first water level value set at each position along the target river channel is calculated through the constant flow dynamic model; based on the inflow flow of the target river channel in the first time period as the upper boundary condition, and the water level in front of the dam of the target river channel in the first time period as the lower boundary condition, the first water level value set at each position is input into the non-constant flow dynamic model to obtain the test results of the first water level changing with time.

[0071] Specifically, the inflow flow of the target river channel in the second time period is used as the upper boundary condition, and the water level in front of the dam of the target river channel in the second time period is used as the lower boundary condition, and the second water level value set at each position along the target river channel is calculated through the constant flow dynamic model; based on the inflow flow of the target river channel in the second time period as the upper boundary condition, and the water level in front of the dam of the target river channel in the second time period as the lower boundary condition, the second water level value set at each position is input into the non-constant flow dynamic model to obtain the test results of the second water level changing with time.

[0072] The embodiment of the present application can take the inflow flow of the target river channel in the first time period (June) as the upper boundary condition and the water level before the dam of the target river channel as the lower boundary condition, so as to calculate a first water level value set of each position along the target river channel by the constant flow dynamic model. The first water level value set of each position along the river channel can be taken as the initial condition of the unsteady flow dynamic model, and the inflow flow of the target river channel in the first time period is taken as the upper boundary condition and the water level before the dam of the target river channel is taken as the lower boundary condition, so as to obtain the first water level change test result corresponding to the first time period by calculation of the unsteady flow dynamic model.

[0073] The embodiment of the present application can take the inflow flow of the target river channel in the second time period (September) as the upper boundary condition and the water level before the dam of the target river channel as the lower boundary condition, so as to calculate a second water level value set of each position along the target river channel by the constant flow dynamic model. The second water level value set of each position along the river channel can be taken as the initial condition of the unsteady flow dynamic model, and the inflow flow of the target river channel in the second time period is taken as the upper boundary condition and the water level before the dam of the target river channel is taken as the lower boundary condition, so as to obtain the second water level change test result corresponding to the second time period by calculation of the unsteady flow dynamic model.

[0074] The present application utilizes the constant flow dynamic model to obtain the steady water level value conforming to the measured boundary, and utilizes the unsteady flow dynamic model to capture the dynamic response of the water level with time, so as to realize the leap from the static working condition analysis to the dynamic process prediction, and to provide more accurate time sequence basis for the scientific regulation and control of the river channel system.

[0075] In step S404, the initial roughness set is adjusted based on the first water level change test result and the actual water level change result of the target river channel in the first time period to obtain a first roughness set, and the initial roughness set is adjusted based on the second water level change test result and the actual water level change result of the target river channel in the second time period to obtain a second roughness set.

[0076] The initial roughness set is composed of the initial roughness set respectively set for each river section.

[0077] Specifically, the initial roughness set respectively set for each river section is adjusted based on the target condition that the first water level change test result is consistent with the actual water level change result of the target river channel in the first time period, and the adjusted roughness of each river section forms the first roughness set, or the initial roughness set respectively set for each river section is adjusted based on the target condition that the second water level change test result is consistent with the actual water level change result of the target river channel in the second time period, and the adjusted roughness of each river section forms the second roughness set.

[0078] In an embodiment of the present invention, the initial roughness value corresponding to each river section can be adjusted so that the test result of the first water level change over time is consistent with the result of the actual water level change over time of the target river channel in the first time period, that is, the adjusted roughness of each river section is combined into a first roughness set. At the same time, the initial roughness value corresponding to each river section can also be adjusted so that the result of the second water level change over time is consistent with the result of the actual water level change over time of the target river channel in the second time period, and the adjusted roughness of each river section is combined into a second roughness set, as shown in Table 1, including each river section, the corresponding initial roughness, the first roughness set, and the second roughness set.

[0079] This method adjusts the initial roughness to ensure that the results of the first water level change test are consistent with the actual water level change over time in the target river during the first time period. This yields an adjusted roughness, achieving a transition from "fixed parameter simulation" to "adaptive dynamic prediction." This method not only improves the simulation accuracy of the hydrodynamic model but also, through time-series comparison of roughness sets, creates a quantitative analysis tool for river channel evolution. This provides both scientific and practical technical support for multi-objective decision-making, such as flood control safety, ecological protection, and water resource management.

[0080] Step S405: extract the first roughness and the second roughness corresponding to the current river section from the first roughness set and the second roughness set respectively; determine whether the river section is in a scouring or silting state based on the first roughness and the second roughness; repeat the steps of extracting the first roughness and the second roughness corresponding to the current river section from the first roughness set and the second roughness set respectively until the scouring or silting state corresponding to all river sections in the target river channel is obtained.

[0081] In the embodiment of the present invention, based on the first roughness and second roughness corresponding to each river channel, it can be determined whether the river section is in a scouring or silting state. For example, as shown in Table 1, the first roughness corresponding to river section D is 23.4, and the second roughness is 22.5. The decrease in roughness indicates that river section D is in a scouring state. As an example only, by judging the changes in the roughness corresponding to all river sections, it can be determined whether all river sections in the target river channel are in a scouring state or a silting state.

[0082] Furthermore, based on the first roughness set and the second roughness set, the following formula is used to determine whether the target river channel is in a scouring or silting state:

[0083]

[0084] Where r represents the rate of change. When r is positive, the target river channel is determined to be in a siltation state. When r is negative, the target river channel is determined to be in a scouring state. K1 represents the second roughness set; K0 represents the first roughness set.

[0085] The calculation formula for river flow capacity is:

[0086]

[0087] Among them, Q is the flow rate of a certain section in the river, m 3 / s; A is the water flow area of ​​a certain section in the river, m 2 ; R is the hydraulic radius, m; i is the slope of the riverbed, and K is the roughness.

[0088] According to the formula for calculating the flow capacity of a river channel in the embodiment of the present invention, when the incoming flow rate and water level remain the same and the terrain does not change, that is, when Q, A, R, and i remain unchanged, the roughness should remain unchanged. When the incoming flow rate and water level remain the same and the roughness changes, in order to maintain the consistency of Q, R, and i, A changes, that is, the water flow area of ​​a certain section in the river channel changes, and the rate of change r is:

[0089]

[0090] When r is positive, it indicates that siltation has occurred in the river section, and when r is negative, it indicates that erosion has occurred in the river section.

[0091] As shown in Table 1, the roughness K0 of river section D is 23.4 and K1 is 22.5. Substituting into the above formula, it can be seen that scour occurred in river section D with a scour rate of 3.8%; the roughness K0 and K1 of other river sections did not change, indicating that no obvious scouring and deposition changes occurred.

[0092] By querying the topographic section data of this large river-type reservoir in different years before and after the construction, it was found that among the sections where scouring and deposition changes occurred, scouring occurred in the D section of the river, with a corresponding scouring rate of about 5%, indicating that the scouring and deposition trend and magnitude of the river channel determined by this method are relatively reliable.

[0093] This embodiment also provides a device for determining river scouring and silting characteristics. This device is used to implement the above-mentioned embodiments and preferred embodiments, and details already described will not be repeated. As used below, the term "module" may refer to a combination of software and / or hardware that implements a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, implementation using hardware, or a combination of software and hardware, is also possible and contemplated.

[0094] This embodiment provides a device for determining the characteristics of river scouring and silting. Figure 5As shown, the method comprises: a dynamic model establishing module 501, configured to acquire flow topographic information of a target river channel, and establish a fluid dynamic model corresponding to the target river channel based on the flow topographic information and a preset initial roughness set corresponding to the target river channel; a water flow information acquiring module 502, configured to acquire water flow information of the target river channel in a first time period and a second time period respectively, wherein there is an abnormal change in water flow between the first time period and the second time period, and the water flow information at least comprises actual water level change results of the target river channel in the first time period and the second time period respectively; a water level change test module 503, configured to input the water flow information in the first time period into the fluid dynamic model corresponding to the target river channel to obtain a first water level change test result, and input the water flow information in the second time period into the fluid dynamic model corresponding to the target river channel to obtain a second water level change test result; a roughness adjusting module 504, configured to adjust the initial roughness set based on the first water level change test result and the actual water level change result of the target river channel in the first time period to obtain a first roughness set, and adjust the initial roughness set based on the second water level change test result and the actual water level change result of the target river channel in the second time period to obtain a second roughness set; and a river channel state determining module 505, configured to determine that the target river channel is in an erosion or deposition state based on the first roughness set and the second roughness set.

[0095] In some optional embodiments, the fluid dynamic model comprises a steady flow dynamic model and an unsteady flow dynamic model, and the water flow information acquiring module 502 comprises: a river section dividing unit, configured to divide the target river channel into a plurality of river sections based on all along-the-way observation stations of the target river channel as section nodes; a roughness setting unit, configured to set respective initial roughnesses for the respective river sections based on flow topographic features corresponding to the respective river sections; and a model establishing unit, configured to establish the steady flow dynamic model and the unsteady flow dynamic model corresponding to the target river channel based on the flow topographic information of the respective river sections and the respective initial roughnesses.

[0096] In some optional embodiments, the water flow information at least comprises inflow flow and dam front water level of the target river channel, and the water level change test module 503 comprises: a first water level value calculating unit, configured to calculate a first water level value set of each position of the target river channel by the steady flow dynamic model, with the inflow flow of the target river channel in the first time period as an upper boundary condition and the dam front water level of the target river channel in the first time period as a lower boundary condition; and a first water level change unit, configured to input the first water level value set of each position into the unsteady flow dynamic model to obtain the first water level change test result, with the inflow flow of the target river channel in the first time period as the upper boundary condition and the dam front water level of the target river channel in the first time period as the lower boundary condition.

[0097] In some optional embodiments, the water flow information includes at least the inflow flow of the target river channel and the water level in front of the dam, and the water level change over time test module 503 includes: a second water level value calculation unit, which is used to use the inflow flow of the target river channel in the second time period as the upper boundary condition and the water level in front of the dam of the target river channel in the second time period as the lower boundary condition, and calculate the second water level value set of each position along the target river channel through a constant flow dynamic model; a second water level change over time unit, which is used to input the second water level value set of each position into the non-constant flow dynamic model based on the inflow flow of the target river channel in the second time period as the upper boundary condition and the water level in front of the dam of the target river channel in the second time period as the lower boundary condition, to obtain the second water level change over time test result.

[0098] In some optional embodiments, the initial roughness set is composed of the initial roughness set for each river section, and the roughness adjustment module 504 includes: a roughness adjustment unit, which is used to adjust the initial roughness set for each river section based on the target condition that the first water level change over time test result is consistent with the actual water level change over time result of the target river channel in the first time period, and form the corresponding adjusted roughness of each river section into a first roughness set, or, based on the target condition that the second water level change over time test result is consistent with the actual water level change over time result of the target river channel in the second time period, adjust the initial roughness set for each river section, and form the corresponding adjusted roughness of each river section into a second roughness set.

[0099] In some optional embodiments, the river channel state determination module 505 includes: a roughness extraction unit, used to extract the first roughness and the second roughness corresponding to the current river section from the first roughness set and the second roughness set respectively; a river section state determination unit, used to determine whether the river section is in a scouring or silting state based on the first roughness and the second roughness; a river channel state determination unit, used to repeatedly execute the steps of extracting the first roughness and the second roughness corresponding to the current river section from the first roughness set and the second roughness set respectively until the scouring or silting state corresponding to all river sections in the target river channel is obtained.

[0100] In some optional embodiments, based on the first roughness set and the second roughness set, it is determined whether the target river channel is in a scouring or silting state by the following formula:

[0101]

[0102] Where r represents the rate of change. When r is positive, the target river channel is determined to be in a siltation state. When r is negative, the target river channel is determined to be in a scouring state. K1 represents the second roughness set; K0 represents the first roughness set.

[0103] The further functional description of each of the above modules and units is the same as that of the above corresponding embodiments and will not be repeated here.

[0104] The river channel scouring and silting characteristics judgment device in this embodiment is presented in the form of a functional unit, where the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and memory that executes one or more software or fixed programs, and / or other devices that can provide the above functions.

[0105] The embodiment of the present invention also provides a computer device having the above Figure 5 The device for judging the scouring and silting characteristics of a river channel is shown.

[0106] See also Figure 6 , Figure 6 is a structural diagram of a computer device provided by an optional embodiment of the present invention, such as Figure 6 As shown, the computer device includes: one or more processors 10, memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. Various components utilize different buses to communicate with each other and can be installed on a common mainboard or installed in other ways as needed. The processor can process the instructions executed in the computer device, including instructions stored in the memory or on the memory to display the graphical information of the GUI on an external input / output device (such as, a display device coupled to the interface). In some optional embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Equally, multiple computer devices can be connected, and each device provides part of the necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Figure 6 A processor 10 is taken as an example.

[0107] The processor 10 may be a central processing unit, a network processor, or a combination thereof. The processor 10 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit, a programmable logic device, or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic, or any combination thereof.

[0108] The memory 20 stores instructions that can be executed by at least one processor 10, so as to enable at least one processor 10 to execute the method shown in the above embodiment.

[0109] The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system and application programs required for at least one function; the data storage area may store data created based on the use of the computer device, etc. In addition, the memory 20 may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some optional embodiments, the memory 20 may optionally include a memory remotely located relative to the processor 10, and these remote memories may be connected to the computer device via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0110] The memory 20 may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a flash memory, a hard disk or a solid-state drive; the memory 20 may also include a combination of the above types of memory.

[0111] The computer device further includes an input device 30 and an output device 40. The processor 10, the memory 20, the input device 30 and the output device 40 may be connected via a bus or other means. Figure 6 The bus connection is taken as an example.

[0112] The input device 30 can receive input digital or character information and generate key signal input related to user settings and function control of the computer device, such as a touch screen, a keypad, a mouse, a trackpad, a touch pad, an indicator stick, one or more mouse buttons, a trackball, a joystick, etc. The output device 40 can include a display device, an auxiliary lighting device (e.g., an LED), and a tactile feedback device (e.g., a vibration motor). The above-mentioned display device includes but is not limited to a liquid crystal display, a light emitting diode, a display, and a plasma display. In some optional embodiments, the display device can be a touch screen.

[0113] The embodiment of the present invention also provides a computer-readable storage medium. The above-mentioned method according to the embodiment of the present invention can be implemented in hardware, firmware, or implemented as a computer code that can be recorded in a storage medium, or implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and downloaded through a network and will be stored in a local storage medium, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state drive, etc.; further, the storage medium can also include a combination of the above-mentioned types of memory. It can be understood that a computer, a processor, a microprocessor controller or 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 a computer, a processor or hardware, the method shown in the above embodiment is implemented.

[0114] A portion of the present invention may be applied as a computer program product, such as a computer program instruction, which, when executed by a computer, can call or provide the method and / or technical solution according to the present invention through the operation of the computer. Those skilled in the art should understand that the form in which the computer program instruction exists in a computer-readable medium includes, but is not limited to, a source file, an executable file, an installation package file, etc. Accordingly, the way in which the computer program instruction is executed by the computer includes, but is not limited to: the computer directly executes the instruction, or the computer compiles the instruction and then executes the corresponding compiled program, or the computer reads and executes the instruction, or the computer reads and installs the instruction and then executes the corresponding installed program. Here, the computer-readable medium may be any available computer-readable storage medium or communication medium that can be accessed by the computer.

[0115] Although the embodiments of the present invention have been described with reference to the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the present invention. Such modifications and variations are all within the scope defined by the appended claims.

Claims

1. A method for determining river scouring and silting characteristics, characterized in that: The method comprises: Acquiring flow topography information of a target river channel, and establishing a fluid dynamics model corresponding to the target river channel based on the flow topography information and a preset initial roughness set corresponding to the target river channel; Obtaining flow information of a target river channel in a first time period and a second time period, respectively, wherein there is an abnormal change in flow between the first time period and the second time period, and the flow information at least includes a result of an actual water level change over time in the target river channel in each of the first time period and the second time period; Inputting the water flow information of the first time period into the fluid dynamics model corresponding to the target river channel to obtain a first water level change over time test result, and inputting the water flow information of the second time period into the fluid dynamics model corresponding to the target river channel to obtain a second water level change over time test result; adjusting the initial roughness set based on the first water level variation over time test result and the actual water level variation over time result of the target river channel in the first time period to obtain a first roughness set; and adjusting the initial roughness set based on the second water level variation over time test result and the actual water level variation over time result of the target river channel in the second time period to obtain a second roughness set; Based on the first roughness set and the second roughness set, it is determined that the target river channel is in a scouring or silting state.

2. The method according to claim 1, characterized in that The fluid dynamic model includes a steady flow dynamic model and a non-steady flow dynamic model. The fluid dynamic model corresponding to the target river channel is established based on the flow topography information and the preset initial roughness, including: The target river is divided into multiple sections based on all the observation stations along the target river as segment nodes; Based on the flow and topographic characteristics of each river section, the corresponding initial roughness is set for each river section; Based on the flow topography information of each river section and the corresponding initial roughness, a steady flow dynamic model and a non-steady flow dynamic model corresponding to the target river channel are established.

3. The method according to claim 2, characterized in that The water flow information includes at least the incoming flow of the target river channel and the water level in front of the dam. The water flow information of the first time period is input into the fluid dynamic model corresponding to the target river channel to obtain a test result of the first water level changing over time, including: Using the inflow flow of the target river channel in the first time period as an upper boundary condition and the water level in front of the dam of the target river channel in the first time period as a lower boundary condition, a set of first water level values ​​at various locations along the target river channel is calculated using the constant flow dynamics model; Based on the inflow flow of the target river channel in the first time period as the upper boundary condition, and the water level in front of the dam of the target river channel in the first time period as the lower boundary condition, the set of first water level values ​​at each position is input into the non-steady flow dynamic model to obtain the test results of the first water level change over time.

4. The method according to claim 2, characterized in that The water flow information includes at least the incoming flow of the target river channel and the water level in front of the dam. The water flow information of the second time period is input into the fluid dynamic model corresponding to the target river channel to obtain a test result of the second water level changing with time, including: Using the inflow flow of the target river channel in the second time period as an upper boundary condition and the water level in front of the dam of the target river channel in the second time period as a lower boundary condition, a set of second water level values ​​at various locations along the target river channel is calculated using the constant flow force model; Based on the inflow flow of the target river channel in the second time period as the upper boundary condition, and the water level in front of the dam of the target river channel in the second time period as the lower boundary condition, the set of second water level values ​​at each position is input into the non-steady flow dynamic model to obtain the test results of the second water level change over time.

5. The method according to claim 3 or 4, characterized in that The initial roughness set is composed of initial roughnesses set for respective river sections. The initial roughness set is adjusted based on the first water level variation over time test result and the actual water level variation over time result of the target river channel in the first time period to obtain the first roughness set, or the initial roughness set is adjusted based on the second water level variation over time test result and the actual water level variation over time result of the target river channel in the second time period to obtain the second roughness set, including: Based on the target condition that the first water level variation over time test result is consistent with the actual water level variation over time result of the target river channel in the first time period, the initial roughness set for each river section is adjusted, and the adjusted roughness corresponding to each river section is combined into a first roughness set, or, Based on the target condition that the second water level change over time test result is consistent with the actual water level change over time result of the target river channel in the second time period, the initial roughness set for each river section is adjusted, and the adjusted roughness corresponding to each river section is combined into a second roughness set.

6. The method according to claim 5, characterized in that The determining, based on the first roughness set and the second roughness set, whether the target river channel is in a scouring or silting state includes: Extracting the first roughness and the second roughness corresponding to the current river section from the first roughness set and the second roughness set respectively; Based on the first roughness and the second roughness, determining that the river section is in a scouring or silting state; The steps of respectively extracting the first roughness and the second roughness corresponding to the current river section from the first roughness set and the second roughness set are repeated until the scouring or silting states corresponding to all river sections in the target river channel are obtained.

7. The method according to claim 1, characterized in that Based on the first roughness set and the second roughness set, the target river channel is determined to be in a scouring or silting state using the following formula: Where r represents the rate of change. When r is positive, the target river channel is determined to be in a siltation state. When r is negative, the target river channel is determined to be in a scouring state. K1 represents the second roughness set; K0 represents the first roughness set.

8. A device for judging the characteristics of river scouring and silting, characterized in that: The device comprises: a dynamic model building module, configured to obtain flow topography information of a target river channel, and to build a fluid dynamic model corresponding to the target river channel based on the flow topography information and a preset initial roughness set corresponding to the target river channel; a water flow information acquisition module, configured to respectively acquire water flow information of a target river channel in a first time period and a second time period, wherein an abnormal change in water flow occurs between the first time period and the second time period, and the water flow information at least includes a result of an actual water level change over time in the target river channel in each of the first time period and the second time period; a water level variation over time test module, configured to input the water flow information of a first time period into the fluid dynamics model corresponding to the target river channel to obtain a first water level variation over time test result, and to input the water flow information of a second time period into the fluid dynamics model corresponding to the target river channel to obtain a second water level variation over time test result; a roughness adjustment module, configured to adjust the initial roughness set based on the first water level variation over time test result and the actual water level variation over time result of the target river channel in the first time period to obtain a first roughness set, and to adjust the initial roughness set based on the second water level variation over time test result and the actual water level variation over time result of the target river channel in the second time period to obtain a second roughness set; The channel state determination module is configured to determine whether the target channel is in a scouring or silting state based on the first roughness set and the second roughness set.

9. A computer device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the river scouring and silting characteristics judgment method according to any one of claims 1 to 7 by executing the computer instructions.

10. A computer-readable storage medium, characterized in that The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the river scouring and silting characteristics judgment method according to any one of claims 1 to 7.

Images