Reservoir siltation control method, device, electronic equipment and storage medium

Abstract

The present invention discloses a reservoir siltation control method, device, electronic device and storage medium, including: obtaining the water-sand movement process and water-sand movement model within a target time period, and calculating the reservoir sediment discharge ratio, variable backwater area sedimentation, and perennial backwater area sedimentation; determining the sediment discharge ratio index, variable backwater area sedimentation ratio index, and perennial backwater area sedimentation ratio index based on the water-sand movement model, reservoir sediment discharge ratio, variable backwater area sedimentation, and perennial backwater area sedimentation, and then determining the monotonicity of the corresponding functions to obtain the avoidance water level corresponding to each index; determining the target avoidance water level corresponding to the target sedimentation control demand from the avoidance water levels corresponding to each index according to the target sedimentation control demand, and determining the target water level interval based on the target avoidance water level to control the reservoir sedimentation. This method ensures that the reservoir can control the sedimentation in the target area to be as small as possible, and will not cause significant sedimentation.

Description

Technical Field

[0001] The present invention relates to the field of water conservancy projects, and in particular to a reservoir siltation control method, device, electronic equipment and storage medium. Background Art

[0002] In order to make full use of the water energy resources generated by the natural terrain drop and to intercept floods, many reservoirs were built to block rivers and form river-type reservoirs. The reservoir area of ​​this type of reservoir is generally divided into a perennial backwater area and a variable backwater area. Among them, the perennial backwater area is close to the dam, is located below the water surface for a long time, and is in a submerged state; the variable backwater area is relatively far from the dam, and switches between submerged and natural river forms due to changes in the reservoir water level. Therefore, the variable backwater area and the perennial backwater area show different sedimentation trends due to different water flow states.

[0003] When a flood carries a large amount of sediment into the reservoir area, the reservoir will intercept the flood and store water, which will inevitably lead to an increase in the reservoir water level. However, the rise in the reservoir water level will inevitably lead to increased sedimentation. However, due to the influence of three main factors: river topography, inflow and inflow sediment, the increase in sedimentation in the reservoir area caused by the reservoir rising to different water levels is not the same, and there are also differences in the impact on the variable backwater area, the perennial backwater area and the reservoir sediment discharge ratio. At present, the measures for sedimentation control and sediment reduction of reservoirs are usually aimed at sedimentation control of the entire reservoir, without considering the various areas of the reservoir, so the effect of sedimentation control is not ideal. Summary of the invention

[0004] In view of this, an embodiment of the present invention provides a reservoir siltation control method, device, electronic device and storage medium to solve the problem of poor reservoir siltation control effect.

[0005] According to a first aspect, an embodiment of the present invention provides a reservoir siltation control method, comprising:

[0006] Obtaining the water-sediment movement process and the water-sediment movement model within the target time period, and determining the reservoir sediment discharge ratio, the variable backwater area sedimentation, and the perennial backwater area sedimentation based on the water-sediment movement process and the water-sediment movement model;

[0007] Determine the sediment discharge ratio index, the sedimentation ratio index of the variable backwater area, and the sedimentation ratio index of the permanent backwater area based on the water-sediment movement model, the reservoir sedimentation ratio, the sedimentation amount of the variable backwater area, and the sedimentation amount of the permanent backwater area;

[0008] Calculate the functions and function monotonicity corresponding to the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index, respectively, and determine the avoidance water level corresponding to each index based on the function monotonicity;

[0009] Obtain a target siltation control demand, and based on the target siltation control demand, determine a target avoidance water level corresponding to the target siltation control demand from the avoidance water levels corresponding to each indicator, and determine a target water level interval based on the target avoidance water level to control reservoir siltation.

[0010] The reservoir siltation control method provided by the embodiment of the present invention determines the reservoir sediment discharge ratio, the variable backwater area sedimentation amount, and the perennial backwater area sedimentation amount, and then establishes the sediment discharge ratio index, the variable backwater area sedimentation ratio index, and the perennial backwater area sedimentation ratio index. The monotonicity of the function corresponding to each index is analyzed to determine the avoidance water level corresponding to each index. Different indicators are determined for different areas and siltation reduction needs. When siltation control is required through water level regulation, the corresponding target avoidance water level is selected and the target water level range is determined according to the target siltation control demand, to ensure that the reservoir can control the sedimentation in the target area as little as possible under the target avoidance water level, and will not cause significant sedimentation.

[0011] In some embodiments, the determining of the reservoir sediment discharge ratio, the amount of siltation in the variable backwater area, and the amount of siltation in the perennial backwater area based on the water-sediment movement process and the water-sediment movement model includes:

[0012] Determine the flow and sediment content of the inlet river section according to the water and sediment movement process, and calculate the inlet sediment volume based on the flow and sediment content of the inlet river section;

[0013] Determine the flow and sediment content of the outflow river section based on the water and sediment movement model, and calculate the outflow sediment volume based on the flow and sediment content of the outflow river section;

[0014] Determine the reservoir sediment discharge ratio according to the ratio of the outgoing sediment volume to the incoming sediment volume;

[0015] The flow and sediment content of the backwater zone boundary section are determined based on the water-sediment movement model, and the siltation amount of the variable backwater zone and the siltation amount of the perennial backwater zone are determined based on the flow and sediment content of the backwater zone boundary section.

[0016] In some embodiments, the determining of the flow and sediment content of the backwater zone boundary section based on the water-sediment movement model, and the determining of the variable backwater zone siltation amount and the perennial backwater zone siltation amount based on the flow and sediment content of the backwater zone boundary section, comprises:

[0017] Calculating the sediment transport capacity of the backwater zone boundary section based on the flow and sediment content of the backwater zone boundary section;

[0018] The amount of sedimentation in the variable backwater area is determined based on the difference between the amount of sand entering the reservoir and the amount of sand discharged at the boundary section of the backwater area.

[0019] The annual siltation volume of the backwater area is determined based on the difference between the sediment transport volume at the boundary section of the backwater area and the sediment discharge volume.

[0020] In some embodiments, the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index are determined according to the following formula:

[0021]

[0022]

[0023] Among them, i represents the i-th level of operating water level, and the reservoir operating water level can be divided into n levels, i=1 represents the flood limit water level, and i=n represents the flood control high water level; P represents the reservoir sediment discharge ratio; B represents the siltation in the variable backwater area; C represents the siltation in the perennial backwater area; I P Indicates the sediment discharge ratio index; I B Indicates the proportion of sedimentation in the changing backwater area; I C It indicates the proportion of siltation in the perennial backwater area; the subscript min indicates that the reservoir is at the flood limit water level; the subscript max indicates that the reservoir is at the flood control high water level; It represents the ratio of the sediment discharge ratio attenuation at the i-th operating water level to the total attenuation from the water level min to max; It indicates the ratio of the sedimentation in the variable backwater area under the i-th operating water level to the total sedimentation from the water level min to max; It represents the ratio of the annual siltation in the backwater area at the i-th operating water level to the total siltation when the water level ranges from min to max.

[0024] In some embodiments, the determining of the avoidance water level corresponding to each indicator based on the monotonicity of the functions corresponding to the sediment discharge ratio indicator, the variable backwater area siltation ratio indicator, and the perennial backwater area siltation ratio indicator, respectively, includes:

[0025] Delimiting at least one operating water level;

[0026] Calculate the sediment discharge ratio corresponding to the operating water level, the variable backwater area sedimentation ratio, and the perennial backwater area sedimentation ratio respectively;

[0027] Draw function graphs based on the calculated values ​​of the sediment discharge ratio index, the calculated values ​​of the variable backwater area siltation ratio index, and the calculated values ​​of the perennial backwater area siltation ratio index, respectively, to determine the sediment discharge ratio index function, the variable backwater area siltation ratio index function, and the perennial backwater area siltation ratio index function;

[0028] The monotonicity of the sediment discharge ratio index function, the variable backwater area siltation ratio index function and the perennial backwater area siltation ratio index function are analyzed respectively, and the avoidance water level corresponding to each index is determined based on the monotonicity.

[0029] In some embodiments, when the target siltation control requirement includes more than one indicator, the method further includes:

[0030] Get the weight coefficient corresponding to each indicator;

[0031] Determine a weighted index for reducing siltation based on the weight coefficients corresponding to the various indicators;

[0032] Dividing at least one target water level, and calculating the calculated values ​​of the siltation reduction weighted index corresponding to different target water levels, so as to determine a function of the target water level with respect to the siltation reduction weighted index;

[0033] Analyzing the monotonicity of the function of the weighted siltation reduction index to determine the avoided water level of the weighted siltation reduction index;

[0034] Based on the target siltation control demand, a target avoidance water level corresponding to the target siltation control demand is determined from the avoidance water level of the siltation reduction weighted index, and a target water level interval is determined based on the target avoidance water level.

[0035] In some embodiments, the weighted index for reducing siltation is determined according to the following formula:

[0036]

[0037] Among them, W i represents the weighted index of silt reduction, i represents the i-th level operating water level, I P Indicates the sediment discharge ratio index; I B Indicates the proportion of sedimentation in the changing backwater area; I C It represents the proportion of siltation in the backwater area throughout the year, a 1 、a 2 、a 3 are weight coefficients, where a 1 +a 2 +a 3 =1.

[0038] According to a second aspect, an embodiment of the present invention provides a reservoir siltation control device, comprising:

[0039] A data acquisition module is used to acquire the water-sediment movement process and the water-sediment movement model within the target time period, and determine the reservoir sediment discharge ratio, the variable backwater area siltation, and the perennial backwater area siltation based on the water-sediment movement process and the water-sediment movement model;

[0040] An index determination module, for determining a sediment discharge ratio index, a sediment accumulation ratio index of a variable backwater area, and a sediment accumulation ratio index of a perennial backwater area based on the water-sediment movement model, the reservoir sediment discharge ratio, the sediment accumulation amount of a variable backwater area, and the sediment accumulation amount of a perennial backwater area;

[0041] A water level calculation module, used to calculate the functions corresponding to the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index and the function monotonicity, and determine the avoidance water level corresponding to each index based on the function monotonicity;

[0042] The target water level determination module is used to obtain the target siltation control demand, and based on the target siltation control demand, determine the target avoidance water level corresponding to the target siltation control demand from the avoidance water levels corresponding to each indicator, and determine the target water level interval based on the target avoidance water level to control the siltation of the reservoir.

[0043] According to the third aspect, an embodiment of the present invention provides an electronic device, comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the reservoir siltation control method described in the first aspect or any one of the embodiments of the first aspect by executing the computer instructions.

[0044] According to a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable the computer to execute the reservoir siltation control method described in the first aspect or any one embodiment of the first aspect. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0046] Figure 1 is a flow chart of a reservoir sedimentation control method according to an embodiment of the present invention;

[0047] Figure 2 is a water-sand movement process prediction diagram according to an embodiment of the present invention;

[0048] Figure 3 is a graph of a sediment discharge ratio index function according to an embodiment of the present invention;

[0049] Figure 4 is a graph of a variable backwater area sedimentation ratio index function according to an embodiment of the present invention;

[0050] Figure 5 is a function graph of the index of the proportion of siltation in the perennial backwater area according to an embodiment of the present invention;

[0051] Figure 6is a weighted index function image of silt reduction according to an embodiment of the present invention;

[0052] Figure 7 is a block diagram of a reservoir sedimentation control device according to an embodiment of the present invention;

[0053] Figure 8 It is a schematic diagram of the hardware structure of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0054] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the 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 creative work are within the scope of protection of the present invention.

[0055] According to an embodiment of the present invention, an embodiment of a reservoir siltation control method 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.

[0056] At present, measures for reservoir siltation control and siltation reduction include the use of mechanical devices and reservoir scheduling. Mechanical means to remove silt in the reservoir area are suitable for short-term and small-scale implementation, otherwise the cost is high. The reservoir scheduling method has low cost and can alleviate siltation in the long term. The existing reservoir scheduling method increases the discharge flow and lowers the reservoir water level by scheduling the siltation reduction means in the reservoir area to achieve sediment discharge. It usually sets the threshold directly, is highly empirical, is applicable to specific reservoirs, and is not universal. In addition, it does not take into account that each backwater area in the river-type reservoir has different water and sand movement characteristics. For example, in a flood process, the operating water level Z has the most serious siltation in the perennial backwater area, but has little effect on the siltation in the variable backwater area; for another flood process entering the reservoir, the same operating water level Z causes the most serious siltation in the variable backwater area, and has little effect on the siltation in the variable backwater area. Therefore, when controlling the siltation of river-type reservoirs, it is necessary to consider the siltation conditions in different areas of the reservoir for targeted regulation.

[0057] In this embodiment, a reservoir sedimentation control method is provided. Figure 1 is a flow chart of a reservoir sedimentation control method according to an embodiment of the present invention. Figure 1 As shown, the process includes the following steps:

[0058] S11, obtaining the water-sediment movement process and the water-sediment movement model within the target time period, and determining the reservoir sediment discharge ratio, the variable backwater area siltation amount, and the perennial backwater area siltation amount based on the water-sediment movement process and the water-sediment movement model.

[0059] The target time period refers to the next few days. According to actual needs, the specific number of days in the next few days is not limited. The water and sand movement process of the next few days is provided on a daily basis by the forecast department. The water and sand movement process includes the movement process of the flood peak and the movement process of the sand peak. A mathematical model of water and sand movement in the reservoir area is established based on the water and sand movement to obtain a water and sand movement model. The water and sand movement model is used to calculate the water and sand process and the change process of the water level, flow rate, and sediment content of each river section corresponding to the reservoir water level over time. The water and sand movement model can be a one-dimensional mathematical model, and the specific establishment method is not limited. The reservoir sediment discharge ratio is obtained by calculating the ratio of the amount of sediment leaving the reservoir to the amount of sediment entering the reservoir. The amount of sediment entering the reservoir can be calculated from the flow rate and sediment content obtained during the water and sand movement process, and the amount of sediment leaving the reservoir can be simulated by the water and sand movement model. The amount of siltation in the variable backwater area is obtained by calculating the difference between the amount of sediment entering the reservoir and the amount of sediment discharged at the boundary section of the backwater area. The amount of siltation in the perennial backwater area is obtained by calculating the difference between the amount of sediment discharged at the boundary section of the backwater area and the amount of sediment discharged from the reservoir. The sediment discharge at the boundary section of the backwater area can be simulated by a water-sediment movement model.

[0060] S12, based on the water-sediment movement model, reservoir sediment discharge ratio, variable backwater area sedimentation amount and perennial backwater area sedimentation amount, determine the sediment discharge ratio index, variable backwater area sedimentation proportion index and perennial backwater area sedimentation proportion index.

[0061] The sediment discharge ratio index can be used to indicate the proportion of sediment discharge ratio attenuation under different operating water levels, the variable backwater area siltation ratio index can be used to indicate the proportion of variable backwater area siltation under different operating water levels, and the perennial backwater area siltation ratio index can be used to indicate the proportion of perennial backwater area siltation under different operating water levels. After obtaining the water-sediment movement model, the changes corresponding to different operating water levels, the perennial backwater area siltation ratio, and the variable backwater area siltation ratio are statistically analyzed to establish the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index.

[0062] S13, respectively calculating the functions corresponding to the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index and the function monotonicity, and determining the avoidance water level corresponding to each index based on the function monotonicity.

[0063] For the three indicators, the calculated values ​​of the indicators corresponding to the step-by-step increase of the operating water level from the flood limit water level to the flood control high water level are statistically analyzed, and the functions corresponding to the three indicators are constructed respectively, namely, the function of the sediment discharge ratio indicator on the water level, the function of the variable backwater area siltation ratio indicator on the water level, and the function of the perennial backwater area siltation ratio indicator on the water level. The three functions are analyzed for monotonicity respectively, and the extreme value points, slope value points, etc. are taken as the avoidance water level. Each function has at least one avoidance water level, and there may be multiple avoidance water levels according to the actual situation. The avoidance water level is the minimum value higher than the current operating water level, and the water level interval represents the interval formed by the operating water level and the avoidance water level.

[0064] When performing water level control, a suitable avoidance water level may be selected according to the actual operating water level. For example, the operating water level may be controlled not to be higher than the adjacent avoidance water level, and may be between the actual operating water level and the avoidance water level.

[0065] S14, obtaining a target siltation control demand, and based on the target siltation control demand, determining a target avoidance water level corresponding to the target siltation control demand from the avoidance water levels corresponding to the indicators, and determining a target water level interval based on the target avoidance water level to control reservoir siltation.

[0066] The target siltation reduction control demand can be one or more. In this embodiment, a single siltation reduction control demand is taken as an example, and the corresponding index is selected for analysis according to the actual target siltation reduction control demand. For example: when it is necessary to control the siltation of the variable backwater area, according to the current reservoir operating water level, the smallest avoidance water level higher than the current reservoir operating water level is selected from the avoidance water levels corresponding to the variable backwater area siltation ratio index. The avoidance water level is the target avoidance water level, and the water level interval formed by the current reservoir operating water level or the next level operating water level of the current operating water level and the target avoidance water level is regarded as the target water level interval. By controlling the water level, it is ensured that the reservoir controls the siltation of the target area within the target water level interval as little as possible. If there are multiple solutions for the target avoidance water level, all of them must be avoided.

[0067] The reservoir siltation control method provided by the embodiment of the present invention determines the reservoir sediment discharge ratio, the variable backwater area sedimentation amount, and the perennial backwater area sedimentation amount, and then establishes the sediment discharge ratio index, the variable backwater area sedimentation ratio index, and the perennial backwater area sedimentation ratio index. The monotonicity of the function corresponding to each index is analyzed to determine the avoidance water level corresponding to each index. Different indicators are determined for different areas and siltation reduction needs. When siltation control is required through water level regulation, the corresponding target avoidance water level is selected and the target water level range is determined according to the target siltation control demand, to ensure that the reservoir can control the sedimentation in the target area as little as possible under the target avoidance water level, and will not cause significant sedimentation.

[0068] In some embodiments, Figure 1 S11 includes the following steps:

[0069] S21, determine the flow and sediment content of the river section entering the reservoir according to the water and sediment movement process, and calculate the amount of sediment entering the reservoir based on the flow and sediment content of the river section entering the reservoir.

[0070] The flow and sediment content of the river section entering the reservoir are counted during the movement of water and sediment. The amount of sediment entering the reservoir is calculated based on the obtained flow and sediment content. The calculation method of the amount of sediment entering the reservoir is as follows:

[0071] Flow rate (m 3 / s)×sand content (kg / m 3 )×time(s)=amount of sand entering the warehouse(kg)

[0072] S22, determining the flow and sediment content of the outlet river section based on the water and sediment movement model, and calculating the amount of sediment in the reservoir based on the flow and sediment content of the outlet river section.

[0073] The water-sediment movement process, reservoir water level and topographic parameters of the river section entering the reservoir are input into the water-sediment movement model to simulate the water-sediment movement. The dynamic change process of the water level, flow rate and sediment content of each section of the river in the reservoir area over time is simulated, and the sediment volume out of the reservoir is calculated based on this. The calculation method of the sediment volume out of the reservoir is as follows:

[0074] Flow rate (m 3 / s)×sand content (kg / m 3 )×time(s)=sand quantity discharged from the warehouse(kg)

[0075] S23, determining the reservoir sediment discharge ratio according to the ratio of the outgoing sediment volume to the incoming sediment volume.

[0076] The calculation method of reservoir sediment discharge ratio P is: P = sediment outflow ÷ sediment inflow.

[0077] S24, determining the flow and sediment content of the backwater zone boundary section based on the water-sediment movement model, and determining the variable backwater zone siltation amount and the perennial backwater zone siltation amount based on the flow and sediment content of the backwater zone boundary section. Specifically, step S24 includes the following steps:

[0078] S241, calculating the sediment transport volume of the backwater zone boundary section based on the flow and sediment content of the backwater zone boundary section.

[0079] The flow and sediment content of the backwater zone boundary section are obtained from the water-sediment movement model, and the sediment transport of the backwater zone boundary section is calculated. The calculation method of the sediment transport of the backwater zone boundary section is as follows:

[0080] Flow rate at the junction section of the backwater area (m 3 / s)×Sediment content at the junction of backwater area (kg / m 3 )×time (s)=sediment transport volume at the junction section of the backwater area (kg)

[0081] S242, determine the amount of sedimentation in the variable backwater area based on the difference between the amount of sediment entering the reservoir and the amount of sediment discharged at the cross section of the backwater area.

[0082] The calculation method for the variable backwater area siltation volume B is: B = amount of sediment entering the reservoir - amount of sediment discharged at the backwater area boundary section.

[0083] S243, determine the annual siltation volume of the backwater area based on the difference between the sediment discharge volume at the backwater area boundary section and the sediment discharge volume.

[0084] The calculation method for the annual siltation volume C in the backwater area is: C = sediment transport volume at the backwater area boundary section - sediment discharge volume.

[0085] In some embodiments, the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index are determined according to the following formula:

[0086]

[0087]

[0088] Among them, i represents the i-th level of operating water level, and the reservoir operating water level can be divided into n levels, i=1 represents the flood limit water level, and i=n represents the flood control high water level; P represents the reservoir sediment discharge ratio; B represents the siltation in the variable backwater area; C represents the siltation in the perennial backwater area; I P Indicates the sediment discharge ratio index; I B Indicates the proportion of sedimentation in the changing backwater area; I C It indicates the proportion of siltation in the perennial backwater area; the subscript min indicates that the reservoir is at the flood limit water level; the subscript max indicates that the reservoir is at the flood control high water level; It represents the ratio of the sediment discharge ratio attenuation at the i-th operating water level to the total attenuation from the water level min to max; It indicates the ratio of the sedimentation in the variable backwater area under the i-th operating water level to the total sedimentation from the water level min to max; It represents the ratio of the annual siltation in the backwater area at the i-th operating water level to the total siltation when the water level ranges from min to max.

[0089] In some embodiments, Figure 1 S13 includes the following steps:

[0090] S31, dividing at least one operating water level.

[0091] The operating water level can be set according to the range from the flood limit water level to the flood control high water level of the reservoir, and one or more operating water levels can be set according to actual needs.

[0092] Taking the flood limit water level as 100 meters and the flood control high water level as 150 meters, 10 operating water levels are set as an example, with an interval of 5 meters between each operating water level.

[0093] S32, respectively calculating the sediment discharge ratio index calculation value corresponding to the operating water level, the variable backwater area siltation ratio index calculation value, and the perennial backwater area siltation ratio index calculation value.

[0094] The operating water level is substituted into the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index to obtain the corresponding calculated values ​​of the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index.

[0095] S33, drawing function graphs based on the calculated values ​​of the sediment discharge ratio index, the calculated values ​​of the variable backwater area siltation ratio index, and the calculated values ​​of the perennial backwater area siltation ratio index, respectively, to determine the sediment discharge ratio index function, the variable backwater area siltation ratio index function, and the perennial backwater area siltation ratio index function.

[0096] According to the calculated values ​​of indicators corresponding to each operating water level, the function graphs of the sediment discharge ratio indicator function, the variable backwater area siltation ratio indicator function and the perennial backwater area siltation ratio indicator function are drawn.

[0097] S34, respectively analyzing the monotonicity of the sediment discharge ratio index function, the variable backwater area siltation ratio index function, and the perennial backwater area siltation ratio index function, and determining the avoidance water level corresponding to each index based on the monotonicity.

[0098] In this embodiment, a monotonicity analysis is performed based on the obtained function graph, and the maximum point, the positive slope value, and the high point can be taken as the avoidance water level.

[0099] In some embodiments, when the target siltation control requirement includes more than one indicator, Figure 1 The reservoir siltation control method further comprises the following steps:

[0100] S41, obtaining the weight coefficient corresponding to each indicator.

[0101] The target silt reduction control demand may involve multiple indicators, taking into account 2 to 3 of the sediment discharge ratio indicator, the variable backwater area sedimentation ratio indicator, and the perennial backwater area sedimentation ratio indicator. Different weight coefficients can be set according to the importance of each indicator. For example, if you do not want the reservoir sediment discharge ratio to be significantly reduced, you can give the sediment discharge ratio indicator a greater weight; if you want to alleviate the problem of siltation at the tail of the reservoir that hinders navigation, you need to minimize the sedimentation in the variable backwater area, and you need to give the variable backwater area sedimentation ratio indicator a greater weight; if you want to minimize the sedimentation problem in front of the dam, you need to minimize the sedimentation in the perennial backwater area, and you need to give the perennial backwater area sedimentation ratio indicator a greater weight.

[0102] S42, determining a weighted index for reducing siltation based on the weight coefficients corresponding to the various indexes.

[0103] In some embodiments, the weighted index for reducing siltation is determined according to the following formula:

[0104]

[0105] Among them, W i represents the weighted index of silt reduction, i represents the i-th level operating water level, I P Indicates the sediment discharge ratio index, I B Indicates the proportion of sedimentation in the changing backwater area, I C It represents the proportion of siltation in the backwater area throughout the year, a 1 , a 2 , a 3 are weight coefficients, where a 1 +a 2 +a 3 =1.

[0106] S43, dividing at least one target water level, and calculating the siltation reduction weighted index calculation values ​​corresponding to different target water levels to determine a function of the target water level with respect to the siltation reduction weighted index.

[0107] The target water level can be set according to the flood limit water level to the flood control high water level of the reservoir, and one or more target water levels can be set according to actual needs. Substitute all target water levels into the siltation reduction weighted index to obtain the calculated values ​​of the siltation reduction weighted index corresponding to all target water levels, and establish the relationship between the siltation reduction weighted index and the water level, that is, the function of the target water level on the siltation reduction weighted index.

[0108] S44, analyzing the monotonicity of the function of the weighted siltation reduction index to determine the avoided water level of the weighted siltation reduction index.

[0109] Monotonicity analysis can be performed by drawing a function graph of the weighted index of silt reduction. The method of analyzing monotonicity is not limited here. The maximum value point and the positive slope value and high point are taken as the avoidance water level.

[0110] S45, determining a target avoidance water level corresponding to the target siltation control demand from the avoidance water level of the siltation reduction weighted index based on the target siltation control demand, and determining a target water level interval based on the target avoidance water level.

[0111] According to the actual operating water level in the target siltation control demand, the smallest avoidance water level that is higher than the current reservoir operating water level is selected from the avoidance water levels of the siltation reduction weighted index. This avoidance water level is the target avoidance water level, and the water level interval formed by the current reservoir operating water level or the next level operating water level of the current operating water level and the target avoidance water level is regarded as the target water level interval.

[0112] In the embodiment of the present invention, a river-type reservoir area is selected as the research object. The total length of the reservoir area is 600km, of which 0-200km from the upstream is the variable backwater area, and 200-600km is the perennial backwater area. The current water level of the reservoir is 121.5m. Facing the upcoming flood process, the reservoir will intercept and store water, resulting in an inevitable rise in the reservoir water level. However, in order to control the increase in siltation in the reservoir area as little as possible, it is necessary to select the avoidance water level and water level interval according to the different regional siltation control requirements. The specific steps are as follows:

[0113] (1) Obtain the change process of the predicted flood peak and sand peak over time in the next week, that is, the water and sand movement process, such as Figure 2 shown.

[0114] (2) Establish a reservoir water and sediment movement model, input the water and sediment movement process in step (1), the different reservoir water levels, and the reservoir cross-sectional terrain parameters. According to the reservoir flood limit water level of 100m to the flood control high water level of 150m, set 10 operating water levels, each with a 5m interval. Under the same reservoir water and sediment process, input 10 water levels into the model to simulate the water flow and sediment movement in the reservoir area, so a total of 10 simulations are performed.

[0115] (3) The amount of sediment entering the reservoir is calculated from the flow and sediment content obtained during the water-sediment movement process, and the amount of sediment leaving the reservoir is simulated by the water-sediment movement model. The ratio of the amount of sediment leaving the reservoir to the amount of sediment entering the reservoir is used to obtain the sediment discharge ratio. The flow and sediment content of the cross-section at the junction of the backwater area are simulated by the water-sediment movement model, and the sediment transport at the cross-section at the junction of the backwater area is calculated. The difference between the amount of sediment entering the reservoir and the amount of sediment transported at the cross-section at the junction of the backwater area is then calculated to obtain the variable amount of sedimentation in the backwater area. The difference between the amount of sediment transported at the cross-section at the junction of the backwater area and the amount of sediment leaving the reservoir is calculated to obtain the perennial amount of sedimentation in the backwater area.

[0116] Establish the sediment discharge ratio index, the variable backwater area sedimentation ratio index, and the perennial backwater area sedimentation ratio index, as follows:

[0117]

[0118]

[0119] Among them, i represents the i-th level of operating water level, and the reservoir operating water level can be divided into n levels, i=1 represents the flood limit water level, and i=n represents the flood control high water level; P represents the reservoir sediment discharge ratio; B represents the siltation in the variable backwater area; C represents the siltation in the perennial backwater area; I P Indicates the sediment discharge ratio index; I B Indicates the proportion of sedimentation in the changing backwater area; I C It indicates the proportion of siltation in the perennial backwater area; the subscript min indicates that the reservoir is at the flood limit water level; the subscript max indicates that the reservoir is at the flood control high water level; It represents the ratio of the sediment discharge ratio attenuation at the i-th operating water level to the total attenuation from the water level min to max; It indicates the ratio of the sedimentation in the variable backwater area under the i-th operating water level to the total sedimentation from the water level min to max; It represents the ratio of the annual siltation in the backwater area at the i-th operating water level to the total siltation when the water level ranges from min to max.

[0120] (4) For I P ,I B ,I C The three indicators are statistically analyzed as the operating water level gradually rises from the flood limit water level to the flood control high water level. The three functional relationships are obtained by drawing: P =(I p ), f B = Z(I B ), f C =(I C ). The functions obtained are as follows Figure 3 , Figure 4 , Figure 5 According to the monotonicity of the function, the maximum point, the positive slope value, and the high point are taken as the avoidance water level.

[0121] The avoidance water levels corresponding to the sediment discharge ratio index are: 100m, 105m, 110m, 125m, and 140m.

[0122] The avoidance water levels corresponding to the changing backwater area siltation ratio indicators are: 100m, 105m, 115m, and 140m.

[0123] The avoidance water levels corresponding to the annual backwater area siltation ratio indicators are: 100m, 120m, 135m, and 145m.

[0124] (5) When the target sedimentation control demand includes one indicator, for a single sedimentation reduction demand, only one indicator needs to be considered. According to the current reservoir water level, the reservoir operating water level is controlled to be no higher than the adjacent avoidance water level. If there are multiple solutions for the avoidance water level, all of them must be avoided. Since the current reservoir water level is 121.5m, if I P To control the index, the water level of the reservoir should be controlled to be no higher than the nearby avoidance water level of 125m when the reservoir is storing flood water. That is, the target avoidance water level is 125m, and the water level interval of the current operating water level of 121.5m to 125m is the target water level interval. B As the control index, the reservoir water level is controlled to be no higher than the adjacent avoidance water level of 140m, that is, the target avoidance water level is 140m, and the water level interval from the current operating water level of 121.5m to 140m is the target water level interval; if I CTo control the index, the reservoir water level must be controlled not higher than 135m, that is, the target avoidance water level is 135m, and the current operating water level interval of 121.5m to 135m is the target water level interval. If flood storage will cause the reservoir to exceed the above water level limit, the reservoir water level can be controlled to exceed the avoidance water level and dispatched within the next water level interval as much as possible.

[0125] (6) When the target sedimentation control requirement includes more than one indicator, a weighted sedimentation reduction indicator W is established:

[0126]

[0127] Among them, W i represents the weighted index of silt reduction, i represents the i-th level operating water level, I P Indicates the sediment discharge ratio index; I B Indicates the proportion of sedimentation in the changing backwater area; I C It represents the proportion of siltation in the backwater area throughout the year, a 1 、a 2 、a 3 are weight coefficients, where a 1 +a 2 +a 3 =1.

[0128] If the sedimentation control requirements of all three indicators are considered at the same time, for example, if a 1 =0.4, a 2 =0.3, a 3 = 0.3, the multi-objective problem can be transformed into a weighted single-objective weighted index W for reducing siltation. P ,I B and I C Substituting the calculation result into the above formula, the relationship between the target water level and the weighted index of silt reduction is finally obtained as follows: Figure 6 As shown, it should be noted that the target water level and the operating water level mentioned above are the same here. The same method as in step (4) above is used to determine the avoidance water level corresponding to the weighted index for silt reduction, and the avoidance water levels are: 100m, 105m, and 140m. When the current reservoir water level is 121.5m, the avoidance water level is 140m. The water level interval of the current operating water level of 121.5m to 140m is the target water level interval. The reservoir can intercept floods within this target water level interval to meet the silt reduction demand under the weighted target.

[0129] In this embodiment, a reservoir siltation control device is also provided, which is used to implement the above-mentioned embodiments and implementation methods, and the descriptions that have been made will not be repeated. As used below, the term "module" can implement a combination of software and / or hardware of a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceived.

[0130] This embodiment provides a reservoir sedimentation control device, such as Figure 7 As shown, including:

[0131] A data acquisition module is used to acquire the water-sediment movement process and the water-sediment movement model within the target time period, and determine the reservoir sediment discharge ratio, the variable backwater area siltation, and the perennial backwater area siltation based on the water-sediment movement process and the water-sediment movement model;

[0132] An index determination module, for determining a sediment discharge ratio index, a sediment accumulation ratio index of a variable backwater area, and a sediment accumulation ratio index of a perennial backwater area based on the water-sediment movement model, the reservoir sediment discharge ratio, the sediment accumulation amount of a variable backwater area, and the sediment accumulation amount of a perennial backwater area;

[0133] A water level calculation module, used to calculate the functions corresponding to the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index and the function monotonicity, and determine the avoidance water level corresponding to each index based on the function monotonicity;

[0134] The target water level determination module is used to obtain the target siltation control demand, and based on the target siltation control demand, determine the target avoidance water level corresponding to the target siltation control demand from the avoidance water levels corresponding to each indicator, and determine the target water level interval based on the target avoidance water level to control the siltation of the reservoir.

[0135] In some embodiments, the data acquisition module includes:

[0136] An inflow sediment volume determination unit, used to determine the flow and sediment content of the inflow river section according to the water and sediment movement process, and calculate the inflow sediment volume based on the flow and sediment content of the inflow river section;

[0137] A discharge sediment quantity determination unit, used to determine the flow and sediment content of the discharge river section based on the water-sediment movement model, and calculate the discharge sediment quantity based on the flow and sediment content of the discharge river section;

[0138] A sediment discharge ratio determination unit, used to determine the sediment discharge ratio of the reservoir according to the ratio of the outgoing sediment volume to the incoming sediment volume;

[0139] The siltation amount calculation unit is used to determine the flow and sediment content of the backwater area boundary section based on the water and sediment movement model, and to determine the siltation amount of the variable backwater area and the siltation amount of the perennial backwater area based on the flow and sediment content of the backwater area boundary section.

[0140] In some embodiments, the sedimentation amount calculation unit includes:

[0141] A sediment transport determination subunit, used to calculate the sediment transport at the junction section of the backwater area based on the flow and sediment content at the junction section of the backwater area;

[0142] The first sedimentation amount determination subunit is used to determine the sedimentation amount of the variable backwater area according to the difference between the amount of sediment entering the reservoir and the sediment transport amount of the cross section of the backwater area;

[0143] The second siltation amount determination subunit is used to determine the annual siltation amount of the backwater area according to the difference between the sediment transport volume of the backwater area boundary section and the sediment discharge volume.

[0144] In some embodiments, the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index are determined according to the following formula:

[0145]

[0146]

[0147] Among them, i represents the i-th level of operating water level, and the reservoir operating water level can be divided into n levels, i=1 represents the flood limit water level, and i=n represents the flood control high water level; P represents the reservoir sediment discharge ratio; B represents the siltation in the variable backwater area; C represents the siltation in the perennial backwater area; I P Indicates the sediment discharge ratio index; I B Indicates the proportion of sedimentation in the changing backwater area; I C It indicates the proportion of siltation in the perennial backwater area; the subscript min indicates that the reservoir is at the flood limit water level; the subscript max indicates that the reservoir is at the flood control high water level; It represents the ratio of the sediment discharge ratio attenuation at the i-th operating water level to the total attenuation from the water level min to max; It indicates the ratio of the sedimentation in the variable backwater area under the i-th operating water level to the total sedimentation from the water level min to max; It represents the ratio of the annual siltation in the backwater area at the i-th operating water level to the total siltation when the water level ranges from min to max.

[0148] In some embodiments, the water level calculation module includes:

[0149] A water level division unit, used for dividing at least one operating water level;

[0150] An index calculation unit is used to calculate the sediment discharge ratio index calculation value corresponding to the operating water level, the variable backwater area siltation ratio index calculation value, and the perennial backwater area siltation ratio index calculation value;

[0151] An index function determination unit is used to draw function images based on the calculated value of the sediment discharge ratio index, the calculated value of the variable backwater area siltation ratio index, and the calculated value of the perennial backwater area siltation ratio index, so as to determine the sediment discharge ratio index function, the variable backwater area siltation ratio index function, and the perennial backwater area siltation ratio index function;

[0152] The water level determination unit is used to analyze the monotonicity of the sediment discharge ratio index function, the variable backwater area siltation ratio index function and the perennial backwater area siltation ratio index function respectively, and determine the avoidance water level corresponding to each index based on the monotonicity.

[0153] In some embodiments, when the target siltation control requirement includes more than one indicator, the apparatus further includes:

[0154] A weight determination module is used to obtain the weight coefficient corresponding to each indicator;

[0155] A weighted index determination module, used to determine a weighted index for reducing siltation based on weight coefficients corresponding to the indicators;

[0156] A weighted function determination module, used to divide at least one target water level, and calculate the siltation reduction weighted index calculation values ​​corresponding to different target water levels, so as to determine the function of the target water level with respect to the siltation reduction weighted index;

[0157] A weighted function analysis module, used for analyzing the monotonicity of the function of the siltation reduction weighted index to determine the avoided water level of the siltation reduction weighted index;

[0158] A weighted water level determination module is used to determine a target avoidance water level corresponding to the target siltation control demand from the avoidance water level of the siltation reduction weighted index based on the target siltation control demand, and determine a target water level interval based on the target avoidance water level.

[0159] In some embodiments, the weighted index for reducing siltation is determined according to the following formula:

[0160]

[0161] Among them, W i represents the weighted index of silt reduction, i represents the i-th level operating water level, I P Indicates the sediment discharge ratio index; I B Indicates the proportion of sedimentation in the changing backwater area; I C It represents the proportion of siltation in the backwater area throughout the year, a 1 、a 2 、a 3are weight coefficients, where a 1 +a 2 +a 3 =1.

[0162] The reservoir siltation control device in this embodiment is presented in the form of a functional unit, where the unit refers to an ASIC circuit, a processor and memory that executes one or more software or fixed programs, and / or other devices that can provide the above functions.

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

[0164] The embodiment of the present invention also provides an electronic device having the above Figure 8 Reservoir siltation control device shown.

[0165] See also Figure 8 , Figure 8 is a schematic diagram of the structure of an electronic device provided by an optional embodiment of the present invention, such as Figure 8 As shown, the electronic device may include: at least one processor 601, such as a CPU (Central Processing Unit), at least one communication interface 603, a memory 604, and at least one communication bus 602. The communication bus 602 is used to realize the connection and communication between these components. The communication interface 603 may include a display screen (Display), a keyboard (Keyboard), and the optional communication interface 603 may also include a standard wired interface and a wireless interface. The memory 604 may be a high-speed RAM memory (Random Access Memory) or a non-volatile memory (non-volatile memory), such as at least one disk storage. The memory 604 may optionally be at least one storage device located away from the aforementioned processor 601. The processor 601 may be combined with Figure 5 In the described device, the memory 604 stores an application program, and the processor 601 calls the program code stored in the memory 604 to execute any of the above method steps.

[0166] The communication bus 602 may be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus. The communication bus 602 may be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 8Only one thick line is used in the diagram, but this does not mean that there is only one bus or only one type of bus.

[0167] Among them, the memory 604 may include a volatile memory (English: volatile memory), such as a random access memory (English: random-access memory, abbreviated: RAM); the memory may also include a non-volatile memory (English: non-volatile memory), such as a flash memory (English: flash memory), a hard disk drive (English: hard disk drive, abbreviated: HDD) or a solid-state drive (English: solid-state drive, abbreviated: SSD); the memory 604 may also include a combination of the above types of memory.

[0168] The processor 601 may be a central processing unit (CPU), a network processor (NP), or a combination of a CPU and a NP.

[0169] The processor 601 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit (ASIC), a programmable logic device (PLD), or a combination thereof. The PLD may be a complex programmable logic device (CPLD), a field-programmable gate array (FPGA), a generic array logic (GAL), or any combination thereof.

[0170] Optionally, the memory 604 is also used to store program instructions. The processor 601 can call the program instructions to implement the reservoir siltation control method shown in the embodiment of the present application.

[0171] The embodiment of the present invention further provides a non-transitory computer storage medium, wherein the computer storage medium stores computer executable instructions, and the computer executable instructions can execute the reservoir siltation control method in any of the above method embodiments. The storage medium may be a disk, an optical disk, a read-only memory (ROM), a random access memory (RAM), a flash memory, a hard disk (HDD) or a solid-state drive (SSD), etc.; the storage medium may also include a combination of the above types of memory.

[0172] Although the embodiments of the present invention have been described in conjunction with 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, and such modifications and variations are all within the scope defined by the appended claims.

Claims

1. A reservoir sedimentation control method, characterized in that: include: The water-sediment movement process and the water-sediment movement model within the target time period are obtained, and based on the water-sediment movement process and the water-sediment movement model, the reservoir sediment discharge ratio, the variable backwater area siltation, and the perennial backwater area siltation are determined, including: determining the flow and sediment content of the inlet river section according to the water-sediment movement process, and calculating the inlet sediment volume based on the flow and sediment content of the inlet river section; determining the flow and sediment content of the outlet river section based on the water-sediment movement model, and calculating the outlet sediment volume based on the flow and sediment content of the outlet river section; determining the reservoir sediment discharge ratio according to the ratio of the outlet sediment volume to the inlet sediment volume; and determining the sediment discharge ratio of the reservoir based on the water-sediment movement model. The method comprises: determining the flow rate and sediment content of the cross-section of the backwater area according to the flow rate and sediment content of the cross-section of the backwater area, and determining the variable backwater area siltation amount and the perennial backwater area siltation amount based on the flow rate and sediment content of the cross-section of the backwater area; determining the variable backwater area siltation amount and the perennial backwater area siltation amount based on the flow rate and sediment content of the cross-section of the backwater area, including: calculating the sediment transport amount of the cross-section of the backwater area according to the flow rate and sediment content of the cross-section of the backwater area; determining the variable backwater area siltation amount according to the difference between the inflow sediment amount and the sediment transport amount of the cross-section of the backwater area; determining the perennial backwater area siltation amount according to the difference between the sediment transport amount of the cross-section of the backwater area and the outflow sediment amount; Based on the water-sediment movement model, the reservoir sediment discharge ratio, the variable backwater area sedimentation amount and the perennial backwater area sedimentation amount, the sediment discharge ratio index, the variable backwater area sedimentation proportion index and the perennial backwater area sedimentation proportion index are determined; the sediment discharge ratio index, the variable backwater area sedimentation proportion index and the perennial backwater area sedimentation proportion index are determined according to the following formula: Among them, i represents the i-th level of operating water level, and the reservoir operating water level is divided into n levels. i=1 represents the flood limit water level, and i=n represents the flood control high water level; P represents the reservoir sediment discharge ratio; B represents the siltation in the variable backwater area; C represents the siltation in the perennial backwater area; I P Indicates the sediment discharge ratio index; I B Indicates the proportion of sedimentation in the changing backwater area; I C It indicates the proportion of siltation in the backwater area throughout the year; the subscript min indicates that the reservoir is at the flood limit water level; the subscript max indicates that the reservoir is at the high flood control water level; It represents the ratio of the sediment discharge ratio attenuation at the i-th operating water level to the total attenuation from the water level min to max; It represents the ratio of the siltation in the variable backwater area under the i-th operating water level to the total siltation from the water level min to max; It represents the ratio of the annual siltation in the backwater area under the i-th operating water level to the total siltation from the water level min to max; The functions and function monotonicity corresponding to the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index are calculated respectively, and the avoidance water level corresponding to each index is determined based on the function monotonicity; the functions and function monotonicity corresponding to the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index are calculated respectively, and the avoidance water level corresponding to each index is determined based on the function monotonicity, including: dividing at least one operating water level; calculating the calculated value of the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index respectively corresponding to the operating water level The calculated value of the proportion of siltation in the backwater area and the calculated value of the proportion of siltation in the perennial backwater area; the function graphs are drawn based on the calculated value of the sediment discharge ratio index, the calculated value of the proportion of siltation in the variable backwater area and the calculated value of the proportion of siltation in the perennial backwater area, so as to determine the sediment discharge ratio index function, the variable backwater area siltation ratio index function and the perennial backwater area siltation ratio index function; the monotonicity of the sediment discharge ratio index function, the variable backwater area siltation ratio index function and the perennial backwater area siltation ratio index function are analyzed respectively, and the avoidance water level corresponding to each index is determined based on the monotonicity; Obtain a target siltation control demand, and based on the target siltation control demand, determine a target avoidance water level corresponding to the target siltation control demand from the avoidance water levels corresponding to each indicator, and determine a target water level interval based on the target avoidance water level to control reservoir siltation.

2. The method according to claim 1, characterized in that When the target siltation control requirement includes more than one indicator, the method further includes: Get the weight coefficient corresponding to each indicator; Determine a weighted index for reducing siltation based on the weight coefficients corresponding to the various indicators; Dividing at least one target water level, and calculating the calculated values ​​of the siltation reduction weighted index corresponding to different target water levels, so as to determine a function of the target water level with respect to the siltation reduction weighted index; Analyzing the monotonicity of the function of the weighted siltation reduction index to determine the avoided water level of the weighted siltation reduction index; Based on the target siltation control demand, a target avoidance water level corresponding to the target siltation control demand is determined from the avoidance water level of the siltation reduction weighted index, and a target water level interval is determined based on the target avoidance water level.

3. The method according to claim 2, characterized in that The weighted index for reducing siltation is determined according to the following formula: Among them, W i represents the weighted index of silt reduction, i represents the i-th level operating water level, I P Indicates the sediment discharge ratio index; I B Indicates the proportion of sedimentation in the changing backwater area; I C It represents the proportion of siltation in the perennial backwater area. a1, a2, and a3 are all weight coefficients, among which a1+a2+a3=1.

4. A reservoir sedimentation control device, characterized in that: include: A data acquisition module is used to acquire the water-sediment movement process and the water-sediment movement model within the target time period, and determine the reservoir sediment discharge ratio, the variable backwater area siltation, and the perennial backwater area siltation based on the water-sediment movement process and the water-sediment movement model; The data acquisition module includes: an inflow sand volume determination unit, which is used to determine the flow and sediment content of the inflow river section according to the water-sand movement process, and calculate the inflow sand volume based on the flow and sediment content of the inflow river section; an outflow sand volume determination unit, which is used to determine the flow and sediment content of the outflow river section based on the water-sand movement model, and calculate the outflow sand volume based on the flow and sediment content of the outflow river section; a sediment discharge ratio determination unit, which is used to determine the reservoir sediment discharge ratio according to the ratio of the outflow sand volume to the inflow sand volume; a siltation volume calculation unit, which is used to determine the flow and sediment content of the backwater zone boundary section based on the water-sand movement model, and determine the siltation volume of the variable backwater zone and the siltation volume of the perennial backwater zone based on the flow and sediment content of the backwater zone boundary section; The sedimentation amount calculation unit includes: a sedimentation amount determination subunit, which is used to calculate the sedimentation amount of the backwater area boundary section based on the flow and sediment content of the backwater area boundary section; a first sedimentation amount determination subunit, which is used to determine the variable backwater area sedimentation amount according to the difference between the inflow sediment amount and the sedimentation amount of the backwater area boundary section; a second sedimentation amount determination subunit, which is used to determine the perennial backwater area sedimentation amount according to the difference between the sedimentation amount of the backwater area boundary section and the outflow sediment amount; An index determination module is used to determine the sediment discharge ratio index, the sedimentation ratio index of the variable backwater area, and the sedimentation ratio index of the perennial backwater area based on the water-sediment movement model, the reservoir sediment discharge ratio, the sedimentation amount of the variable backwater area, and the sedimentation amount of the perennial backwater area; the sediment discharge ratio index, the sedimentation ratio index of the variable backwater area, and the sedimentation ratio index of the perennial backwater area are determined according to the following formula: Among them, i represents the i-th level of operating water level, and the reservoir operating water level is divided into n levels. i=1 represents the flood limit water level, and i=n represents the flood control high water level; P represents the reservoir sediment discharge ratio; B represents the siltation in the variable backwater area; C represents the siltation in the perennial backwater area; I P Indicates the sediment discharge ratio index; I B Indicates the proportion of sedimentation in the changing backwater area; I C It indicates the proportion of siltation in the backwater area throughout the year; the subscript min indicates that the reservoir is at the flood limit water level; the subscript max indicates that the reservoir is at the high flood control water level; It represents the ratio of the sediment discharge ratio attenuation at the i-th operating water level to the total attenuation from the water level min to max; It represents the ratio of the siltation in the variable backwater area under the i-th operating water level to the total siltation from the water level min to max; It represents the ratio of the annual siltation in the backwater area under the i-th operating water level to the total siltation from the water level min to max; A water level calculation module, used to calculate the functions corresponding to the sediment discharge ratio index, the variable backwater area siltation ratio index, and the perennial backwater area siltation ratio index and the function monotonicity, and determine the avoidance water level corresponding to each index based on the function monotonicity; The water level calculation module includes: a water level division unit, used to divide at least one operating water level; An index calculation unit is used to respectively calculate the sediment discharge ratio index calculation value, the variable backwater area siltation ratio index calculation value, and the perennial backwater area siltation ratio index calculation value corresponding to the operating water level; an index function determination unit is used to draw function graphs based on the sediment discharge ratio index calculation value, the variable backwater area siltation ratio index calculation value, and the perennial backwater area siltation ratio index calculation value, so as to determine the sediment discharge ratio index function, the variable backwater area siltation ratio index function, and the perennial backwater area siltation ratio index function; a water level determination unit is used to respectively analyze the monotonicity of the sediment discharge ratio index function, the variable backwater area siltation ratio index function, and the perennial backwater area siltation ratio index function, and determine the avoidance water level corresponding to each index based on the monotonicity; a target water level determination module is used to obtain a target siltation control demand, and determine the target avoidance water level corresponding to the target siltation control demand from the avoidance water levels corresponding to each index based on the target siltation control demand, and determine the target water level interval based on the target avoidance water level to control reservoir siltation.

5. An electronic 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 reservoir siltation control method according to any one of claims 1 to 3 by executing the computer instructions.

6. 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 reservoir siltation control method according to any one of claims 1 to 3.

Images